Jump to content

Volume of fluid method

From Wikipedia, the free encyclopedia
An illustration of fluid simulation using VOF method.

In computational fluid dynamics, the volume of fluid (VOF) method is a family of free-surface modelling techniques, i.e. numerical techniques for tracking and locating the free surface (or fluid–fluid interface). They belong to the class of Eulerian methods which are characterized by a mesh that is either stationary or is moving in a certain prescribed manner to accommodate the evolving shape of the interface. As such, VOF methods are advection schemes capturing the shape and position of the interface, but are not standalone flow solving algorithms. The Navier–Stokes equations describing the motion of the flow have to be solved separately.

History

[edit]

The volume of fluid method is based on earlier Marker-and-cell (MAC) methods[1][2] developed at Los Alamos National Laboratory. MAC used Lagrangian marker particles to track the distribution of fluid in a fixed Eulerian grid. The use of marker particles was computationally expensive because it required many marker particles per grid cell, to reduce numerical noise when discrete marker particles move across grid cells. The original idea of the VOF method was to replace marker particles with a single scalar variable per grid cell representing the volume fraction of fluid in it. Thereby, the volume of fluid is governed by an advection equation. This idea arose from studies of two-phase mixture (water and steam) problems where it was customary to use a volume of steam variable.[3] The VOF approach was first demonstrated in a 1975 publication “Methods for Calculating Multi-Dimensional, Transient Free Surface Flows Past Bodies” by Nichols and Hirt.[4] This publication described how to advect the fluid fraction with a Donor-Acceptor scheme, how to estimate the orientation and position of the free surface inside surface cells, and how to prescribe appropriate boundary conditions (continuity and zero shear stress) at the free surface. This approach was much simpler than other techniques tracking the surface of fluid, yet more versatile as it could model the coalescence and breakup of fluid regions. In 1976, Noh & Woodward[5] presented the Simple Line Interface Calculation (SLIC), a technique to approximate fluid interfaces based on volume fractions, designed for directional-split advection scheme of volume fractions. SLIC could also handle an arbitrary number of immiscible fluid phases per grid cells. Thereby, SLIC was well suited to the VOF approach, although the two methods were initially independent and remained separate till the 90s. The term “Volume of Fluid method” and it acronym “VOF” method were coined in the 1980 Los Alamos Scientific Laboratory report, “SOLA-VOF: A Solution Algorithm for Transient Fluid Flow with Multiple Free Boundaries,” by Nichols, Hirt and Hotchkiss[6] and in the journal publication “Volume of Fluid (VOF) Method for the Dynamics of Free Boundaries” by Hirt and Nichols in 1981.[7] These two publications provided more details about the specific procedures used to approximate the position of the free surface (locally represented by an inclined line in surface cells) and apply the free surface boundary conditions on it. Since VOF method surpassed MAC by lowering computer storage requirements, it quickly became popular. Early applications of the SOLA-VOF program developed at Los Alamos include light-water-reactor safety studies.[8][9] A variant of the SOLA-VOF code was also adopted by NASA.[10][11] In 1982, Youngs developed the Piecewise-Linear Interface Calculation (PLIC) scheme,[12] which improved accuracy of interface reconstruction upon the SLIC and early VOF methods.[13]

Overview

[edit]

The method is based on the idea of a so-called fraction function . It is a scalar function, defined as the integral of a fluid's characteristic function in the control volume, namely the volume of a computational grid cell. The volume fraction of each fluid is tracked through every cell in the computational grid, while all fluids share a single set of momentum equations, i.e. one for each spatial direction. From a cell-volume averaged perspective, when a cell is empty of the tracked phase, the value of is zero; when the cell is full of tracked phase, ; and when the cell contains an interface between the tracked and non-tracked volumes, . From a perspective of a local point that contains no volume, is a discontinuous function insofar as its value jumps from 0 to 1 when the local point moves from the non-tracked to the tracked phase. The normal direction of the fluid interface is found where the value of changes most rapidly. With this method, the free-surface is not defined sharply, instead it is distributed over the height of a cell. Thus, in order to attain accurate results, local grid refinements have to be done. The refinement criterion is simple, cells with have to be refined. A method for this, known as the marker and micro-cell method, has been developed by Raad and his colleagues in 1997.[14]

The evolution of the -th fluid in a system on fluids is governed by the transport equation (actually the same equation that has to be fulfilled by the level-set method distance function ):

with the following constraint

,

i.e., the volume of the fluids is constant. For each cell, properties such as density are calculated by a volume fraction average of all fluids in the cell

These properties are then used to solve a single momentum equation through the domain, and the attained velocity field is shared among the fluids.

The VOF method is computationally friendly, as it introduces only one additional equation and thus requires minimal storage. The method is also characterized by its capability of dealing with highly non-linear problems in which the free-surface experiences sharp topological changes. By using the VOF method, one also evades the use of complicated mesh deformation algorithms used by surface-tracking methods. The major difficulty associated with the method is the smearing of the free-surface. This problem originates from excessive diffusion of the transport equation.

Discretization

[edit]

