TY - JOUR
T1 - Adaptive Osher-type scheme for the Euler equations with highly nonlinear equations of state
AU - Lee, Bok Jik
AU - Toro, Eleuterio F.
AU - Castro, Cristóbal E.
AU - Nikiforakis, Nikolaos
PY - 2013/8/1
Y1 - 2013/8/1
N2 - For the numerical simulation of detonation of condensed phase explosives, a complex equation of state (EOS), such as the Jones-Wilkins-Lee (JWL) EOS or the Cochran-Chan (C-C) EOS, are widely used. However, when a conservative scheme is used for solving the Euler equations with such equations of state, a spurious solution across the contact discontinuity, a well known phenomenon in multi-fluid systems, arises even for single materials. In this work, we develop a generalised Osher-type scheme in an adaptive primitive-conservative framework to overcome the aforementioned difficulties. Resulting numerical solutions are compared with the exact solutions and with the numerical solutions from the Godunov method in conjunction with the exact Riemann solver for the Euler equations with Mie-Grüneisen form of equations of state, such as the JWL and the C-C equations of state. The adaptive scheme is extended to second order and its empirical convergence rates are presented, verifying second order accuracy for smooth solutions. Through a suite of several tests problems in one and two space dimensions we illustrate the failure of conservative schemes and the capability of the methods of this paper to overcome the difficulties.
AB - For the numerical simulation of detonation of condensed phase explosives, a complex equation of state (EOS), such as the Jones-Wilkins-Lee (JWL) EOS or the Cochran-Chan (C-C) EOS, are widely used. However, when a conservative scheme is used for solving the Euler equations with such equations of state, a spurious solution across the contact discontinuity, a well known phenomenon in multi-fluid systems, arises even for single materials. In this work, we develop a generalised Osher-type scheme in an adaptive primitive-conservative framework to overcome the aforementioned difficulties. Resulting numerical solutions are compared with the exact solutions and with the numerical solutions from the Godunov method in conjunction with the exact Riemann solver for the Euler equations with Mie-Grüneisen form of equations of state, such as the JWL and the C-C equations of state. The adaptive scheme is extended to second order and its empirical convergence rates are presented, verifying second order accuracy for smooth solutions. Through a suite of several tests problems in one and two space dimensions we illustrate the failure of conservative schemes and the capability of the methods of this paper to overcome the difficulties.
KW - Dumbser-Osher-Toro solver
KW - Equation of state
KW - Euler equations
KW - Exact Riemann solver
KW - Godunov method
KW - Mie-Grüneisen
KW - Osher solver
KW - Primitive and conservative scheme
UR - https://www.scopus.com/pages/publications/84877315239
U2 - 10.1016/j.jcp.2013.03.046
DO - 10.1016/j.jcp.2013.03.046
M3 - Article
AN - SCOPUS:84877315239
SN - 0021-9991
VL - 246
SP - 165
EP - 183
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -