竖直环形通道内液氮流动沸腾的MUSIG模型分析 第二部分:径向气泡直径和界面面积浓度的预测

由于考虑了气泡的破裂和聚合,同两流体模型相比,MUSIG模型(多尺寸组模型)能更准确地描述流场内气泡直径。采用MUSIG模型详细分析了不同壁面热流量,液体入口速度,过冷度以及不同管道高度时通道内气泡相界面面积、当地气泡直径、空泡系数等参数沿径向的分布。分析结果表明,MUSIG模型可用来预测泡弹状流型转变区的流动参数,也即该模型拓展了两流体模型的使用范围。     

竖直环形通道内液氮流动沸腾的MUSIG模型分析 Part 1:轴向气泡直径和界面面积浓度的预测

近年来群体平衡模型受到了极大的关注。MUSIG模型(多尺寸组模型),为群体平衡模型和双流体模型的结合提供了一个框架。文中介绍了MUSIG模型,并采用MUSIG模型分析了环形通道中的液氮流动沸腾。详细分析了不同壁面热流量、液体入口速度和过冷度时通道内气泡相界面面积、当地气泡直径、空泡系数等参数沿轴向的分布。     

数值模拟两相汽蚀流动的新模型和算法

提出了数值模拟两相汽蚀流动现象的新汽蚀模型和算法.提出的汽蚀模型和算法耦合了考虑紊流粘性效应的Reynolds-Averaged Navier-Stokes方程求解方法,可以自动模拟空泡起始点、空泡长度和汽蚀空泡形状.在流场计算和界面修正的迭代计算过程中,跟踪并得到液相/气相界面.数值模拟了圆锥形圆柱体和半球形圆柱体在不同汽蚀系数下的汽蚀流动现象,得到了与实验值完全吻合的数值模拟结果,并且与已发表的数值研究结果进行了比较.计算结果表明提出的汽蚀模型和算法能够有效地模拟汽蚀流动中的气泡界面和气泡长度,汽蚀模型和算法的正确性和实用性得到了相应的验证.     

方截面鼓泡床气液两相瞬态数值研究

在双流体模型中引入界面浓度输运方程,利用界面浓度和气泡平均Sauter直径模化各相间作用力。引入一附加 湍动能输运方程模化气泡诱导引起的液相湍流。利用该模型对方截面鼓泡床内气液两相流进行三维瞬态数值模拟。计算结 果表明该模型能较好得模拟方截面鼓泡床内气液两相流时均和瞬态流动特征。     

数值模拟颗粒旋转对流化床内气固两相流动特性的影响

流化床内颗粒自旋转将影响颗粒相的流动特性.本文运用基于颗粒动理学理论的欧拉-欧拉气固多相流模型,考虑颗粒自旋转流动对颗粒碰撞能量交换和耗散的影响,数值模拟流化床内气体颗粒两相流动特性.计算结果表明颗粒的自旋转使得床内更容易形成气泡,颗粒浓度分布变化增大.颗粒自旋转运动将导致床内非均匀结构更明显.     

竖直矩形通道内液体流动

通过对沸腾气泡在液体中的受力分析,建立了沸腾气泡长大过程的动力学方程;进而获得了沸腾气泡的生长速率与脱离直径的计算方法.采用图象捕集与处理系统,对竖直矩形通道内液体流动沸腾气泡长大与脱离行为进行实验测定,结合模型求解,获得了气泡生长速率、气泡脱离直径、气泡与加热壁面的接触角等参数随操作条件的变化;由模型计算所得的气泡脱离直径与实测值较为符合.     

双大涡模拟方法模拟提升管内气固两相流动

基于Smagorinsky涡黏模型以及颗粒动理学理论,建立了气固两相流双大涡模拟模型。考虑大涡模拟中过滤尺度的影响,给出颗粒相亚格子压力和热传导系数计算模型。考虑颗粒聚团对两相作用的影响,给出了考虑颗粒聚团作用的气固两相多尺度曳力系数模型。数值模拟了提升管内气固两相流动特性,合理地预测出了提升管内气固两相环-核流动结构。模拟结果与Knowlton等实测结果相吻合。     

高温堆球流模拟在堆物理中的应用

介绍了高温气冷球床反应堆物理计算中燃料元件流动特性模拟的方法,对10MW高温堆进行了计算,与未考虑堆芯中燃料曲线流动的简化计算结果进行了比较。     

声致发光气泡内水蒸气的影响

