Rayleigh Taylor Instability

  Rayleigh Taylor Instability

AIM:

• Study of some practical CFD models that have been based on the mathematical analysis of Rayleigh Taylor waves
• Perform the Rayleigh Taylor instability simulation for 2 different mesh sizes with the base mesh being 0.5 mm.
• explain why a steady-state approach might not be suitable for these types of simulation.[ Choose Air and Water for this simulation]
• Run one more simulation with water and user-defined material(density = 400 kg/m3, viscosity = 0.001 kg/m-s) for refined mesh. 
• Define the Atwood Number. Find out the Atwood number for both the cases and explain how the variation in Atwood number in the above two cases affects the behavior of the instability. 

 Rayleigh Taylor Instability :

Introduction:

Rayleigh Taylor Instability is an instability of the interface betweeen the two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid. It is clear from Newtons law of viscosity,fluid cannot resist the slightest shear force and Density is one of the main property of the fluid to incorporate the shear force. Based on the density, when lighter and heavier fluids are mixed together, denser fluid comes to top of the lighter fluid, due to the gravity and the density variation at the interface of the liquid is not stable. This instability of the interface between the different densities of fluid when the lighter fluid is accelerated towards the heavier fluid is called  RAYLEIGH TAYLOR INSTABILITY.

Examples:

• Formation of high luminosity twin exhaust jets in rotating gas clouds in an external gravitational potential,
• Instabilities in plasma fusion reactors and intertial confinement fusion,
• Electromagnetic implosion of metal liner
• Behaviour of water that is suspended above oil in gravity of Earth

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

Model	Description
Ritchmyer-Meshkov instability:	The Richtmyer-Meshkov instability arises  to the instability of an interface between the two fluids that is impulsively accelerated, generally by a shock wave. When a shock wave interacts with an interface separating two different fluids hydrodynamic instability occurs, including nonlinear growth and subsequent transition to turbulence, across a wide range of Mach numbers. Example: Supersonic combustion in a scramjet, Combustion, supernova, hypersonic air-breathing engines In combustion, Instability begins with a very small amplitude perturbations which initially grow linearly with respect to time. The instability can be considered the implusive-acceleration limit of the Rayleigh Taylor Instability. During the implosion of an inertia confinement fusion target , the hot shell material surrounding the cold D-T fuel layer is shock accelerated. Such instability is visualised in Magnetic target fusion.
Turbulence model	It is required to predict the average mixing behaviour in flows that are on average one- or two-dimensional. The governing equations of turbulent flows can be solved directly for simple cases of flow. For turbulent flows, CFD simulation uses turbulent models to predict the evolution of turbulence and are simplified constitutive equations that predict the statistical evolution of turbulent flows.
Plateau–Rayleigh instability:	It is also known as Rayleigh instability, and explains why the falling stream of fluid breaks up into smaller packets with the same volume but less surface area. This fluid instability is exploited in the design of a particular type of ink jet technology whereby a jet of liquid is perturbed into a steady stream of droplets.
Kelvin–Helmholtz instability:	This instability is observed when there is velocity shear in a single continuous fluid or additionally where there is a velocity difference across the interface between the two fluids. Kelvin-Helmholtz  studied the dynamics of two fluids of different densities subjected to a small disturbance such as wave was introduced at the boundary connecting both the fluids. This phenomena is mostly responsible for natures basic structures. It can occur when there is velocity shear in a single continuos fluid, or where there is a velocity difference across the interface between two fluids. Examples: wind blowing over water, clouds, the ocean, Saturn's bands, and the sun's corona
Smoother particle hydrodynaics (SPH) method	It is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and fluid flows. SPH is  increasingly used to model fluid motion and it has several benefits over traditional grid-based techniques. It guarantees conservation of mass without extra computation as the particles themselves represent mass. SPH also computes pressure from weighted contributions of neighboring particles rather than by solving linear systems of equations. Furthermore, unlike grid-based techniques, which must track fluid boundaries, SPH creates a free surface for two-phase interacting fluids directly since the particles represent the denser fluid  and empty space represents the lighter fluid.
Single Fluid Model	The analysis of two-phase flows (mixture model) is done using single fluid model. In this model individual fluid phases are assumed to behave as a flowing mixture described in terms of the mixture properties. It is a five-equation model consisting of the mass, momentum and energy equations for a vapor/liquid mixture, and two equations describing the formation and growth of the liquid phase.
Two Fluid Model	This model separates the sets of governing equation for the vapor and liquid phases but the interaction between the droplets and  heat exchange between the liquid phase and the solid boundary are not modelled. In addition the velocity slip between vapor and the liquid phase is in this model taken into account.
Descrete Boltzman Model	It is used for premixied, non-premixed or partially pri-mixed n0n equilibrium reactive flows. It is suitable for both subsonic and supersonic flows with/without chemical reaction and external force. It is developed to investigate the hydrodynamic and thermodynamic non-equilibrium effects in phase seperation processes. Based on the thermodynamic non-equilibrium features, its estimates the deviation amplitude from the thermodynamic  equilibrium.

Performing Rayleigh Taylor instability simulation:

1.For 2 different mesh sizes with the base mesh being 0.5 mm(Choose Air and Water for this simulation)

2.One more simulation with water and user-defined material(density = 400 kg/m3, viscosity = 0.001 kg/m-s) for refined mesh. 

Explanation:

1.To perform this operation we have to make a Geometry in ansys and then proceed.

Geometry:

First in geometry choose DEsign-Mode-Sketch Mode from 3D Mode.Then ctrl+6 to choose front face . Then in Sketch mode choose Rectangle and make a rectangle of 20mm*20mm (After selecting rectangle choose sketch mode at bottom side then tap at any point on the plane ).Then in Sketch-Edit choose pull tool to select the sketch as a plane and then escape.

Then again in Design change the mode into Sketch mode and in Sketch select rectangle and click on the previous rectangle to select the plane and draw same rectangle above the previously drawn rectangle.Then again pull tool and select the sketch.

 

Then in the upper left corner rename the two surfaces , rename upper surface as water and lower as Air.Then go to workbench and select overlap bodies and select the geometry and then select the share option to share common part of the two geometries,then click on the green tick and then click escape. We can see at the left bottom the topology has been shared and we will get a message one common edge found.

Geometry has been created so then close it open meshing.

MESH:

In mesh click on mesh - right click -generate mesh . Keep the element size as 0.5mm .

Then close the mesh and open set up.

SETUP:

In Domain we choose pressure-based transient solver and choose gravity and in y axis choose -9.81 as gravity works downward.Then in physics choose viscous medium as laminar,choose Muliphase choose Model-volume of fluid,Number of Eulerian phases-2,Formulation-impliit,choose Phases-air as first phase and water as second phase. Before choosing phases choose material and then in fluent database select water liquid copy it and close it and then in mulyiphase choose phases.

Then in solution first initialize then we have to patch the solution as we have choosed a mixture then we have to choose an alpha value that will indicate how much percent air or water in the mixture we have.In patch choose phases-water,variable-volume fraction,value-1,zone-rayleigh_taylor_model-water,then patch. we have to remember that after every initialization we have to patch.

Then choose timestep as 0.01 and no of iteration 400

we have to make an animation of the solution so in results section choose contour-new and then choose the following below and then in solution-create and choose solution animation.

Results:

case-1

element size 0.5mm

Frame-8

Frame 70:

 

Frame 110:

Frame 400:

 

Residual:

 

ANIMATION LINK:

https://www.youtube.com/watch?v=KARlg7VB1nU

Explanation:

it is seen from the contour that lower density fluid (Air) is setteled on top and the higher density fluid (bottom) is settled at the bottom. This is formed due to the instability. This growing instabilities are known as Rayleigh Taylor Instability. We could see the formation of bubbles whose position is later replaced by water. Simulation begins with fluids of two phases at hydrostatic equilibrium state and then Rayleigh Taylor instability is seen at the interference when the lower density fluid pushes a higher density fluid due to which formation of shock waves occur at interface.The disturbance at the interface ends with lighter fluid settling up and the heavier fluid settles down. 

 

element size 0.4mm

Frame-8:

 

Frame-70:

 

Frame-110:

Frame-400:

Residual:

 

Animation link:https://www.youtube.com/watch?v=G4BKoddFpB8

 

element size 0.3mm

Frame-8:

 

Frame-70:

Frame-110:

 

Frame-400:

 

Residual:

Animation link:

https://www.youtube.com/watch?v=nqw6xQEkOiI

Explanation:

It is clear from the results that instabilities that are occuring at the interface ie, formation of bubbles and vortices are explicitly seen as we refine the mesh quality or finer mesh.

Steady State Approach: It is used when we are only interested in end result and the intermittent behaviour is not important. Whereas in this case we are interested in the real-time changes at the interface as well as to see the changes in the fluid region.  This real  time observations  can only be studied in Transient state Analysis. We are interested in the behaviour at interface when the fluids of various densities collapses as a function of time and this time step is too small to capture explicit behaviour the Instabilitty effect known as Rayleigh Taylor Instability. Thus Transient approach was selected.

Case-2

(water and user_defined material)

Frame-8

Frame-70:

Frame-110:

Frame-400:

Residual:

Animation Link:

https://www.youtube.com/watch?v=CRXhLhxdBpo

Explanation:

It is seen from the results that simulation begins with fluids of two phases at hydrostatic equilibrium state and then Rayleigh Taylor instability is seen at the interference when the lower density fluid pushes a higher density fluid due to which formation of shock waves occur at interface. This in turn lead to formation of bubble like but later it becomes assymetric, which compresses the heavy fluid around it due to shock waves that pushes heavy fluid updwards. Later disturbance at the interface ends with lighter fluid settling up and the heavier fluid settles down. 

 

Atwood Number:  It 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:

Iti is paramount in the study of Rayleigh-Taylor instability. 

Atwood number close to 0, RT instability flows take the form asymmetric “finger” of fluid; 

Atwood number close to 1, the much lighter fluid “below” the heavier fluid takes the form of larger bubble-like plumes.

Explanation:

It is observed from the simulations that air forms larger bubbles when water is burst down the air and forms large bubbles. The Atwood number calculated in first case is 1 (approx.) and other case is 0.42. Rayleigh Taylor instability shows three distinct stages namely  Linear stability ,  bubble(rising) formation and long term evolution due to bubble merging and mixing.

For low Atwood number, bubble and mushroom head shape is relatively symmetrical, the disturbances are linear for a larger time. In contrast, for large Atwood number, bubble and mushroom head shape is assymetrical and occurs rapidly at an accelerating rate leading to non linear growth rate.

Hence simulation results and Atwood calculation results are similar,  thus validating the results.

Conclusion:

Rayleigh Taylor Instability simulation was studied at interface of both fluids having different densities under the influence of gravity and the following observations were made:

• Finer mesh shows formation of bubbles and vortices explicity
• Transient state analysis is preffered instead of steady state analysis as we are interested at real-time changes at the interface
• The Atwood number depends on the density of the fluid and ranges from 0 to 1.
• Both cases have larger Atwood number showing the Rayleigh Taylor Instability.  Thus matching the simulation results.