利用格子Boltzmann方法模拟单个气泡在复杂流道内的运动特性

利用格子Boltzmann方法模拟了单个气泡在具有三个半圆形喉部的复杂流道内的上升过程.通过分析气泡运动过程中的形态、运动轨迹及运动速度的变化,研究复杂流道对气泡运动特性的影响.研究结果表明,在上升过程中,由于壁面的影响,气泡的形状发生严重的变形,运动轨迹也发生相应的偏转.通过实验结果的对比,验证了模拟结果的正确性.结果表明格子Boltzmann方法可以用于模拟具有复杂边界的两相流问题.     

气泡在液体中运动过程的数值模拟

本文用数值方法预测气泡在液体中的百定常运动。运用位标函数进行界面的隐含跟踪并且与有限体积法相结合构成一种可行的计算方法。本文方法允许在界面处存在很大的物性差，而且较容易将表面张力引入控制方程。我们对气液两相流中单个气泡的运动进行了计算，得到了与实验结果符合很好的数值结果。     

水平布置等径气泡上升运动数值模拟分析

采用FTM(Front Tracking Method)方法,使用Fortran语言自编程序,对气泡上升运动进行数值模拟,并且选用小气泡作为研究对象,忽略气泡变形对气泡运动的影响。通过对比分析1~4个气泡的流场与运动速度变化,发现气泡间相互作用主要受流场变化的影响,由于双气泡的漩涡抵消和尾迹合并作用,双气泡的终极速度要比单气泡小,而三个气泡和四个气泡向上运动时中间气泡超越两侧气泡之后逐渐表现为单气泡和双气泡的运动特性,但发展到终极速度的时间更长。     

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

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

基于室内实验的海底渗漏甲烷气泡运动行为及吸附回收研究

海底甲烷渗漏可能会影响海洋环境,乃至影响全球气候和碳循环,但目前人们对渗漏甲烷在海水中的运动行为了解有限。本文基于水体甲烷释放室内实验,研究了不同甲烷渗漏条件下的甲烷气泡行为及其运动特征。结果显示,甲烷气泡在上浮过程中发生合并、分离和破碎,运动轨迹呈“S”形;气泡尺寸及上浮速度随释放气体流量增大而增大;气泡运动还受水中障碍物影响,具体作用与释放气体流量、障碍物表面粗糙度及形态有关。采用日本水合物试采数据开展甲烷喷发式释放模拟实验发现,短期释放大量甲烷会引起水流速度和动压力显著增大。此外,水体甲烷吸附实验表明,经白炭黑疏水处理后,活性炭吸附甲烷能力可提高5%。     

浮力气泡对水平壁面的回弹动力学特性

黏性液体中的气泡浮升运动有趣而又复杂,而气泡与固壁边界的相互作用更是广泛存在于实际工程中.基于轴对称数值计算,模拟了浮力驱动下气泡在黏性液体中上升并与顶部水平固壁面碰撞、回弹的过程.采用考虑表面张力的不可压、变密度Navier-Stokes方程来描述气液两相流流动,并通过基于分级八叉树的有限体积法进行数值求解.为准确捕捉气泡在回弹过程中局部而迅速的拓扑变化,采用了动态自适应网格技术耦合流体体积法(volume of fluid,VOF)来重构气泡的形状. 从气泡对壁面的碰撞和回弹的基本现象入手,研究了伽利略数 Ga和接触速度$U_{a}$对气泡回弹动力学特性的影响, 分析了气泡碰撞过程中涡结构的变化.用回弹高度$H$、回弹周期$T$、长宽比{$A_{r}$}、浮升速度$U$、轴向位置$z$和回复系数$C_{r}$等参数来表征不同条件时气泡的运动和形状特性. 研究结果表明,气泡的回弹运动特性对 Ga十分敏感. Ga的增大可加剧气泡形变, 促进气泡的回弹运动, 增多回弹次数,增大回弹参数($T$和$H)$, 提升回复系数. 然而,接触速度并非决定气泡回弹动力学的控制参数, $U_{a}$的改变并不会改变回复系数.      

