Abstract:In this paper, we introduce a framework for the discretization of a class of constrained Hamilton-Jacobi equations, a system coupling a Hamilton-Jacobi equation with a Lagrange multiplier determined by the constraint. The equation is non-local, and the constraint has bounded variations. We show that, under a set of general hypothesis, the approximation obtained with a finite-differences monotonic scheme, converges towards the viscosity solution of the constrained Hamilton-Jacobi equation.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of numerical approximation of a class of constrained Hamilton - Jacobi equations (HJ). Specifically, the paper proposes a framework to discretize such equations, which couple a Hamilton - Jacobi equation and Lagrange multipliers determined by constraint conditions. Such equations often appear in population models in quantitative genetics in the long - time and small - mutation asymptotic cases.
#### Main problems and background
1. **Mathematical model**:
- The research object is a class of constrained Hamilton - Jacobi equations:
\[
\begin{cases}
\partial_t u(t, x) + b(x, I(t)) H(\nabla_x u(t, x)) + R(t, x, I(t)) = 0, & x \in \mathbb{R}^d, t > 0 \\
\min_{x \in \mathbb{R}^d} u(t, x) = 0, & t \geq 0
\end{cases}
\]
where \( H(p) \) is a strictly convex Hamiltonian, \( b(x, I(t)) \) is the birth rate, \( R(t, x, I(t)) \) is the net growth rate, and \( I(t) \) is the weighted population size.
2. **Biological background**:
- Such equations are important tools for describing the long - term evolutionary behavior of populations with phenotypic characteristics, especially in the long - time and small - mutation asymptotic limits. For example, the Lotka - Volterra integral equation:
\[
\varepsilon \partial_t n^\varepsilon(t, x) - \frac{1}{\varepsilon} \int_{\mathbb{R}^d} K\left(\frac{y - x}{\varepsilon}\right) b(y, J^\varepsilon(t)) n^\varepsilon(t, y) dy - R(t, x, J^\varepsilon(t)) n^\varepsilon(t, x) = 0
\]
where \( J^\varepsilon(t)=\int_{\mathbb{R}^d} \psi(x) n^\varepsilon(t, x) dx \), representing the weighted population size.
3. **Numerical approximation**:
- The paper proposes a finite - difference monotone scheme to approximate the viscous solutions of the above equations and proves that this scheme converges to the viscous solutions under certain assumptions.
- In particular, the paper also constructs and analyzes an asymptotic - preserving (AP) scheme for handling the model stability in non - asymptotic cases.
4. **Main challenges**:
- Regularity problem of equation solutions: The solution \( u \) is expected to have only Lipschitz regularity, while \( I \) may have jump discontinuities.
- Selection of appropriate function spaces: To ensure the well - definition of the equation, it is necessary to select an appropriate function space \( W^{1,\infty}_{\text{loc}}(\mathbb{R}^d)\times BV_{\text{loc}} \).
- Selection of discretization parameters: It is necessary to select appropriate discretization parameters to ensure the stability and convergence of the numerical scheme.
#### Summary
By introducing the finite - difference monotone scheme and the asymptotic - preserving scheme, this paper solves the problem of numerical approximation of constrained Hamilton - Jacobi equations, which is especially suitable for describing population dynamics models in the long - time and small - mutation asymptotic limits.