二元喷管跨音速流场的数值计算

本采用有限体积法进行空间离散，Runge-Kutta法进行时间推进的办法，通过求解Euler方程来模拟二元收敛－扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡，在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解，计算结果与试验数据和其他数值方法吻合良好。     

二元矢量喷管流场的N－S方程计算

对二元矢量喷管在不同工况下的流场进行了数值模拟。求解的控制方程为三维非定常Ｎ－Ｓ方程。空间方向用有限体积离散，时间方向用Ｒｕｎｇｅ－Ｋｕｔｔａ法推进。粘性项的计算上采用了一种有限体积框架下新的离散化方法。湍流模型采用了广泛使用的Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模型。计算表明，本文方法能准确地计算二元矢量喷管、非矢量喷管在设计与大设计落压比下的流场。     

用可变扩散湍流模型研究轴对称排气系统红外辐射特性

为研究可变扩散湍流模型在排气系统红外辐射特性计算中的适用性,计算了轴对称亚声速喷管在3 ～5μm波段的空腔-喷流组合红外辐射特性.喷管流场及温度场采用有限体积法求解N-S方程,湍流模型采用标准k-ε模型和可变扩散模型.红外特性计算采用有限体积法求解吸收-发射性介质条件下的三维辐射传输方程,并考虑大气衰减作用.计算结果与试验对比表明,与标准k-ε模型相比,可变扩散湍流模型能够有效改善亚声速热喷流的核心区长度,从而获得与试验结果一致的红外辐射强度,说明可变扩散湍流模型适用于精确模拟轴对称排气系统的红外辐射特性.     

二元喷管跨音速流场的数值计算

本文采用有限体积法进行空间离散,Runge-Kutta法进行时间推进的办法,通过求解Euler方程来模拟二元收敛—扩散喷管跨音速流场。运用隐式残差平均技术和当地时间步长法加速计算收敛。为了消除气流参数的波动和激波前后的振荡,在方程中添加了自适应耗散项。在100个时间步以内便可获得稳态解。计算结果与试验数据和其他数值方法吻合良好。      

一种高效的多重网格三维N-S方程计算方法及其应用

介绍一种基于 Jameson中心差分和 Runge-Kutta时间推进的三维 N-S方程计算的有限体积方法 ,计算点选取体积单元中心 ,为了提高收敛效率和程序稳定性 ,采用了当地时间步长和多重网格方法 ,特别是针对粘性计算中计算网格形状比高的特点 ,隐式残差光顺采用当地变系数的方法 ,并专门设计了考虑计算网格形状比的人工粘性模型。应用本文方法对收扩喷管的三维流场和NACA0 0 1 2翼型二维流场进行了计算 ,结果表明程序具有良好的收敛能力和较高的稳定性。     

轴对称收敛——扩散喷管几何参数优化设计

本文通过对轴对称收敛－扩散喷管的四个主要几何参数（面积比AR，扩张半角β，收敛半角α，喉部园角半径与喉道半径之比Rt）采用二次回归方程，建立了在设计工况时轴对称收敛－扩散喷管推力系数CF与几何参数AR，β,α.Rt的数学模型，用随机射线法进行优化设计。本方法为排气系统轴对称收敛－扩散喷管方案的选择提供了优化方案。     

多重网格加速的LUSGS算法用于喷管流场数值模拟

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用独特的显式算法克服了两方程湍流模型刚性对流场计算收敛速度的影响,并给出了该方法的理论依据。通过轴对称喷管内流场、二元喷管内外流场的数值模拟验证了采用以上方法在计算效率上的提高。最后将RNG模型计算结果与SA模型计算结果及实验数据做对比分析。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

次流喷射控制推力矢量喷管的流场数值模拟

为了探讨次流喷射控制二维推力矢量喷管的性能,用时间推进求解Ｎ－Ｓ方程的方法数值模拟了喷管的流场,数值格式用中心有限体积格式,借助于当地时间步长和隐式残差光顺技术加速收敛,得到了在不同二次喷流压力下推力矢量喷管内流场的变化特征。计算结果与实验结果比较有满意的一致性。研究表明,应用次流喷射控制主流流动可以实现较大的推力矢量转折,但是,二次喷流必须具有足够的压力值。     

火箭喷管推力偏心高精度计算软件