To avoid smearing of the free-surface, the transport equation has to be solved without excessive diffusion. Thus, the success of a VOF method depends heavily on the scheme used for the advection of the field. Any chosen scheme needs to cope with the fact that is discontinuous, unlike e.g. the distance function used in the Level-Set method.

Whereas a first order upwind scheme smears the interface, a downwind scheme of the same order will cause a false distribution problem which will cause erratic behavior in case of the flow is not oriented along a grid line. As these lower-order schemes are inaccurate, and higher-order schemes are unstable and induce oscillations, it has been necessary to develop schemes which keep the free-surface sharp while also producing monotonic profiles for .[15] Over the years, a multitude of different methods for treating the advection have been developed. In the original VOF-article by Hirt, a donor-acceptor scheme was employed. This scheme formed a basis for the compressive differencing schemes.

The different methods for treating VOF can be roughly divided into three categories, namely the donor-acceptor formulation, higher order differencing schemes and line techniques.

The Donor-Acceptor Schemes

[edit]

The donor-acceptor scheme is based on two fundamental criteria, namely the boundedness criterion and the availability criterion. The first one states that the value of has to be bounded between zero and one. The latter criterion ensures that the amount of fluid convected over a face during a time step is less than or equal to the amount available in the donor cell, i.e., the cell from which the fluid is flowing to the acceptor cell. In his original work, Hirt treated this with a blended scheme consisting of controlled downwinding and upwind differencing.

Higher Order Differencing Schemes

[edit]

In the higher order differencing schemes, as the name suggests, the convective transport equation is discretized with higher order or blended differencing schemes. Such methods include the Compressive Interface Capturing Scheme for Arbitrary Meshes (CICSAM) [16] and High Resolution Interface Capturing (HRIC) [17] scheme, which are both based on the Normalized Variable Diagram (NVD) by Leonard.[18]

Geometrical Reconstruction Techniques

[edit]
A spherical droplet represented by PLIC [19](Piecewise Linear Interface Calculation) geometrical reconstruction technique in a VOF simulation; (a) general view, (b) zoom into the cavity region. Reconstruction yields a planar segment in each of the control volumes; the segments are generally discontinuous, which is visible especially in the under-resolved regions.Obtained using the Basilisk code [1].

Line techniques circumvent the problems associated with the discretization of the transport equation by not tracking the interface in a cell explicitly. Instead, the fluid distribution in a cell an interface is obtained by using the volume fraction distribution of neighbouring cells. The Simple Line Interface Calculation (SLIC) by Noh and Woodward from 1976[5] uses a simple geometry to reconstruct the interface. In each cell the interface is approximated as a line parallel to one of the coordinate axes and assumes different fluid configurations for the horizontal and vertical movements respectively. A widely used technique today is the Piecewise Linear Interface Calculation by Youngs.[20] PLIC is based on the idea that the interface can be represented as a line in R2 or a plane in R3; in the latter case we may describe the interface by:

where is a vector normal to the interface. Components of the normal are found e.g. by using the finite difference method or its combination with least squares optimization. The free term is then found (analytically or by approximation) by enforcing mass conservation within computational cell. Once the description of the interface is established, the advection equation of is solved using geometrical techniques such as finding the flux of between grid cells, or advecting the endpoints of interface using discrete values of fluid velocity.

Interface capture issues

[edit]

In two-phase flows in which the properties of the two phases are vastly different, errors in the computation of the surface tension force at the interface cause Front-Capturing methods such as Volume of Fluid (VOF) and Level-Set method (LS) to develop interfacial spurious currents. To better solve such flows, special treatment is required to reduce such spurious currents. A few studies have looked at improving interface tracking by combining Level-set method and Volume of fluid methods while a few others have looked at improving the numerical solving algorithm by adding smoothening loops or improving property averaging techniques.[21]

See also

[edit]

References

