二元矢量喷管内流特性的数值模拟

应用多重网格ＭａｃＣｏｒｍａｃｋ显示格式对不同矢量角下二元喷管的流动进行了模拟分析。计算时,在细网络上采用了雷诺平均Ｎ－Ｓ方程和Ｋ－ε湍流模型方程;在粗网格上则采用Ｅｕｌｅｒ方程。数值计算给出了较清晰的流谱。壁面压力分布的计算值也与实验结果吻合。通过计算结果分析了喷管几何参数对流场的影响。所得结果显示出多重网格算法高效可靠,可用于开展参数研究     

轴对称矢量喷管内流特性的数值模拟研究

利用矢通量分裂有限体积格式和B-L湍流模型,采用TTM网格生成和局部网格加密技术,求解三维N-S方程,对轴对称矢量喷管的内流流场进行数值模拟,并将计算结果与模型试验结果进行了分析比较。研究结果表明,在矢量状态下,计算结果和试验结果相比,壁面静压分布的相对误差不大于13%,气动矢量角的相对误差不大于10%,推力系数的相对误差不大于2%,所开发的计算方法和计算程序可用于轴对称矢量喷管的工程设计。      

轴对称矢量喷管外流对内流干扰研究

应用有限体积法在同位网格中对控制方程进行离散化,对某一轴对称矢量喷管的内外流场进行了数值模拟,数值模拟和试验结果进行了比较,分析了该喷管在矢量偏角20°下,不同的落压比、不同的外流速度对轴对称矢量喷管内流的影响。研究结果表明,喷管内的分离激波是影响轴对称矢量喷管性能的重要因素。      

机械式轴对称矢量喷管动态偏转的矢量特性研究

采用限体积法在同位网格中对满足几何守恒定律的控制方程进行离散化,对某一轴对称矢量喷管在几何偏角速度大小为40°/s下从0°偏转到14.2°过程中的内外流场采用压力隐式分离算法(PISO)进行了数值模拟,分析了该喷管在不同的落压比、不同的外流速度下偏转时,喷管内外流场的动态变化特性,几何偏角和气动矢量角的关系,当喷管偏转后在内壁面产生分离时,喷管的气动矢量角振荡。研究结果为矢量喷管的设计和发动机控制技术的研究提供了参考。      

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

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

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

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

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

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

高超声速喷管三维热化学非平衡流场的数值模拟

利用数值求解三维热化学非平衡Navier-Stokes方程的方法和多块网格技术,数值模拟了典型高焓激波风洞锥型喷管的三维非平衡流场,分析了热化学非平衡效应对喷管流动的影响,给出了平动温度、组分N2和O2振动温度、马赫数和组分质量分数在轴对称线上的分布,也给出了喷管出口附近的速度和温度沿径向分布等结果,获得了喷管流场的详细信息,并与已经存在的采用轴对称Navier-Stokes方程计算的结果,进行了比较,二者吻合得很好.计算结果表明:目前的代码是可以模拟多块结构网格下的热化学非平衡流动的.     

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

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

激波诱导二元矢量喷管内流特性数值研究

采用三种不同湍流模型分别对激波诱导二元矢量喷管内流场进行计算,并将计算结果与模型试验结果进行了分析比较,得出重整化群(RNG)k-ε二方程模型为文中数值模拟中最有效的湍流模型。对该型喷管的计算结果分析表明,喷管矢量角的计算值和实验值基本吻合;在NPR=4,ω=9.6%情况下,矢量角最大可以达到16.61°;注气位置、注气压力和喷管落压比对矢量角的影响比较大。      

A 3D hybrid grid generation technique and a multigrid/parallel algorithm based on anisotropic agglomeration approach

A hybrid grid generation technique and a multigrid/parallel algorithm are presented in this paper for turbulence flow simulations over three-dimensional (3D) complex geometries. The hybrid grid generation technique is based on an agglomeration method of anisotropic tetrahedrons. Firstly, the complex computational domain is covered by pure tetrahedral grids, in which anisotropic tetrahedrons are adopted to discrete the boundary layer and isotropic tetrahedrons in the outer field. Then, the anisotropic tetrahedrons in the boundary layer are agglomerated to generate prismatic grids. The agglomeration method can improve the grid quality in boundary layer and reduce the grid quantity to enhance the numerical accuracy and efficiency. In order to accelerate the convergence history, a multigrid/parallel algorithm is developed also based on anisotropic agglomeration approach. The numerical results demonstrate the excellent accelerating capability of this multigrid method.     

