定常轴对称超音速喷管和超音速射流的计算

本文叙述轴对称超音速喷管和超音速射流的数值计算。内点用MacCormack有限差分格式,分步加粘性;边界用从特征关系列出的二阶精度格式;喷管口用Prandtl-Meyer流公式;射流边界用一个迭代过程确定。对有激波喷管和无激波喷管做了大量的数值实验,并作了比较和讨论,给出了三个喷管的计算结果。     

计算矢量喷管流场的有限元法

本文通过采用扫描有限元法及多步龙格－库塔时间推进格式求解非定常平均Ｎ－Ｓ方程来模拟二元矢量喷管的内部跨速流场,成功地计算了二维有激波的跨声速喷管内部流场及仰 仰角３０°的三维矢量喷管内部流场。     

气动塞式喷管的流场数值模拟和设计参数探讨

针对气动塞式喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象 ,本文采用扩展的压强校正法求解用紊流模型封闭的可压缩湍流平均 N-S方程组 ,利用交错网格系统抹平迭代过程中的数值振荡 ,并通过高分辨率的 TVD格式进一步改善压强校正法的激波捕获效果 ,对气动塞式喷管的流场进行了详细的数值模拟。并进一步探讨了喷管倾角、塞锥长度和二次流流量比对喷管流场结构和性能的影响。     

N2O/HTPB固液火箭发动机喷管两相流计算

利用二维轴对称N-S方程对选用氧化亚氮/丁羟基燃料推进剂的固液混合火箭发动机的喷管两相流进行了计算.计算采用MacCormack时间推进预报校正二步格式,采用了Baldwin-Lomax代数湍流模型和两相平衡流模型.计算了三种氧燃比下4个不同喷管的喷管流场参数,并计算了喷管性能,通过比较两相流和气相流的计算结果,分析了不同氧燃比和喷管形状对喷管性能的影响,认为固液火箭发动机的性能主要受氧燃比的影响,为固液混合火箭发动机的设计提供了依据.      

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

高精度ENO格式在喷管流动模拟中的应用

本文分析了ＥＮＯ格式的特点,并应用于Ｅｕｌｅｒ方程和全ＮＳ方程的迁移项和压力项,模拟了各种喷管流动。首先计算了ＪＰＬ喷管流动,取得了和实验一致的结果,分析了喷管内产生弱激波的原因;其次计算了非定常二维喷管流动,给出了相应的计算结果,并分析了壁面附近流场的拓扑结构;最后计算了摆动喷管的流动及推力矢量计算,取得了推力随摆角变化的规律,和实验结果表明基本一致。     

三维可压缩流欧拉方程组的隐式时间推进法

本文提出了求解三维定常无粘可压缩流场问题的一种隐式时间推进法。在计算空间以逆变速度分量为未知变量的欧拉方程组为控制方程,从而简化壁面边界的处理。在Beam-Warming的AF-方法基础上,用对角化及矢通量分裂等方法提高了计算速度及解的精度和稳定性。作为算例,对伴有激波的方形收放喷管跨声速流动进行了计算,以显示方法的有效性。     

塞锥型面简化与截短对塞式喷管性能的影响

用NND差分格式求解N－S方程，对截短直锥塞式喷管和用特征线法设计的截短型面锥塞式喷管的流场和性能进行了比较。成功地捕获了流场激波、膨胀波及反射压缩波，计算结果与实验数据吻合较好，对不同型面塞式喷管在不同高度下的推力性能和塞锥截短对发动机性能的影响。研究表明，相同长度的截短直锥发动机总推力比型面锥塞式喷管发动机低1％－1.5%；相对30％截短型面塞锥发动机，50％截短型面塞锥发动机的总推力可提高约2.5%。     

环流喷管的轴对称湍流流场计算

对火箭发动机环流喷管进行了二维数值仿真,用数值方法求解Navier－Stokes方程,数值格式为显格式时间推进格式,湍流模型选择可靠性较好的两层代数湍流模型。由于所研究的流场边界形状复杂,计算中采用了分区的方法。通过多个算例的对比,确定了使环形喷管推力效率较高的几何形状。画出了环流喷管内部的计算网格图、速度矢量图、等马赫线图、等压线图及喷管上壁面及轴线的压力分布图。     

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

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

固体火箭发动机有约束流场推力特性计算

为计算带约束流场的固体火箭发动机推力特性，建立了一维准定常变截面的内弹道计算和约束流场流动的数学模型。计算结果表明，约束流场使推力增加14%，冲波进入喷管扩张段。与试验结果很接近，表明计算模型正确，可用于其它型号带约束流场的推力特性计算。     

