Comparison between level set and phase field method for simulating bubble movement behavior under electric field

In this study, based on different numberical simulation methods, the gas-liquid two-phase flow is taken as the research object. By coupling the continuity equation of incompressible fluid, Navier-Stokes equation, electric field equation and other control equations, a multi-field coupling model for rising bubbles in viscous fluids is established, and numerical simulations are carried out. The two-phase popularity of coupled electric field is studied, and the effect of electric field on bubble motion is analyzed.The Level-set and phase field method are used to track the changes of deformation and rupture during the rising of the bubble. The accuracy and validity of the two methods are verified by mass conservation. At the same time, the calculation area is determined for the accuracy of calculation, and the optimal mesh size is calculated by using mesh independence test. Compared with the level set method, the phase field method has a certain improvement in the calculation efficiency and accuracy. Among them, the calculation efficiency of the phase field calculation method in the same grid is increased by 5 times, and by 3 times in the vertical electric field environment. Moreover, using the phase field method is easier to capture the bubbles slight changes while they are rising, and the quality of the simulation results is better.The simulation analysis of bubble rising process under coupled electric field by two methods shows that under the interaction of electrostatic force, buoyancy and surface tension, the bubble is stretched into an ellipsoid along the direction of the electric field line, and the ratio of the length to the short axis is proportional to the applied electric field strength. In addition, the bubble rising velocity is affected by the electric field, and the vertical electric field accelerates the rising of the bubble.     

Multi-bubble motion behavior of uniform magnetic field based on phase field model

Aiming at the interaction and coalescence of bubbles in gas–liquid two-phase flow, a multi-field coupling model was established to simulate deformation and dynamics of multi-bubble in gas–liquid two-phase flow by coupling magnetic field, phase field, continuity equation, and momentum equation. Using the phase field method to capture the interface of two phases, the geometric deformation and dynamics of a pair of coaxial vertical rising bubbles under the applied uniform magnetic field in the vertical direction were investigated. The correctness of results is verified by mass conservation method and the comparison of the existing results. The results show that the applied uniform magnetic field can effectively shorten the distance between the leading bubble and the trailing bubble, the time of bubbles coalescence, and increase the velocity of bubbles coalescence. Within a certain range, as the intensity of the applied uniform magnetic field increases, the velocity of bubbles coalescence is proportional to the intensity of the magnetic field, and the time of bubbles coalescence is inversely proportional to the intensity of the magnetic field.     

声场与电场作用下空化泡的动力学特性

以液体为工作介质, 利用空化泡的RP控制方程, 模拟分析了无量纲化的电场频率、场强的幅值以及无量纲化的声波频率、 声压幅值的变化对空化泡运动特性的影响. 结果表明: 声场和电场联合作用时, 空化泡运动处于混沌区域范围远高于两者单独作用下空化泡的混沌区域范围. 这不仅对声空化的进一步研究具有重要的理论意义, 而且对于提高和改进空化降解有机污染物的技术也具有指导意义.     

液晶在电场和剪切耦合作用下的流变学行为

本文通过理论和实验对液晶 5CB在剪切和电场耦合作用下流变行为进行了研究. 采用液晶连续理论, 建立了包括界面锚定能, 弹性自由能, 介电自由能和流动能在内的系统 Gibbs自由能公式, 通过最小化系统自由能的方法求解液晶在剪切和电场耦合作用下的取向分布及其黏度变化, 从分子基础模型上揭示了液晶在耦合作用下的流变行为、微观机理及其影响规律, 并通过流变测试对此进行验证. 对比分析了理论和试验结果的误差和原因, 发现界面锚定效应对于液晶分子的取向和黏度具有重要影响. 理论和试验结果均表明, 液晶在电场作用下具有明显的电黏效应, 表现出非牛顿流变行为, 其黏度值由剪切和电场的竞争和耦合作用共同决定. 在外电场作用下液晶的黏度可以增加到初始值的 4倍左右, 液晶这种其自身黏度可随着外场 (例如运动速度) 改变的特性在一定的条件下可以自适应地满足不同工况对黏度的要求, 这对实现智能摩擦润滑具有重要的意义. 关键词： 液晶 流变行为 电黏效应 耦合作用     

Experimental investigation on flow boiling bubble motion under ultrasonic field in vertical minichannel by using bubble tracking algorithm

Bubble dynamics is important in flow boiling of minichannel, and ultrasonic field effects bubble behaviors. However, flow boiling bubble movements in minichannels under ultrasonic field have received little research attention and are still poorly understood. In this paper, the effects of ultrasonic field on bubble dynamics are experimentally studied by capturing the bubble motion behaviors of the flow boiling bubbles. The ultrasonic frequencies are set to 23, 28, 32, and 40 kHz. Bubble tracking algorithm, which studies the growth, trajectories, velocities, and traveled distances for bubbles, is created to qualitatively describe bubble motion behavior of flow boiling in minichannel. It is found that after the application of ultrasound, the detachment frequency, velocity, and travel distance of the bubbles significantly increases, and the growth behavior and trajectory are extremely complex, the two-phase gas-liquid flow is extremely unstable. The bubbles gain kinetic energy as the ultrasound frequency increases. Finally, numerical simulations are used to quantitatively investigate the mechanism of bubble motion in microchannels under ultrasonic fields.     

相同尺度下气泡与复杂壁面的耦合特性研究

采用格子Boltzmann方法(LBM)建立了气液固三相耦合的动力学模型,研究了相同尺度下上浮气泡与复杂壁面的相互耦合作用.首先,基于黏性流体理论,通过构建一组格子Boltzmann(LB)方程来描述气液两相的运动,并以LB离散体积力的形式计入了黏性力、表面张力和重力.同时,采用LBM中的Half-way反弹模型与有限差分格式相结合的方式进行固壁边界的处理.然后,利用本文建立的模型,对不同特征尺寸比条件下,气泡与考虑边缘效应的平面固壁和曲面固壁的耦合特性进行了研究.研究发现固壁边界条件以及特征尺寸比对气泡的运动和拓扑结构的变化都具有明显的非线性影响.最后,研究了流体属性对气泡与复杂壁面耦合规律的影响.     

静止水中单个上升气泡的直接数值模拟

本文发展了基于Front Tracking的直接数值模拟方法研究气液两相界面的迁移特性,该方法对气液两相采用半隐式的分步法直接求解N-S方程,耦合Front Tracking Method获得两相界面的三维变形。针对无边界以及垂直壁面附近静止水中的单个气泡上升过程进行模拟,研究气泡运动的机理以及气泡与壁面的相互作用。数值模拟准确再现了气泡的上升过程和变形,不同Re数下气泡的上升速度计算结果同经验关联式非常吻合,验证了该方法的有效性。随后分析了气泡周围流场的结构,发现壁面对气泡周围流场的抑制是壁面对气泡作用力的主要原因,将导致气泡逐渐偏离垂直壁面。     

II. Recirculation flow induced by a bubble stream rising in a viscous liquid

The recirculation flow induced by the rising motion of a bubble stream in a viscous fluid within an open-top rectangular enclosure is studied. The three-dimensional volume averaged conservation equations are solved by a control-volume method using a hybrid finite differencing scheme to describe the liquid phase hydrodynamics. The momentum exhange between the bubbles and the liquid phase is modeled with a source term equals to the volumetric buoyancy force acting on the gas in the bubble stream. The volumetric buoyancy force accounts for in line interactions between bubbles through the average gas volume fraction in the gas liquid column which depends on the size and the rising velocity of bubbles. The fluid flow within an open-top rectangular enclosure is further investigated by particle image velocimetry for a bubble stream rising in a water-glycerol solution. The measured fluid velocities in a vertical plane are compared with the predictions of the numerical model over a wide range of fluid viscosity (43 mPa s-800 mPa s) and gas flow rates. Finally, the recirculation flows resulting from the interaction of two neighbouring vertical bubble streams are studied. Received: 23 July 1997 / Revised: 19 December 1997 / Accepted: 11 May 1998     

强电场下气泡形变对液体绝缘的影响

 基于液体气泡击穿的椭球模型,推导了气泡形变的流体动力学方程;利用软件Comsol模拟了气泡受力后的形变;根据模拟结果,结合气体击穿的帕邢定律,讨论了气泡形变对液体绝缘的影响。结果表明:气泡在静电引力和表面张力的作用下,沿电场线方向拉伸成椭球,椭球长短轴之比与外加电场强度和气泡初始半径成正比;形变导致电场方向气体通道延长,气泡更容易击穿,外加气压是避免气泡击穿的有效途径之一,外加气压的大小与电场强度和气泡初始半径成正比。     

On the parametric excitation of electrothermal instability in a dielectric liquid layer using an alternating electric field

