5,767 research outputs found
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
Entropy Stable Finite Volume Approximations for Ideal Magnetohydrodynamics
This article serves as a summary outlining the mathematical entropy analysis
of the ideal magnetohydrodynamic (MHD) equations. We select the ideal MHD
equations as they are particularly useful for mathematically modeling a wide
variety of magnetized fluids. In order to be self-contained we first motivate
the physical properties of a magnetic fluid and how it should behave under the
laws of thermodynamics. Next, we introduce a mathematical model built from
hyperbolic partial differential equations (PDEs) that translate physical laws
into mathematical equations. After an overview of the continuous analysis, we
thoroughly describe the derivation of a numerical approximation of the ideal
MHD system that remains consistent to the continuous thermodynamic principles.
The derivation of the method and the theorems contained within serve as the
bulk of the review article. We demonstrate that the derived numerical
approximation retains the correct entropic properties of the continuous model
and show its applicability to a variety of standard numerical test cases for
MHD schemes. We close with our conclusions and a brief discussion on future
work in the area of entropy consistent numerical methods and the modeling of
plasmas
Entropy stable DGSEM for nonlinear hyperbolic systems in nonconservative form with application to two-phase flows
In this work, we consider the discretization of nonlinear hyperbolic systems
in nonconservative form with the high-order discontinuous Galerkin spectral
element method (DGSEM) based on collocation of quadrature and interpolation
points (Kopriva and Gassner, J. Sci. Comput., 44 (2010), pp.136--155; Carpenter
et al., SIAM J. Sci. Comput., 36 (2014), pp.~B835-B867). We present a general
framework for the design of such schemes that satisfy a semi-discrete entropy
inequality for a given convex entropy function at any approximation order. The
framework is closely related to the one introduced for conservation laws by
Chen and Shu (J. Comput. Phys., 345 (2017), pp.~427--461) and relies on the
modification of the integral over discretization elements where we replace the
physical fluxes by entropy conservative numerical fluxes from Castro et al.
(SIAM J. Numer. Anal., 51 (2013), pp.~1371--1391), while entropy stable
numerical fluxes are used at element interfaces. Time discretization is
performed with strong-stability preserving Runge-Kutta schemes. We use this
framework for the discretization of two systems in one space-dimension: a
system with a nonconservative product associated to a
linearly-degenerate field for which the DGSEM fails to capture the physically
relevant solution, and the isentropic Baer-Nunziato model. For the latter, we
derive conditions on the numerical parameters of the discrete scheme to further
keep positivity of the partial densities and a maximum principle on the void
fractions. Numerical experiments support the conclusions of the present
analysis and highlight stability and robustness of the present schemes
An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry
We design an arbitrary high-order accurate nodal discontinuous Galerkin
spectral element approximation for the nonlinear two dimensional shallow water
equations with non-constant, possibly discontinuous, bathymetry on
unstructured, possibly curved, quadrilateral meshes. The scheme is derived from
an equivalent flux differencing formulation of the split form of the equations.
We prove that this discretisation exactly preserves the local mass and
momentum. Furthermore, combined with a special numerical interface flux
function, the method exactly preserves the mathematical entropy, which is the
total energy for the shallow water equations. By adding a specific form of
interface dissipation to the baseline entropy conserving scheme we create a
provably entropy stable scheme. That is, the numerical scheme discretely
satisfies the second law of thermodynamics. Finally, with a particular
discretisation of the bathymetry source term we prove that the numerical
approximation is well-balanced. We provide numerical examples that verify the
theoretical findings and furthermore provide an application of the scheme for a
partial break of a curved dam test problem
- …