85,866 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
High order direct Arbitrary-Lagrangian-Eulerian schemes on moving Voronoi meshes with topology changes
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
2D granular flows with the rheology and side walls friction: a well balanced multilayer discretization
We present here numerical modelling of granular flows with the
rheology in confined channels. The contribution is twofold: (i) a model to
approximate the Navier-Stokes equations with the rheology through an
asymptotic analysis. Under the hypothesis of a one-dimensional flow, this model
takes into account side walls friction; (ii) a multilayer discretization
following Fern\'andez-Nieto et al. (J. Fluid Mech., vol. 798, 2016, pp.
643-681). In this new numerical scheme, we propose an appropriate treatment of
the rheological terms through a hydrostatic reconstruction which allows this
scheme to be well-balanced and therefore to deal with dry areas. Based on
academic tests, we first evaluate the influence of the width of the channel on
the normal profiles of the downslope velocity thanks to the multilayer approach
that is intrinsically able to describe changes from Bagnold to S-shaped (and
vice versa) velocity profiles. We also check the well balance property of the
proposed numerical scheme. We show that approximating side walls friction using
single-layer models may lead to strong errors. Secondly, we compare the
numerical results with experimental data on granular collapses. We show that
the proposed scheme allows us to qualitatively reproduce the deposit in the
case of a rigid bed (i. e. dry area) and that the error made by replacing the
dry area by a small layer of material may be large if this layer is not thin
enough. The proposed model is also able to reproduce the time evolution of the
free surface and of the flow/no-flow interface. In addition, it reproduces the
effect of erosion for granular flows over initially static material lying on
the bed. This is possible when using a variable friction coefficient
but not with a constant friction coefficient
- …