1,125 research outputs found

    Discrete Lie Advection of Differential Forms

    Get PDF
    In this paper, we present a numerical technique for performing Lie advection of arbitrary differential forms. Leveraging advances in high-resolution finite volume methods for scalar hyperbolic conservation laws, we first discretize the interior product (also called contraction) through integrals over Eulerian approximations of extrusions. This, along with Cartan's homotopy formula and a discrete exterior derivative, can then be used to derive a discrete Lie derivative. The usefulness of this operator is demonstrated through the numerical advection of scalar fields and 1-forms on regular grids.Comment: Accepted version; to be published in J. FoC

    High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes

    Full text link
    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

    Fully computable a posteriori error bounds for hybridizable discontinuous Galerkin finite element approximations

    Get PDF
    We derive a posteriori error estimates for the hybridizable discontinuous Galerkin (HDG) methods, including both the primal and mixed formulations, for the approximation of a linear second-order elliptic problem on conforming simplicial meshes in two and three dimensions. We obtain fully computable, constant free, a posteriori error bounds on the broken energy seminorm and the HDG energy (semi)norm of the error. The estimators are also shown to provide local lower bounds for the HDG energy (semi)norm of the error up to a constant and a higher-order data oscillation term. For the primal HDG methods and mixed HDG methods with an appropriate choice of stabilization parameter, the estimators are also shown to provide a lower bound for the broken energy seminorm of the error up to a constant and a higher-order data oscillation term. Numerical examples are given illustrating the theoretical results

    ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement

    Full text link
    We present the first high order one-step ADER-WENO finite volume scheme with Adaptive Mesh Refinement (AMR) in multiple space dimensions. High order spatial accuracy is obtained through a WENO reconstruction, while a high order one-step time discretization is achieved using a local space-time discontinuous Galerkin predictor method. Due to the one-step nature of the underlying scheme, the resulting algorithm is particularly well suited for an AMR strategy on space-time adaptive meshes, i.e.with time-accurate local time stepping. The AMR property has been implemented 'cell-by-cell', with a standard tree-type algorithm, while the scheme has been parallelized via the Message Passing Interface (MPI) paradigm. The new scheme has been tested over a wide range of examples for nonlinear systems of hyperbolic conservation laws, including the classical Euler equations of compressible gas dynamics and the equations of magnetohydrodynamics (MHD). High order in space and time have been confirmed via a numerical convergence study and a detailed analysis of the computational speed-up with respect to highly refined uniform meshes is also presented. We also show test problems where the presented high order AMR scheme behaves clearly better than traditional second order AMR methods. The proposed scheme that combines for the first time high order ADER methods with space--time adaptive grids in two and three space dimensions is likely to become a useful tool in several fields of computational physics, applied mathematics and mechanics.Comment: With updated bibliography informatio
    • …
    corecore