复杂前机身—进气道组合体的流场计算

本文提出了一种能够对复杂前机身-进气道组合体进行流场计算与气动分析的方法。通过求解Euler方程来模拟前机身-进气道组合体在亚、跨、超声速,不同攻角和发动机流量条件下的各类流场。     

超音速来流中侧向喷流干扰流场的数值模拟

通过求解三维可压缩N S方程组,对尖锥/圆柱组合体模型在超音速绕流中的侧向喷流干扰流场进行了模拟,数值方法采用有限体积法,差分格式采用NND格式。分析了干扰流场的结构,研究了干扰效应对气动性能的影响,得到来流Ma∞=3 3和4 5时力放大因子和干扰力矩系数随攻角变化的规律,计算结果与实验一致。     

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

亚超声速位流场中复杂组合体的压力分布与气动特性的数值计算

本文对F.A.伍德华得的亚超声速位流场中组合体的气动分析方法作了进一步的研究和改进。改进表面边界条件方案中前缘涡面的处理;在翼身交接区内,采用改变机身内部过渡涡位置和纵向分块数的方法,改善了超声速流中组合体上的压力分布;采用了实用的物体表面几何处理方法,扩大了计算程序的几何描述能力。本文简要描述了机翼-机身-尾翼(包括垂尾短舱等部件)复杂组合体的压力分布和气动特性的数值计算方法,介绍了多种简单组合体和复杂组合体在亚超声速位流场中的计算结果,并与实验数据作对比。理论计算和实验数据比较一致,复杂组合体的计算结果表明,计算精度基本上能够满足初步选型设计的要求,可供实际使用。     

倾转过渡状态旋翼-机翼气动干扰特性

针对倾转旋翼机倾转过渡状态旋翼-机翼之间的气动干扰开展了高分辨率数值模拟。计算采用有限体积法对积分形式的非定常Navier-Stokes方程进行离散,使用重叠网格处理旋翼桨叶与机翼之间的相对运动。背景网格采用基于八叉树结构的笛卡尔结构化网格,并使用基于无量纲Q准则的网格自适应技术来获得高分辨率尾迹流场。首先,计算了孤立旋翼和旋翼-机翼组合体悬停状态,根据计算结果证明高分辨率求解器对于旋翼气动力预测的准确性和旋翼尾迹结构演化的高分辨率特性,对比了旋翼-机翼组合体悬停状态改进延迟分离涡模拟(IDDES)和RANS计算结果,IDDES所得流场更为精细且气动力结果与RANS结果也有所差异。随后,对倾转过渡状态不同倾转角下旋翼-机翼组合体飞行状态进行了仿真计算,得到并分析了旋翼与机翼之间气动干扰流场。结果表明,在倾转过渡状态中期旋翼产生诱导滑流对机翼有着一定的升力增益效果,而倾转初期与末期则无明显升力增益,初期甚至为负增益;旋翼与机翼气动变化情况证明倾转过渡状态旋翼-机翼气动干扰对倾转旋翼气动性能有着重要影响。     

高超声速安定面／舵面干扰特性研究

为了改善高超声速乘波体飞行器的稳定性及操作性,针对安定面／舵面组合体的绕流特性进行研究。通过数值求解RANS方程,对Ma=6时安定面／全动舵组合体关键参数对舵效的影响规律进行研究,并对流动机理进行分析。计算结果表明,安定面与全动舵间水平距离、垂直距离是影响舵效的关键参数。分析干扰流场可知,全动舵主要受到安定面尾迹、激波系及翼尖涡的干扰。其中,安定面尾迹对舵面干扰最大,故全动舵应尽量避开安定面尾迹。     

三维复杂无粘流场的自适应八叉树结构直角网格算法模拟

运用自适应八叉树结构直角网格对三维复杂流场进行了划分 ,并采用点面判别法进行外法向量的判断 ,使与之关联的与物面相交网格的体积计算更为简便准确。在此基础上 ,采用以中心差分为基础的 Jameson的有限体积法进行 Euler方程求解 ,在通量计算过程中 ,采用以面为基础的算法 ,减少了计算工作量。对多个外挂、翼身组合体及全机等复杂外形绕流流场进行了模拟 ,均取得了较好的效果     

横向喷流对鸭式导弹滚转特性影响研究

通过数值方法求解三维可压缩雷诺平均N-S方程,对导弹横向喷流的干扰流场进行了数值模拟,计算结果与实验数据吻合较好,基本验证了计算方法在横向喷流复杂流场数值模拟方面的有效性。在此基础上,对固定尾翼鸭式布局导弹亚、跨、超声速流场进行了数值模拟,并且计算分析了横向喷流对鸭式布局导弹滚转控制特性的影响,计算结果表明,横向喷流可以有效提高鸭式布局导弹的滚转控制能力。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

