Week 5 - Rayleigh Taylor Instability
Rayleigh Taylor Instability
Objective:
- 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.
- Perform the Rayleigh Taylor instability simulation for 2 different mesh sizes with the base mesh being 0.5 mm. Compare the results by showing the animations( Attach them from your google drive). Observe the differences in the animation and explain the effect of mesh size on simulation. Also, 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.
the mathematical analysis of Rayleigh Taylor waves
- Diesel Spray Simulation and KH-RT Wave Model 2003-01-3231 :This study presents diesel spray breakup regimes and the wave model basic theory from literature. The RD wave model and the KH-RT wave model are explained. The implementation of the KH-RT wave model in a commercial CFD code is briefly presented. This study relies on experimental data from non-evaporating sprays that have earlier been measured at Helsinki University of Technology.The simulated fuel spray in a medium-speed diesel engine had a satisfactory match with the experimental data. The KH-RT wave model resulted in a much faster drop breakup than with the RD wave model. This resulted in a thin spray core with the KH-RT model. The fuel viscosity effect on drop sizes was well predicted by the KH-RT wave model.
- 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
-
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 ink jet technology whereby a jet of liquid is perturbed into a steady stream of droplets.
The driving force of the Plateau–Rayleigh instability is that liquids, by virtue of their surface tensions, tend to minimize their surface area. A considerable amount of work has been done recently on the final pinching profile by attacking it with self-similar solutions.[1][
- The Kelvin–Helmholtz instability (after Lord Kelvin and Hermann von Helmholtz) is a fluid instability that occurs when there is velocity shear in a single continuous fluid or a velocity difference across the interface between two fluids. Kelvin-Helmholtz instabilities are visible in the atmospheres of planets and moons, such as in cloud formations on Earth or the Red Spot on Jupiter, and the atmospheres of the Sun and other star
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
- A = ρ 1 − ρ 2 ρ 1 + ρ 2 {\displaystyle \mathrm {A} ={\frac {\rho _{1}-\rho _{2}}{\rho _{1}+\rho _{2}}}}

where
- ρ 1 {\displaystyle \rho _{1}}
= density of heavier fluid - ρ 2 {\displaystyle \rho _{2}}
= density of lighter fluid
stimulation of Rayleigh Taylor Instability
case 1: base line stimulation
steps invalving in analysis are shown below
- geometry
- meshing
- setup
- solution
- result
geometry: we have to create the two rectangular surfaces as shown below
- for air and water of demensions 20mmX20mm
- after taht we have to share the information between two layes as shown below

meshing: here first we have to create the basic mesh as show below

setup:here we are using transient stete and presure solver as shown below

here we are using k-omega based viscouse model




solution:
residuals:

valume fraction:

animation link:https://1drv.ms/u/s!AvvLa2FhdyA4tjVpaC6kuT0QW6QS
case 2: meshing 0.5mm
residuals:

valume fraction:

animation link:https://1drv.ms/u/s!AvvLa2FhdyA4tj_LbblYWguK1-ve
case 3: meshing 0.25mm also we have definesd a material
material properties shown below

residuals:

valume fraction:

animation link:https://1drv.ms/u/s!AvvLa2FhdyA4tj_LbblYWguK1-ve