High speed cine observations of cavitating flow in a duct

The dynamics of cavities produced in cavitating flow confined in a duct was studied. The ultimate purpose of the work is to develop models of the flow to assist in predicting cavitation erosion and noise. Observations of the cavitating flow using high speed cine photography allowed confirmation to be made of the shedding mechanism originally described by Knapp, and measurements of the cavity dimensions to be determined as a function of time. It was found that the time for a cavity to collapse was found three times greater than expected from Rayleigh's classical theory.     

Explosive vaporization in microenclosures

The explosive vaporization of a liquid above planar microheaters induces a fast increase of pressure that is exploited in many thermally driven actuators in MEMS components such as ink jet printer cartridges, pumps, valves and optical switches. Some of these components need to enclose the working fluid as it is the case of valves in which the heated liquid is separated from the flow that it regulates by a flexible membrane. To achieve a better insight into the thermodynamic processes involved, the present work investigates experimentally an enclosed microsystem designed and fabricated for this purpose, composed of a small liquid volume (8 nL) heated by a electric pulse for 2 μs supplied to a planar microfabricated heater. During the heating, the temperature-induced change in resistance can be determined by imposing a defined current and measuring the voltage drop over the heater. While the chip is based on a silicon substrate with integrated platinum heaters and sensors, the structure enclosing the fluid (cavity and fluidic access to it) is made of a silicone elastomer, poly(dimethylsiloxane) (PDMS). This transparent material is widely used in microfluidics and allows for flexible and transparent walls that can be deflected by increasing the pressure inside the cavity. To seal the system the inlet and the outlet were closed by blocking them with a metallic stab. In the present work we visualize vaporization of isopropanol in contact with a suddenly heated planar resistor for two different cavity heights, 150 μm and 16 μm. The rate of temperature rise of the thin liquid layer in contact with the heater is of the order of 10⁷ K s⁻¹ for a pulse duration of 2 μs. We compare bubble growth and collapse for the open and closed systems. Compared to the open system, the bubble growth in the closed system is considerably damped.     

The motion of a small bubble in waves

The motion of a single spherical small bubble due to buoyancy in the ideal fluid with waves is investigated theoretically and experimentally in this article. Assuming that the bubble has no effect on the wave field, equations of a bubble motion are obtained and solved. It is found that the nonlinear effect increases with the increase of the bubble radius and the rising time. The rising time and the motion orbit are given by calculations and experiments. When the radius of a bubble is smaller than 0.5mm and the distance from the free surface is greater than the wave height, the results of the present theory are in close agreement with measurements. The project supported by the National Natural Science Foundation of China     

Cavitation Research and Ship Propeller Design

The role of cavitation research in the design of ship propellers and the influence of research on propeller design is reviewed. The historical development of research on bubble cavitation is an example of a lack of communication between research and design. Research on sheet cavitation is starting now and simplifications such as two dimensional cavitation are being made. It is argued from observations on propellers that the use of two-dimensional cavitaties is not a proper simplification to investigate sheet cavitation. An illustration is also given of the gap between the assessment of the risk of erosion on propeller models and research on erosion. Finally, the simplifications of tip vortex inception and the problems of the inception speed of propellers are discussed.     

缓变主流中三维气泡的非线性振动

空化现象和水下噪声机制与液体中气泡的动力学行为密切相关．在无粘势流的假定下，采用多参数摄动分析，研究了缓变主流中三维气泡的非线性体积模态振动．推导了关于缓变泡形展开的各阶扰动方程，获得了一阶振动的演化方程和一些特殊情况下的解析解；并采用高阶有限元离散的边界积分方程方法，对平面固壁和自由面附近三维气泡的固有频率进行了数值计算     

球状泡群内气泡的耦合振动

