12,939 research outputs found
Discontinuous Galerkin approximations in computational mechanics: hybridization, exact geometry and degree adaptivity
Discontinuous Galerkin (DG) discretizations with exact representation of the geometry and local polynomial degree adaptivity are revisited. Hybridization techniques are employed to reduce the computational cost of DG approximations and devise the hybridizable discontinuous Galerkin (HDG) method. Exact geometry described by non-uniform rational B-splines (NURBS) is integrated into HDG using the framework of the NURBS-enhanced finite element method (NEFEM). Moreover, optimal convergence and superconvergence properties of HDG-Voigt formulation in presence of symmetric second-order tensors are exploited to construct inexpensive error indicators and drive degree adaptive procedures. Applications involving the numerical simulation of problems in electrostatics, linear elasticity and incompressible viscous flows are presented. Moreover, this is done for both high-order HDG approximations and the lowest-order framework of face-centered finite volumes (FCFV).Peer ReviewedPostprint (author's final draft
On the stability of projection methods for the incompressible Navier-Stokes equations based on high-order discontinuous Galerkin discretizations
The present paper deals with the numerical solution of the incompressible
Navier-Stokes equations using high-order discontinuous Galerkin (DG) methods
for discretization in space. For DG methods applied to the dual splitting
projection method, instabilities have recently been reported that occur for
coarse spatial resolutions and small time step sizes. By means of numerical
investigation we give evidence that these instabilities are related to the
discontinuous Galerkin formulation of the velocity divergence term and the
pressure gradient term that couple velocity and pressure. Integration by parts
of these terms with a suitable definition of boundary conditions is required in
order to obtain a stable and robust method. Since the intermediate velocity
field does not fulfill the boundary conditions prescribed for the velocity, a
consistent boundary condition is derived from the convective step of the dual
splitting scheme to ensure high-order accuracy with respect to the temporal
discretization. This new formulation is stable in the limit of small time steps
for both equal-order and mixed-order polynomial approximations. Although the
dual splitting scheme itself includes inf-sup stabilizing contributions, we
demonstrate that spurious pressure oscillations appear for equal-order
polynomials and small time steps highlighting the necessity to consider inf-sup
stability explicitly.Comment: 31 page
Stabilised hybrid discontinuous Galerkin methods for the Stokes problem with non-standard boundary conditions
In several studies it has been observed that, when using stabilised
elements for both velocity and pressure,
the error for the pressure is smaller, or even of a higher order in some cases,
than the one obtained when using inf-sup stable
(although no formal proof of either
of these facts has been given). This increase in polynomial order requires the
introduction of stabilising terms, since the finite element pairs used do not
stability the inf-sup condition. With this motivation, we apply the
stabilisation approach to the hybrid discontinuous Galerkin discretisation for
the Stokes problem with non-standard boundary conditions
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
An optimal order interior penalty discontinuous Galerkin discretization of the compressible Navier-Stokes equations
In this article we propose a new symmetric version of the
interior penalty discontinuous Galerkin finite element
method for the numerical approximation
of the compressible Navier-Stokes equations. Here, particular
emphasis is devoted to the construction of an optimal numerical
method for the evaluation of certain target functionals of practical
interest, such as the lift and drag coefficients of a body immersed in a
viscous fluid. With this in mind, the key ingredients in the
construction of the method include: (i) An adjoint consistent imposition of
the boundary conditions; (ii) An adjoint consistent reformulation of the
underlying target functional of practical interest; (iii) Design of appropriate
interior--penalty stabilization terms. Numerical experiments presented
within this article clearly indicate the optimality of the proposed
method when the error is measured in terms of both the L2-norm, as well as
for certain target functionals. Computational comparisons with
other discontinuous Galerkin schemes proposed in the literature,
including the second scheme
of Bassi and Rebay, the standard SIPG method
outlined in [Hartmann,Houston-2006], and an NIPG variant of the new scheme
will be undertaken
An advection-robust Hybrid High-Order method for the Oseen problem
In this work, we study advection-robust Hybrid High-Order discretizations of
the Oseen equations. For a given integer , the discrete velocity
unknowns are vector-valued polynomials of total degree on mesh elements
and faces, while the pressure unknowns are discontinuous polynomials of total
degree on the mesh. From the discrete unknowns, three relevant
quantities are reconstructed inside each element: a velocity of total degree
, 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 of diameter contributes to the discretization error with
an -term in the diffusion-dominated regime, an
-term in the advection-dominated regime, and
scales with intermediate powers of in between. Numerical results complete
the exposition
- …