skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Semi-implicit method to solve compressible multiphase fluid flows without acoustic time step restrictions

Journal Article · · Computers and Fluids

A new 3D multiphase numerical capability is presented here for simulating multiphase flow regimes at all Mach numbers (M). The new method is a semi-implicit extension of the finite volume discrete equation method (DEM) of Chinnaya et al. (2004), which originally used explicit time-stepping. The capability is also developed to work with another extension of the DEM to moving grids for arbitrary Lagrangian-Eularian (ALE) methods detailed in Dunn (2011). Rather than solving all phase equations simultaneously, the DEM reduces the equations to a system of single-phase Riemann solves, where each phase has its own velocity and thermodynamic state. Exchanges between the phases are modeled through source terms accounting for the phase interactions. Furthermore, since the original multiphase scheme uses an explicit time-advancement scheme, it has time step restrictions dictated by the speed of sound, which limits the model's ability to simulate weakly compressible flows which typically need to be integrated for longer time periods. Here, we extend the current multiphase formulation by implementing a pressure-correcting step to enable implicit calculations and remove acoustic time step limitations. The new semi-implicit algorithm allows use of relatively large time steps compared to an explicit method. Validation and benefits of the new implicit time-step method are illustrated using several examples including weakly compressible flows and strong shock waves. The scheme presented here is general and may be used for a variety of applications which require capabilities for handling multiphase flow at a wide range of Mach numbers. However, the goal of this research is to simulate all stages of high energy explosions, including the shock formation (high Mach numbers) and evolution of the buoyant cloud (low Mach numbers).

Research Organization:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Organization:
USDOE National Nuclear Security Administration (NNSA); USDOE Laboratory Directed Research and Development (LDRD) Program
Grant/Contract Number:
AC52-07NA27344; ASC52-07NA27344
OSTI ID:
1765291
Alternate ID(s):
OSTI ID: 1682507
Report Number(s):
LLNL-JRNL-797139; 998624; TRN: US2206232
Journal Information:
Computers and Fluids, Vol. 210, Issue no. 15; ISSN 0045-7930
Publisher:
ElsevierCopyright Statement
Country of Publication:
United States
Language:
English

References (19)

Numerical solution of the Navier-Stokes equations journal January 1968
An improved reconstruction method for compressible flows with low Mach number features journal May 2008
Adaptive solution techniques for simulating underwater explosions and implosions journal January 2008
An Euler-Lagrange particle approach for modeling fragments accelerated by explosive detonation: EULER-LAGRANGE PARTICLE APPROACH FOR MODELING FRAGMENTS
  • Price, Matthew A.; Nguyen, Vinh-Tan; Hassan, Oubay
  • International Journal for Numerical Methods in Engineering, Vol. 106, Issue 11 https://doi.org/10.1002/nme.5155
journal December 2015
A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows journal April 1999
Chapter 7 Pyroclastic surges and compressible two-phase flow journal January 1998
Compressible, multiphase semi-implicit method with moment of fluid interface representation journal December 2014
Advanced multi-domain method for multiphase flow interaction with Lagrangian structural meshes journal July 2019
An analysis and comparison of the time accuracy of fractional-step methods for the Navier-Stokes equations on staggered grids journal January 2002
The Compressibility of Media under Extreme Pressures journal September 1944
Explicit Gibbs free energy equation of state applied to the carbon phase diagram journal April 2000
Algorithm for non-hydrostatic dynamics in the Regional Oceanic Modeling System journal January 2007
Numerical simulation of cloud rise phenomena associated with nuclear bursts journal October 2009
An artificially upstream flux vector splitting scheme for the Euler equations journal July 2003
Approximate factorization as a high order splitting for the implicit incompressible flow equations journal October 1992
A method for avoiding the acoustic time step restriction in compressible flow journal June 2009
Modelling detonation waves in heterogeneous energetic materials journal May 2004
A second-order projection method for the incompressible navier-stokes equations journal December 1989
Multifluid flow calculations at all mach numbers journal September 1974