A high‐resolution method for compressible two‐fluid flows and simulation of three‐dimensional shock–bubble interactions

Abstract

SUMMARYA high‐resolution method is developed to capture the material interfaces of compressible two‐fluid flows in multiple dimensions. A fluid mixture model system with single velocity and pressure is used, and viscous effect can also be taken into account. A consistent thermodynamic law based on the assumption of pressure equilibrium is employed to describe the thermodynamic behaviors of the pure fluids and mixture of two components. The splitting and unsplit Eulerian formulations of piecewise parabolic method are extended to numerically integrate the hyperbolic part of the model system, whereas the system of diffusion equations is solved using an explicit, central difference scheme. The block‐structured adaptive mesh refinement (AMR) capability is built in the hydrodynamic code to locally improve grid resolution. The resulting method is verified to be at least second‐order accurate in space. Numerical results show that the discontinuities, particularly contact discontinuities, can be resolved sharply. The use of AMR allows flow features at disparate scales to be resolved sufficiently. In addition, three‐dimensional shock–bubble interactions are simulated to investigate effects of Mach number on bubble evolution. The flow structures including those peculiar to three‐dimensional bubble are resolved correctly, and some physical phenomena with increasing Mach number are reported. Copyright © 2012 John Wiley & Sons, Ltd.

The piecewise parabolic method for elastic-plastic flow in solids

“…The single-component regions are then advanced in time with high-order Godunov methods i.e. piecewise parabolic method (PPM) 33 – 35 . The approach we use for elastic-plastic solid is basically identical the approach described in Zhang 25 .…”

The Effect of Shock Wave Strength on Shock Bubble Interaction

“…To solve the Equations (1)-(4), the inviscid fluxes are discretized by PPM method, while the viscous and heat fluxes are discretized by second-order central difference scheme. Using the dimension-splitting technique, the 3D problem can be simplified to 1D problem, and then two-step Lagrange-Remap algorithm is applied for time advancing [11]. About 37 million cells are adopted in the computation, and at least 120 cells are distributed in the bubble in each direction.…”

A Direct Ghost Fluid Method for Modeling Explosive Gas and Water Flows

“…e location of marker points is decided via the solution of fluid equations [15][16][17]. As a cutcell type of method, it comes with a problem that the elements cut by interfaces may contain small subcells, resulting in extremely small time step size in an explicit time integration scheme in order to meet the Courant-Friedrichs-Levy (CFL) requirement [18]. Also, the complicated front-tracking method is not always effective especially when a large interface deformation exists [18].…”