Abstract
This paper addresses entropy generation and the corresponding dissipation of kinetic energy associated with high-resolution, shock-capturing (Godunov) methods. Analytical formulae are derived for the rate of increase of entropy given arbitrary jumps in primitive variables at a cell interface. It is demonstrated that for general continuously varying flows the inherent numerical entropy increase of Godunov methods is not proportional to the velocity jump cubed as is commonly assumed, but it is proportional to the velocity jump squared. Furthermore, the dissipation of kinetic energy is directly linked to temperature multiplied by change in entropy at low Mach numbers. The kinetic energy dissipation rate is shown to be proportional to the velocity jump squared and the speed of sound. The leading order dissipation rate associated with jumps in pressure, density and shear waves is detailed and further shown that at low Mach number it is the dissipation due to the perpendicular velocity jumps which dominates. This explains directly the poor performance of Godunov methods at low Mach numbers. The analysis is also applied to high-order accurate methods in space and time and all analytical results are validated with simple numerical experiments.
Original language | English |
---|---|
Pages (from-to) | 4853-4872 |
Number of pages | 20 |
Journal | Journal of Computational Physics |
Volume | 227 |
Issue number | 10 |
DOIs | |
Publication status | Published - 1 May 2008 |
Keywords
- high-resolution methods
- Godunov methods
- dissipation
- kinetic energy
- entropy
- large eddy simulation
- low Mach number