激光烧蚀RT不稳定性线性增长率和非线性行为的数值研究

 给出了不同情况下多种波长的烧蚀瑞利 泰勒不稳定性线性增长率的二维计算结果，并与Takabe公式和Sanz公式进行了比较，最后给出了单模激光烧蚀RTI非线性发展行为的数值结果。线性增长率的二维计算结果很好验证了Sanz公式，表明β值应是2，不是3。当烧蚀速度较大或烧蚀面较宽时，用Takabe公式估计的烧蚀RTI线性增长率与二维计算值明显不符。     

Nonlinear theory of the ablative Rayleigh-Taylor instability

A fully nonlinear sharp-boundary model of the ablative Rayleigh-Taylor instability is derived and closed in a similar way to the self-consistent closure of the linear theory. It contains the stabilizing effect of ablation and accurately reproduces the results of 2D DRACO simulations. The single-mode saturation amplitude, bubble and spike evolutions in the nonlinear regimes, and the seeding of long-wavelength modes via mode coupling are determined and compared with the classical theory without ablation. Nonlinear stability above the linear cutoff is also predicted.     

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.     

Lattice Boltzmann study of three-dimensional immiscible Rayleigh–Taylor instability in turbulent mixing stage

In this paper, we numerically studied the late-time evolutional mechanism of three-dimensional (3D) single-mode immiscible Rayleigh–Taylor instability (RTI) by using an improved lattice Boltzmann multiphase method implemented on graphics processing units. The influences of extensive dimensionless Reynolds numbers and Atwood numbers on phase interfacial dynamics, spike and bubble growth were investigated in details. The longtime numerical experiments indicate that the development of 3D singlemode RTI with a high Reynolds number can be summarized into four different stages: linear growth stage, saturated velocity growth stage, reacceleration stage and turbulent mixing stage. A series of complex interfacial structures with large topological changes can be observed at the turbulent mixing stage, which always preserve the symmetries with respect to the middle axis for a low Atwood number, and the lines of symmetry within spike and bubble are broken as the Atwood number is increased. Five statistical methods for computing the spike and bubble growth rates were then analyzed to reveal the growth law of 3D single-mode RTI in turbulent mixing stage. It is found that the spike late-time growth rate shows an overall increase with the Atwood number, while the bubble growth rate experiences a slight decrease with the Atwood number at first and then basically maintains a steady value of around 0.1. When the Reynolds number decreases, the later stages cannot be reached gradually and the evolution of phase interface presents a laminar flow state.     

Analytical model of nonlinear, single-mode, classical Rayleigh-Taylor instability at arbitrary Atwood numbers

An analytical model of the nonlinear bubble evolution of single-mode, classical Rayleigh-Taylor instability at arbitrary Atwood numbers (A(T)) is presented. The model is based on an extension of Layzer's theory [Astrophys. J. 122, 1 (1955)] previously applied only to the fluid-vacuum interfaces (A(T) = 1). The model provides a continuous bubble evolution from the earlier exponential growth to the nonlinear regime when the bubble velocity saturates at U(b) = square root of [2A(T)/(1+A(T)) (g/C(g)k)], where k is the perturbation wave number, g is the interface acceleration, and C(g) = 3 and C(g) = 1 for the two-dimensional and three-dimensional geometries, respectively.     

等离子体中的消融激波及其时间发展

本文研究了等离子体消融所造成的激波，找出了位于激波波前及消融波前的各流体动力学参量间的相互依赖关系，对消融激波的生成与时间发展过程进行了描述。     

Theoretical and simulation research of hydrodynamic instabilities in inertial-confinement fusion implosions

