A numerical model for cloud cavitation based on bubble cluster

The cavitation cloud of different internal structures results in different collapse pressures owing to the interaction among bubbles. The internal structure of cloud cavitation is required to accurately predict collapse pressure. A cavitation model was developed through dimensional analysis and direct numerical simulation of collapse of bubble cluster. Bubble number density was included in proposed model to characterize the internal structure of bubble cloud. Implemented on flows over a projectile, the proposed model predicts a higher collapse pressure compared with Singhal model. Results indicate that the collapse pressure of detached cavitation cloud is affected by bubble number density.     

Observations on the dynamics of bubble cluster in an ultrasonic field

The dynamics driven interaction between the bubbles in a cavitation cluster is known to be a complex phenomenon indicative of a highly active nonlinear as well as chaotic behavior in ultrasonic fields. By considering the cluster of encapsulated microbubble with a thin elastic shell in ultrasonic fields, in this paper, the dynamics of microbubbles has been studied via applying the methods of chaos physics. Bifurcation, Lyapunov exponent, and time series are plotted with respect to variables such as amplitude, initial bubble radius, frequency and viscosity. The findings of the study indicate that a bubble cluster undergoes a chaotic unstable region as the amplitude and frequency of ultrasonic pulse are increased mainly due to the period doubling phenomenon. The results of the present study are supported by findings of previous studies.     

Elongated bubble shape in inclined air-water slug flow

The shape of elongated bubbles in upward inclined air-water slug flow was studied experimentally by quantitative measurements of the cross sectional distribution of the phases within the pipe, using a wire mesh sensor. Ensemble-averaged shapes of elongated bubbles were determined for a wide range of gas and liquid flow rates, as well as for different pipe inclination angles. The elongated bubble nose can be characterized by an annular domain where liquid is present above the gas. The effect of gas and liquid flow rates, as well as of the pipe inclination angle on the bubble shape (front and tail) is studied. A simplified theoretical model is proposed to determine the bubble front shape. The model predictions compare favorably with the experimental results.     

Interactions of multiple spark-generated bubbles with phase differences

This paper aims to study the complex interaction between multiple bubbles, and to provide a summary and physical explanation of the phenomena observed during the interaction of two bubbles. High-speed photography is utilized to observe the experiments involving multiple spark-generated bubbles. Numerical simulations corresponding to the experiments are performed using the Boundary Element Method (BEM). The bubbles are typically between 3 and 5 mm in radius and are generated either in-phase (at the same time) or with phase differences. Complex phenomena are observed such as bubble splitting, and high-speed jetting inside a bubble caused by another collapsing bubble nearby (termed the ‘catapult’ effect). The two-bubble interactions are broadly classified in a graph according to two parameters: the relative inter-bubble distance and the phase difference (a new parameter introduced). The BEM simulations provide insight into the physics, such as bubble shape changes in detail, and jet velocities. Also presented in this paper are the experimental results of three bubble interactions. The interesting and complex observations of multiple bubble interaction are important for a better understanding of real life applications in medical ultrasonic treatment and ultrasonic cleaning. Many of the three bubble interactions can be explained by isolating bubble pairs and classifying their interaction according to the graph for the two bubble case. This graph can be a useful tool to predict the behavior of multiple bubble interactions.     

Effect of orifice shape on bubble formation mechanism

The effect of orifice shape on the mechanism of bubble formation in gas–liquid two-phase flow is investigated experimentally with three different orifice geometries regarding a circle, a square, and a triangle with same cross-sectional areas. The liquid and gas phases are purified water at 20 °C and air at room temperature, respectively. Gas is injected at the rate of 50–1200 mlph into a stagnant pool of liquid in distances of 5, 10, and 15 cm below the liquid surface. The position, velocity, and acceleration of bubbles are measured at bubbles’ centers of mass (CM) and the effects of these parameters on the bubble volume are investigated. Moreover, the forces acting on a bubble are balanced and the effects of geometry and gas flow rate on each force are presented. In addition, the changes of the acting forces versus time are plotted and discussed for a specific condition. Results show the bubbles formed with the square and circular orifice cross-sectional areas have the most and least volumes at detachment, respectively.     

基于气泡形态影响的水下气幕对冲击波衰减效果分析

