Numerical simulations of granular dam break: Comparison between discrete element, Navier-Stokes, and thin-layer models.

Granular flows occur in various contexts, including laboratory experiments, industrial processes, and natural geophysical flows. To investigate their dynamics, different kinds of physically based models have been developed. These models can be characterized by the length scale at which dynamic processes are described. Discrete models use a microscopic scale to individually model each grain, Navier-Stokes models use a mesoscopic scale to consider elementary volumes of grains, and thin-layer models use a macroscopic scale to model the dynamics of elementary columns of fluids. In each case, the derivation of the associated equations is well-known. However, few studies focus on the extent to which these modeling solutions yield mutually coherent results. In this article, we compare the simulations of a granular dam break on a horizontal or inclined planes for the discrete model convex optimization contact dynamics (COCD), the Navier-Stokes model Basilisk, and the thin-layer depth-averaged model SHALTOP. We show that, although all three models allow reproducing the temporal evolution of the free surface in the horizontal case (except for SHALTOP at the initiation), the modeled flow dynamics are significantly different, and, in particular, during the stopping phase. The stresses measured at the flow's bottom, reflecting the flow dynamics, are in relatively good agreement, but significant variations are obtained with the COCD model due to complex and fast-varying granular lattices. Similar conclusions are drawn using the same rheological parameters to model a granular dam break on an inclined plane. This comparison exercise is essential for assessing the limits and uncertainties of granular flow modeling.

Experimental assessment of the effective friction at the base of granular chute flows on a smooth incline.

We report on direct measurements of the basal force components for granular material flowing down a smooth incline. We investigate granular flows for a large range of inclination angles from θ=13.4^{∘} to 83.6° and various gate openings of the chute. We find that the effective basal friction coefficient µ_{B}, obtained from the ratio of the longitudinal force to the normal one, exhibits a systematic increase with increasing slope angle and a significant weakening with increasing particle holdup H (the depth-integrated particle volume fraction). At low angles, the basal friction is slightly less than or equal to tanθ. The deviation from tanθ can be interpreted as a contribution from the sidewall to the overall friction. At larger angles, the basal friction µ_{B} saturates at an asymptotic value that is dependent on the gate opening of the chute. Importantly, our data confirm the outcomes of recent discrete numerical simulations. First, for steady and fully developed flows as well as for moderately accelerated ones, the variation of the basal friction can be captured through a unique dimensionless number, the Froude number Fr, defined as Fr=U[over ¯]/(gHcosθ)^{1/2}, where U[over ¯] is the mean flow velocity. Second, the mean velocity scales with the particle holdup H with a power exponent close to 1/4, contrasting with the Bagnold scaling (U[over ¯]∼H^{3/2}).

Modeling of partial dome collapse of La Soufrière of Guadeloupe volcano: implications for hazard assessment and monitoring.

Over the past 9,150 years, at least 9 flank collapses have been identified in the history of La Soufrière of Guadeloupe volcano. On account of the volcano's current unrest, the possibility of such a flank collapse should not be dismissed in assessing hazards for future eruptive magmatic as well as non-magmatic scenarios. We combine morphological and geophysical data to identify seven unstable structures (volumes ranging from 1 × 106 m3 to 100 × 106 m3), including one that has a volume compatible with the last recorded flank collapse in 1530 CE. We model their dynamics and emplacement with the SHALTOP numerical model and a simple Coulomb friction law. The best-fit friction coefficient to reproduce the 1530 CE event is tan(7°) = 0.13, suggesting the transformation of the debris avalanche into a debris flow, which is confirmed by the texture of mapped deposits. Various friction angles are tested to investigate less water-rich and less mobile avalanches. The most densely populated areas of Saint-Claude and Basse-Terre, and an area of Gourbeyre south of the Palmiste ridge, are primarily exposed in the case of the more voluminous and mobile flank collapse scenarios considered. However, topography has a prominent role in controlling flow dynamics, with barrier effects and multiple channels. Classical mobility indicators, such as the Heim's ratio, are thus not adequate for a comprehensive hazard analysis.

Empirical investigation of friction weakening of terrestrial and Martian landslides using discrete element models.

Understanding what controls the travelling distance of large landslides has been the topic of considerable debate. By combining observation and experimental data with depth-averaged continuum modelling of landslides and generated seismic waves, it was empirically observed that lower effective friction had to be taken into account in the models to reproduce the dynamics and runout distance of larger volume landslides. Moreover, such simulation and observation results are compatible with a friction weakening with velocity as observed in earthquake mechanics. We investigate here as to whether similar empirical reduced friction should be put into discrete element models (DEM) to reproduce observed runout of large landslides on Earth and on Mars. First we show that, in the investigated parameter range and for a given volume, the runout distance simulated by 3D DEM is not much affected by the number (i.e. size) of grains once this number attains ~ 8000. We then calibrate the model on laboratory experiments and simulate other experiments of granular flows on inclined planes, making it possible for the first time to reproduce the observed effect of initial volume and aspect ratio on runout distances. In particular, the normalised runout distance starts to depend on the volume involved only above a critical slope angle > 16-19°, as observed experimentally. Finally, based on field data (volume, topography, deposit), we simulate a series of landslides on simplified inclined topography. The empirical friction coefficient, calibrated to reproduce the observed runout for each landslide, is shown to decrease with increasing landslide volume (or velocity), going down to values as low as 0.1-0.2. No distinguishable difference is observed between the behaviour of terrestrial and Martian landslides.

Frictional velocity-weakening in landslides on Earth and on other planetary bodies.

One of the ultimate goals in landslide hazard assessment is to predict maximum landslide extension and velocity. Despite much work, the physical processes governing energy dissipation during these natural granular flows remain uncertain. Field observations show that large landslides travel over unexpectedly long distances, suggesting low dissipation. Numerical simulations of landslides require a small friction coefficient to reproduce the extension of their deposits. Here, based on analytical and numerical solutions for granular flows constrained by remote-sensing observations, we develop a consistent method to estimate the effective friction coefficient of landslides. This method uses a constant basal friction coefficient that reproduces the first-order landslide properties. We show that friction decreases with increasing volume or, more fundamentally, with increasing sliding velocity. Inspired by frictional weakening mechanisms thought to operate during earthquakes, we propose an empirical velocity-weakening friction law under a unifying phenomenological framework applicable to small and large landslides observed on Earth and beyond.