Inertial fusion energy(IFE)has been considered a promising,nearly inexhaustible source of sustainable carbon-free power for the world's energy future.It has long been recognized that the control of hydrodynamic instabilities is of critical importance for ignition and high-gain in the inertial-confinement fusion(ICF)hot-spot ignition scheme.In this mini-review,we summarize the progress of theoretical and simulation research of hydrodynamic instabilities in the ICF central hot-spot implosion in our group over the past decade.In order to obtain sufficient understanding of the growth of hydrodynamic instabilities in ICF,we first decompose the problem into different stages according to the implosion physics processes.The decomposed essential physics processes that are associated with ICF implosions,such as Rayleigh-Taylor instability(RTI),Richtmyer-Meshkov instability(RMI),Kelvin-Helmholtz instability(KHI),convergent geometry effects,as well as perturbation feed-through are reviewed.Analytical models in planar,cylindrical,and spherical geometries have been established to study different physical aspects,including density-gradient,interface-coupling,geometry,and convergent effects.The influence of ablation in the presence of preheating on the RTI has been extensively studied by numerical simulations.The KHI considering the ablation effect has been discussed in detail for the first time.A series of single-mode ablative RTI experiments has been performed on the Shenguang-Ⅱlaser facility.The theoretical and simulation research provides us the physical insights of linear and weakly nonlinear growths,and nonlinear evolutions of the hydrodynamic instabilities in ICF implosions,which has directly supported the research of ICF ignition target design.The ICF hot-spot ignition implosion design that uses several controlling features,based on our current understanding of hydrodynamic instabilities,to address shell implosion stability,has been briefly described,several of which are novel.     

Effect of initial phase on the Rayleigh—Taylor instability of a finite-thickness fluid shell

Rayleigh—Taylor instability (RTI) of finite-thickness shell plays an important role in deep understanding the characteristics of shell deformation and material mixing. The RTI of a finite-thickness fluid layer is studied analytically considering an arbitrary perturbation phase difference on the two interfaces of the shell. The third-order weakly nonlinear (WN) solutions for RTI are derived. It is found the main feature (bubble-spike structure) of the interface is not affected by phase difference. However, the positions of bubble and spike are sensitive to the initial phase difference, especially for a thin shell (kd<1), which will be detrimental to the integrity of the shell. Furthermore, the larger phase difference results in much more serious RTI growth, significant shell deformation can be obtained in the WN stage for perturbations with large phase difference. Therefore, it should be considered in applications where the interface coupling and perturbation phase effects are important, such as inertial confinement fusion.     

Pulse chirping effect on controlling the transverse cavity oscillations in nonlinear bubble regime

The propagation of an intense laser pulse in an under-dense plasma induces a plasma wake that is suitable for the acceleration of electrons to relativistic energies. For an ultra-intense laser pulse which has a longitudinal size shorter than the plasma wavelength, λp, instead of a periodic plasma wave, a cavity free from cold plasma electrons, called a bubble, is formed behind the laser pulse. An intense charge separation electric field inside the moving bubble can capture the electrons at the base of the bubble and accelerate them with a narrow energy spread. In the nonlinear bubble regime, due to localized depletion at the front of the pulse during its propagation through the plasma, the phase shift between carrier waves and pulse envelope plays an important role in plasma response. The carrier–envelope phase(CEP) breaks down the symmetric transverse ponderomotive force of the laser pulse that makes the bubble structure unstable. Our studies using a series of two-dimensional(2D) particle-in-cell(PIC) simulations show that the frequency-chirped laser pulses are more effective in controlling the pulse depletion rate and consequently the effect of the CEP in the bubble regime. The results indicate that the utilization of a positively chirped laser pulse leads to an increase in rate of erosion of the leading edge of the pulse that rapidly results in the formation of a steep intensity gradient at the front of the pulse. A more unstable bubble structure, the self-injections in different positions, and high dark current are the results of using a positively chirped laser pulse. For a negatively chirped laser pulse, the pulse depletion process is compensated during the propagation of the pulse in plasma in such a way that results in a more stable bubble shape and therefore, a localized electron bunch is produced during the acceleration process. As a result, by the proper choice of chirping, one can tune the number of self-injected electrons, the size of accelerated bunch and its energy spectrum to the values required for practical applications.     

