A wave-propagation based volume tracking method for compressible multicomponent flow in two space dimensions

We present a simple volume-of-fluid approach to interface tracking for inviscid compressible multicomponent flow problems in two space dimensions. The algorithm uses a uniform Cartesian grid with some grid cells subdivided by tracked interfaces, approximately aligned with the material interfaces in the flow field. A standard volume-moving procedure that consists of two basic steps: (1) the update of a discrete set of volume fractions from the current time to the next and (2) the reconstruction of interfaces from the resulting volume fractions, is employed to find the new location of the tracked interfaces in piecewise linear form at the end of a time step. As in the previous work by LeVeque and the author to front tracking based on a surface-moving procedure (R.J. LeVeque, K.-M. Shyue, Two-dimensional front tracking based on high-resolution wave propagation methods, J. Comput. Phys. 123 (1996) 354-368), a conservative high-resolution wave propagation method is applied on the resulting slightly non-uniform grid to update all the cell values, while the stability of the method is maintained by using a large time step idea even in the presence of small cells and the use of a time step with respect to the uniform grid cells. We validate our algorithm by performing the simulation of a Mach 1.22 shock wave in air over a circular R22 gas bubble, where sensible agreement of some key flow features of the computed solutions are observed when direct comparison of our results are made with the existing experimental and numerical ones that appear in the literature. Other computations are also presented that show the feasibility of the algorithm together with a mixture type of the model equations developed by the author (K.-M. Shyue, A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Gruneisen equation of state, J. Comput. Phys. 171 (2001) 678-707) for practical multicomponent problems with general compressible materials characterized by a Mie-Gruneisen form of the equation of state.

[1]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[2]  Ralph Menikoff,et al.  Anomalous reflection of a shock wave at a fluid interface , 1990, Journal of Fluid Mechanics.

[3]  R. LeVeque Wave Propagation Algorithms for Multidimensional Hyperbolic Systems , 1997 .

[4]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[5]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

[6]  Randall J. LeVeque,et al.  A wave propagation method for three-dimensional hyperbolic conservation laws , 2000 .

[7]  J. Craggs Applied Mathematical Sciences , 1973 .

[8]  E. Puckett,et al.  Second-Order Accurate Volume-of-Fluid Algorithms for Tracking Material Interfaces , 2013 .

[9]  M. Rudman INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

[10]  G. J. Ball,et al.  A free-Lagrange augmented Godunov method for the simulation of elastic-plastic solids , 2002 .

[11]  David H. Sharp,et al.  Numerical investigation of Richtmyer–Meshkov instability using front tracking , 1995, Journal of Fluid Mechanics.

[12]  Marius Paraschivoiu,et al.  Second order accurate volume tracking based on remapping for triangular meshes , 2003 .

[13]  P. Colella,et al.  A conservative three-dimensional Eulerian method for coupled solid-fluid shock capturing , 2002 .

[14]  Keh-Ming Shyue,et al.  Regular Article: A Fluid-Mixture Type Algorithm for Compressible Multicomponent Flow with van der Waals Equation of State , 1999 .

[15]  Randall J. LeVeque,et al.  Large time step shock-capturing techniques for scalar conservation laws , 1981 .

[16]  J. Brock,et al.  Volume tracking of interfaces having surface tension in two and three dimensions , 1996 .

[17]  Richard Saurel,et al.  A compressible flow model with capillary effects , 2005 .

[18]  Randall J. LeVeque,et al.  One-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods , 1995, SIAM J. Sci. Comput..

[19]  D. B. Kothe,et al.  A SECOND-ORDER ACCURATE, LINEARITY-PRESERVING VOLUME TRACKING ALGORITHM FOR FREE SURFACE FLOWS ON 3-D UNSTRUCTURED MESHES , 1999 .

[20]  Rainald Löhner,et al.  On the computation of multi-material flows using ALE formulation , 2004 .

[21]  Randall J. LeVeque Hyperbolic conservation laws and numerical methods , 1990 .

[22]  R. Brent Table errata: Algorithms for minimization without derivatives (Prentice-Hall, Englewood Cliffs, N. J., 1973) , 1975 .

[23]  Gretar Tryggvason,et al.  Bifurcation of tracked scalar waves , 1986 .

[24]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

[25]  J. O´Rourke,et al.  Computational Geometry in C: Arrangements , 1998 .

