A parallel multiphase flow code for the 3D simulation of explosive volcanic eruptions

A new parallel code for the simulation of the transient, 3D dispersal of volcanic particles in the atmosphere is presented. The model equations, describing the multiphase flow dynamics of gas and solid pyroclasts ejected from the volcanic vent during explosive eruptions, are solved by a finite-volume discretization scheme and a pressure-based iterative non-linear solver suited to compressible multiphase flows. The solution of the multiphase equation set is computationally so demanding that the simulation of the transient 3D dynamics of eruptive columns would not be cost-effective on a single workstation. The new code has been parallelized by adopting an ad hoc domain partitioning scheme that enforces the load balancing in the presence of a large number of topographic blocking-cells. An optimized communication layer has been built over the Message-Passing Interface. It is shown that the present code has a remarkable efficiency on several high-performance platforms and makes it possible, for the first time, to simulate fully 3D eruptive scenarios on realistic volcano topography.

[1]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[2]  Peter S. Pacheco.  Parallel programming with MPI , 1996 .

[3]  B. Voight,et al.  Computational modelling of the transient dynamics of the August 1997 Vulcanian explosions at Soufrière Hills Volcano, Montserrat: influence of initial conduit conditions on near-vent pyroclastic dispersal , 2002, Geological Society, London, Memoirs.

[4]  Augusto Neri,et al.  Numerical simulation of collapsing volcanic columns with particles of two sizes , 1996 .

[5]  M. Rosi,et al.  The 1631 Vesuvius eruption. A reconstruction based on historical and stratigraphical data , 1993 .

[6]  A. A. Amsden,et al.  A numerical fluid dynamics calculation method for all flow speeds , 1971 .

[7]  B. Voight,et al.  The eruption of Soufrière Hills Volcano, Montserrat from 1995 to 1999 , 2002 .

[8]  Andrew W. Woods,et al.  Explosive Volcanic Eruptions , 2001 .

[9]  Maria Vittoria Salvetti,et al.  An immersed boundary method for compressible multiphase flows: application to the dynamics of pyroclastic density currents , 2007 .

[10]  L. Wilson,et al.  THEORETICAL MODELING OF THE GENERATION, MOVEMENT, AND EMPLACEMENT , 1978 .

[11]  Cornelis Vuik,et al.  Segregated solution methods for compressible flow , 2000 .

[12]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[13]  Herbert E. Huppert,et al.  Particle-driven gravity currents: asymptotic and box model solutions , 2000 .

[14]  Cornelis Vuik,et al.  EFFICIENT COMPUTATION OF FLOW WITH CAVITATION BY COMPRESSIBLE PRESSURE CORRECTION , 2000 .

[15]  A. Neri,et al.  Pyroclastic flow hazard assessment at Vesuvius (Italy) by using numerical modeling. II. Analysis of flow variables , 2002 .

[16]  Loyce Adams,et al.  New Parallel SOR Method by Domain Partitioning , 1999, SIAM J. Sci. Comput..

[17]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[18]  Chia-Jung Hsu.  Numerical Heat Transfer and Fluid Flow , 1981 .

[19]  M. Syamlal,et al.  MFIX documentation theory guide , 1993 .

[20]  William I. Rose,et al.  Numerical modeling of geophysical granular flows: 2. Computer simulations of plinian clouds and pyroclastic flows and surges , 2004 .

[21]  S. Dartevelle,et al.  Numerical modeling of geophysical granular flows: 1. A comprehensive approach to granular rheologies and geophysical multiphase flows , 2004 .

[22]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[23]  G.E. Moore,et al.  Cramming More Components Onto Integrated Circuits , 1998, Proceedings of the IEEE.

[24]  emontmej,et al.  High Performance Computing , 2003, Lecture Notes in Computer Science.

[25]  Pieter Wesseling,et al.  Unified methods for computing incompressible and compressible flow , 2009 .

[26]  H. Huppert,et al.  Emplacement of the Taupo ignimbrite by a dilute turbulent flow , 1996, Nature.

[27]  Dimitri Gidaspow,et al.  Multiparticle simulation of collapsing volcanic columns and pyroclastic flow , 2003 .

[28]  Paolo Papale,et al.  The role of magma composition and water content in explosive eruptions: 1. Conduit ascent dynamics , 1998 .

[29]  Greg A. Valentine,et al.  Numerical models of Plinian eruption columns and pyroclastic flows , 1989 .

[30]  Torrey,et al.  K-FIX: a computer program for transient, two-dimensional, two-fluid flow. THREED: an extension of the K-FIX code for three-dimensional calculations. [BWR; PWR] , 1976 .

[31]  D. Gidaspow.  Multiphase Flow and Fluidization , 1994 .

[32]  Hans-Ulrich Schmincke,et al.  Explosive Volcanic Eruptions , 1984 .

[33]  C. Merkle,et al.  Computation of Multiphase Mixture Flows with Compressibility Effects , 2002 .

[34]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .

[35]  Paolo Papale,et al.  The role of magma composition and water content in explosive eruptions: 2. Pyroclastic dispersion dynamics , 1998 .

[36]  R. B. Mason,et al.  THE DEVELOPMENT AND DECAY OF THE 100-MB. SUMMERTIME ANTICYCLONE OVER SOUTHERN ASIA , 1963 .

[37]  Augusto Neri,et al.  NUMERICAL SIMULATION OF COLLAPSING VOLCANIC COLUMNS , 1993 .

[38]  M. Rosi,et al.  From Magma to Tephra: Modelling Physical Processes of Explosive Volcanic Eruptions , 1998 .

[39]  B. Voight,et al.  A new high-performance 3D multiphase flow code to simulate volcanic blasts and pyroclastic density currents: example from the Boxing Day event, Montserrat , 2005 .

[40]  A. Woods,et al.  The dynamics and thermodynamics of large ash flows , 1996 .

[41]  Charles Meneveau,et al.  Generalized Smagorinsky model for anisotropic grids , 1993 .

[42]  Hester Bijl,et al.  A Unified Method for Computing Incompressible and Compressible Flows in Boundary-Fitted Coordinates , 1998 .

[43]  Jungwoo Kim,et al.  An immersed-boundary finite-volume method for simulations of flow in complex geometries , 2001 .

[44]  Michael Herzog,et al.  Volcanic plume simulation on large scales , 1998 .

[45]  Lionel Wilson,et al.  Explosive Volcanic Eruptions–III. Plinian Eruption Columns , 1958 .

[46]  K. Wohletz.  Chapter 7 Pyroclastic surges and compressible two-phase flow , 1998 .

[47]  Izumi Hachisu,et al.  A numerical study of turbulent mixing in eruption clouds using a three-dimensional fluid dynamics model , 2005 .

[48]  A. A. Amsden,et al.  Numerical calculation of multiphase fluid flow , 1975 .

[49]  A. Neri,et al.  Pyroclastic flow hazard assessment at Vesuvius (Italy) by using numerical modeling. I. Large-scale dynamics , 2002 .

[50]  W. Shyy,et al.  A Pressure-Based Method for Turbulent Cavitating Flow Computations , 2001 .

[51]  R. Santacroce,et al.  Un GIS per il Vesuvio , 2000 .

[52]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[53]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[54]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .