预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

非平衡流计算方法及其模拟激波诱导振荡燃烧

给出一种数值模拟高超声速化学非平衡流动的计算方法，采用流动项与化学反应生成源项解耦处理，即可以有效地解决非平衡流动计算中遇到的刚性问题，还具有公式推导简单和应用方便的优点。从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超音速H2／Air预混合气体钝体绕流场，较好分辨出流场中激波诱导振荡燃烧非定常流动现象，表明计算方法具有较高的精度。对计算方法的效率和网格数的影响也进行研究。     

热与化学耦合非平衡对热流和电子密度的影响

采用单、双温度模型两种粘性激波层方程组对弹头驻点区非平衡流场作了数值模拟。计算结果表明，在算例条件下热非平衡与化学非平衡的耦合效应对壁面热流和电子密度有明显影响。文中给出了完全催化壁和非催化壁两种情况下以及采用不同振动松弛时间公式时上述影响的变化规律。     

磁场对高超声速弱电离气体流动的影响

对偶极子磁场作用下的三维钝头体高超声速黏性绕流的化学非平衡流动进行了数值模拟。电导率利用组分公式计算,化学模型为7组元、6反应模型。应用诱导磁场方法将磁场与流场耦合,控制方程的空间离散采用有限体积法,扩散项用中心差分格式,对流项采用二阶迎风格式,时间推进采用3步Runge-Kutta法。计算结果表明,外加偶极子磁场使激波脱体距离增加,壁面摩擦系数和表面热流密度增大。与无磁场作用时相比,在0.153T外加磁场作用下,冻结流中的激波脱体距离增加约3倍,局部壁面摩擦系数最大增加63%,局部表面热流密度最大增加61%;而化学非平衡流中的激波脱体距离增加约0.5倍,局部壁面摩擦系数最大增加47%,局部表面热流密度最大增加31%,且化学非平衡流中激波层内温度的最大值约为冻结流中的64%。     

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

烧蚀对再入体绕流电子数密度影响的数值研究

从化学非平衡NS方程或其简化形式(粘性激波层方程、PNS方程)出发,采用分区求解的方法数值模拟含泰氟隆烧蚀产物的高超声速再入体流场.首先计算分析了RAM-C飞行试验状态和一个考核算例,对计算方法和程序作了初步验证.其次耦合烧蚀壁面边界条件和流体控制方程,采用19组分28个反应的化学反应模型,对有、无烧蚀情况下高超声速钝锥体再入绕流及底部流场进行了数值模拟,深入分析了烧蚀产物对再入流场电子数密度及温度分布的影响.     

化学动力学模型对H2/AIR超燃模拟的影响

为了研究超燃流动中不同化学反应模型的影响，提出流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似N－S方程出发，采用ENO差分格式，数值模拟了超声速H2／Air 预混合气体中钝体激波诱导的振荡燃烧流场，针对飞行速度大于和接近爆轰速度两种不同流动状态，比较Evans-Schexnayder提出的7组分8反应机理（ES模型），Jachimowski提出的13组分32反应机理（J模型）及其简化为9组分19反应机理（J－1模型）对流场参数的影响，计算表明反应机理引起压力场变化较小，对温度等参数影响较大，对于飞行马赫数低于7的问题，取不包括N原子反应的Jachimowski的9组分19反应机理较合适。     

高速高温流场电子能非平衡的数值模拟

采用描述电子能非平衡的三温度模型,结合11组分空气的化学反应模型,对多种高速高温热化学非平衡流场开展数值模拟,并与描述电子能平衡的两温度模型结果进行对比,研究电子能非平衡对高超声速流场特性的影响。圆球弹道靶试验算例表明电子能非平衡不影响激波脱体距离。RAM-C II飞行器的4个飞行工况算例表明,尽管两温度和三温度模型结果存在差异,但二者电子数密度分布的趋势和量级接近,均可与飞行试验数据保持一致,其中三温度模型的预测效果更好。FIRE II飞行器极高温流场模拟结果显示,电子能非平衡几乎不影响飞行器表面的对流传热。     

甲烷/空气超声速燃烧流动数值模拟

用流动项与化学反应生成源项解耦处理的化学非平衡流动计算方法，从薄层近似三维N-S方程出发，采用ENO差分格式数值模拟了超声速冲压加速器简化模型中，高速甲烷气流从后体多个喷口射入高超声速空气流形成流场，研究了高温异质气体效应和19组分65反应模型的非平衡效应对冲压加速器表面压力分布的影响。计算表明新型解耦方法适合反应机理复杂的碳氢燃料超燃冲压发动机内部流动模拟，为开发应用软件系统打下基础。     

带座舱飞船高超声速再入非平衡流场的数值研究

