非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

一类逆风ENO线方法的构造及其对Euler方程的应用

本文利用ENO(essentially-non-oscillatory)插值思想,对非线性双曲型守恒律=0提出了一类逆风ENO线方法,使用这种格式其计算量比C.Shu在文[6]中提出的通量分裂方法明显减少,进而同Runge-kutta时间离散相结合,提出了时空高精度差分格式。本文还把这类方法推广到方程组并用于计算Euler方程组的黎曼问题和拟一维喷管问题,得到了比较满意的数值结果。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

基于非定常低速预处理和DES的三角翼数值模拟

发展出一套基于低速预处理技术采用SA-DES和SST-DES模型,对三维低速非定常流动进行数值模拟的方法.采用非结构混合网格有限体积方法求解非定常流场,时间离散采用基于LU-SGS隐式格式的双时间步长方法.对非定常低速预处理方法进行了推导和分析,并应用该方法对大攻角三角翼的非定常流场进行了数值模拟,对不同站位压强分布与试验值和参考文献值进行了对比,吻合较好.     

多管脉冲爆震发动机共用喷管特性的数值模拟

基于有限速率基元化学反应模型和二维迎风型TVD格式的数值模拟方法,研究了3种不同形式的共用喷管设计对多管PDE非定常流动的影响,化学反应模型采用H2/O2/N2的11组分23步反应.数值研究的结果表明:共用喷管对多管PDE中流场产生了重要影响;单管点火的情况下,爆震波传播过程中由于喷管内表面作用产生的反射激波回传有利于旁通爆震管的起爆,喷管可以改变点火爆震管封闭端的压力变化历程,收敛扩张喷管的多管PDE封闭端的压力高于收敛喷管和扩张喷管;共用喷管使得各个爆震管的起爆和运转过程中存在相互影响,设计合理的共用喷管对优化多管PDE的性能和保证PDE协调工作具有重要作用.     

嵌套网格粘性流动数值模拟用于风洞洞壁干扰研究

运用嵌套网格技术和Navier Stokes数值模拟对机翼半模和翼身组合体试验时风洞的四壁干扰进行综合模拟、评估和修正。计算格式在空间上采用中心有限体积离散 ,在时间上采用多步Runge Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

基于WENO-分段线性格式的旋翼流场数值模拟

为了降低二阶JST(Jameson-schmidt-turkel)格式导致的数值耗散,发展了一套适合于格心格式的基于加权本质无振荡(Weighted essentially non-oscillatory,WENO)-分段线性格式的旋翼流场数值模拟方法。采用运动嵌套网格方法生成围绕旋翼的网格系统,主控方程选择Navier-Stokes(N-S)方程。为了更加有效地对桨尖涡等流动细节进行捕捉,在笛卡尔背景网格上采用七阶Roe-WENO格式计算对流通量;在桨叶贴体非结构网格上采用二阶精度的Roe-分段线性格式计算对流通量。时间离散采用了高效的双时间隐式LU-SGS(Lower upper symmetric Gauss-Seidel)方法进行时间推进。最后,应用上述方法对悬停状态的C-T(Caradonna-tung)旋翼和Helishape 7A旋翼进行了数值模拟,将数值计算结果与实验数据进行了对比,计算值与实验值吻合较好;并将桨尖涡模拟效果与二阶JST格式的模拟效果进行了对比。对比结果表明:本文方法能有效对旋翼流场进行计算,且在相同计算条件下,WENO-分段线性格式能够更有效地捕捉旋翼涡流场的流动特性,表明在计算旋翼涡流场时WENO-分段线性格式相比传统二阶JST格式具有更低的数值耗散。     

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。     

辐射及化学非平衡流耦合场计算方法

把一种新的三维热辐射计算方法应用到高超声速非平衡流场的数值模拟过程中:在非结构网格上对高超声速化学非平衡流场的数值模拟过程中,耦合了辐射输运方程来求解辐射热流。控制方程为含有化学反应源项和辐射源项的三维N-S方程,数值离散格式采用了Jam eson有限体积法,辐射输运方程采用有限体积法求解。化学反应模型为7组元7反应模型,并采用了化学反应特征时间步长。本文分别对RAM C-Ⅱ和MU SES-C模型进行了计算,计算结果与参考文献的结果吻合。     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

马赫数连续可变跨声速湿蒸汽风洞的研制

用于跨声速气动测量的探针须从亚声速到超声速范围进行标定。变质量槽式喷管通过扩张段壁面上槽缝流出部分气流的自适应特性可在不同背压下得到不同出口马赫数,从而使标定气动探针的风洞实现马赫数从0到超声速的连续变化。为了研究采用湿蒸汽为工质的变质量槽式喷管的性能及优化其结构,采用三维犖-犛方程以及可实现犽-ε湍流模型对其进行了详细的数值仿真。结果表明收缩段型线、扩张段长度及壁槽尺寸等对喷管流场特性有重要影响,喷管进出口压比在一定范围内,槽式喷管有最优的收缩段型线、扩张段长度和开槽尺寸。根据数值仿真结果研制了马赫数从0到1.6连续可变的跨声速湿蒸汽风洞,对此风洞性能进行验证,表明该风洞在马赫数从零到超声速范围内可获得均匀、稳定的出口气流,满足跨声速湿蒸汽气动探针的标定要求。     

NUMERICAL SIMULATION OF THREE DIMENSIONAL INTERACTING TURBULENCE FLOW FIELD OVER PROJECTILE WITH LATERAL JETS

INTRODUCTIONIn recent years,the study of lateral jets hasaroused more interest of aerodynamicist all overthe world.On the one hand,the interacting flowfield caused by lateral jet includes complicatedshock wave patterns and has important values tostuding flow phenomena.It has been not veryclear to this day and it is very difficult for nume-rical simulation.On the other hand,lateral jethas important applications for the field of aero-nautics,astronautics and weapons:for example,lateral jet …     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

大涡模拟在加力燃烧室冷态流场中的应用

对带V型稳定器模型加力燃烧室二维冷态流动进行了大涡模拟。采用了代数和k方程两种亚网格尺度模型，并用SIMPLE算法和混合差分格式求解离散方程，用壁面函数处理固壁边界条件，得到了模型加力燃烧室二维瞬态流场的模拟结果。通过对不同亚网格尺度模型和不同入口速度下的瞬态流场的比较，揭示了不同亚网格尺度模型和不同入口速度分布对带V型稳定器燃烧室二维冷态瞬态流场的影响。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

Model for Asymmetry of Shock/Boundary Layer Interactions in Nozzle Flows

The reason for the asymmetry phenomenon of shock/boundary layer interactions(SBLI)in a completely symmetric nozzle with symmetric flow conditions is still an open question.A model for the asymmetry of nozzle flows was proposed based on the properties of fluid entrainment in the mixing layer and momentum conservation.The asymmetry model is deduced based on the nozzle flow with restricted shock separation,and is still applicable for free shock separation.Flow deflection angle at nozzle exit is deduced from this model.Steady numerical simulations are conducted to model the asymmetry of the SBLIs in a planar convergent-divergent nozzle tested by previous researchers.The obtained values of deflection angle based on the numerical results of forced symmetric nozzle flows can judge the asymmetry of flows in a nozzle at some operations.It shows that the entrainment of shear layer on the separation induced by SBLTs is one of the reasons for the asymmetry in the confined SBLIs.