一个用于界面不稳定性实验研究的竖式激波管

 研制了一个可用于气-气界面Rayleigh-Taylor（RT）不稳定性、气-液界面Richtmyer-Meshkov（RM）不稳定性、气-液界面RT和RM耦合不稳定性实验的竖式激波管装置,配备了相应的光学电学测量系统,给出了装置的主要技术参数,以及典型实验的压力-时间曲线。利用这套实验装置,可以精确地控制不稳定性的初始扰动界面,观测记录界面演化发展全过程。     

斜界面Rayleigh-Taylor不稳定性混合实验研究

 实验研究了不相溶流体斜界面Rayleigh-Taylor（R-T）不稳定性湍流混合区的混合不对称性特征。利用高压气体加速装有不同液体的箱体,加速度方向由轻液体指向重液体,此时界面是R-T不稳定性的。利用阴影测试技术,研究了初始倾角9°的ZnCl₂溶液/正己烷斜界面的演化规律,得出混合区内混合不对称、斜界面倾角渐增的规律。     

Suppression of Rayleigh-Taylor instability by gyroviscosity and sheared axial flow in imploding plasma pinches

The synergistic stabilizing effect of gyroviscosity and sheared axial flow on the Rayleigh-Taylor instability in Z-pinch implosions is studied by means of the incompressible viscid magneto-hydrodynamic equations.The gyroviscosity(or finite Larmor radius) effects are introduced in the momentum equation through an anisotropic ion stress tensor.Dispersion relation with the effect of a density discontinuity is derived.The results indicate that the short-wavelength modes of the Rayleigh-Taylor instability are easily stabilized by the gyroviscosity effects.The long wavelength modes are stabilized by the sufficient sheared axial flow.However,the synergistic effects of the finite Larmor radius and sheared axial flow can heavily mitigate the Rayleigh-Taylor instability.This synergistic effect can compress the Rayleigh-Taylor instability to a narrow wave number region.Even with a sufficient gyroviscosity and large enough flow velocity,the synergistic effect can completely suppressed the Rayleigh-Taylor instability in whole wave number region.     

轴向剪切流对Z箍缩等离子体瑞利-泰勒不稳定性的抑制

利用理想磁流体力学模型对有轴向剪切流的Z箍缩等离子体不稳定性进行了分析。给出了可压缩模型的色散关系,分别对可压缩及不可压缩模型中轴向剪切流对Z箍缩等离子体瑞利-泰勒不稳定性的抑制作用进行了比较,讨论了可压缩性对含有轴向剪切流系统不稳定性的影响。结果表明,可压缩性能够减缓瑞利-泰勒P凯尔文-亥姆霍兹(RTPKH)模扰动的增长,因而使得轴向剪切流对系统不稳定性的抑制作用表现得更为突出。计算结果还说明,在RT不稳定性线性增长阶段,等离子体温度较低,使用可压缩模型能够更真实地描述系统的状态。     

SYNERGETIC MITIGATION OF THE RAYLEIGH—TAYLOR INSTABILITY BY SHEAR AXIAL FLOWS AND FINITE LARMOR RADIUS

1　IntroductionImplodingZ pinchesemployinggaspuffs,cy lindricalfoils,andwirearraysareverycompactandefficientdevicesforcouplingelectromagneticenergyfromapulseforminglineintoadenseplas macolumn.Morerecently,aluminumwirearrayimplosionshaveproducedthehighestX raypowerof40TWonthe20TWSaturnacceleratoratSandiaNationalLaboratory[1].Thesedeviceshavepoten tialapplicationstocontrolledfusionaswellastointensesourcesofkeVX raysfornuclearweaponseffects.Unfortunately,theimplosionsarehighlysusceptibletoR…     

微小卫星贮箱内两相流界面变形定位及控制管理

液化气推进技术是微小伴随卫星在轨飞行采用的一种新型推进技术。微小伴星空间调姿、变轨过程需精确的推力控制,因此必须了解卫星推进系统中推进剂的形位分布。本文理论分析了空间微重力环境下液化气推进剂气/液界面的形位分布及变化,并通过落塔实验验证了微重力环境下有效控制液体推进剂的管理方法。     

