Japan Geoscience Union Meeting 2015

Presentation information

Oral

Symbol S (Solid Earth Sciences) » S-VC Volcanology

[S-VC46] Dynamics of volcanic eruptions and igneous activities

Mon. May 25, 2015 9:00 AM - 10:45 AM 304 (3F)

Convener:*Tomofumi Kozono(Department of Geophysics, Graduate School of Science, Tohoku University), Yujiro Suzuki(Earthquake Research Institute, The University of Tokyo), Satoshi Okumura(Division of Earth and Planetary Materials Science, Department of Earth Science, Graduate School of Science, Tohoku University), Chair:Satoshi Okumura(Division of Earth and Planetary Materials Science, Department of Earth Science, Graduate School of Science, Tohoku University)

9:45 AM - 10:00 AM

[SVC46-11] A numerical study of pyroclastic flow dynamics: Development of a two-layer model based on Shallow-Water equations

*Hiroyuki SHIMIZU1, Takehiro KOYAGUCHI1, Yujiro SUZUKI1 (1.Earthquake Research Institute, The University of Tokyo)

Keywords:pyroclastic flows, density currents, gravity currents, Shallow-Water equations, two-layer model

Pyroclastic flows are density currents of mixtures of volcanic particles and gas that flow on the ground surface. They are characterized by wide ranges of density ratio ρ/ρa(〜100-103), where ρ and ρa are the densities of current and ambient, respectively. The dynamics of pyroclastic flows are affected by many physical factors such as settling of particles and entrainment of ambient air. In order to understand these effects on the dynamics of pyroclastic flows with wide ranges of ρ/ρa, we have developed a new numerical model based on shallow-water equations. The model has been verified by applying it to the dam-break problem, for which analytical solutions of different fluid dynamical stages are available.

In order to calculate the dynamics of pyroclastic flows with wide ranges of ρ/ρa using the shallow-water equations, the balance between the driving force and the dynamic reaction at the flow front (i.e., front condition) should be correctly taken into account. In previous works, two types of numerical models have been proposed to express the front condition: Boundary-Condition (BC) type model and Artificial-Bed (AB) type model. BC type model is a numerical method in which the front condition is applied to the boundary condition at the flow front. In AB type model, the front condition is calculated by setting a thin artificial bed ahead of the front. We have revealed that AB type model is applicable to the currents of ρ/ρa≥100, whereas BC type model should be used for the currents of ρ/ρa≤100. We have also developed a rigorous algorithm calculating the front condition using BC type model.

Our numerical method enables us to investigate the effects of the settling of particles and the entrainment of ambient air on the flow dynamics for wide ranges of ρ/ρa. Particle settling decelerates the front speed and enhances the formation of "jump" that separates the current into head and tail parts. Entrainment of ambient air also decelerates the front speed and suppresses the formation of the jump between the head and the tail of the current. Our results also suggested that the rate of entrainment tends to increase as ρ/ρa increases. When both the effects of particle settling and entrainment are present, the competition between these two effects is considered to result in diverse flow patterns.

Pyroclastic flows with large density gradients are generally divided into a dilute overriding part (100≤ρ/ρa≤101) and a dense basal part (ρ/ρa〜103). In order to simulate the dynamics of such realistic pyroclastic flows, we have developed a two-layer model where the dynamics of the dilute overriding part is solved by BC type model and that of the dense basal part is solved by AB type model. In the dilute overriding part, the effects of the settling of particles and the entrainment of ambient air are taken into account. In the dense basal part, the effects of the settling of particles and the basal friction are taken into account. We present some preliminary results showing temporal evolutions of the dilute overriding part, the dense basal part and the deposits.