Search

[Departement_IRSTEA]Eaux [TR1_IRSTEA]RIVAGE [ADD1_IRSTEA]Hydrosystèmes et risques naturels; International audience; The ability of viscoplastic fluids to self-select their flow geometry through the formation of unyielded dead zones has important consequences for flows in wavy channels, flows over obstacles, etc. Yet, the mechanisms controlling the formation and dimensions of the dead zones remain poorly understood. We present a detailed cross-comparison of experimental and numerical results concerning channel flows of a viscoplastic fluid over a rectangular cavity filled by the same material. In all the configurations studied, which correspond to moderate values of the Bingham number, a continuous dead zone forms inside the cavity. Both numerical and experimental data reveal that, unlike at high Bingham numbers, the shear-rate profiles above the dead zone display an asymmetric shape. Accordingly, two different flow zones can be distinguished: a Poiseuille-like zone, in which the velocity profile is similar to that over a rigid wall, and a boundary layer ensuring the transition with the dead zone below. It is shown that the effective boundary condition felt by the Poiseuille-like layer is essentially controlled by incoming flow characteristics, such that the thickness of this layer does not obey simple relations with cavity length. Interestingly, however, the thickness of the boundary layer appears to follow a generalized Oldroyd’s scaling with cavity length.

Accepted Sept 25, 2017. 31 pages, 21 figures.; International audience; This paper deals with the numerical resolution of a shallow water viscoplastic flow model. Viscoplastic materials are characterized by the existence of a yield stress: below a certain critical threshold in the imposed stress, there is no deformation and the material behaves like a rigid solid, but when that yield value is exceeded, the material flows like a fluid. In the context of avalanches, it means that after going down a slope, the material can stop and its free surface has a non-trivial shape, as opposed to the case of water (Newtonian fluid). The model involves variational inequalities associated with the yield threshold: finite volume schemes are used together with duality methods (namely Augmented Lagrangian and Bermudez-Moreno) to discretize the problem. To be able to accurately simulate the stopping behaviour of the avalanche, new schemes need to be designed, involving the classical notion of well-balancing. In the present context, it needs to be extended to take into account the viscoplastic nature of the material as well as general bottoms with wet/dry fronts which are encountered in geophysical geometries. Here we derive such schemes in 2D as the follow up of the companion paper treating the 1D case. Numerical tests include in particular a generalized 2D benchmark for Bingham codes (the Bingham-Couette flow with two non-zero boundary conditions on the velocity) and a simulation of the avalanche path of Taconnaz in Chamonix - Mont-Blanc to show the usability of these schemes on real topographies from digital elevation models (DEM).