导弹多喷口侧向喷流流动数值模拟

针对典型导弹外形,通过求解三维Navier-Stoke方程开展了多喷口喷流与超声速来流的复杂干扰流动数值模拟,研究了来流马赫数、来流攻角、喷口数量及位置、喷流压力比对流场结构和弹体表面流动的影响.结果表明,相邻喷口间会产生明显的干扰,来流马赫数、攻角、喷口数量及位置和喷流压力都会对喷流与主流之间的干扰以及导弹表面流动产生明显影响,且在沿流向布置的多喷口中,最上游喷口决定了喷口前的流场结构.     

自然超空泡与尾喷流相互作用的数值模拟

认识带尾喷流和自然超空泡的水下高速航行体流体动力特性并发展其预报与控制方法仍是水动力学领域极具挑战性的课题．本研究采用CFD方法对尾喷流和自然超空泡之间的相互作用进行了数值研究．针对发动机欠膨胀超音速喷流,采用现有实验结果验证了基于两方程湍流模型的二维轴对称流动数值模型的可靠性．尾喷流在空气和蒸汽环境中流动的数值计算结果表明,由于蒸汽环境中背压较低,欠膨胀尾喷流无法及时形成压缩波,使得蒸汽环境中尾喷流的过膨胀区和气相扩散区的体积比空气中大;尾喷流很难形成规则的激波格栅,波系结构相对简单．针对携尾喷流的平头航行体超空泡流状态的数值模拟结果表明,尾喷流注入超空泡后迅速充满航行体周围的空腔区域;尾喷流与超空泡尾迹区域形成的回射流相互作用最终导致超空泡断裂,断裂过程中伴随着燃气泡的下泄现象;受空泡壁面约束,尾喷流难以在狭窄的超空泡空腔内完全膨胀,尾喷流的激波波系结构有显著的变化：在喷嘴附近受到压缩,在远离喷嘴区域受到超空泡水汽掺混的破坏;空泡内压强基本维持在饱和蒸汽压附近,没有显著增加．     

二维振荡叶栅非定常粘性流动数值模拟

采用显式四步Runge-Kutta格式，结合Baldwon-Lomax紊流模型求解Navier-Stokes方程，借助运动网格技术，完成了对二维振荡叶栅非定常粘性流动的数值模拟。为了加速求解过程，引入了变系数隐式残差光顺方法，取得了较好效果。数值结果与已公布的数据有很好的一致性。     

超声速钝体逆向喷流减阻的数值模拟研究

为研究逆向喷流技术对超声速钝体减阻的影响,采用标准k-ε湍流模型,通过求解二维Navier-Stokes方程对超声速球头体逆向冷喷流流场进行了数值模拟,并着重分析了喷口总压、喷口尺寸对流场模态和减阻效果的影响。计算结果显示:随着喷流总压的变化,流场可出现两种流动模态,即长射流穿透模态和短射流穿透模态;喷流能使球头体受到的阻力明显减小;存在最大减阻临界喷流总压值(在所研究参数范围内最大减阻可达51.1%);在其它喷流物理参数不变时,随着喷口尺寸的增大,同一流动模态下的减阻效果下降。本文的研究对超声速钝体减阻技术在工程上的应用具有一定的参考价值。     

微差爆破震动叠加起始位置数值模拟

采用ANSYS/LS-DYNA3D程序对岩体中双孔微差爆破近源场爆破振动进行了模拟。结果表明,随着离爆源距离的增加,质点振动波形的持续时间变长,同时二炮孔爆破振动波形中质点处于静止状态的间隔变小并产生叠加。因此通过数值模拟可确定微差爆破近区振动叠加起始位置。     

非定常IV型激波－激波干扰数值模拟研究

对IV型激波-激波干扰非定常流动进行了数值模拟，采用有限体积法，结合空间半隐的二阶OC-TVD格式与时间二阶显式Runge-Kutta法求解三维全N-S方程，并且使用了Baldwin-Lomax代数湍流模型．得到了周期性流场变化结果，其中包括周期变化的双涡结构．壁面压强峰值的大小和位置均呈周期性振动，壁面压力系数和Stanton数的时均分布与定常实验结果符合得较好．并从一周期内流场结构的扰动传播出发，分析了结构变化的相位，说明了IV型激波一激波干扰内在的非定常性机理与影响因素。     

加权基本无振荡格式研究进展

加权型基本无振荡WENO格式是近十年发展起来的一类高阶、高精度格式,它是在ENO格式的基础上采用加权思想构造的,对流场内的间断和细致结构具有较高的分辨率,适于求解包含激波、膨胀波以及接触间断等复杂结构的流场,目前已发展成为计算流体力学中一种重要的方法。本文针对加权型基本无振荡格式近年来的进展作一简要介绍。     

