竖直流道宽度对气泡运动行为影响的数值模拟

用数值方法模拟了竖直通道宽度对气泡在液体中的非定常运动、变形以及传热特性的影响。在这个模拟中,界面跟踪采用了VOF方法,并采用PL IC进行界面重构。主流场计算采用有限容积方法将控制方程离散,其中扩散项采用中心差分格式,对流项采用一阶迎风格式。用成熟的S IM PLE算法求解N-S方程的速度与压力的耦合问题。引入CSF模型处理运动界面的表面张力。利用所编制的程序计算了竖直流道中的单个气泡的形状、运动特性以及气泡内外流场与传热特性,并对竖直通道宽度在不同情况下,对气泡的形状、运动特性以及传热特性进行了进一步的研究。得到了一系列有价值的结果,并与实验结果比较。表明数值模拟结果与实验结果吻合的较好。     

A piecewise linear approach to volume tracking a triple point

An approach to volume tracking three materials is presented that, in contrast with the so‐called ‘onion‐skin’ methodology, assumes the existence of a ‘triple point’ at which two interfaces between three materials intersect. The reconstruction of any cell that contains three materials is iterative: the approach is to locate a point of intersection between two interfaces that minimizes a given error expression. The advantages and limitations of the algorithm are presented via a series of advection tests that demonstrate that triple points can be reconstructed and advected just as well as simpler interfaces in typical applications. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of sharp gas–liquid interface using VOF method and adaptive grid local refinement around the interface

The volume of fluid (VOF) method is used to perform two‐phase simulations (gas–liquid). The governing Navier–Stokes conservation equations of the flow field are numerically solved on two‐dimensional axisymmetric or three‐dimensional unstructured grids, using Cartesian velocity components, following the finite volume approximation and a pressure correction method. A new method of adaptive grid local refinement is developed in order to enhance the accuracy of the predictions, to capture the sharp gas–liquid interface and to speed up the calculations. Results are compared with experimental measurements in order to assess the efficiency of the method. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical simulation of a solitary wave interaction with submerged multi-bodies

The problems of a solitary wave passing over rectangular cylinders have been analysed. The numerical simulation is based on the full nonlinear two-dimensional Navier-Stokes equations which are solved by the finite difference method. The free surface is dealt with by the Volume of Fluid method (VOF). Results for a solitary wave passing over a single cylinder are compared with the experimental data of Seabra-Santos, Penouard and Temperville^[2] and better agreement is obtained than those obtained from the long wave equation based on the potential flow theory. Results are also given for two cylinders with different gaps. The project supported by the National Natural Science Foundation of China and the Development Foundation of Science and Technology of Shanghai Education Committee and the Royal Society.     

Interactions of three-dimensional vicous axisymmetric vortex rings with a free surface

The unsteady nonlinear interaction of three-dimensional vortices with a free surface is a great challenge in fluid mechanics, which has deep theoretical significance and important practical background. Applying the three-dimensional VOF method, the interactions of three-dimensional axisymmetric vortex rings with a free surface in an incompressible viscous fluid are numerically simulated. The influence of the Froude number and the surface tension are studied and the evolution of the vorticity, the trajectories of the vortex rings and the baroclinic vorticity on the surface are obtained. The results agreed well with the experiments reported in the literature. The project supported by the National Natural Science Foundation of China     

不同界面计算方法对液滴在电场作用下变形数值模拟精度的影响

液滴在电场作用下的变形是电流体动力学的基础课题之一,表面张力的计算精度对液滴变形量的模拟结果有重要影响。本文以开源计算流体动力学平台OpenFOAM的VOF模型为框架,研究了MULES和isoAdvector两类界面更新算法与相分数梯度和RDF函数两类曲率算法对电场作用下液滴变形模拟精度的影响。研究表明,isoAdvector算法相比MULES算法对网格密度的要求更低,但其耦合相分数梯度算法计算表面张力的误差较高。isoAdvector算法耦合RDF函数算法计算误差较低,并且在使用轴对称网格时,只有该算法能够同时处理液滴平行于电场和垂直于电场方向的变形,得到的数值结果与解析解吻合较好。     

Numerical and experimental investigation of flow over a semicircular weir

The water flow over a semicircular weir is investigated numerically and experimentally in this paper. The numerical model solves the Reynolds equation for a mean flow field with thek-ε-turbulent model. To trace the motion of the free surface, the COF method with geometric reconstruction is employed. The velocity of the flow is measured by means of LDV technique. Four types of flow patterns, the position of the separation and reattachment point, the distribution of shear stress on the bed at downstream of the weir are presented and discussed. The numerical results agree well with the experiment data.     

An improved three-dimensional model for interface pressure calculations in free-surface flows

A three-dimensional method for the calculation of interface pressure in the computational modeling of free surfaces and interfaces is developed. The methodology is based on the calculation of the pressure force at the interfacial cell faces and is mainly designed for volume of fluid (VOF) interface capturing approach. The pressure forces at the interfacial cell faces are calculated according to the pressure imposed by each fluid on the portion of the cell face that is occupied by that fluid. Special formulations for the pressure in the interfacial cells are derived for different orientations of an interface. The present method, referred to as pressure calculation based on the interface location (PCIL), is applied to both static and dynamic cases. First, a three-dimensional motionless drop of liquid in an initially stagnant fluid with no gravity force is simulated as the static case and then two different small air bubbles in water are simulated as dynamic cases. A two-fluid, piecewise linear interface calculation VOF method is used for numerical simulation of the interfacial flow. For the static case, both the continuum surface force (CSF) and the continuum surface stress (CSS) methods are used for surface tension calculations. A wide range of Ohnesorge numbers and density and viscosity ratios of the two fluids are tested. It is shown that the presence of spurious currents (artificial velocities present in case of considerable capillary forces) is mainly due to the inaccurate calculation of pressure forces in the interfacial computational cells. The PCIL model reduces the spurious currents up to more than two orders of magnitude for the cases tested.

Also for the dynamic bubble rise case, it is shown that using the numerical solver employed here, without PCIL, the magnitude of spurious currents is so high that it is not possible to simulate this type of surface tension dominated flows, while using PCIL, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.     

Multiphase modelling of pouring glass over the spout lip of an industrial float in the flat glass forming process

This paper describes the numerical modelling of the three‐phase flow that is established when pouring molten glass onto a tin bath (float) under a non‐reactive atmosphere. Discharged from the furnace, the glass is spouted over a lip onto the bath, where it floats and spreads out in the form of a long ribbon, leading to the flat glass forming process. Numerically, the physical modelling must consider a multiphase problem with three immiscible phases: glass, tin and air. The simulation of this complex free‐surface flow involves the solution of the Navier–Stokes set of equations for all the phases simultaneously, using a volume‐of‐fluid formulation that introduces a marker function convected by the flow to identify each phase. The evolution of the interphases is tracked over time with the implementation of a continuous surface force algorithm. A general purpose, well‐tested commercial code, FLUENT, is employed for the computations. Firstly, a two‐dimensional model considering the symmetry plane of the float is developed to fit accurate physical and numerical parameters. The high complexity of the interphases in addition to great differences between the physical properties of the phases has required extensive tests to ensure the consistency and accuracy of the solutions. Afterwards, a complete three‐dimensional model is built to simulate the pouring process in a real geometry and predict the behaviour of the industrial facility when the operating conditions are modified. Typical flow phenomena inside the tin bath, such as tin currents, wet back flow or even the equilibrium thickness, are also obtained as a promising result of the numerical modelling. Copyright © 2008 John Wiley & Sons, Ltd.