The Direct Numerical Simulation of the Rising Gas Bubble With the Volume of Fluid (VOF) Method

Author(s):  
Shuji Hironaka ◽  
Saki Manabe ◽  
Yuki Fujisawa ◽  
Gen Inoue ◽  
Yosuke Matsukuma ◽  
...  

A gas-liquid two phase flow is complicated and it has not been understood well thus far, in spite of extensive investigation. Numerical simulation is a potential approach to understand this phenomenon. Although a number of studies have been conducted to understand the behavior of bubbles on the basis of computational fluid dynamics (CFD), it is difficult to completely simulate a complicated three-phase flow, including coalescence and breakup of bubbles. Although the two-fluid model based on the semi-empirical model can well estimate the actual behavior of the system in which the equations are derived, the estimation over the applicable region of equations does not always agree with the actual result. Since the 1960s, various procedures have been proposed to directly track the free surface between two phases, for example, the adaptive mesh method and the particle method. Although each of these methods has certain advantages and disadvantages, the volume of fluid (VOF) method is the most acceptable method for capturing the free surface accurately and clearly. However, a concern related to this method is the maintenance of a constant volume of the fluid. In this study, a simulation code using the VOF method is developed in order to estimate the behavior of bubbles in a vertical pipe. Further, an offset of the volume fraction is introduced to stably calculate and minimize the volume fluctuation. The effect of the surface tension is also built into the program in order to estimate the behavior of the bubbles rising through the liquid medium. The simulations of the collapsed water column and a single rising bubble are conducted with the proposed simulation code. Consequently, we confirm that these results fairly agree with the experimental ones.

Author(s):  
Y. G. Chen ◽  
W. G. Price ◽  
P. Temarel

This investigation continues the development of an anti-diffusive volume of fluid method [1] by improving accuracy through the addition of an artificial diffusion term, with a negative diffusion coefficient, to the original advection equation describing the evolution of the fluid volume fraction. The advection and diffusion processes are split into a set of two partial differential equations (PDEs). The improved anti-diffusive Volume of Fluid (VOF) method is coupled with a two-fluid flow solver to predict free surface flows and illustrated by examples given in two-dimensional flows. The first numerical example is a solitary wave travelling in a tank. The second example is a plunging wave generated by flow over a submerged obstacle of prescribed shape on a horizontal floor. The computational results are validated against available experimental data.


2008 ◽  
Vol 14 (2) ◽  
pp. 72-87 ◽  
Author(s):  
Koustuv Debnath ◽  
Amartya Kumar Bhattacharya ◽  
Biswanath Mahato ◽  
Agnimitro Chakrabarti

Author(s):  
Marco Pellegrini ◽  
Giulia Agostinelli ◽  
Hidetoshi Okada ◽  
Masanori Naitoh

Steam condensation is characterized by a relatively large interfacial region between gas and liquid which, in computational fluid dynamic (CFD) analyses, allows the creation of a discretized domain whose average cell size is larger than the interface itself. For this reason generally one fluid model with interface tracking (e.g. volume of fluid method, VOF) is employed for its solution in CFD, since the solution of the interface requires a reasonable amount of cells, reducing the modeling efforts. However, for some particular condensation applications, requiring the computation of long transients or the steam ejected through a large number of holes, one-fluid model becomes computationally too expensive for providing engineering information, and a two-fluid model (i.e. Eulerian two-phase flow) is preferable. Eulerian two-phase flow requires the introduction of closure terms representing the interactions between the two fluids in particular, in the condensation case, drag and heat transfer. Both terms involve the description of the interaction area whose definition is different from the typical one adopted in the boiling analyses. In the present work a simple but effective formulation for the interaction area is given based on the volume fraction gradient and then applied to a validation test case of steam bubbling in various subcooling conditions. It has been shown that this method gives realistic values of bubble detachment time, bubble penetration for the cases of interest in the nuclear application and in the particular application to the Fukushima Daiichi accident.


2014 ◽  
Vol 700 ◽  
pp. 643-646
Author(s):  
Dong Wang ◽  
Si Qing Zhang ◽  
Yun Long Zhang

In order to investigate the silt abrasion of modified trailing edge of stay vane in Francis turbine, the numerical simulation of trailing edge with different geometries were carried out based on the solid-liquid two-phase flow by means of Computation Fluid Dynamics. The results show that low solid volume fraction distributes on the chamfered surface of trailing edge, and high solid volume fraction distributes on the end of oblique surface. The smaller the modified angle is, the larger the distribution area of high solid volume fraction is, which show the trailing edge with smaller oblique angle may suffer from silt abrasion. Therefore, in order to solve the vibration caused by Karman vortex the trailing edge has to be sharpened, the oblique angle of trailing edge should not be too small. At end of trailing edge needs to ensure a certain thickness, especially the trailing edge near the lower ring can be thicker, which can meet the anti-abrasion requirements.


Author(s):  
Rudolf Schilling ◽  
Moritz Frobenius

The numerical simulations of three types of two-phase flow in centrifugal pump impellers are described. First, the liquid-solid particle flow is modeled by an Euler-Lagrangeian approach assuming a mass concentration less than 5% and particle diameters being less than 1000 microns. The empirical erosion model to predict the local and total wear is calibrated by measurements. Second, the influence of the relative air contents on the head-drop is simulated assuming a relatively small volume fraction and applying a simple one-fluid model. The mixture is characterized by a common density depending on the flow field. Finally, the cavitating flow is studied by implementing the Rayleigh equation into the numerical procedure describing the transient process of bubble growth and collapse. The developed simulation tools are applied to predict the three types of two-phase flows in impellers. Within the defined ranges of application the simulation results agree fairly well with the experimental data.


Author(s):  
Stamatis Kalogerakos ◽  
Mustapha Gourma ◽  
Chris Thompson

Severe limitations of the use of three-dimensional computational fluid dynamics codes (CFD) arise when trying to simulate multiphase flow in long pipes due to time constraints. 1D codes for two-phase flow, based on two-fluid models, are fast but are known to be accurate only when the velocities are within the Kelvin-Helmholtz inviscid limit [1]. An alternative is to carry out a two-dimensional CFD simulation of a channel based on the Volume of Fluid (VOF) model. 2D CFD has a wider applicability range compared to 1D, it does not have the issue of ill-posedness and it also has better turbulence models built in. Again compared to 1D the 2D VOF model has a better interface description and wall treatment. In this paper a novel method is introduced that allows swift simulations of pipeline two-phase flow in the stratified and slug flow regime, by approximating the pipe as a channel and with a methodology that solves the problem of the interfacial velocity differences, inherent in the volume of fluid model. An initial validation using the wave growth problem has already been carried out [2]. Here a set consisting of 92 experimental cases in the slug flow regime has been simulated with 2D CFD, and the simulation results showed a good agreement with experimental results. Discussions in the paper include also the question of the range of applicability for 2D CFD, and the advantages and disadvantages compared to 3D CFD and also to 1D code based on the two-fluid model. Shear stresses are then extracted from the 2D CFD simulations and used to recalibrate the friction factors [3] used in the 1D code.


2019 ◽  
Vol 195 ◽  
pp. 820-827 ◽  
Author(s):  
X. Yin ◽  
I. Zarikos ◽  
N.K. Karadimitriou ◽  
A. Raoof ◽  
S.M. Hassanizadeh

Sign in / Sign up

Export Citation Format

Share Document