跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

喷管一维流场化学非平衡通用计算软件研制

研制的火箭发动机一维喷管非平衡流通用计算软件，建立在一个统一的化学反应库之上。该通用计算软件可对推进剂元素为碳、氢、氧、氮、氟、氯的火箭发动机喷管进行一维化学非平衡流数值模拟。通过计算推进剂元素为氢、氧和碳、氢、氧、氮的喷管一维化学非平衡流，得到了喷管流场的速度、温度、密度、压力及各组分的质量百分比等参数，为工程设计提供了依据。     

喷管碳基喉衬的烧蚀

一、计算方法和结果 喷管材料的烧蚀主要是热和化学作用的结果。热量来自碳与燃气中的氧化元素之间的异相化学反应,以及对流和幅射传热。本文中的解析模型用来预测喷管温度分布及烧蚀率。在处理联结核心流及喷管壁面的边界层流动时,把传热传质数值作为壁内导热计算的输入数据。边界层的边界条件是由一维无粘定常流的解来确定的。在材料响应的分析中介入了幅射边界条件,在UTC数据上引入了一个有效辐射系数。用有限元计算方法求解瞬态热传导方程来确定喷管壁内温度分布随时间的变化关系。     

基于CFD数值模拟的拉瓦尔喷管流场分析

结合一维等熵管流理论,采用CFD方法对某型拉瓦尔喷管典型工作的流场进行整场数值模拟,阐述了网格生成,流场初始化,以及RNG湍流模型使用方法,并对拉瓦尔喷管工作特性进行了气动分析。研究表明,CFD计算结合一维管流理论能明显提高计算效率,数值计算结果与理论相符较好。     

化学非平衡喷管流计算方法的研究

本文提供了考虑化学反应非平衡和振动非平衡的喷管流动计算方法。在忽略扩散、粘性和热传导及辐射损失的假设下，得到完全气体一维定常流动方程，再加上振动能松驰方程和化学反应速率方程即可数值积分求解。采用扰动法得到接近平衡流的一组积分初值。采用“逆积分方法”求解通过喉道。计算结果是令人满意的。     

轴对称拉瓦尔喷管流场分析

采用一维管流理论分析了轴对称拉瓦尔喷管典型工作状态,结合理论分析对喷管流场进行CFD仿真计算。经研究表明：（1）结合理论分析的流场初始化、以及FMG和外推计算方法,CFD计算效率明显提高;（2）不同拉瓦尔喷管工作状态,需采取不同的壁面函数法来处理近壁区;（3）一维管流理论分析中,将管内强激波简化为正激波的常规做法,计算误差较大,不尽合理。     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

一维非平衡喷管流动计算

引言 在电弧风洞的运行中,喷管前室的气流温度很高,空气各组元之间会发生化学反应,双原子组元还要出现振动激发。混合气体通过喷管加速,化学反应和双原子分子的振动都达不到平衡,结果使得喷管各处流场参数与平衡流的计算数据有些差异,这就是喷管流动的非平衡效应。     

C/C复合材料喷管烧蚀壁面退移流固耦合仿真研究

为了保障固体火箭发动机C/C喷管的可靠性,建立了一套正确反映发动机喷管烧蚀过程的流固耦合计算模型,以实现对喷管烧蚀率的高精度预估。依据热化学烧蚀理论以及喷管内燃气与喷管结构体界面的质量平衡和能量平衡关系,建立并验证了考虑壁面退移的C/C喷管流固耦合方法,实现了燃气流动、异相化学反应、结构体传热三者间的耦合。通过实验发动机喷管的烧蚀计算,论证了模型的正确性,并分析了不同金属铝含量对烧蚀率的影响,计算所得的烧蚀率与实验值最大相对误差为4.3%,与不考虑壁面退移的耦合算法计算结果对比,计算精度最高可提升46%。计算结果表明：C/C喷管在喉部附近烧蚀最为严重;推进剂中Al含量的增加导致燃气中氧化组分浓度降低,进而减少了烧蚀速率,这些结论与C/C喷管烧蚀相关研究结果一致。     

固体火箭发动机喷管结构缝隙设计

通过热和结构计算，对固体火箭发动机喷管缝隙设计进行了分析。将燃气简化为一维等熵流以确定喷管内型面所承受的温度和压力载荷．基于三维有限元模型．计算了喷管的瞬时温度场。导入温度场分析结果，并用点一点接触单元模拟不同材料之间的接触状态，计算了温度和压力载荷联合作用下喷管的位移场和应力场。综合多个时刻的计算结果得到发动机工作过程中喷管结构缝隙和接触应力的变化趋势，并据此对设计缝隙提出了修改意见，结果可为发动机喷管设计提供参考。     

