再冲击载荷作用下流动混合的数值模拟

基于多介质的流体体积分数VOF(volume of fluid)方法和PPM(piecewise parabolic method)方法,给出和发展了可用于多介质粘性流体动力学的数值计算方法和计算代码MVPPM(multi-viscosity-fluidpiecewise parabolic method)。为了检验和验证此计算代码,对某激波管实验——再冲击载荷作用下的流体动力学不稳定性及其导致的混合过程进行了数值模拟,计算结果与实验结果一致。同时还研究了激波反射冲击作用下流体混合区的演化情况,在反射激波和混合区相互作用的瞬间,混合区的宽度明显减小,之后又迅速增大;另外,混合增长率与初始扰动的频谱有很大关联。通过对有粘性(分子动力学粘性)和无粘性结果的对比,发现粘性对混合区的影响很小。

单模态Richtmyer-Meshkov不稳定性数值模拟

采用VOF(Volume of Fluid)方法和PPM(Piecewise Parabolic Method)方法,发展了可用于可压缩多介质粘性流体动力学问题的数值模拟方法MVPPM(Multi-Viscous-Fluid Piecewise Parabolic Method)。利用MVPPM对多个具有不同初始扰动振幅的二维和三维单模态RM(Richtmyer-Meshkov)不稳定性模型进行了数值计算,并与理论模型的计算结果进行了比较。结果表明,无论二维还是三维情况,当初始扰动振幅相对于波长较小的时候,计算的扰动振幅和增长率与理论模型的计算结果一致。当初始扰动波长不变而振幅逐渐增大时,界面振幅和增长率也逐渐增大。对于具有相同初始扰动的情况,三维计算结果在线性段与二维计算结果相同,但是在非线性段比二维结果大,说明非线性和三维效应在RM不稳定性发展过程中起着重要作用。

激波管实验和果冻实验界面不稳定性数值计算

运用算子分裂技术,增加二阶空间中心差方法和两步Rung-Kutta时间推进方法计算动力学粘性以及热流部分对流场的影响,将可压缩多介质流体动力学高精度欧拉计算方法MFPPM发展到适用于NS(Navier Stokes)方程的可压缩多介质粘性流体计算方法MVPPM。通过与界面不稳定性实验结果的比较,来检验计算方法的正确性,并验证计算程序的有效性。主要包括一个激波管实验和两个果冻实验,即英国AWE(Atomic Weapons Establishment)激波管实验和LLNL实验室的果冻环实验以及中国工程物理研究院流体物理研究所冲击波物理与爆轰物理实验室进行的爆轰驱动下单层果冻界面不稳定性实验。研究结果表明:数值模拟结果与实验测量结果以及对应时刻的实验图象均吻合较好。

Experimental and numerical investigation of inclined air/SF_6 interface instability under shock wave

The shock tube experiments of inclined air/SF6 interface instability under the shock wave with the Mach numbers 1.23 and 1.41 are conducted. The numerical simulation is done with the parallel algorithm and the multi-viscous-fluid and turbulence （MVFT） code of the large-eddy simulation （LES）. The developing process of the interface accelerated by the shock wave is reproduced by the simulations. The complex wave structures, e.g., the propagation, refraction, and reflection of the shock wave, are clearly revealed in the flows. The simulated evolving images of the interface are consistent with the experimental ones. The simulated width of the turbulent mixing zone （TMZ） and the displacements of the bubble and the spike also agree well with the experimental data. Also, the reliability and effectiveness of the MVFT in simulating the problem of interface instability are validated. The more energies are injected into the TMZ when the shock wave has a larger Mach number. Therefore, the perturbed interface develops faster.

Evolution of Hydrodynamic Instability on Planar Jelly Interface Driven by Explosion

A high precision numerical algorithm MVPPM (multi-viscous-fluid piecewise parabolic method) is proposed and applied to study the multi-viscous-fluid dynamics problems. Three planar jelly experiments with periodic cosine perturbation on the initial interface have been conducted and numerically simulated by MVPPM. Good agreement between experimental and numerical results has been achieved, including the shape of jelly interface, the displacements of front face of jelly layer, bubble top and spike head. The effects of initial conditions (including amplitude and wave length of perturbation, thickness of jelly layer, etc.) on the evolution of the jelly interface have been numerically analyzed. It is found that the key affecting factors are the perturbation amplitude and thickness of jelly layer. The hydrodynamic instability on double planar jelly layers driven by explosion has been investigated numerically to examine their laws of evolution, and an interesting phenomenon is observed.

单模态RM不稳定性的二维和三维数值计算研究(英文)

发展了可用于可压缩多介质粘性流体动力学问题的数值模拟方法和代码MVPPM(Multi-Viscous-Fluid Piecewise Parabolic Method)。利用MVPPM方法对多个具有不同初始扰动振幅和波长的二维和三维单模态Richtmyer-Meshkov(RM)不稳定性模型进行了数值计算,并和理论模型的计算结果进行了比较。结果显示,扰动界面的发展与扰动的初始条件密切相关。无论二维还是三维情况,当初始扰动强度较小的时候,数值计算的扰动振幅及增长率和理论模型的计算结果一致。对于具有相同初始扰动的情况,三维数值计算结果在线性段与二维计算结果相同,但是在非线性段比二维结果大,说明非线性和三维效应在RM不稳定性发展过程中起着重要作用。