for treating surface tension under the framework of volume of fluid(VOF)
the exact location of fluid interface cannot be easily determined. In order to resolve such a problem
a new sub-grid interface locating algorithm is proposed in this work. In the new algorithm
a variable curvature curve
used to approximate the interface in sub-grid
is constructed based on the curvature of the two grid cells on both sides of the interface. As a result
the intersection point of the connection line of centers of interfacial cells and the interface is calculated accurately. The influence of the curvature on the interface location is taken into account. Moreover
this algorithm is independent on mesh type and space dimension. Numerical tests
including capillary rise
static droplet
bubble rise in two-dimensions and three-dimensions
show that interface curvature has a significant influence on the sub-grid interface location. Calculating the sub-grid interface location and the curvature at such a location with high accuracy can improve the numerical prediction of VOF method significantly
and by 3.5 times at the best. This method can provide a new tool for predicting physical behaviors of high surface tension driven flow in complex geometrical structures.
关键词
Keywords
references
UNVERDI S O, TRYGGVASON G. A front-tracking method for viscous, incompressible, multi-fluid flows [J]. Journal of Computational Physics, 1992, 100(1): 25-37.
SUSSMAN M, SMEREKA P, OSHER S. A Level Set approach for computing solutions to incompressible two-phase flow [J]. Journal of Computational Physics, 1994, 114(1): 146-159.
HIRTC W, NICHOLSB D. Volume of fluid(VOF)method for the dynamics of free boundaries [J]. Journal of Computational Physics, 1981, 39(1): 201-225.
COQUERELLE M, GLOCKNER S. A fourth-order accurate curvature computation in a level set framework of two-phase flows subjected to surface tension forces [J]. Journal of Computational Physics, 2016, 305(1): 838-876.
YOUNGS D L. An interface tracking method for a 3D Eulerian hydrodynamics code: technical report 44/92/35 [R]. Berks, UK: AWRE, 1984: 3-9.
EDWARD P J J, GERRY E P. Second-order accurate volume-of-fluid algorithms for tracking material interfaces [J]. Journal of Computational Physics, 2004, 199(2): 465-502.
SCARDOVELLI R, ZALESKI S. Interface reconstruction with least-square fit and split Eulerian-Lagrangian advection [J]. Int J Numer Meth Fluids, 2003, 41(3): 251-274.
SUSSMAN M, GERRY E P. A coupled level set and volume-of-fluid method for computing 3d and axisymmetric incompressible two-phase flows [J]. Journal of Computational Physics, 2000, 162(2): 301-337.
BALCÁZAR N, ORIOL L, JOFRE L, et al. A coupled volume-of-fluid/level-set method for simulation of two-phase flows on unstructured meshes [J]. Computers Fluids, 2016, 124(1): 12-29.
SCARDOVELLI R, ZALESKI R. Direct numerical simulation of free-surface and interfacial flow [J]. Annual Review of Fluid Mechanics, 1999, 31(1): 567-603.
DIOT S, FRANCOIS M M. An interface reconstruction method based on an analytical formula for 3D arbitrary convex cells [J]. Journal of Computational Physics, 2016, 305(15): 63-74.
LPEZ J, HERNÁNDEZ J, GMEZ P, et al. A new volume conservation enforcement method for PLIC reconstruction in general convex grids [J]. Journal of Computational Physics, 2016, 316(1): 338-359.
FRANCOIS M M, CUMMINS S J, DENDY E D, et al. A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework [J]. Journal of Computational Physics, 2006, 213(1): 141-173.
ISSA R I. Solution of the implicitly discretized fluid flow equations by operator-splitting [J]. Journal of Computational Physics, 1985, 62(1): 40-65.
SU J, Gu Z L, CHEN C G, et al. A two-layer mesh method for discrete element simulation of gas-particle systems with arbitrarily polyhedral mesh [J]. International Journal for Numerical Methods in Engineering, 2015, 103(10): 759-780.
ADAMI S, HU X Y, ADAMS N A. A new surface-tension formulation for multi-phase SPH using a reproducing divergence approximation [J]. Journal of Computational Physics, 2010, 229(13): 5011-5021.
SUN Pengnan, LI Yunbo, MING Furen. Numerical simulation on the motion characteristics of freely rising bubbles using smoothed particle hydrodynamics method [J]. Acta Phys Sin, 2015, 64(17): 1747011-1747015.
WALTERS J K, DAVIDSON J F. The initial motion of a gas bubble formed in an inviscid liquid: Part 1 The two-dimensional bubble [J]. Journal of Fluid Mechanics, 1962, 12(3): 408-416.
BHAGA D, WEBER M E. Bubbles in viscous liquids: shapes, wakes, and velocities [J]. Journal of Fluid Mechanics, 1981, 105(1): 61-85.
Study of the Mechanism of Electrostatic Accumulation and Current Density Distribution Characteristics During Liquid Hydrogen Pipeline Transportation
Numerical Simulation of Space Charge Distribution Characteristics of Polyimide Materials under X-Ray Irradiation
Numerical Simulation Study on Geothermal Extraction Characteristics of Supercritical CO2 in Heterogeneous Fractures
Numerical Study on the Effect of Relative Wettability on Binary Droplet Collision
Action Mechanism of Heat Transfer between the Working Medium and the Impeller on Aerodynamic Performance and Flow Characteristics of Radial Inflow Turbines with Supercritical Carbon Dioxide
Related Author
SUN Wenhao
LIU Hongbao
WANG Lei
QING Ziyou
LI Zhuolun
MA Yuan
LI Yanzhong
ZHONG Hui
Related Institution
School of Energy and Power Engineering,Xi’an Jiaotong University
Aerospace System Engineering Shanghai
School of Electrical Engineering,Xi’an Jiaotong University
Northwest Institute of Nuclear Technology
Institute of Electronic Engineering,China Academy of Engineering Physics