An application of parallel computing to the simulation of volcanic eruptions

Abstract

A parallel code for the simulation of the transient 3D dispersal of volcanic particles produced by explosive eruptions is presented. The model transport equations, based on the multiphase flow theory, describe the atmospheric dynamics of the gas-particle mixture ejected through the volcanic crater. The numerics is based on a finite-volume discretization scheme and a pressure-based iterative non-linear solver suited to compressible multiphase flows. The code has been parallelized by adopting an ad hoc domain partitioning scheme that enforces the load balancing. An optimized communication layer has been built over the Message-Passing Interface. The code proved to be remarkably efficient on several high-performance platforms and makes it possible to simulate fully 3D eruptive scenarios on realistic volcano topography.