ABSTRACT
Two-dimensional numerical simulations for the Rayleigh-Taylor instability in an elastic-plastic medium are presented. Recent predictions of the theory regarding the asymmetric growth of peaks and valleys during the linear phase of the instability evolution are confirmed. Extension to the nonlinear regime reveals singular features, such as the long delay in achieving the nonlinear saturation and an intermediate phase with growth rate larger than the classical one.
ABSTRACT
The study of the linear stage of the incompressible Rayleigh-Taylor instability in elastic-plastic solids is performed by considering thick plates under a constant acceleration that is also uniform except for a small sinusoidal ripple in the horizontal plane. The analysis is carried out by using an analytical model based on the Newton second law and it is complemented with extensive two-dimensional numerical simulations. The conditions for marginal stability that determine the instability threshold are derived. Besides, the boundary for the transition from the elastic to the plastic regime is obtained and it is demonstrated that such a transition is not a sufficient condition for instability. The model yields complete analytical solutions for the perturbation amplitude evolution and reveals the main physical process that governs the instability. The theory is in general agreement with the numerical simulations and provides useful quantitative results. Implications for high-energy-density-physics experiments are also discussed.
Subject(s)
Linear Models , Rheology/methods , Computer Simulation , Phase TransitionABSTRACT
We present an analytical model for the Rayleigh-Taylor instability that allows for an approximate but still very accurate and appealing description of the instability physics in the linear regime. The model is based on the second law of Newton and it has been developed with the aim of dealing with the instability of accelerated elastic solids. It yields the asymptotic instability growth rate but also describes the initial transient phase determined by the initial conditions. We have applied the model to solid/solid and solid/fluid interfaces with arbitrary Atwood numbers. The results are in excellent agreement with previous models that yield exact solutions but which are of more limited validity. Our model allows for including more complex physics. In particular, the present approach is expected to lead to a more general theory of the instability that would allow for describing the transition to the plastic regime.
