6,449 research outputs found
Parallelization of a relaxation scheme modelling the bedload transport of sediments in shallow water flow
In this work we are interested in numerical simulations for bedload erosion
processes. We present a relaxation solver that we apply to moving dunes test
cases in one and two dimensions. In particular we retrieve the so-called
anti-dune process that is well described in the experiments. In order to be
able to run 2D test cases with reasonable CPU time, we also describe and apply
a parallelization procedure by using domain decomposition based on the
classical MPI library.Comment: 19 page
An unstructured finite-volume method for coupled models of suspended sediment and bed load transport in shallow-water flows
The aim of this work is to develop a well-balanced finite-volume method for the accurate numerical solution of the equations governing suspended sediment and bed load transport in two-dimensional shallow-water flows. The modelling system consists of three coupled model components: (i) the shallow-water equations for the hydrodynamical model; (ii) a transport equation for the dispersion of suspended sediments; and (iii) an Exner equation for the morphodynamics. These coupled models form a hyperbolic system of conservation laws with source terms. The proposed finite-volume method consists of a predictor stage for the discretization of gradient terms and a corrector stage for the treatment of source terms. The gradient fluxes are discretized using a modified Roe's scheme using the sign of the Jacobian matrix in the coupled system. A well-balanced discretization is used for the treatment of source terms. In this paper, we also employ an adaptive procedure in the finite-volume method by monitoring the concentration of suspended sediments in the computational domain during its transport process. The method uses unstructured meshes and incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep sediment concentrations and bed load gradients that may form in the approximate solutions. Details are given on the implementation of the method, and numerical results are presented for two idealized test cases, which demonstrate the accuracy and robustness of the method and its applicability in predicting dam-break flows over erodible sediment beds. The method is also applied to a sediment transport problem in the Nador lagoon
On the influence of the thickness of the sediment moving layer in the definition of the bedload transport formula in Exner systems
In this paper we study Exner system and introduce a modified general definition for bedload transport flux. The new formulation has the advantage of taking into account the thickness of the sediment layer which avoids mass conservation problems in certain situations. Moreover, it reduces to a classical solid transport discharge formula in the case of quasi-uniform regime. We also present several numerical tests where we compare the proposed sediment transport formula with the classical formulation and we show the behavior of the new model in different configurations
Resonance and morphological stability of tidal basins
The paper describes the concept of a network model for the morphological behaviour of a near-resonant multiple-inlet tidal basin, as part of a model system which includes the barrier island coasts and the outer deltas. It addresses the question whether a small interference somewhere in such a basin can have major effects on sediment transport and morphology elsewhere in the system.\ud
\ud
In order to explain the basic ideas of the model, only the main tidal constituent (M2) and the associated topography-induced residual current are considered, not the overtides. Furthermore, the model concerns only non-cohesive sediment (sand).\ud
\ud
In spite of these simplifications, the model concept is shown to be effective, in that it results in a morphological evolution equation for each branch of the network and a picture of the influence of each branch on the resonance-sensitivity of the system as a whole
Derivation of a multilayer approach to model suspended sediment transport: application to hyperpycnal and hypopycnal plumes
We propose a multi-layer approach to simulate hyperpycnal and hypopycnal
plumes in flows with free surface. The model allows to compute the vertical
profile of the horizontal and the vertical components of the velocity of the
fluid flow. The model can describe as well the vertical profile of the sediment
concentration and the velocity components of each one of the sediment species
that form the turbidity current. To do so, it takes into account the settling
velocity of the particles and their interaction with the fluid. This allows to
better describe the phenomena than a single layer approach. It is in better
agreement with the physics of the problem and gives promising results. The
numerical simulation is carried out by rewriting the multi-layer approach in a
compact formulation, which corresponds to a system with non-conservative
products, and using path-conservative numerical scheme. Numerical results are
presented in order to show the potential of the model
Modeling bed erosion in free surface flows by the particle finite element method
We present a general formulation for modeling bed erosion in free surface flows using the particle finite element method (PFEM). The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes (particles) in domains containing fluid and solid subdomains. Nodes are viewed as material points (called particles) which can freely move and even separate from the fluid and solid subdomains representing, for instance, the effect of water drops or soil/rock particles. A mesh connects the nodes defining the discretized domain in the fluid and solid regions where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the nonlinear transient coupled fluid-structure problem is described. The erosion mechanism is modeled by releasing the material adjacent to the bed surface according to the frictional work generated by the fluid shear stresses. The released bed material is subsequently transported by the fluid flow. Examples of application of the PFEM to solve a number of bed erosion problems involving large motions of the free surface and splashing of waves are presented
Modeling bed erosion in free surface flows by the particle finite element method
We present a general formulation for modeling bed erosion in free surface flows using the particle finite element method (PFEM). The key feature of the PFEM is the use of an updated Lagrangian description to model the motion of nodes (particles) in domains containing fluid and solid subdomains. Nodes are viewed as material points (called particles) which can freely move and even separate from the fluid and solid subdomains representing, for instance, the effect of water drops or soil/rock particles. A mesh connects the nodes defining the discretized domain in the fluid and solid regions where the governing equations, expressed in an integral form, are solved as in the standard FEM. The necessary stabilization for dealing with the incompressibility of the fluid is introduced via the finite calculus (FIC) method. An incremental iterative scheme for the solution of the nonlinear transient coupled fluid-structure problem is described. The erosion mechanism is modeled by releasing the material adjacent to the bed surface according to the frictional work generated by the fluid shear stresses. The released bed material is subsequently transported by the fluid flow. Examples of application of the PFEM to solve a number of bed erosion problems involving large motions of the free surface and splashing of waves are presented
Shallow Water Moment models for bedload transport problems
In this work a simple but accurate shallow model for bedload sediment
transport is proposed. The model is based on applying the moment approach to
the Shallow Water Exner model, making it possible to recover the vertical
structure of the flow. This approach allows us to obtain a better approximation
of the fluid velocity close to the bottom, which is the relevant velocity for
the sediment transport. A general Shallow Water Exner moment model allowing for
polynomial velocity profiles of arbitrary order is obtained. A regularization
ensures hyperbolicity and easy computation of the eigenvalues. The system is
solved by means of an adapted IFCP scheme proposed here. The improvement of
this IFCP type scheme is based on the approximation of the eigenvalue
associated to the sediment transport. Numerical tests are presented which deal
with large and short time scales. The proposed model allows to obtain the
vertical structure of the fluid, which results in a better description on the
bedload transport of the sediment layer
Modeling Shallow Water Flows on General Terrains
A formulation of the shallow water equations adapted to general complex
terrains is proposed. Its derivation starts from the observation that the
typical approach of depth integrating the Navier-Stokes equations along the
direction of gravity forces is not exact in the general case of a tilted curved
bottom. We claim that an integration path that better adapts to the shallow
water hypotheses follows the "cross-flow" surface, i.e., a surface that is
normal to the velocity field at any point of the domain. Because of the
implicitness of this definition, we approximate this "cross-flow" path by
performing depth integration along a local direction normal to the bottom
surface, and propose a rigorous derivation of this approximation and its
numerical solution as an essential step for the future development of the full
"cross-flow" integration procedure. We start by defining a local coordinate
system, anchored on the bottom surface to derive a covariant form of the
Navier-Stokes equations. Depth integration along the local normals yields a
covariant version of the shallow water equations, which is characterized by
flux functions and source terms that vary in space because of the surface
metric coefficients and related derivatives. The proposed model is discretized
with a first order FORCE-type Godunov Finite Volume scheme that allows
implementation of spatially variable fluxes. We investigate the validity of our
SW model and the effects of the bottom geometry by means of three synthetic
test cases that exhibit non negligible slopes and surface curvatures. The
results show the importance of taking into consideration bottom geometry even
for relatively mild and slowly varying curvatures
- …