Advanced

High-order surface tension VOF-model for 3D bubble flows with high density ratio

Lörstad, Daniel LU and Fuchs, Laszlo LU (2004) In Journal of Computational Physics 200(1). p.153-176
Abstract
An improved Volume of Fluid (VOF) method is presented which is applicable to high density ratio 3D flows for a large range of bubble Reynolds number (Re). The method is based on the Navier-Stokes equations for incompressible multi-phase flows which are discretized on a Cartesian staggered grid. The multi-grid technique together with the pressure-velocity coupling scheme for multi-phase flows have resulted in an efficient solver which nearly exponentially converge with the number of iterations. The convergence speed also shows negligible dependence on density ratio, viscosity ratio and Re. A second-order accurate, non-diffusive, mass conservative phase transport model is presented which does not suffer from unphysical over- or under-shoots... (More)
An improved Volume of Fluid (VOF) method is presented which is applicable to high density ratio 3D flows for a large range of bubble Reynolds number (Re). The method is based on the Navier-Stokes equations for incompressible multi-phase flows which are discretized on a Cartesian staggered grid. The multi-grid technique together with the pressure-velocity coupling scheme for multi-phase flows have resulted in an efficient solver which nearly exponentially converge with the number of iterations. The convergence speed also shows negligible dependence on density ratio, viscosity ratio and Re. A second-order accurate, non-diffusive, mass conservative phase transport model is presented which does not suffer from unphysical over- or under-shoots of the phase variable. The high accurate normal, curvature and surface tension force model in combination with the high-order defect-correction scheme for multi-phase flows shows second-order global accuracy when applied to the transient bubble rise where the viscosity ratio is equal to one. In contrast, the commonly used viscosity model for VOF introduces a first order error for the same problem. The VOF method has been tested for different types of bubble flows at low Re and for path-oscillating and wobbling air bubbles (in water) with a diameter range of 1.82 < D < 6 mm. The numerical results agree quantitatively with the available experimental data. The investigations show that the proposed high accurate surface tension model can be used successfully for wobbling flows with bubble deformation while maintaining the mass of the phases. The error in mass conservation is directly proportional to the residual in solving the discrete problem. (C) 2004 Elsevier Inc. All rights reserved. (Less)
Please use this url to cite or link to this publication:
author
organization
publishing date
type
Contribution to journal
publication status
published
subject
keywords
surface tension, volume of fluid, bubble two-phase flows, curvature
in
Journal of Computational Physics
volume
200
issue
1
pages
153 - 176
publisher
Elsevier
external identifiers
  • wos:000224284500008
  • scopus:4544285126
ISSN
0021-9991
DOI
10.1016/j.jcp.2004.04.001
language
English
LU publication?
yes
id
19c48bf5-e7e7-45ca-9b9d-b0f19f83715d (old id 265353)
date added to LUP
2007-10-24 20:47:03
date last changed
2017-08-20 03:39:15
@article{19c48bf5-e7e7-45ca-9b9d-b0f19f83715d,
  abstract     = {An improved Volume of Fluid (VOF) method is presented which is applicable to high density ratio 3D flows for a large range of bubble Reynolds number (Re). The method is based on the Navier-Stokes equations for incompressible multi-phase flows which are discretized on a Cartesian staggered grid. The multi-grid technique together with the pressure-velocity coupling scheme for multi-phase flows have resulted in an efficient solver which nearly exponentially converge with the number of iterations. The convergence speed also shows negligible dependence on density ratio, viscosity ratio and Re. A second-order accurate, non-diffusive, mass conservative phase transport model is presented which does not suffer from unphysical over- or under-shoots of the phase variable. The high accurate normal, curvature and surface tension force model in combination with the high-order defect-correction scheme for multi-phase flows shows second-order global accuracy when applied to the transient bubble rise where the viscosity ratio is equal to one. In contrast, the commonly used viscosity model for VOF introduces a first order error for the same problem. The VOF method has been tested for different types of bubble flows at low Re and for path-oscillating and wobbling air bubbles (in water) with a diameter range of 1.82 &lt; D &lt; 6 mm. The numerical results agree quantitatively with the available experimental data. The investigations show that the proposed high accurate surface tension model can be used successfully for wobbling flows with bubble deformation while maintaining the mass of the phases. The error in mass conservation is directly proportional to the residual in solving the discrete problem. (C) 2004 Elsevier Inc. All rights reserved.},
  author       = {Lörstad, Daniel and Fuchs, Laszlo},
  issn         = {0021-9991},
  keyword      = {surface tension,volume of fluid,bubble two-phase flows,curvature},
  language     = {eng},
  number       = {1},
  pages        = {153--176},
  publisher    = {Elsevier},
  series       = {Journal of Computational Physics},
  title        = {High-order surface tension VOF-model for 3D bubble flows with high density ratio},
  url          = {http://dx.doi.org/10.1016/j.jcp.2004.04.001},
  volume       = {200},
  year         = {2004},
}