40 research outputs found
Hybridisable Discontinuous Galerkin Formulation of Compressible Flows
This work presents a review of high-order hybridisable discontinuous Galerkin
(HDG) methods in the context of compressible flows. Moreover, an original
unified framework for the derivation of Riemann solvers in hybridised
formulations is proposed. This framework includes, for the first time in an HDG
context, the HLL and HLLEM Riemann solvers as well as the traditional
Lax-Friedrichs and Roe solvers. HLL-type Riemann solvers demonstrate their
superiority with respect to Roe in supersonic cases due to their positivity
preserving properties. In addition, HLLEM specifically outstands in the
approximation of boundary layers because of its shear preservation, which
confers it an increased accuracy with respect to HLL and Lax-Friedrichs. A
comprehensive set of relevant numerical benchmarks of viscous and inviscid
compressible flows is presented. The test cases are used to evaluate the
competitiveness of the resulting high-order HDG scheme with the aforementioned
Riemann solvers and equipped with a shock treatment technique based on
artificial viscosity.Comment: 60 pages, 31 figures. arXiv admin note: substantial text overlap with
arXiv:1912.0004
High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
We present a new family of very high order accurate direct
Arbitrary-Lagrangian-Eulerian (ALE) Finite Volume (FV) and Discontinuous
Galerkin (DG) schemes for the solution of nonlinear hyperbolic PDE systems on
moving 2D Voronoi meshes that are regenerated at each time step and which
explicitly allow topology changes in time.
The Voronoi tessellations are obtained from a set of generator points that
move with the local fluid velocity. We employ an AREPO-type approach, which
rapidly rebuilds a new high quality mesh rearranging the element shapes and
neighbors in order to guarantee a robust mesh evolution even for vortex flows
and very long simulation times. The old and new Voronoi elements associated to
the same generator are connected to construct closed space--time control
volumes, whose bottom and top faces may be polygons with a different number of
sides. We also incorporate degenerate space--time sliver elements, needed to
fill the space--time holes that arise because of topology changes. The final
ALE FV-DG scheme is obtained by a redesign of the fully discrete direct ALE
schemes of Boscheri and Dumbser, extended here to moving Voronoi meshes and
space--time sliver elements. Our new numerical scheme is based on the
integration over arbitrary shaped closed space--time control volumes combined
with a fully-discrete space--time conservation formulation of the governing PDE
system. In this way the discrete solution is conservative and satisfies the GCL
by construction.
Numerical convergence studies as well as a large set of benchmarks for
hydrodynamics and magnetohydrodynamics (MHD) demonstrate the accuracy and
robustness of the proposed method. Our numerical results clearly show that the
new combination of very high order schemes with regenerated meshes with
topology changes lead to substantial improvements compared to direct ALE
methods on conforming meshes
A staggered semi-implicit hybrid finite volume / finite element scheme for the shallow water equations at all Froude numbers
We present a novel staggered semi-implicit hybrid FV/FE method for the
numerical solution of the shallow water equations at all Froude numbers on
unstructured meshes. A semi-discretization in time of the conservative
Saint-Venant equations with bottom friction terms leads to its decomposition
into a first order hyperbolic subsystem containing the nonlinear convective
term and a second order wave equation for the pressure. For the spatial
discretization of the free surface elevation an unstructured mesh of triangular
simplex elements is considered, whereas a dual grid of the edge-type is
employed for the computation of the depth-averaged momentum vector. The first
stage of the proposed algorithm consists in the solution of the nonlinear
convective subsystem using an explicit Godunov-type FV method on the staggered
grid. Next, a classical continuous FE scheme provides the free surface
elevation at the vertex of the primal mesh. The semi-implicit strategy followed
circumvents the contribution of the surface wave celerity to the CFL-type time
step restriction making the proposed algorithm well-suited for low Froude
number flows. The conservative formulation of the governing equations also
allows the discretization of high Froude number flows with shock waves. As
such, the new hybrid FV/FE scheme is able to deal simultaneously with both,
subcritical as well as supercritical flows. Besides, the algorithm is well
balanced by construction. The accuracy of the overall methodology is studied
numerically and the C-property is proven theoretically and validated via
numerical experiments. The solution of several Riemann problems attests the
robustness of the new method to deal also with flows containing bores and
discontinuities. Finally, a 3D dam break problem over a dry bottom is studied
and our numerical results are successfully compared with numerical reference
solutions and experimental data