Week 5 - Rayleigh Taylor Instability

Rayleigh Taylor Instability :

                                     The Rayleigh–Taylor instability is an instability  of an interface between two fluids of different

densities which occurs  when the lighter fluid is pushing the heavier fluid. Examples include the behaviour of water

suspended above oil in the gravity of Earth  ,mushroom clouds like those from volcanic eruptions and atmospheric nuclear

explosion, supernova explosions in which expanding core gas is accelerated into denser shell gas, instabilities in plasma

fusion reactors and inertial confinement fusion.

 Some practical CFD models that have been based on the mathematical analysis of Rayleigh Taylor waves are:

Kelvin Helmholtz instability Model:

                                                 This instability can occur when there is velocity shear in a continuous  fluid, or where there

is a velocity difference across the interface between two fluids. An example is wind blowing over the water. other examples

showing this instability  are clouds,oceans,saturns bands.

    2. Richtmeyer  meshkow instability :

                                                  This instability occurs when a shock wave interacts with an interface separating two

different fluids or in other words when two fluids of different densities are impulsively accelerated.The development of

instability begins with small amplitude pertubations which initially grow linear with time.This instability can be considered the

impulsive acceleration limit of the Rayleigh Taylor Instability.

    3. Plateau Rayleigh instability :

                                               This instability explains why and how a falling stream of fluid breaks up into smaller packets

with the same volume but different or small surface area . It is often called as Rayleigh instability.The driving force of this

instability is that the liquids by virtue of their surface tensions ,tends to minimize their surface area .       

   4.Rayleigh Taylor Instability :

                                            The Rayleigh–Taylor instability is an instability  of an interface between two fluids of different

densities which occurs  when the lighter fluid is pushing the heavier fluid.

Case 1) Rayleigh Taylor instability simulation for base mesh 0.5 mm :

Step 1) Geometry creation :

   This is 2D simulation so we create the two square of size 20X20mm one is above the other and we give the name of top

square is water and bottom square is air. and use share topology due to that the interface edge is shared .The created

geometry is shown below:

 

Step 2) Mesh : Use mesh size 0.5 mm 

The mesh is shown below :

 

Step 3) Setup and solution:

 

Viscous model: Laminar 

Material : Selected Air and water 

Phase selections :

 

 

Method of initialization: Standard 

After initialization we required to patch the two phases of the fluid :

Create a contour :

 

By using contour we create the animation of every time step :

Time step size (sec)= 0.005

No of Time step : 1000 

Result:

Residual plot:

 

Case 2) Mesh size: 0.3 mm

Time step size (sec)= 0.005

No of Time steps : 800 

Result:

Residual plot:

 

Case 3) Mesh size: 0.2 mm

 

 Time step size (sec)= 0.005

 No of Time steps : 800

Result:

Residual plot:

 

Case 4) Run the simulation with water and user-defined material(density = 400 kg/m³, viscosity = 0.001 kg/m-s) for refined

mesh. ie Mesh size =0.2mm 

 Time step size (sec)= 0.005

 No of Time steps : 800

Result :

Residual plot:

Observations :

The finer mesh captures the instability better than the coarse mesh.but it takes more time for the heavier fluid to settle

down.

The computational time is also increased for fine mesh but the result obtained having better accuracy.

The steady state simulation is performed if we are concerned more about the final state results or the equilibrium state,In RT

instability CFD model we are more concerned to learn about the transition of the irregularities that start developing when we

pour the high densiy fluid upon the low density fluid under the effect of gravity.and we want to know how the interface

collapses as a function of time so In this case Transient state is more suitable over the steady state.

Atwood Number:

It is a dimensionless number that is used to study hydrodynamic instabilities  in density stratified flows.

It is denoted by letter A

  $A=\frac{{\rho }_1-{\rho }_2}{{\rho }_1+{\rho }_2}$

      $\to 1$

 

where,

${\rho }_1$

=Density of heavier fluid `

${\rho }_2$

=Density of lighter fluid

Calculate the Atwood Number for Air-water RT instability

${\rho }_1$

=998.2 $\frac{kg}{m^3}$

(water density)

${\rho }_2$

=1.225  $\frac{kg}{m^3}$

(Air density)

put this value in equation 1 we get A=0.9975 

Calculate the Atwood Number for Water-UD material RT instability

${\rho }_1=$

998.2 $\frac{kg}{m^3}$

(Water density)

${\rho }_2$

=400 $\frac{kg}{m^3}$

(User defined material density)

put this value in equation 1 we get A=0.4278 

The Atwood number close to 1 then the lighter fluid below the heavier fluid which forms a larger bubbles which is less

symmetrical and also it occurs very rapidly.

The Atwood number close to 0 then the lighter fluid below the heavier fluid which forms a  bubbles which is relatively

symmetrical and the pertubations are linear for larger time.