[edit]
  1. ^ Harlow, Francis H.; Welch, J. Eddie (1965-12-01). "Numerical Calculation of Time-Dependent Viscous Incompressible Flow of Fluid with Free Surface". The Physics of Fluids. 8 (12): 2182–2189. doi:10.1063/1.1761178. ISSN 0031-9171.
  2. ^ Welch, J. E.; Harlow, F. H.; Shannon, J. P.; Daly, B. J. (1965-11-01). THE MAC METHOD-A COMPUTING TECHNIQUE FOR SOLVING VISCOUS, INCOMPRESSIBLE, TRANSIENT FLUID-FLOW PROBLEMS INVOLVING FREE SURFACES (Report). Los Alamos National Lab. (LANL), Los Alamos, NM (United States). doi:10.2172/4563173.
  3. ^ Hirt, C. W. "Volume of Fluid (VOF) History". FLOW-3D. Retrieved 2024-07-05.
  4. ^ Nichols, B. D.; Hirt, C. W. (1975). "Methods for calculating multidimensional, transient free surface flows past bodies". Conference: 1. international conference on numerical ship hydrodynamics, Gaithersburg, MD, USA, 20 Oct 1975. Bibcode:1975STIN...7625526N. OSTI 5044626.
  5. ^ a b Noh, W.F.; Woodward, P. (1976). van de Vooren, A.I.; Zandbergen, P.J. (eds.). SLIC (Simple Line Interface Calculation). proceedings of 5th International Conference of Fluid Dynamics. Lecture Notes in Physics. Vol. 59. pp. 330–340. doi:10.1007/3-540-08004-x_336. ISBN 3-540-08004-X.
  6. ^ Nichols, B. D.; Hirt, C. W.; Hotchkiss, R. S. (1980-08-01). SOLA-VOF: a solution algorithm for transient fluid flow with multiple free boundaries (Report). Los Alamos National Lab. (LANL), Los Alamos, NM (United States). doi:10.2172/5122053.
  7. ^ Hirt, C.W.; Nichols, B.D. (1981). "Volume of fluid (VOF) method for the dynamics of free boundaries". Journal of Computational Physics. 39 (1): 201–225. Bibcode:1981JCoPh..39..201H. doi:10.1016/0021-9991(81)90145-5.
  8. ^ Nichols, B. D.; Hirt, C. W. (February 1980). "Numerical Simulation of Boiling Water Reactor Vent-Clearing Hydrodynamics". Nuclear Science and Engineering. 73 (2): 196–209. doi:10.13182/NSE80-A18699. ISSN 0029-5639.
  9. ^ Hirt, C. W.; Nichols, B. D. (1981-05-01). "A Computational Method for Free Surface Hydrodynamics". Journal of Pressure Vessel Technology. 103 (2): 136–141. doi:10.1115/1.3263378. ISSN 0094-9930.
  10. ^ Torrey, M.; Cloutman, L. (1985). NASA-VOF2D: a computer program for incompressible with free surfaces (Technical report). Los Alamos National Laboratory. Bibcode:1985STIN...8630116T. LA-10612-MS.
  11. ^ Torrey, M. D.; Mjolsness, R. C.; Stein, L. R. (July 1987). "NASA-VOF3D: A three-dimensional computer program for incompressible flows with free surfaces". Los Alamos National Laboratory report LA-11009-MS. Bibcode:1987STIN...8810288T.
  12. ^ Youngs, D. (January 1982). Morton, K. W.; Baines, M. J. (eds.). "Time-Dependent Multi-material Flow with Large Fluid Distortion". In book: Numerical Methods in Fluid Dynamics. Academic Press.
  13. ^ Rider, William J.; Kothe, Douglas B. (1998-04-10). "Reconstructing Volume Tracking". Journal of Computational Physics. 141 (2): 112–152. doi:10.1006/jcph.1998.5906. ISSN 0021-9991.
  14. ^ Chen, S.; Raad, D.B. (1997). "The surface marker and micro-cell method". International Journal for Numerical Methods in Fluids. 25 (7): 749–778. Bibcode:1997IJNMF..25..749C. doi:10.1002/(SICI)1097-0363(19971015)25:7<749::AID-FLD584>3.3.CO;2-F.
  15. ^ Darwish, M.; Moukalled, F. (2006). "Convective Schemes for Capturing Interfaces of Free-Surface Flows on Unstructured Grids". Numerical Heat Transfer Part B. 49 (1): 19–42. Bibcode:2006NHTB...49...19D. doi:10.1080/10407790500272137. S2CID 121067159.
  16. ^ Ubbink, O.; Issa, R.I. (1999). "Method for Capturing Sharp Fluid Interfaces on ArbitraryMeshes". J. Comput. Phys. 153 (1): 26–50. Bibcode:1999JCoPh.153...26U. doi:10.1006/jcph.1999.6276.
  17. ^ Muzaferija, S.; Peric, M.; Sames, P; Schelin, T. (1998). "A two-fluid Navier-Stokes solver to simulate water entry". Twenty-Second Symposium on Naval Hydrodynamics. ISBN 978-0-309-18453-3.
  18. ^ Leonard, B.P. (1991). "The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection". Computer Methods in Applied Mechanics and Engineering. 88 (1): 17–74. Bibcode:1991CMAME..88...17L. doi:10.1016/0045-7825(91)90232-U.
  19. ^ Aniszewski, Wojciech (2014). "Volume of Fluid (VOF) type advection methods in two-phase flow: A comparative study". Computers & Fluids. 97: 52–73. arXiv:1405.5140. Bibcode:2014arXiv1405.5140A. doi:10.1016/j.compfluid.2014.03.027. S2CID 119661007.
  20. ^ Youngs, D.L. (1982). "Time-dependent multi-material flow with large fluid distortion". Numerical Methods for Fluid Dynamics. Academic Press. pp. 273–285. ISBN 978-0-12-508360-7. OCLC 9918216.
  21. ^ Rajendran, Sucharitha; Manglik, Raj M.; Jog, Milind A. (2022-06-01). "New Property Averaging Scheme for Volume of Fluid Method for Two-Phase Flows With Large Viscosity Ratios". Journal of Fluids Engineering. 144 (6): 061101. arXiv:2301.01638. doi:10.1115/1.4053548. ISSN 0098-2202. S2CID 246024850.
  • Pilliod, J.E. (1992). An analysis of Piecewise Linear Interface Reconstruction Algorithms for Volume of Fluid Methods (Thesis). University of California, Davis. OCLC 1012402545.