No Arabic abstract
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
Overland flow on agricultural fields may have some undesirable effects such as soil erosion, flood and pollutant transport. To better understand this phenomenon and limit its consequences, we developed a code using state-of-the-art numerical methods: FullSWOF (Full Shallow Water equations for Overland Flow), an object oriented code written in C++. It has been made open-source and can be downloaded from http://www.univ-orleans.fr/mapmo/soft/FullSWOF/. The model is based on the classical system of Shallow Water (SW) (or Saint-Venant system). Numerical difficulties come from the numerous dry/wet transitions and the highly-variable topography encountered inside a field. It includes runon and rainfall inputs, infiltration (modified Green-Ampt equation), friction (Darcy-Weisbach and Manning formulas). First we present the numerical method for the resolution of the Shallow Water equations integrated in FullSWOF_2D (the two-dimension version). This method is based on hydrostatic reconstruction scheme, coupled with a semi-implicit friction term treatment. FullSWOF_2D has been previously validated using analytical solutions from the SWASHES library (Shallow Water Analytic Solutions for Hydraulic and Environmental Studies). Finally, FullSWOF_2D is run on a real topography measured on a runoff plot located in Thies (Senegal). Simulation results are compared with measured data. This experimental benchmark demonstrate the capabilities of FullSWOF to simulate adequately overland flow. FullSWOF could also be used for other environmental issues, such as river floods and dam-breaks.
We propose a variational form of the BDF2 method as an alternative to the commonly used minimizing movement scheme for the time-discrete approximation of gradient flows in abstract metric spaces. Assuming uniform semi-convexity --- but no smoothness --- of the augmented energy functional, we prove well-posedness of the method and convergence of the discrete approximations to a curve of steepest descent. In a smooth Hilbertian setting, classical theory would predict a convergence order of two in time, we prove convergence order of one-half in the general metric setting and under our weak hypotheses. Further, we illustrate these results with numerical experiments for gradient flows on a compact Riemannian manifold, in a Hilbert space, and in the $L^2$-Wasserstein metric.
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.
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.
Many geophysical flow or wave propagation problems can be modeled with two-dimensional depth-averaged equations, of which the shallow water equations are the simplest example. We describe the GeoClaw software that has been designed to solve problems of this nature, consisting of open source Fortran programs together with Python tools for the user interface and flow visualization. This software uses high-resolution shock-capturing finite volume methods on logically rectangular grids, including latitude--longitude grids on the sphere. Dry states are handled automatically to model inundation. The code incorporates adaptive mesh refinement to allow the efficient solution of large-scale geophysical problems. Examples are given illustrating its use for modeling tsunamis, dam break problems, and storm surge. Documentation and download information is available at www.clawpack.org/geoclaw