Rayleigh-Taylor instability in variable density swirling flows

Using linear instability theory and nonlinear dynamics, the Rayleigh-Taylor instability of variable density swirling flows is studied. It is found that the flow topology could be predicted, when the instability sets in, using a function χ dependent on density and axial and azimuthal velocities. It is shown that even when the inner axial-flow is heavier than the outer one (a favorable case for the development of the Rayleigh-Taylor instability thanks to the centrifugal force) the instability is not necessarily Rayleigh-Taylor-dominated. It is also shown that when the Rayleigh-Taylor instability develops, it is helical.     

液滴超声速气动破碎初期界面不稳定性

针对液滴破碎问题,获得并揭示两相界面演化特征机理.采用数值模拟方法,观察了超声速条件下的液滴气动破碎初期的界面不稳定性.基于数值模拟结果和线性稳定性理论,综合分析表明,Rayleigh-Taylor不稳定性和Kelvin-Helmholtz不稳定性均对源于驻点和外环之间中段附近处的主导扰动产生作用.保持其他流动特性不变,降低K-H不稳定性的影响,对数值模拟进行了专门改进,进一步验证了前述结论.      

Direct Numerical Simulation of the Rayleigh-Taylor Instability with the Spectral Element Method

A novel method is proposed to simulate Rayleigh-Taylor instabilities using a specially-developed unsteady three-dimensional high-order spectral element method code. The numerical model used consists of Navier-Stokes equations and a transport-diffusive equation. The code is first validated with the results of linear stability perturbation theory. Then several characteristics of the Rayleigh-Taylor instabilities are studied using this three-dimensional unsteady code, including instantaneous turbulent structures and statistical turbulent mixing heights under different initial wave numbers. These results indicate that turbulent structures of Rayleigh-Taylor instabilities are strongly dependent on the initial conditions. The results also suggest that a high-order numerical method should provide the capability of simulating small scale fluctuations of Rayleigh-Taylor instabilities of turbulent flows.     

Weak Nonlinearity of Ablative Rayleigh--Taylor Instability

The weakly nonlinear regime of single mode ablative Rayleigh-Taylor instability is studied, with consideration of preheat effect and the width of the ablation front. The Rayleigh-Taylor linear growth rate agrees well with the direct numerical simulation. For the density perturbation, the amplitude distribution of the fundamental mode has one peak value whereas those of the second and third harmonics have two and three peak values, respectively. Harmonics generation versus wave number is also given and it is close to the result of direct numerical simulation.     

Three-dimensional array structures associated with Richtmyer-Meshkov and Rayleigh-Taylor instability

A boundary separating adjacent gas or liquid media is frequently unstable. Richtmyer-Meshkov and Rayleigh-Taylor instability cause the growth of intricate structures on such boundaries. All the lattice symmetries [rectangular (pmm2), square (p4mm), hexagonal (p6mm), and triangular (p3m1) lattices] which are of interest in connection with the instability of the surface of a fluid are studied for the first time. They are obtained from initial disturbances consisting of one (planar case, two-dimensional flow), two (rectangular cells), or three (hexagons and triangles) harmonic waves. It is shown that the dynamic system undergoes a transition during development from an initial, weakly disturbed state to a limiting or asymptotic stationary state (stationary point). The stability of these points (stationary states) is investigated. It is shown that the stationary states are stable toward large-scale disturbances both in the case of Richtmyer-Meshkov instability and in the case of Rayleigh-Taylor instability. It is discovered that the symmetry increases as the system evolves in certain cases. In one example the initial Richtmyer-Meshkov or Rayleigh-Taylor disturbance is a sum of two waves perpendicular to one another with equal wave numbers, but unequal amplitudes: a ₁(t=0)≠a ₂(t=0). Then, during evolution, the flow has p2 symmetry (rotation relative to the vertical axis by 180°), which goes over to p4 symmetry (rotation by 90°) at t→∞, since the amplitudes equalize in the stationary state: a ₁(t=∞)=a ₂(t=∞). It is shown that the hexagonal and triangular arrays are complementary. Upon time inversion (t→−t), “rephasing” occurs, and the bubbles of a hexagonal array transform into jets of a triangular array and vice versa. Zh. éksp. Teor. Fiz. 116, 908–939 (September 1999)     

