Abstract
Methods are presented for integrating differential equations in conservation form in an Eulerian mesh in the presence of discontinuities and shocks. We make use of the monotonicity property, namely, that a solution is assumed to be monotone between mesh points. Monotonicity is preserved by convective equations and can be rigorously implemented in a code. In more general systems of conservative equations, monotonicity is enforced when interpolating a curve, given on an integer mesh, onto a half-integer mesh. This can be interpreted as introducing a local diffusive term whenever a discontinuity or shock occurs. The integration of a Riemann problem shows that with this method shocks are represented within two zones and contact discontinuities exhibit little diffusion.
Original language | English |
---|---|
Pages (from-to) | 212-226 |
Number of pages | 15 |
Journal | Journal of Computational Physics |
Volume | 38 |
Issue number | 2 |
DOIs | |
State | Published - 1 Jan 1980 |
Externally published | Yes |
ASJC Scopus subject areas
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- Physics and Astronomy (all)
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics