A high-order discontinuous Galerkin method with Lagrange multipliers for advection-diffusion problems

Abstract A high-order Discontinuous Galerkin method with Lagrange Multipliers (DGLM) is presented for the solution of advection–diffusion problems on unstructured adaptive meshes. Unlike other hybrid discretization methods for transport problems, it operates directly on the second-order form the advection–diffusion equation. Like the Discontinuous Enrichment Method (DEM), it chooses the basis functions among the free-space solutions of the homogeneous form of the governing partial differential equation, and relies on Lagrange multipliers for enforcing a weak continuity of the approximated solution across the element interface boundaries. However unlike DEM, the proposed hybrid discontinuous Galerkin method approximates the Lagrange multipliers in a space of polynomials, instead of a space of traces on the element boundaries of the normal derivatives of a subset of the basis functions. For a homogeneous problem, the design of arbitrarily high-order elements based on this DGLM method is supported by a detailed mathematical analysis. For a non-homogeneous one, the approximated solution is locally decomposed into its homogeneous and particular parts. The homogeneous part is captured by the DGLM elements designed for a homogenous problem. The particular part is obtained analytically after the source term is projected onto an appropriate polynomial space. This decoupling between the two parts of the solution is another differentiator between DGLM and DEM with attractive computational advantages. An a posteriori error estimator for the proposed method is also derived to enable adaptive mesh refinement. All theoretical results are illustrated by numerical simulations that furthermore highlight the potential of the proposed high-order hybrid DG method for transport problems with steep gradients in the high Peclet number regime.

[1]  Song Wang,et al.  Multidimensional exponentially fitted simplicial finite elements for convection-diffusion equations with tensor-valued diffusion , 2005 .

[2]  C. Farhat,et al.  Bubble Functions Prompt Unusual Stabilized Finite Element Methods , 1994 .

[3]  Charbel Farhat,et al.  A discontinuous enrichment method for capturing evanescent waves in multiscale fluid and fluid/solid problems , 2008 .

[4]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[5]  Charbel Farhat,et al.  Unusual Stabilized Finite Element Methods and Residual-Free-Bubbles , 1996 .

[6]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[7]  P. Grisvard.  Elliptic Problems in Nonsmooth Domains , 1985 .

[8]  C. Farhat,et al.  The Discontinuous Enrichment Method , 2000 .

[9]  Charbel Farhat,et al.  A higher‐order discontinuous enrichment method for the solution of high péclet advection–diffusion problems on unstructured meshes , 2010 .

[10]  I. Babuska,et al.  The partition of unity finite element method: Basic theory and applications , 1996 .

[11]  Charbel Farhat,et al.  A discontinuous enrichment method for variable‐coefficient advection–diffusion at high Péclet number , 2011 .

[12]  Jim Douglas,et al.  An absolutely stabilized finite element method for the stokes problem , 1989 .

[13]  C. Farhat,et al.  The discontinuous enrichment method for elastic wave propagation in the medium‐frequency regime , 2006 .

[14]  T. Hughes.  Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[15]  Charbel Farhat,et al.  A discontinuous enrichment method for the finite element solution of high Péclet advection-diffusion problems , 2009 .

[16]  S. P. Oliveira,et al.  Streamline design of stability parameters for advection—diffusion problems , 2001 .

[17]  Michel Fortin,et al.  Variational Formulations and Finite Element Methods , 1991 .

[18]  Charbel Farhat,et al.  Higher‐order extensions of a discontinuous Galerkin method for mid‐frequency Helmholtz problems , 2004 .

[19]  Charbel Farhat,et al.  The discontinuous enrichment method for multiscale analysis , 2003 .

[20]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[21]  I. Babuska,et al.  The Partition of Unity Method , 1997 .

[22]  Alessandro Russo,et al.  CHOOSING BUBBLES FOR ADVECTION-DIFFUSION PROBLEMS , 1994 .

[23]  Charbel Farhat,et al.  Convergence Analysis of a Discontinuous Galerkin Method with Plane Waves and Lagrange Multipliers for the Solution of Helmholtz Problems , 2009, SIAM J. Numer. Anal..

[24]  Charbel Farhat,et al.  Residual-Free Bubbles for the Helmholtz Equation , 1996 .

[25]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

[26]  Charbel Farhat,et al.  A discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of short wave exterior Helmholtz problems on unstructured meshes , 2004 .

[27]  Thomas J. R. Hughes,et al.  Galerkin/least-squares finite element methods for the reduced wave equation with non-reflecting boundary conditions in unbounded domains , 1992 .

[28]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[29]  Charbel Farhat,et al.  A Variational Multiscale Method for the Large Eddy Simulation of Compressible Turbulent Flows on Unstructured Meshes - Application to vortex shedding , 2004 .

[30]  C. Farhat,et al.  Overview of the discontinuous enrichment method, the ultra‐weak variational formulation, and the partition of unity method for acoustic scattering in the medium frequency regime and performance comparisons , 2012 .

[31]  Charbel Farhat,et al.  Three‐dimensional discontinuous Galerkin elements with plane waves and Lagrange multipliers for the solution of mid‐frequency Helmholtz problems , 2006 .

[32]  Alessandro Russo,et al.  Applications of the pseudo residual-free bubbles to the stabilization of convection-diffusion problems , 1998 .