固体火箭发动机内流场分区耦合数值计算

采用分区耦合方法计算固体火箭轴对称燃烧室与喷管流场。对于低速的燃烧室流场，选用不可压流的Ｎ－Ｓ方程描述并用ＳＩＭＰＬＥＣ方法数值求解；对于高Ｒｅ数的喷管流场，则采用Ｅｕ－ｌｅｒ方程描述并用ＳＣＭ方法求解。计算时用燃烧室出流为喷管流场提供入口参数，同时用喷管流场压强分布反馈影响燃烧室流动状况。对耦合边界条件处理方法进行了探讨。对典型的侧壁加质燃烧室与喷管流场进行了计算，计算结果揭示了单独喷管流场计算难以反映的喷管收敛段近壁区的低速区域，与已有的燃烧室流场实验结果一致并反映了燃烧室与喷管流场之间的联系，较好地模拟了流动中的物理现象。     

轴对称超临界粘性喷流的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式中的矢通量分裂技术，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了收敛喷管的内流场和喷射流场，通过对趋临界工况下喷管跨音速流场的数值计算，表明本文的结果与试验结果是吻合的。     

固体火箭喷管两相粘性跨音速流场计算

本文进行了固体火箭喷管两相粘性湍流跨音速流场计算.粘性的气相控制方程用隐式近似因子分解法求解,粒子方程采用跟踪粒子轨迹的特征线法求解,粘性湍流选用代数模型,气相和凝相充分地偶合.CFL数可以取500左右,收敛速度快,使粘性两相跨音速喷管流场计算耗费的机时达到工程计算可以接受的程度.通过计算,获得了在粘性和粒子同时作用下的流场参数分布.这对固体火箭发动机流场需要同时考虑粒子和粘性作用的专题研究大有裨益.     

粘性跨音速喷管流场计算

本文继用隐式近似因子分解法成功地计算无粘跨音速喷管流场之后,用同样的方法,结合任意曲线坐标系,通过非定常方程在相当长时间后的定常解,针对矩形截面喷管和轴对称喷管粘性跨音速流场,求解了可压缩层流薄层N-S方程,获得了核心流,特别是边界层中的流动参数.对于矩形截面和轴对称两种喷管进行了计算,其结果和实验数据相当一致.将本文的方法应用到粘性两相流动计算,可望较多的节省机时.     

轴对称喷管流场分析的有限体积法

本应用有限体积法对轴对称喷管的跨音速流场进行了分析，采用显式预校两步MacCormack格式进行计算。运用隐式残差平均技术及当地时间步长来加速计算的收敛，提高格式的稳定性。对三个不同几何参数的轴对称喷管的数值分析表明，计算结果与试验数据吻合良好。     

超声速侧向多喷流干扰特性数值模拟

为了研究多个喷流喷管对导弹控制力的干扰影响,本文通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格。通过对多种喷管控制组合的超声速横向喷流干扰流场的数值模拟,研究和分析了喷口附近流场的涡系结构和波系结构,并将喷管几种排列组合对导弹喷流干扰力放大因子的影响进行了分析研究,得出一些多喷流干扰的结论。     

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

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

并联TBCC排气系统流场结构数值模拟及实验研究

为探索并联TBCC排气系统在宽马赫数飞行范围内的性能变化规律,针对某TBCC排气系统构型开展了飞行包线内的数值模拟分析。为验证数值模拟方法的可靠性,对该实验模型完成了典型工作点静特性风洞实验,获得了相应的流场结构和壁面压力分布规律,并将实验结果与数值模拟结果进行了对比。研究结果表明:并联TBCC排气系统双通道气流在涡轮下壁面出口处发生严重干扰导致出口流场复杂。优化设计TBCC排气喷管,有利于降低气流干扰对整个排气系统流场结构和推力性能的影响;飞行包线内,排气系统总的推力系数均高于0.9。气流大面积分离导致冲压发动机低马赫数时性能较低,但其对排气系统整体的推力性能影响很小;数值模拟得到的壁面沿程压力分布、流场波系结构等均与实验结果非常吻合,因此,数值计算结果用来评价TBCC排气系统性能的可靠性得到实验验证。     

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

密切内锥乘波前体进气道一体化设计和性能分析

采用特征线方法设计了具有直线初始激波、内收缩段消除激波反射、出口参数均匀可控的基准内锥流场。基于密切内锥(Osculating Inward turning Cone,OIC)乘波体设计方法,发展了密切内锥乘波前体进气道(Os-culating Inward turning Cone Waverider Inlet,OICWI)一体化设计技术。基于基准内锥流场和前体进气道一体化设计技术,设计了密切内锥乘波前体进气道。采用数值方法对设计的密切内锥乘波前体进气道进行了计算分析,结果表明无粘流场结构和基准内锥流场吻合,无粘模拟结果和理论设计结果吻合。粘性数值模拟结果显示一体化进气道具有较高的流量捕获率及总压恢复特性,进气道出口流场分布均匀。     

Godunov方法在复杂流场数值模拟中的应用

本文根据Ｖａｎ Ｌｅｅｒ提出的方法给出了求解三维非定常Ｅｕｌｅｒ方程和三维定常Ｅｕｌｅｒ方程的二阶Ｇｏｄｕｎｏｖ方法。方法要点为：在各体积离散单元中,用流动变量的分段线性分布代替一阶方法中的分段常量分布,在确定线性分布的斜率中引入单调性限制条件;在离散单元界面构造Ｒｉｅｍａｎｎ问题,由Ｒｉｅｍａｎｎ间断解确定单元侧面的通量;用预测－修正两步法对离散方程推进求解,用上述方法结合单区网格或多区重迭网格     

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

型面旋转变马赫数风洞喷管的优化设计