气泡帷幕是水下爆炸冲击波防护的重要手段,对其作用机理及技术参数的深入研究对水下爆破安全与应用具有重要意义。采用高速摄影技术对室内小型水下气泡帷幕模型拍摄发现气幕在形成过程和与水下爆炸冲击波相互作用过程中均具有高度非连续性和非均匀性,且气幕区域内气体与液体混杂,界面轮廓复杂多样。在此基础上,考虑气泡形状及界面影响下,通过LS-DYNA有限元软件自带的APDL语言进行编程,实现了在设定的气幕区域内,通过设定气泡直径变化范围及气泡直径之间的最小差异值随机投放一定数量不同直径的气泡来模拟真实气幕中气泡的分布,并通过改变固定区域内气泡个数来模拟不同气压值工况下的气幕效果。分析发现该方法能够更加真实反映气幕在冲击波防护过程中的防护机理,随着单位区域内气泡数量的增大,防护效果越明显,但当气泡数量达到一定数量后气幕整体连续性及稳定性基本固定,防护效果也趋于稳定。     

Numerical simulation of bubble breakup phenomena in a narrow flow field

Based on the boundary integral method, a 3D bubble breakup model in a narrow flow field is established, and a corresponding computation program is developed to simulate the symmetrical and asymmetrical bubble breakup. The calculated results are compared with the experimental results and agree with them very well, indicating that the numerical model is valid. Based on the basic behavior of bubbles in a narrow flow field, the symmetrical and asymmetrical bubble breakup is studied systematically using the developed program. A feasibility rule of 3D bubble breakup is presented. The dynamics of sub-bubbles after splitting is studied. The influences of characteristic parameters on bubble breakup and sub-bubble dynamics are analyzed.     

Physical and transport properties governing bubble column operations

In this work, physical and transport properties governing fluid dynamics in bubble columns were studied in the heterogeneous liquid circulation regime of operation based on several published experimental data. It was confirmed that the axial dispersion phenomenon of the liquid phase could be defined by the constant Peclet number predicated on the two-phase slip velocity, which comprised the time-averaged bubbly flow velocity superimposed by the velocity of large bubbles rising in the fluid. In addition, the slip velocity of mean bubbles in the bubbly flow not only indicated a constant rising velocity of a particular liquid but also apparently reduced to that for the Higbie model with respect to the bubble diameter and its terminal rising velocity, when applied to Akita and Yoshida's dimensionless correlation for mass transfer operation (Akita and Yoshida, 1974). In the course of demonstration, it was confirmed that a semi-theoretical correlation for the gas holdup proposed by Mersmann (1978) is similar to the empirical correlation proposed by Akita and Yoshida (1973) with respect to the expression and the physical perspective.     

Marangoni效应下气泡上升过程的FTM模拟

采用界面跟踪法FTM(Front-Tracking Method),并结合CSF(continuum surface force)模型,研究了在垂直方向上温度分布不均匀的对称流场中由Marangoni效应引起的气泡上升运动问题。模拟了在不同的M a数和Pr数下单气泡及同轴双气泡的运动。研究结果表明,在不同的M a数下气泡的运动差异较大,M a数越大,气泡运动至稳态时的速度越大,且气泡运动的最大速度值与M a数呈正相关关系;增大Pr数所造成的粘度增大或热扩散率减小将削弱气泡的迁移运动;Marangoni对流中双气泡的局部运动证实了温度梯度和气泡运动速度紧密相关。     

Strong interaction between a buoyancy bubble and a free surface

The growth and collapse of buoyant vapor bubbles close to a free surface in an inviscid incompressible fluid is investigated in this paper. The strong interaction between the deforming bubble and the free surface is simulated numerically by a boundary-integral method (Taib 1985; Blake et al., 1987). Improvements are made in the calculation of the singular integrals, the use of nonuniform boundary elements, and the choice of time-step size. The present numerical results agree better with the experimental observations of Blake and Gibson (1981) than previous numerical predictions for bubbles initiated at one maximum radius from the free surface. There is also concurrence of flow features with the experiments for a bubble initiated as close as half maximum radius from the free surface, where other numerical efforts have failed. The effects of buoyancy on bubbles initiated close to a free surface are also investigated. Vastly different features, depending on the distance of the bubble to the free surface and the buoyancy-force parameter, have been observed.     

Motion of a gas bubble in fluid under vibration

