非结构网格催化壁面的隐式处理方法

推进技术 ›› 2018, Vol. 39 ›› Issue (2): 251-260.

非结构网格催化壁面的隐式处理方法

原志超,黄世璋,高效伟

大连理工大学工业装备结构分析国家重点实验室，航空航天学院，辽宁大连 116024,大连理工大学工业装备结构分析国家重点实验室，航空航天学院，辽宁大连 116024,大连理工大学工业装备结构分析国家重点实验室，航空航天学院，辽宁大连 116024

发布日期:2021-08-15
作者简介:原志超，男，博士生，研究领域为高超声速气动热数值模拟。E-mail: zhichaoy@126.com 通讯作者:高效伟，男，博士，教授，研究领域为边界单元法，飞行器热防护系统优化设计。

An Implicit Scheme with Unstructured Gird for Catalytic Surface

School of Aeronautics and Astronautics，State Key Laboratory of Structural Analysis for Industrial Equipment，Dalian University of Technology，Dalian 116024，China,School of Aeronautics and Astronautics，State Key Laboratory of Structural Analysis for Industrial Equipment，Dalian University of Technology，Dalian 116024，China and School of Aeronautics and Astronautics，State Key Laboratory of Structural Analysis for Industrial Equipment，Dalian University of Technology，Dalian 116024，China

Published:2021-08-15

摘要/Abstract

摘要： 为了解决格点型非结构有限体积法隐式处理完全催化壁面时边界上出现非物理解的情况，提出了壁面边界隐式矩阵“松弛”处理方法，使得格点型非结构有限体积法可以用于高超声速非平衡条件下气动热的预测。提出了考虑等温壁面和辐射平衡壁面下，处理完全催化壁面的“松弛”方法，并给出了边界控制体上隐式矩阵的具体形式。针对高超声速化学非平衡流动，使用非结构网格对典型的圆柱绕流以及球头绕流问题进行模拟，并与结构网格计算结果进行比较分析。验证了所提出的强制边界条件的“松弛”格式可以正确模拟完全催化壁面边界条件。在此基础上，针对不同质量分数配比的自由来流进行计算，分析壁面处组元的质量分数以及热流，计算结果显示即使在自由来流中两个组元质量分数各占50%的情况下，被“松弛”的组元质量方程的残差仍然能够到10-8量级，壁面上热流计算正确并且组元质量分数与自由来流保持一致。这表明本文方法在自由来流中质量分数最大组元不占优的情况下仍然有效。

关键词: 非结构网格；高超声速；化学非平衡；壁面催化；隐式处理方法

Abstract: In order to solve the nonphysical phenomenon when the unstructured finite volume solver based on cell-vertex scheme is employed for fully catalytic surface， the ‘relaxation’ method for implicit matrix on boundary is developed. Thus， the unstructured finite volume solver based on cell-vertex scheme can be used to predict aero-heating for hypersonic nonequilibrium. The strong implicit scheme with ‘relaxation’ considering fully catalytic surface under isothermal and radiative equilibrium wall is developed and the implicit matrix for the control volume on boundary is given. This scheme based on unstructured grid is implemented for the simulation of hypersonic chemical nonequilibrium flows around a cylinder and a sphere. Numerical results are obtained and compared with that of structured grid solver. It indicates that the ‘relaxation’ scheme for strong boundary condition can simulate fully catalytic surface accurately. Furthermore， the simulation for different mass fraction compositions of freestream is performed and the mass fraction and heat flux on the wall are analyzed. The results show that the residual order of mass conservation equation for the ‘relaxation’ species can reach 10-8 and the mass fraction on the wall is the same to that of freestream even if the mass fraction of two species both equal to 50%. It shows that the developed method is still effective when the maximum mass fraction species does not dominate.

Key words: Unstructured grid；Hypersonic；Chemical nonequilibrium；Surface catalysis；Implicit scheme

原志超,黄世璋,高效伟. 非结构网格催化壁面的隐式处理方法[J]. 推进技术, 2018, 39(2): 251-260.

YUAN Zhi-chao,HUANG Shi-zhang,GAO Xiao-wei. An Implicit Scheme with Unstructured Gird for Catalytic Surface[J]. Journal of Propulsion Technology, 2018, 39(2): 251-260.

参考文献

