轴对称欠膨胀羽流流场的分离系数矩阵解法和实验研究

本文对轴对称超音速羽流流场进行了数值计算和实验研究。数值计算方法是空间步进的分离系数矩阵差分法,并用粘性和无粘分裂法处理粘性问题。同时对同轴超音速喷流进行了流场观察和测量,所得纹影照片清晰地显示了尾喷流场的详细结构。计算结果在近场与文献[1]一致,在远场与本文的实验结果基本一致。计算实践表明,与中心差分格式相比,分离系数矩阵方法的稳定性好,精度高,捕捉的激波耗散小,且无数值振荡。     

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

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

复杂几何形状喷管内外三维流场的数值模拟

研究了球面收敛矩形扩张的二元喷管（ＳＣＦＮ）流场的数值模拟,为了有效地捕获激波与滑流面及较好地模拟分离,采用了Ｒｏｅ格式、ＭＵＳＣＬ方法及Ｖｅｎｄａｔａｋｒｉｓｈｎａｎ限制器、Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ方程紊流模型。研究中给出了两个喷管模型四个流动状态的数值模拟结果。结果表明,除了有严重分离的流动情况以外,该算法能得到相当满意的计算结果。     

轴对称矢量喷管内流特性的不同湍流模型计算

利用二阶迎风格式，在采用三种不同湍流模型分别对轴对称矢量喷管三维内流场进行对比计算后，S-A-方程模型成为文中数值模拟中最有效的湍流模型，并将此模型计算的结果与模型试验结果进行了分析比较。研究结果表明，在矢量状态下，计算结果和试验结果相比，气动矢量角的相对误差不大于6％，流量系数的相对误差不大于l％，推力系数的相对误差不大于l％，所开发的计算方法可用于轴对称矢量喷管的工程设计。     

固体火箭发动机喷管喉衬动边界条件下的温度场计算

本文探讨了固体火箭发动机喷管喉村在烧蚀退移、燃气中固体微粒的沉积以及沉积物吹除所引起的动边界条件下的温度场数值计算模型。模型的建立以有限体积能量平衡的基本原理为依据,采用有限差分的数值计算方法及稳定性较好的交替方向隐式格式求解。为确保数值计算结果在边界退移及网格脱落时的稳定性,采用了插值与迭代相结合的方法,以全隐式格式求解表层网格节点和受热壁面的温度。模型还考虑了多层背村材料结构、变物性材料特性等诸种因素,使之具有一定的普遍性,并可应用于工程设计和研究中。数次发动机实验证明,本文提出的温度场计算模型能够较真实地反映喉衬温度响应,具有较为满意的计算精度。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

高超声速再入飞行器气动轴研究

根据根据力学中的力螺旋概念，论述了再入飞行器的气动轴存在的一般性，分析了再入飞行器气动轴的三种类型，提供了确定气动轴的方法，建立了计算非对称再入飞行器飞动轴的近似方法和数值方法，近似方法基于修正的牛顿流理论和推广的内伏牛顿流理论，数值方法基于Ｊ．Ｅ．Ｄａｙｗｉｔｔ教授发展的分离系数矩阵欧拉方程时间相关方法和张鲁民研究员发展的三维定常分离系数矩阵欧拉方程的空间推进方法。     

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

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

入口温度剖面对喷管流场结构的影响

应用质量平均的Navier-Stokes方程和B－L代数湍流模型，对超燃冲压发动机尾喷管的流场进行了数值模拟研究，在计算过程中，对方程中的对流项采用了空间为二阶精度的TVD格式，扩散项则采用了二阶中心差分离散，通过数值模拟，对比研究了温度非均匀性，三维效应对尾喷管的流场结构的影响。     

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

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

倒置喷管的湍流流场计算

数值模拟采用基于 Van Leer的矢通量分裂格式 ;应用线 Gauss- Seidel迭代法求解 ;并结合局部时间步长技术对收敛过程加速。对 Re=1 .66× 1 0 6 的倒置喷管流场进行数值计算的结果与实验符合     

A low-diffusion robust flux splitting scheme towards wide-ranging Mach number flows

This paper develops a low-diffusion robust flux splitting scheme termed TVAP to achieve the simulation of wide-ranging Mach number flows. Based on Toro-Vázquez splitting approach, the new scheme splits inviscid flux into convective system and pressure system. This method introduces Mach number splitting function and numerical sound speed to evaluate advection system. Meanwhile, pressure diffusion term, pressure momentum flux, interface pressure and interface velocity are constructed to measure pressure system. Then, typical test problems are utilized to systematically assess the robustness and accuracy of the resulting scheme. Matrix stability analysis and a series of numerical cases, such as double shear-layer problem and hypersonic viscous flow over blunt cone, demonstrate that TVAP scheme achieves excellent low diffusion, shock stability, contact discontinuity and low-speed resolution, and is potentially a good candidate for wide-ranging Mach number flows.     

分开排气系统的1种设计方法及其性能研究

针对大涵道比涡扇发动机,开展了其排气系统气动型面参数化设计方法和气动性能的研究。通过控制外涵、内涵喷管流道的中心线形状和流通面积,设计了外涵、内涵分开排气喷管气动型面。采用数值模拟的方法模拟了喷管的流场结构,并分析了外涵落压比和自由流马赫数对喷管推力系数的影响。在自由流马赫数为0.050时,推力系数随外涵落压比的增大而增大,在自由流马赫数为0.785时,推力系数随外涵落压比的增大而减小;在外涵落压比相同时,自由流马赫数越大,核心机舱外罩b段与尾锥b段所受轴向力在喷管总推力中所占比例越大,对喷管推力系数的影响也增大。     

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

钝体绕流流场分离涡的数值模拟

 采用近似因式分解和特殊Jacobin系数矩阵分裂法,求解可压缩Navier-Stokes方程,数值模拟矩形柱体绕流流场中分离涡的形成和发展。给出了Ma_∞=0.3,0.8;Re=10~4,10~6的计算结果。从结果中可清楚地看到柱体非对称涡形成的非定常过程。     

多点集中力下高超进气道弹性变几何研究

针对二元高超声速进气道需在宽马赫数Ma为4.5～6.0的范围内工作的要求,探索了一种在多点集中力作用下进气道曲面压缩面弹性可调的变几何方案.通过数值计算,首先对集中力数目的选取进行了研究,发现3个集中力对型面变形更为有效;然后对在3个集中力作用下产生变形的进气道进行了流固耦合分析,结果表明变形型面的静压分布与理想型面的静压分布基本吻合;对其性能进行分析,结果表明:Ma为6.0设计的进气道,变形后在Ma为4.5,5.0流量系数分别达到0.98和1.00,且出口总压恢复系数与未变形的进气道相比基本保持不变.由此说明,多点集中力下弹性可调的变几何方案是可行的,并有助于提高非设计点下进气道的性能.