Suppression of the Rayleigh-Taylor instability due to self-radiation in a multiablation target

A scheme to suppress the Rayleigh-Taylor instability has been investigated for a direct-drive inertial fusion target. In a high-Z doped-plastic target, two ablation surfaces are formed separately-one driven by thermal radiation and the other driven by electron conduction. The growth of the Rayleigh-Taylor instability is significantly suppressed on the radiation-driven ablation surface inside the target due to the large ablation velocity and long density scale length. A significant reduction of the growth rate was observed in simulations and experiments using a brominated plastic target. A new direct-drive pellet was designed using this scheme.     

A two-phase diffuse-interface model for Hele-Shaw flows with large property contrasts

A novel two-phase diffuse-interface model is used to simulate flows inside a Hele-Shaw cell. The model assumes that the two phases coexist inside the diffuse interface, with different velocities and properties. A separate equation is used to calculate the slip velocity between the two phases inside the diffuse interface. It is shown that for one-dimensional flows parallel to the diffuse interface, the results are independent of the diffuse-interface width, regardless of the magnitude of the density and viscosity contrasts between the phases. This two-phase approach is coupled with a phase-field equation for calculating the interface motion. The model is applied to a buoyancy-driven two-phase flow involving a Rayleigh-Taylor instability and validated through a comparison with available sharp-interface results. The flows and interface topology changes are investigated for large density and viscosity contrasts between the phases. The convergence of the results with respect to the interface width is examined in detail. It is shown that the two-phase model converges better than a standard diffuse-interface model that assumes the presence of a single velocity inside the diffuse interface. Remaining interface width dependencies can be attributed to the capillary stress term in the momentum equation.     

瑞利-泰勒不稳定性中尺度效应的数值模拟研究

 瑞利-泰勒不稳定性是一种由于密度梯度引起的界面不稳定性，在惯性约束聚变中具有重要的意义。利用被动标量输运模型对包含不同尺度初始扰动的界面演化过程进行数值模拟。计算结果表明界面的初始形状对不稳定性的发展具有很大的影响，狭长型扰动比正方型扰动发展慢。另外，不同尺度扰动的相互作用一般会减小沿界面发展方向运动的动能，使能量更多地用于平行于界面方向的运动。     

R-T不稳定性发展的二维和三维对比研究

 应用摄动方法研究不同密度流体交界面上的Rayleigh-Taylor（R-T）不稳定性问题，分别研究了二维（2D）与三维（3D）R-T不稳定性的线性与非线性增长，并进行了对比。得出的结论为：（1） 三维不稳定性的线性增长比二维情况要快，且同维度上相同的初始扰动波长，其时间增长指数受另一维度上的初始波长影响。（2） R-T不稳定性发展在考虑了非线性增长的影响之后，二维与三维扰动增长的图像也不一样。     

Numerical simulation of the hydrodynamic instability experiments and flow mixing

Based on the numerical methods of volume of fluid (VOF) and piecewise parabolic method (PPM) and parallel circumstance of Message Passing Interface (MPI), a parallel multi-viscosity-fluid hydrodynamic code MVPPM (Multi-Viscosity-Fluid Piecewise Parabolic Method) is developed and performed to study the hydrodynamic instability and flow mixing. Firstly, the MVPPM code is verified and validated by simulating three instability cases: The first one is a Riemann problem of viscous flow on the shock tube; the second one is the hydrodynamic instability and mixing of gaseous flows under re-shocks; the third one is a half height experiment of interfacial instability, which is conducted on the AWE’s shock tube. By comparing the numerical results with experimental data, good agreement is achieved. Then the MVPPM code is applied to simulate the two cases of the interfacial instabilities of jelly models accelerated by explosion products of a gaseous explosive mixture (GEM), which are adopted in our experiments. The first is implosive dynamic interfacial instability of cylindrical symmetry and mixing. The evolving process of inner and outer interfaces, and the late distribution of mixing mass caused by Rayleigh-Taylor (RT) instability in the center of different radius are given. The second is jelly layer experiment which is initialized with one periodic perturbation with different amplitude and wave length. It reveals the complex processes of evolution of interface, and presents the displacement of front face of jelly layer, bubble head and top of spike relative to initial equilibrium position vs. time. The numerical results are in excellent agreement with that experimental images, and show that the amplitude of initial perturbations affects the evolvement of fluid mixing zone (FMZ) growth rate extremely, especially at late times.     