[1] 杨勇, 张辉, 郑宏涛. 有翼高超声速再入飞行器气动设计难点问题[J]. 航空学报, 2015, 36(1): 49-57.
[2] 王江峰, 伍贻兆, 季卫栋, 等. 高超声速复杂气动问题数值方法研究进展[J]. 航空学报, 2015, 36(1): 159-175.
[3] Meng Z, Fan H, Peng K, et al. A Hypersonic Aeroheating Calculation Method Based on Inviscid Outer Edge of Boundary Layer Parameters[J]. Acta Astronautica, 2016, 129: 429-437.
[4] 刘健, 原志超, 杨恺, 等. 高超声速飞行器多层复杂热防护结构气-固耦合快速热分析方法[J]. 推进技术, 2016, 37(2): 227-234. (LIU Jian, YUAN Zhi-chao, YANG Kai, et al. Fast Algorithm of Coupled Flow-Thermal for Multi-Layered Complex TPS of Hypersonic Aircraft[J]. Journal of Propulsion Technology, 2016, 37(2): 227-234.)
[5] 杨恺, 高效伟. 高超声速飞行器关键部位气动热计算[J]. 计算力学学报, 2012, 29(1): 13-18.
[6] 潘静, 阎超, 耿云飞, 等. 高超声速飞行器气动防热新概念研究[J]. 力学学报, 2010, 42(3): 383-388.
[7] 罗祖分, 宋保银. 航天器再入大气层热力分析[J]. 航空动力学报, 2016, 31(10): 2507-2514.
[8] Mazaheri A, Gnoffo P A, Johnston C O, et al. LAURA Users Manual: 5.3-48528[R]. NASA TM-2010-216836.
[9] Wright M, Candler G, Bose D. Data-Parallel Line Relaxation Method for the Navier-Stokes Equations[J]. AIAA Journal, 1998, 36(9): 1603-1609.
[10] 董维中, 高铁锁, 丁明松, 等. 高超声速飞行器表面温度分布与气动热耦合数值研究[J]. 航空学报, 2015, 36(1): 311-324.
[11] 李海燕, 唐志共, 杨彦广, 等. 高超声速飞行器高温流场数值模拟面临的问题[J]. 航空学报, 2015, 36(1): 176-191.
[12] 李海燕. 高超声速高温气体流场的数值模拟[D]. 绵阳:中国空气动力研究与发展中心, 2007.
[13] 李桦, 潘沙, 田正雨. 高超声速滑移流气动热并行计算数值模拟[M]. 北京:国防工业出版社, 2013.
[14] 柳军. 热化学非平衡流及其辐射现象的实验和数值计算研究[D]. 长沙:国防科学技术大学, 2004.
[15] Gnoffo P. Simulation of Stagnation Region Heating in Hypersonic Flow on Tetrahedral Grids[R]. AIAA 2007-3960.
[16] Gnoffo P, White J. Computational Aerothermodynamic Simulation Issues on Unstructured Grids[R]. AIAA 2004-2371.
[17] 唐静, 郑鸣, 邓有奇, 等. 网格自适应技术在复杂外形流场模拟中的应用[J]. 计算力学学报, 2015, 32(6): 752-757.
[18] Anderson W K, Bonhaus D L. An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids[J]. Computers & Fluids, 1994, 23(1): 1-21.
[19] 王江峰, 余奇华, 伍贻兆. 高超声速热化学非平衡绕流分布式并行计算[J]. 中国科学技术大学学报, 2008, 38 (5): 529-533.
[20] 王江峰, 伍贻兆, 余奇华. 混合网格化学非平衡绕流通量分裂格式及并行算法[J]. 中国科学技术大学学报, 2006, 36(10): 1034-1039.
[21] Diskin B, Thomas J L. Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Inviscid Fluxes[J]. AIAA Journal, 2011, 49(4): 836-854.
[22] Diskin B, Thomas J L, Nielsen E J, et al. Comparison of Node-Centered and Cell-Centered Unstructured Finite-Volume Discretizations: Viscous Fluxes[J]. AIAA Journal, 2010, 48(7): 1326-1338.
[23] Blazek J. Computational Fluid Dynamics: Principles and Applications(3nd)[M]. Oxford: Butterworth-Heinemann, 2015.
[24] 欧阳水吾, 谢中强. 高温非平衡空气绕流[M]. 北京: 国防工业出版社, 2001.
[25] 王羽, 屈崑, 蔡晋生. 高超声速化学反应源项雅克比矩阵的对角化[J]. 航空学报, 2016, 37(5):1419-1427.
[26] 王羽, 蔡晋生, 屈崑, 等. 增强高超声速化学反应流数值计算的稳定性方法[J]. 宇航学报, 2016, 37(9): 1135-1141.
[27] Gnoffo P A, Gupta R N, Shinn J L. Conservation Equations and Physical Models for Hypersonic Air Flows in Thermal and Chemical Nonequilibrium[R]. NASA-TP-2867.
[28] Harten A. High Resolution Schemes for Hyperbolic Conservation Laws[J]. Journal of Computational Physics, 1983, 49(3): 357-393.
[29] 王兰. 超燃冲压发动机整机非结构网格并行数值模拟研究[D]. 绵阳:中国空气动力研究与发展中心, 2007.
[30] Yuan Z, Huang S, Gao X, et al. Effects of Surface-Catalysis Efficiency on Aeroheating Characteristics in Hypersonic Flow[R]. ASCE/AS.1943-5525, 2016.
[31] 杨恺, 原志超, 朱强华, 等. 高超声速热化学非平衡钝体绕流数值模拟[J]. 推进技术, 2014, 35(12): 1585-1591. (YANG Kai, YUAN Zhi-chao, ZHU Qiang-hua, et al. Numerial Simulation of Hypersonic Thermochemical Nonequilibrium Blunt Body Flows[J]. Journal of Propulsion Technology, 2014, 35(12): 1585-1591.)
[32] 程晓丽, 艾邦成, 王强. 基于分子平均自由程的热流计算壁面网格准则[J]. 力学学报, 2010, 42(6): 1083-1089.
[33] 黎作武. 近似黎曼解对高超声速气动热计算的影响研究[J]. 力学学报, 2008, 40(1): 19-25.(编辑:梅瑛)　　　　　　　　　　　　　*　收稿日期:2016-12-14；修订日期:2017-02-23。作者简介:原志超,男,博士生,研究领域为高超声速气动热数值模拟。E-mail: zhichaoy@126.com通讯作者:高效伟,男,博士,教授,研究领域为边界单元法,飞行器热防护系统优化设计。 E-mail: xwgao@dlut.edu.cn