高温化学非平衡计算中源项处理方法.pdfVIP

2004年 7月 北 京 航 空 航 天 大 学 学 报 July 2004 第30卷 第7期 JournalofBeijingUniversityofAeronauticsandAs~onautics Vo1．30 No．7 高温化学非平衡计算 中源项处理方法 马 辉 王发民 (中国科学院力学研究所 ，北京 100080) 摘 要 ：高超声速三维化学非平衡绕流流场的数值模拟中，在计算初期容 易发生组元密度 出现负值 的非物理现象，另外源项的刚性是影响计算稳定性和收敛 速度 的主要原因。根据源项线性化理论 ，采用 2种源项线性化方法处理化学反应源 项．一种为根据流场 内化学反应物理规律构造的线性化方法，该方法能抑制在计算过 程 中组元密度 出现负值现象，提高 了计算稳定性，加速 了收敛速度 ．同时采用 了时间 预处理矩 阵的线性化方法，较好地解决 了非平衡化学反应与流场耦合 的刚性 问题 。 证 明了所构造源项线性化方法相容性 ．数值实验表 明该方法有效地避免了源项计算 中密度 出负的问题 ，加快了计算的收敛速度 ，从而提高了计算效率 ． 关 键 词 ：化学非平衡 ；三维计算 ；高超音速再入 ；源项线性化方法 中图分类号 ：V211．24 文献标识码 ：A 文 章 编 号 ：1001．5965(2004)07—0679—03 Calculationmethodofsourcetermsinhightemperature chemicaInonequilibrium flows MaHui WangFamin (InstituteofMechanics，ChineseAcademyofSciences，Beijing100080，China) Abstract：In thenumericalsimulationofhypersonicchemicalnonequilibrium flows，speciesdensitymaybe— comenegativeatthebeginningofcalculation，chemicalstiffnessisthemainreasonthataffectsthecalculation’Ssta— bilityandconvergence．Based on thebasictheory ofsource term ’Slinearazition，twolinearazition methodswere chosen．Oneisbasedonthephysical rulesofchemicalreactionflows．Themethod’Sadvantageisthatcanimprove thecalculation’Sstabilityandcanavoidthespecies’densitytobecomenegativeduringthe calculation． Th eother methodisthetimepreconditionermethod，which hasthefunction ofdiminishingthestiffnessofthesourcetemr ． Th emethodwasprovedtobewellapproximatesthepaaial differential equations．Th ecalculation resultshowsthat thismethodcaneffect