-
Multigrid on unstructured meshes with regions of low quality cells
Authors:
Yuxuan Chen,
Garth N. Wells
Abstract:
The convergence of multigrid methods degrades significantly if a small number of low quality cells are present in a finite element mesh, and this can be a barrier to the efficient and robust application of multigrid on complicated geometric domains. The degraded performance is observed also if intermediate levels in a non-nested geometric multigrid problem have low quality cells, even when the fin…
▽ More
The convergence of multigrid methods degrades significantly if a small number of low quality cells are present in a finite element mesh, and this can be a barrier to the efficient and robust application of multigrid on complicated geometric domains. The degraded performance is observed also if intermediate levels in a non-nested geometric multigrid problem have low quality cells, even when the fine grid is high quality. It is demonstrated for geometric multigrid methods that the poor convergence is due to the local failure of smoothers to eliminate parts of error around cells of low quality. To overcome this, a global--local combined smoother is developed to maintain effective relaxation in the presence of a small number of poor quality cells. The smoother involves the application of a standard smoother on the whole domain, followed by local corrections for small subdomains with low quality cells. Two- and three-dimensional numerical experiments demonstrate that the degraded convergence of multigrid for low quality meshes can be restored to the high quality mesh reference case using the proposed smoother. The effect is particularly pronounced for higher-order finite elements. The results provide a basis for developing efficient, non-nested geometric multigrid methods for complicated engineering geometries.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
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.
-
Accelerating frequency-domain numerical methods for weakly nonlinear focused ultrasound using nested meshes
Authors:
Samuel P. Groth,
Pierre Gélat,
Seyyed R. Haqshenas,
Nader Saffari,
Elwin van 't Wout,
Timo Betcke,
Garth N. Wells
Abstract:
The numerical simulation of weakly nonlinear ultrasound is important in treatment planning for focused ultrasound (FUS) therapies. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough to resolve the highest harmonic present in the problem, leading to a v…
▽ More
The numerical simulation of weakly nonlinear ultrasound is important in treatment planning for focused ultrasound (FUS) therapies. However, the large domain sizes and generation of higher harmonics at the focus make these problems extremely computationally demanding. Numerical methods typically employ a uniform mesh fine enough to resolve the highest harmonic present in the problem, leading to a very large number of degrees of freedom. This paper proposes a more efficient strategy in which each harmonic is approximated on a separate mesh, the size of which is proportional to the wavelength of the harmonic. The increase in resolution required to resolve a smaller wavelength is balanced by a reduction in the domain size. This nested meshing is feasible owing to the increasingly localised nature of higher harmonics near the focus.
Numerical experiments are performed for FUS transducers in homogeneous media in order to determine the size of the meshes required to accurately represent the harmonics. In particular, a fast \emph{volume potential} approach is proposed and employed to perform convergence experiments as the computation domain size is modified. This approach allows each harmonic to be computed via the evaluation of an integral over the domain. Discretising this integral using the midpoint rule allows the computations to be performed rapidly with the FFT. It is shown that at least an order of magnitude reduction in memory consumption and computation time can be achieved with nested meshing. Finally, it is demonstrated how to generalise this approach to inhomogeneous propagation domains.
△ Less
Submitted 22 July, 2021; v1 submitted 5 November, 2020;
originally announced November 2020.
-
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.
-
Scalable computation of thermomechanical turbomachinery problems
Authors:
Chris N. Richardson,
Nathan Sime,
Garth N. Wells
Abstract:
A commonly held view in the turbomachinery community is that finite element methods are not well-suited for very large-scale thermomechanical simulations. We seek to dispel this notion by presenting performance data for a collection of realistic, large-scale thermomechanical simulations. We describe the necessary technology to compute problems with $O(10^7)$ to $O(10^9)$ degrees-of-freedom, and em…
▽ More
A commonly held view in the turbomachinery community is that finite element methods are not well-suited for very large-scale thermomechanical simulations. We seek to dispel this notion by presenting performance data for a collection of realistic, large-scale thermomechanical simulations. We describe the necessary technology to compute problems with $O(10^7)$ to $O(10^9)$ degrees-of-freedom, and emphasise what is required to achieve near linear computational complexity with good parallel scaling. Performance data is presented for turbomachinery components with up to 3.3 billion degrees-of-freedom. The software libraries used to perform the simulations are freely available under open source licenses. The performance demonstrated in this work opens up the possibility of system-level thermomechanical modelling, and lays the foundation for further research into high-performance formulations for even larger problems and for other physical processes, such as contact, that are important in turbomachinery analysis.
△ Less
Submitted 15 November, 2018; v1 submitted 26 April, 2018;
originally announced April 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.
-
Containers for portable, productive and performant scientific computing
Authors:
Jack S. Hale,
Lizao Li,
Chris N. Richardson,
Garth N. Wells
Abstract:
Containers are an emerging technology that hold promise for improving productivity and code portability in scientific computing. We examine Linux container technology for the distribution of a non-trivial scientific computing software stack and its execution on a spectrum of platforms from laptop computers through to high performance computing (HPC) systems. We show on a workstation and a leadersh…
▽ More
Containers are an emerging technology that hold promise for improving productivity and code portability in scientific computing. We examine Linux container technology for the distribution of a non-trivial scientific computing software stack and its execution on a spectrum of platforms from laptop computers through to high performance computing (HPC) systems. We show on a workstation and a leadership-class HPC system that when deployed appropriately there are no performance penalties running scientific programs inside containers. For Python code run on large parallel computers, the run time is reduced inside a container due to faster library imports. The software distribution approach and data that we present will help developers and users decide on whether container technology is appropriate for them. We also provide guidance for the vendors of HPC systems that rely on proprietary libraries for performance on what they can do to make containers work seamlessly and without performance penalty.
△ Less
Submitted 3 November, 2016; v1 submitted 26 August, 2016;
originally announced August 2016.
-
Automatic calibration of damping layers in finite element time domain simulations
Authors:
Steven Vandekerckhove,
Garth N. Wells,
Herbert De Gersem,
Koen Van Den Abeele
Abstract:
Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the domain of interest and the absorbing region. Selecting parameters in the attenuation functions is non-trivial. In this work, an optimisation procedure for autom…
▽ More
Matched layers are commonly used in numerical simulations of wave propagation to model (semi-)infinite domains. Attenuation functions describe the damping in layers, and provide a matching of the wave impedance at the interface between the domain of interest and the absorbing region. Selecting parameters in the attenuation functions is non-trivial. In this work, an optimisation procedure for automatically calibrating matched layers is presented. The procedure is based on solving optimisation problems constrained by partial differential equations with polynomial and piecewise-constant attenuation functions. We show experimentally that, for finite element time domain simulations, piecewise-constant attenuation function are at least as efficient as quadratic attenuation functions. This observation leads us to introduce consecutive matched layers as an alternative to perfectly matched layers, which can easily be employed for problems with arbitrary geometries. Moreover, the use of consecutive matched layers leads to a reduction in computational cost compared to perfectly matched layers. Examples are presented for acoustic, elastodynamic and electromagnetic problems. Numerical simulations are performed with the libraries FEniCS/DOLFIN and dolfin-adjoint, and the computer code to reproduce all numerical examples is made freely available.
△ Less
Submitted 28 January, 2016;
originally announced January 2016.
-
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.
-
Unified Form Language: A domain-specific language for weak formulations of partial differential equations
Authors:
Martin S. Alnaes,
Anders Logg,
Kristian B. Oelgaard,
Marie E. Rognes,
Garth N. Wells
Abstract:
We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential oper…
▽ More
We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.
△ Less
Submitted 25 April, 2013; v1 submitted 16 November, 2012;
originally announced November 2012.
-
Automated code generation for discontinuous Galerkin methods
Authors:
Kristian B. Ølgaard,
Anders Logg,
Garth N. Wells
Abstract:
A compiler approach for generating low-level computer code from high-level input for discontinuous Galerkin finite element forms is presented. The input language mirrors conventional mathematical notation, and the compiler generates efficient code in a standard programming language. This facilitates the rapid generation of efficient code for general equations in varying spatial dimensions. Key con…
▽ More
A compiler approach for generating low-level computer code from high-level input for discontinuous Galerkin finite element forms is presented. The input language mirrors conventional mathematical notation, and the compiler generates efficient code in a standard programming language. This facilitates the rapid generation of efficient code for general equations in varying spatial dimensions. Key concepts underlying the compiler approach and the automated generation of computer code are elaborated. The approach is demonstrated for a range of common problems, including the Poisson, biharmonic, advection--diffusion and Stokes equations.
△ Less
Submitted 4 April, 2011;
originally announced April 2011.
-
Optimisations for quadrature representations of finite element tensors through automated code generation
Authors:
Kristian B. Ølgaard,
Garth N. Wells
Abstract:
We examine aspects of the computation of finite element matrices and vectors which are made possible by automated code generation. Given a variational form in a syntax which resembles standard mathematical notation, the low-level computer code for building finite element tensors, typically matrices, vectors and scalars, can be generated automatically via a form compiler. In particular, the generat…
▽ More
We examine aspects of the computation of finite element matrices and vectors which are made possible by automated code generation. Given a variational form in a syntax which resembles standard mathematical notation, the low-level computer code for building finite element tensors, typically matrices, vectors and scalars, can be generated automatically via a form compiler. In particular, the generation of code for computing finite element matrices using a quadrature approach is addressed. For quadrature representations, a number of optimisation strategies which are made possible by automated code generation are presented. The relative performance of two different automatically generated representations of finite element matrices is examined, with a particular emphasis on complicated variational forms. It is shown that approaches which perform best for simple forms are not tractable for more complicated problems in terms of run time performance, the time required to generate the code or the size of the generated code. The approach and optimisations elaborated here are effective for a range of variational forms.
△ Less
Submitted 1 April, 2011;
originally announced April 2011.
-
DOLFIN: Automated Finite Element Computing
Authors:
Anders Logg,
Garth N. Wells
Abstract:
We describe here a library aimed at automating the solution of partial differential equations using the finite element method. By employing novel techniques for automated code generation, the library combines a high level of expressiveness with efficient computation. Finite element variational forms may be expressed in near mathematical notation, from which low-level code is automatically generate…
▽ More
We describe here a library aimed at automating the solution of partial differential equations using the finite element method. By employing novel techniques for automated code generation, the library combines a high level of expressiveness with efficient computation. Finite element variational forms may be expressed in near mathematical notation, from which low-level code is automatically generated, compiled and seamlessly integrated with efficient implementations of computational meshes and high-performance linear algebra. Easy-to-use object-oriented interfaces to the library are provided in the form of a C++ library and a Python module. This paper discusses the mathematical abstractions and methods used in the design of the library and its implementation. A number of examples are presented to demonstrate the use of the library in application code.
△ Less
Submitted 31 March, 2011;
originally announced March 2011.
-
A FEniCS-Based Programming Framework for Modeling Turbulent Flow by the Reynolds-Averaged Navier-Stokes Equations
Authors:
Mikael Mortensen,
Hans Petter Langtangen,
Garth N. Wells
Abstract:
Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by fi…
▽ More
Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model.
△ Less
Submitted 31 March, 2011; v1 submitted 14 February, 2011;
originally announced February 2011.
-
Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations
Authors:
Robert Jan Labeur,
Garth N. Wells
Abstract:
A hybrid method for the incompressible Navier--Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux conti…
▽ More
A hybrid method for the incompressible Navier--Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier--Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a point-wise solenoidal velocity field. Mass conservation, momentum conservation and global energy stability are proved for the time-continuous case, and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations.
△ Less
Submitted 15 January, 2012; v1 submitted 16 December, 2010;
originally announced December 2010.