382 research outputs found
ADER-WENO Finite Volume Schemes with Space-Time Adaptive Mesh Refinement
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
Lagrangian ADER-WENO Finite Volume Schemes on Unstructured Triangular Meshes Based On Genuinely Multidimensional HLL Riemann Solvers
In this paper we use the genuinely multidimensional HLL Riemann solvers
recently developed by Balsara et al. to construct a new class of
computationally efficient high order Lagrangian ADER-WENO one-step ALE finite
volume schemes on unstructured triangular meshes. A nonlinear WENO
reconstruction operator allows the algorithm to achieve high order of accuracy
in space, while high order of accuracy in time is obtained by the use of an
ADER time-stepping technique based on a local space-time Galerkin predictor.
The multidimensional HLL and HLLC Riemann solvers operate at each vertex of the
grid, considering the entire Voronoi neighborhood of each node and allows for
larger time steps than conventional one-dimensional Riemann solvers. The
results produced by the multidimensional Riemann solver are then used twice in
our one-step ALE algorithm: first, as a node solver that assigns a unique
velocity vector to each vertex, in order to preserve the continuity of the
computational mesh; second, as a building block for genuinely multidimensional
numerical flux evaluation that allows the scheme to run with larger time steps
compared to conventional finite volume schemes that use classical
one-dimensional Riemann solvers in normal direction. A rezoning step may be
necessary in order to overcome element overlapping or crossing-over. We apply
the method presented in this article to two systems of hyperbolic conservation
laws, namely the Euler equations of compressible gas dynamics and the equations
of ideal classical magneto-hydrodynamics (MHD). Convergence studies up to
fourth order of accuracy in space and time have been carried out. Several
numerical test problems have been solved to validate the new approach
A multidimensional grid-adaptive relativistic magnetofluid code
A robust second order, shock-capturing numerical scheme for multi-dimensional
special relativistic magnetohydrodynamics on computational domains with
adaptive mesh refinement is presented. The base solver is a total variation
diminishing Lax-Friedrichs scheme in a finite volume setting and is combined
with a diffusive approach for controlling magnetic monopole errors. The
consistency between the primitive and conservative variables is ensured at all
limited reconstructions and the spatial part of the four velocity is used as a
primitive variable. Demonstrative relativistic examples are shown to validate
the implementation. We recover known exact solutions to relativistic MHD
Riemann problems, and simulate the shock-dominated long term evolution of
Lorentz factor 7 vortical flows distorting magnetic island chains.Comment: accepted for publication in Computer Physics Communication
High-Resolution Mathematical and Numerical Analysis of Involution-Constrained PDEs
Partial differential equations constrained by involutions provide the highest fidelity mathematical models for a large number of complex physical systems of fundamental interest in critical scientific and technological disciplines. The applications described by these models include electromagnetics, continuum dynamics of solid media, and general relativity. This workshop brought together pure and applied mathematicians to discuss current research that cuts across these various disciplines’ boundaries. The presented material illuminated fundamental issues as well as evolving theoretical and algorithmic approaches for PDEs with involutions. The scope of the material covered was broad, and the discussions conducted during the workshop were lively and far-reaching
A Second-Order Unsplit Godunov Scheme for Cell-Centered MHD: the CTU-GLM scheme
We assess the validity of a single step Godunov scheme for the solution of
the magneto-hydrodynamics equations in more than one dimension. The scheme is
second-order accurate and the temporal discretization is based on the
dimensionally unsplit Corner Transport Upwind (CTU) method of Colella. The
proposed scheme employs a cell-centered representation of the primary fluid
variables (including magnetic field) and conserves mass, momentum, magnetic
induction and energy. A variant of the scheme, which breaks momentum and energy
conservation, is also considered. Divergence errors are transported out of the
domain and damped using the mixed hyperbolic/parabolic divergence cleaning
technique by Dedner et al. (J. Comput. Phys., 175, 2002). The strength and
accuracy of the scheme are verified by a direct comparison with the eight-wave
formulation (also employing a cell-centered representation) and with the
popular constrained transport method, where magnetic field components retain a
staggered collocation inside the computational cell. Results obtained from two-
and three-dimensional test problems indicate that the newly proposed scheme is
robust, accurate and competitive with recent implementations of the constrained
transport method while being considerably easier to implement in existing hydro
codes.Comment: 31 Pages, 16 Figures Accepted for publication in Journal of
Computational Physic
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
- …