Rayleigh-Taylor Instability

Rayleigh-Taylor Instability

Numerical Simulation Initial & Boundary Condition

	Euler-equations for two ideal, polytropic gases (air and helium)
	Wave propagation scheme with Monotinized-Centered wave-limiter and transverse wave propagation, approximate Riemann-solver for multiple calorically perfect gases
	Quasi-steady wave propagation method used to incorporate gravitational force
	Transport equation for adiabatic coefficient used for pressure correction at material interface
	Computation with CFL-No. 0.95 to tend = 0.08s

	3D computation, one quarter of computational domain - AMR computation with a coarse grid of 20 x 20 x 60 cells. 2 levels with refinement factor 2 are used. Finest level corresponds to 80 x 80 x 240 grid (1.5M cells). - 1750 time steps calculated to tend - 80h real time on 16 nodes SP2 (nonadaptive 150 h)

	2D reference computation, cylindrical symmetric - Source terms due to symmetry treated with Strang splitting using a two-step Runge-Kutta method - Uniform 240 x 720 grid - 6870 time steps calculated to tend - 43h real time on 5 nodes SP2

A detailed paper including the computations is available:

Accurate Simulation of Rayleigh-Taylor instabilities
In P. Jonas and V. Uruba, editors, Proc. of Colloquium on Fluid Dynamics,
Institute of Thermodynamics, Academy of Sciences of Czech Republic, Prague, pages 37 - 44, Oct 19 - 20, 1999.

Postscript (gzipped, 568 kB)

Mpeg-Movies

Isosurface Plot in 3D (361k)

Reference computation in 2D, color plot of cross section (209k)

Reference computation in 2D, schlieren plot of cross cection in 2D (112k)

The computational domain is doubled for visualization.

Rayleigh-Taylor Instability

last update: 11/14/06