振动气泡形成辐射场影响其他气泡的运动, 故多气泡体系中气泡处于耦合振动状态. 本文在气泡群振动模型的基础上, 考虑气泡间耦合振动的影响, 得到了均匀球状泡群内振动气泡的动力学方程, 以此为基础分析了气泡的非线性声响应特征. 气泡间的耦合振动增加了系统对每个气泡的约束, 降低了气泡的自然共振频率, 增强了气泡的非线性声响应. 随着气泡数密度的增加, 振动气泡受到的抑制增强; 增加液体静压力同样可抑制泡群内气泡的振动, 且存在静压力敏感区(1–2 atm, 1 atm=1.01325×10⁵ Pa); 驱动声波对气泡振动影响很大, 随着声波频率的增加, 能够形成空化影响的气泡尺度范围变窄. 在同样的声条件、泡群尺寸以及气泡内外环境下, 初始半径小于5 μm 的气泡具有较强的声响应. 气泡耦合振动会削弱单个气泡的空化影响, 但可延长多气泡系统空化泡崩溃发生的时间间隔和增大作用范围, 整体空化效应增强.     

单液滴撞击倾斜液膜飞溅过程的耦合Level Set-VOF模拟

采用耦合水平集--体积分数法(CLSVOF)对液滴撞击倾斜表面液膜后液膜的形态演化及飞溅过程进行数值模拟, 并对液滴撞击液膜过程中形成的空气卷吸现象进行研究并探讨了撞击角对此的影响, 分析了液滴撞击后液体内部的压力和速度分布, 对液滴撞击倾斜表面液膜的飞溅过程进行讨论, 并与实验结果进行了对比, 验证了CLSVOF方法研究液滴撞击倾斜液膜的可行性. 结果表明, 液滴撞击倾斜液膜时前后两部分飞溅现象产生的机理不同, 前半部分飞溅是由于压差引起的颈部射流, 而后半部分则是由液膜径向流动产生的飞溅现象. 随着撞击角的增大, 空气卷吸气泡数量减少.     

一个新的非常规各向异性元的收敛性分析

构造了一个新的非常规各向异性Hermite型矩形单元并据此对二阶椭圆问题提出了一个混合元格式,同时给出了该格式的收敛性分析.     

The effect of the VOF–CSF static contact angle boundary condition on the dynamics of sliding and bouncing ellipsoidal bubbles

The static contact angle is the only empiricism introduced in a Volume of Fluid–Continuum Surface Force (VOF–CSF) model of bubbly flow. Although it has previously been shown to have a relatively limited effect on the accuracy of velocity and shape predictions in the case of large gas bubbles sliding under inclined walls (e.g. Cook and Behnia, 2001), it may have a more determining influence on the numerical prediction of the dynamics of smaller ellipsoidal bubbles which were shown by Tsao and Koch (1997) to bounce repeatedly when sliding under inclined walls at certain wall inclinations. The present paper reports on the influence of surface tension and the static contact angle on the dynamics of an ellipsoidal air bubble of equivalent diameter D_e = 3.4 mm. The bubble Eötvös and Morton numbers are Eo = 1.56 and Mo = 2 × 10⁻¹¹ respectively. The computational results are achieved with a Piecewise Linear Construction (PLIC) of the interface and are reviewed with reference to experimental measurements of bubble velocity and interface shape oscillations recorded using a high speed digital camera. Tests are performed at plate inclination angles θ ∈ {10°, 20°, 30°, 45°} to the horizontal and computational models consider three static contact angles θ_c ∈ {10°, 20°, 30°}. The static contact angle has been found to have a significant effect on the bubble dynamics but to varying degree depending on the plate inclination. It is shown to promote lift off and bouncing when the plate inclination angle reaches 30° in a way that does not necessarily reflect experimental observations.     

Computational approach for a pair of bubble coalescence process

The coalescence of bubbles has great value in mineral recovery and oil industry. In this paper, two co-axial bubbles rising in a cylinder is modelled to study the coalescence of bubbles for four computational experimental test cases. The Reynolds’ (Re) number is chosen in between 8.50 and 10, Bond number, Bo ∼4.25-50, Morton number, M 0.0125-14.7. The viscosity ratio (μ_r) and density ratio (ρ_r) of liquid to bubble are kept constant (100 and 850 respectively). It was found that the Bo number has significant effect on the coalescence process for constant Re, μ_r and ρ_r. The bubble-bubble distance over time was validated against published experimental data. The results show that VOF approach can be used to model these phenomena accurately. The surface tension was changed to alter the Bo and density of the fluids to alter the Re and M, keeping the μ_r and ρ_r the same. It was found that for lower Bo, the bubble coalesce is slower and the pocket at the lower part of the leading bubble is less concave (towards downward) which is supported by the experimental data.