Rayleigh Taylor Instability

Objective:

Q1. What are some practical CFD models that have been based on the mathematical analysis of Rayleigh Taylor waves? In your own words, explain how these mathematical models have been adapted for CFD calculations

Q2. Perform the Rayleigh Taylor instability simulation for 3 different mesh sizes with the base mesh being 0.5 mm. Compare the results by showing the animations. Also, explain why a steady-state approach might not be suitable for these types of simulation

Q3. Define the Atwood Number. Find out the Atwood number for this case and explain how the variation in Atwood number affects the behavior of the instability. How can we use the Atwood number to validate our simulation result

Introduction:

Rayleigh Taylor Instability:

RT instability is the instability of an interface between two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid.

Q1. What are some practical CFD models that have been based on the mathematical analysis of Rayleigh Taylor waves?

• Plateau-Rayleigh Taylor Instability:

The Plateau–Rayleigh instability, often just called the Rayleigh instability, explains why and how a falling stream of fluid breaks up into smaller packets with the same volume but less surface area. It is related to the Rayleigh–Taylor instability and is part of a greater branch of fluid dynamics concerned with fluid thread breakup. This fluid instability is exploited in the design of a particular type of inkjet technology whereby a jet of liquid is perturbed into a steady stream of droplets.

• Kelvin-Helmholtz Instability:

The Kelvin–Helmholtz instability (after Lord Kelvin and Hermann von Helmholtz) can occur when there is velocity shear in a single continuous fluid, or where there is a velocity difference across the interface between two fluids. An example is a wind blowing over water: The instability manifests in waves on the water surface. More generally, clouds, the ocean, Saturn\'s bands, Jupiter\'s Red Spot, and the sun\'s corona show this instability.

• Richtmyer-Meshkov Instability:

The Richtmyer–Meshkov instability (RMI) occurs when two fluids of different density are impulsively accelerated. Normally this is by the passage of a shock wave. The development of the instability begins with small amplitude perturbations which initially grow linearly with time. This is followed by a nonlinear regime with bubbles appearing in the case of a light fluid penetrating a heavy fluid, and with spikes appearing in the case of a heavy fluid penetrating a light fluid. A chaotic regime eventually is reached and the two fluids mix. This instability can be considered the impulsive-acceleration limit of the Rayleigh–Taylor instability

The Richtmyer-Meshkov instability arises when a shock wave interacts with an interface separating two different fluids. It combines compressible phenomena, such as shock interaction and refraction, with hydrodynamic instability, including nonlinear growth and subsequent transition to turbulence, across a wide range of Mach numbers. This review focuses on the basic physical processes underlying the onset and development of the Richtmyer-Meshkov instability in simple geometries. It examines the principal theoretical results along with their experimental and numerical validation. It also discusses the different experimental approaches and techniques and how they can be used to resolve outstanding issues in this field.

• Varicose Instabilities:

An investigation of a model of turbulence generation in the wall region of a turbulent boundary layer is made through direct numerical simulations. The model is based on the varicose instability of a streak. First, a laminar boundary layer disturbed by a continuous blowing through a slot is simulated in order to reproduce and further investigate the results reported from the experiments of Acarlar and Smith [J. Fluid Mech. 175, 43 (1987)]. An isolated streak with an inflectional profile is generated that becomes unstable, resulting in a train of horseshoe vortices. The frequency of the vortex generation is equal to the experimental results. Comparison of the instability characteristics to those predicted through an Orr–Sommerfeld analysis is in good agreement. Second, direct numerical simulation of a turbulent boundary layer is performed to point out the similarities between the horseshoe vortices in a turbulent and a laminar boundary layer. The characteristics of streaks and the vortical structures surrounding them in a turbulent boundary layer compare well with the model streak. The results of the present study show that one mechanism for the generation of horseshoe vortices in turbulent boundary layers is related to a normal inflectional instability of the streaks.

• Richtmyer-Meshkov Instabilities:

The Richtmyer–Meshkov instability (RMI) occurs when two fluids of different density are impulsively accelerated. Normally this is by the passage of a shock wave. The development of the instability begins with small amplitude perturbations which initially grow linearly with time. This is followed by a nonlinear regime with bubbles appearing in the case of a light fluid penetrating a heavy fluid, and with spikes appearing in the case of a heavy fluid penetrating a light fluid. A chaotic regime eventually is reached and the two fluids mix. This instability can be considered the impulsive-acceleration limit of the Rayleigh–Taylor instability

