Computing granular avalanches and landslides

Geophysical mass flows—debris flows, volcanic avalanches, landslides—are often initiated by volcanic activity. These flows can contain O(10⁶-10⁷) m³ or more of material, typically soil and rock fragments that might range from centimeters to meters in size, are typically O(10 m) deep, and can run out over distances of tens of kilometers. This vast range of scales, the rheology of the geological material under consideration, and the presence of interstitial fluid in the moving mass, all make for a complicated modeling and computing problem. Although we lack a full understanding of how mass flows are initiated, there is a growing body of computational and modeling research whose goal is to understand the flow processes, once the motion of a geologic mass of material is initiated. This paper describes one effort to develop a tool set for simulations of geophysical mass flows. We present a computing environment that incorporates topographical data in order to generate a numerical grid on which a parallel, adaptive mesh Godunov solver can simulate model systems of equations that contain no interstitial fluid. The computational solver is flexible, and can be changed to allow for more complex material models, as warranted.