1,100 research outputs found
The VOLNA code for the numerical modelling of tsunami waves: generation, propagation and inundation
A novel tool for tsunami wave modelling is presented. This tool has the
potential of being used for operational purposes: indeed, the numerical code
\VOLNA is able to handle the complete life-cycle of a tsunami (generation,
propagation and run-up along the coast). The algorithm works on unstructured
triangular meshes and thus can be run in arbitrary complex domains. This paper
contains the detailed description of the finite volume scheme implemented in
the code. The numerical treatment of the wet/dry transition is explained. This
point is crucial for accurate run-up/run-down computations. Most existing
tsunami codes use semi-empirical techniques at this stage, which are not always
sufficient for tsunami hazard mitigation. Indeed the decision to evacuate
inhabitants is based on inundation maps which are produced with this type of
numerical tools. We present several realistic test cases that partially
validate our algorithm. Comparisons with analytical solutions and experimental
data are performed. Finally the main conclusions are outlined and the
perspectives for future research presented.Comment: 47 pages, 27 figures. Other author's papers can be downloaded at
http://www.lama.univ-savoie.fr/~dutykh
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
IllinoisGRMHD: An Open-Source, User-Friendly GRMHD Code for Dynamical Spacetimes
In the extreme violence of merger and mass accretion, compact objects like
black holes and neutron stars are thought to launch some of the most luminous
outbursts of electromagnetic and gravitational wave energy in the Universe.
Modeling these systems realistically is a central problem in theoretical
astrophysics, but has proven extremely challenging, requiring the development
of numerical relativity codes that solve Einstein's equations for the
spacetime, coupled to the equations of general relativistic (ideal)
magnetohydrodynamics (GRMHD) for the magnetized fluids. Over the past decade,
the Illinois Numerical Relativity (ILNR) Group's dynamical spacetime GRMHD code
has proven itself as a robust and reliable tool for theoretical modeling of
such GRMHD phenomena. However, the code was written "by experts and for
experts" of the code, with a steep learning curve that would severely hinder
community adoption if it were open-sourced. Here we present IllinoisGRMHD,
which is an open-source, highly-extensible rewrite of the original
closed-source GRMHD code of the ILNR Group. Reducing the learning curve was the
primary focus of this rewrite, with the goal of facilitating community
involvement in the code's use and development, as well as the minimization of
human effort in generating new science. IllinoisGRMHD also saves computer time,
generating roundoff-precision identical output to the original code on
adaptive-mesh grids, but nearly twice as fast at scales of hundreds to
thousands of cores.Comment: 37 pages, 6 figures, single column. Matches published versio
Athena: A New Code for Astrophysical MHD
A new code for astrophysical magnetohydrodynamics (MHD) is described. The
code has been designed to be easily extensible for use with static and adaptive
mesh refinement. It combines higher-order Godunov methods with the constrained
transport (CT) technique to enforce the divergence-free constraint on the
magnetic field. Discretization is based on cell-centered volume-averages for
mass, momentum, and energy, and face-centered area-averages for the magnetic
field. Novel features of the algorithm include (1) a consistent framework for
computing the time- and edge-averaged electric fields used by CT to evolve the
magnetic field from the time- and area-averaged Godunov fluxes, (2) the
extension to MHD of spatial reconstruction schemes that involve a
dimensionally-split time advance, and (3) the extension to MHD of two different
dimensionally-unsplit integration methods. Implementation of the algorithm in
both C and Fortran95 is detailed, including strategies for parallelization
using domain decomposition. Results from a test suite which includes problems
in one-, two-, and three-dimensions for both hydrodynamics and MHD are given,
not only to demonstrate the fidelity of the algorithms, but also to enable
comparisons to other methods. The source code is freely available for download
on the web.Comment: 61 pages, 36 figures. accepted by ApJ
3D free surface flow simulations based on the integral form of the equations of motion
This work deals with a novel three-dimensional finite-volume non-hydrostatic shock-capturing model for the simulation of wave transformation processes and wave-structure interaction. The model is based on an integral formulation of the Navier-Stokes equations solved on a coordinate system in which the vertical coordinate is varying in time. A finite-volume shock-capturing numerical technique based on high order WENO reconstructions is adopted in order to discretize the fluid motion equations
ECHO: an Eulerian Conservative High Order scheme for general relativistic magnetohydrodynamics and magnetodynamics
We present a new numerical code, ECHO, based on an Eulerian Conservative High
Order scheme for time dependent three-dimensional general relativistic
magnetohydrodynamics (GRMHD) and magnetodynamics (GRMD). ECHO is aimed at
providing a shock-capturing conservative method able to work at an arbitrary
level of formal accuracy (for smooth flows), where the other existing GRMHD and
GRMD schemes yield an overall second order at most. Moreover, our goal is to
present a general framework, based on the 3+1 Eulerian formalism, allowing for
different sets of equations, different algorithms, and working in a generic
space-time metric, so that ECHO may be easily coupled to any solver for
Einstein's equations. Various high order reconstruction methods are implemented
and a two-wave approximate Riemann solver is used. The induction equation is
treated by adopting the Upwind Constrained Transport (UCT) procedures,
appropriate to preserve the divergence-free condition of the magnetic field in
shock-capturing methods. The limiting case of magnetodynamics (also known as
force-free degenerate electrodynamics) is implemented by simply replacing the
fluid velocity with the electromagnetic drift velocity and by neglecting the
matter contribution to the stress tensor. ECHO is particularly accurate,
efficient, versatile, and robust. It has been tested against several
astrophysical applications, including a novel test on the propagation of large
amplitude circularly polarized Alfven waves. In particular, we show that
reconstruction based on a Monotonicity Preserving filter applied to a fixed
5-point stencil gives highly accurate results for smooth solutions, both in
flat and curved metric (up to the nominal fifth order), while at the same time
providing sharp profiles in tests involving discontinuities.Comment: 20 pages, revised version submitted to A&
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
- …