提出了一个单气泡声致发光的简单计算模型.这个模型是在均匀压强近似下，考虑质量和温度在气泡内的非均匀分布，同时考虑了水蒸气在气泡壁上的凝结与蒸发以及水蒸气在气泡内相对惰性气体的质量扩散.通过Saha方程估算气体电离密度，利用电子与离子、电子与中性粒子的轫致辐射，电子与离子的复合辐射公式估算气泡的辐射强度.不考虑化学反应，计算了不同水温时的气泡发光强度，发现当水温在0 ℃时轫致辐射发光模型比较符合实验结果，水温升高时，如水温为20 ℃或以上，轫致辐射发光模型的计算与实验结果出现数量级差别.考虑化学反应，轫致辐射发光模型的计算则总是比实验结果低2个数量级.     

油滴撞击油膜层内气泡的变形与破裂过程的数值模拟

旋转工作的机械零部件和机械设备的润滑系统工作过程中普遍存在着油滴和油膜的碰撞行为,这一行为易引起气泡夹带现象.气泡将对油滴撞击油膜时的运动过程和附壁油膜层的形成质量造成不可忽视的影响.基于耦合的水平集-体积分数方法,对油滴撞击含气泡油膜的行为进行数值模拟研究,考察油膜层内气泡的变形运动过程,分析气泡大小和位置等因素对撞击过程中气泡变形特征参数的影响规律,并探讨气泡破裂的动力学机制.研究表明,随着气泡直径的增大,油滴撞击含气泡油膜后气泡会依次出现自由表面破裂、稳定变形以及油膜内部破裂等现象;直径d=20μm的气泡能较稳定地存在于油膜层内,同时该值也是气泡发生自由表面破裂和油膜内部破裂的临界值.此外,气泡所在位置同样对气泡变形历程有一定影响,气泡越接近油膜表面,其变形量越大;位于油膜底层的气泡会附着在壁面上.在自由表面破裂和油膜内部破裂过程中,气泡破裂是由气-液界面不稳定引起的,表面张力对这两种现象起重要作用;而黏性剪切力对油膜内部破裂现象也有着不可忽视的影响.     

Simulation of turbulent non-isothermal polydisperse bubbly flow behind a sudden tube expansion

The results of numerical simulation of the structure of non-isothermal polydisperse bubbly turbulent flow and heat transfer behind a sudden tube expansion are presented. The study was carried out at a change in the initial diameter of the air bubbles within d _m1 = 1–5 mm and their volumetric void fraction β = 0–10 %. Small bubbles are available in almost the entire cross section of the tube, while the large bubbles pass mainly through the flow core. An increase in the size of dispersed phase causes the growth of turbulence in the liquid phase due to flow turbulization, when there is a separated flow of liquid past the large bubbles. Adding the air bubbles causes a significant reduction in the length of the separation zone and heat transfer enhancement, and these effects increase with increasing bubble size and their gas volumetric flow rate ratio.     

枝晶生长和气泡形成的数值模拟

本文建立了二维的格子玻尔兹曼方法-元胞自动机(lattice Boltzmann method-cellular automaton, LBM-CA)耦合模型, 对凝固过程中枝晶生长和气泡形成进行模拟研究. 本模型采用CA方法模拟枝晶的生长, 根据界面溶质平衡法计算枝晶生长的驱动力. 采用基于Shan-Chen多相流的LBM模拟气泡在液相中的生长和运动. 在LBM-CA的耦合模型中包含了固-液-气三相之间的相互作用. 应用Laplace定理和模拟气-液-固三相之间的润湿现象对模型进行了验证. 应用所建立的LBM-CA耦合模型模拟研究了气-液相互作用系数对单气泡生长的影响. 发现单气泡的生长速度和平衡半径随气-液相互作用系数的增大而增大. 定向凝固过程中枝晶和气泡生长的模拟结果再现了枝晶的择优生长、 气泡的优先形核位置、气泡的长大、合并、在枝晶间受挤变形以及在液相通道中的运动等物理现象, 与实验结果符合良好. 此外, 初始气体含量越高, 凝固结束时气泡的体积分数也相对较高. 本模型的模拟结果可以揭示在凝固过程中气泡形核、 生长和运动演化以及与枝晶生长相互作用的物理机理.     

泡状流运动的湍流模型及液相速度和湍动能的分布预测

１前言在核能、航天等技术领域和能源、动力、石油化工等工业过程都存在泡状流现象,其中绝大多数的泡状流流动状态为湍流,因此对湍流泡状流进行深入研究极为必要。近年来Ｌｅｅ［１］、Ｌｏｐｅｚ　ｄｅ　Ｂｅｒｔｏｄａｎｏ［２,３］等提出了一些多维的湍流输运的双流...     

气泡-液体闭式多股射流的大涡模拟研究