Rayleigh-Taylor不稳定性弱非线性阶段界面效应对谐波的影响

为了更好地理解不同空间坐标系下流体界面对Rayleigh-Taylor（RT）不稳定性弱非线性阶段谐波的影响,文章采用3阶小扰动展开法,解析研究了球坐标空间经典RT不稳定性弱非线性阶段谐波的演化规律,并和柱坐标空间以及直角坐标空间相应结果进行了对比研究.当球坐标系和直角坐标系中RT不稳定性界面扰动波长相同,球坐标系中初始扰动半径为无穷大时（即球坐标下RT不稳定性初始扰动半径相对于扰动波长为无穷大时）,球坐标下RT不稳定性前4次谐波的结果和直角坐标系下的相应结果相同.研究表明:由初始界面曲率引起的Bell-Plesset（BP）效应和空间效应（直角坐标空间、柱坐标空间和球坐标空间）对谐波发展有较大的影响.即在不同正交曲线坐标系下,不同曲率的流体界面效应对RT不稳定性谐波发展有较大的影响.对于柱坐标空间和球坐标空间,2阶对0次谐波的反馈加强了界面向内收缩.研究还表明:界面效应增加了2次谐波的负反馈,然而,对于基模和3次谐波却有不同的影响.      

金属锡Rayleigh-Taylor不稳定性对模型参数敏感性的数值分析

利用自研的爆轰与冲击动力学欧拉计算程序和Steinberg-Guinan(SG)本构模型,数值模拟分析了样品初始参数(初始振幅、初始波长、样品初始厚度)和SG本构模型初始参数对爆轰驱动锡Rayleigh-Taylor(RT)不稳定性增长的影响。结果表明金属锡样品的初始参数对其RT不稳定性增长有很大的影响。RT不稳定性增长随着初始振幅的减小而减小,且存在一个截止初始振幅;存在一个最不稳定的模态(波长),当初始波长大于该波长时,RT不稳定性增长随着初始波长的减小而增大,反之,RT不稳定性增长随着初始波长的减小而减小;样品厚度的增大可以抑制RT不稳定性增长,而且存在一个样品截止厚度。金属锡的RT不稳定性增长对其SG本构模型应变硬化系数和应变硬化指数的变化不敏感,而对压力硬化系数和热软化系数比较敏感。从采用扰动增长法预估材料强度的角度来说,修正压力硬化系数以获得锡合理的材料强度是合理的途径。     

Effects of compressibility on the Rayleigh-Taylor instability in Z-pinch implosions

1 IntroductionZ-pinch imploding plasma thin shells formed bygas puffs,cylindrical foils,and wire arrays can radiatehundreds of kilojoules of soft X-rays.They are highefficient and compact devices which have wideapplied prospects,specially,as intense X-ray…     

界面张力对Rayleigh-Taylor不稳定性的影响

将具有简单速度势的Layzer模型和Zufiria模型推广至非理想流体情况, 并分别利用这两种模型研究了界面张力对Rayleigh-Taylor不稳定性的影响. 首先得到了两种模型下气泡的渐近速度和渐近曲率的解析表达式; 其次系统研究了界面张力对气泡的渐近速度和渐近曲率的影响; 最后将两种模型进行了比较, 并将气泡的渐近速度和数值模拟进行了比较. 研究表明: 界面张力压低了气泡的速度, 但对曲率没有影响; 利用简单速度势的Layzer模型所得的气泡的渐近速度比复杂速度势的Layzer模型的值小, 但是比Zufiria模型的值大; 当阿特伍德数等于1时, 简单速度势的Layzer模型和复杂速度势的Layzer模型给出的结果一致. 关键词： Rayleigh-Taylor不稳定性 界面张力 Layzer模型 Zufiria模型