差分干涉CT技术诊断激波管界面不稳定性

对差分干涉计算机层析技术诊断激波加载下气柱界面不稳定性的可行性进行初步理论和实验研究。针对有限投影数的联合代数重建算法,提出了一种新的投影矩阵计算方法,采用特定投影方向,各投影方向具有不同采样间隔。投影方向的对称性和投影线的规律性使计算过程大为减化。用该重建算法对气柱密度场进行数值模拟,结果表明4方向投影能够重建激波作用下气柱演化初期的体密度场,8方向投影基本能够重建演化中期的体密度场。用差分干涉系统对激波管气柱进行3方向诊断的实验结果表明,差分干涉计算机层析技术用于激波冲击下的气柱不稳定性密度场诊断是可行的。     

差分干涉CT技术诊断激波管界面不稳定性

对差分干涉计算机层析技术诊断激波加载下气柱界面不稳定性的可行性进行初步理论和实验研究。针对有限投影数的联合代数重建算法,提出了一种新的投影矩阵计算方法,采用特定投影方向,各投影方向具有不同采样间隔。投影方向的对称性和投影线的规律性使计算过程大为减化。用该重建算法对气柱密度场进行数值模拟,结果表明4方向投影能够重建激波作用下气柱演化初期的体密度场,8方向投影基本能够重建演化中期的体密度场。用差分干涉系统对激波管气柱进行3方向诊断的实验结果表明,差分干涉计算机层析技术用于激波冲击下的气柱不稳定性密度场诊断是可行的。     

Large-eddy simulations of the Richtmyer-Meshkov instability of rectangular interfaces accelerated by shock waves

A parallel algorithm and code MVFT(multi-viscous-fluid and turbulence) of large-eddy simulation(LES) is developed from our MVPPM(multi-viscous-fluid piecewise parabolic method),and performed to solve the multi compressible Navier-Stokes(N-S) equations.The effect of the unresolved subgrid-scale(SGS) motions on the large scales is represented by different SGS stress models in LES.A Richtmyer-Meshkov instability experiment of the evolution of a rectangular block of SF6,which occupies half of the height of the ...     

Molecular dynamics simulation of cylindrical Richtmyer-Meshkov instability

The microscopic-scale Richtmyer-Meshkov(RM) instability of a single-mode Cu-He interface subjected to a cylindrically converging shock is studied through the classical molecular dynamics simulation. An unperturbed interface is first considered to examine the flow features in the convergent geometry, and notable distortions at the circular inhomogeneity are observed due to the atomic fluctuation. Detailed processes of the shock propagation and interface deformation for the single-mode interface impacted by a converging shock are clearly captured. Different from the macroscopic-scale situation, the intense molecular thermal motions in the present microscale flow introduce massive small wavelength perturbations at the single-mode interface, which later significantly impede the formation of the roll-up structure. Influences of the initial conditions including the initial amplitude,wave number and density ratio on the instability growth are carefully analyzed. It is found that the late-stage instability development for interfaces with a large perturbation does not depend on its initial amplitude any more. Surprisingly, as the wave number increases from 8 to 12, the growth rate after the reshock drops gradually. The distinct behaviors induced by the amplitude and wave number increments indicate that the present microscopic RM instability cannot be simply characterized by the amplitude over wavelength ratio(η). The pressure history at the convergence center shows that the first pressure peak caused by the shock focusing is insensitive to η, while the second one depends heavily on it.     

激波与SF_6梯形气柱相互作用的数值模拟

基于大涡模拟,结合五阶加权基本无振荡格式以及沉浸边界法对平面入射激波与两种SF_6梯形重气柱的相互作用过程进行了数值模拟,数值结果清晰地显示了激波诱导Richtmyer-Meshkov不稳定性所导致的两种梯形重气柱的变形过程,详细分析了入射激波在两种梯形重气柱界面发生反射、折射、绕射以及折射激波与透射激波在气柱内部来回反射的过程,并研究了该过程中所产生的复杂波系结构,对两种梯形气柱变形过程中与周围空气的混合过程进行了分析;通过记录气柱界面四个特征尺寸随时间的变化对两种梯形气柱界面的不同演化过程进行了定量分析。     

Direct visualization of laser-driven focusing shock waves

Direct real-time visualization and measurement of laser-driven shock generation, propagation, and 2D focusing in a sample are demonstrated. A substantial increase of the pressure at the convergence of the cylindrical acoustic shock front is observed experimentally and simulated numerically. Single-shot acquisitions using a streak camera reveal that at the convergence of the shock wave in water the supersonic speed reaches Mach 6, corresponding to the multiple gigapascal pressure range ～30 GPa.     

Plane and cylindrical shock waves in magnetogasdynamics

Summary The Chisnell-Chester-Whitham method has been used to study the propagation of diverging plane and cylindrical shock waves through an ideal gas in the presence of a magnetic field having only constant axial and azimuthal components, simultaneously for both weak and strong cases. By assuming an initial density distribution ϱ₀ = ϱ′r ^−W , where ϱ′ is the density at the plane/axis of symmetry andw is a constant, the analytical expressions for shock velocity and shock strength have been obtained. The expressions for the pressure, the density and the particle velocity immediately behind the shock have been derived for both the cases.

