Space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting for multidimensional detonation waves simulation

TL;DR

This paper introduces a novel space-time adaptive ADER-DG finite element method with a specialized predictor and a posteriori sub-cell limiting, enabling accurate and efficient simulation of multidimensional detonation waves without splitting methods.

Contribution

The paper develops a new modification of the LST-DG predictor based on local time step partitioning, improving simulation of reactive flows with detonation waves.

Findings

Successfully simulates classical detonation problems with high accuracy.

Resolves detonation wave structures even on coarse meshes.

Avoids non-physical artifacts common in detonation simulations.

Abstract

The space-time adaptive ADER-DG finite element method with LST-DG predictor and a posteriori sub-cell ADER-WENO finite-volume limiting was used for simulation of multidimensional reacting flows with detonation waves. The presented numerical method does not use any ideas of splitting or fractional time steps methods. The modification of the LST-DG predictor has been developed, based on a local partition of the time step in cells in which strong reactivity of the medium is observed. This approach made it possible to obtain solutions to classical problems of flows with detonation waves and strong stiffness, without significantly decreasing the time step. The results obtained show the very high applicability and efficiency of using the ADER-DG-PN method with a posteriori sub-cell limiting for simulating reactive flows with detonation waves. The numerical solution shows the correct formation and propagation of ZND detonation waves. The structure of detonation waves is resolved by this numerical method with subcell resolution even on coarse spatial meshes. The smooth components of the numerical solution are correctly and very accurately reproduced by the numerical method. Non-physical artifacts of the numerical solution, typical for problems with detonation waves, such as the propagation of non-physical shock waves and weak detonation fronts ahead of the main detonation front, did not arise in the results obtained. The results of simulating rather complex problems associated with the propagation of detonation waves in significantly inhomogeneous domains are presented, which show that all the main features of detonation flows are correctly reproduced by this numerical method. It can be concluded that the space-time adaptive ADER-DG-PN method with a posteriori sub-cell ADER-WENO finite-volume limiting is perfectly applicable to simulating complex reacting flows with detonation waves.