Modeling cavitation in a rapidly changing pressure field – Application to a small ultrasonic horn

Ultrasonic horn transducers are frequently used in applications of acoustic cavitation in liquids. It has been observed that if the horn tip is sufficiently small and driven at high amplitude, cavitation is very strong, and the tip can be covered entirely by the gas/vapor phase for longer time intervals. A peculiar dynamics of the attached cavity can emerge with expansion and collapse at a self-generated frequency in the subharmonic range, i.e. below the acoustic driving frequency. The term “acoustic supercavitation” was proposed for this type of cavitation Žnidarčič et al. (2014) [1].We tested several established hydrodynamic cavitation models on this problem, but none of them was able to correctly predict the flow features. As a specific characteristic of such acoustic cavitation problems lies in the rapidly changing driving pressures, we present an improved approach to cavitation modeling, which does not neglect the second derivatives in the Rayleigh–Plesset equation. Comparison with measurements of acoustic supercavitation at an ultrasonic horn of 20 kHz frequency revealed a good agreement in terms of cavity dynamics, cavity volume and emitted pressure pulsations.The newly developed cavitation model is particularly suited for simulation of cavitating flow in highly fluctuating driving pressure fields.     

Characterization of cavitation under ultrasonic horn tip – Proposition of an acoustic cavitation parameter

Acoustic cavitation, generated by a piezo-driven transducer, is a commonly used technique in a variety of processes, from homogenization, emulsification, and intensification of chemical reactions to surface cleaning and wastewater treatment. An ultrasonic horn, the most commonly used acoustic cavitation device, creates unique cavitation conditions under the horn tip that depend on various parameters such as the tip diameter, the driving frequency of the horn, its amplitude, and fluid properties. Unlike for hydrodynamic cavitation, the scaling laws for acoustic cavitation are poorly understood. Empirical relationships between cavitation dynamics, ultrasonic horn operating conditions, and fluid properties were found through systematic characterization of cavitation under the tip. Experiments were conducted in distilled water with various sodium chloride salt concentrations under different horn amplitudes, tip geometries, and ambient pressures. Cavitation characteristics were monitored by high-speed (200,000 fps) imaging, and numerous relations were found between operating conditions and cavitation dynamics. The compared results are discussed along with a proposal of a novel acoustic cavitation parameter and its relationship to the size of the cavitation cloud under the horn tip. Similar to the classical hydrodynamic cavitation number, the authors propose for the first time an acoustic cavitation parameter based on experimental results.     

超声空化及其声流结构实验研究*

超声空化及其声流效应在医学、化工和能源等领域得到广泛应用。本文采用高速摄像和粒子图像测速系统分别研究了超声场下的空化形态和声流场结构的时空演化规律。实验研究了50W,100W,200W和250W等四种不同输入功率对18kHz的超声变幅杆附近空化及其声流场的影响。研究结果表明：（1）在变幅杆下端面处观察到由大量空化气泡均匀分布组成的倒置锥形空泡结构,并且锥形空泡结构为稳态流动结构。（2）在超声变幅杆附近产生了两种不同的声流形式,第一种是变幅杆底端的射流型声流,第二种是变幅杆两侧的回旋流。此外,通过研究空泡与声流场中最大速度点之间的空间对应关系,发现声流是因为空泡流动带动而产生的。（3）空间位置和输入功率能显著影响射流型声流的流场结构,但是对回旋流的影响十分微弱。     

Combining COMSOL modeling with acoustic pressure maps to design sono-reactors

