113 research outputs found
Vertex-based Compatible Discrete Operator schemes on polyhedral meshes for advection-diffusion equations
International audienceWe devise and analyze vertex-based, Péclet-robust, lowest-order schemes for advection-diffusion equations that support polyhedral meshes. The schemes are formulated using Compatible Discrete Operators (CDO), namely primal and dual discrete differential operators, a discrete contraction operator for advection, and a discrete Hodge operator for diffusion. Moreover, discrete boundary operators are devised to weakly enforce Dirichlet boundary conditions. The analysis sheds new light on the theory of Friedrichs' operators at the purely algebraic level. Moreover, an extension of the stability analysis hinging on inf-sup conditions is presented to incorporate divergence-free velocity fields under some assumptions. Error bounds and convergence rates for smooth solutions are derived, and numerical results are presented on three-dimensional polyhedral meshes
An HMM--ELLAM scheme on generic polygonal meshes for miscible incompressible flows in porous media
We design a numerical approximation of a system of partial differential
equations modelling the miscible displacement of a fluid by another in a porous
medium. The advective part of the system is discretised using a characteristic
method, and the diffusive parts by a finite volume method. The scheme is
applicable on generic (possibly non-conforming) meshes as encountered in
applications. The main features of our work are the reconstruction of a Darcy
velocity, from the discrete pressure fluxes, that enjoys a local consistency
property, an analysis of implementation issues faced when tracking, via the
characteristic method, distorted cells, and a new treatment of cells near the
injection well that accounts better for the conservativity of the injected
fluid
SUPG-stabilized stabilization-free VEM: a numerical investigation
We numerically investigate the possibility of defining stabilization-free
Virtual Element (VEM) discretizations of advection-diffusion problems in the
advection-dominated regime. To this end, we consider a SUPG stabilized
formulation of the scheme. Numerical tests comparing the proposed method with
standard VEM show that the lack of an additional arbitrary stabilization term,
typical of VEM schemes, that adds artificial diffusion to the discrete
solution, allows to better approximate boundary layers, in particular in the
case of a low order scheme.Comment: 15 page
Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains
The numerical approximation of partial differential equations (PDEs) posed on complicated geometries, which include a large number of small geometrical features or microstructures, represents a challenging computational problem. Indeed, the use of standard mesh generators, employing simplices or tensor product elements, for example, naturally leads to very fine finite element meshes, and hence the computational effort required to numerically approximate the underlying PDE problem may be prohibitively expensive. As an alternative approach, in this article we present a review of composite/agglomerated discontinuous Galerkin finite element methods (DGFEMs) which employ general polytopic elements. Here, the elements are typically constructed as the union of standard element shapes; in this way, the minimal dimension of the underlying composite finite element space is independent of the number of geometrical features. In particular, we provide an overview of hp-version inverse estimates and approximation results for general polytopic elements, which are sharp with respect to element facet degeneration. On the basis of these results, a priori error bounds for the hp-DGFEM approximation of both second-order elliptic and first-order hyperbolic PDEs will be derived. Finally, we present numerical experiments which highlight the practical application of DGFEMs on meshes consisting of general polytopic elements
Review of discontinuous Galerkin finite element methods for partial differential equations on complicated domains
The numerical approximation of partial differential equations (PDEs) posed on complicated geometries, which include a large number of small geometrical features or microstructures, represents a challenging computational problem. Indeed, the use of standard mesh generators, employing simplices or tensor product elements, for example, naturally leads to very fine finite element meshes, and hence the computational effort required to numerically approximate the underlying PDE problem may be prohibitively expensive. As an alternative approach, in this article we present a review of composite/agglomerated discontinuous Galerkin finite element methods (DGFEMs) which employ general polytopic elements. Here, the elements are typically constructed as the union of standard element shapes; in this way, the minimal dimension of the underlying composite finite element space is independent of the number of geometrical features. In particular, we provide an overview of hp-version inverse estimates and approximation results for general polytopic elements, which are sharp with respect to element facet degeneration. On the basis of these results, a priori error bounds for the hp-DGFEM approximation of both second-order elliptic and first-order hyperbolic PDEs will be derived. Finally, we present numerical experiments which highlight the practical application of DGFEMs on meshes consisting of general polytopic elements
Arbitrary-Lagrangian-Eulerian discontinuous Galerkin schemes with a posteriori subcell finite volume limiting on moving unstructured meshes
We present a new family of high order accurate fully discrete one-step
Discontinuous Galerkin (DG) finite element schemes on moving unstructured
meshes for the solution of nonlinear hyperbolic PDE in multiple space
dimensions, which may also include parabolic terms in order to model
dissipative transport processes. High order piecewise polynomials are adopted
to represent the discrete solution at each time level and within each spatial
control volume of the computational grid, while high order of accuracy in time
is achieved by the ADER approach. In our algorithm the spatial mesh
configuration can be defined in two different ways: either by an isoparametric
approach that generates curved control volumes, or by a piecewise linear
decomposition of each spatial control volume into simplex sub-elements. Our
numerical method belongs to the category of direct
Arbitrary-Lagrangian-Eulerian (ALE) schemes, where a space-time conservation
formulation of the governing PDE system is considered and which already takes
into account the new grid geometry directly during the computation of the
numerical fluxes. Our new Lagrangian-type DG scheme adopts the novel a
posteriori sub-cell finite volume limiter method, in which the validity of the
candidate solution produced in each cell by an unlimited ADER-DG scheme is
verified against a set of physical and numerical detection criteria. Those
cells which do not satisfy all of the above criteria are flagged as troubled
cells and are recomputed with a second order TVD finite volume scheme. The
numerical convergence rates of the new ALE ADER-DG schemes are studied up to
fourth order in space and time and several test problems are simulated.
Finally, an application inspired by Inertial Confinement Fusion (ICF) type
flows is considered by solving the Euler equations and the PDE of viscous and
resistive magnetohydrodynamics (VRMHD).Comment: 39 pages, 21 figure
The conforming virtual element method for polyharmonic and elastodynamics problems: a review
In this paper, we review recent results on the conforming virtual element
approximation of polyharmonic and elastodynamics problems. The structure and
the content of this review is motivated by three paradigmatic examples of
applications: classical and anisotropic Cahn-Hilliard equation and phase field
models for brittle fracture, that are briefly discussed in the first part of
the paper. We present and discuss the mathematical details of the conforming
virtual element approximation of linear polyharmonic problems, the classical
Cahn-Hilliard equation and linear elastodynamics problems.Comment: 30 pages, 7 figures. arXiv admin note: text overlap with
arXiv:1912.0712
- …