Rarefactiondriven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime
Abstract
Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Singlemode twodimensional, and singlemode threedimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a nonconstant acceleration, and a time decreasing Atwood number,$$A=(\unicode[STIX]{x1D70C}_{2}\unicode[STIX]{x1D70C}_{1})/(\unicode[STIX]{x1D70C}_{2}+\unicode[STIX]{x1D70C}_{1})$$, where$$\unicode[STIX]{x1D70C}_{2}$$and$$\unicode[STIX]{x1D70C}_{1}$$are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of$A=0.49$$,$$A=0.63$$,$$A=0.82$$and$$A=0.94$$. Nominally twodimensional (2D) experiments (initiated with nearly 2D perturbations) and 2D simulations are observed to approach an intermediatetime velocity plateau that is in disagreement with the latetime velocity obtained from the incompressible model of Goncharov (Phys. Rev. Lett., vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2D bubbles in large wavenumber,$$k=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D706}=0.247~\text{mm}^{1}$$, experiments and simulations, where$$\unicode[STIX]{x1D706}$is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These latetime velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian (Phys. Fluids., vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a freefall like behaviour. Finally, experiments initiated with threedimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.
 Authors:
 Univ. of Arizona, Tucson, AZ (United States)
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Publication Date:
 Research Org.:
 Univ. of Arizona, Tucson, AZ (United States)
 Sponsoring Org.:
 USDOE National Nuclear Security Administration (NNSA)
 OSTI Identifier:
 1419851
 Grant/Contract Number:
 NA0002929
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Journal of Fluid Mechanics
 Additional Journal Information:
 Journal Volume: 838; Journal ID: ISSN 00221120
 Publisher:
 Cambridge University Press
 Country of Publication:
 United States
 Language:
 English
 Subject:
 42 ENGINEERING; gas dynamics; nonlinear instability; turbulent mixing
Citation Formats
Morgan, R. V., Cabot, W. H., Greenough, J. A., and Jacobs, J. W. Rarefactiondriven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime. United States: N. p., 2018.
Web. doi:10.1017/jfm.2017.893.
Morgan, R. V., Cabot, W. H., Greenough, J. A., & Jacobs, J. W. Rarefactiondriven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime. United States. doi:10.1017/jfm.2017.893.
Morgan, R. V., Cabot, W. H., Greenough, J. A., and Jacobs, J. W. 2018.
"Rarefactiondriven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime". United States.
doi:10.1017/jfm.2017.893.
@article{osti_1419851,
title = {Rarefactiondriven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime},
author = {Morgan, R. V. and Cabot, W. H. and Greenough, J. A. and Jacobs, J. W.},
abstractNote = {Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Singlemode twodimensional, and singlemode threedimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a nonconstant acceleration, and a time decreasing Atwood number,$A=(\unicode[STIX]{x1D70C}_{2}\unicode[STIX]{x1D70C}_{1})/(\unicode[STIX]{x1D70C}_{2}+\unicode[STIX]{x1D70C}_{1})$, where$\unicode[STIX]{x1D70C}_{2}$and$\unicode[STIX]{x1D70C}_{1}$are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of$A=0.49$,$A=0.63$,$A=0.82$and$A=0.94$. Nominally twodimensional (2D) experiments (initiated with nearly 2D perturbations) and 2D simulations are observed to approach an intermediatetime velocity plateau that is in disagreement with the latetime velocity obtained from the incompressible model of Goncharov (Phys. Rev. Lett., vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2D bubbles in large wavenumber,$k=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D706}=0.247~\text{mm}^{1}$, experiments and simulations, where$\unicode[STIX]{x1D706}$is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These latetime velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian (Phys. Fluids., vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a freefall like behaviour. Finally, experiments initiated with threedimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.},
doi = {10.1017/jfm.2017.893},
journal = {Journal of Fluid Mechanics},
number = ,
volume = 838,
place = {United States},
year = 2018,
month = 1
}

Evidence for a BubbleCompetition Regime in Indirectly Driven Ablative RayleighTaylor Instability Experiments on the NIF
In this paper, we investigate on the National Ignition Facility the ablative RayleighTaylor instability in the transition from weakly nonlinear to highly nonlinear regimes. A planar plastic package with preimposed twodimensional broadband modulations is accelerated for up to 12 ns by the xray drive of a gasfilled Au radiation cavity with a radiative temperature plateau at 175 eV. This extended tailored drive allows a distance traveled in excess of 1 mm for a 130 μm thick foil. Measurements of the modulation optical density performed by xray radiography show that a bubblemerger regime for the RayleighTaylor instability at an ablation frontmore »Cited by 10 
Interface width effect on the classical RayleighTaylor instability in the weakly nonlinear regime
In this paper, the interface width effects (i.e., the density gradient effects or the density transition layer effects) on the RayleighTaylor instability (RTI) in the weakly nonlinear (WN) regime are investigated by numerical simulation (NS). It is found that the interface width effects dramatically influence the linear growth rate in the linear growth regime and the mode coupling process in the WN growth regime. First, the interface width effects decrease the linear growth rate of the RTI, particularly for the short perturbation wavelengths. Second, the interface width effects suppress (reduce) the thirdorder feedback to the fundamental mode, which induces themore » 
Preheating ablation effects on the RayleighTaylor instability in the weakly nonlinear regime
The twodimensional RayleighTaylor instability (RTI) with and without thermal conduction is investigated by numerical simulation in the weakly nonlinear regime. A preheat model {kappa}(T)={kappa}{sub SH}[1+f(T)] is introduced for the thermal conduction [W. H. Ye, W. Y. Zhang, and X. T. He, Phys. Rev. E 65, 057401 (2002)], where {kappa}{sub SH} is the SpitzerHaerm electron thermal conductivity coefficient and f(T) models the preheating tongue effect in the cold plasma ahead of the ablation front. The preheating ablation effects on the RTI are studied by comparing the RTI with and without thermal conduction with identical density profile relevant to inertial confinement fusionmore » 
Harmonic growth of spherical RayleighTaylor instability in weakly nonlinear regime
Harmonic growth in classical RayleighTaylor instability (RTI) on a spherical interface is analytically investigated using the method of the parameter expansion up to the third order. Our results show that the amplitudes of the first four harmonics will recover those in planar RTI as the interface radius tends to infinity compared against the initial perturbation wavelength. The initial radius dramatically influences the harmonic development. The appearance of the secondorder feedback to the initial unperturbed interface (i.e., the zeroth harmonic) makes the interface move towards the spherical center. For these four harmonics, the smaller the initial radius is, the faster theymore »