Abstract:In this paper, we propose a new MUSCL scheme by combining the ideas of the Kurganov and Tadmor scheme and the so-called Deviation method which results in a well-balanced finite volume method for the hyperbolic balance laws, by evolving the difference between the exact solution and a given stationary solution. After that, we derive a semi-discrete scheme from this new scheme and it can be shown to be essentially TVD when applied to a scalar conservation law. In the end, we apply and validate the developed methods by numerical experiments and solve classical problems featuring Euler equations with gravitational source term.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a new numerical method for solving hyperbolic conservation laws with source terms (i.e., non - homogeneous conservation laws). Specifically, the author proposes a new MUSCL scheme. By combining the Kurganov and Tadmor scheme and the so - called Deviation method, an equilibrium scheme applicable to the Euler equations with gravitational source terms in one - dimensional space is constructed. This new method aims to preserve discrete steady - state solutions and can guarantee the total variation diminishing (TVD) property when applied to scalar conservation laws.
The core contributions of the paper are as follows:
1. **Propose a new central - upwind scheme**: This scheme is based on the KT scheme and the Deviation method and can preserve steady - state solutions without using a Riemann solver.
2. **Construct a semi - discrete scheme**: It is proved that this scheme has the TVD property when dealing with the modified homogeneous scalar conservation law.
3. **Conduct numerical experiment verification**: Through a series of classic test cases (such as isothermal equilibrium, isothermal equilibrium with perturbation, moving equilibrium and shock - tube problems), the effectiveness and accuracy of the proposed numerical method are verified.
### Summary of Mathematical Formulas
- **Conservation Laws and Source Terms**:
\[
\frac{\partial q(x, t)}{\partial t}+\nabla \cdot f(q(x, t)) = S(q(x, t)), \quad x\in\Omega\subset\mathbb{R}, \quad t > 0
\]
where \(q(x, t)\) is the vector of conserved variables, \(f(q)\) is the flux function, and \(S(q)\) is the source term.
- **Deviation Definition**:
\[
\Delta q=q - \tilde{q}
\]
where \(\tilde{q}\) is the given steady - state solution.
- **Modified Conservation Law**:
\[
(\Delta q)_t+[f(\Delta q+\tilde{q}) - f(\tilde{q})]_x = S(\Delta q+\tilde{q}, x)-S(\tilde{q}, x)
\]
- **Proof of TVD Property**:
\[
\frac{d}{dt}(\Delta q)_j(t)=-\frac{H_{j + 1/2}(t)-H_{j - 1/2}(t)}{\Delta x}+S_j(t)
\]
where \(H_{j+1/2}(t)\) is the numerical flux and \(S_j(t)\) is the source term.
Through these improvements, this method can accurately capture the steady - state solutions in physical phenomena in numerical simulations and avoid the numerical oscillation problems existing in traditional methods.