952 research outputs found
FullSWOF_Paral: Comparison of two parallelization strategies (MPI and SKELGIS) on a software designed for hydrology applications
In this paper, we perform a comparison of two approaches for the
parallelization of an existing, free software, FullSWOF 2D (http://www.
univ-orleans.fr/mapmo/soft/FullSWOF/ that solves shallow water equations for
applications in hydrology) based on a domain decomposition strategy. The first
approach is based on the classical MPI library while the second approach uses
Parallel Algorithmic Skeletons and more precisely a library named SkelGIS
(Skeletons for Geographical Information Systems). The first results presented
in this article show that the two approaches are similar in terms of
performance and scalability. The two implementation strategies are however very
different and we discuss the advantages of each one.Comment: 27 page
A limitation of the hydrostatic reconstruction technique for Shallow Water equations
Because of their capability to preserve steady-states, well-balanced schemes
for Shallow Water equations are becoming popular. Among them, the hydrostatic
reconstruction proposed in Audusse et al. (2004), coupled with a positive
numerical flux, allows to verify important mathematical and physical properties
like the positivity of the water height and, thus, to avoid unstabilities when
dealing with dry zones. In this note, we prove that this method exhibits an
abnormal behavior for some combinations of slope, mesh size and water height.Comment: 7 page
A "well-balanced" finite volume scheme for blood flow simulation
We are interested in simulating blood flow in arteries with a one dimensional
model. Thanks to recent developments in the analysis of hyperbolic system of
conservation laws (in the Saint-Venant/ shallow water equations context) we
will perform a simple finite volume scheme. We focus on conservation properties
of this scheme which were not previously considered. To emphasize the necessity
of this scheme, we present how a too simple numerical scheme may induce
spurious flows when the basic static shape of the radius changes. On contrary,
the proposed scheme is "well-balanced": it preserves equilibria of Q = 0. Then
examples of analytical or linearized solutions with and without viscous damping
are presented to validate the calculations. The influence of abrupt change of
basic radius is emphasized in the case of an aneurism.Comment: 36 page
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
A large time-step and well-balanced Lagrange-Projection type scheme for the shallow-water equations
This work focuses on the numerical approximation of the Shallow Water
Equations (SWE) using a Lagrange-Projection type approach. We propose to extend
to this context recent implicit-explicit schemes developed in the framework of
compressibleflows, with or without stiff source terms. These methods enable the
use of time steps that are no longer constrained by the sound velocity thanks
to an implicit treatment of the acoustic waves, and maintain accuracy in the
subsonic regime thanks to an explicit treatment of the material waves. In the
present setting, a particular attention will be also given to the
discretization of the non-conservative terms in SWE and more specifically to
the well-known well-balanced property. We prove that the proposed numerical
strategy enjoys important non linear stability properties and we illustrate its
behaviour past several relevant test cases
- …