$L^\infty$ bounds for numerical solutions of noncoercive convection-diffusion equations

Claire Chainais-Hillairet,Maxime Herda
DOI: https://doi.org/10.1007/978-3-030-43651-3_12
2019-12-11
Abstract:In this work, we apply an iterative energy method à la de Giorgi in order to establish $L^{\infty}$ bounds for numerical solutions of noncoercive convection-diffusion equations with mixed Dirichlet-Neumann boundary conditions.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to establish \( L^\infty \) estimates for the numerical solutions of non - coercive convection - diffusion equations. Specifically, the author applies an iterative energy method (similar to the De Giorgi method) to prove that under the mixed Dirichlet - Neumann boundary conditions, the numerical solutions of non - coercive convection - diffusion equations have explicit \( L^\infty \) upper and lower bounds. ### Main problem description Consider the following convection - diffusion equation with mixed boundary conditions: \[ \begin{cases} \text{div}(-\nabla v + Uv)+bv = f & \text{in } \Omega, \\ (-\nabla v + Uv)\cdot n = 0 & \text{on } \Gamma_N, \\ v = v_D & \text{on } \Gamma_D, \end{cases} \] where: - \(\Omega\) is a bounded polygonal region, - \(\Gamma_D\) and \(\Gamma_N\) are the Dirichlet and Neumann boundaries respectively, - \(U\) is the velocity field, \(b\) is a non - negative function, \(f\) is the source term, and \(v_D\) is the boundary condition. ### Research background Although non - coercive elliptic linear problems have been widely studied, as far as the author knows, no literature has yet derived explicit \( L^\infty \) estimates. Therefore, the main contribution of this paper is to establish \( L^\infty \) estimates for numerical solutions through the iterative energy method. ### Numerical scheme To numerically solve the above problem, the author adopts a finite - volume scheme and assumes that the mesh satisfies certain regularity conditions. Specifically, the numerical flux is defined as: \[ F_{K,\sigma}= \begin{cases} 0, & \forall K\in T, \forall \sigma\in E^N_K, \\ \tau_\sigma\left(B(-U_{K,\sigma}d_\sigma)v_K - B(U_{K,\sigma}d_\sigma)v_{K,\sigma}\right), & \forall K\in T, \forall \sigma\in E_K\setminus E^N_K, \end{cases} \] where \(B(s)\) is a Lipschitz continuous function satisfying specific conditions. ### Main results Theorem 1 shows that under appropriate assumptions, the numerical solution \(v\) satisfies: \[ - M\leq v_K\leq M, \quad \forall K\in T, \] where \(M\) and \(M\) are non - negative constants depending only on \(\Omega\), \(\xi\), the function \(B\), \(\|U\|_{L^\infty}\), \(\|f^+\|_{L^\infty}\) and \(\|(v_D)^+\|_{L^\infty}\). ### Conclusion By introducing the De Giorgi iterative method, the author successfully establishes \( L^\infty \) estimates for the numerical solutions of non - coercive convection - diffusion equations. This result not only fills the gap in the existing literature but also provides a theoretical basis for further research.