A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework

• Published in: Journal of computational physics Vol. 213; no. 1; pp. 141 - 173
• Main Authors: Francois, Marianne M., Cummins, Sharen J., Dendy, Edward D., Kothe, Douglas B., Sicilian, James M., Williams, Matthew W.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier Inc 01.03.2006; Elsevier
• Subjects: Algorithms; Balancing; Computational techniques; Continuum surface force; Curvature; Curvatures; Drop dynamics; Exact sciences and technology; Ghost fluid method; Incompressible flow; Mathematical analysis; Mathematical methods in physics; Mathematical models; Multiphase flow; Physics; Representations; Surface tension; Volume fraction; Volume-of-fluid; Curvatures; Incompressible flow; Multiphase flow; Volume-of-fluid; Ghost fluid method; Surface tension; Continuum surface force; Drop dynamics; Pressure gradients; Calculation methods; Continuum; Algorithms; Dynamics; Models; Modelling; Calculation; Curvature; Corrections
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2005.08.004

A new balanced-force algorithm is presented for modeling interfacial flow with surface tension. The algorithm is characterized by a pressure-correction method with the interfaces represented by volume fractions. Within this flow algorithm, we devise a continuous (e.g., continuum surface tension model) and a sharp (e.g., a ghost fluid method) interface representation of the surface-tension-induced interfacial pressure jump condition. The sharp interface representation is achieved by temporarily reconstructing distance functions from volume fractions. We demonstrate that a flow algorithm designed to legislate force balance retains an exact balance between surface tension forces and the resulting pressure gradients. This balance holds for both continuous and sharp representations of interfacial surface tension. The algorithm design eliminates one of the elusive impediments to more accurate models of surface tension-driven flow, the remaining of which is accurate curvature estimation. To validate our formulation, we present results for an equilibrium (static) drop in two and three dimensions having an arbitrary density jump across the interface. We find that the sharp surface tension method yields an abrupt pressure jump across the interface, whereas the continuous surface tension method results in a smoother transition. Both methods, however, yield spurious velocities of the same order, the origin of which is due solely to errors in curvature. Dynamic results are also presented to illustrate the versatility of the method.