Scaled-up and economically viable sonochemical systems are critical for increased use of ultrasound in environmental and chemical processing applications. In this study, computational simulations and acoustic pressure maps were used to design a larger-scale sono-reactor containing a multi-stepped ultrasonic horn. Simulations in COMSOL Multiphysics showed ultrasonic waves emitted from the horn neck and tip, generating multiple regions of high acoustic pressure. The volume of these regions surrounding the horn neck were larger compared with those below the horn tip. The simulated acoustic field was verified by acoustic pressure contour maps generated from hydrophone measurements in a plexiglass box filled with water. These acoustic pressure contour maps revealed an asymmetric and discrete distribution of acoustic pressure due to acoustic cavitation, wave interaction, and water movement by ultrasonic irradiation. The acoustic pressure contour maps were consistent with simulation results in terms of the effective scale of cavitation zones (∼10 cm and <5 cm above and below horn tip, respectively). With the mapped acoustic field and identified cavitation location, a cylindrically-shaped sono-reactor with a conical bottom was designed to evaluate the treatment capacity (∼5 L) for the multi-stepped horn using COMSOL simulations. In this study, verification of simulation results with experiments demonstrates that coupling of COMSOL simulations with hydrophone measurements is a simple, effective and reliable scientific method to evaluate reactor designs of ultrasonic systems.     

PIV quantification of the flow induced by an ultrasonic horn and numerical modeling of the flow and related processing times

The flow in a confined container induced by an ultrasonic horn is measured by Particle Image Velocimetry (PIV). This flow is caused by acoustic streaming and highly influenced by the presence of cavitation. The jet-like experimentally observed flow is compared with the available theoretical solution for a turbulent free round jet. The similarity between both flows enables a simplified numerical model to be made, whilst the phenomenon is very difficult to simulate otherwise. The numerical model requires only two parameters, i.e. the flow momentum and turbulent kinetic energy at the position of the horn tip. The simulated flow is used as a basis for the calculation of the time required for the entire liquid volume to pass through the active cavitation region.     

Designing and characterizing a multi-stepped ultrasonic horn for enhanced sonochemical performance

The commonly used ultrasonic horn generates localized cavitation below its converging tip resulting in a dense bubble cloud near the tip and limiting diffusion of reactive components into the bubble cloud or reactive radicals out of the bubble cloud. To improve contact between reactive components, a novel ultrasonic horn design was developed based on the principles of the dynamic wave equation. The horn, driven at 20 kHz, has a multi-stepped design with a cone-shaped tip increasing the energy-emitting surface areas and creating multiple reactive zones. Through different physical and chemical experiments, performance of the horn was compared to a typical horn driven at 20 kHz. Hydrophone measurements showed high acoustic pressure areas around the horn neck and tip. Sonochemiluminescence experiments verified multiple cavitation zones consistent with hydrophone readings. Calorimetry and dosimetry results demonstrated a higher energy efficiency (31.3%) and a larger hydroxyl radical formation rate constant (0.36 μM min⁻¹) compared to typical horns. In addition, the new horn degraded naphthalene faster than the typical horn tested. The characterization results demonstrate that the multi-stepped horn configuration has the potential to improve the performance of ultrasound as an advanced oxidation technology by increasing the cavitation zone in the solution.     

Modeling acoustic cavitation with inhomogeneous polydisperse bubble population on a large scale

A model for acoustic cavitation flows able to depict large geometries and time scales is proposed. It is based on the Euler–Lagrange approach incorporating a novel Helmholtz solver with a non-linear acoustic attenuation model. The method is able to depict a polydisperse bubble population, which may vary locally. The model is verified and analyzed in a setup with a large sonotrode. Influences of the initial void fraction and the population type are studied. The results show that the velocity is strongly influenced by these parameters. Furthermore, the largest bubbles determine the highest pressure amplitude reached in the domain, which corresponds to the Blake threshold of these bubbles. Additionally, a validation is performed with a small sonotrode. The model reproduces most of the experimentally observed phenomena. In the experiments, neighboring bubbles are found which move in different directions depending on their size. The numerical results show that the responsible mechanism here is the reversal of the primary Bjerknes force at a certain pressure amplitude.     

Understanding the ultrasound field of high viscosity mixtures: Experimental and numerical investigation of a lab scale batch reactor

