Volume of fluid (VOF) method for the dynamics of free boundaries

doi:10.1016/0021-9991(81)90145-5

Journal of Computational Physics

Volume 39, Issue 1, January 1981, Pages 201-225

https://doi.org/10.1016/0021-9991(81)90145-5 Get rights and content

Abstract

Several methods have been previously used to approximate free boundaries in finite-difference numerical simulations. A simple, but powerful, method is described that is based on the concept of a fractional volume of fluid (VOF). This method is shown to be more flexible and efficient than other methods for treating complicated free boundary configurations. To illustrate the method, a description is given for an incompressible hydrodynamics code, SOLA-VOF, that uses the VOF technique to track free fluid surfaces.

References (22)

C.W. Hirt et al.
J. Comput. Phys.
(1974)
B.D. Nichols et al.
J. Comput. Phys.
(1973)
B.D. Nichols et al.
J. Comput. Phys.
(1971)
W.E. Johnson
Development and Application of Computer Programs Related to Hypervelocity Impact
Systems Science, and Software report 3SR-353
(1970)
SOLA Codes may be obtained from the National Energy Software Center, Argonne National Laboratory 9700 South Cass...
J.D. Kershner et al.
2DE: A Two-Dimensional Continuous Eulerian Hydrodynamic Code for Computing Multicomponent Reactive Hydrodynamic Problems
Los Alamos Scientific Laboratory report LA-4846
(1972)
J.C. Martin et al.
Philos. Trans. Roy. Soc. London Ser. A
(1952)
C.W. Hirt et al.
SOLA-A Numerical Solution Algorithm for Transient Fluid Flows
Los Alamos Scientific Laboratory report LA-5852
(1975)
B.D. Nichols et al.
F.H. Harlow et al.
Phys. Fluids
(1965)
J.E. Welch et al.
The MAC Method: A Computing Technique for Solving Viscous, Incompressible, Transient Fluid Flow Problems Involving Free Surfaces
Los Alamos Scientific Laboratory report LA-3425
(1966)

F.H. Harlow et al.

J. Comput. Phys.

(1976)

Cited by (13946)

