-
Design and analysis of an exactly divergence-free hybridized discontinuous Galerkin method for incompressible flows on meshes with quadrilateral cells
Authors:
Joseph P. Dean,
Sander Rhebergen,
Garth N. Wells
Abstract:
We generalise a hybridized discontinuous Galerkin method for incompressible flow problems to non-affine cells, showing that with a suitable element mapping the generalised method preserves a key invariance property that eludes most methods, namely that any irrotational component of the prescribed force is exactly balanced by the pressure gradient and does not affect the velocity field. This invari…
▽ More
We generalise a hybridized discontinuous Galerkin method for incompressible flow problems to non-affine cells, showing that with a suitable element mapping the generalised method preserves a key invariance property that eludes most methods, namely that any irrotational component of the prescribed force is exactly balanced by the pressure gradient and does not affect the velocity field. This invariance property can be preserved in the discrete problem if the incompressibility constraint is satisfied in a sufficiently strong sense. We derive sufficient conditions to guarantee discretely divergence-free functions are exactly divergence-free and give examples of divergence-free finite elements on meshes with triangular, quadrilateral, tetrahedral, or hexahedral cells generated by a (possibly non-affine) map from their respective reference cells. In the case of quadrilateral cells, we prove an optimal error estimate for the velocity field that does not depend on the pressure approximation. Our analysis is supported by numerical results.
△ Less
Submitted 26 September, 2023; v1 submitted 8 June, 2023;
originally announced June 2023.
-
Preconditioning for a pressure-robust HDG discretization of the Stokes equations
Authors:
Sander Rhebergen,
Garth N. Wells
Abstract:
We introduce a new preconditioner for a recently developed pressure-robust hybridized discontinuous Galerkin (HDG) finite element discretization of the Stokes equations. A feature of HDG methods is the straightforward elimination of degrees-of-freedom defined on the interior of an element. In our previous work (J. Sci. Comput., 77(3):1936--1952, 2018) we introduced a preconditioner for the case in…
▽ More
We introduce a new preconditioner for a recently developed pressure-robust hybridized discontinuous Galerkin (HDG) finite element discretization of the Stokes equations. A feature of HDG methods is the straightforward elimination of degrees-of-freedom defined on the interior of an element. In our previous work (J. Sci. Comput., 77(3):1936--1952, 2018) we introduced a preconditioner for the case in which only the degrees-of-freedom associated with the element velocity were eliminated via static condensation. In this work we introduce a preconditioner for the statically condensed system in which the element pressure degrees-of-freedom are also eliminated. In doing so the number of globally coupled degrees-of-freedom are reduced, but at the expense of a more difficult problem to analyse. We will show, however, that the Schur complement of the statically condensed system is spectrally equivalent to a simple trace pressure mass matrix. This result is used to formulate a new, provably optimal preconditioner. Through numerical examples in two- and three-dimensions we show that the new preconditioned iterative method converges in fewer iterations, has superior conservation properties for inexact solves, and is faster in CPU time when compared to our previous preconditioner.
△ Less
Submitted 18 October, 2021; v1 submitted 19 May, 2021;
originally announced May 2021.
-
An embedded--hybridized discontinuous Galerkin finite element method for the Stokes equations
Authors:
Sander Rhebergen,
Garth N. Wells
Abstract:
We present and analyze a new embedded--hybridized discontinuous Galerkin finite element method for the Stokes problem. The method has the attractive properties of full hybridized methods, namely an $H({\rm div})$-conforming velocity field, pointwise satisfaction of the continuity equation and \emph{a priori} error estimates for the velocity that are independent of the pressure. The embedded--hybri…
▽ More
We present and analyze a new embedded--hybridized discontinuous Galerkin finite element method for the Stokes problem. The method has the attractive properties of full hybridized methods, namely an $H({\rm div})$-conforming velocity field, pointwise satisfaction of the continuity equation and \emph{a priori} error estimates for the velocity that are independent of the pressure. The embedded--hybridized formulation has advantages over a full hybridized formulation in that it has fewer global degrees-of-freedom for a given mesh and the algebraic structure of the resulting linear system is better suited to fast iterative solvers. The analysis results are supported by a range of numerical examples that demonstrate rates of convergence, and which show computational efficiency gains over a full hybridized formulation.
△ Less
Submitted 18 September, 2019; v1 submitted 22 November, 2018;
originally announced November 2018.
-
A hybridizable discontinuous Galerkin method for the Navier--Stokes equations with pointwise divergence-free velocity field
Authors:
Sander Rhebergen,
Garth N. Wells
Abstract:
We introduce a hybridizable discontinuous Galerkin method for the incompressible Navier--Stokes equations for which the approximate velocity field is pointwise divergence-free. The method builds on the method presented by Labeur and Wells [SIAM J. Sci. Comput., vol. 34 (2012), pp. A889--A913]. We show that with modifications of the function spaces in the method of Labeur and Wells it is possible t…
▽ More
We introduce a hybridizable discontinuous Galerkin method for the incompressible Navier--Stokes equations for which the approximate velocity field is pointwise divergence-free. The method builds on the method presented by Labeur and Wells [SIAM J. Sci. Comput., vol. 34 (2012), pp. A889--A913]. We show that with modifications of the function spaces in the method of Labeur and Wells it is possible to formulate a simple method with pointwise divergence-free velocity fields which is momentum conserving, energy stable, and pressure-robust. Theoretical results are supported by two- and three-dimensional numerical examples and for different orders of polynomial approximation.
△ Less
Submitted 31 January, 2018; v1 submitted 25 April, 2017;
originally announced April 2017.
-
Three-field block-preconditioners for models of coupled magma/mantle dynamics
Authors:
Sander Rhebergen,
Garth N. Wells,
Andrew J. Wathen,
Richard F. Katz
Abstract:
For a prescribed porosity, the coupled magma/mantle flow equations can be formulated as a two-field system of equations with velocity and pressure as unknowns. Previous work has shown that while optimal preconditioners for the two-field formulation can be obtained, the construction of preconditioners that are uniform with respect to model parameters is difficult. This limits the applicability of t…
▽ More
For a prescribed porosity, the coupled magma/mantle flow equations can be formulated as a two-field system of equations with velocity and pressure as unknowns. Previous work has shown that while optimal preconditioners for the two-field formulation can be obtained, the construction of preconditioners that are uniform with respect to model parameters is difficult. This limits the applicability of two-field preconditioners in certain regimes of practical interest. We address this issue by reformulating the governing equations as a three-field problem, which removes a term that was problematic in the two-field formulation in favour of an additional equation for a pressure-like field. For the three-field problem, we develop and analyse new preconditioners and we show numerically that they are optimal in terms of problem size and less sensitive to model parameters, compared to the two-field preconditioner. This extends the applicability of optimal preconditioners for coupled mantle/magma dynamics into parameter regimes of physical interest.
△ Less
Submitted 4 June, 2015; v1 submitted 19 November, 2014;
originally announced November 2014.
-
Analysis of block-preconditioners for models of coupled magma/mantle dynamics
Authors:
Sander Rhebergen,
Garth N. Wells,
Richard F. Katz,
Andrew J. Wathen
Abstract:
This article considers the iterative solution of a finite element discretisation of the magma dynamics equations. In simplified form, the magma dynamics equations share some features of the Stokes equations. We therefore formulate, analyse and numerically test a Elman, Silvester and Wathen-type block preconditioner for magma dynamics. We prove analytically and demonstrate numerically the optimalit…
▽ More
This article considers the iterative solution of a finite element discretisation of the magma dynamics equations. In simplified form, the magma dynamics equations share some features of the Stokes equations. We therefore formulate, analyse and numerically test a Elman, Silvester and Wathen-type block preconditioner for magma dynamics. We prove analytically and demonstrate numerically the optimality of the preconditioner. The presented analysis highlights the dependence of the preconditioner on parameters in the magma dynamics equations that can affect convergence of iterative linear solvers. The analysis is verified through a range of two- and three-dimensional numerical examples on unstructured grids, from simple illustrative problems through to large problems on subduction zone-like geometries. The computer code to reproduce all numerical examples is freely available as supporting material.
△ Less
Submitted 29 May, 2014; v1 submitted 25 November, 2013;
originally announced November 2013.