带控制片方柱绕流的非定常数值模拟

对高雷诺数下带控制片的方体(在方柱之前放置一小尺度的柱形薄片)非定常绕流进行了数值模拟。数值模拟方法采用RNG重整化群紊流模型，SIMPLE算法，在非定常计算中引入双重时间步方法，将方体和控制片区域作为求解域的一部分作整体求解。结果表明，控制片与方柱之间距离的不同，使得控制片所形成的涡被抑止或产生卡门涡街，从而对柱体侧面边界层的分离发生影响，产生三种典型的流态，柱体时均阻力系数相对于无控制片的情况迅速减小，当控制片偏离柱体轴心时，柱体升力系数迅速增大，偏心至柱体外壁面附近达到极大值。阻力系数达到最小值后将产生跳跃式的增大，而升力系数在达到最大值后也将产生跳跃式的减小。     

高超声速稀薄流中横向喷流干扰特性实验研究

喷流干扰是高超声速飞行高精度控制的一种有效手段,研究者们以往大部分都主要集中于连续流条件下喷流干扰效应的机理研究,并给出了喷流干扰流场的典型结构,而稀薄流条件下喷流干扰特性的实验数据还十分匮乏.本文利用JFX爆轰激波风洞产生高超声速稀薄自由流,基于平板模型开展不同喷流压力和自由来流参数对横向喷流干扰特性影响的实验研究,采用高速纹影成像及图像处理技术,获得稀薄流条件下喷流干扰流场演化过程及流场结构的变化规律.相比于无喷流条件形成的流场,横向喷流与稀薄自由流相互作用形成的流场结构更为复杂,喷流压力由于受到稀薄来流的扰动,斜激波会短暂穿透喷流干扰流场并延伸至楔形体上部.喷流干扰流场内桶状激波的影响范围随着喷流压力的升高而逐渐变宽,位于三波点上游的斜激波空间位置不会随喷流压力的变化而改变,而位于三波点下游的弓形激波则向上游移动,当喷流压力过低时,桶状激波不会与其他两种激波交汇形成三波点.高超声速稀薄来流压力的降低同样会使桶状激波的影响范围变宽,弓形激波同样也会向上游移动,但基本不会对斜激波空间位置产生任何影响.     

用基于M-SST模型的DES数值模拟喷流流场

脱体涡数值模拟方法(dettached eddy simulation,DES)是把雷诺平均Navier-Stokes方程(RANS)方法及大涡模拟方法(LES)结合起来模拟有脱体涡的湍流流场的数值模拟方法,其主要思想是在物面附近解雷诺平均Navier-Stokes方程、在其他区域采用Smagorinski大涡模拟方法。本文在剪切应力传输(SST)湍流模型的基础上用DES及混合非结构网格数值模拟具有横向喷流的湍流流场,算法采用Osher逆风格式,利用该套程序(包括网格生成及算法),对导弹在不同马赫数下的喷流流场进行了数值模拟,并与同时开展的实验研究的结果进行了对比,结果表明用该方法处理这类问题是较准确的。     

多级轴流压气机叶栅内三维紊流流场的数值模拟

对多级轴流压气机叶栅内三维紊流流场进行数值模拟,采用高精度高分辨率的三阶ENN格式以保证对激波的捕捉和对紊流特征的正确模拟,利用LU-SGS隐式解法提高了计算速度,从而构成了一种既准确又高效的多级跨声速轴流压气机紊流流动数值求解系统。重点研究了动静网格的交接方法及相应的动静交接面处理模型。计算了某带进口导叶的三级轴流压气机,计算结果与实验数据吻合较好。     

爆炸与冲击问题的大规模高精度计算

爆炸与冲击问题的数值模拟在国防和民用安全领域具有重要的工程实用价值.由于爆炸与冲击问题是一个多物质在高应变率、高温及高压条件下的强非线性的瞬态动力学问题,给数值模拟带来了很多的困难,为此,针对爆炸与冲击问题数值模拟中的一些关键和难点问题开展了研究.提出了三维非线性双曲守恒系统的伪弧长自适应网格算法,分析了算法的实现过程,数值结果表明该算法有效地提高了冲击波强间断处的分辨率.发展了针对气相爆轰数值模拟的附加龙格-库塔方法,对非线性对流项进行显示计算,化学反应源项进行半隐式计算,有效地解决了源项引起的刚性问题,计算结果表明该算法可以准确地捕捉和描述爆轰波的复杂结构和典型特征.针对三维工程实际物理问题中的大规模计算需求,给出了三维多物质流体动力学欧拉数值方法的并行化方法,开发了三维爆炸与冲击问题并行计算程序,并给出了针对该并行程序的测试方法.上述工作有利地解决了爆炸与冲击问题大规模、高精度计算中的一些难题.最后,开展了大口径聚能射流侵彻混凝土靶问题的数值模拟和实验研究,通过典型爆炸与冲击工程问题的计算验证了所研究数值方法的有效性.      