三维气泡与自由表面相互作用的直接数值模拟

采用VOF中的PLIC界面重构方法数值模拟了三维气泡在液流中上升并与自由表面相互作用的运动.分别考察了不同初始高度,有无来流及有无再生气泡对气泡上升高度、上升速度、压力及与自由表面相互作用等的影响.结果表明:气泡初始位置越低,顶端上升的高度越大,自由面隆起的范围更广.越靠近自由表面,底部射流横向发展越窄,而向上的压力梯度,气泡上升速度,底部射流上升高度越大,反之则反;但如果底部射流均在接近自由表面以前已横向发展充分,则差别不大.气泡外形、上升高度、破裂时间以及上升速度与来流无关.产生再生气泡后,原生气泡与再生气泡相吸,相互加速对方的上升;自由表面抬升的高度增幅较大,范围拓宽,上升速度也大大增加,且再生气泡越多,自由表面隆起的范围越大.     

水下爆炸气泡运动的理论研究

在适当深度的无黏、无旋的流体中对水下爆炸气泡运动特性进行理论研究。综合运用势流理论、能量方程以及拉格朗日方程建立气泡在不可压缩流体中的运动方程。并以此为基础,考虑重力、浮力以及阻力等多种因素对气泡运动特性的影响,通过引入新的边界积分方程,结合分析力学中完整非保守系统的Hamilton原理建立气泡在可压缩流体中的运动微分方程,并对微分方程进行求解。将方程的数值解与MSC.DYTRAN非线性有限元软件的计算结果以及经验公式进行对比,方程数值解与二者都具有较好的一致性。结果表明,基于非保守系统可压缩流体建立的气泡运动方程正确、可行,相关的理论研究和计算具有一定参考价值。     

单个上升气泡与异质液膜相互作用的数值模拟

为了探讨在流体区域内气泡上升然后与异质液膜接触的过程中,气泡速度、形状随时间的变化,以及区域内液膜的运动特性。本文采用FTM方法模拟研究了在2D×D的空间范围内单个气泡从静止开始上升到与异质液膜接触过程中的运动特性。研究结果表明:气泡上升到与液膜一定距离后,会使液膜中间凸起变薄,并将液膜挤向两侧,气泡的竖直速度先增大后减小;液膜厚度不同时,气泡上升初期相同时刻的纵横比?基本相同;气泡上升后期相同时刻的纵横比随着液膜厚度的增厚而增大;不同表面张力下,相同时刻气泡纵横比随着表面张力的减小而减小;表面张力越小,气泡顶部和底端的压力差会越小,上升速度越小。     

气泡夹带对壁面油膜流动特性的影响

油-气润滑系统工作过程中,润滑油膜受微油滴冲击和压缩空气扰动影响易形成气泡夹带现象,气泡夹带行为将对壁面润滑油膜层的形成及流动过程产生重要影响。基于VOF数值模拟方法,对含气泡油膜沿倾斜壁面的流动行为进行研究,考察了气泡的存在对油膜形态和流动速度的影响规律,以及气泡破裂阶段空腔邻域内流体压力变化特性。研究表明,油膜夹带气泡的形变和迁移诱发气泡周围微流场的速度扰动现象,导致气液界面处产生非均匀速度梯度分布,进而引发油膜表面的形态波动。气泡发生破裂时,油膜空穴部位发生明显的正负压力波动现象,气泡附近壁面将承受一定的交变载荷作用。     

Numerical simulation of a single rising bubble by VOF with surface compression

