Finite element approximation of high-dimensional transport-dominated diffusion problems

High-dimensional partial differential equations with nonnegative characteristic form arise in numerous mathematical models in science. In problems of this kind, the computational challenge of beating the exponential growth of complexity as a function of dimension is exacerbated by the fact that the problem may be transport-dominated. We develop the analysis of stabilised sparse finite element methods for such high-dimensional, non-self-adjoint and possibly degenerate partial differential equations.\ud \ud (Presented as an invited lecture under the title "Computational multiscale modelling: Fokker-Planck equations and their numerical analysis" at the Foundations of Computational Mathematics conference in Santander, Spain, 30 June - 9 July, 2005.

An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior--penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi and Rebay, the standard SIPG method outlined in [Hartmann,Houston-2006], and an NIPG variant of the new scheme will be undertaken

hp-version discontinuous Galerkin methods for advection-diffusion-reaction problems on polytopic meshes

We consider the hp-version interior penalty discontinuous Galerkin finite element method (DGFEM) for the numerical approximation of the advection-diffusion-reaction equation on general computational meshes consisting of polygonal/polyhedral (polytopic) elements. In particular, new hp-version a priori error bounds are derived based on a specific choice of the interior penalty parameter which allows for edge/face-degeneration. The proposed method employs elemental polynomial bases of total degree p (P_p-basis) defined in the physical coordinate system, without requiring the mapping from a given reference or canonical frame. Numerical experiments highlighting the performance of the proposed DGFEM are presented. In particular, we study the competitiveness of the p-version DGFEM employing a P_p-basis on both polytopic and tensor-product elements with a (standard) DGFEM employing a (mapped) Q_p-basis. Moreover, a computational example is also presented which demonstrates the performance of the proposed hp-version DGFEM on general agglomerated meshes

Recovered finite element methods on polygonal and polyhedral meshes

Recovered Finite Element Methods (R-FEM) have been recently introduced in Georgoulis and Pryer [Comput. Methods Appl. Mech. Eng. 332 (2018) 303–324]. for meshes consisting of simplicial and/or box-type elements. Here, utilising the flexibility of the R-FEM framework, we extend their definition to polygonal and polyhedral meshes in two and three spatial dimensions, respectively. An attractive feature of this framework is its ability to produce arbitrary order polynomial conforming discretizations, yet involving only as many degrees of freedom as discontinuous Galerkin methods over general polygonal/polyhedral meshes with potentially many faces per element. A priori error bounds are shown for general linear, possibly degenerate, second order advection-diffusion-reaction boundary value problems. A series of numerical experiments highlight the good practical performance of the proposed numerical framework

An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken

An Optimal Order Interior Penalty Discontinuous Galerkin Discretization of the Compressible Navier-Stokes Equations

In this article we propose a new symmetric version of the interior penalty discontinuous Galerkin finite element method for the numerical approximation of the compressible Navier-Stokes equations. Here, particular emphasis is devoted to the construction of an optimal numerical method for the evaluation of certain target functionals of practical interest, such as the lift and drag coefficients of a body immersed in a viscous fluid. With this in mind, the key ingredients in the construction of the method include: (i) An adjoint consistent imposition of the boundary conditions; (ii) An adjoint consistent reformulation of the underlying target functional of practical interest; (iii) Design of appropriate interior-penalty stabilization terms. Numerical experiments presented within this article clearly indicate the optimality of the proposed method when the error is measured in terms of both the L_2-norm, as well as for certain target functionals. Computational comparisons with other discontinuous Galerkin schemes proposed in the literature, including the second scheme of Bassi & Rebay, cf. [11], the standard SIPG method outlined in [25], and an NIPG variant of the new scheme will be undertaken