用高精度有限差分法模拟烧蚀瑞利-泰勒不稳定性

求解烧蚀面附近流场的定常解,并以此作为基本流实现用高精度的WENO格式对烧蚀瑞利-泰勒不稳定性的数值模拟.线性增长率与Lindl公式以及线性稳定性分析给出的结果相符合,证明该数值模拟方法的准确性与精度,该方法还具有较好的界面变形捕捉能力.     

Laser-induced bubble formation and collapse in a nonlinear Rayleigh-Taylor unstable interface in a thin layer of molten metal

Bubble formation and bubble collapse zones in a laser-induced nonlinear Rayleigh-Taylor (R-T) instability were shown to follow the contours of a series of R-T shock-fronts generated on an irregular planar (terrace-like) target surface. Bubble formation occurs in the regime of target planar vaporization. Spatial variation of the bubble density distribution, the bubble size, and the bubble-bubble distance, as a function of distance from the shock front envelope, were determined. Bubble collapse occurs in the regime of planar-to-volume boiling transition and proceeds by the so-called “chain reaction” collapse mechanism inside a 2D bubble array. The contours of bubble generation and of bubble collapse were simulated by using the analytical model of Ott with variable phase and variable amplitude of R-T modes.     

任意Atwood数Rayleigh-Taylor和 Richtmyer-Meshkov 不稳定性气泡速度研究

本文将Layzer气泡模型推广到任意界面Atwood数情形,得到了自洽的微分方程组.该模型描述了气泡从早期的指数增长阶段到气 泡以渐近速度上升的非线性阶段的发展过程,给出了Rayleigh-Taylor(RT)和Richtmyer-Meshkov(RM)不稳定性的二维和 三维气泡速度渐近解,还求出了二维和三维RT不稳定性气泡顶点附近速度的解析解.     

表面张力对Rayleigh-Taylor不稳定性谐波的影响

Using the method of the parameter expansion up to the third order, explicitly investigates surface tension effect on harmonics at weakly nonlinear stage in Rayleigh-Taylor instability (RTI) for arbitrary Atwood numbers and compares the results with those of classical RTI within the framework of the third-order weakly nonlinear theory. It is found that surface tension strongly reduces the linear growth rate of time, resulting in mild growth of the amplitude of the fundamental mode, and changes amplitudes of the second and third harmonics, as is expressed as a tension factor coupling in amplitudes of the harmonics. On the one hand, surface tension can either decrease or increase the space amplitude; on the other hand, surface tension can also change their phases for some conditions which are explicitly determined.     

Three-dimensional effects in directional solidification in hele-shaw cells: nonlinear evolution and pattern selection

Directional solidification of a dilute binary alloy in a Hele-Shaw cell is modeled by a long-wave nonlinear evolution equation with zero flux and contact-angle conditions at the walls. The basic steady-state solution and its linear stability criteria are found analytically, and the nonlinear system is solved numerically. Concave-down (toward the solid) interfaces under physically realistic conditions are found to be more unstable than the planar front. Weakly nonlinear analysis indicates that subcritical bifurcation is promoted, the domain of modulational instability is expanded and transition to three-dimensional patterns is delayed due to the contact-angle condition. In the strongly nonlinear regime fully three-dimensional steady-state solutions are found whose characteristic amplitude is larger than that for the two-dimensional problem. In the subcritical regime secondary bifurcation to stable solutions is promoted.     

Phase-field modeling of epitaxial growth: Applications to step trains and island dynamics