直升机旋翼悬停流场的粘性数值模拟

在旋转坐标系下求解可压缩NS方程,对Caradona-Tung旋翼的跨声速悬停流场进行了数值模拟,计算中采用嵌套网格方法,并且应用了多重网格技术加速流场收敛。计算结果表明,多重网格技术显著改善了流场的收敛效率,计算量仅为单重网格的四分之一。文中同时还对无粘和粘性流场进行了对比,NS方程考虑了物面的粘性作用,能正确模拟出激波和边界层的干扰现象,计算得到的桨叶表面压力分布与实验值更加吻合。     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

聚合多重网格法在旋翼前飞流场计算中的应用

基于非结构动态嵌套网格, 将聚合多重网格方法引入到旋翼的前飞流场计算中.从初始的整体非结构网格出发, 快速建立相互重叠的两个子域网格, 同时无需建立背景网格便可以确定初始洞边界和初始插值贡献单元, 且动态嵌套时无需更新洞边界, 并采用了一种无需背景网格便可快速查找新的贡献单元的方法.分区利用聚合法生成多重网格, 利用双时间方法分区运用三层网格V循环求解非定常Euler方程, 两区的流场参数在重叠区实现耦合, 对旋翼的前飞流场进行了计算.结果表明, 多重网格方法能够有效地加快旋翼前飞流场收敛.      

Development process of muzzle flows including a gun-launched missile

Numerical investigations on the launch process of a gun-launched missile from the muzzle of a cannon to the free-flight stage have been performed in this paper. The dynamic overlapped grids approach are applied to dealing with the problems of a moving gun-launched missile. The high-resolution upwind scheme(AUSMPW+) and the detailed reaction kinetics model are adopted to solve the chemical non-equilibrium Euler equations for dynamic grids. The development process and flow field structure of muzzle flows including a gun-launched missile are discussed in detail.This present numerical study confirms that complicated transient phenomena exist in the shortly launching stages when the gun-launched missile moves from the muzzle of a cannon to the freeflight stage. The propellant gas flows, the initial environmental ambient air flows and the moving missile mutually couple and interact. A complete structure of flow field is formed at the launching stages, including the blast wave, base shock, reflected shock, incident shock, shear layer, primary vortex ring and triple point.     

旋翼/机身气动干扰的数值模拟

通过求解粘性非定常可压缩N-S方程,发展了旋翼/机身气动干扰的数值模拟方法,主要技术措施包括:(1)通过ADT数据结构,建立了高效的多块对接网格动态重叠技术;(2)在可压缩方程求解中引入低速预处理方法,克服低速流场求解时的刚性问题;(3)加入了多重网格方法,提高收敛效率;(4)湍流模型方程求解需要反复计算壁面最近距离,采用Wigton的优化算法很大程度上提高了计算速度。利用建立的计算方法对佐治亚大学(GIT)旋翼/机身干扰模型进行了数值模拟,分别比较了机身表面的平均压力分布、旋翼的下洗流场以及非定常压力脉动,对干扰现象及其流动机理进行了分析。     

跨声速轴流风扇转子流场三维粘性数值分析

采用一种快速求解三维粘性流场的计算方法求解跨声速压气机风扇转子内部流场。该方法以ＬＵ型隐式格式和ＭＵＳＣＬ ＴＶＤ迎风格式为基础,结合壁面函数方法和简单的混合长度湍流模型,使用多重网格迭代加速收敛技术,对三维可压缩雷诺平均Ｎ－Ｓ方程进行求解。计算得到了ＮＡＳＡ Ｌｅｗｉｓ ６７跨声速、低展弦比轴流风扇转子的性能曲线,并重点分析了最高效率点工况和近失速工况下的内部流场。计算与实验结果的对比表明,此方     

二维结构/非结构混合网格的生成及其流场模拟

采用分区求解方法在二维结构/非结构混合网格中对粘性绕流作了数值模拟。利用层面推进法生成绕物体的O形结构网格,利用推进阵面法在结构网格外生成非结构网格。在结构网格区内,采用B-L湍流模型,求解雷诺平均N-S方程。在非结构网格区内求解Euler方程。两区的解通过在交界面上满足通量守恒条件进行耦合。数值算例表明了方法的有效性和实用性。