** Progress in Earth and Planetary Science is the official journal of the Japan Geoscience Union, published in collaboration with its 50 society members.

>>Japan Geoscience Union

>>Links to 50 society members

• • • • • Progress in Earth and Planetary Science

## Gallery View of PEPS Articles

Research

Solid earth sciences

Session convener-recommended article JpGU Meeting 2015

201703201703

### A numerical shallow-water model for gravity currents for a wide range of density differences

#### Shimizu H A, Koyaguchi T, Suzuki Y J

Gravity currents, Numerical model, Shallow-water model, Front condition Numerical results of a two-layer pyroclastic density current model. The non-dimensional high-density current height (hH) and low-density current height (hL + hH) at non-dimensional time t = 1 (a) and 10 (b) are shown.

Gravity currents with various contrasting densities play a role in mass transport in a number of geophysical situations. The ratio of the density of the current, ρc, to the density of the ambient fluid, ρa, can vary between 100 and 103. In this paper, we present a numerical method of simulating gravity currents for a wide range of ρc / ρa using a shallow-water model. In the model, the effects of varying ρc / ρa are taken into account via the front condition (i.e., factors describing the balance between the driving pressure and the ambient resistance pressure at the flow front). Previously, two types of numerical models have been proposed to solve the front condition. These are referred to here as the Boundary Condition (BC) model and the Artificial Bed (AB) model. The front condition is calculated as a boundary condition at each time step in the BC model, whereas it is calculated by setting a thin artificial bed ahead of the front in the AB model. We assessed the BC and AB models by comparing their numerical results with the analytical results for a simple case of homogeneous currents. The results from the BC model agree well with the analytical results when ρc / ρa ≲ 102, but the model tends to overestimate the speed of the front position when ρc / ρa ≳ 102. In contrast, the AB model generates good approximations of the analytical results for ρc / ρa ≳ 102, given a sufficiently small artificial bed thickness, but fails to reproduce the analytical results when ρc / ρa ≲ 102. Therefore, we propose a numerical method in which the BC model is used for currents with ρc / ρa ≲ 102 and the AB model is used for currents with ρc / ρa ≳ 102.