This paper is concerned with the analysis of motion of a gas bubble in a uniformly oscillating incompressible fluid. A theoretical model explaining the effect of sinking of gas bubbles in the absence of a standing pressure wave is validated experimentally. The conditions under which this effect occurs are determined, and a simple formula is derived for the average velocity of a gas bubble in the fluid.     

近自由面的多个水下爆炸气泡相互作用研究

将气泡运动阶段周围的流场假设为无黏、无旋、不可压缩的理想流体,运用边界积分法模拟流场中气泡的运动,并开发了三维计算程序,计算值与实验值吻合较好.用该方法模拟了近自由面多气泡之间的相互作用,包括同相气泡和异相气泡. 通过计算发现,气泡的周期随两气泡中心的距离减小而增大,这是由于多气泡之间存在抑制作用,特别是对异相气泡,这种抑制作用更加明显,称之为多气泡之间的抑制效应. 无论有、无自由面存在,多气泡之间均存在抑制效应,由于抑制效应导致同相与异相气泡相互耦合作用的动态特性存在巨大的差异,这些现象可为将来研究多个同时或延时产生的水下爆炸气泡的威力提供参考.      

Robust numerical analysis of the dynamic bubble formation process in a viscous liquid

The dynamics of bubble formation from a submerged nozzle in a highly viscous liquid with relatively fast inflow gas velocity is studied numerically. The numerical simulations are carried out using a sharp interface coupled level set/volume-of-fluid (CLSVOF) method and the governing equations are solved through a hydrodynamic scheme with formal second-order accuracy. Numerical results agree well with experimental results and it is shown that the sharp interface CLSVOF method enables one to reproduce the bubble formation process for a wide range of inflow gas velocities. From numerical results, one can improve their understanding of the mechanisms regarding the dynamics of bubble formation. For example, it is found that for some sets of parameters that the bubble formation process reaches steady state after several bubbles are released from the nozzle. At steady state, bubbles uniformly rise freely in the viscous liquid. It is observed that the fluid flow around a formed bubble has a significant role in determining the overall dynamic process of bubble formation; e.g. the effect of the fluid flow from the preceding bubble can be seen on newly formed bubbles.     

On the modeling and simulation of a laser‐induced cavitation bubble

Single cavitation bubbles exhibit severe modeling and simulation difficulties. This is due to the small scales of time and space as well as due to the involvement of different phenomena in the dynamics of the bubble. For example, the compressibility, phase transition, and the existence of a noncondensable gas inside the bubble have strong effects on the dynamics of the bubble. Moreover, the collapse of the bubble involves the occurrence of critical conditions for the pressure and temperature. This adds extra difficulties to the choice of equations of state. Even though several models and simulations have been used to study the dynamics of the cavitation bubbles, many details are still not clearly accounted for. Here, we present a numerical investigation for the collapse and rebound of a laser‐induced cavitation bubble in liquid water. The compressibility of the liquid and vapor are involved. In addition, great focus is devoted to study the effects of phase transition and the existence of a noncondensable gas on the dynamics of the collapsing bubble. If the bubble contains vapor only, we use the six‐equation model for two‐phase flows that was modified in our previous work [A. Zein, M. Hantke, and G. Warnecke, J. Comput. Phys., 229(8):2964‐2998, 2010]. This model is an extension to the six‐equation model with a single velocity of Kapila et al. (Phys. Fluid, 13:3002‐3024, 2001) taking into account the heat and mass transfer. To study the effect of a noncondensable gas inside the bubble, we add a third phase to the original model. In this case, the phase transition is considered only at interfaces that separate the liquid and its vapor. The stiffened gas equations of state are used as closure relations. We use our own method to determine the parameters to obtain reasonable equations of state for a wide range of temperatures and make them suitable for the phase transition effects. We compare our results with experimental ones. Also our results confirm some expected physical phenomena. Copyright © 2013 John Wiley & Sons, Ltd.     

Numerical simulation of a single bubble by compressible two‐phase fluids

The present work deals with the numerical investigation of a collapsing bubble in a liquid–gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non‐conservative transport equations of the pressure law coefficients. The original first‐order discretization is extended to higher order applying second‐order ENO reconstruction to the primitive variables. The derivation of the non‐conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale‐based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2‐D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser‐induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd.     

FCC辐照金属铜材料中氦泡平衡内压的分子模拟及其尺度效应

