Advection of axisymmetric interfaces by the volume‐of‐fluid method

SUMMARY A criterion is proposed for the advection of axisymmetric interfaces. The location of an interface is followed by a volume-tracking technique wherein a volume fraction parameter is assigned to each of the cells in a Eulerian grid system. The interface is discretized into a set of line segments fitted at the boundary of every pair of neighburing computational cells. The orientation of a line segment is obtained by inspecting the volume fractions of two neighbouring cells. The volume fractions are then advected using the velocity components at the boundary of the two cells. The following advection criterion is proposed: for advection in the axial direction the axial velocity I( is assumed constant in the vicinity of each cell face; for advection in the radial direction the radial velocity v times the radial distance r is assumed constant in the vicinity of each cell face, i.e. gv = const., where fi = 0 for Cartesian and fi = 1 for axisymmetric systems. The above criterion is used to develop an algorithm for the advection of axisymmetric interfaces which is referred to as the ‘axisymmetric flux line segment model for advection and interface reconstruction’ or A-FLAIR.

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

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

[3]  Ahmed F. Ghoniem,et al.  Formation and inflammation of a turbulent jet , 1984 .

[4]  E. Puckett A Numerical Study of Shock Wave Refraction at a CO2/CH4 Interface , 1989 .

[5]  John A. Trapp,et al.  A numerical technique for low-speed homogeneous two-phase flow with sharp interfaces , 1976 .

[6]  Jerzy M. Floryan,et al.  Numerical Methods for Viscous Flows With Moving Boundaries , 1989 .

[7]  Alexandre J. Chorin,et al.  Curvature and solidification , 1985 .

[8]  Heinz W. Engl,et al.  Distribution of event times in time-resolved fluorescence: the exponential series approach—algorithm, regularization, analysis , 1991 .

[9]  P. Woodward,et al.  SLIC (Simple Line Interface Calculation) , 1976 .

[10]  C. W. Hirt,et al.  SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries , 1980 .

[11]  P. K. Barr,et al.  Interface scheme for turbulent flame propagation , 1984 .

[12]  J. Sethian Curvature and the evolution of fronts , 1985 .

[13]  A. K. Oppenheim,et al.  Numerical modeling of turbulent combustion in premixed gases , 1981 .

[14]  A. Chorin Flame advection and propagation algorithms , 1980 .

[15]  L. F. Henderson,et al.  On the refraction of shock waves at a slow–fast gas interface , 1991, Journal of Fluid Mechanics.

[16]  K. Hain The partial donor cell method , 1978 .