利用火箭喷管从亚声速、跨声速到超声速的非对称流动数学模型建立了计算喷管推力偏心的理论模型，并采用高精度有限体积TVD格式求解数值方程，可以对火箭发动机喷管的推力偏心特性进行评估与设计，并集成为Nvc软件。该软件完成了从喷管网格生成、数值计算到推力偏心分析全过程，界面友好，使用方便，大大提高了设计效率，通过大量的计算与试验结果的比较，计算结果可以满足工需要，减少试验次数，节约成本，具有广泛的用途。     

轴对称推力矢量喷管的静态内部性能分析

采用Jamson提出的有限体积法对4种结构的轴对称推力矢量喷管的内流场进行了数值分析，预测了其内流场性能参数，计算结果与参考文献符合的很好，表明该方法对建立轴对称推力矢量喷管的数学模型具有很高的模拟精度和较快的收敛速度。     

二元矢量喷管流场及性能数值模拟

用有限差分法对二元矢量喷管的内流场及推力性能进行数值模拟。流场计算采用ＭａｃＣｏｒ－ｍａｃｋ两步格式，运用两层代数紊流模型，对初场、网格进行处理，采用加快收敛的当地时间步长，使计算精度和时间都能较好地满足工程应用的要求，得到了与实验数据较一致的喷管内特性。     

二元喷管流体矢量控制方案数值研究

采用基于雷诺平均的三维N-S方程和RNGk-ε湍流模型对两种流体注入方式下的某型二元收敛-扩散喷管全流场进行了数值模拟.计算结果表明:在激波诱导和喉道倾斜矢量控制方案下, 扩散段射流位置对喷管矢量效率的影响趋势是一致的, 射流位置越靠近喷口, 喷管获得的矢量角越大.相同射流流量在同一位置注入时, 由于喉道倾斜方案下的喷管主流可以实现亚声速偏转, 所以其总压损失较激波诱导方案要小.      

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

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

用温度修正的k-ε模型研究排气系统红外辐射特性

为研究温度修正的k-ε湍流模型对排气系统红外辐射特性计算精度的改善程度,计算了某轴对称收敛喷管在3~5μm波段的空腔-喷流组合红外辐射特性,并与试验结果对比.喷管流场及温度场采用有限体积法(FVM)求解Navier-Stokes(N-S)方程得到,湍流模型分别采用标准k-ε模型和温度修正的k-ε模型.红外辐射特性的计算采用有限体积法求解吸收-发射性介质条件下的三维辐射传输方程,并考虑大气衰减作用.计算结果表明:温度修正的k-ε模型可有效预测排气系统热喷流的核心区长度,从而获得与试验结果一致的红外辐射强度.与标准k-ε模型相比,采用温度修正的k-ε模型计算得到的红外辐射强度的最大相对误差由64%降至了4.3%,说明温度修正的k-ε模型适用于精确模拟排气系统的红外辐射特性.      

长尾喷管两相流流场计算

为了给长尾喷管热防护结构的烧蚀和温度场计算提供准确的热边界条件，从二维轴对称N－S方程出发，采用颗粒轨道模型和有限体积的Jameson格式计算了长层喷管纯气相流场和两相流流场，分析和比较表明：在长尾喷管喉部附近以及扩散段的相同轴向位置，两相流情况下轴线上的燃气马赫数小于纯气相情况下的燃气温度；在进行长尾喷管热防护结构设计时应以两相流流场作为热边界条件。     

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

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

姿控发动机喷管内流场及推力计算研究

运用时间相关法数值求解定常Ｎ－Ｓ方程和采用有限体积方法离散控制方程,对姿控发动机喷管内流场及其产生的推力进行了计算,得到了较为精细的流场图象。计算结果表明,不同壁面温度条件可以显著改变喷管壁面附近的流场结构及其产生的推力。     

喷管激波贴口压比计算

本根据试验数据，采用多元线性回归分析方法建立了轴对称收敛－扩散喷管激波贴口时落压比与面积比，收敛半角，扩张半角，无量纲喉道曲率半径之间关系的数学模型，为喷管研究提供依据。     

收敛-扩散喷管几何参数对内特性的影响

航空发动机广泛采用机械可调式收敛-扩散喷管,正确地选择几何参数是喷管设计的主要任务之一。本文通过模型试验研究具有不同几何参数的收敛-扩散喷管在不同工况下的推力特性和壁面压力分布。试验设备及测量系统见文献。模型喷管有轴对称收敛-扩散型和两侧壁为有机玻璃的二元收敛-扩散型两种,其几何参数为:面积比C_A为1.0～2.4扩张半角α为0～24°,收敛半角β为8～24°,喉部圆角半径R_0与喉道半径R_1之比为0.31～2.0。