型面旋转的变马赫数可调风洞喷管因其能够实现出口马赫数的连续快速变化、调节容易和流场品质较好等优点,成为目前地面变马赫数风洞研究的热点。将出口马赫数范围为2.0~4.0的变马赫数风洞喷管作为研究对象,基于iSIGHT和Fluent等软件搭建了适于该风洞喷管的优化设计平台,研究了型线设计中设计点马赫数、初始膨胀线和附面层修正量等因素对变马赫数出口流场的影响,应用NSGAⅡ算法在依据试验设计点所构建的Kriging近似模型上搜索,得到最优的风洞喷管基准型线,并对该基准型线进行二维、三维数值仿真校核,结果显示：通过旋转优选得到的基准型线能够得到马赫数范围为2.0~4.0的变马赫数出口流场,且流场均匀性基本满足固定型面、单设计点风洞喷管的国军标（GJB-1179-91）合格指标,为变马赫数可调风洞喷管在变马赫数风洞中的实际应用奠定了良好的基础。     

超声速/高超声速双拐点喷管设计

为实现直连式试验台、高温风洞等试验设备的多马赫数运行,提出了双拐点喷管设计方法.喷管分2段设计,第1段共用,采用3次B-Spline函数描述喷管轴线马赫数分布.首先采用特征线方法求解Euler方程,得到无黏的理想喷管型面.其次采用参考温度方法求解边界层位移厚度,对无黏壁面进行修正得到实际壁面.共用段喷管出口的平行均匀流作为第2段喷管设计的初值.为验证设计方法的可行性,设计了中间马赫数为3.0,出口马赫数分别为4.0,4.5和5.0的双拐点喷管,并采用雷诺平均的Navier-Stokes方程对设计的喷管流场进行数值模拟.计算结果表明:喷管出口流场均匀,试验菱形区的马赫数误差小于1.2%.该方法提高了喷管设计精度,保证消波干净,为直连式试验台、高温风洞等设备的多个喷管共用一套动力系统提供了基础.      

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

双弯曲入射激波的可控中心体内收缩基准流场设计

采用有旋特征线法设计了一种双弯曲入射激波的可控中心体内收缩基准流场,两道入射激波交于中心体起始点,入射激波和反射激波通过给定激波径向总压恢复系数分布进行反设计,壁面通过给定轴向马赫数分布规律进行反设计.该基准流场分为“三波四区”且压缩效率较高.基于该基准流场设计了圆形进口内收缩进气道并进行了黏性修正,数值计算结果表明:内收缩进气道设计点核心区的流场特征和激波形状与基准流场基本一致;在来流马赫数为4.0~7.0时进气道具有较高的压缩效率和良好的流量捕获能力,设计点喉道截面增压比和总压恢复系数分别为17.7和0.729;来流马赫数为5.0~7.0时内部总阻力系数变化平缓,从0.23下降为0.22.      

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

高超声速内收缩进气道轴对称基准流场改进

针对轴对称基准流场中前缘曲激波靠近中心体的部分激波强度过大现象,基于马赫数分布可控反设计方法,将这道前缘曲激波分解为一道较弱弯曲激波和部分等熵压缩波,改进的基准流场存在“四波四区”结构且压缩效率明显提高.基于该改进的基准流场和常规“两波三区”基准流场分别设计了圆形进口的内收缩进气道并对其流场特点和性能进行数值研究.结果表明:改进的进气道的流场能较好保持基准流场的特点;在来流马赫数为4.0~7.0范围内具有较高的压缩效率和良好的流量捕获能力,设计点的出口压比和总压恢复系数分别为17.56和0.540;改进的进气道采用来流马赫数从高到低前缘弯曲激波和汇集的等熵压缩波依次封口的设计概念,在提高流量捕获能力的同时减小了总压损失,总体性能优于常规进气道,来流马赫数为7.0时总压恢复系数相对提高了23.6%,来流马赫数为4.0时流量系数相对提高了5.7%.      

1

Ｄｉｒｅｃｔｏｐｅｒａｔｉｎｇｃｏｓｔｓｆｏｒｔｒａｎｓｐｏｒｔａｉｒｃｒａｆｔｈａｖｅｉｎｃｒｅａｓｅｄｒａｐｉｄｌｙｓｉｎｃｅｔｈｅｏｉｌｃｒｉｓｉｓｏｆｔｈｅ１９７０ｓ．Ｔｈｅｐｏｔｅｎｔｉａｌｆｏｒｌａｒｇｅｆｕｅｌｓａｖｉｎｇｓａｎｄｔｈｅ...     

一种适合于旋翼前飞非定常流场计算的新型运动嵌套网格方法

针对旋翼流场计算中关键的运动嵌套网格技术,提出了一个新的网格单元关系快速判断方法。该方法通过计算特殊网格点在背景网格上的位置来建立桨叶网格在背景网格中所占据范围的关系,并用于嵌套网格洞单元的识别,又进一步将该方法推广用于Inverse Map的建立,在此基础上,建立了一套用于前飞流场计算的运动嵌套方法。通过多种前飞旋翼测试算例,分别评估了背景网格、桨叶网格尺寸改变对计算时间的影响,同时验证了该方法的可靠性,结果表明构建网格嵌套关系所用时间随桨叶网格和背景网格尺寸的变化幅度较小,且具有良好的鲁棒性,能够满足旋翼前飞状态计算的需要。最后,在基于前飞非定常N-S/Euler流场求解器中,采用该方法模拟了UH-60A和7A旋翼悬停/前飞状态的气动特性,通过计算时间的对比体现出该方法的高效性,并且流场计算结果与试验数据对比表明该方法可以有效地用于旋翼前飞非定常流场的计算。