** Progress in Earth and Planetary Science is the official journal of Japan Geoscience Union (JpGU)
** Progress in Earth and Planetary Science is partly financially supported by a Grant-in-Aid for Publication of Scientific Research Results to enhance dissemination of information of scientific research.
Gallery View of PEPS Articles
Solid earth sciences
Session convener-recommended article JpGU Meeting 2015
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.