In this work, mixtures of increasing viscosity (from 0.9 to ≈720 mPas) are sonicated directly using an ultrasonic horn at 30 kHz to investigate the effect of viscosity on the ultrasound field both from an experimental and numerical point of view. The viscosity of the mixtures is modified by preparing water-polyethylene glycol solutions. The impact of the higher viscosity on the acoustic pressure distribution is studied qualitatively and semi-quantitatively using sonochemiluminescence. The velocity of light scattering particles added in the mixtures is also explored to quantify acoustic streaming effects using Particle Image Velocimetry (PIV). A numerical model is developed that is able to predict cavitationally active zones accounting for both thermoviscous and cavitation based attenuation. The results show that two cavitation zones exist: one directly under the horn tip and one around the part of the horn body that is immersed in the liquid. The erosion patterns on aluminum foil confirm the existence of both zones. The intensity of the cavitationally active zones decreases considerably with increasing viscosity of the solutions. A similar reduction trend is observed for the velocity of the particles contained in the jet directly under the tip of the horn. Less erratic flow patterns relate to the high viscosity mixtures tested. Finally, two numerical models were made combining different boundary conditions related to the ultrasonic horn. Only the model that includes the radial horn movements is able to qualitatively predict well the location of the cavitation zones and the decrease of the zones intensity, for the highest viscosities studied. The current findings should be taken into consideration in the design and modelling phase of horn based sonochemical reactors.     

Studies of a novel sensor for assessing the spatial distribution of cavitation activity within ultrasonic cleaning vessels

This paper describes investigations of the spatial distribution of cavitation activity generated within an ultrasonic cleaning vessel, undertaken using a novel cavitation sensor concept. The new sensor monitors high frequency acoustic emissions (>1 MHz) generated by micron-sized bubbles driven into acoustic cavitation by the applied acoustic field. Novel design features of the sensor, including its hollow, cylindrical shape, provide the sensor with spatial resolution, enabling it to associate the megahertz acoustic emissions produced by the cavitating bubbles with specific regions of space within the vessel. The performance of the new sensor has been tested using a 40 kHz ultrasonic cleaner employing four transducers and operating at a nominal electrical power of 140 W under controlled conditions. The results demonstrate the ability of the sensors to identify 'hot-spots' and 'cold-spots' in cavitation activity within the vessel, and show good qualitative agreement with an assessment of the spatial distribution of cavitation determined through erosion monitoring of thin sheets of aluminium foil. The implications of the studies for the development of reliable methods of quantifying the performance of cleaning vessels are discussed in detail.     

Comparison of frequency domain and time domain methods for the numerical simulation of contactless ultrasonic cavitation

The use of a top-mounted electromagnetic induction coil has been demonstrated as a contactless alternative to traditional ultrasonic treatment (UST) techniques that use an immersed mechanical sonotrode for the treatment of metals in the liquid state. This method offers similar benefits to existing UST approaches, including degassing, grain refinement, and dispersion of nanoparticles, while also preventing contact contamination due to erosion of the sonotrode. Contactless treatment potentially extends UST to high temperature or reactive melts. Generally, the method relies on acoustic resonance to reach pressure levels suitable for inertial cavitation and as a result the active cavitation volume tends to lie deep in the melt rather than in the small volume surrounding the immersed sonotrode probe. Consequently, (i) with suitable tuning of the coil supply frequency for resonance, the treatment volume can be made arbitrarily large, (ii) the problem of shielding and pressure wave attenuation suffered by the immersed sonotrode is avoided. However, relying on acoustic resonance presents problems: (i) the emergence of bubbles alters the speed of sound, resonance is momentarily lost, and cavitation becomes intermittent, (ii) as sound waves travel through and reflect on all the materials surrounding the melt, the sound characteristics of the crucible and supporting structures need to be carefully considered. The physics of cavitation coupled with this intermittent behaviour poses a challenge to sonotrode modelling orthodoxy, a problem we are trying to address in this publication. Two alternative approaches will be discussed, one of which is in the time domain and one in the frequency domain, which couple the solution of a bubble dynamics solver with that of an acoustics solver, to give an accurate prediction of the acoustic pressure generated by the induction coil. The time domain solver uses a novel algorithm to improve simulation time, by detecting an imminent bubble collapse and prescribing its subsequent behaviour, rather than directly solving a region that would normally require extremely small time steps. This way, it is shown to predict intermittent cavitation. The frequency domain solver for the first time couples the nonlinear Helmholtz model used for studying cavitation, with a background source term for the contribution of Lorentz forces. It predicts comparable RMS pressures to the time domain solver, but not the intermittent behaviour due to the underlying harmonic assumption. As further validation, the frequency domain method is also used to compare the generated acoustic pressure with that of traditional UST using a mechanical sonotrode.     

