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

Fuente: arXiv
Guardado en:
Detalles Bibliográficos
Autores principales: Bermejo, Rodolfo, Colera, Manuel
Formato: Preprint
Publicado: 2024
Materias:
Acceso en línea:
Etiquetas: Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
_version_ 1866916080935501824
author Bermejo, Rodolfo
Colera, Manuel
author_facet Bermejo, Rodolfo
Colera, Manuel
contents 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 $μ$ is large, the error is $O(h^{k+1}+Δt^{q})$, whereas in the advective regime, $μ\ll 1$, the convergence is $O(\min (h^{k},\frac{h^{k+1} }{Δt})+Δt^{q})$. It is worth remarking that the error constant does not have exponential $μ^{-1}$ dependence.
format Preprint
id arxiv_https___arxiv_org_abs_2401_02249
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle High order Lagrange-Galerkin methods for the conservative formulation of the advection-diffusion equation
Bermejo, Rodolfo
Colera, Manuel
Numerical Analysis
65M12, 65M25, 65M60, 65M50
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 $μ$ is large, the error is $O(h^{k+1}+Δt^{q})$, whereas in the advective regime, $μ\ll 1$, the convergence is $O(\min (h^{k},\frac{h^{k+1} }{Δt})+Δt^{q})$. It is worth remarking that the error constant does not have exponential $μ^{-1}$ dependence.
title High order Lagrange-Galerkin methods for the conservative formulation of the advection-diffusion equation
topic Numerical Analysis
65M12, 65M25, 65M60, 65M50
url https://arxiv.org/abs/2401.02249