An Asymptotic Preserving and Energy Stable Scheme for the Euler System with Congestion Constraint

K.R. Arun,Amogh Krishnamurthy,Harihara Maharana

2024-06-20

Abstract:In this work, we design and analyze an asymptotic preserving (AP), semi-implicit finite volume scheme for the scaled compressible isentropic Euler system with a singular pressure law known as the congestion pressure law. The congestion pressure law imposes a maximal density constraint of the form $0\leq \varrho <1$, and the scaling introduces a small parameter $\varepsilon$ in order to control the stiffness of the density constraint. As $\varepsilon\to 0$, the solutions of the compressible system converge to solutions of the so-called free-congested Euler equations that couples compressible and incompressible dynamics. We show that the proposed scheme is positivity preserving and energy stable. In addition, we also show that the numerical densities satisfy a discrete variant of the constraint. By means of extensive numerical case studies, we verify the efficacy of the scheme and show that the scheme is able to capture the two dynamics in the limiting regime, thereby proving the AP property.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to design and analyze an Asymptotic Preserving (AP), semi - implicit finite - volume scheme for solving the compressible isentropic Euler system with a choking pressure law. Specifically, this problem involves the following aspects:

1. **System Description**:
- This research focuses on the one - dimensional isentropic Euler system with a special choking pressure law. The equations of this system are:
\[
\begin{cases}
\partial_t \varrho^\varepsilon + \partial_x (\varrho^\varepsilon u^\varepsilon) = 0, \\
\partial_t (\varrho^\varepsilon u^\varepsilon) + \partial_x (\varrho^\varepsilon (u^\varepsilon)^2) + \partial_x (p^\varepsilon(\varrho^\varepsilon)) = 0,
\end{cases}
\]
where $\varrho^\varepsilon$ and $u^\varepsilon$ represent density and velocity respectively, and the choking pressure $p^\varepsilon(\varrho^\varepsilon)$ is defined as:
\[
p^\varepsilon(\varrho^\varepsilon) = \frac{\varepsilon \varrho^\varepsilon}{(1 - \varrho^\varepsilon)^\gamma}, \quad \gamma\in(1,3].
\]

2. **Physical Background**:
- The choking phenomenon occurs in many macroscopic transport models, such as traffic flow, pedestrian flow, and concentration limitations in multiphase flow. These phenomena introduce the maximum density constraint $\varrho^\varepsilon < 1$, and when $\varepsilon\rightarrow0$, the solution of the system will converge to the so - called free - choked Euler equations, coupling compressible and incompressible dynamics.

3. **Numerical Challenges**:
- Design a numerical scheme that can accurately approximate the solution for $\varepsilon > 0$ and can capture two different dynamic behaviors (i.e., the free dynamics in the low - density region and the choked, incompressible dynamics in the high - density region) in the limit case of $\varepsilon\rightarrow0$.
- Ensure the stability of the numerical scheme when transitioning between different density regions, especially when the boundary between the low - density region and the high - density region changes.

4. **Objectives**:
- The proposed scheme needs to satisfy the following properties:
- Asymptotic Preserving (AP): When $\varepsilon\rightarrow0$, the scheme should be able to automatically adapt and stably approximate the limit system.
- Positivity - preserving: Ensure that the numerical density is always non - negative.
- Energy stability: Ensure that the total energy remains stable under a discrete version of the energy balance condition.

By designing this asymptotically - preserving and energy - stable numerical scheme, researchers hope to be able to effectively handle complex fluid - dynamic problems with a choking pressure law in numerical simulations, especially to provide reliable solutions in the transition region from free flow to choked flow.