采用气泡-液体两相流动的欧拉-拉氏大涡模拟,研究了矩形通道内多股射流形成的气泡-液体两相湍流流动,得到了气液两相湍流瞬态结构,产生和发展过程。研究结果发现气泡和液体都有瞬态大涡结构,气泡脉动比液体的强。大涡模拟统计结果给出了有无气泡两种情况下的液体湍流脉动速度均方根值分布。瞬态和统计结果都表明,气泡增大了液体湍流,液体湍流来源于其自身的剪切产生和气泡的作用。这与二阶矩模型模拟结果定性一致。     

Constructing ultrasonic images of soft spherical scatterers

The paper considers specific features of ultrasonic visualization of gas bubbles in a liquid or a medium of like soft biological tissue type under conditions when the size of scatterers is comparable to the acoustic wavelength. It was proposed to use styrofoam specimens as the experimental model of stationary gas bubbles. Patterns of ultrasound scattering by a styrofoam sphere in water were obtained experimentally. It was shown that the measurement results agree well with the prediction of the classical theoretical model of scattering of a plane wave by a perfectly soft sphere. Several experiments were performed illustrating the specific features of visualizing millimeter-sized bubbles. A Terason commercial ultrasonic scanner was used; gelatin specimens with embedded styrofoam spheres served as the objects of study. The simulation and experimental results showed that when bubbles with diameters of <1 mm are visualized, it is impossible to measure the diameter of scatterers because bubbles of different diameters are imaged as bright spots of identical diameter, which is equal to the scanner resolution. To eliminate this difficulty, it is recommended to use the results of theoretical simulation performed in this study, which revealed a monotonic increase in the backscattered signal intensity with an increase in bubble radius. An ultrasonic visualization mode is proposed in which the brightness of scattered signals is used to differentiate between bubbles of different size.     

An Euler–Lagrange method considering bubble radial dynamics for modeling sonochemical reactors

Unsteady numerical computations are performed to investigate the flow field, wave propagation and the structure of bubbles in sonochemical reactors. The turbulent flow field is simulated using a two-equation Reynolds-Averaged Navier–Stokes (RANS) model. The distribution of the acoustic pressure is solved based on the Helmholtz equation using a finite volume method (FVM). The radial dynamics of a single bubble are considered by applying the Keller–Miksis equation to consider the compressibility of the liquid to the first order of acoustical Mach number. To investigate the structure of bubbles, a one-way coupling Euler–Lagrange approach is used to simulate the bulk medium and the bubbles as the dispersed phase. Drag, gravity, buoyancy, added mass, volume change and first Bjerknes forces are considered and their orders of magnitude are compared. To verify the implemented numerical algorithms, results for one- and two-dimensional simplified test cases are compared with analytical solutions. The results show good agreement with experimental results for the relationship between the acoustic pressure amplitude and the volume fraction of the bubbles. The two-dimensional axi-symmetric results are in good agreement with experimentally observed structure of bubbles close to sonotrode.     

Numerical simulation of 3D bubbles rising in viscous liquids using a front tracking method

The rise of bubbles in viscous liquids is not only a very common process in many industrial applications, but also an important fundamental problem in fluid physics. An improved numerical algorithm based on the front tracking method, originally proposed by Tryggvason and his co-workers, has been validated against experiments over a wide range of intermediate Reynolds and Bond numbers using an axisymmetric model [J. Hua, J. Lou, Numerical simulation of bubble rising in viscous liquid, J. Comput. Phys. 22 (2007) 769–795]. In the current paper, this numerical algorithm is further extended to simulate 3D bubbles rising in viscous liquids with high Reynolds and Bond numbers and with large density and viscosity ratios representative of the common air–water two-phase flow system. To facilitate the 3D front tracking simulation, mesh adaptation is implemented for both the front mesh on the bubble surface and the background mesh. On the latter mesh, the governing Navier–Stokes equations for incompressible, Newtonian flow are solved in a moving reference frame attached to the rising bubble. Specifically, the equations are solved using a finite volume scheme based on the Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm, and it appears to be robust even for high Reynolds numbers and high density and viscosity ratios. The 3D bubble surface is tracked explicitly using an adaptive, unstructured triangular mesh. The numerical model is integrated with the software package PARAMESH, a block-based adaptive mesh refinement (AMR) tool developed for parallel computing. PARAMESH allows background mesh adaptation as well as the solution of the governing equations in parallel on a supercomputer. Further, Peskin distribution function is applied to interpolate the variable values between the front and the background meshes. Detailed sensitivity analysis about the numerical modeling algorithm has been performed. The current model has also been applied to simulate a number of cases of 3D gas bubbles rising in viscous liquids, e.g. air bubbles rising in water. Simulation results are compared with experimental observations both in aspect of terminal bubble shapes and terminal bubble velocities. In addition, we applied this model to simulate the interaction between two bubbles rising in a liquid, which illustrated the model’s capability in predicting the interaction dynamics of rising bubbles.