6,449 research outputs found

    Parallelization of a relaxation scheme modelling the bedload transport of sediments in shallow water flow

    Get PDF
    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

    Get PDF
    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

    Get PDF
    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

    Get PDF
    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

    Full text link
    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

    Get PDF
    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

    Get PDF
    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

    Get PDF
    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

    Full text link
    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
    • …
    corecore