The capability of the direct volume of fluid method for describing the surface dynamics of a free two‐dimensional rising bubble is evaluated using quantities of a recently published benchmark. The model equations are implemented in the open source computational fluid dynamics library OpenFOAM®. Here, a main ingredient of the numerical method is the so‐called surface compression that corrects the fluxes near the interface between two phases. The application of this method with respect to two test cases of a benchmark is considered in the main part. The test cases differ in physical properties, thus in different surface tension effects. The quantities centre of mass position, circularity and rise velocity are tracked over time and compared with the ones given in the benchmark. For test case one, where surface tension effects are more pre‐eminent, deviations from the benchmark results become more obvious. However, the flow features are still within reasonable range. Nevertheless, for test case two, which has higher density and viscosity ratios and above all a lower influence of the surface tension force, good agreement compared with the benchmark reference results is achieved. This paper demonstrates the good capabilities of the direct volume of fluid method with surface compression with regard to the preservation of sharp interfaces, boundedness, mass conservation and low computational time. Some limitation regarding the occurrence of parasitic currents, bad pressure jump prediction and bad grid convergence have been observed. With these restrictions in mind, the method is suitable for the simulation of similar two‐phase flow configurations. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical investigation of electrohydrodynamic (EHD) effects on bubble deformation under pseudo-nucleate boiling conditions

In this article, the electrohydrodynamic (EHD) effects on nucleate boiling are studied by developing a numerical modelling of EHD effect on bubble deformation in pseudo-nucleate boiling conditions. The volume of fluid (VOF) method is employed to track the interface between the gas–liquid two phases; the user-defined code is written and added to the commercial software FLUENT to solve the electric field and the corresponding electric body force. On this basis, the model is applied to study the EHD effects on heat transfer and fluid flows. An initial air bubble surrounded by liquid CCl₄ and attached to a horizontal superheated wall under the action of electric field is studied. The results of the EHD effect on bubble shape evolution are compared with those of available experiments showing good agreement. The mechanism of EHD enhancement of heat transfer and the EHD induced phenomena including bubble elongation and detachment are analyzed in detail.     

Mass and heat transfer by natural convection in a vertical cavity

This paper reports a fundamental study of laminar natural convection in a rectangular enclosure with heat and mass transfer from the side, when the bouyancy effect is due to density variations caused by either temperature or concentration variations. In the first part of the study scale analysis is used to determine the scales of the flow, temperature and concentration fields in boundary layer flow for all values of Prandtl and Lewis numbers. In particular, scale analysis shows that in the extreme case where the flow is driven by bouyancy due to temperature variations, the ratio of mass transfer rate divided by heat transfer rate scales as $\text{Le}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ only if (Pr > 1, Le < 1) or (Pr < 1, Sc < 1), and as $\text{Le}{}^{\mathrm{<mtext>1</mtext><mtext>3</mtext>}}$ if (Pr > 1, Le > 1) or (Pr < 1, Sc > 1). In the second part of the study, the boundary layer scales derived in the first part are used to determine the heat and mass transport characteristics of a vertical slot filled with fluid. Criteria for the existence of distinct thermal and concentration boundary layers in the slot are determined. Numerical solutions for the flow and concentration fields in a slot without distinct thermal boundary layers are reported. These solutions support further the method of scale analysis employed in the first part of the study     

Surface condition effects on flame impingement heat transfer

Flame impingement heating is used in many industrial applications, including the heating and melting of both glass and metal. This heating process usually comprises multiple heat transfer mechanisms, such as forced convection, thermal radiation, and thermochemical heat release. However, little experimental data are available that can be used to determine the importance of each mechanism. This information would be useful for optimizing the heating process and for developing computer models. The objective of this study was to determine the relative importance of thermal radiation and to determine how the thermochemical heat release is affected by the surface properties of the target.This study investigated the heat transfer from oxygen-enhanced, natural gas flames (15 kW) impinging normal to a water-cooled metal disk (d_b = 135 mm) segmented into concentric calorimetric rings. Polished, untreated, and blackened surfaces were used to study emissivity effects. The heat flux to the blackened and polished surfaces was the highest and lowest, respectively. The flux to untreated surfaces was between the highest and lowest fluxes. The largest difference in the flux, between the polished and blackened surfaces, was only 9.8%. Catalyticity effects were investigated by using alumina-coated (nearly noncatalytic), untreated, and platinum-coated (highly catalytic) surfaces. The heat flux to platinum-coated surfaces was the highest. The fluxes to untreated surfaces were similar to those for alumina-coated surfaces. The largest difference in the flux, between the platinum-coated and the alumina-coated surfaces, was only 12%. Therefore, both nonluminous flame radiation and the thermochemical heat release from surface catalytic reactions were relatively small fractions of the total heat flux.     

