An Interface Tracking Method Based on Volume Tracking in Embedded Micro Cells

Author(s):  
Akio Tomiyama ◽  
Yusuke Nakahara ◽  
Satoru Abe

An interface tracking method based on volume tracking in micro cells embedded in regular cells is proposed. The volume tracking is based on a CIP scheme, which yields extremely sharp and jaggy interface in micro cells. However, in regular cells, the jagged interface is smoothed and the interface thickness becomes comparable to a regular cell size Δx. Single air bubbles in stagnant water and single large air bubbles in stagnant water filled in a vertical pipe are simulated with the proposed method. As a result, it is confirmed that (1) the proposed method conserves bubble volumes, (2) the interface thickness around a bubble is kept within 2Δx, and thereby the surface tension force is well evaluated by the CSF model, and (3) the method can give good predictions for bubbles in a high density-ratio system.

Author(s):  
Mark Ho ◽  
Guan Heng Yeoh ◽  
John Arthur Reizes ◽  
Victoria Timchenko

Purpose Interface distinct two-phase computational fluid dynamics (CFD) simulations require accurate tracking in surface curvature, surface area and volume fraction data to precisely calculate effects such as surface tension, interphase momentum and interphase heat and mass transfer exchanges. To attain a higher level of accuracy in two-phase flow CFD simulations, the intersection marker (ISM) method was developed. The ISM method has cell-by-cell remeshing capability that is volume conservative, maintains surface continuity and is suited for the tracking of interface deformation in transient two-phase flow simulations. Studies of isothermal single bubbles rising in quiescent water were carried out to test the ISM method for two-phase flow simulations. Design/methodology/approach The ISM method is a hybrid Lagrangian–Eulerian front tracking algorithm which can model an arbitrary three-dimensional surface within an array of cubic control volumes. Fortran95 was used to implement the ISM method, which resulted in approximately 25,000+ lines of written code and comments. To demonstrate the feasibility of the ISM algorithm for two-phase flow simulations, the ISM algorithm was coupled with an in-house CFD code, which was modified to simulate two-phase flows using a single fluid formulation. The constitutional equations incorporated terms of variable density and viscosity. In addition, body force source terms were included in the momentum equation to account for surface tension and buoyancy effects. Findings The performance of two-phase flow simulations was benchmarked against experimental data for four air/water bubbles with 1, 2.5, 5 and 10 mm of diameter rising in quiescent fluid. A variety of bubble sizes were tested to demonstrate the accuracy of the ISM interface tracking method. The results attained were in close agreement with experimental observations. Practical implications The results obtained show that the ISM method is a viable means for interface tracking of two-phase flow CFD simulations. Other applications of the ISM method include simulations of solid–fluid interaction and other immersed boundary flow problems. Originality/value The ISM method is a novel approach to front tracking, and the results shown are original in content.


2007 ◽  
Author(s):  
Djemai Merrouche ◽  
Kamal Mohammedi ◽  
Idir Belaidi ◽  
Bachir Mabrouki

Author(s):  
Syed Sharif ◽  
Mark Ho ◽  
Victoria Timchenko ◽  
Guan Yeoh

Abstract In this paper, the growth of a rising vapour bubble in superheated water was numerically studied using an advanced interface tracking method, called the InterSection Marker (ISM) method. The ISM method is a hybrid Lagrangian-Eulerian Front Tracking algorithm that can model an arbitrary Three-Dimensional (3D) surface within an array of cubic control-volumes. The ISM method has cell-by-cell remeshing capability that is volume conservative, maintains surface continuity and is suited for tracking interface deformation in multiphase flow simulations. This method was previously used in adiabatic bubble rise simulation with no heat and mass transfers to or from the bubble were considered. This present work will extend the ISM method's application to simulate vapour bubble growth in superheated water with the inclusion of additional physics, such as the convective heat transfer mechanism and the phase change. Coupled with an in-house variable-density and variable-viscosity single-fluid flow solver, the method was used to simulate vapour bubble growth due to the convective action. The forces such as the surface tension and the buoyancy were included in the momentum equation. The source terms for the mass transfer were also modelled in the CFD governing equations to simulate the growth. Bubble properties such as size, shape, velocity, drag coefficient, and convective heat transfer coefficient were predicted. Effects of surface tension and temperature on the bubble characteristic were also discussed. Obtained numerical results were compared against the analytical and past works and found to be in good agreement.


Author(s):  
Geoffroy Vaudor ◽  
Alain Berlemont ◽  
Thibaut Ménard ◽  
Mathieu Doring

In this paper, a computational method is presented that addresses the problem of multiphase flow characterized by phases with significant density ratio accompanied by strong shearing. The Coupled Level-Set Volume-of-Fluid (CLSVOF) technique is used for interface tracking, while the momentum transfer is coupled to that of mass by means of momentum fluxes computed using a sub-grid. This is an extended adaptation of Rudman’s volume tracking technique [1]. The new method is shown to conserve kinetic energy when applied to cases otherwise unfeasible, such as shear layer with high density ratio.


Author(s):  
Kosuke Hayashi ◽  
Akio Tomiyama

An interface tracking method for predicting bubble dissolution process is proposed. A non-diffusive scheme for advecting species concentrations is adopted to accurately compute the volume change due to mass transfer. The applicability of the proposed method is examined through several benchmark tests, i.e. mass transfer from a static bubble and that from free rising bubbles. Predicted species concentration distributions and mass transfer coefficients agree well with theoretical and empirical models. Dissolution of single carbon dioxide bubbles in a vertical pipe filled with water is also simulated. The bubbles consist only of carbon dioxide, and nitrogen and oxygen are initially dissolved in water. The volume change due to dissolution of carbon dioxide from the bubbles and evaporation of nitrogen and oxygen from water are well predicted.


Author(s):  
Marianne M. Francois ◽  
Douglas B. Kothe ◽  
Edward D. Denby ◽  
James M. Sicilian ◽  
Matthew W. Williams

A consistent formulation is presented for modeling surface tension driven flow with the continuum surface tension force (CSF) model within a volume of fluid (VOF) method using a pressure-correction projection method. We show that a flow algorithm whose inherent design is motivated by legislating force balance gives an exact (to round off) balance between surface tension forces and pressure gradients that arise as a result. This design eliminates one of the elusive impediments to more accurate CSF-based surface tension models, the remaining of which is curvature estimation accuracy. To validate our formulation, we present results for an equilibrium (static) drop in two and three dimensions having an arbitrary density ratio and demonstrate in the process that scaling effects within the CSF framework are insignificant.


Sign in / Sign up

Export Citation Format

Share Document