A new stabilized linear finite element method for solving reaction-convection-diffusion equations

Po Wen Hsieh, Suh Yuh Yang

Research output: Contribution to journalArticlepeer-review

8 Scopus citations

Abstract

In this paper, we propose a new stabilized linear finite element method for solving reaction-convection-diffusion equations with arbitrary magnitudes of reaction and diffusion. The key feature of the new method is that the test function in the stabilization term is taken in the adjoint-operator-like form -εδv-(a·∇v)/γ+σv, where the parameter γ is appropriately designed to adjust the convection strength to achieve high accuracy and stability. We derive the stability estimates for the finite element solutions and establish the explicit dependence of L2 and H1 error bounds on the diffusivity, modulus of the convection field, reaction coefficient and the mesh size. The analysis shows that the proposed method is suitable for a wide range of mesh Péclet numbers and mesh Damköhler numbers. More specifically, if the diffusivity ε is sufficiently small with ε<∥a∥h and the reaction coefficient σ is large enough such that ∥a∥<σh, then the method exhibits optimal convergence rates in both L2 and H1 norms. However, for a small reaction coefficient satisfying ∥a∥≥σh, the method behaves like the well-known streamline upwind/Petrov-Galerkin formulation of Brooks and Hughes. Several numerical examples exhibiting boundary or interior layers are given to demonstrate the high performance of the proposed method. Moreover, we apply the developed method to time-dependent reaction-convection-diffusion problems and simulation results show the efficiency of the approach.

Original languageEnglish
Pages (from-to)362-382
Number of pages21
JournalComputer Methods in Applied Mechanics and Engineering
Volume307
DOIs
StatePublished - 1 Aug 2016

Keywords

  • Boundary layer
  • Interior layer
  • Reaction-convection-diffusion equation
  • Stabilization parameter
  • Stabilized finite element method

Fingerprint

Dive into the research topics of 'A new stabilized linear finite element method for solving reaction-convection-diffusion equations'. Together they form a unique fingerprint.

Cite this