A powder-scale multiphysics framework for powder bed fusion of fiber-reinforced polymer composites
2024, Advanced Powder Materials
Additive manufacturing of fiber-reinforced polymer composites has garnered great interest due to its potential in fabricating functional products with lightweight characteristics and unique material properties. However, the major concern in polymer composites remains the presence of pore defects, as a thorough understanding of pore formation is insufficient. In this study, a powder-scale multiphysics framework has been developed to simulate the printing process of fiber-reinforced polymer composites in powder bed fusion additive manufacturing. This numerical framework involves various multiphysics phenomena such as particle flow dynamics of fiber-reinforced polymer composite powder, infrared laser–particle interaction, heat transfer, and multiphase fluid flow dynamics. The melt depths of one-layer glass fiber–reinforced polyamide 12 composite parts fabricated by selective laser sintering are measured to validate modelling predictions. The numerical framework is employed to conduct an in-depth investigation of pore formation mechanisms within printed composites. Our simulation results suggest that an increasing fiber weight fraction would lead to a lower densification rate, larger porosity, and lower pore sphericity in the composites.
Interfacial analysis and development of condensation models with OpenFOAM
2024, International Journal of Thermal Sciences
Interfacial mass and heat transfer are crucial factors in phase-change models. In this study, based on the open source toolbox OpenFOAM, an extended volume-of-fluid (VOF) solver with phase transition inserted is utilized to study three typical phase-change models (Lee model, energy balance model, and Rattner model). Through simulations and analysis of phase-change problems, we found that these three models exhibit different sensitivities to gas thermal conductivity, particularly through the interfacial heat conduction equation. The temperature and heat source distributions at the vapor-liquid interface in each model were analyzed using analytical problems, such as the one-dimensional Stefan problem and the Nusselt film condensation problem. Based on the analysis, we proposed a new energy balance model that considers the ratio of vapor-liquid thermal conductivity to correct the impact of vapor thermal conductivity. Furthermore, we introduced a temperature correction term to stabilize the interfacial temperature at the saturation temperature. The simulation results of the proposed model demonstrated excellent agreement with the analytical solutions of quasi-static condensation processes, including the one-dimensional Stefan problem and the Nusselt film condensation problem. Throughout the dynamic condensation process with rising bubbles, different phase-change models exhibited significant differences in condensing power distribution over space and time.
Temporal evolution of demarcations between the unsteady and steady regions and between the film and meniscus during dip coating with vertical withdrawal of Newtonian fluids
2024, Chemical Engineering Science
Thus far studies on dip coating focus on the asymptotic layer thickness which is reached only in the lower part of the layer when the withdrawal time is sufficiently long. This work attempts to offer a complete description of the evolution of layer thickness profile with withdrawal time for vertical withdrawal of Newtonian fluids using simulation and the Buckingham π theorem. The results show that over time, the withdrawal process successively undergoes two regimes: (I) an unsteady film region and an unsteady meniscus region and (II) an unsteady film region connected to a steady meniscus region through an unsteady meniscus region. The controlling factors for the two regimes are identified at different values of capillary number and Reynolds number. Models are proposed to predict the heights of the demarcations between different regions. The models agree well with experimental results when the controlling forces are viscous force and gravity.
Numerical investigation of the water-entry impact performance of a bionic unmanned aerial-underwater vehicle
2024, Ocean Engineering
Numerical investigation of the water-entry impact performance of a bionic unmanned aerial–underwater vehicle (UAUV) is performed. According to the modified SST turbulent model method and volume of fluid (VOF) model, a robust solution procedure for a numerical method is described, and a verification study is conducted on water-entry cases of a cone, hemisphere, and gannet. Biological observation of flying fish has led to the proposal of this three-dimensional modelling approach for water entry processes of bionic UAUVs. The influence of the water-entry velocity and angle on the impact characteristics of a bionic UAUV is numerically investigated by this proposed numerical method. The impact pattern for different water-entry states is used to evaluate the correlation between impact acceleration and the water-entry state. The results indicate the feasibility of investigating the water-entry impact performance of bionic UAUVs in detail. The maximum axial impact acceleration is less than 2.58 g when the water-entry angle is θ ≤ 30° and the velocity is V ≤ 10 m/s. A quadratic relationship exists between the maximum impact acceleration and the initial water-entry velocity, and an approximate cubic relationship exists between the maximum impact acceleration and water-entry angle. The UAUV water entry at too large an angle and too high a velocity should be avoided to avoid a significant impact load.The impact rule investigation of numerical simulation can provide a feasible analysis strategy and a bionic efficient design configuration for UAUV.
On the influence of flow-front orientation on stringer stiffened composite panels in water impacts
2024, Ocean Engineering
Water impacts form the critical load case for high-performance carbon fibre reinforced polymer (CFRP) racing craft. Such events produce a peaked, non-uniform pressure distribution that travels along a hull panel as it is immersed. Current design standards are based on static, uniform pressure loads that do not account for the directional nature of water impacts. With recent trends towards the use of directionally stiffened hull structures in the form of stringer stiffened composite panels (SSCPs), such simplifications of the load case may no longer be valid. In this study, a marine-based SSCP was tested experimentally and numerically to investigate the effects of flow-front orientation on high-performance hull panels. Parallel and perpendicular impacts at constant velocity were carried out using the novel Servo-hydraulic Slam Testing System (SSTS) and the results were used to validate a one-way coupled computational fluid dynamics — finite element analysis (CFD-FEA) Fluent/Abaqus solution. The highest strains in the monolithic skin and stringer were observed for perpendicular impacts. A parameter sweep across a range of impact orientations between parallel and perpendicular impacts was carried out. An approximately linear relationship between flow orientation angle and key structural strains was observed, with the highest strains reported at 70–90°. Results indicate that the critical load case for SSCPs occurs at 75°orientation angles, where strains in the stringer capping are maximum.
Experimental and numerical study on the adjustment of twin hydrofoils of a variable-structure SWATH
2024, Ocean Engineering
The form-changeable unmanned surface vehicles (USVs) have garnered considerable research interest due to their wide-ranging applicability and robust combat capabilities. This paper specifically addresses the longitudinal motion stability of a variable-structure Small Waterplane Area Twin-hull Ship (SWATH) equipped with twin hydrofoils, aiming to evaluate the impact of hydrofoils on sailing attitudes. Initially, model tests were systematically conducted under different form-states, both with and without hydrofoils, in calm water conditions. Subsequently, verification and validation of computational fluid dynamics (CFD) simulations were performed against the experimental results. Finally, strategies to mitigate bow-diving during high-speed tests were proposed based on the validated CFD simulations, focusing on adjusting hydrofoil attack angles. The test results indicate that with outwardly spreading struts, the bow tends to dip down more as speed increases. Model tests, combined with CFD simulations, serve as a reference for stable sailing at full speeds, presenting strategies for adjusting attack angles under varying form-states and speeds. Moreover, effective foil adjustment modes are identified to suppress bow-diving, considering thrust. Strategies for adjusting attack angles during the acceleration process are summarized, verifying the stability of the optimal 10° form-state at high speeds.

View all citing articles on Scopus

^☆: The U. S. Government's right to retain a nonexclusive royalty-free license in and to the copyright covering this paper, for governmental purposes, is acknowledged. This work was performed under the auspices of the U. S. Department of Energy.

View full text

Volume of fluid (VOF) method for the dynamics of free boundaries☆

Abstract

J. Comput. Phys.

J. Comput. Phys.

J. Comput. Phys.

Systems Science, and Software report 3SR-353

Los Alamos Scientific Laboratory report LA-4846

Philos. Trans. Roy. Soc. London Ser. A

SOLA-A Numerical Solution Algorithm for Transient Fluid Flows

Los Alamos Scientific Laboratory report LA-5852

Phys. Fluids

The MAC Method: A Computing Technique for Solving Viscous, Incompressible, Transient Fluid Flow Problems Involving Free Surfaces

Los Alamos Scientific Laboratory report LA-3425

J. Comput. Phys.