Numerical Simulation of a Supersonic Flow with Transverse Injection of Jets

A plane supersonic flow with symmetric perpendicular injection of jets through slots in the walls is numerically simulated with the use of Navier–Stokes equations. The effect of the jet pressure ratio and Mach number on the flow structure is considered. The angle of inclination of the shock wave and the separationregion length are found as functions of the jet pressure ratio. The influence of the jet pressure ratio on the increase in the lift force arising owing to interaction of the flow with the injected jet is found.     

基于混合型WENO格式数值模拟四类Riemann问题

本文基于三阶WENO格式和三阶WENO-Z格式,利用有限体积法研究了同一格式在不同方向使用不同的求解器,以及在不同方向采用相同精度的不同格式数值模拟4类Riemann问题,分析各向异性对数值计算结果的影响．数值结果表明,无论是使用不同的格式,还是使用不同的求解器进行数值模拟,都会导致数值结果不同程度地失去图像对称性．     

半封闭狭缝湍流冲击射流的数值模拟

将Yakhot和Orszag提出的RNGk-ε模型推广应用于半封闭狭缝冲击射流场的数值模拟，以评价该模型对这种复杂湍流的预测能力。将计算得到的流场平均速度分布、湍流强度分布和流函数分布与标准k—ε模型的预测结果以及相应的实验数据进行了比较，结果表明：RNGk—ε模型的预测结果总体上要好于标准k—ε模型，但与实验值相比，所有预测结果都还存在不同程度的误差，尤其是近壁区和滞止点较远下游处的湍流强度分布。说明RNG模型虽然已在某些湍流的预测中取得了一定的成功，但要定量准确地预测冲击射流场，还必须针对其流动特征对模型加以改进。     

FORCED OSCILLATIONS OF BOUNDARY VALUE PROBLEMS OF HIGHER ORDER FUNCTIONAL PARTIAL DIFFERENTIAL EQUATIONS

ＦＯＲＣＥＤＯＳＣＩＬＬＡＴＩＯＮＳＯＦＢＯＵＮＤＡＲＹＶＡＬＵＥＰＲＯＢＬＥＭＳＯＦＨＩＧＨＥＲＯＲＤＥＲＦＵＮＣＴＩＯＮＡＬＰＡＲＴＩＡＬＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＪｉｎＭｉｎｇＺｈｏｎｇ（靳明忠），ＤｏｎｇＹｉｎｇ（董莹），Ｌｉ...     

Numerical Study of a Forced Pendulum with Friction

We first describe the model of a forced pendulum with viscousdamping and Coulomb friction. Then we show that a unique local solutionof the mathematically well-posed problem exists. An adapted numericalscheme is built. Attention is devoted to the study of the nonlinearbehaviour of a pendulum via a numerical scheme with small constant timesteps. We describe the global behaviour of the free and forcedoscillations of the pendulum due to friction. We show that chaoticbehaviour occurs when friction is not too large. Lyapunov exponents arecomputed and a Melnikov relation is obtained as a limit of regularisedCoulomb friction. For larger friction, asymptotic behaviour correspondsto equilibrium.     

Numerical Simulation of Plasma Dynamics in a Nonuniform Magnetic Field

An efficient algorithm is proposed enabling numerical simulations of plasma dynamics in a nonuniform magnetic field. The present numerical data are in good agreement with experimental data obtained in a GOL-3 setup and with previous simulations. The experimentally observed effect of fast transfer of energy to ions is confirmed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 1, pp. 35–45, January–February, 2006.     

飞机大迎角俯仰运动下的横侧向气动特性及被动控制

在南航3m低速风洞内,利用一套两自由度动态试验机构,通过测力实验研究了某飞机模型静态和俯仰动态过程中大迎角下的横侧向气动特性,分析比较了在模型头部加上扰动片后,对横侧向气动特性产生的影响.研究结果表明,模型在静态大迎角下会产生较大的侧向力和偏航力矩,而模型的快速上仰过程则进一步加剧了模型头部流动的非对称性,在大迎角下产生较大的偏航力矩迟滞环;当在模型头部加扰动片后,不论是静态过程还是动态过程,都使得模型的侧向力和偏航力矩减小,从而改善了俯仰运动过程中大迎角下的横侧向气动特性.