Free nonlinear oscillations of a thermally relaxing spherical gas bubble in an incompressible liquid

The quasi-adiabatic regime of free oscillation of a bubble in the presence of irreversible interphase heat transfer between the bubble and the ambient liquid is studied. On the basis of simplified model equations of a rarefield bubble mixture, a nonlinear-oscillation equation of the relaxation type is obtained. In constructing an exact particular solution of this equation, the heat transfer law associated with bubble compression is established. For studying the harmonic oscillations, the Krylov-Bogolyubov-Mitropol’skii asymptotic method is used. It is shown that, for a small bubble, the viscosity and heat transfer effects are of the same order. For a small bubble, the influence of these effects on the formation of the natural-oscillation frequency, which is small in the linear approximation, may be significant in the nonlinear formulation. For a large bubble, the influence of these effects is negligible in both approximations. For the approximate solution of the nonlinear equation, a uniformly valid second-order expansion is constructed.     

Shear-opposed mixed-convection flow and heat transfer in a narrow, vertical cavity

A combined experimental and theoretical study has been carried out to investigate mixed-convection heat transfer in a narrow, vertical cavity. The shear force is produced at the left side of a cavity by a belt moving upward that constitutes the sixth wall of the cavity. The left wall of the cavity was cooled and the opposite (right wall) was heated. Hence, the buoyancy force tries to bring the fluid down, and the shear force tends to induce upward fluid flow. The test cell was equipped with two heat exchangers and three thermocouple racks for measuring the temperature distributions at 12 different positions. The temperature field was scanned in the cavity for various flow and temperature boundary conditions. Three-and-two-dimensional laminar models were used to analyze the problem theoretically. The experimental measurements were found to be in good agreement with the numerical predictions.     

Simultaneous effects of induced magnetic field and heat and mass transfer on the peristaltic motion of second‐order fluid in a channel

In this article, we investigate the influence of heat and mass transfer on the peristaltic flow of magnetohydrodynamic second‐order fluid in a channel when the induced magnetic field effects are present. Problem formulation in a wave frame of reference is presented. The governing nonlinear analysis is carried out under the assumption of small wave number. Explicit expressions of the pressure gradient, the stream function, the magnetic force function, the axial induced magnetic field, the current density distribution, the temperature, and the concentration distribution are derived. The effects of embedded parameters are also examined. Copyright © 2011 John Wiley & Sons, Ltd.     

Direct numerical simulation of turbulent heat transfer on the Reynolds analogy over irregular rough surfaces

The effect of rough surface topography on heat and momentum transfer is studied by direct numerical simulations of turbulent heat transfer over uniformly heated three-dimensional irregular rough surfaces, where the effective slope and skewness values are systematically varied while maintaining a fixed root-mean-square roughness. The friction Reynolds number is fixed at 450, and the temperature is treated as a passive scalar with a Prandtl number of unity. Both the skin friction coefficient and Stanton number are enhanced by the wall roughness. However, the Reynolds analogy factor for the rough surface is lower than that for the smooth surface. The semi-analytical expression for the Reynolds analogy factor suggests that the Reynolds analogy factor is related to the skin friction coefficient and the difference between the temperature and velocity roughness functions, and the Reynolds analogy factor for the present rough surfaces is found to be predicted solely based on the equivalent sand-grain roughness. This suggests that the relationship between the Reynolds analogy factor and the equivalent sand-grain roughness is not affected by the effective slope and skewness values. Analysis of the heat and momentum transfer mechanisms based on the spatial- and time-averaged equations suggests that two factors decrease the Reynolds analogy factor. One is the increased effective Prandtl number within the rough surface in which the momentum diffusivity due to the combined effects of turbulence and dispersion is larger than the corresponding thermal diffusivity. The other is the significant increase in the pressure drag force term above the mean roughness height.