---

Riassunto Il metodo di Chisnell-Chester-Whitham è stato usato per studiare la propagazione di onde d’urto divergenti cilindriche e piane attraverso un gas ideale in presenza di un campo magnetico che ha solo componenti costanti assiali e azimutali, contemporaneamente per casi deboli e forti. Assumendo una distribuzione di densità iniziale ϱ₀ = ϱ′r ^−W , dove ϱ′ è la densità al piano/asse di simmetria ew è una costante, si sono ottenute le espressioni analitiche per la velocità d’urto e la forza d’urto. Le espressioni per la pressione, la densità e la velocità delle particelle immediatamente dopo l’urto sono state determinate per entrambi i casi.

     

High-amplitude single-mode perturbation evolution at the Richtmyer-Meshkov instability

The single-mode Richtmyer-Meshkov hydrodynamic instability at light/heavy (air/SF6 and air/CO2), close density (air/N2), and heavy/light (air/He) interfaces has been experimentally studied for a low incident shock wave Mach number. Two identical 2D half sinusoidal initial perturbations, with a relatively high initial amplitude, were considered in order to rapidly reach the nonlinear regime and check the reduction of the initial growth rate compared to that predicted by the small-amplitude theory. The growth rate measurements for the air/SF6 and air/CO2 cases are in excellent agreement with the nonlinear model of Sadot et al. coupled with a reduction factor suggested by Rikanati et al. In the air/N2 case, the reversal phase can be precisely described by the linear theory. Finally, the heavy/light experiment is well described by the Vandenboomgaerde model also coupled with a smaller reduction factor.     

Investigation of the Richtmyer-Meshkov instability with stereolithographed interfaces

A novel method to set highly accurate initial conditions has been designed in the context of shock tube experiments for the Richtmyer-Meshkov instability study. Stereolithography has been used to design the membrane supports which initially materialize the gaseous interface. The visualizations of both heavy-light and light-heavy sinusoidal interfaces were carried out with laser sheet diagnostics. Experiments are in very good agreement with theory and simulations for the heavy-light case, but probably due to the membrane effects, quickly deviate from them in the light-heavy configuration.     

激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

Surface tension and instability of shock waves in gases

The pressure anisotropy within a shock layer results in an additional property of the shock wave: the surface tension. Its value amounts to 10³ dyn/cm for Mach number 10. The surface tension term in the characteristic equation causes an absolute shock wave instability in the spontaneous sound emission regime.     

The D'yakov-kontorovich instability of shock waves in real gases

In the 1950s, D'yakov and Kontorovich predicted that under certain conditions perturbed shock waves in nonideal gases can become unstable by emitting undamped sound and entropy-vortex waves. For the last 45 years, though, little progress has been made in the identification and numerical modeling of physical conditions for which this phenomenon might occur. Using a van der Waals equation of state, we present for the first time a dynamical simulation of a D'yakov-Kontorovich instability. The two-dimensional emission pattern of acoustic waves appearing in the simulation agrees with the prediction of a linearized theory.     

Focusing of shock waves induced by optical breakdown in water

The focusing of laser-generated shock waves by a truncated ellipsoidal reflector was experimentally and numerically investigated. Pressure waveform and distribution around the first (F(1)) and second foci (F(2)) of the ellipsoidal reflector were measured. A neodymium doped yttrium aluminum garnet laser of 1046 nm wavelength and 5 ns pulse duration was used to create an optical breakdown at F(1), which generates a spherically diverging shock wave with a peak pressure of 2.1-5.9 MPa at 1.1 mm stand-off distance and a pulse width at half maximum of 36-65 ns. Upon reflection, a converging shock wave is produced which, upon arriving at F(2), has a leading compressive wave with a peak pressure of 26 MPa and a zero-crossing pulse duration of 0.1 mus, followed by a trailing tensile wave of -3.3 MPa peak pressure and 0.2 mus pulse duration. The -6 dB beam size of the focused shock wave field is 1.6 x 0.2 mm(2) along and transverse to the shock wave propagation direction. Formation of elongated plasmas at high laser energy levels limits the increase in the peak pressure at F(2). General features in the waveform profile of the converging shock wave are in qualitative agreement with numerical simulations based on the Hamilton model.     

Study for particular solutions of cylindrical shock waves in magnetogasdynamics

In this research article, a problem of shock waves propagation in magnetogasdynamics is considered in ideal gas. The problem can be represented by hyperbolic non-linear system of partial differential equations (PDEs). It is solved by similarity method under invariant surface conditions. I reduced non-linear system of magnetogasdynamics into first order system of ordinary differential equations (ODEs) and it is non-linear. The ratio of heat capacity at constant pressure to heat capacity at constant volume is also known as the adiabatic constant and their values for monatomic, diatomic gasses are more than one in ideal gas. So, I have obtained exact solution of magnetogasdynamics system in ideal gas for particular value of adiabatic index$\gamma =2$. The effect of the ambient density exponent $\theta$ on flow quantities density, velocity, pressure, and magnetic field are determined before the shock.