Large-eddy simulation of pyroclastic density currents

Abstract

We investigate the dynamics of turbulent pyroclastic density currents (PDCs) by adopting a 3D, Eulerian-Eulerian multiphase flow model, in which solid particles are treated as a continuum and the grain-size distribution is simplified by assuming two particulate phases. The turbulent sub-grid stress of the gas phase is modelled within the framework of Large-Eddy Simulation (LES) by means of a eddy-viscosity model together with a wall closure. Despite the significant numerical diffusion associated to the upwind method adopted for the Finite-Volume discretization, numerical simulations demonstrate the need of adopting a Sub-Grid Scale (SGS) model, while revealing the complex interplay between the grid and the SGS filter sizes. We also analyse the relationship between the averaged flow dynamic pressure and the action exerted by the PDC on a cubic obstacle, to evaluate the impact of a PDC on a building. Numerical results suggest that the average flow dynamic pressure can be used as a proxy for the force per unit surface acting on the building envelope (Fig. 5), even for such steeply stratified flows. However, it is not possible to express such proportionality as a constant coefficient such as the drag coefficient in a steady-state current. The present results indeed indicate that the large epistemic and aleatory uncertainty on initial and boundary conditions has an impact on the numerical predictions which is comparable to that of grid resolution.