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

Authors: Afshin Ahmadi Nadooshan

Abstract:

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.

Keywords: Interfacial flow, Incompressible flow, surface tension, Volume-Of-Fluid, Kelvin-Helmholtz.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1081003

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.