In this paper, we present a new phase-field model including combined effects of edge diffusion, the Ehrlich-Schwoebel barrier, deposition and desorption to simulate epitaxial growth. A new free energy function together with a correction to the initial phase variable profile is used to efficiently capture the morphological evolution when a large deposition flux is imposed. A formal matched asymptotic analysis is performed to show the reduction of the phase-field model to the classical sharp interface Burton-Cabrera-Frank model for step flow when the interfacial thickness vanishes. The phase-field model is solved by a semi-implicit finite difference scheme, and adaptive block-structured Cartesian meshes are used to dramatically increase the efficiency of the solver. The numerical scheme is used to investigate the evolution of perturbed circularly shaped small islands. The effect of edge diffusion is investigated together with the Ehrlich-Schwoebel barrier. We also investigate the linear and nonlinear regimes of a step meandering instability. We reproduce the predicted scaling law for the growth of the meander amplitude, which was based on an analysis of a long wavelength regime. New nonlinear behavior is observed when the meander wavelength is comparable to the terrace width. In particular, a previously unobserved regime of coarsening dynamics is found to occur when the meander wavelength is comparable to the terrace width.     

三维激光烧蚀瑞利-泰勒不稳定性并行计算

给出了用三维激光烧蚀不稳定性程序求解的物理方程组,简要介绍了所用的数值方法.根据所模拟问题的物理性质和数值模拟方法探索了在多处理机上实行并行计算的技术和方法,并行效率达75%,在此基础上进行了三维多模激光烧蚀不稳定性的数值模拟,深化了对有关物理规律的认识.     

Linear perturbation amplification in self-similar ablation flows of inertial confinement fusion

Exact similarity solutions for inviscid compressible ablative flows in slab symmetry with nonlinear heat conduction are proposed for studying unsteadiness and compressibility effects on the hydrodynamic stability of ablation fronts relevant to inertial confinement fusion. Both the similarity solutions and their linear perturbations are numerically computed with a dynamical multidomain Chebyshev pseudospectral method. Herewith the first analysis of laser-imprinting based on a dynamic solution is presented, showing that maximum perturbation amplification occurs for a laser-intensity modulation of zero transverse wave number, with growth dominated by the mean flow stretching.     

Acoustic phase conjugation in a three-phase medium

Acoustic phase conjugation is studied in a sandy marine sediment that contains air bubbles in its fluid fraction. The considered phase conjugation is a four-wave nonlinear parametric sound interaction caused by nonlinear bubble oscillations which are known to be dominant in acoustic nonlinear interactions in three-phase marine sediments. Two various mechanisms of phase conjugation are studied. One of them is based on the stimulated Raman-type sound scattering on resonance bubble oscillations. The other is associated with sound interactions with bubble oscillations whose frequencies are far from resonance bubble frequencies. Nonlinear equations to solve the phase conjugation problem are derived, expressions for acoustic wave amplitudes with a conjugate wave front are obtained and compared for various frequencies of the excited bubble oscillations.     

Contact and non-contact ultrasonic diagnostics of shock waves driven by single-shot femtosecond laser ablation of titanium

Initial pressures and expansion speeds of air shock waves, driven by single-shot femtosecond laser ablation of titanium surface at fluences up to 6 J/cm², were acquired by means of front-side non-contact time-of-flight ultrasonic measurements. Similar contact ultrasonic measurements on the rear side of the titanium target demonstrate general correlation of ultrasonic amplitudes the measured in the contact and non-contact modes at high ablative pressures (above 100 GPa, superelastic propagation regime), and strong deviation of these dependences at lower ablative pressures (strong dissipation regime) without any indications of new titanium phases, besides the initial α-Ti phase.     

Effect of picosecond laser induced cavitation bubbles generated on Au targets in a nanoparticle production set-up

This work is aimed at an analysis of the influence on the efficiency of nanoparticle production of a cavitation bubble (CB), which forms during the laser ablation process in high-fluence regime. The CB is produced on an Au metal target immersed in water by 1064 nm ps Nd:YAG laser pulses at different fluences. Its time–space evolution is monitored by a shadowgraphic set-up, while the Au nanoparticles production rate is tagged by the growth of the plasmon resonance, which is detected by measuring shot-by-shot the UV-Vis absorbance. We analyze the dependence of bubble size on the experimental parameters. Our results appear of interest to enhance the nanoparticle production efficiency in a liquid medium.