氦泡内的氦密度及其内压对含氦泡辐照材料的力学性能具有重要影响，本文采用分子模拟方法系统地研究了面心立方金属铜材料内氦泡平衡内压及其尺度效应、温度效应和引起的应力场。基于能量最小原理与应力平衡准则，提出了确定平衡内压的方法，研究发现：1）采用能量最低原理与应力平衡准则所获得的平衡内压是自洽的；2）存在一个临界尺度，当氦泡孔径小于该临界尺度时，平衡内压出现反常尺度效应，即当氦泡孔径小于3nm时，平衡内压并不随孔径减小而明显增大，甚至出现减小的异常现象；3）传统采用基于球形孔洞的理论公式估算严重高估了氦泡平衡内压，引起的误差随氦泡孔径减小而显著增大，如孔径为3nm时，误差超过63%；4）由于纳米孔洞多面体特性和材料的各向异性，即使在平衡内压条件下，基体铜材料内仍存在一定的应力分布，该应力体现出显著的局域特征，即随着距孔洞中心距离的增加而快速下降。此外，本文还提出采用时间平均叠加区域平均来减小热涨落对应力场的影响，取得了较好的效果。     

FCC辐照金属铜材料中氦泡平衡内压的分子模拟及其尺度效应

氦泡内的氦密度及其内压对含氦泡辐照材料的力学性能具有重要影响,本文采用分子模拟方法系统地研究了面心立方金属铜材料内氦泡平衡内压及其尺度效应、温度效应和引起的应力场。基于能量最小原理与应力平衡准则,提出了确定平衡内压的方法,研究发现：1）采用能量最低原理与应力平衡准则所获得的平衡内压是自洽的;2）存在一个临界尺度,当氦泡孔径小于该临界尺度时,平衡内压出现反常尺度效应,即当氦泡孔径小于3nm时,平衡内压并不随孔径减小而明显增大,甚至出现减小的异常现象;3）传统采用基于球形孔洞的理论公式估算严重高估了氦泡平衡内压,引起的误差随氦泡孔径减小而显著增大,如孔径为3nm时,误差超过63%;4）由于纳米孔洞多面体特性和材料的各向异性,即使在平衡内压条件下,基体铜材料内仍存在一定的应力分布,该应力体现出显著的局域特征,即随着距孔洞中心距离的增加而快速下降。此外,本文还提出采用时间平均叠加区域平均来减小热涨落对应力场的影响,取得了较好的效果。     

A bench top experimental model of bubble transport in multiple arteriole bifurcations

Motivated by a novel gas embolotherapy technique, a bench top vascular bifurcation model is used to investigate the splitting of long bubbles in a series of liquid-filled bifurcations. The developmental gas embolotherapy technique aims to treat cancer by infarcting tumors with gas emboli that are formed by selective acoustic vaporization of 6 μm, intravascular, perfluorcarbon droplets. The resulting gas bubbles are large enough to extend through several vessel bifurcations. The current bench top experiments examine the effects of gravity and flow on bubble transport through multiple bifurcations. The effect of gravity is varied by changing the roll angle of the bifurcating network about its parent tube. Splitting at each bifurcation is nearly even when the roll angle is zero. It is demonstrated that bubbles can either stick at one of the second bifurcations or in the second generation daughter tubes, even though the flow rate in the parent tube is constant. The findings of this work indicate that both gravity and flow are important in determining the bubble transport, and suggest that a treatment strategy that includes multiple doses may be effective in delivering emboli to vessels not occluded by the initial dose.     

Quasi-static bubble in a yield stress fluid: elasto-plastic model

The influence of the yield stress of Carbopol® gel dispersions on the behaviour of quasi-static bubbles was investigated. Many fluids, from many different industrial fields, have yield stress behaviour. Most of them contain gas bubbles. To study bubble behaviour in such suspensions, a transparent model fluid (dispersion of Carbopol® in water) was used. The experimental device allowed to quasi-statically increase bubble internal pressure with small pressure step to reach a maximum target internal pressure and the pressure setpoint was inverted to return to the initial pressure. Hysterical behaviour of the bubbles was highlighted as they did not regain their initial shape because of yield stress. We show that the rheological behaviour is related to the internal pressure, bubble geometry and yield stress in quasi-static conditions. A modification of the Laplace law depending on the yield stress of the fluid and bubble sphericity was proposed and validated.