Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA multilayer Saint-Venant system with mass exchanges for Shallow Water flows. Derivation and numerical validation
164
0
0.0
(
0
)
Added by
Jacques Sainte-Marie
Publication date
2009
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The standard multilayer Saint-Venant system consists in introducing fluid layers that are advected by the interfacial velocities. As a consequence there is no mass exchanges between these layers and each layer is described by its height and its average velocity. Here we introduce another multilayer system with mass exchanges between the neighborhing layers where the unknowns are a total height of water and an average velocity per layer. We derive it from Navier-Stokes system with an hydrostatic pressure and prove energy and hyperbolicity properties of the model. We also give a kinetic interpretation leading to effective numerical schemes with positivity and energy properties. Numerical tests show the versatility of the approach and its ability to compute recirculation cases with wind forcing.
rate research
Read More
We pursue here the development of models for complex (viscoelastic) fluids in shallow free-surface gravity flows which was initiated by [Bouchut-Boyaval, M3AS (23) 2013] for 1D (translation invariant) cases. The models we propose are hyperbolic quasilinear systems that generalize Saint-Venant shallow-water equations to incompressible Maxwell fluids. The models are compatible with a formulation of the thermo-dynamics second principle. In comparison with Saint-Venant standard shallow-water model, the momentum balance includes extra-stresses associated with an elastic potential energy in addition to a hydrostatic pressure. The extra-stresses are determined by an additional tensor variable solution to a differential equation with various possible time rates. For the numerical evaluation of solutions to Cauchy problems, we also propose explicit schemes discretizing our generalized Saint-Venant systems with Finite-Volume approximations that are entropy-consistent (under a CFL constraint) in addition to satisfy exact (discrete) mass and momentum conservation laws. In comparison with most standard viscoelastic numerical models, our discrete models can be used for any retardation-time values (i.e. in the vanishing solvent-viscosity limit). We finally illustrate our hyperbolic viscoelastic flow models numerically using computer simulations in benchmark test cases. On extending to Maxwell fluids some free-shear flow testcases that are standard benchmarks for Newtonian fluids, we first show that our (numerical) models reproduce well the viscoelastic physics, phenomenologically at least, with zero retardation-time. Moreover, with a view to quantitative evaluations, numerical results in the lid-driven cavity testcase show that, in fact, our models can be compared with standard viscoelastic flow models in sheared-flow benchmarks on adequately choosing the physical parameters of our models. Analyzing our models asymptotics should therefore shed new light on the famous High-Weissenberg Number Problem (HWNP), which is a limit for all the existing viscoelastic numerical models.
We develop well-balanced central schemes on overlapping cells for the Saint-Venant shallow water system and its variants. The main challenge in deriving the schemes is related to the fact that the Saint-Venant system contains a geometric source term due to nonflat bottom topography and therefore a delicate balance between the flux gradients and source terms has to be preserved. We propose a constant subtraction technique, which helps one to ensure a well-balanced property of the schemes, while maintaining arbitrary high-order of accuracy. Hierarchical reconstruction limiting procedure is applied to eliminate spurious oscillations without using characteristic decomposition. Extensive one- and two-dimensional numerical simulations are conducted to verify the well-balanced property, high-order of accuracy, and non-oscillatory high-resolution for both smooth and nonsmooth solutions.
The purpose of this paper is to derive rigorously the so called viscous shallow water equations given for instance page 958-959 in [A. Oron, S.H. Davis, S.G. Bankoff, Rev. Mod. Phys, 69 (1997), 931?980]. Such a system of equations is similar to compressible Navier-Stokes equations for a barotropic fluid with a non-constant viscosity. To do that, we consider a layer of incompressible and Newtonian fluid which is relatively thin, assuming no surface tension at the free surface. The motion of the fluid is described by 3d Navier-Stokes equations with constant viscosity and free surface. We prove that for a set of suitable initial data (asymptotically close to shallow water initial data), the Cauchy problem for these equations is well-posed, and the solution converges to the solution of viscous shallow water equations. More precisely, we build the solution of the full problem as a perturbation of the strong solution to the viscous shallow water equations. The method of proof is based on a Lagrangian change of variable that fixes the fluid domain and we have to prove the well-posedness in thin domains: we have to pay a special attention to constants in classical Sobolev inequalities and regularity in Stokes problem.
Numerical simulations of flows are required for numerous applications, and are usually carried out using shallow water equations. We describe the FullSWOF software which is based on up-to-date finite volume methods and well-balanced schemes to solve this kind of equations. It consists of a set of open source C++ codes, freely available to the community, easy to use, and open for further development. Several features make FullSWOF particularly suitable for applications in hydrology: small water heights and wet-dry transitions are robustly handled, rainfall and infiltration are incorporated, and data from grid-based digital topographies can be used directly. A detailed mathematical description is given here, and the capabilities of FullSWOF are illustrated based on analytic solutions and datasets of real cases. The codes, available in 1D and
In this paper, we investigate the formation of singularity for general two dimensional and radially symmetric solutions for rotating shallow water system from different aspects. First, the formation of singularity is proved via the study for the associated moments for two dimensional solutions. For the radial symmetric solutions, the formation of singularity is established for the initial data with compact support. Finally, the global existence or formation of singularity for the radial symmetric solutions of the rotating shallow water system are analyzed in detail when the solutions are of the form with separated variables.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