The onset of electrothermal convective instability of a liquid dielectric subjected to an unsteady electric field is studied in the EHD approximation, when charge formation is produced only due to dielectrophoresis. Convective thresholds are found in two different cases: (i) instability of the liquid equilibrium in a horizontal layer, and (ii) instability of the liquid flow in a vertical layer. The stability boundaries are obtained when there is interaction of dielectrophoretic and gravitational forces. Stability plots of electrical Rayleigh number versus thermal Rayleigh number are given. We show that only synchronous response to variations of the external electric field of finite frequency exists when heating a horizontal layer from above. Quasiperiodic response to the external alternating action is possible in the case of a vertical layer. The influence of the Prandtl number on the stability thresholds is also examined. The asymptotic behavior of the critical parameters in the limiting case of low-frequency modulation is studied using the Wentzel–Kramers–Brillouin method.     

气泡在内置交错矩形肋管道中自由上升的FTM直接数值模拟

采用直接数值模拟的FTM（Front-Tracking Method）方法研究在重力作用下单气泡在竖直内置矩形肋管内的自由上升.选择矩形肋的肋高、肋距与方管宽度的比值作为管道几何特征值,用莫顿数作为流体特征参数,分析不同肋高、肋距及莫顿数情况下的气泡自由上升,研究不同几何特征及流体参数对气泡运动轨迹的影响.研究表明：在矩形肋片的影响下气泡上升时形状不对称,存在水平方向位移,表现为"蛇形"上升.这种现象与管道几何结构及气泡周围流体性质有关,肋高越大,现象越明显;气泡周围流体的粘度越小,水平方向位移越大.     

匀强磁场对水中气泡运动的影响

基于Rayleigh-Plesset方程, 考虑极性水分子在均匀磁场运动受到磁场力作用, 根据能量守恒建立了外磁场作用下单气泡运动的控制方程, 并对附加压强的大小、性质及对气泡运动的影响进行了计算和分析. 结果表明: 随磁场强度的增强, 附加压强线性增大, 气泡膨胀率降低, 最大半径减小, 气泡坍缩速度下降; 外加磁场引起的气泡振动变化规律与增大静态压具有相似的效果.     

电场作用下的气泡受力分析

本文在蠕动流近似的基础上分析计算了具有不同介电常数工质的气泡在外电场力作用下的受力状况。计算结果表明电场力的水平分力使气泡沿场强方向伸长,并且液体和气体相对介电常数之比越大的工质所受的电场力越大。从而,沸腾换热过程中外加电场对液体和气体相对介电常数比值较大的工质能够产生更大的影响。     

水下火箭水平射流初期特征研究

水中工作固体火箭发动机处于重浮力同时作用环境下, 与工作于大气环境下的固体火箭发动机具有不同的工作特性. 为进一步掌握水下固体火箭发动机的工作特性, 对具有重浮力特征的水下射流进行研究, 重点分析重浮力作用下水平喷射射流结构及推力振荡情况, 采用VOF模型对水平喷射且具有重浮力特征的三维发动机模型进行仿真模拟, 对比有/无重浮力下射流气泡的上浮特征, 并采用动量原理对发动机工作初期的射流扰动进行分析, 获得了重浮力下水下固体火箭发动机的推力振荡特征. 研究结果表明: 由于重浮力逐渐占据主导地位, 射流气泡具有明显的上浮特征, 推力与重浮力耦合后在竖直方向产生的翻转力矩更大, 通过与文献中实验对比可见, 采用VOF模型并考虑重浮力后仿真所得射流结构与实验结果更吻合.     

复杂微通道内气泡在浮力作用下上升行为的格子Boltzmann方法模拟

本文采用气-液两相流格子Boltzmann方法模拟了复杂微通道内气泡在浮力作用下的上升过程,主要研究障碍物表面润湿性、浮力大小、障碍物尺寸和气泡初始位置对气泡变形、分裂、合并的动力学行为以及对气泡上升速度、终端速度和气泡剩余质量的运动特性的影响.研究发现,障碍物表面接触角较小时气泡能够完整地通过障碍物通道,随着障碍物表面接触角增加,气泡通过障碍物通道时严重变形,并会发生分裂行为,使得部分气泡黏附在障碍物表面,从而导致气泡到达终端时质量减少.相应地,气泡上升速度以及终端速度也随着微通道表面接触角的增加而减小.另一方面,随着浮力的增加,气泡在上升过程中更容易发生分裂和合并现象,且气泡剩余质量和终端速度随着浮力的增加呈对数形式增加.此外,随着微通道障碍物半径增加,气泡剩余质量首先缓慢减小然后快速减小,而气泡终端速度近似呈线性减小.最后,数值结果还表明当气泡初始位置偏离管道中间时,其上升速度、气泡剩余质量以及气泡终端速度都与初始位置在管道中间时的变化趋势一致,然而对应的数值均减小,且气泡在上升过程中变形更严重.     