Pre-Processing:

Geometry:

Size: 20*20

Case: 1

• Base Mesh
• Nodes: 3321
• Elements: 3200

Residual:

Volume Fraction of Water:

Case: 2

• Base Mesh: 0.3
• Nodes: 9180
• Elements: 8978

Residual:

Volume Fraction of Water:

Case: 3

• Base Mesh: 0.1
• Nodes: 80601
• Elements: 80000

Residual:

Volume Fraction of Water:

Conclusion:

• Rayleigh-Taylor instability The fluid with a higher density is located above the fluid with a smaller density. Due to the gravity in the system fluid with higher density moves below the fluid with lower density.
• Where there is lighter fluid pushes to the heavier fluid because of that the water comes down and air goes up.
• When During this processes air bubbles start taking place, On account of their lower density with respect to the outer environment, the bubbles rise buoyantly and sweep up some of the intracluster medium (ICM), as evident by the enhanced emission at their rims.
• It would seem that the top of these bubbles is prone to the Rayleigh–Taylor instability. Indeed, the gravitational acceleration of the air pushes downwards, and the dense swept-up gas shell lies on the top of the thinner one inside the bubble itself.
• Because the shell density ρ_s is much higher than the density ρ_b inside the bubble, and since the gravitational acceleration is directed downward, the flow at the bubble\'s top would be RT unstable if only gravity is considered. This simple picture, however, overlooks the bubble\'s inflation process.
• Since the mesh refined from Base mesh, the simulation results become very smother at case 1 and case 2 compared to baseline mesh. From visualizing the simulation results become very smothering.
• The water molecule itself, H₂O, is actually less massive than a Nitrogen molecule, N₂, which makes up 80% of our air. It’s just that there are more water molecules in one liter than there are air molecules…even when the air is pressurized. Think of it this way: if you pump up the bottle to 140 pounds per square inch (PSI) of air pressure, then for every air molecule expelled out of the bottle, you could have expelled roughly 83 water molecules. That’s a lot more mass!
• this processes initially, the air has lied at the bottom during the simulation the air traveled to an upper region and water comes down due to gravity of the fluid. At the day the results divided into three phases. Below having water and above the air in between two fluids some diffusivity forming due to lost volume fraction of air and water. Where the air is contained some water droplets often we called water vapor or humidity.

Q3. Define the Atwood Number?

The Atwood number (A) is a dimensionless number in fluid dynamics used in the study of hydrodynamic instabilities in density stratified flows. It is a dimensionless density ratio defined as

`At = (rho_1-rho_2)/(rho_1+rho_2)`

where

`rho_1`= density of the heavier fluid

`rho_2 =` density of the lighter fluid

`At = (998-1.2)/(998+1.2)`

At = 0.997

Atwood number is an important parameter in the study of Rayleigh–Taylor instability and Richtmyer–Meshkov instability. In Rayleigh–Taylor instability, the penetration distance of heavy fluid bubbles into the lighter fluid is a function of acceleration time scale.

• A small density differences, At < 0.1, up-welling (‘bubble’) and down-welling
  (‘spike’) disturbances form a symmetrical mixing region that expands away
  from the initial position of the density interface and is based on the reviewers’ experience rather than on any extensive reported study; indeed, asymmetries in the penetration of rising bubbles and falling spikes are not readily apparent until At > 0.5
• As the Atwood number increases, the spike penetration is greater than that of the bubble owing to a narrowing of the spike, and a consequent reduction of drag on the head of the spike.
• In this case, The Atwood number is 0.997 nearly 1, Atwood regime above A > 0.90, where density differences may be significantly closer to conditions seen in many engineering and scientific applications.
• The first simulations of RTI at Atwood numbers A>0.90 are then used to explore the effects of varying the density ratio on mixing dynamics at very high Atwood numbers. Spike heights and mixing layer growth rates are shown to be strongly affected by the initial density ratio.
• At high Atwood number, however, spike stems are shown to produce forward transfer due to compressive straining generated from the larger spike penetration velocity and the differing spike morphology produced at the higher Atwood number.