Flux-conservative Hermite methods for simulation of nonlinear conservation laws

Adeline Kornelus,Daniel Appelö
DOI: https://doi.org/10.48550/arXiv.1703.06848
2017-03-21
Abstract:A new class of Hermite methods for solving nonlinear conservation laws is presented. While preserving the high order spatial accuracy for smooth solutions in the existing Hermite methods, the new methods come with better stability properties. Artificial viscosity in the form of the entropy viscosity method is added to capture shocks.
Numerical Analysis
What problem does this paper attempt to address?
The problems that this paper attempts to solve are as follows: When simulating nonlinear conservation laws, existing methods face challenges in dealing with the co - existence of shock waves and small smooth scales. Specifically: 1. **Shock wave capture**: Shock waves often appear in the solutions of nonlinear conservation laws. These shock waves are characterized by rapid changes in the solution within a very thin region. Accurately approximating shock waves and smooth waves is a challenge because it is necessary to solve both small - scale and large - scale problems simultaneously. 2. **Numerical stability**: Although the existing Hermite method has high - order spatial accuracy in the case of smooth solutions, it has poor stability when dealing with shock waves. 3. **Discrete conservation**: The traditional Hermite method leads to discontinuities when calculating fluxes, thus losing the conservation property. This is especially obvious when dealing with nonlinear problems. To solve these problems, the author proposes a new class of flux - conservative Hermite methods. The main improvements of this method include: - **Better stability**: While maintaining high - order spatial accuracy, the new method improves stability, especially in capturing shock waves, by introducing artificial viscosity (entropy viscosity method). - **Flux continuity**: The new method ensures the continuity of the numerical flux at the cell interface by calculating the numerical flux at the cell boundary and interpolating it to the cell center, thus solving the discrete conservation problem in traditional methods. In summary, this paper aims to develop a new numerical method that can efficiently and stably simulate nonlinear conservation laws, especially for complex physical systems containing shock waves and small - scale smooth waves.