激波绕拐角运动流场的数值模拟

本文采用时间和空间方向均为二阶精度的ＴＶＤ有限差分格式来数值模拟冲击波绕三角形、半菱形障碍物的非定常复杂流场，算出了与理论分析一致的流场结构，成功地捕获了激波，膨胀波，二次弱反射激波等物理特征。     

运动激波绕射的数值模拟

用高分辨率的ＴＶＤ格式对复杂的激波反射过程做数值模拟。激波斜壁反射计算结果与实验干涉照片、压力分布、密度分布均符合得很好。通过对运动激波绕射凸起物的数值模拟,得到非定常反射流场发展过程及复杂流场间断的清晰图谱。     

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

激波绕射方台的非定常流场研究

本文以数值方法研究了运动激波绕射矩形突台的非定常流场发展过程。计算中采用ＴＶＤ格式，得到了清晰的流场变化图谱。从流场图中可分辨出旋涡的形式过程及分离域形态。研究结果表明，本方法可以很好地描述激波绕射方台形成的复杂流动过程，对流场间断有较高的捕获力。     

粘性高超声速激波—激波碰撞的计算和分析

本文用特片变量型的隐式ＮＮＤ格式计算了五种类型的高超声速激波－激波干扰流场，得到了清晰的流场结构图画。激波和剪切层等分辨得很清楚。计算的物面压力与实验甚为吻合。借助于激波极曲线，对第４类激波干扰结构进行了分析，解释了干扰引起物面压力增高的原因。     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

压力恢复系统扩压器激波串现象的数值模拟

通过求解由S-A一方程湍流模式封闭的三维雷诺平均Navier-Stokes方程,采用有限体积方法对激光器压力恢复系统变截面扩压器常温状态下的三维流场进行了数值模拟.数值计算和实验符合良好,较好模拟了扩压器中由激波/边界层干扰诱导的复杂流场的流场特性并再现了流场中的"激波串"和"伪激波"现象.本文对激波串的形成机理进行了分析,发现了几种变截面扩压器特有的激波串结构,并研究了扩压器扩张角对扩压效率的影响.     

可压缩球涡和平面激波相互作用的数值研究

行进了可压缩球涡和平面激波相互作用的数值研究。一方面，用非定常Ｅｕｌｅｒ方程数值模拟了一种可压缩球涡的自身演化，同时，利用Ｒａｎｋｉｎｅ－Ｈｕｇｏｎｉｏｔ关系，在流场中嵌入运动激波，求解了激波－球涡相干的流动过程。     

Shock Wave Interaction with the Magnetopause

The magnetopause is in continuous motion and shock waves and impulsive acceleration events can occur. As an example, we show that the interaction of an interplanetary shock with the bow shock can generate a shock wave that after passing through the magnetosheath can interact with the magnetopause. In fluid dynamics, when a shock wave encounters a fluid discontinuity, the interface may become unstable and form bubbles and spikes. We consider this Richtmyer-Meshkov instability in magnetohydrodynamics. At the dayside magnetopause, the instability tends to be stabilized by the magnetic field. However, the shock wave interaction can initiate magnetic field reconnection for the southward IMF, which may be important in strong interplanetary shock events. This revised version was published online in August 2006 with corrections to the Cover Date.     

可压缩涡环和平面激波相互作用的数值研究

本文从非定常的Ｅｕｌｅｒ方程出发,进行了可压缩涡环和平面激波相互作用的数值研究。     

Initial ISEE magnetometer results: Shock observation

ISEE-1 and -2 magnetic field profiles across 6 terrestrial bow shock and one interplanetary shock are examined. The interplanetary shock illustrates the behavior of a low Mach number shock. It had an upstream whistler wave precursor with an apparent wavelength of 180 km. The shock thickness was about 90 km for the thickness of the final field jump or 270 km for the exponential growth of the precursor wave packet. The ion inertial length was 50 km, upstream of the shock.Three examples of low or moderate , high Mach number, quasiperpendicular shocks are examined. These did not have upstream waves, but rather had waves growing in the field gradient. The growth length for these waves and the shock profile was of the order of the ion inertial length.Two examples of high shocks showed little coherence in field variation even though the two vehicles were only a few hundred kilometers apart. Thus, we cannot gauge their velocity and turn the time profiles into distances. The final crossing examined shows clearly the effect of changing the orientation of the interplanetary magnetic field. Initially the upstream magnetic field made an angle of about 80° to the shock normal and the shock position remained fairly steady. Then the field rotated to 45° to the normal and the field profiles became very irregular and the shock position very unstable. Discrete wave packets appeared.Finally, we present the joint behavior of wave, particle and field data across some of these shocks to show some of the myriad of shock features whose behavior we are now beginning to investigate.     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

受限空间内激波与火焰作用的三维计算

激波与火焰的相互作用常发生在超声速燃烧和燃烧转爆轰过程中,为深入了解这一现象,采用带化学反应的三维Navier-Stokes方程,对平面入射激波及其反射激波与火焰的作用进行了计算研究,其中燃烧过程使用单步反应模型描述。研究结果显示:在受限空间内模拟激波与火焰作用,能更好地符合实验结果,从而体现出受限空间的三维效应;火焰在入射激波的作用下主要经Richtmyer-Meshkov不稳定而发生变形,此时火焰变形以物理作用为主,燃烧膨胀效果相对不明显;当反射激波与变形火焰再次作用后,火焰迅速膨胀变形,放热率维持在较高水平,此时化学反应过程起主要作用;在反射激波的作用下,变形火焰复杂三维涡结构的形成能强化热量与质量的输送,提高燃烧速率。     

Shock wave configurations and reflection hysteresis outside a planar Laval nozzle

When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection(RR) pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection(MR) pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited condition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these conditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RRMMR transition in the dual-solution domain.     

动态微涡流发生器对激波/边界层干扰的影响研究

本文基于OPENFOAM数值仿真平台,采用动态网格技术和湍流离散涡(DES)模型,研究了微涡流发生器以一定速度向下游移动时,激波/边界层干扰(SWBLI)流场特性的变化,重点关注干扰区域内的流向和展向的流场特性。来流马赫数为4,微涡流发生器向下游移动速度为0m/s,20m/s和40m/s。研究表明：当MVG向下游移动时,SWBLI区域的“弓”形高压区会演化成“双弓”形；入射激波形成高压区的压力明显降低,同时,入射激波和反射激波形成高压区的峰值位置均会向下游移动；流场下游 “双圆弧”状高压区的高度逐渐降低；SWBLI区域边界层的高度逐渐降低,同时边界层底部的速度也有所降低；随着MVG移动速度的增加,对SWBLI流场的控制效果更加明显；动态MVG对流场的控制是通过尾迹涡和波系结构实现的。     

关于建立高阶差分格式的问题

为了能在不太密的网格上捕捉到流场的细致结构。通常采用高阶精度的差分 式进行数值模拟。为能抑制计算中和在激波附近产生的虚假波动,本文从物理构思出发,提出了建立高阶格式的两个基本原则,作者称之为抑制波动的原则和稳定性原则。