Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 5.240
IF 5-year value: 5.768
IF 5-year
CiteScore value: 8.9
SNIP value: 1.713
IPP value: 5.53
SJR value: 3.18
Scimago H <br class='widget-line-break'>index value: 71
Scimago H
h5-index value: 51
Volume 6, issue 6
Geosci. Model Dev., 6, 1905–1924, 2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.
Geosci. Model Dev., 6, 1905–1924, 2013
© Author(s) 2013. This work is distributed under
the Creative Commons Attribution 3.0 License.

Model evaluation paper 04 Nov 2013

Model evaluation paper | 04 Nov 2013

A semi-implicit, second-order-accurate numerical model for multiphase underexpanded volcanic jets

S. Carcano1, L. Bonaventura1, T. Esposti Ongaro2, and A. Neri2 S. Carcano et al.
  • 1MOX – Modelling and Scientific Computing, Dipartimento di Matematica "F. Brioschi", Politecnico di Milano, Via Bonardi 9, 20133 Milano, Italy
  • 2Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Pisa, Via della Faggiola 32, 56126 Pisa, Italy

Abstract. An improved version of the PDAC (Pyroclastic Dispersal Analysis Code, Esposti Ongaro et al., 2007) numerical model for the simulation of multiphase volcanic flows is presented and validated for the simulation of multiphase volcanic jets in supersonic regimes. The present version of PDAC includes second-order time- and space discretizations and fully multidimensional advection discretizations in order to reduce numerical diffusion and enhance the accuracy of the original model. The model is tested on the problem of jet decompression in both two and three dimensions. For homogeneous jets, numerical results are consistent with experimental results at the laboratory scale (Lewis and Carlson, 1964). For nonequilibrium gas–particle jets, we consider monodisperse and bidisperse mixtures, and we quantify nonequilibrium effects in terms of the ratio between the particle relaxation time and a characteristic jet timescale. For coarse particles and low particle load, numerical simulations well reproduce laboratory experiments and numerical simulations carried out with an Eulerian–Lagrangian model (Sommerfeld, 1993). At the volcanic scale, we consider steady-state conditions associated with the development of Vulcanian and sub-Plinian eruptions. For the finest particles produced in these regimes, we demonstrate that the solid phase is in mechanical and thermal equilibrium with the gas phase and that the jet decompression structure is well described by a pseudogas model (Ogden et al., 2008). Coarse particles, on the other hand, display significant nonequilibrium effects, which associated with their larger relaxation time. Deviations from the equilibrium regime, with maximum velocity and temperature differences on the order of 150 m s−1 and 80 K across shock waves, occur especially during the rapid acceleration phases, and are able to modify substantially the jet dynamics with respect to the homogeneous case.

Publications Copernicus