New developments in the analysis of column-collapse pyroclastic density currents through numerical simulations of multiphase flows

Abstract

Abstract.A granular multiphase model has been used to evaluate the action of differently sized particles on the dynamics of fountains and associated pyroclastic density currents. The model takes into account the overall disequilibrium conditions between a gas phase and several solid phases, each characterized by its own physical properties. The dynamics of the granular flows (fountains and pyroclastic density currents) has been simulated by adopting a Reynolds-averaged Navier-Stokes model for describing the turbulence effects. Numerical simulations have been carried out by using different values for the eruptive column temperature at the vent, solid particle frictional concentration, turbulent kinetic energy, and dissipation. The results obtained provide evidence of the multiphase nature of the model and describe several disequilibrium effects. The low concentration (≤5 × 10 −4 ) zones lie in the upper part of the granular flow, above the fountain, and above the tail and body of pyroclastic density current as thermal plumes. The high concentration zones, on the contrary, lie in the fountain and at the base of the current. Hence, pyroclastic density currents are assimilated to granular flows constituted by a low concentration suspension flowing above a high concentration basal layer (boundary layer), from the proximal regions to the distal ones. Interactions among the solid particles in the boundary layer of the granular flow are controlled by collisions between particles, whereas the dispersal of particles in the suspension is determined by the dragging of the gas phase. The simulations describe well the dynamics of a tractive boundary layer leading to the formation of stratified facies during Strombolian to Plinian eruptions.

Density Stratification and Buoyancy Evolution in Pyroclastic Density Currents

“…PDCs that result from collapsing eruption columns are thought to be highly density stratified (Breard & Lube, 2017; Breard et al., 2016; Brown & Branney, 2013; Lepore & Scarpati, 2012; Valentine, 1987). Density stratification is thought to hinder turbulence, mixing, and even particle settling within PDCs and other sediment‐laden gravity currents (Calder et al., 1997; Cantero et al., 2009; Davarpanah Jazi & Wells, 2016; Ellison & Turner, 1959; Felix et al., 2005; Fisher, 1979; Hanenkamp, 2011; Jacobson & Testik, 2014; Kneller & Buckee, 2000; Meiburg & Kneller, 2010; Migeon et al., 2012; Riley et al., 2008; Ruddick et al., 1989; Sulpizio et al., 2007), which impedes lift‐off relative to less stratified or non‐stratified currents.…”