过热液体中蒸汽泡上升过程的格子Boltzmann三维数值模拟

应用格子Boltzmann相变模型,在三维空间研究蒸汽泡在过热液体中生长、上升和变形等动力学行为.为研究传热传质对蒸汽泡运动的影响,对比模拟相同条件下气泡在等温环境中上升的物理过程.结果表明：蒸汽泡在过热液体中上升发生的变形程度较小,意味着相变对蒸汽泡的影响和表面张力一样使汽泡保持初始的形状.蒸汽泡在过热液体中的上升速度较小,说明随着汽泡生长拖拽力的影响比浮力大.蒸汽泡生长率在初始阶段达到最大值,随后会趋于一个恒定的值.随着汽泡体积增大和上升速度的增加,其对流场的扰动也越来越剧烈.蒸汽泡生长和上升引起的对流运动对温度场的演化造成很大的影响.     

Flow field around growing and rising vapour bubble by PIV measurement

A study on flow field measurement around growing and rising vapour bubbles by use of PIV technique is presented. Bubbles were generated from single artificial cavities. Experiments have been conducted with saturated boiling of distilled water at atmospheric pressure. In the experiment fluid velocity field surrounding the bubbles was visualized by use of polyamide tracer particles and a sheet of a YAG pulse laser beam. The images were recorded with a cross-correlation CCD-camera. It has been shown that for lower heat flux density bubble growths in an almost quiescent bulk of liquid. For higher heat flux density the train of bubbles creates a vapour column with strong wake effect. Maximum liquid velocity recorded is approximately equal to the terminal velocity of bubble rising in a stagnant liquid.     

金纳米环双体尺寸和耦合效应对表面等离子体共振特性的影响

采用离散偶极子近似方法系统地研究了金纳米环双体的消光光谱及其电场分布. 计算结果表明, 金纳米环双体在耦合作用下的共振消光峰对应着不同振动模式, 改变金纳米环双体的排列方式、 间距和尺寸大小, 其表面等离子体共振消光峰发生红移或蓝移. 因此可以通过对金纳米环双体结构参数和排列方式的设定, 调节其表面等离子体共振消光峰的位置. 电场分布表明, 水平排列的金纳米环双体较单个金纳米环产生更强的局部表面增强电场. 适当的小间距, 较大的内外半径的金纳米环水平阵列更适合做表面增强拉曼散射的衬底, 在生物分子检测等领域具有潜在的应用.     

The Armstrong experiment revisited

When a high-voltage direct-current is applied to two beakers filled with water or polar liquid dielectrica, a horizontal bridge forms between the two beakers. This experiment was first carried out by Lord Armstrong in 1893 and then forgotten until recently. Such bridges are stable by the action of electrohydrodynamic (EHD) forces caused by electric field gradients counteracting gravity. Due to these gradients a permanent pumping of liquid from one beaker into the other is observed. At macroscopic scale several of the properties of a horizontal water bridge can be explained by modern electrohydrodynamics, analyzing the motion of fluids in electric fields. Whereas on the molecular scale water can be described by quantum mechanics, there is a conceptual gap at mesoscopic scale which is bridged by a number of theories including quantum mechanical entanglement and coherent structures in water – theories that we discuss here. Much of the phenomenon is already understood, but even more can still be learned from it, since such “floating” liquid bridges resemble a small high voltage laboratory of their own: The physics of liquids in electric fields of some kV/cm can be studied, even long time experiments like neutron or light scattering are feasible since the bridge is in a steady-state equilibrium and can be kept stable for hours. It is also an electro-chemical reactor where compounds are transported through by the EHD flow, enabling the study of electrochemical reactions under potentials which are otherwise not easily accessible. Last but not least the bridge provides the experimental biologist with the opportunity to expose living organisms such as bacteria to electric fields without killing them, but with a significant influence on their behavior, and possibly, even on their genome.