A pressure-based, compressible, two-phase flow finite volume method for underwater explosions

Underwater explosions present many challenges for numerical modeling. The density ratio between explosion gases and the surrounding liquid are typically O(103), and the pressure ratios can be just as high. We develop a cell-centered finite volume method to solve the governing equations of a two-phas...

Full description

Saved in:
Bibliographic Details
Published in:Computers & fluids 2013-10, Vol.87, p.132-143
Main Authors: Miller, S.T., Jasak, H., Boger, D.A., Paterson, E.G., Nedungadi, A.
Format: Article
Language:English
Subjects:
Computer simulation
Density
Explosions
Finite volume
Fluid flow
Mathematical analysis
Mathematical models
Multiphase flow
Pressure-based
Segregated solution
Underwater explosions
UNDEX
Volume fraction
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Underwater explosions present many challenges for numerical modeling. The density ratio between explosion gases and the surrounding liquid are typically O(103), and the pressure ratios can be just as high. We develop a cell-centered finite volume method to solve the governing equations of a two-phase homogeneous fluid. Isentropic and isothermal equations of state are used to relate densities and pressure, in lieu of solving an energy equation. A volume fraction is used to distinguish between the disparate phases. A pressure-based, segregated algebraic solution procedure is used to solve for the primitive variables representing volume fraction, velocity, and piezometric pressure. Numerical examples include verification and validation for a number of canonical test cases. In particular, we examine an advecting material interface and show the absence of pressure oscillations across the contact discontinuity. The non-conservative form of our equations does not guarantee exact mass balance, but numerical experiments indicate that this dissipative mechanism has a stabilizing affect on the method. Shock tube and axisymmetric underwater explosion problems are presented to demonstrate the robustness of the algorithm when very large density and pressure discontinuities are present. We simulate a shallow water explosion in three dimensions and examine the effects on the free surface and explosion bubble topology.
ISSN:0045-7930
1879-0747
DOI:10.1016/j.compfluid.2013.04.002