Abstract:We present and analyze a new finite volume scheme of Gudonov-type for a nonlinear scalar conservation law whose flux function has a discontinuous coefficient due to time-dependent changes in its sign along a Lipschitz continuous curve.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to propose and analyze a new Godunov - type finite - volume scheme for solving nonlinear scalar conservation law equations with spatio - temporal flux discontinuities. Specifically, the paper focuses on the situation where the flux function changes sign on a Lipschitz continuous curve due to time - dependent changes.
### Problem Background
The considered conservation law equation is:
\[
\begin{cases}
\partial_t \rho+\partial_x F(t, x, \rho)=0, \\
\rho(0, x)=\rho_0(x), \quad \rho(t, \pm 1)=0,
\end{cases}
\]
where:
- \(t\in\mathbb{R}^+ = [0,\infty)\) is the time variable,
- \(x\in\Omega=(-1, 1)\) is the space variable,
- \((t, x)\mapsto\rho(t, x)\) is the unknown function, and \(\rho\in L^\infty(\mathbb{R}^+\times\Omega; [0, 1])\),
- The flux function \(F(t, x, \rho)\) is defined as:
\[
F(t, x, \rho)=\text{sign}(x-\xi(t)) f(\rho),
\]
where \(f: [0, 1]\to\mathbb{R}^+\) is a concave function, satisfying \(f(0) = 0=f(1)\), \(\max_{\rho\in[0, 1]} f(\rho)>0\), and \(\lim_{\rho\to0^+}\frac{f(\rho)}{\rho}\in\mathbb{R}\);
- \(\xi:\mathbb{R}^+\to\Omega\) is Lipschitz continuous, and \(\inf_{\mathbb{R}^+}\xi,\sup_{\mathbb{R}^+}\xi\subset(-1, 1)\).
### Main Challenges
This type of equations has applications in many problems in physical and engineering sciences, such as two - phase flow in porous media, traffic flow models, ion etching in semiconductor device manufacturing, etc. However, due to the discontinuity of the flux function, traditional numerical methods are difficult to handle these problems. In particular, for problems with spatio - temporal discontinuous fluxes, the classical Cauchy problem may not have a global classical solution, so the concept of weak solutions needs to be introduced.
### Solutions
The paper proposes a new Godunov - type finite - volume scheme, which deals with flux discontinuities by introducing a moving grid. The specific steps include:
1. **Define the moving grid**: In the discrete space - time domain, adjust the grid according to the position of the flux interface \(\xi\) so that \(\xi\) is included in the computational grid.
2. **Modify the numerical flux**: Use a modified numerical flux function based on the non - classical Riemann solver (\(R_\alpha\)) in the two cells near the flux interface, and use the standard Godunov flux formula in the places far from \(\xi\).
3. **Time evolution**: Ensure that in each time step, the change of the grid complies with the CFL condition to ensure numerical stability.
### Theoretical Contributions
The paper also discusses the existence of entropy - weak solutions and proves their convergence by constructing approximate solutions. In addition, the paper introduces the concept of "admissibility germs" to describe the entropy conditions on the flux - discontinuous interface, thus ensuring the uniqueness and well - posedness of the solution.
In summary, this paper aims to provide an effective method for dealing with scalar conservation law equations with spatio - temporal flux discontinuities, filling some gaps in the existing literature.