类升力体外形俯仰阻尼特性数值研究

采用有限差分方法求解薄层近似的非定常Navier-Stokes方程，定常流场采用交替方向隐式分解的NND格式，非定常流场采用四步Runge-Kutta方法，在保证时、空二阶精度前提下引入变系数残值光顺技术提高非定常流场计算效率，复杂带翼外形的空间网格通过求解抛物化的椭圆型方程生成，最后在Etkin理论下给出球锥及类升力体外形的俯仰阻尼导数计算结果。     

飞机亚跨超绕流的数值模拟研究

本文以Euler方程为数学模型,采用一种高精度的TVD(Total Variation Dimishing)离散格式及一种含近似因式分解的推进迭代方法,求解亚跨超绕流’流场。通过若干算例的试算,证明方法是可行的,可以用来模拟飞机的复杂流场。     

多段翼构型的积冰数值模拟

基于非结构网格提出一种针对多段翼构型的积冰数值模拟方法.通过求解气体欧拉方程得到多段翼的绕流流场,基于空气绕流流场的气流速度分布,采用迭代法对水滴控制方程进行数值求解,得到多段翼各翼面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对某典型三段翼霜冰和明冰的结冰情况进行了预测研究,根据多段翼的流场特点提出一种明冰冻结模型的求解思路,并分析了积冰对多段翼性能的影响.      

低超声速来流中机翼跨声速非定常气动力数值解法

本文给出一种低超声速来流中机翼跨声速非定常流的计算方法。该方法使用跨声速小扰动理论,正确地计入了沿流向平面流动特征线的压力波的运动边界条件。流动方程是使用时间精确交替方向隐式(ADI)有限差分解法求解。试算了M=1.1的矩形翼和M=1.1、1.34的F5机翼的定常和非定常气动力,并与国外的风洞试验和计算结果作了比较,符合很好。表明本方法是有效的。     

旋翼三维非定常跨音速绕流的数值模拟

 利用小扰动位势方程，用有限差分方法求解了旋翼跨音速流场。导出了严格的三维非定常位势方程远场无反射边界条件，并成功地用于数值计算中，缩小了物理求解域的空间尺度。采用了贴体坐标和网格加密技术，以得到合理的非均匀网格分布。应用较少的网格点数成功地计算了无升力情况下旋翼的三维非定常亚、跨音速绕流流场，得到了与实验结果符合良好的数值解     

平面大头叶栅跨声速流动的一种有效解法

准确数值预估叶片头部表面的压力分布,一直是叶片气膜冷却设计中十分重要而又尚未圆满解决的问题。本文针对这个问题提出了一个高准确度、高效率的计算方法。 在文献[1]的基础上,本文从任意正交曲线坐标系中的基本方程出发,引入Von Misses变换,导出了相应的流线控制方程,并具体地在拟边界层坐标系中计算了RA叶栅流场,结果与实验符合甚好,据作者所知,本文方法是已见到的计算同类叶栅方法中准确度最高、花费机时最少的。     

三维湍流室内气流数值解

本文基于有限差分法,k-ε模型和p-v法提出了计算三维湍流室内气流问题的数学模型,且用它求解了某水电厂厂房的流场。同时还做了相应的模型试验,并把二者进行了比较,获得了满意的结果。     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

环形回流燃烧室两相反应流场的数值研究

本文在任意曲线坐标系下,采用区域数据即时传递法,对包括冷却通道在内的整个环形回流燃烧室的三维两相反应流进行了数值模拟。气相采用Euler方法处理,并采用标准k-ε双方程紊流模型,EBU-Arrhenius紊流燃烧模型,六通量热辐射模型,液相采用Lagrange法处理。在非交错网格体系下,气相用SIMPLE法求解,液相采用颗粒群轨道模型,并用PSIC算法对其进行数值求解。      

二维涡轮级的定常/非定常数值模拟

为了解定常与非定常计算差别对工程应用参考价值的影响，本文通过求争二维Navier-Stokes方程，针对一级涡轮叶栅分别进行了定常／非常数值模拟。结果发现，对于叶片表面压力分布，第一排导叶定常计算与非定常计算得到一致的结果，但是在动叶上差别较大，尤其在叶排小间距时。当叶排间距拉大后，定常计算结果趋近于非定常计算结果时间平均值。