An Improved Volume of Fluid Method for Two-Phase Flow Computations on Collocated Grid System

2011 ◽  
Vol 133 (4) ◽  
Author(s):  
Dong-Liang Sun ◽  
Yong-Ping Yang ◽  
Jin-Liang Xu ◽  
Wen-Quan Tao

An improved volume of fluid method called the accurate density and viscosity volume of fluid (ADV-VOF) method is proposed to solve two-phase flow problems. The method has the following features: (1) All operations are performed on a collocated grid system. (2) The piecewise linear interface calculation is used to capture interfaces and perform accurate estimations of cell-edged density and viscosity. (3) The conservative Navier–Stokes equations are solved with the convective term discretized by a second and third order interpolation for convection scheme. (4) A fractional-step method is applied to solve the conservative Navier–Stokes equations, and the BiCGSTAB algorithm is used to solve the algebraic equations by discretizing the pressure-correction equation. The above features guarantee a simple, stable, efficient, and accurate simulation of two-phase flow problems. The effectiveness of the ADV-VOF method is verified by comparing it with the conventional volume of fluid method with rough treatment of cell-edged density and viscosity. It is found that the ADV-VOF method could successfully model the two-phase problems with large density ratio and viscosity ratio between two phases and is better than the conventional volume of fluid method in this respect.

2013 ◽  
Vol 71 ◽  
pp. 283-296 ◽  
Author(s):  
Sudarshan Tiwari ◽  
Axel Klar ◽  
Steffen Hardt ◽  
Alexander Donkov

2011 ◽  
Vol 1 (32) ◽  
pp. 17 ◽  
Author(s):  
Ivo Wenneker ◽  
Peter Wellens ◽  
Reynald Gervelas

ComFLOW is a 3D Volume-of-Fluid (VOF) model to solve the incompressible Navier-Stokes equations including free surface, or to solve the Navier-Stokes equations for two-phase flow problems (two-phase flow: both an incompressible viscous fluid (e.g., water) and a compressible viscous fluid (e.g., air) are present). The problem statement of the present study reads: ‘Is ComFLOW capable of accurate prediction of wave impacts on (impermeable) coastal structures such as dikes? And, if so, what are the preferred model settings and associated computing times?’. In this paper, ComFLOW is validated for this purpose by comparison against pressure data as measured in the Delta flume by pressure sensors at dikes. We have selected three different experiments, with typical dike geometries (slope 1:3.5, with and without berm) at which more than 20 pressure sensors were installed. The results can be summarized as follows. The pressure measurements are reproduced well in the simulations. A grid with about 170 grid cells per wave length in the horizontal, and between 4 and 6 grid cells per wave height in the vertical, proves to be sufficiently fine. At such a grid resolution and with about 450 by 35 grid cells in the computational domain, a typical CPU time is 35 minutes for simulations with a model time of 10 wave periods. For the present application, it is preferable to use the one-phase flow model rather than the two-phase flow model, since the former gives better results in the lower located pressure sensors and consumes less CPU time.


Water ◽  
2021 ◽  
Vol 13 (21) ◽  
pp. 3092
Author(s):  
Lourenço Sassetti Mendes ◽  
Javier L. Lara ◽  
Maria Teresa Viseu

Spillway design is key to the effective and safe operation of dams. Typically, the flow is characterized by high velocity, high levels of turbulence, and aeration. In the last two decades, advances in computational fluid dynamics (CFD) made available several numerical tools to aid hydraulic structures engineers. The most frequent approach is to solve the Reynolds-averaged Navier–Stokes equations using an Euler type model combined with the volume-of-fluid (VoF) method. Regardless of a few applications, the complete two-phase Euler is still considered to demand exorbitant computational resources. An assessment is performed in a spillway offset aerator, comparing the two-phase volume-of-fluid (TPVoF) with the complete two-phase Euler (CTPE). Both models are included in the OpenFOAM® toolbox. As expected, the TPVoF results depend highly on the mesh, not showing convergence in the maximum chute bottom pressure and the lower-nappe aeration, tending to null aeration as resolution increases. The CTPE combined with the k–ω SST Sato turbulence model exhibits the most accurate results and mesh convergence in the lower-nappe aeration. Surprisingly, intermediate mesh resolutions are sufficient to surpass the TPVoF performance with reasonable calculation efforts. Moreover, compressibility, flow bulking, and several entrained air effects in the flow are comprehended. Despite not reproducing all aspects of the flow with acceptable accuracy, the complete two-phase Euler demonstrated an efficient cost-benefit performance and high value in spillway aerated flows. Nonetheless, further developments are expected to enhance the efficiency and stability of this model.


2018 ◽  
Vol 52 (6) ◽  
pp. 2357-2408 ◽  
Author(s):  
Stefan Metzger

We construct a Galerkin finite element method for the numerical approximation of weak solutions to a recent micro-macro bead-spring model for two-phase flow of dilute polymeric solutions derived by methods from nonequilibrium thermodynamics ([Grün, Metzger, M3AS 26 (2016) 823–866]). The model consists of Cahn-Hilliard type equations describing the evolution of the fluids and the unsteady incompressible Navier-Stokes equations in a bounded domain in two or three spatial dimensions for the velocity and the pressure of the fluids with an elastic extra-stress tensor on the right-hand side in the momentum equation which originates from the presence of dissolved polymer chains. The polymers are modeled by dumbbells subjected to a finitely extensible, nonlinear elastic (FENE) spring-force potential. Their density and orientation are described by a Fokker-Planck type parabolic equation with a center-of-mass diffusion term. We perform a rigorous passage to the limit as the spatial and temporal discretization parameters simultaneously tend to zero, and show that a subsequence of these finite element approximations converges towards a weak solution of the coupled Cahn-Hilliard-Navier-Stokes-Fokker-Planck system. To underline the practicality of the presented scheme, we provide simulations of oscillating dilute polymeric droplets and compare their oscillatory behaviour to the one of Newtonian droplets.


2020 ◽  
Vol 10 (23) ◽  
pp. 8739
Author(s):  
Vitaly Sergeev ◽  
Nikolai Vatin ◽  
Evgeny Kotov ◽  
Darya Nemova ◽  
Svyatoslav Khorobrov

The main objective of the study is to propose a technical solution integrated into the pipeline for the transition of the flow regime from slug to bubbly two-phase flow. The object of research is isothermal two-phase gas–Newtonian-liquid flow in a horizontal circular pipeline. There is local resistance in the pipe in the form of a streamlined transverse mesh partition. The mesh partition ensures the transition of the flow from the slug regime to the bubbly regime. The purpose of the study is to propose a technical solution integrated into the pipeline for changing the flow regime of a two-phase flow from slug to bubbly flow. The method of research is a simulation using computational fluid dynamics (CFD) numerical simulation. The Navier–Stokes equations averaged by Reynolds describes the fluid motion. The k-ε models were used to close the Reynolds-averaged Navier–Stokes (RANS) equations. The computing cluster «Polytechnic—RSK Tornado» was used to solve the tasks. The results of simulation show that pressure drop on the grid did not exceed 10% of the pressure drop along the length of the pipeline. The mesh partition transits the flow regime from slug to layered one, which will help to increase the service life and operational safety of a real pipeline at insignificant energy costs to overcome the additional resistance integrated into the pipeline.


Sign in / Sign up

Export Citation Format

Share Document