超声速和高超声速横向喷流数值模拟

为了解横向喷流的干扰影响和喷口附近的流场结构,采用分块对接网格和“0”型网格技术,数值求解N—S方程来模拟超声速和高超声速流场中横向喷流的干扰流场。对两种尖拱弹身外形的超声速和高超声速喷流干扰流场进行了数值计算,计算结果与风洞实验数据吻合一致。在此基础上,开展了某型导弹多喷构型干扰流场的数值模拟,得出一些横向喷流数值模拟的结论。     

任意回转面亚、跨音流场计算方法

本文所设计的任意回转面亚、跨音流场的计算软件 ,能根据叶型和进出口条件分析流动特性 ,并自动选择精确和节省机时的合适的数值解法。用户在给出叶型坐标和定解条件后 ,可以得到叶面马赫数分布或压力分布 ,全流场的等马赫数或等压图     

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

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

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

Numerical simulation of multi-dimensional inviscid compressible flows by using TV-HLL scheme

This paper deals with the numerical solution of inviscid compressible flows. The threedimensional Euler equations describing the mentioned problem are presented and solved numerically with the finite volume method. The evaluation of the numerical flux at the interfaces is performed by using the Toro Vazquez-Harten Lax Leer(TV-HLL) scheme. An essential feature of the proposed scheme is to associate two systems of differential equations, called the advection system and the pressure system. It can be implemented with a very simple manner in the standard finite volume Euler and Navier–Stokes codes as extremely simple task. The scheme is applied to some test problems covering a wide spectrum of Mach numbers, including hypersonic, low speed flow and three-dimensional aerodynamics applications.     

二维超燃冲压发动机磁控进气道的数值模拟

针对飞行马赫数大于设计马赫数的情况,采用二维磁流体动力学方程对磁控进气道进行了数值模拟研究。数值模拟结果表明磁流体装置的电磁作用可以使非设计马赫数下进气道激波满足SOL(shock on lip)条件,并使出口处的流动变得均匀。分析了磁流体作用位置、宽度和深度等关键控制参数对该类进气道性能的影响,计算结果表明,磁流体作用区域越靠近飞行器前缘,而且在纵向越深,则进气道出口处的流动越均匀,但流率会有所下降;若磁流体作用区域较宽,则需较小的磁场就能使非设计马赫数下进气道的激波结构满足SOL条件。     

高超声速弯曲激波压缩侧压式进气道数值研究

采用壁面马赫数呈线性分布的曲面压缩系统改进参考侧压式进气道的顶板,得到弯曲激波压缩侧压式进气道,并与参考侧压式进气道进行了比较.数值研究结果表明:设计状态无黏时曲面压缩顶板壁面马赫数分布与给定的马赫数分布基本一致,并且有黏时其壁面压力分布也与二维曲面的基本相同;同参考侧压式进气道相比,顶板采用曲面压缩能够一定程度地改善壁面压力分布,使其末端压力梯度变化平缓;并且非设计状态下的性能也得到有效地改善,特别是来流马赫数为4时,其流量系数提高6.0%、达到0.799,喉道截面总压恢复系数提高1.9%;来流马赫数为5时,其流量系数提高5.2%、达到0.909,喉道截面总压恢复系数提高3.2%.随着攻角增大,该进气道流量捕获能力增强、隔离段出口截面流场畸变减小,但喉道截面总压恢复系数下降剧烈.      

Simulation of underexpanded supersonic jet flows with chemical reactions

To achieve a detailed understanding of underexpanded supersonic jet structures influenced by afterburning and other flow conditions, the underexpanded turbulent supersonic jet with and without combustions are investigated by computational fluid dynamics(CFD) method.A program based on a total variation diminishing(TVD) methodology capable of predicting complex shocks is created to solve the axisymmetric expanded Navier–Stokes equations containing transport equations of species. The finite-rate ratio model is employed to handle species sources in chemical reactions. CFD solutions indicate that the structure of underexpanded jet is typically influenced by the pressure ratio and afterburning. The shock reflection distance and maximum value of Mach number in the first shock cell increase with pressure ratio. Chemical reactions for the rocket exhaust mostly exist in the mixing layer of supersonic jet flows. This tends to reduce the intensity of shocks existing in the jet, responding to the variation of thermal parameters.     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

弹体穿越平面激波三维非定常流场的数值模拟

用数值方法求解三维非定常可压缩Navier-Stokes方程，计算了复杂弹体在相同马赫数、不同迎角下的流场和气动力特征，计算结果与实验结果吻合，证明本文采用的方法能够准确计算此类复杂流场；然后采用该方法模拟了弹体穿越正激波的非定常过程，得到了弹体周围流场及气动力特征随穿越过程的非定常变化，阻力特征穿越后急剧下降。