液体材料超声处理过程中声场和流场的分布规律研究

针对合金熔体等液体材料的超声处理过程,选取水作为透明模型材料,采用数值模拟计算和示踪粒子实验方法,研究了20和490 kHz两种频率超声作用下水中的声场和流场分布.结果表明,增大变幅杆半径能够提高水中声压水平,扩大空化效应的发生区域.当超声频率为20 kHz时,水中声压最大值出现在超声变幅杆下端面处,且声压沿传播距离的增大而显著减小.如果超声频率增加至490 kHz,水中的声压级相比于20 kHz时明显提高,且声压沿着超声传播方向呈现出周期性振荡特征.两种频率超声作用下水中的流场呈现相似的分布特征,且平均流速均随着变幅杆半径增大表现出先升高后降低的趋势.变幅杆半径相同时,20 kHz频率超声作用下水中的平均流速高于490 kHz频率超声.采用示踪粒子图像测速技术实时观察和测定了水中的流速分布,发现其与计算结果基本一致.     

High-speed imaging of ultrasonic emulsification using a water-gallium system

Aiming at elucidating ultrasonic emulsification mechanisms, the interaction between a single or multiple acoustic cavitation bubbles and gallium droplet interface was investigated using an high-speed imaging technique. To our best knowledge, the moment of emulsification and formation of fine droplets during ultrasound irradiation were observed for the first time. It was found that the detachment of fine gallium droplets occurs from the water-gallium interface during collapse of big cavitation bubbles. The results suggest that the maximum size of cavitation bubble before collapsing is of prime importance for emulsification phenomena. Previous numerical simulation revealed that the collapse of big cavitation bubble is followed by generation of high-velocity liquid jet directed toward the water-gallium interface. Such a jet is assumed to be the prime cause of liquid emulsification. The distance between cavitation bubbles and water-gallium interface was found to slightly affect the emulsification onset. The droplet fragmentation conditions are also discussed in terms of the balance between (1) interfacial and kinetic energies and (2) dynamic and Laplace pressure during droplet formation.     

High-frequency acoustic emissions generated by a 20 kHz sonochemical horn processor detected using a novel broadband acoustic sensor: a preliminary study

This paper describes the application of a novel broadband acoustic sensor to evaluating the acoustic emissions from cavitation produced by a typical commercial 20 kHz sonochemical horn processor. Investigations of the reproducibility of the processor, and of the variation in cavitation emissions as a function of output setting and sensor location are described, and resulting trends discussed in terms of the broadband integrated power in the megahertz frequency range. Companion studies with a conventional membrane hydrophone have illustrated for the first time that cavitation emissions produced by a sonochemical horn processor can extend to frequencies beyond 20 MHz, and the sensor shows that significant nonlinearity can be seen in measured cavitation activity with increasing nominal output power.     

声空化工程初步*

声空化气泡内的高温、高压和高密度是声空化工程的机理和基础。该文简要回顾了国内外声空化理论和实验研究的进展,针对当前在液体中进行工业规模声处理等声能应用方面的现状和存在的问题,提出了扩大声能应用的一种可能途径——声空化空间分布控制并在实验室内进行了实验研究。结果表明该方法具有可行性。在此基础上,文中给出了基于空化空间分布控制方法在稠油井口辅助降粘和高固污泥预处理方面取得声空化工程初步应用的两个例子。现场的试验结果表明,这两种声空化工程化样机在实际现场工况条件下,均取得了较好的应用效果。最后,对当前声空化工程应用前景进行了初步探讨。     

