Search

Rayleigh-Taylor instability (RTI) under radiation background is commonly found in both engineering applications and natural phenomena. In the optically thin and incompressible limit, the corresponding problem can be simplified as an interface discontinuous acceleration (IDA) RTI problem, but to date has only been studied in the linear stage. In this paper, the entire IDA-RTI evolution was studied numerically and theoretically, particularly for the stages beyond the linear stage. The results show that the IDA-RTI problem is equivalent to the classical RTI with the effective acceleration g_{eff}^{*} that is introduced in this work. Moreover, our studies further show that IDA-RTI can occur if and only if g_{eff}^{*}>0 (from heavy fluid to light fluid). This criterion means that IDA-RTI can occur when (i) heavy fluid supports (or accelerates) the light fluid or (ii) the two fluids have the same density, in contrast to the classical RTI problem. Moreover, the quasisteady bubble and spike velocities are theoretically predicted with quantitative accuracy, showing good agreement with the results of numerical simulations in a wide range of density ratios and acceleration configurations.

Implicit large eddy simulations of two-dimensional Rayleigh-Taylor instability at different density ratios (i.e., Atwood number A=0.05, 0.5, and 0.9) are conducted to investigate the late-time dynamics of bubbles. To produce a flow field full of bounded, semibounded, and chaotic bubbles, three problems with distinct perturbations are simulated: (I) periodic sinusoidal perturbation, (II) isolated W-shaped perturbation, and (III) random short-wave perturbations. The evolution of height h, velocity v, and diameter D of the (dominant) bubble with time t are formulated and analyzed. In problem I, during the quasisteady stage, the simulations confirm Goncharov's prediction of the terminal speed v_{∞}=Frsqrt[Agλ/(1+A)], where Fr=1/sqrt[3π]. Moreover, the diameter D at this stage is found to be proportional to the initial perturbation wavelength λ as D≈λ. This differed from Daly's simulation result of D=λ(1+A)/2. In problem II, a W-shaped perturbation is designed to produce a bubble environment similar to that of chaotic bubbles in problem III. We obtain a similar terminal speed relationship as above, but Fr is replaced by Fr_{w}≈0.63. In problem III, the simulations show that h grows quadratically with the bubble acceleration constant α≡h/(Agt^{2})≈0.05, and D expands self-similarly with a steady aspect ratio β≡D/h≈(1+A)/2, which differs from existing theories. Therefore, following the mechanism of self-similar growth, we derive a relationship of β=4α(1+A)/Fr_{w}^{2} to relate the evolution of chaotic bubbles in problem III to that of semibounded bubbles in problem II. The validity of this relationship highlights the fact that the dynamics of chaotic bubbles in problem III are similar to the semibounded isolated bubbles in problem II, but not to that of bounded periodic bubbles in problem I.

Turbulent mixing induced by Rayleigh-Taylor (RT) instability occurs ubiquitously in many natural phenomena and engineering applications. As the simplest and primary descriptor of the mixing process, the evolution of mixing width of the mixing zone plays a notable role in the flows. The flows generally involve complex varying acceleration histories and widely varying density ratios, two dominant factors affecting the evolution of mixing width. However, no satisfactory theory for predicting the evolution has yet been established. Here a theory determining the evolution of mixing width in general RT flows is established to reproduce, first, all of the documented experiments conducted for diverse (i.e., constant, impulsive, oscillating, decreasing, increasing, and complex) acceleration histories and all density ratios. The theory is established in terms of the conservation principle, with special consideration given to the asymmetry of the volume-averaged density fields occurring in actual flows. The results reveal the sensitivity or insensitivity of the evolution of a mixing front of a neighboring light or heavy fluid to the degree of asymmetry and thus explain the distinct evolutions in two experiments with the same configurations.