A 3D cell-centered Lagrangian scheme applied to the simulation of 3D non-stationary Rayleigh–Taylor Instability in supernova remnants

Several astronomical flows can be studied thanks to the gas dynamics equations under the Lagrangian formalism. Here we propose to study the plerion dynamic, i.e. supernova remnants blown-up by a central pulsar as well as the Rayleigh–Taylor Instability (RTI) development at the inner interface of this gas shell.The scheme used is a multi-dimensional second order cell-centered Lagrangian scheme. It satisfies the Geometric Conservation Law (GCL), a semi-discrete entropy inequality and it conserves globally the momentum and the total energy. The convergence of the scheme towards the analytical solution is tested for the plerion test case and in the case of an axi-symmetric perturbation. Finally, the scheme is used to study the perturbation growth on the shell inner face perturbed with spherical harmonics.