Motivation
The Rayleigh-Taylor Instability (RTI) is a fundamental fluid instability occurring when a heavy fluid is supported by a light fluid against gravity. In astrophysics, this is critical for understanding:
- Supernova Remnants: The "fingers" seen in the Crab Nebula are classic RTI structures formed as pulsar wind nebulae accelerate ejecta.
- Solar Atmosphere: Magnetic flux tubes rising through the convection zone are subject to RTI, affecting solar flare stability.
Methodology
1. Linear Stability Analysis
We utilized Ideal Magneto-Hydrodynamics (MHD) to linearize the equations of motion. For a horizontal magnetic field $\mathbf{B} = B \hat{x}$, the dispersion relation is derived as:
Significance: The term with $B^2$ is negative, indicating stabilization. High wavenumbers $k$ (small structures) are stabilized first by magnetic tension.
2. Finite Volume Simulation
To study the non-linear "mushroom cap" evolution, we solved the compressible Euler equations using a Python-based Finite Volume Method (FVM).
Implications
Our results confirm that a critical magnetic field $B_c$ exists. If $B > B_c$, the instability is completely suppressed. This explains why certain astrophysical boundaries remain sharp despite strong density contrasts—magnetic tension acts as a "skin" holding the fluid interface together.
Detailed Project Report
Access the full detailed derivation of effect of magnetic field in supressing RTI including stability region plots and 2-D finite volume analysis code.
Read Full ReportNon-Linear Evolution (Simulation)
Finite Volume MethodSimulation showing the evolution of density. Note the formation of characteristic Kelvin-Helmholtz mushroom structures on the spikes as the instability grows.
Magnetic Field Stabilization
Growth Rate Squared ($\gamma^2$) vs Wavenumber ($k$) and Magnetic Field ($B$). The "valley" indicates stability.
Stability Threshold Surface