激波诱导矢量喷管流场的数值模拟及试验

为研究固体火箭发动机激波诱导矢量控制效率的影响因素,及验证数值模拟方法的准确性,对激波诱导轴对称推力矢量喷管模型进行了壁面测压试验,采用二阶精度Roe格式和k-ωSST两方程湍流模型求解强守恒型Navier-Stokes方程对矢量喷管复杂干扰内流场进行数值模拟。根据试验和数值模拟结果分析了喷管内主流和次流相互作用产生的复杂流场结构,比较了在不同喷管落压比NPR和次主流压力比SPR下喷管壁面静压的分布情况。结果表明,数值计算和试验结果基本吻合,验证了计算方法的准确性;在一定范围内减小喷管落压比,增大次主流压力比可以增大喷管周向壁面静压差,提高喷管的推力矢量偏角。     

螺旋式发射管内火箭燃气冲击流场数值模拟

用改进的MUSCL格式解三维、可压缩平均雷诺纳维尔-斯托克斯方程组,湍流模型为Spalart-Allmaras代数模型,计算和分析了火箭发射管内燃气冲击流场.在建立差分算法时,将有限体积离散和黎曼解算器相结合,简化了计算工作.首先,以小喉截面超声速喷管流动问题开展数值实验,数值结果与实验数据吻合很好;其次,计算了火箭发射管内燃气流场,获得发射管内气流压强和温度等物理参数分布.     

三维数值模拟再生冷却喷管的换热

为了解液体火箭发动机喷管再生冷却的换热特点，采用数值模拟的方法，对内喷管燃气、壁面和冷却液建立不同的三维控制方程，进行流动和传热的耦合计算。在计算中，假定喷管流动为冻结流动，考虑燃气向壁面的对流换热和辐射换热；采用二阶迎风格式离散控制方程，采用DO模型离散求解辐射换热方程，水蒸气的吸收系数根据Leckner公式计算。计算模型采用缩比热试车发动机，数值计算结果与实验结果吻合较好，较准确地模拟出了喷管的壁面热流密度，得到了喷管燃气和冷却液的流场和温度场，对高压再生冷却喷管的设计具有指导意义。     

跨音速喷管流场波系分布的数值仿真

应用分离系数矩阵方法对跨音速喷管流场进行数值仿真, 同时计算喷管内马赫波分布。采用5151个网格点的计算结果清晰反映出喷管内膨胀波和压缩波的生成和反射。数值仿真与实验结果在宽广马赫数范围内吻合一致, 并揭示了波系分区分布, 反映出波系影响下跨音速喷管流场参数分布的规律性。      

激波诱导控制推力矢量喷管实验及数值计算

采用实验方法,通过在二元收敛-扩张喷管扩张段引入二次流喷射,开展了激波诱导控制的流体推力矢量技术研究.实验过程通过喷管上、下壁面压力测量及出口射流纹影观测,研究了主流压力、二次流喷射压力以及二次流喷嘴几何(缝或孔)对推力矢量喷管性能的影响.同时,结合数值计算方法,对各实验工况下的喷管流场进行数值模拟,获得了实验手段难以得到的流场数据和性能,对实验结果进行了辅助分析.初步研究结果表明:在给定的实验条件下,主流压力越高,喷管推力矢量角越小,同时推力系数越大;二次流压力越高,喷管推力矢量角越大,同时推力系数减小;同孔喷射相比,采用喷缝几何下的上壁面激波诱导分离点更趋于向上游移动,分离点后压升显著,射流穿透能力强,对主流的扰动强烈.      

超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

固体火箭喷管颗粒尺寸分级两相跨音速流场计算

本文中对固体火箭喷管颗粒尺寸分级的两相跨音速流场作了计算.气相控制方程采用隐式近似因子分解法求解,尺寸分级的颗粒控制方程采用特征线求解,然后,二者进行充分的耦合,可以获得固体火箭发动机含有任意颗粒质量分数和不同颗粒尺寸时轴对称喷管跨音速流场的参数分布.文中讨论了不同颗粒半径和质量分数对流场的影响,对单一颗粒尺寸和颗粒尺寸分级的参数进行了比较.两相耦合计算的迭代收敛速度取决于气相,本文中气相方程求解的格式除部分边界外是隐式的,CFL数可取至6左右,收敛速度快.特别是对颗粒尺寸分级的计算,得益更大,其得益的倍数为颗粒的分级数.