Numerical Simulation of Interfacial Flow with Volume-Of-Fluid Method

In this article, various models of surface tension force (CSF, CSS and PCIL) for interfacial flows have been applied to dynamic case and the results were compared. We studied the Kelvin- Helmholtz instabilities, which are produced by shear at the interface between two fluids with different physical properties. The velocity inlet is defined as a sinusoidal perturbation. When gravity and surface tension are taking into account, we observe the development of the Instability for a critic value of the difference of velocity of the both fluids. The VOF Model enables to simulate Kelvin-Helmholtz Instability as dynamic case.





References:
[1] D. Gerlach, G. Tomar, G. Biswas and F. Durst ,"Comparison of volumeof-
fluid methods for surface tension-dominant two-phase flows",
International Journal of Heat and Mass Transfer, 49, 740-754, 2006
[2] Osher, S., and Sethian, J., "A fronts propagating with curvaturedependent
speed: algorithms based on Hamilton-Jacobi formulations", J.
Comp. Phys. 79(1), 12(1988)
[3] Noh, W.F. and Woodward, P.R., "Slic (simple line interface method)", in
Lecture Notes in Physics, 59, 1976.
[4] Hirt, C.W. and Nichols, B.D., "Volume of fluid (VOF) method for the
dynamics of free boundaries", J. Comp. Phys., 39, 201-225, 1981.
[5] Youngs, D.L. "Time-dependent multi-material flow with large fluid
distribution", in Numerical methods for fluid dynamics, Morton and
Norman,Editor,187-221,1996
[6] Ashgriz, N and Poo, J.Y., "FLAIR: Flux Line-segment model for
advection and interface reconstruction". J. Comp. Phys. , 93,449-468,
1991.
[7] Rider, W.J. and Kothe D.B., "Reconstruction volume tracking", J.
Comp. Phys., 14, 112, 1998.
[8] Pilliod J.E. and E.G. Puckett, "Second-order accurate volume-of-fluid
algorithms for tracking material interfaces", Lawrence Berkley Lab.
Tech. Report, No.LBNL-40744, 1997.
[9] Welch S.W.J., T. Rachidi, Numerical computation of .lm boiling
including conjugated heat transfer, Numer. Heat Transfer, Part B 42
(2002) 35-53.
[10] Agarwal D.K., S.W.J. Welch, G. Biswas, F. Durst, Planar simulation of
bubble growth in .lm boiling in near-critical water using a variant of the
VOF method, J. Heat Transfer (ASME) 126 (2004) 329-338.
[11] Renardy Y., M. Renardy, "PROST: a parabolic reconstruction of surface
tension for the volume-of-fluid method", J. Comp. Phys. 183 (2002)
400-421.
[12] Jamet D., D. Torres, J.U. Brackbill, On the theory and computation of
surface tension: the elimination of parasitic currents through energy
conservation in the second-gradient method, J. Comp. Phys. 182 (2002)
262-276.
[13] Brackbill J.U., D.B. Kothe, C. Zemach, A continuum method for
modeling surface tension, J. Comp. Phys. 100 (1992) 335-354.
[14] Kothe, D.B., W.J. Rider, S.J. Mosso, and J.S. Brock, "Volume tracking
of interfaces having surface tension in two and three dimensions" AIAA
96-0859, 1996.
[15] Lafaurie B., C. Nardone, R. Scardovelli, S. Zaleski, G. Zanetti, Modeling
merging and fragmentation in multiphase .ows with SURFER, J. Comp.
Phys. 113 (1994) 134-147.
[16] Shirani, E., Ashgriz, N. and Mostaghimi, J., "Interface pressure
calculation based on conservative of momentum for front tracking
methods", J. Comp. Phys., 203, 153-175, 2005.