High order Lagrange-Galerkin methods for the conservative formulation of the advection-diffusion equation

Rodolfo Bermejo,Manuel Colera
2024-01-04
Abstract:We introduce in this paper the numerical analysis of high order both in time and space Lagrange-Galerkin methods for the conservative formulation of the advection-diffusion equation. As time discretization scheme we consider the Backward Differentiation Formulas up to order $q=5$. The development and analysis of the methods are performed in the framework of time evolving finite elements presented in C. M. Elliot and T. Ranner, IMA Journal of Numerical Analysis \textbf{41}, 1696-1845 (2021). The error estimates show through their dependence on the parameters of the equation the existence of different regimes in the behavior of the numerical solution; namely, in the diffusive regime, that is, when the diffusion parameter $\mu$ is large, the error is $O(h^{k+1}+\Delta t^{q})$, whereas in the advective regime, $\mu \ll 1$, the convergence is $O(\min (h^{k},\frac{h^{k+1} }{\Delta t})+\Delta t^{q})$. It is worth remarking that the error constant does not have exponential $\mu ^{-1}$ dependence.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop and analyze high - order Lagrange - Galerkin methods for solving the advection - diffusion equation in conservative form. Specifically, the paper focuses on how to achieve high - order accuracy simultaneously in time and space, and particularly considers the influence of the advection term and the diffusion term in the equation under different conditions. ### Main Problem Description 1. **Numerical Solution of the Advection - Diffusion Equation**: - The paper studies the numerical solution methods of the advection - diffusion equation, especially for its conservative form. The advection - diffusion equation describes the transport process of substances in fluids, which contains an advection term and a diffusion term. - The form of the equation is: \[ \begin{cases} \frac{\partial c}{\partial t}+\nabla\cdot(uc)-\mu\Delta c + a_0c = f(x,t)&\text{in }Q_T,\\ (\mu\nabla c)\cdot\nu = 0&\text{on }\partial\Omega\times(0,T),\\ c(x,0)=c_0(x)&\text{in }\Omega \end{cases} \] where \(c\) is the concentration function, \(u\) is the flow velocity vector, \(\mu\) is the diffusion coefficient, \(a_0\) is the reaction coefficient, and \(f\) is the source term. 2. **High - Order Time Discretization**: - The paper uses high - order backward difference formulas (BDF) for time discretization, with an accuracy of up to fifth - order (\(q = 5\)). - This high - order time discretization method can improve the accuracy of the numerical solution, especially in long - time simulations. 3. **Space Discretization**: - Space discretization adopts the finite element method and is analyzed within the framework of time - evolving finite elements. - In this way, regions with complex geometric shapes can be processed, especially when the boundary layer is very thin. 4. **Error Estimation**: - The paper analyzes in detail the error estimates under different parameter conditions (advection - dominated and diffusion - dominated). - In the diffusion - dominated case (\(\mu\) is large), the error is \(O(h^{k + 1}+\Delta t^q)\). - In the advection - dominated case (\(\mu\ll1\)), the error is \(O(\min(h^k,h^{k + 1}\Delta t)+\Delta t^q)\). - Importantly, the error constant does not depend on \(\mu^{- 1}\), which makes the method more reliable in advection - dominated problems. 5. **Mass Conservation Property**: - The paper emphasizes that the developed numerical method has the mass conservation property, that is, global mass balance. - This property is very important for many physical and engineering applications because it ensures the physical meaning of the numerical solution. ### Summary The main objective of this paper is to develop a high - order Lagrange - Galerkin method for solving the conservative form of the advection - diffusion equation. By combining high - order time discretization and finite - element space discretization, this method can provide high - accuracy numerical solutions under different parameter conditions and maintain the mass conservation property. This is of great significance for dealing with complex problems that are advection - dominated and diffusion - dominated.