In this paper, we propose a new idea of tracking material interface. Since the regions filled with different materials at the initial time are merely transported by the velocity field, the material type at present is determined by its original location. We introduce the advection equation of the base coordinates to specify the material type, and solve this equation in the Euler framework. Thanks to the initial linear distribution, this method is free of numerical diffusion for the problems with a constant or a rigid body rotation velocity field and can produce accurate results for the general case. Moreover, it is applicable to the advection function of arbitrary distribution, for example, problems with more than two types of fluids. The new method is incorporated into a newly developed flow solver employing the semi-Lagrangian model to successfully solve the flow problems with multiple types of fluids. [S0098-2202(00)02001-0]

Issue Section:
Technical Papers
Keywords:
transport processes, flow simulation, numerical analysis, interface phenomena
Topics:
Density, Diffusion (Physics), Flow (Dynamics), Fluids, Numerical analysis, Pressure, Surface tension, Dynamics (Mechanics), Errors, Deformation, Space
1.
Hirt
,
C. W.
, and
Nicholls
,
B. D.
,
1981
, “
Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries
,”
J. Comput. Phys.
,
39
, pp.
201
–
225
.
2.
Osher
,
S.
, and
Sethian
,
J. A.
,
1988
, “
Fronts Propagating with Curvature-Dependent Speed: Algorithms Based on Hamilton-Jacobi Formulations
,”
J. Comput. Phys.
,
79
, pp.
12
–
49
.
3.
Sethian, J. A., 1996, Level Set Methods: Evolving Interfaces in Geometry, Fluid Mechanics, Computer Vision, and Material Science, Cambridge University Press.
4.
Brackbill
,
J. U.
,
Kothe
,
D. B.
, and
Zemach
,
C.
,
1992
, “
A Continuum Method for Modeling Surface Tension
,”
J. Comput. Phys.
,
100
, pp.
335
–
354
.
5.
Aleinov, I., et al., 1995, “Computing Surface Tension with High-Order Kernels,” Proceedings of the 6th International Symposium on Computational Fluid Dynamics, pp. 13–18.
6.
Puckett
,
E. G.
,
Almgren
,
A. S.
,
Bell
,
J. B.
,
Marcus
,
D. L.
, and
Rider
,
W. J.
,
1997
, “
A Second-Order Projection Method for Tracking Fluid Interfaces in Variable Density Incompressible Flows
,”
J. Comput. Phys.
,
130
, pp.
269
–
282
.
7.
Rudman
,
M.
,
1997
, “
Volume Tracking Methods for Interfacial Flow Calculations
,”
IJNMF
,
24
, pp.
671
–
691
.
8.
Rider
,
W. J.
, and
Kothe
,
D. B.
,
1998
, “
Reconstructing Volume Tracking
,”
J. Comput. Phys.
,
141
, pp.
112
–
152
.
9.
Jia
,
W.
,
1998
, “
An Accurate Semi-Lagrangian Scheme Designed for Incompressible Navier-Stokes Equations Written in Generalized Coordinates
,”
Trans. Japan Soc. Aero. Space Sci.
,
41
, No.
133
, pp.
105
–
117
.
10.
Staniforth
,
A.
, and
Cote
,
J.
,
1991
, “
Semi-Lagrangian Integration Schemes for Atmospheric Models—A Review
,”
Mon. Weather Rev.
,
119
, pp.
2206
–
2223
.
11.
Hirt
,
C. W.
,
Amsden
,
A. A.
, and
Cook
,
J. L.
,
1974
, “
An Arbitrary Lagrangian Eulerian Computing Method for all Flow Speeds
,”
J. Comput. Phys.
,
14
, pp.
227
–
253
.
12.
Rider, W. J., et al., 1995, “Stretching and Tearing Interfaces Tracking Methods,” LANL Preprint.
13.
Lafaurie
,
B.
,
Nardone
,
C.
,
Scardovelli
,
R.
,
Zaleski
,
S.
, and
Zanetti
,
G.
,
1994
, “
Modelling Merging and Fragmentation in Multiphase Flows with SURFER
,”
J. Comput. Phys.
,
113
, pp.
134
–
147
.
Copyright © 2000
 by ASME
You do not currently have access to this content.