Numerical Simulation of Rayleigh-Taylor Instability with Large Density Ratios Based on RKDG Method

Abstract

Rayleigh-Taylor instability problem with large density ratios is simulated by RKDG method which is developed for Euler equations with an additional body force corresponding to the gravity. The interface capturing ability of RKDG method is testified, while the density ratio (heavy to light) ranges from 3 to 20. Numerical results show that RKDG method has capability to pursue contact discontinuity in Rayleigh-Taylor instability with large density ratio. In the late stage of Rayleigh-Taylor instability problem, the contact line begins to crash, but the numerical solution is still smooth near the interface and has high resolution.