本文针对带座舱飞船高超声速再入大气层过程中存在的严重气动加热现象 ,利用简单隐式TVD差分格式和激波捕捉法 ,数值求解三维化学非平衡Navier Stokes方程 ,其中化学模型是 7组分 1 5个化学反应的空气化学模型 ,对带座舱飞船再入高度为 70km和速度为6 7km/s的化学非平衡流场进行了数值模拟 ,给出了带座舱飞船再入迎角为 0°,1 0°和 2 0°情况下的壁面热流、表面压力和电子数密度等参数分布     

11组元轴对称热化学非平衡流场的数值研究

本文从轴对称热化学非平衡Ｎ－Ｓ方程出发,采用双温度、十一组元反应气体模型,利用隐式ＮＮＤ有限差分格式的时间预处理技术数值求解了ＦＩＲＥＩＩ飞船热化学非平衡流场。文中同时给出了七组元模型的计算结果,并进行了结果比较。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

超音速钝体无粘绕流的计算

研究用三维特征线相容性关系的迎风差分格式,数值模拟超音速和高超音速钝体无粘绕流流场。用钝体反方法计算头激波及其波后的流场,用轴对称流特征线法确定初值面上的流动参数。算例的计算结果与实验数据比较吻合     

非结构网格下含冷却孔的涡轮转子三维流场计算

在转动坐标系里,讨论非结构网格下,含多排冷却孔的涡轮转子三维流场雷诺平均Navier-Stokes方程组的求解。计算表明:含冷却孔的转子流场计算要比静子流场复杂的多。为了得到较为满意的流场数值解,本文对叶片压力面与吸力面采用了不同的湍流模式、并且还对流场迭代中初始场选取方式以及转子叶尖间隙计算站的设置等问题进行了一系列的改进。几个典型流场算例表明:初始场选取、湍流模式的合理选用对转动坐标系下获得一个收敛的三维流场数值解至关重要;另外,转子叶尖间隙站的设置直接影响着间隙附近的流场,尤其对现代高负荷涡轮转子,间隙流动的研究就更为重要。      

高超声速喷管非平衡尺度效应的数值分析

采用轴对称热化学非平衡全N—S方程，数值分析了高超声速喷管流场中非平衡引起的尺度效应。并给出计算条件下喷管轴线上的物理量和化学组元分布。计算结果表明，在喉部下游不远处流动趋近于化学冻结和热力冻结状态。流场(特别是核心无粘流区)的特性不仅取决于流动尺度的相对值，也取决于其绝对值。     

高超声速喷管非平衡黏性流动的数值研究

分别采用平衡气体模型、化学非平衡气体模型和热化学非平衡气体模型,通过求解轴对称Navier-Stokes方程组,数值模拟了法国Marseille高焓激波风洞锥型喷管(H₀=10.3MJ/kg),分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、振动温度、马赫数和组分质量数在轴对称线上的分布与喷管出口附近的速度和温度沿径向分布等结果。计算结果表明:化学反应速度加快,会导致喷管出口温度增加,振动能的冻结会导致喷管出口温度降低。     

基于Busemann压升规律的可控消波内转基准流场设计

通过将经典Busemann设计方法和特征线反设计方法相结合，实现了对基准Busemann流场的气动截短，构建了一种具有基准Busemann流场截短压升规律的可控消波内转基准流场.通过数值模拟对可控消波内转基准流场及其追踪得到的“糖勺”型进气道进行了无黏验证分析.结果表明:特征线和CFD计算结果相吻合，可控消波内转基准流场设计合理可行.该基准流场继承了Busemann设计方法的高效压缩特性，且反射激波得到有效控制，基本实现消波，性能优于传统的截短Busemann流场.在设计点马赫数为7条件下，喉部截面参数均匀，增压比为18.32，总压恢复系数为0.878，压缩效率为0.936，隔离段内几乎无损失,出口气流匀直，气流角均在±0.4°以内.流线追踪得到的“糖勺”型进气道出口形状更加饱满，流动特征与可控消波内转基准流场基本一致.      

Nonequilibrium recombination after a curved shock wave

The effect of nonequilibrium recombination after a curved two-dimensional shock wave in a hypervelocity dissociating flow of an inviscid Lighthill–Freeman gas is considered. An analytical solution is obtained with the effective shock values derived by Hornung (1976) [5] and the assumption that the flow is ‘quasi-frozen’ after a thin dissociating layer near the shock. The solution gives the expression of dissociation fraction as a function of temperature on a streamline. A rule of thumb can then be provided to check the validity of binary scaling for experimental conditions and a tool to determine the limiting streamline that delineates the validity zone of binary scaling. The effects on the nonequilibrium chemical reaction of the large difference in free stream temperature between free-piston shock tunnel and equivalent flight conditions are discussed. Numerical examples are presented and the results are compared with solutions obtained with two-dimensional Euler equations using the code of Candler (1988) [10].