-
Stability, convergence, and pressure-robustness of numerical schemes for incompressible flows with hybrid velocity and pressure
Authors:
Lorenzo Botti,
Michele Botti,
Daniele Antonio Di Pietro,
Francesco Carlo Massa
Abstract:
In this work we study the stability, convergence, and pressure-robustness of discretization methods for incompressible flows with hybrid velocity and pressure. Specifically, focusing on the Stokes problem, we identify a set of assumptions that yield inf-sup stability as well as error estimates which distinguish the velocity- and pressure-related contributions to the error. We additionally identify…
▽ More
In this work we study the stability, convergence, and pressure-robustness of discretization methods for incompressible flows with hybrid velocity and pressure. Specifically, focusing on the Stokes problem, we identify a set of assumptions that yield inf-sup stability as well as error estimates which distinguish the velocity- and pressure-related contributions to the error. We additionally identify the key properties under which the pressure-related contributions vanish in the estimate of the velocity, thus leading to pressure-robustness. Several examples of existing and new schemes that fit into the framework are provided, and extensive numerical validation of the theoretical properties is provided.
△ Less
Submitted 19 April, 2024;
originally announced April 2024.
-
Homogeneous multigrid for hybrid discretizations: application to HHO methods
Authors:
Daniele A. Di Pietro,
Zhaonan Dong,
Guido Kanschat,
Pierre Matalon,
Andreas Rupp
Abstract:
We prove the uniform convergence of the geometric multigrid V-cycle for hybrid high-order (HHO) and other discontinuous skeletal methods. Our results generalize previously established results for HDG methods, and our multigrid method uses standard smoothers and local solvers that are bounded, convergent, and consistent. We use a weak version of elliptic regularity in our proofs. Numerical experime…
▽ More
We prove the uniform convergence of the geometric multigrid V-cycle for hybrid high-order (HHO) and other discontinuous skeletal methods. Our results generalize previously established results for HDG methods, and our multigrid method uses standard smoothers and local solvers that are bounded, convergent, and consistent. We use a weak version of elliptic regularity in our proofs. Numerical experiments confirm our theoretical results.
△ Less
Submitted 10 April, 2024; v1 submitted 23 March, 2024;
originally announced March 2024.
-
A Péclet-robust discontinuous Galerkin method for nonlinear diffusion with advection
Authors:
Lourenço Beirão da Veiga,
Daniele A. Di Pietro,
Kirubell B. Haile
Abstract:
We analyze a Discontinuous Galerkin method for a problem with linear advection-reaction and $p$-type diffusion, with Sobolev indices $p\in (1, \infty)$. The discretization of the diffusion term is based on the full gradient including jump liftings and interior-penalty stabilization while, for the advective contribution, we consider a strengthened version of the classical upwind scheme. The develop…
▽ More
We analyze a Discontinuous Galerkin method for a problem with linear advection-reaction and $p$-type diffusion, with Sobolev indices $p\in (1, \infty)$. The discretization of the diffusion term is based on the full gradient including jump liftings and interior-penalty stabilization while, for the advective contribution, we consider a strengthened version of the classical upwind scheme. The developed error estimates track the dependence of the local contributions to the error on local Péclet numbers. A set of numerical tests supports the theoretical derivations.
△ Less
Submitted 15 February, 2024;
originally announced February 2024.
-
A pressure-robust Discrete de Rham scheme for the Navier-Stokes equations
Authors:
Daniele A. Di Pietro,
Jerome Droniou,
Jia Jia Qian
Abstract:
In this work we design and analyse a Discrete de Rham (DDR) method for the incompressible Navier-Stokes equations. Our focus is, more specifically, on the SDDR variant, where a reduction in the number of unknowns is obtained using serendipity techniques. The main features of the DDR approach are the support of general meshes and arbitrary approximation orders. The method we develop is based on the…
▽ More
In this work we design and analyse a Discrete de Rham (DDR) method for the incompressible Navier-Stokes equations. Our focus is, more specifically, on the SDDR variant, where a reduction in the number of unknowns is obtained using serendipity techniques. The main features of the DDR approach are the support of general meshes and arbitrary approximation orders. The method we develop is based on the curl-curl formulation of the momentum equation and, through compatibility with the Helmholtz-Hodge decomposition, delivers pressure-robust error estimates for the velocity. It also enables non-standard boundary conditions, such as imposing the value of the pressure on the boundary. In-depth numerical validation on a complete panel of tests including general polyhedral meshes is provided. The paper also contains an appendix where bounds on DDR potential reconstructions and differential operators are proved in the more general framework of Polytopal Exterior Calculus.
△ Less
Submitted 9 January, 2024;
originally announced January 2024.
-
An a posteriori error analysis based on equilibrated stresses for finite element approximations of frictional contact
Authors:
Ilaria Fontana,
Daniele A. Di Pietro
Abstract:
We consider the unilateral contact problem between an elastic body and a rigid foundation in a description that includes both Tresca and Coulomb friction conditions. For this problem, we present an a posteriori error analysis based on an equilibrated stress reconstruction in the Arnold--Falk--Winther space that includes a guaranteed upper bound distinguishing the different components of the error.…
▽ More
We consider the unilateral contact problem between an elastic body and a rigid foundation in a description that includes both Tresca and Coulomb friction conditions. For this problem, we present an a posteriori error analysis based on an equilibrated stress reconstruction in the Arnold--Falk--Winther space that includes a guaranteed upper bound distinguishing the different components of the error. This analysis is the starting point for the development of an adaptive algorithm including a stopping criterion for the generalized Newton method. This algorithm is then used to perform numerical simulations that validate the theoretical results.
△ Less
Submitted 5 January, 2024;
originally announced January 2024.
-
Uniform Poincaré inequalities for the Discrete de Rham complex on general domains
Authors:
Daniele A. Di Pietro,
Marien-Lorenzo Hanot
Abstract:
In this paper we prove Poincaré inequalities for the Discrete de Rham (DDR) sequence on a general connected polyhedral domain $Ω$ of $\mathbb{R}^3$. We unify the ideas behind the inequalities for all three operators in the sequence, deriving new proofs for the Poincaré inequalities for the gradient and the divergence, and extending the available Poincaré inequality for the curl to domains with arb…
▽ More
In this paper we prove Poincaré inequalities for the Discrete de Rham (DDR) sequence on a general connected polyhedral domain $Ω$ of $\mathbb{R}^3$. We unify the ideas behind the inequalities for all three operators in the sequence, deriving new proofs for the Poincaré inequalities for the gradient and the divergence, and extending the available Poincaré inequality for the curl to domains with arbitrary second Betti numbers. A key preliminary step consists in deriving "mimetic" Poincaré inequalities giving the existence and stability of the solutions to topological balance problems useful in general discrete geometric settings. As an example of application, we study the stability of a novel DDR scheme for the magnetostatics problem on domains with general topology.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
A discrete three-dimensional divdiv complex on polyhedral meshes with application to a mixed formulation of the biharmonic problem
Authors:
Daniele A. Di Pietro,
Marien-Lorenzo Hanot
Abstract:
In this work, following the Discrete de Rham (DDR) paradigm, we develop an arbitrary-order discrete divdiv complex on general polyhedral meshes. The construction rests 1) on discrete spaces that are spanned by vectors of polynomials whose components are attached to mesh entities and 2) on discrete operators obtained mimicking integration by parts formulas. We provide an in-depth study of the algeb…
▽ More
In this work, following the Discrete de Rham (DDR) paradigm, we develop an arbitrary-order discrete divdiv complex on general polyhedral meshes. The construction rests 1) on discrete spaces that are spanned by vectors of polynomials whose components are attached to mesh entities and 2) on discrete operators obtained mimicking integration by parts formulas. We provide an in-depth study of the algebraic properties of the local complex, showing that it is exact on mesh elements with trivial topology. The new DDR complex is used to design a numerical scheme for the approximation of biharmonic problems, for which we provide detailed stability and convergence analyses.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
An exterior calculus framework for polytopal methods
Authors:
Francesco Bonaldi,
Daniele A. Di Pietro,
Jerome Droniou,
Kaibo Hu
Abstract:
We develop in this work the first polytopal complexes of differential forms. These complexes, inspired by the Discrete De Rham and the Virtual Element approaches, are discrete versions of the de Rham complex of differential forms built on meshes made of general polytopal elements. Both constructions benefit from the high-level approach of polytopal methods, which leads, on certain meshes, to leane…
▽ More
We develop in this work the first polytopal complexes of differential forms. These complexes, inspired by the Discrete De Rham and the Virtual Element approaches, are discrete versions of the de Rham complex of differential forms built on meshes made of general polytopal elements. Both constructions benefit from the high-level approach of polytopal methods, which leads, on certain meshes, to leaner constructions than the finite element method. We establish commutation properties between the interpolators and the discrete and continuous exterior derivatives, prove key polynomial consistency results for the complexes, and show that their cohomologies are isomorphic to the cohomology of the continuous de Rham complex.
△ Less
Submitted 6 July, 2023; v1 submitted 20 March, 2023;
originally announced March 2023.
-
A polytopal method for the Brinkman problem robust in all regimes
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work we develop a discretisation method for the Brinkman problem that is uniformly well-behaved in all regimes (as identified by a local dimensionless number with the meaning of a friction coefficient) and supports general meshes as well as arbitrary approximation orders. The method is obtained combining ideas from the Hybrid High-Order and Discrete de Rham methods, and its robustness rest…
▽ More
In this work we develop a discretisation method for the Brinkman problem that is uniformly well-behaved in all regimes (as identified by a local dimensionless number with the meaning of a friction coefficient) and supports general meshes as well as arbitrary approximation orders. The method is obtained combining ideas from the Hybrid High-Order and Discrete de Rham methods, and its robustness rests on a potential reconstruction and stabilisation terms that change in nature according to the value of the local friction coefficient. We derive error estimates that, thanks to the presence of cut-off factors, are valid across the all regimes and provide extensive numerical validation.
△ Less
Submitted 4 March, 2023; v1 submitted 9 January, 2023;
originally announced January 2023.
-
An arbitrary-order discrete rot-rot complex on polygonal meshes with application to a quad-rot problem
Authors:
Daniele A. Di Pietro
Abstract:
In this work, following the discrete de Rham (DDR) approach, we develop a discrete counterpart of a two-dimensional de Rham complex with enhanced regularity. The proposed construction supports general polygonal meshes and arbitrary approximation orders. We establish exactness on a contractible domain for both the versions of the complex with and without boundary conditions and, for the former, pro…
▽ More
In this work, following the discrete de Rham (DDR) approach, we develop a discrete counterpart of a two-dimensional de Rham complex with enhanced regularity. The proposed construction supports general polygonal meshes and arbitrary approximation orders. We establish exactness on a contractible domain for both the versions of the complex with and without boundary conditions and, for the former, prove a complete set of Poincaré-type inequalities. The discrete complex is then used to derive a novel discretisation method for a quad-rot problem which, unlike other schemes in the literature, does not require the forcing term to be prepared. We carry out complete stability and convergence analyses for the proposed scheme and provide numerical validation of the results.
△ Less
Submitted 27 October, 2022;
originally announced October 2022.
-
Cohomology of the discrete de Rham complex on domains of general topology
Authors:
Daniele A. Di Pietro,
Jérôme Droniou,
Silvano Pitassi
Abstract:
In this work we prove that, for a general polyhedral domain of $\mathbb{R}^3$, the cohomology spaces of the discrete de Rham complex of [Di Pietro and Droniou, An arbitrary-order discrete de Rham complex on polyhedral meshes: Exactness, Poincaré inequalities, and consistency, Found. Comput. Math., 2021, DOI: 10.1007/s10208-021-09542-8] are isomorphic to those of the continuous de Rham complex. Thi…
▽ More
In this work we prove that, for a general polyhedral domain of $\mathbb{R}^3$, the cohomology spaces of the discrete de Rham complex of [Di Pietro and Droniou, An arbitrary-order discrete de Rham complex on polyhedral meshes: Exactness, Poincaré inequalities, and consistency, Found. Comput. Math., 2021, DOI: 10.1007/s10208-021-09542-8] are isomorphic to those of the continuous de Rham complex. This is, to the best of our knowledge, the first result of this kind for an arbitrary-order complex built from a general polyhedral mesh.
△ Less
Submitted 24 May, 2023; v1 submitted 2 September, 2022;
originally announced September 2022.
-
A serendipity fully discrete div-div complex on polygonal meshes
Authors:
Michele Botti,
Daniele A. Di Pietro,
Marwa Salah
Abstract:
In this work we address the reduction of face degrees of freedom (DOFs) for discrete elasticity complexes. Specifically, using serendipity techniques, we develop a reduced version of a recently introduced two-dimensional complex arising from traces of the three-dimensional elasticity complex. The keystone of the reduction process is a new estimate of symmetric tensor-valued polynomial fields in te…
▽ More
In this work we address the reduction of face degrees of freedom (DOFs) for discrete elasticity complexes. Specifically, using serendipity techniques, we develop a reduced version of a recently introduced two-dimensional complex arising from traces of the three-dimensional elasticity complex. The keystone of the reduction process is a new estimate of symmetric tensor-valued polynomial fields in terms of boundary values, completed with suitable projections of internal values for higher degrees. We prove an extensive set of new results for the original complex and show that the reduced complex has the same homological and analytical properties as the original one. This paper also contains an appendix with proofs of general Poincaré--Korn-type inequalities for hybrid fields.
△ Less
Submitted 14 July, 2022;
originally announced July 2022.
-
A pressure-robust HHO method for the solution of the incompressible Navier-Stokes equations on general meshes
Authors:
Daniel Castanon Quiroz,
Daniele A. Di Pietro
Abstract:
In a recent work [10], we have introduced a pressure-robust Hybrid High-Order method for the numerical solution of the incompressible Navier-Stokes equations on matching simplicial meshes. Pressure-robust methods are characterized by error estimates for the velocity that are fully independent of the pressure. A crucial question was left open in that work, namely whether the proposed construction c…
▽ More
In a recent work [10], we have introduced a pressure-robust Hybrid High-Order method for the numerical solution of the incompressible Navier-Stokes equations on matching simplicial meshes. Pressure-robust methods are characterized by error estimates for the velocity that are fully independent of the pressure. A crucial question was left open in that work, namely whether the proposed construction could be extended to general polytopal meshes. In this paper we provide a positive answer to this question. Specifically, we introduce a novel divergence-preserving velocity reconstruction that hinges on the solution inside each element of a mixed problem on a subtriangulation, then use it to design discretizations of the body force and convective terms that lead to pressure robustness. An in-depth theoretical study of the properties of this velocity reconstruction, and their reverberation on the scheme, is carried out for polynomial degrees $k \geq 0$ and meshes composed of general polytopes. The theoretical convergence estimates and the pressure robustness of the method are confirmed by an extensive panel of numerical examples.
△ Less
Submitted 11 October, 2022; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Homological- and analytical-preserving serendipity framework for polytopal complexes, with application to the DDR method
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work we investigate from a broad perspective the reduction of degrees of freedom through serendipity techniques for polytopal methods compatible with Hilbert complexes. We first establish an abstract framework that, given two complexes connected by graded maps, identifies a set of properties enabling the transfer of the homological and analytical properties from one complex to the other. T…
▽ More
In this work we investigate from a broad perspective the reduction of degrees of freedom through serendipity techniques for polytopal methods compatible with Hilbert complexes. We first establish an abstract framework that, given two complexes connected by graded maps, identifies a set of properties enabling the transfer of the homological and analytical properties from one complex to the other. This abstract framework is designed having in mind discrete complexes, with one of them being a reduced version of the other, such as occurring when applying serendipity techniques to numerical methods. We then use this framework as an overarching blueprint to design a serendipity DDR complex. Thanks to the combined use of higher-order reconstructions and serendipity, this complex compares favorably in terms of degrees of freedom (DOF) count to all the other polytopal methods previously introduced and also to finite elements on certain element geometries. The gain resulting from such a reduction in the number of DOFs is numerically evaluated on two model problems: a magnetostatic model, and the Stokes equations.
△ Less
Submitted 28 September, 2022; v1 submitted 6 March, 2022;
originally announced March 2022.
-
A fully discrete plates complex on polygonal meshes with application to the Kirchhoff-Love problem
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work we develop a novel fully discrete version of the plates complex, an exact Hilbert complex relevant for the mixed formulation of fourth-order problems. The derivation of the discrete complex follows the discrete de Rham paradigm, leading to an arbitrary-order construction that applies to meshes composed of general polygonal elements. The discrete plates complex is then used to derive a…
▽ More
In this work we develop a novel fully discrete version of the plates complex, an exact Hilbert complex relevant for the mixed formulation of fourth-order problems. The derivation of the discrete complex follows the discrete de Rham paradigm, leading to an arbitrary-order construction that applies to meshes composed of general polygonal elements. The discrete plates complex is then used to derive a novel numerical scheme for Kirchhoff--Love plates, for which a full stability and convergence analysis are performed. Extensive numerical tests complete the exposition.
△ Less
Submitted 29 December, 2021;
originally announced December 2021.
-
Arbitrary-order pressure-robust DDR and VEM methods for the Stokes problem on polyhedral meshes
Authors:
Lourenço Beirão da Veiga,
Franco Dassi,
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
This paper contains two major contributions. First we derive, following the discrete de Rham (DDR) and Virtual Element (VEM) paradigms, pressure-robust methods for the Stokes equations that support arbitrary orders and polyhedral meshes. Unlike other methods presented in the literature, pressure-robustness is achieved here without resorting to an $\boldsymbol{H}({\rm div})$-conforming construction…
▽ More
This paper contains two major contributions. First we derive, following the discrete de Rham (DDR) and Virtual Element (VEM) paradigms, pressure-robust methods for the Stokes equations that support arbitrary orders and polyhedral meshes. Unlike other methods presented in the literature, pressure-robustness is achieved here without resorting to an $\boldsymbol{H}({\rm div})$-conforming construction on a submesh, but rather projecting the volumetric force onto the discrete $\boldsymbol{H}({\bf curl})$ space. The cancellation of the pressure error contribution stems from key commutation properties of the underlying DDR and VEM complexes. The pressure-robust error estimates in $h^{k+1}$ (with $h$ denoting the meshsize and $k\ge 0$ the polynomial degree of the DDR or VEM complex) are proven theoretically and supported by a panel of three-dimensional numerical tests. The second major contribution of the paper is an in-depth study of the relations between the DDR and VEM approaches. We show, in particular, that a complex developed following one paradigm admits a reformulation in the other, and that couples of related DDR and VEM complexes satisfy commuting diagram properties with the degrees of freedom maps.
△ Less
Submitted 27 June, 2024; v1 submitted 17 December, 2021;
originally announced December 2021.
-
A posteriori error estimates via equilibrated stress reconstructions for contact problems approximated by Nitsche's method
Authors:
Daniele Antonio Di Pietro,
Ilaria Fontana,
Kyrylo Kazymyrenko
Abstract:
We present an a posteriori error estimate based on equilibrated stress reconstructions for the finite element approximation of a unilateral contact problem with weak enforcement of the contact conditions. We start by proving a guaranteed upper bound for the dual norm of the residual. This norm is shown to control the natural energy norm up to a boundary term, which can be removed under a saturatio…
▽ More
We present an a posteriori error estimate based on equilibrated stress reconstructions for the finite element approximation of a unilateral contact problem with weak enforcement of the contact conditions. We start by proving a guaranteed upper bound for the dual norm of the residual. This norm is shown to control the natural energy norm up to a boundary term, which can be removed under a saturation assumption. The basic estimate is then refined to distinguish the different components of the error, and is used as a starting point to design an algorithm including adaptive stopping criteria for the nonlinear solver and automatic tuning of a regularization parameter. We then discuss an actual way of computing the stress reconstruction based on the Arnold-Falk-Winther finite elements. Finally, after briefly discussing the efficiency of our estimators, we showcase their performance on a panel of numerical tests.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
A Hybrid High-Order method for incompressible flows of non-Newtonian fluids with power-like convective behaviour
Authors:
Daniel Castanon Quiroz,
Daniele Antonio Di Pietro,
André Harnist
Abstract:
In this work, we design and analyze a Hybrid High-Order (HHO) discretization method for incompressible flows of non-Newtonian fluids with power-like convective behaviour. We work under general assumptions on the viscosity and convection laws, that are associated with possibly different Sobolev exponents r > 1 and s > 1. After providing a novel weak formulation of the continuous problem, we study i…
▽ More
In this work, we design and analyze a Hybrid High-Order (HHO) discretization method for incompressible flows of non-Newtonian fluids with power-like convective behaviour. We work under general assumptions on the viscosity and convection laws, that are associated with possibly different Sobolev exponents r > 1 and s > 1. After providing a novel weak formulation of the continuous problem, we study its well-posedness highlighting how a subtle interplay between the exponents r and s determines the existence and uniqueness of a solution. We next design an HHO scheme based on this weak formulation and perform a comprehensive stability and convergence analysis, including convergence for general data and error estimates for shear-thinning fluids and small data. The HHO scheme is validated on a complete panel of model problems.
△ Less
Submitted 9 January, 2022; v1 submitted 28 June, 2021;
originally announced June 2021.
-
A discrete de Rham method for the Reissner-Mindlin plate bending problem on polygonal meshes
Authors:
Daniele A. Di Pietro,
Jerome Droniou
Abstract:
In this work we propose a discretisation method for the Reissner--Mindlin plate bending problem in primitive variables that supports general polygonal meshes and arbitrary order. The method is inspired by a two-dimensional discrete de Rham complex for which key commutation properties hold that enable the cancellation of the contribution to the error linked to the enforcement of the Kirchhoff const…
▽ More
In this work we propose a discretisation method for the Reissner--Mindlin plate bending problem in primitive variables that supports general polygonal meshes and arbitrary order. The method is inspired by a two-dimensional discrete de Rham complex for which key commutation properties hold that enable the cancellation of the contribution to the error linked to the enforcement of the Kirchhoff constraint. Denoting by $k\ge 0$ the polynomial degree for the discrete spaces and by $h$ the meshsize, we derive for the proposed method an error estimate in $h^{k+1}$ for general $k$, as well as a locking-free error estimate for the lowest-order case $k=0$. The theoretical results are validated on a complete panel of numerical tests.
△ Less
Submitted 2 September, 2022; v1 submitted 25 May, 2021;
originally announced May 2021.
-
An arbitrary-order discrete de Rham complex on polyhedral meshes. Part II: Consistency
Authors:
Daniele Antonio Di Pietro,
Jérôme Droniou
Abstract:
In this paper we prove a complete panel of consistency results for the discrete de Rham (DDR) complex introduced in the companion paper [D. A. Di Pietro and J. Droniou, An arbitrary-order discrete de Rham complex on polyhedral meshes. Part I: Exactness and Poincaré inequalities, 2021, submitted], including primal and adjoint consistency for the discrete vector calculus operators, and consistency o…
▽ More
In this paper we prove a complete panel of consistency results for the discrete de Rham (DDR) complex introduced in the companion paper [D. A. Di Pietro and J. Droniou, An arbitrary-order discrete de Rham complex on polyhedral meshes. Part I: Exactness and Poincaré inequalities, 2021, submitted], including primal and adjoint consistency for the discrete vector calculus operators, and consistency of the corresponding potentials. The theoretical results are showcased by performing a full convergence analysis for a DDR approximation of a magnetostatics model. Numerical results on three-dimensional polyhedral meshes complete the exposition.
△ Less
Submitted 1 February, 2021; v1 submitted 13 January, 2021;
originally announced January 2021.
-
An arbitrary-order discrete de Rham complex on polyhedral meshes: Exactness, Poincaré inequalities, and consistency
Authors:
Daniele Antonio Di Pietro,
Jérôme Droniou
Abstract:
In this paper we present a novel arbitrary-order discrete de Rham (DDR) complex on general polyhedral meshes based on the decomposition of polynomial spaces into ranges of vector calculus operators and complements linked to the spaces in the Koszul complex. The DDR complex is fully discrete, meaning that both the spaces and discrete calculus operators are replaced by discrete counterparts, and sat…
▽ More
In this paper we present a novel arbitrary-order discrete de Rham (DDR) complex on general polyhedral meshes based on the decomposition of polynomial spaces into ranges of vector calculus operators and complements linked to the spaces in the Koszul complex. The DDR complex is fully discrete, meaning that both the spaces and discrete calculus operators are replaced by discrete counterparts, and satisfies suitable exactness properties depending on the topology of the domain. In conjunction with bespoke discrete counterparts of $L^2$-products, it can be used to design schemes for partial differential equations that benefit from the exactness of the sequence but, unlike classical (e.g., Raviart--Thomas--Nédélec) finite elements, are nonconforming. We prove a complete panel of results for the analysis of such schemes: exactness properties, uniform Poincaré inequalities, as well as primal and adjoint consistency. We also show how this DDR complex enables the design of a numerical scheme for a magnetostatics problem, and use the aforementioned results to prove stability and optimal error estimates for this scheme.
△ Less
Submitted 3 November, 2021; v1 submitted 13 January, 2021;
originally announced January 2021.
-
Improved error estimates for Hybrid High-Order discretizations of Leray-Lions problems
Authors:
Daniele Antonio Di Pietro,
Jérôme Droniou,
André Harnist
Abstract:
We derive novel error estimates for Hybrid High-Order (HHO) discretizations of Leray-Lions problems set in W^(1,p) with p in (1,2]. Specifically, we prove that, depending on the degeneracy of the problem, the convergence rate may vary between (k+1)(p-1) and (k+1), with k denoting the degree of the HHO approximation. These regime-dependent error estimates are illustrated by a complete panel of nume…
▽ More
We derive novel error estimates for Hybrid High-Order (HHO) discretizations of Leray-Lions problems set in W^(1,p) with p in (1,2]. Specifically, we prove that, depending on the degeneracy of the problem, the convergence rate may vary between (k+1)(p-1) and (k+1), with k denoting the degree of the HHO approximation. These regime-dependent error estimates are illustrated by a complete panel of numerical experiments.
△ Less
Submitted 9 January, 2022; v1 submitted 9 December, 2020;
originally announced December 2020.
-
$p$-Multilevel preconditioners for HHO discretizations of the Stokes equations with static condensation
Authors:
Lorenzo Botti,
Daniele Antonio Di Pietro
Abstract:
We propose a $p$-multilevel preconditioner for Hybrid High-Order discretizations (HHO) of the Stokes equation, numerically assess its performance on two variants of the method, and compare with a classical Discontinuous Galerkin scheme. We specifically investigate how the combination of $p$-coarsening and static condensation influences the performance of the $V$-cycle iteration for HHO. Two differ…
▽ More
We propose a $p$-multilevel preconditioner for Hybrid High-Order discretizations (HHO) of the Stokes equation, numerically assess its performance on two variants of the method, and compare with a classical Discontinuous Galerkin scheme. We specifically investigate how the combination of $p$-coarsening and static condensation influences the performance of the $V$-cycle iteration for HHO. Two different static condensation procedures are considered, resulting in global linear systems with a different number of unknowns and non-zero elements. An efficient implementation is proposed where coarse level operators are inherited using $L^2$-orthogonal projections defined over mesh faces and the restriction of the fine grid operators is performed recursively and matrix-free. The various resolution strategies are thoroughly validated on two- and three-dimensional problems.
△ Less
Submitted 29 September, 2020;
originally announced September 2020.
-
A discrete Weber inequality on three-dimensional hybrid spaces with application to the HHO approximation of magnetostatics
Authors:
Florent Chave,
Daniele A. Di Pietro,
Simon Lemaire
Abstract:
We prove a discrete version of the first Weber inequality on three-dimensional hybrid spaces spanned by vectors of polynomials attached to the elements and faces of a polyhedral mesh. We then introduce two Hybrid High-Order methods for the approximation of the magnetostatics model, in both its (first-order) field and (second-order) vector potential formulations. These methods are applicable on gen…
▽ More
We prove a discrete version of the first Weber inequality on three-dimensional hybrid spaces spanned by vectors of polynomials attached to the elements and faces of a polyhedral mesh. We then introduce two Hybrid High-Order methods for the approximation of the magnetostatics model, in both its (first-order) field and (second-order) vector potential formulations. These methods are applicable on general polyhedral meshes, and allow for arbitrary orders of approximation. Leveraging the previously established discrete Weber inequality, we perform a comprehensive analysis of the two methods. We finally validate them on a set of test-cases.
△ Less
Submitted 25 October, 2021; v1 submitted 7 July, 2020;
originally announced July 2020.
-
An arbitrary-order method for magnetostatics on polyhedral meshes based on a discrete de Rham sequence
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work, we develop a discretisation method for the mixed formulation of the magnetostatic problem supporting arbitrary orders and polyhedral meshes. The method is based on a global discrete de Rham (DDR) sequence, obtained by patching the local spaces constructed in [Di Pietro, Droniou, Rapetti, Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra, arXiv:…
▽ More
In this work, we develop a discretisation method for the mixed formulation of the magnetostatic problem supporting arbitrary orders and polyhedral meshes. The method is based on a global discrete de Rham (DDR) sequence, obtained by patching the local spaces constructed in [Di Pietro, Droniou, Rapetti, Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra, arXiv:1911.03616] by enforcing the single-valuedness of the components attached to the boundary of each element. The first main contribution of this paper is a proof of exactness relations for this global DDR sequence, obtained leveraging the exactness of the corresponding local sequence and a topological assembly of the mesh valid for domains that do not enclose any void. The second main contribution is the formulation and well-posedness analysis of the method, which includes the proof of uniform Poincaré inequalities for the discrete divergence and curl operators. The convergence rate in the natural energy norm is numerically evaluated on standard and polyhedral meshes. When the DDR sequence of degree $k\ge 0$ is used, the error converges as $h^{k+1}$, with $h$ denoting the meshsize.
△ Less
Submitted 11 November, 2020; v1 submitted 14 May, 2020;
originally announced May 2020.
-
A Hybrid High-Order method for creeping flows of non-Newtonian fluids
Authors:
Michele Botti,
Daniel Castanon Quiroz,
Daniele A. Di Pietro,
André Harnist
Abstract:
In this paper, we design and analyze a Hybrid High-Order discretization method for the steady motion of non-Newtonian, incompressible fluids in the Stokes approximation of small velocities. The proposed method has several appealing features including the support of general meshes and high-order, unconditional inf-sup stability, and orders of convergence that match those obtained for scalar Leray-L…
▽ More
In this paper, we design and analyze a Hybrid High-Order discretization method for the steady motion of non-Newtonian, incompressible fluids in the Stokes approximation of small velocities. The proposed method has several appealing features including the support of general meshes and high-order, unconditional inf-sup stability, and orders of convergence that match those obtained for scalar Leray-Lions problems. A complete well-posedness and convergence analysis of the method is carried out under new, general assumptions on the strain rate-shear stress law, which encompass several common examples such as the power-law and Carreau-Yasuda models. Numerical examples complete the exposition.
△ Less
Submitted 9 January, 2022; v1 submitted 30 March, 2020;
originally announced March 2020.
-
An abstract analysis framework for monolithic discretisations of poroelasticity with application to Hybrid High-Order methods
Authors:
Lorenzo Botti,
Michele Botti,
Daniele A. Di Pietro
Abstract:
In this work, we introduce a novel abstract framework for the stability and convergence analysis of fully coupled discretisations of the poroelasticity problem and apply it to the analysis of Hybrid High-Order (HHO) schemes. A relevant feature of the proposed framework is that it rests on mild time regularity assumptions that can be derived from an appropriate weak formulation of the continuous pr…
▽ More
In this work, we introduce a novel abstract framework for the stability and convergence analysis of fully coupled discretisations of the poroelasticity problem and apply it to the analysis of Hybrid High-Order (HHO) schemes. A relevant feature of the proposed framework is that it rests on mild time regularity assumptions that can be derived from an appropriate weak formulation of the continuous problem. To the best of our knowledge, these regularity results for the Biot problem are new. A novel family of HHO discretisation schemes is proposed and analysed, and their performance numerically evaluated.
△ Less
Submitted 8 December, 2019;
originally announced December 2019.
-
Fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra
Authors:
Daniele A. Di Pietro,
Jérôme Droniou,
Francesca Rapetti
Abstract:
In this work, merging ideas from compatible discretisations and polyhedral methods, we construct novel fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra. The spaces and operators that appear in these sequences are directly amenable to computer implementation. Besides proving exactness, we show that the usual three-dimensional sequence of trimmed Finite Eleme…
▽ More
In this work, merging ideas from compatible discretisations and polyhedral methods, we construct novel fully discrete polynomial de Rham sequences of arbitrary degree on polygons and polyhedra. The spaces and operators that appear in these sequences are directly amenable to computer implementation. Besides proving exactness, we show that the usual three-dimensional sequence of trimmed Finite Element spaces forms, through appropriate interpolation operators, a commutative diagram with our sequence, which ensures suitable approximation properties. A discussion on reconstructions of potentials and discrete $L^2$-products completes the exposition.
△ Less
Submitted 17 May, 2021; v1 submitted 9 November, 2019;
originally announced November 2019.
-
A high-order discretization of nonlinear poroelasticity
Authors:
Michele Botti,
Daniele A. Di Pietro,
Pierre Sochala
Abstract:
In this work we construct and analyze a nonconforming high-order discretization method for the quasi-static single-phase nonlinear poroelasticity problem describing Darcean flow in a deformable porous medium saturated by a slightly compressible fluid. The nonlinear elasticity operator is discretized using a Hybrid High-Order method, while the Darcy operator relies on a Symmetric Weighted Interior…
▽ More
In this work we construct and analyze a nonconforming high-order discretization method for the quasi-static single-phase nonlinear poroelasticity problem describing Darcean flow in a deformable porous medium saturated by a slightly compressible fluid. The nonlinear elasticity operator is discretized using a Hybrid High-Order method, while the Darcy operator relies on a Symmetric Weighted Interior Penalty discontinuous Galerkin scheme. The method is valid in two and three space dimensions, delivers an inf-sup stable discretization on general meshes including polyhedral elements and nonmatching interfaces, supports arbitrary approximation orders, and has a reduced cost thanks to the possibility of statically condensing a large subset of the unknowns for linearized versions of the problem. Moreover, the proposed construction can handle both nonzero and vanishing specific storage coefficients.
△ Less
Submitted 3 June, 2019;
originally announced June 2019.
-
Numerical approximation of poroelasticity with random coefficients using Polynomial Chaos and Hybrid High-Order methods
Authors:
Michele Botti,
Daniele A. Di Pietro,
Olivier Le Maître,
Pierre Sochala
Abstract:
In this work, we consider the Biot problem with uncertain poroelastic coefficients. The uncertainty is modelled using a finite set of parameters with prescribed probability distribution. We present the variational formulation of the stochastic partial differential system and establish its well-posedness. We then discuss the approximation of the parameter-dependent problem by non-intrusive techniqu…
▽ More
In this work, we consider the Biot problem with uncertain poroelastic coefficients. The uncertainty is modelled using a finite set of parameters with prescribed probability distribution. We present the variational formulation of the stochastic partial differential system and establish its well-posedness. We then discuss the approximation of the parameter-dependent problem by non-intrusive techniques based on Polynomial Chaos decompositions. We specifically focus on sparse spectral projection methods, which essentially amount to performing an ensemble of deterministic model simulations to estimate the expansion coefficients. The deterministic solver is based on a Hybrid High-Order discretization supporting general polyhedral meshes and arbitrary approximation orders. We numerically investigate the convergence of the probability error of the Polynomial Chaos approximation with respect to the level of the sparse grid. Finally, we assess the propagation of the input uncertainty onto the solution considering an injection-extraction problem.
△ Less
Submitted 11 November, 2019; v1 submitted 28 March, 2019;
originally announced March 2019.
-
A low-order nonconforming method for linear elasticity on general meshes
Authors:
Michele Botti,
Daniele A. Di Pietro,
Alessandra Guglielmana
Abstract:
In this work we construct a low-order nonconforming approximation method for linear elasticity problems supporting general meshes and valid in two and three space dimensions. The method is obtained by hacking the Hybrid High-Order method, that requires the use of polynomials of degree $k\ge1$ for stability. Specifically, we show that coercivity can be recovered for $k=0$ by introducing a novel ter…
▽ More
In this work we construct a low-order nonconforming approximation method for linear elasticity problems supporting general meshes and valid in two and three space dimensions. The method is obtained by hacking the Hybrid High-Order method, that requires the use of polynomials of degree $k\ge1$ for stability. Specifically, we show that coercivity can be recovered for $k=0$ by introducing a novel term that penalises the jumps of the displacement reconstruction across mesh faces. This term plays a key role in the fulfillment of a discrete Korn inequality on broken polynomial spaces, for which a novel proof valid for general polyhedral meshes is provided. Locking-free error estimates are derived for both the energy- and the $L^2$-norms of the error, that are shown to convergence, for smooth solutions, as $h$ and $h^2$, respectively (here, $h$ denotes the meshsize). A thorough numerical validation on a complete panel of two- and three-dimensional test cases is provided.
△ Less
Submitted 6 February, 2019;
originally announced February 2019.
-
A third Strang lemma and an Aubin-Nitsche trick for schemes in fully discrete formulation
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work, we present an abstract error analysis framework for the approximation of linear partial differential equation (PDE) problems in weak formulation. We consider approximation methods in fully discrete formulation, where the discrete and continuous spaces are possibly not embedded in a common space. A proper notion of consistency is designed, and, under a classical inf-sup condition, it…
▽ More
In this work, we present an abstract error analysis framework for the approximation of linear partial differential equation (PDE) problems in weak formulation. We consider approximation methods in fully discrete formulation, where the discrete and continuous spaces are possibly not embedded in a common space. A proper notion of consistency is designed, and, under a classical inf-sup condition, it is shown to bound the approximation error. This error estimate result is in the spirit of Strang's first and second lemmas, but applicable in situations not covered by these lemmas (because of a fully discrete approximation space). An improved estimate is also established in a weaker norm, using the Aubin--Nitsche trick.
We then apply these abstract estimates to an anisotropic heterogeneous diffusion model and two classical families of schemes for this model: Virtual Element and Finite Volume methods. For each of these methods, we show that the abstract results yield new error estimates with a precise and mild dependency on the local anisotropy ratio. A key intermediate step to derive such estimates for Virtual Element Methods is proving optimal approximation properties of the oblique elliptic projector in weighted Sobolev seminorms. This is a result whose interest goes beyond the specific model and methods considered here. We also obtain, to our knowledge, the first clear notion of consistency for Finite Volume methods, which leads to a generic error estimate involving the fluxes and valid for a wide range of Finite Volume schemes. An important application is the first error estimate for Multi-Point Flux Approximation L and G methods.
In the appendix, not included in the published version of this work, we show that classical estimates for discontinuous Galerkin methods can be obtained with simplified arguments using the abstract framework.
△ Less
Submitted 13 November, 2018; v1 submitted 25 April, 2018;
originally announced April 2018.
-
A Hybrid High-Order discretisation of the Brinkman problem robust in the Darcy and Stokes limits
Authors:
Lorenzo Botti,
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work, we develop and analyse a novel Hybrid High-Order discretisation of the Brinkman problem. The method hinges on hybrid discrete velocity unknowns at faces and elements and on discontinuous pressures. Based on the discrete unknowns, we reconstruct inside each element a Stokes velocity one degree higher than face unknowns, and a Darcy velocity in the Raviart-Thomas-Nédélec space. These r…
▽ More
In this work, we develop and analyse a novel Hybrid High-Order discretisation of the Brinkman problem. The method hinges on hybrid discrete velocity unknowns at faces and elements and on discontinuous pressures. Based on the discrete unknowns, we reconstruct inside each element a Stokes velocity one degree higher than face unknowns, and a Darcy velocity in the Raviart-Thomas-Nédélec space. These reconstructed velocities are respectively used to formulate the discrete versions of the Stokes and Darcy terms in the momentum equation, along with suitably designed penalty contributions. The proposed construction is tailored to yield optimal error estimates that are robust throughout the entire spectrum of local (Stokes- or Darcy-dominated) regimes, as identified by a dimensionless number which can be interpreted as a friction coefficient. The singular limit corresponding to the Darcy equation is also fully supported by the method. Numerical examples corroborate the theoretical results. This paper also contains two contributions whose interest goes beyond the specific method and application treated in this work: an investigation of the dependence of the constant in the second Korn inequality on star-shaped domains and its application to the study of the approximation properties of the strain projector in general Sobolev seminorms.
△ Less
Submitted 23 July, 2018; v1 submitted 29 March, 2018;
originally announced March 2018.
-
An advection-robust Hybrid High-Order method for the Oseen problem
Authors:
Joubine Aghili,
Daniele A. Di Pietro
Abstract:
In this work, we study advection-robust Hybrid High-Order discretizations of the Oseen equations. For a given integer $k\ge 0$, the discrete velocity unknowns are vector-valued polynomials of total degree $\le k$ on mesh elements and faces, while the pressure unknowns are discontinuous polynomials of total degree $\le k$ on the mesh. From the discrete unknowns, three relevant quantities are recons…
▽ More
In this work, we study advection-robust Hybrid High-Order discretizations of the Oseen equations. For a given integer $k\ge 0$, the discrete velocity unknowns are vector-valued polynomials of total degree $\le k$ on mesh elements and faces, while the pressure unknowns are discontinuous polynomials of total degree $\le k$ on the mesh. From the discrete unknowns, three relevant quantities are reconstructed inside each element: a velocity of total degree $\le(k+1)$, a discrete advective derivative, and a discrete divergence. These reconstructions are used to formulate the discretizations of the viscous, advective, and velocity-pressure coupling terms, respectively. Well-posedness is ensured through appropriate high-order stabilization terms. We prove energy error estimates that are advection-robust for the velocity, and show that each mesh element $T$ of diameter $h_T$ contributes to the discretization error with an $\mathcal{O}(h_T^{k+1})$-term in the diffusion-dominated regime, an $\mathcal{O}(h_T^{k+\frac12})$-term in the advection-dominated regime, and scales with intermediate powers of $h_T$ in between. Numerical results complete the exposition.
△ Less
Submitted 18 February, 2018; v1 submitted 7 December, 2017;
originally announced December 2017.
-
Discontinuous Skeletal Gradient Discretisation Methods on polytopal meshes
Authors:
Daniele A. Di Pietro,
Jérôme Droniou,
Gianmarco Manzini
Abstract:
In this work we develop arbitrary-order Discontinuous Skeletal Gradient Discretisations (DSGD) on general polytopal meshes. Discontinuous Skeletal refers to the fact that the globally coupled unknowns are broken polynomial on the mesh skeleton. The key ingredient is a high-order gradient reconstruction composed of two terms: (i) a consistent contribution obtained mimicking an integration by parts…
▽ More
In this work we develop arbitrary-order Discontinuous Skeletal Gradient Discretisations (DSGD) on general polytopal meshes. Discontinuous Skeletal refers to the fact that the globally coupled unknowns are broken polynomial on the mesh skeleton. The key ingredient is a high-order gradient reconstruction composed of two terms: (i) a consistent contribution obtained mimicking an integration by parts formula inside each element and (ii) a stabilising term for which sufficient design conditions are provided. An example of stabilisation that satisfies the design conditions is proposed based on a local lifting of high-order residuals on a Raviart-Thomas-Nédélec subspace. We prove that the novel DSGDs satisfy coercivity, consistency, limit-conformity, and compactness requirements that ensure convergence for a variety of elliptic and parabolic problems. Links with Hybrid High-Order, non-conforming Mimetic Finite Difference and non-conforming Virtual Element methods are also studied. Numerical examples complete the exposition.
△ Less
Submitted 6 December, 2017; v1 submitted 29 June, 2017;
originally announced June 2017.
-
A Hybrid High-Order method for Kirchhoff-Love plate bending problems
Authors:
Francesco Bonaldi,
Daniele A. Di Pietro,
Giuseppe Geymonat,
Françoise Krasucki
Abstract:
We present a novel Hybrid High-Order (HHO) discretization of fourth-order elliptic problems arising from the mechanical modeling of the bending behavior of Kirchhoff-Love plates, including the biharmonic equation as a particular case. The proposed HHO method supports arbitrary approximation orders on general polygonal meshes, and reproduces the key mechanical equilibrium relations locally inside e…
▽ More
We present a novel Hybrid High-Order (HHO) discretization of fourth-order elliptic problems arising from the mechanical modeling of the bending behavior of Kirchhoff-Love plates, including the biharmonic equation as a particular case. The proposed HHO method supports arbitrary approximation orders on general polygonal meshes, and reproduces the key mechanical equilibrium relations locally inside each element. When polynomials of degree $k \ge 1$ are used as unknowns, we prove convergence in $h^{k+1}$ (with $h$ denoting, as usual, the meshsize) in an energy-like norm. A key ingredient in the proof are novel approximation results for the energy projector on local polynomial spaces. Under biharmonic regularity assumptions, a sharp estimate in $h^{k+3}$ is also derived for the $L^2$-norm of the error on the deflection. The theoretical results are supported by numerical experiments, which additionally show the robustness of the method with respect to the choice of the stabilization.
△ Less
Submitted 24 January, 2018; v1 submitted 21 June, 2017;
originally announced June 2017.
-
An introduction to Hybrid High-Order methods
Authors:
Daniele A. Di Pietro,
Roberta Tittarelli
Abstract:
This chapter provides an introduction to Hybrid High-Order (HHO) methods. These are new generation numerical methods for PDEs with several advantageous features: the support of arbitrary approximation orders on general polyhedral meshes, the reproduction at the discrete level of relevant continuous properties, and a reduced computational cost thanks to static condensation and compact stencil. Afte…
▽ More
This chapter provides an introduction to Hybrid High-Order (HHO) methods. These are new generation numerical methods for PDEs with several advantageous features: the support of arbitrary approximation orders on general polyhedral meshes, the reproduction at the discrete level of relevant continuous properties, and a reduced computational cost thanks to static condensation and compact stencil. After establishing the discrete setting, we introduce the basics of HHO methods using as a model problem the Poisson equation. We describe in detail the construction, and prove a priori convergence results for various norms of the error as well as a posteriori estimates for the energy norm. We then consider two applications: the discretization of the nonlinear $p$-Laplace equation and of scalar diffusion-advection-reaction problems. The former application is used to introduce compactness analysis techniques to study the convergence to minimal regularity solution. The latter is used to introduce the discretization of first-order operators and the weak enforcement of boundary conditions. Numerical examples accompany the exposition.
△ Less
Submitted 20 April, 2017; v1 submitted 15 March, 2017;
originally announced March 2017.
-
A Hybrid High-Order method for the steady incompressible Navier--Stokes problem
Authors:
Daniele A. Di Pietro,
Stella Krell
Abstract:
In this work we introduce and analyze a novel Hybrid High-Order method for the steady incompressible Navier-Stokes equations. The proposed method is inf-sup stable on general polyhedral meshes, supports arbitrary approximation orders, and is (relatively) inexpensive thanks to the possibility of statically condensing a subset of the unknowns at each nonlinear iteration. We show under general assump…
▽ More
In this work we introduce and analyze a novel Hybrid High-Order method for the steady incompressible Navier-Stokes equations. The proposed method is inf-sup stable on general polyhedral meshes, supports arbitrary approximation orders, and is (relatively) inexpensive thanks to the possibility of statically condensing a subset of the unknowns at each nonlinear iteration. We show under general assumptions the existence of a discrete solution, which is also unique provided a data smallness condition is verified. Using a compactness argument, we prove convergence of the sequence of discrete solutions to minimal regularity exact solutions for general data. For more regular solutions, we prove optimal convergence rates for the energy-norm of the velocity and the $L^2$-norm of the pressure under a standard data smallness assumption. More precisely, when polynomials of degree $k\ge 0$ at mesh elements and faces are used, both quantities are proved to converge as $h^{k+1}$ (with $h$ denoting the meshsize).
△ Less
Submitted 27 July, 2016;
originally announced July 2016.
-
A Hybrid High-Order method for the Cahn-Hilliard problem in mixed form
Authors:
Florent Chave,
Daniele A. Di Pietro,
Fabien Marche,
Franck Pigeonneau
Abstract:
In this work, we develop a fully implicit Hybrid High-Order algorithm for the Cahn-Hilliard problem in mixed form. The space discretization hinges on local reconstruction operators from hybrid polynomial unknowns at elements and faces. The proposed method has several assets: (i) It supports fairly general meshes possibly containing polygonal elements and nonmatching interfaces, (ii) it allows arbi…
▽ More
In this work, we develop a fully implicit Hybrid High-Order algorithm for the Cahn-Hilliard problem in mixed form. The space discretization hinges on local reconstruction operators from hybrid polynomial unknowns at elements and faces. The proposed method has several assets: (i) It supports fairly general meshes possibly containing polygonal elements and nonmatching interfaces, (ii) it allows arbitrary approximation orders, (iii) it has a moderate computational cost thanks to the possibility of locally eliminating element-based unknowns by static condensation. We perform a detailed stability and convergence study, proving optimal convergence rates in energy-like norms. Numerical validation is also provided using some of the most common tests in the literature.
△ Less
Submitted 8 April, 2016; v1 submitted 24 September, 2015;
originally announced September 2015.
-
A Hybrid High-Order method for Leray-Lions elliptic equations on general meshes
Authors:
Daniele A. Di Pietro,
Jérôme Droniou
Abstract:
In this work, we develop and analyze a Hybrid High-Order (HHO) method for steady non-linear Leray-Lions problems. The proposed method has several assets, including the support for arbitrary approximation orders and general polytopal meshes. This is achieved by combining two key ingredients devised at the local level: a gradient reconstruction and a high-order stabilization term that generalizes th…
▽ More
In this work, we develop and analyze a Hybrid High-Order (HHO) method for steady non-linear Leray-Lions problems. The proposed method has several assets, including the support for arbitrary approximation orders and general polytopal meshes. This is achieved by combining two key ingredients devised at the local level: a gradient reconstruction and a high-order stabilization term that generalizes the one originally introduced in the linear case. The convergence analysis is carried out using a compactness technique. Extending this technique to HHO methods has prompted us to develop a set of discrete functional analysis tools whose interest goes beyond the specific problem and method addressed in this work: (direct and) reverse Lebesgue and Sobolev embeddings for local polynomial spaces, $L^{p}$-stability and $W^{s,p}$-approximation properties for $L^{2}$-projectors on such spaces, and Sobolev embeddings for hybrid polynomial spaces. Numerical tests are presented to validate the theoretical results for the original method and variants thereof.
△ Less
Submitted 7 December, 2016; v1 submitted 8 August, 2015;
originally announced August 2015.
-
A nonconforming high-order method for the Biot problem on general meshes
Authors:
Daniele Boffi,
Michele Botti,
Daniele A. Di Pietro
Abstract:
In this work, we introduce a novel algorithm for the Biot problem based on a Hybrid High-Order discretization of the mechanics and a Symmetric Weighted Interior Penalty discretization of the flow. The method has several assets, including, in particular, the support of general polyhedral meshes and arbitrary space approximation order. Our analysis delivers stability and error estimates that hold al…
▽ More
In this work, we introduce a novel algorithm for the Biot problem based on a Hybrid High-Order discretization of the mechanics and a Symmetric Weighted Interior Penalty discretization of the flow. The method has several assets, including, in particular, the support of general polyhedral meshes and arbitrary space approximation order. Our analysis delivers stability and error estimates that hold also when the specific storage coefficient vanishes, and shows that the constants have only a mild dependence on the heterogeneity of the permeability coefficient. Numerical tests demonstrating the performance of the method are provided.
△ Less
Submitted 24 February, 2016; v1 submitted 11 June, 2015;
originally announced June 2015.
-
A discontinuous-skeletal method for advection-diffusion-reaction on general meshes
Authors:
Daniele A. Di Pietro,
Jérôme Droniou,
Alexandre Ern
Abstract:
We design and analyze an approximation method for advection-diffusion-reaction equations where the (generalized) degrees of freedom are polynomials of order $k\ge0$ at mesh faces. The method hinges on local discrete reconstruction operators for the diffusive and advective derivatives and a weak enforcement of boundary conditions. Fairly general meshes with polytopal and nonmatching cells are suppo…
▽ More
We design and analyze an approximation method for advection-diffusion-reaction equations where the (generalized) degrees of freedom are polynomials of order $k\ge0$ at mesh faces. The method hinges on local discrete reconstruction operators for the diffusive and advective derivatives and a weak enforcement of boundary conditions. Fairly general meshes with polytopal and nonmatching cells are supported. Arbitrary polynomial orders can be considered, including the case $k=0$ which is closely related to Mimetic Finite Difference/Mixed-Hybrid Finite Volume methods. The error analysis covers the full range of Péclet numbers, including the delicate case of local degeneracy where diffusion vanishes on a strict subset of the domain. Computational costs remain moderate since the use of face unknowns leads to a compact stencil with reduced communications. Numerical results are presented.
△ Less
Submitted 27 May, 2018; v1 submitted 1 November, 2014;
originally announced November 2014.
-
Equilibrated tractions for the Hybrid High-Order method
Authors:
Daniele A. Di Pietro,
Alexandre Ern
Abstract:
We show how to recover equilibrated face tractions for the hybrid high-order method for linear elasticity recently introduced in [D. A. Di Pietro and A. Ern, A hybrid high-order locking-free method for linear elasticity on general meshes, Comput. Meth. Appl. Mech. Engrg., 2015, 283:1-21], and prove that these tractions are optimally convergent.
We show how to recover equilibrated face tractions for the hybrid high-order method for linear elasticity recently introduced in [D. A. Di Pietro and A. Ern, A hybrid high-order locking-free method for linear elasticity on general meshes, Comput. Meth. Appl. Mech. Engrg., 2015, 283:1-21], and prove that these tractions are optimally convergent.
△ Less
Submitted 1 December, 2014; v1 submitted 1 November, 2014;
originally announced November 2014.