English

Adaptive least-squares space-time finite element methods for convection-diffusion problems

Numerical Analysis 2025-09-16 v1 Numerical Analysis

Abstract

In this paper we formulate and analyse adaptive (space-time) least-squares finite element methods for the solution of convection-diffusion equations. The convective derivative vu\mathbf{v} \cdot \nabla u is considered as part of the total time derivative ddtu=tu+vu\frac{d}{dt}u = \partial_t u + \mathbf{v} \cdot \nabla u, and therefore we can use a rather standard stability and error analysis for related space-time finite element methods. For stationary problems we restrict the ansatz space H01(Ω)H^1_0(\Omega) such that the convective derivative is considered as an element of the dual H1(Ω)H^{-1}(\Omega) of the test space H01(Ω)H^1_0(\Omega), which also allows unbounded velocities v\mathbf{v}. While the discrete finite element schemes are always unique solvable, the numerical solutions may suffer from a bad approximation property of the finite element space when considering convection dominated problems, i.e., small diffusion coefficients. Instead of adding suitable stabilization terms, we aim to resolve the solutions by using adaptive (space-time) finite element methods. For this we introduce a least-squares approach where the discrete adjoint defines local a posteriori error indicators to drive an adaptive scheme. Numerical examples illustrate the theoretical considerations.

Keywords

Cite

@article{arxiv.2509.11955,
  title  = {Adaptive least-squares space-time finite element methods for convection-diffusion problems},
  author = {Christian Köthe and Olaf Steinbach},
  journal= {arXiv preprint arXiv:2509.11955},
  year   = {2025}
}
R2 v1 2026-07-01T05:36:56.336Z