Normalized to: Fambri, F.
[1]
oai:arXiv.org:1801.02839 [pdf] - 1689150
ADER discontinuous Galerkin schemes for general-relativistic ideal
magnetohydrodynamics
Submitted: 2018-01-09, last modified: 2018-05-25
We present a new class of high-order accurate numerical algorithms for
solving the equations of general-relativistic ideal magnetohydrodynamics in
curved spacetimes. In this paper we assume the background spacetime to be given
and static, i.e., we make use of the Cowling approximation. The governing
partial differential equations are solved via a new family of fully-discrete
and arbitrary high-order accurate path-conservative discontinuous Galerkin (DG)
finite-element methods combined with adaptive mesh refinement and time accurate
local timestepping. In order to deal with shock waves and other
discontinuities, the highorder DG schemes are supplemented with a novel
a-posteriori subcell finite-volume limiter, which makes the new algorithms as
robust as classical second-order total-variation diminishing finite-volume
methods at shocks and discontinuities, but also as accurate as unlimited
high-order DG schemes in smooth regions of the flow. We show the advantages of
this new approach by means of various classical two- and three-dimensional
benchmark problems on fixed spacetimes. Finally, we present a performance and
accuracy comparisons between Runge-Kutta DG schemes and ADER high-order
finite-volume schemes, showing the higher efficiency of DG schemes.
[2]
oai:arXiv.org:1612.01410 [pdf] - 1580805
Space-time adaptive ADER-DG schemes for dissipative flows: compressible
Navier-Stokes and resistive MHD equations
Submitted: 2016-12-05
This paper presents an arbitrary h.o. accurate ADER DG method on space-time
adaptive meshes (AMR) for the solution of two important families of non-linear
time dependent PDE for compr. dissipative flows: the compr. Navier-Stokes
equations and the equations of visc. and res. MHD in 2 and 3 space-dimensions.
The work continues a recent series of papers concerning the development and
application of a proper a posteriori subcell FV limiting procedure suitable for
DG methods. It is a well known fact that a major weakness of h.o. DG methods
lies in the difficulty of limiting discontinuous solutions, which generate
spurious oscillations, namely the so-called 'Gibbs phenomenon'. In the present
work the main benefits of the MOOD paradigm, i.e. the computational robustness
even in the presence of strong shocks, are preserved and the numerical
diffusion is considerably reduced also for the limited cells by resorting to a
proper sub-grid. An important feature of our new scheme is its ability to cure
even floating point errors that may occur during a simulation, for example when
taking real roots of negative numbers or after divisions by zero. We apply the
whole approach for the first time to the equations of compr. gas dynamics and
MHD in the presence of viscosity, thermal conductivity and magnetic
resistivity, therefore extending our family of adaptive ADER-DG schemes to
cases for which the numerical fluxes also depend on the gradient of the state
vector. The distinguished high-resolution properties of the presented numerical
scheme stands out against a wide number of non-trivial test cases both for the
compr. Navier-Stokes and the viscous and resistive MHD equations. The present
results show clearly that the shock-capturing capability of the news schemes
are significantly enhanced within a cell-by-cell Adaptive Mesh Refinement
implementation together with time accurate local time stepping (LTS).
[3]
oai:arXiv.org:1412.0081 [pdf] - 1270604
Space-time adaptive ADER discontinuous Galerkin finite element schemes
with a posteriori sub-cell finite volume limiting
Submitted: 2014-11-29, last modified: 2015-09-01
In this paper we present a novel arbitrary high order accurate discontinuous
Galerkin (DG) finite element method on space-time adaptive Cartesian meshes
(AMR) for hyperbolic conservation laws in multiple space dimensions, using a
high order \aposteriori sub-cell ADER-WENO finite volume \emph{limiter}.
Notoriously, the original DG method produces strong oscillations in the
presence of discontinuous solutions and several types of limiters have been
introduced over the years to cope with this problem. Following the innovative
idea recently proposed in \cite{Dumbser2014}, the discrete solution within the
troubled cells is \textit{recomputed} by scattering the DG polynomial at the
previous time step onto a suitable number of sub-cells along each direction.
Relying on the robustness of classical finite volume WENO schemes, the sub-cell
averages are recomputed and then gathered back into the DG polynomials over the
main grid. In this paper this approach is implemented for the first time within
a space-time adaptive AMR framework in two and three space dimensions, after
assuring the proper averaging and projection between sub-cells that belong to
different levels of refinement. The combination of the sub-cell resolution with
the advantages of AMR allows for an unprecedented ability in resolving even the
finest details in the dynamics of the fluid. The spectacular resolution
properties of the new scheme have been shown through a wide number of test
cases performed in two and in three space dimensions, both for the Euler
equations of compressible gas dynamics and for the magnetohydrodynamics (MHD)
equations.
[4]
oai:arXiv.org:1504.07458 [pdf] - 1260782
Solving the relativistic magnetohydrodynamics equations with ADER
discontinuous Galerkin methods, a posteriori subcell limiting and adaptive
mesh refinement
Submitted: 2015-04-28, last modified: 2015-08-10
We present a new numerical tool for solving the special relativistic ideal
MHD equations that is based on the combination of the following three key
features: (i) a one-step ADER discontinuous Galerkin (DG) scheme that allows
for an arbitrary order of accuracy in both space and time, (ii) an a posteriori
subcell finite volume limiter that is activated to avoid spurious oscillations
at discontinuities without destroying the natural subcell resolution
capabilities of the DG finite element framework and finally (iii) a space-time
adaptive mesh refinement (AMR) framework with time-accurate local
time-stepping. The divergence-free character of the magnetic field is instead
taken into account through the so-called "divergence-cleaning" approach. The
convergence of the new scheme is verified up to 5th order in space and time and
the results for a set of significant numerical tests including shock tube
problems, the RMHD rotor and blast wave problems, as well as the Orszag-Tang
vortex system are shown. We also consider a simple case of the relativistic
Kelvin-Helmholtz instability with a magnetic field, emphasizing the potential
of the new method for studying turbulent RMHD flows. We discuss the advantages
of our new approach when the equations of relativistic MHD need to be solved
with high accuracy within various astrophysical systems.