[26]  E. Puckett,et al.  A High-Order Godunov Method for Multiple Condensed Phases , 1996 .

[27]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[28]  A. J. Cable,et al.  High-velocity impact phenomena , 1970 .

[29]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[30]  R. LeVeque Approximate Riemann Solvers , 1992 .

[31]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[32]  S. Welch,et al.  A Volume of Fluid Based Method for Fluid Flows with Phase Change , 2000 .

[33]  Elaine S. Oran,et al.  Numerical Simulation of Reactive Flow , 1987 .

[34]  D. Fletcher,et al.  A New Volume of Fluid Advection Algorithm , 2000 .

[35]  Fotis Sotiropoulos,et al.  A Second-Order Godunov Method for Wave Problems in Coupled Solid-Water-Gas Systems , 1999 .

[36]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[37]  S. Zaleski,et al.  Volume-of-Fluid Interface Tracking with Smoothed Surface Stress Methods for Three-Dimensional Flows , 1999 .

[38]  Randall J. LeVeque,et al.  A large Time Step Generalization of Godunov’s Method for Systems of Conservation Laws , 1985 .

[39]  Mark L. Wilkins,et al.  Computer Simulation of Dynamic Phenomena , 1999 .

[40]  Keh-Ming Shyue A Simple Volume Tracking Method For Compressible Two-Phase Flow , 2001 .

[41]  Lévêque,et al.  High resolution finite volume methods on arbitrary grids via wave propagation. Final report , 1988 .

[42]  Randall J. LeVeque,et al.  H-Box Methods for the Approximation of Hyperbolic Conservation Laws on Irregular Grids , 2003, SIAM J. Numer. Anal..

[43]  Keh-Ming Shyue,et al.  A fluid-mixture type algorithm for compressible multicomponent flow with Mie-Grüneisen equation of state , 2001 .

[44]  Randall J. LeVeque,et al.  Two-Dimensional Front Tracking Based on High Resolution Wave Propagation Methods , 1996 .

[45]  Phillip Colella,et al.  A conservative front tracking method for hyperbolic conservation laws , 1987 .

[46]  Phillip Colella,et al.  A high-order Eulerian Godunov method for elastic-plastic flow in solids , 2001 .

[47]  Gert Vegter,et al.  In handbook of discrete and computational geometry , 1997 .

[48]  H. S. Udaykumar,et al.  An Eulerian method for computation of multimaterial impact with ENO shock-capturing and sharp interfaces , 2003 .

[49]  J. N. Fritz,et al.  CHAPTER VII – THE EQUATION OF STATE OF SOLIDS FROM SHOCK WAVE STUDIES , 1970 .

[50]  R. LeVeque Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[51]  John K. Dukowicz Efficient volume computation for three-dimensional hexahedral cells , 1988 .

[52]  Andrew B. Wardlaw,et al.  Spherical solutions of an underwater explosion bubble , 1998 .

[53]  A. Ruoff,et al.  Linear Shock‐Velocity‐Particle‐Velocity Relationship , 1967 .

[54]  L. Hsiao,et al.  The Riemann problem and interaction of waves in gas dynamics , 1989 .

[55]  P. Mulet,et al.  A flux-split algorithm applied to conservative models for multicomponent compressible flows , 2003 .

[56]  Phillip Colella,et al.  A Cartesian grid embedded boundary method for hyperbolic conservation laws , 2006 .

[57]  Keh-Ming Shyue Front Tracking Methods Based on Wave Propagation. , 1993 .

[58]  Hansong Tang,et al.  A Second-Order Accurate Capturing Scheme for 1D Inviscid Flows of Gas and Water with Vacuum Zones , 1996 .

[59]  M. Rudman International Journal for Numerical Methods in Fluids a Volume-tracking Method for Incompressible Multifluid Flows with Large Density Variations , 2022 .

[60]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[61]  J. Haas,et al.  Interaction of weak shock waves with cylindrical and spherical gas inhomogeneities , 1987, Journal of Fluid Mechanics.

[62]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[63]  D. Anderson,et al.  Algorithms for minimization without derivatives , 1974 .

[64]  Keh-Ming Shyue,et al.  An Efficient Shock-Capturing Algorithm for Compressible Multicomponent Problems , 1998 .

[65]  James M. Hyman,et al.  Numerical methods for tracking interfaces , 1984 .

[66]  Phillip Colella,et al.  Conservative front-tracking for inviscid compressible flow , 1991 .