碟形超声变幅杆的设计

优化超声变幅杆的形状结构可有效地提高水域声场分布和空化区域,提升对水域超声空化效果。通过模拟分析发现传统超声变幅杆在水域中具有声场分布均匀性差、变幅杆端部声压高等特征,不利于声波在水域中传播。基于此,提出并优化设计了一种具有碟形结构的变幅杆,位于变幅杆的最大振幅处的碟形结构,有更大的振动位移;模拟表明其水域声场和声压均衡度显著优于传统变幅杆,铝箔空化腐蚀实验进一步证实了其水域中的声压分布均匀性。同时,实验通过铝箔的空化腐蚀、KI剂量测定及工件表面油渍去除对比了传统变幅杆和碟形变幅杆,分析表明碟形变幅杆所在水域中有较大的空化腐蚀区域,腐蚀速率明显提升,声化学反应速率提高,油渍去除程度增强,说明了设计的碟形变幅杆能够促进空化泡的产生,增加水域空化区域。     

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

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

Radial oscillation and translational motion of a gas bubble in a micro-cavity

According to classical nucleation theory, a gas nucleus can grow into a cavitation bubble when the ambient pressure is negative. Here, the growth process of a gas nucleus in a micro-cavity was simplified to two “events”, and the full confinement effect of the surrounding medium of the cavity was considered by including the bulk modulus in the equation of state. The Rayleigh–Plesset-like equation of the cavitation bubble in the cavity was derived to model the radial oscillation and translational motion of the cavitation bubble in the local acoustic field. The numerical results show that the nucleation time of the cavitation bubble is sensitive to the initial position of the gas nucleus. The cavity size affects the duration of the radial oscillation of the cavitation bubble, where the duration is shorter for smaller cavities. The equilibrium radius of a cavitation bubble grown from a gas nucleus increases with increasing size of the cavity. There are two possible types of translational motion: reciprocal motion around the center of the cavity and motion toward the cavity wall. The growth process of gas nuclei into cavitation bubbles is also dependent on the compressibility of the surrounding medium and the magnitude of the negative pressure. Therefore, gas nuclei in a liquid cavity can be excited by acoustic waves to form cavitation bubbles, and the translational motion of the cavitation bubbles can be easily observed owing to the confining influence of the medium outside the cavity.     

Memory effect and redistribution of cavitation nuclei in a thin liquid layer

Temporal evolution and spatial distribution of acoustic cavitation structures in a thin liquid layer were investigated experimentally with high-speed photography. The inception and disappearance processes of cavitation bubble cloud revealed that the metastable cavitaton structures formed in the thin liquid layer caused a long-term “memory effect”. A factor which weakens the memory effect was identified. The distribution of cavitation nuclei was investigated by changing the temporal decay of the memory effect.     

Numerical simulation of the nonlinear ultrasonic pressure wave propagation in a cavitating bubbly liquid inside a sonochemical reactor

We investigate the acoustic wave propagation in bubbly liquid inside a pilot sonochemical reactor which aims to produce antibacterial medical textile fabrics by coating the textile with ZnO or CuO nanoparticles. Computational models on acoustic propagation are developed in order to aid the design procedures. The acoustic pressure wave propagation in the sonoreactor is simulated by solving the Helmholtz equation using a meshless numerical method. The paper implements both the state-of-the-art linear model and a nonlinear wave propagation model recently introduced by Louisnard (2012), and presents a novel iterative solution procedure for the nonlinear propagation model which can be implemented using any numerical method and/or programming tool. Comparative results regarding both the linear and the nonlinear wave propagation are shown. Effects of bubble size distribution and bubble volume fraction on the acoustic wave propagation are discussed in detail. The simulations demonstrate that the nonlinear model successfully captures the realistic spatial distribution of the cavitation zones and the associated acoustic pressure amplitudes.