Jet dynamics of a cavitation bubble near unequal-sized dual particles.

The jet dynamics of a cavitation bubble near unequal-sized dual particles is investigated employing OpenFOAM, and the effects of the jets on the particles are quantitatively analyzed in terms of their pressure impacts. Different from single-particle cases, the necks that evolve between dual particles are closely linked to the formation mechanism of the jets. Based on the simulation results, the jet dynamics can be divided into five scenarios: (1) the contraction of the annular depression produced by the collision of the two necks causes the bubble to split into two daughter bubbles and generates a single jet inside each daughter bubble; (2) the annular depression impacts the particle, leading to the bubble to fracture and producing a single jet inside a daughter bubble; (3) the bubble is split by a single neck constriction and produces a single jet; (4) the bubble is split by a single neck constriction and generates two jets; and (5) the bubble is split by the contraction of two necks and produces four jets together with three daughter bubbles. As the bubble-particle distance or the radius ratio of the dual particles increases, the maximum force on the small particle generated by the bubble decreases.

Investigation of mechanisms of shock wave generation by collapse of cavitation bubbles near particles.

The mechanism of generation of shock waves by the collapse of a cavitation bubble near a single particle or dual particles is numerically investigated using OpenFOAM. For the single-particle case, shock waves during bubble inception and jet impacting on the particle surface are revealed in detail. The pressure induced on the particle by the inception shock wave of the bubble decreases with increasing bubble-particle distance, and the pressure is proportional to 1/r1.26 (r being the distance from the center of the shock wave). For the dual particles, the evolution of the neck structure is closely related to the generation mechanism of the shock wave. At extremely close particle-bubble distances, two shock waves propagating in opposite directions are emitted outside and inside the bubble after two necks collide. At long particle-bubble distances, a shock wave is emitted after the neck contracts, and simultaneously the bubble splits into two daughter bubbles. The energy of the shock wave generated by the neck constriction (i.e., the pressure at its generation location) first increases and then decreases with increasing bubble-particle distance. For unequal-sized double particles, the size of the daughter bubble depends on the bubble-particle distance and the particle size. These findings provide new perspectives for understanding the damage sustained by hydro-mechanical components operating in sand-laden water flows.

A new model for bubble cluster dynamics in a viscoelastic media.

Bubble cluster dynamics in viscoelastic media is instructive for ultrasound diagnosis and therapy. In this paper, we propose a statistical model for bubble cluster dynamics in viscoelastic media considering the radius distribution of bubble nuclei. By investigating and comparing the response for a bubble in three conditions: single bubble; multi-bubble with the same radius; multi-bubble with different radius, the following rules are found: The promotion or suppression of the bubble cluster on the bubble vibration is not monotonous with the increase of the number of bubbles. The promotion or suppression of the bubble cluster on the bubble vibration varies alternately with the frequency. The effect of bubble cluster on bubble vibration is mostly suppressed when the driving acoustic pressure amplitude pa is high (5000 kPa). Usually, the bubble cluster promotes the vibration of the large bubbles (R0 = 10 µm) more, or suppresses it less.

Investigation on the effects of an elliptical wall on the dynamic behaviors of a bubble restricted by two parallel plates.

The present paper investigates the dynamic behaviors of a bubble restricted by two parallel plates near an elliptical wall. The typical experimental phenomena of the bubble are recorded employing the high-speed photography and a theoretical Kelvin impulse model is established. The impacts of the spatial position and the curvature of the wall on the bubble collapse behaviors are quantitatively investigated through the theoretical model and verified against the experimental results. The Kelvin impulse intensity and the direction during the bubble collapse process are compared and discussed for different elliptical-shaped walls. The main conclusions include: (1) During the bubble collapse process, the phenomenon of the bubble uneven splitting is discovered. (2) At different spatial positions and wall curvatures, the bubble collapse jet angle, movement distance, and velocity are in good agreement with the theoretical Kelvin impulse predictions. (3) As the short-to-long axis ratio increases, the differences in the distributions of the Kelvin impulse intensity and the direction near the elliptical wall gradually become larger, and the range of the influence of the impulse intensity expands.

Microbubble-enhanced transcranial MR-guided focused ultrasound brain hyperthermia: heating mechanism investigation using finite element method.

Recently, our group developed a synergistic brain drug delivery method to achieve simultaneous transcranial hyperthermia and localized blood-brain barrier opening via MR-guided focused ultrasound (MRgFUS). In a rodent model, we demonstrated that the ultrasound power required for transcranial MRgFUS hyperthermia was significantly reduced by injecting microbubbles (MBs). However, the specific mechanisms underlying the power reduction caused by MBs remain unclear. The present study aims to elucidate the mechanisms of MB-enhanced transcranial MRgFUS hyperthermia through numerical studies using the finite element method. The microbubble acoustic emission (MAE) and the viscous dissipation (VD) were hypothesized to be the specific mechanisms. Acoustic wave propagation was used to model the FUS propagation in the brain tissue, and a bubble dynamics equation for describing the dynamics of MBs with small shell thickness was used to model the MB oscillation under FUS exposures. A modified bioheat transfer equation was used to model the temperature in the rodent brain with different heat sources. A theoretical model was used to estimate the bubble shell's surface tension, elasticity, and viscosity losses. The simulation reveals that MAE and VD caused a 40.5% and 52.3% additional temperature rise, respectively. Compared with FUS only, MBs caused a 64.0% temperature increase, which is consistent with our previous animal experiments. Our investigation showed that MAE and VD are the main mechanisms of MB-enhanced transcranial MRgFUS hyperthermia.

Ammonia production by microbubbles: A theoretical analysis of achievable energy intensity.

The present paper studies the energy intensity of ammonia production by a freely oscillating microbubble placed in an infinite domain of liquid. The initial content of the bubble is a mixture of hydrogen and nitrogen. The bubble is expanded isothermically to a maximum radius, then it is "released" and oscillates freely. The input energy is composed of the potential energy of the bubble at the maximum radius, the energy required to produce hydrogen, and the pumping work in case a vacuum is employed. The chemical yield is computed by solving the underlying governing equations: the Keller-Miksis equation for the radial dynamics, the first law of thermodynamics for the internal temperature and the reaction mechanism for the evolution of the concentration of the chemical species. The control parameters during the simulations are the equilibrium bubble size, initial expansion ratio, ambient pressure, the initial concentration ratio of hydrogen and the material properties of the liquid. At the optimal parameter setup, the energy intensity is 90.17GJ/t that is 2.31 times higher than the best available technology, the Haber-Bosch process. In both cases, the hydrogen is generated via water electrolysis.

Corrosive effect on saturated pool boiling heat transfer characteristics of metallic surfaces with hierarchical micro/nano structures.

Surface modification is of critical interest to enhance boiling heat transfer in terms of heat transfer coefficient or critical heat flux (CHF), which is significantly affected by distinct surface morphology and wettability and it can improve the efficiency and safety of equipment. Furthermore, actual service environment may cause severe corrosion to the processed structured surfaces while its consequence on boiling heat transfer is still obscure. In this article, comprehensive researches are conducted to unravel corrosive effect on metallic samples made of stainless steel (SS) and Inconel materials with microstructures. Different constructions (i.e., microgroove, microcavity and micropillar array) and characteristic dimensions (∼20, 50 µm) of microstructure, various duration time (up to 300 days) and pH values (∼7.0-8.5) of corrosive environment are compared thoroughly. Conclusions can be drawn that not all microstructures can enhance pool boiling heat transfer characteristics, especially in terms of CHF values. More importantly, CHF value of SS microgroove sample firstly increases from 60.94 to 94.09 W·cm-2 in 50 days, then decreases to 47.77 W·cm-2 in 300 days, which can be attributed to competition result between formation of hierarchical micro/nano structure with enhancing wicking capability and chemistry condition with increasing contact angle. In addition, distinct bubble dynamics during pool boiling is also analyzed. The insights obtained from this article can be used to guide surface modification method and to reveal evolvement rule of engineered metallic surface in highly corrosive and harsh boiling heating transfer environment.

Numerical investigation of the effect of bubble properties on the linear resonance frequency shift due to inter-bubble interactions in ultrasonically excited lipid coated microbubbles.

Ultrasonically excited microbubbles (MBs) have numerous applications in various fields, such as drug delivery, and imaging. Ultrasonically excited MBs are known to be nonlinear oscillators that generate secondary acoustic emissions in the media when excited by a primary ultrasound wave. The propagation of acoustic waves in the liquid is limited to the speed of sound, resulting in each MB receiving the primary and secondary waves at different times depending on their distance from the ultrasound source and the distance between MBs. These delays are referred to as primary and secondary delays, respectively. A previous study demonstrated that the inclusion of secondary delays in a model describing the interactions between MBs exposed to ultrasound results in an increase in the linear resonance frequency of MBs as they approach each other. This work investigates the impact of various MB properties on the change in linear resonance frequency resulting from changes in inter-bubble distances. The effects of shell properties, including the initial surface tension, surface dilatational viscosity of the shell monolayer, elastic compression modulus of the shell, and the initial radius of the MBs, are examined. MB size is a significant factor influencing the rate of linear resonance frequency increase with increasing concentration. Moreover, it is found that the shell properties of MBs play a negligible role in the rate of change in linear resonance frequency of MBs as the inter-bubble distances change.The findings of this study have implications for various applications of MBs in the biomedical field. By understanding the impact of inter-bubble distances and shell properties on the linear resonance frequency of MBs, the utilization of MBs in applications reliant on their resonant behavior can be optimized.

On efficient modelling of radical production in cavitation assisted reactors.

Process intensification by cavitation is gaining widespread attention due to the benefits that the intense bubble collapse conditions can provide, yet, several knowledge gaps exist in the modelling of such systems. This work studies the numerical prediction of single bubble dynamics and the various approaches that can be employed to estimate the changes in the chemical composition of cavitating bubbles. Specific emphasis is placed on the prediction of the radical production rates during bubble collapse and the computational performance, with the aim of coupling the single bubble dynamics to flow models for reactor hydrodynamics. The results reveal that the choice of chemical reaction approach has virtually no effect on the bubble dynamics, whereas the predicted radical production rates can differ substantially. It is found that evaluating the radical production only on temperature peaks, an approach commonly followed in literature, may result in the most erroneous estimations (on average 12.8 times larger than those of the full kinetic model), while a simplified kinetic model yields more accurate predictions (2.3 times larger) at the expense of increased computational times. Continuous evaluation of the bubble content by assuming equilibrium when the bubble temperature is above a certain threshold (≈1500K) is shown to be capable of predicting total radical production values close to those estimated by solving the kinetics of a detailed reaction model (19.8% difference), as well as requiring only 22.2% more computational costs compared to simulations without chemical reaction modelling. Such an equilibrium approach is therefore recommended for future studies aiming to couple flow simulations with single bubble dynamics to accurately predict radical production rates in cavitation devices, involving numerous bubbles following different flow trajectories. Furthermore, an algebraic expression that successfully approximates the full kinetic simulation results is proposed as a function of the initial nucleus size and the time integral of the liquid pressure when it is under vapor pressure. Such a model can be applied in modelling efforts that do not require local instantaneous radical concentrations, and paves the way for efficient closure modelling of radical production in CFD simulations of hydrodynamic reactors.

Dynamics of cavitation bubbles in viscous liquids in a tube during a transient process.

We experimentally, numerically, and theoretically investigate the dynamics of cavitation bubbles in viscous liquids in a tube during a transient process. In experiments, cavitation bubbles are generated by a modified tube-arrest setup, and the bubble evolution is captured with high-speed imaging. Numerical simulations using OpenFOAM are employed to validate our quasi-one-dimensional theoretical model, which effectively characterizes the bubble dynamics. We find that cavitation onset is minimally affected by the liquid viscosity. However, once cavitation occurs, various aspects of bubble dynamics, such as the maximum bubble length, bubble lifetime, collapse time, and collapse speed, are closely related to the liquid viscosity. We further establish that normalized bubble dynamics are solely determined by the combination of the Reynolds number and the Euler number. Moreover, we also propose a new dimensionless number, Ca2, to predict the maximum bubble length, a critical factor in determining the occurrence of liquid column separation.

Ultra-high-speed dynamics of acoustic droplet vaporization in soft biomaterials: Effects of viscoelasticity, frequency, and bulk boiling point.

Phase-shift droplets are a highly adaptable platform for biomedical applications of ultrasound. The spatiotemporal response of phase-shift droplets to focused ultrasound above a certain pressure threshold, termed acoustic droplet vaporization (ADV), is influenced by intrinsic features (e.g., bulk boiling point) and extrinsic factors (e.g., driving frequency and surrounding media). A deep understanding of ADV dynamics is critical to ensure the robustness and repeatability of an ADV-assisted application. Here, we integrated ultra-high-speed imaging, at 10 million frames per second, and confocal microscopy for a full-scale (i.e., from nanoseconds to seconds) characterization of ADV. Experiments were conducted in fibrin-based hydrogels to mimic soft tissue environments. Effects of fibrin concentration (0.2 to 8 % (w/v)), excitation frequency (1, 2.5, and 9.4 MHz), and perfluorocarbon core (perfluoropentane, perfluorohexane, and perfluorooctane) on ADV dynamics were studied. Several fundamental parameters related to ADV dynamics, such as expansion ratio, expansion velocity, collapse radius, collapse time, radius of secondary rebound, resting radius, and equilibrium radius of the generated bubbles were extracted from the radius vs time curves. Diffusion-driven ADV-bubble growth was fit to a modified Epstein-Plesset equation, adding a material stress term, to estimate the growth rate. Our results indicated that ADV dynamics were significantly impacted by fibrin concentration, frequency, and perfluorocarbon liquid core. This is the first study to combine ultra-high-speed and confocal microscopy techniques to provide insights into ADV bubble dynamics in tissue-mimicking hydrogels.

Ultrasound for the degradation of endocrine disrupting compounds in aqueous solution: A review on mechanisms, influence of operating parameters and cost estimation.

Availability of drinking water is one of the basic humanitarian goals but remains as a grand challenge that the world is facing today. Currently, water bodies are contaminated not only with conventional pollutants but also with numerous recalcitrant pollutants, such as PPCPs, endocrine disrupting compounds, etc. These emerging pollutants require special attention because of their toxicity to living organisms, bio-resistant and can sustain even after primary and secondary treatments of wastewater. Among different treatment technologies, sonolysis is found to be an innovative and promising technique for the treatment of emerging pollutants present in aqueous solution. Sonolysis is the use of ultrasound to enhance or alter chemical reactions by the formation of free radicals and shock waves which ultimately helps in degradation of pollutants. This review summarizes several studies in the sonochemical literature, including mechanisms of sonochemical process, physical and chemical effects of ultrasound, and the influence of several process variables such as ultrasound frequency, power density, temperature and pH of the medium on degradation performance for endocrine disrupting compounds. In addition, this review highlighted techno-economic perspectives focusing on the total cost required for translating the ultrasound-based processes on a large scale. Overall, the objective of this study is to exhibit a critical review of information available in the literature to encourage and promote future research on sonolysis for the degradation of Endocrine Disrupting Compounds (EDCs).

Investigation of the ultrasound-induced collapse of air bubbles near soft materials.

A numerical investigation into the ultrasound-induced collapse of air bubbles near soft materials, utilizing a novel multi-material diffuse interface method (DIM) model with block-structured adaptive mesh refinement is presented. The present work expands from a previous five-equation DIM by incorporating Eulerian hyperelasticity. The model is applicable to any arbitrary number of interacting fluid and solid material. A single conservation law for the elastic stretch tensor enables tracking the deformations for all the solid materials. A series of benchmark cases are conducted, and the solution is found to be in excellent agreement against theoretical data. Subsequently, the ultrasound-induced bubble-tissue flow interactions are examined. The bubble radius was found to play a crucial role in dictating the stresses experienced by the tissue, underscoring its significance in medical applications. The results reveal that soft tissues primarily experience tensile forces during these interactions, suggesting potential tensile-driven injuries that may occur in relevant treatments. Moreover, regions of maximal tensile forces align with tissue elongation areas. It is documented that while early bubble dynamics remain relatively unaffected by changes in shear modulus, at later stages of the penetration processes and the deformation shapes, exhibit notable variations. Lastly, it is demonstrated that decreasing standoff distances enhances the interaction between bubbles and tissue, thereby increasing the stress levels in the tissue, although the behavior of the bubble dynamics remains largely unchanged.

An experimentally validated cavitation inception model for spring-driven autoinjectors.

Cavitation, the formation and collapse of vapor-filled bubbles, poses a problem in spring-driven autoinjectors (AIs). It occurs when the syringe accelerates abruptly during activation, causing pressure fluctuations within the liquid. These bubbles expand and then collapse, generating shock waves that can harm both the device and the drug molecules. This issue stems from the syringe's sudden acceleration when the driving rod hits the plunger. To better understand cavitation in AIs, we explore how design factors like drive spring force, air gap size, and fluid viscosity affect its likelihood and severity. We use a dynamic model for spring-driven autoinjectors to predict and analyze the factors contributing to cavitation initiation and severity. This model predicts the motion of AI components, such as the displacement and velocity of the syringe barrel, and allows us to investigate pressure wave propagation and the subsequent dynamics of cavitation under various operating conditions. We investigated different air gap heights (from 1 to 4 mm), drive spring forces (from 8 to 30 N), and drug solution viscosities (from 1 to 18 cp) to assess cavitation inception based on operational parameters. Results reveal that AI dynamics and cavitation onset and severity strongly depend upon AI operating parameters, namely drive spring force and air gap height. The maximum syringe acceleration increases with spring stiffness and decreases with air gap height; increases in air gap height prolong the time interval of syringe acceleration but diminish the maximum syringe acceleration. From actuation to injection, air gap pressure peaks twice, first due to impact with the rod/plunger and secondly due to the deacceleration event upon injection. The maximum air gap pressure increases with spring stiffness and decreases with air gap height. Results show that maximum cavitation bubble radii and collapse-driven extension rates occur with higher driver spring forces, smaller air gap heights, and less viscous solutions. A cavitation criterion is developed for cavitation in autoinjectors that concludes that cavitation in autoinjectors depends on the peak syringe acceleration.

Review on Bubble Dynamics in Proton Exchange Membrane Water Electrolysis: Towards Optimal Green Hydrogen Yield.

Water electrolysis using a proton exchange membrane (PEM) holds substantial promise to produce green hydrogen with zero carbon discharge. Although various techniques are available to produce hydrogen gas, the water electrolysis process tends to be more cost-effective with greater advantages for energy storage devices. However, one of the challenges associated with PEM water electrolysis is the accumulation of gas bubbles, which can impair cell performance and result in lower hydrogen output. Achieving an in-depth knowledge of bubble dynamics during electrolysis is essential for optimal cell performance. This review paper discusses bubble behaviors, measuring techniques, and other aspects of bubble dynamics in PEM water electrolysis. It also examines bubble behavior under different operating conditions, as well as the system geometry. The current review paper will further improve the understanding of bubble dynamics in PEM water electrolysis, facilitating more competent, inexpensive, and feasible green hydrogen production.

Experimental study on influence of particle shape on shockwave from collapse of cavitation bubble.

The bubble dynamics under the influence of particles is an unavoidable issue in many cavitation applications, with a fundamental aspect being the shockwave affected by particles during bubble collapse. In our experiments, the method of spark-induced bubbles was used, while a high-speed camera and a piezoresistive pressure sensor were utilized to investigate how particle shape affects the evolution of shockwaves. Through the high-speed photography, we found that the presence of the particle altered the consistency of the liquid medium around the bubble, which result in the emitting of water hammer shockwave and implosion shockwave respectively during the collapse of the bubble. This stratification effect was closely related to the bubble-particle relative distance φ and particle shape δ. Specifically, when the bubble-particle relative distance φ < 1.34 e-0.10δ, particles disrupted the medium consistency around the bubbles and led to a nonspherical collapse and the consequent stratification of the shockwave. By measuring the stratified shockwave intensity affected by different particle shapes, we found that the stratified shockwave intensity experienced varying degrees of attenuation. Furthermore, as the particle shape δ increased, the attenuation of the particle on shockwave intensity gradually reduced. These new findings hold significant theoretical implications for elucidating cavitation erosion mechanisms in liquid-solid two-phase flows and applications and prevention strategies in liquid-solid two-phase cavitation fields.

Desynchronization of neuronal firing in multiparameter ultrasound stimulation.

Low-intensity transcranial ultrasound stimulation, a novel neuromodulation technique, that possesses the advantages of non-invasiveness, high penetration depth, and high spatial resolution, has achieved positive neuromodulation effects in animal studies. But the regulatory mechanism remains controversial. The intramembrane cavitation effect is considered one of the mechanisms for ultrasound neuromodulation. In this study, the modified equations of ultrasonic cavitation bubble dynamics were coupled with the dual-coupled neuron Hindmarsh-Rose model, small-world neural network model, and the Jansen-Rit neural mass model, which simulate simple coupled neurons, complex neuronal networks, and discharge signals in epileptic disorders respectively. The results demonstrated that ultrasound stimulation has an appreciable modulatory effect on neuronal firing desynchronization in Hindmarsh-Rose model and small-world neural network model. The desynchronization effect is related to the stimulation frequency and intensity. Furthermore, ultrasound stimulation has an inhibitory effect on epileptic seizures, and the effect is enhanced by increasing ultrasound frequency from 0.1-1.0 MHz. This is the first combination of ultrasonic intramembrane cavitation effect theory with neurons and neural network firing desynchronization, which can provide guidance of parametric and theories support for the studies of neurological diseases such as epilepsy and Parkinson's disease.

Investigation of interaction effects on dual-frequency driven cavitation dynamics in a two-bubble system.

The cavitation dynamics of a two-bubble system in viscoelastic media excited by dual-frequency ultrasound is studied numerically with a focus on the effects of inter-bubble interactions. Compared to the isolated bubble cases, the enhancement or suppression effects can be exerted on the amplitude and nonlinearity of the bubble oscillations to different degrees. Moreover, the interaction effects are found to be highly sensitive to multiple paramount parameters related to the two-bubble system, the dual-frequency ultrasound and the medium viscoelasticity. Specifically, the larger bubble of a two-bubble system shows a stronger effect on the smaller one, and this effect becomes more pronounced when the larger bubble undergoes harmonic and/or subharmonic resonances as well as the two bubbles get closer (e.g., d0 < 100 µm). For the influences of the dual-frequency excitation, the results show that the bubbles can achieve enhanced harmonic and/or subharmonic oscillations as the frequency combinations with small frequency differences (e.g., Δf < 0.2 MHz) close to the corresponding resonance frequencies of bubbles, and the interaction effects are consequently intensified. Similarly, the bubble oscillations and the interaction effects can also be enhanced as the acoustic pressure amplitude of each frequency component is equal and the pressure amplitude pA increases. Above a pressure threshold (pA = 215 kPa), a larger bubble undergoes period 2 (P2) oscillations, which can force a smaller bubble to change its oscillation pattern from period 1 (P1) into P2 oscillations. In addition, it is found that the medium viscosity dampens the bubble oscillations while the medium elasticity affects the bubble resonances, accordingly exhibiting stronger interaction effects at smaller viscosities (e.g., µ < 4 mPa·s) or certain elasticities (approximately G = 70-120 kPa, G = 160-200 kPa and G = 640-780 kPa) at which the bubble resonances occur. The study can contribute to a better understanding of the complex dynamic behaviors of interacting cavitation bubbles in viscoelastic tissues for high efficient cavitation-mediated biomedical applications using dual-frequency ultrasound.

Numerical and experimental investigations on the jet and shock wave dynamics during the cavitation bubble collapsing near spherical particles based on OpenFOAM.

The interaction between cavitation bubbles and particles is essential for the operational performance many kinds of fluid machineries. In the present paper, jet dynamics and shock waves induced by the cavitation bubble collapsing near two spherical particles are numerically investigated based on OpenFOAM. The numerical scheme is validated by the experimental data obtained based on our high-speed camera cavitation system. Our results reveal that bubble split induced by annular jet is the primary feature during bubble collapsing with four typical cases defined. For the jet formation, the localized high pressure produced at the bubble split point is the main reason and the split point also serves as the source of the shock waves shown by the numerical schlieren. Furthermore, the nondimensional bubble-particle distance is the most paramount parameter influencing the jet phenomenon (e.g. jet velocity).

Bubble dynamics in a pressure gradient with reentrant jet break through and energy loss.

The dynamics of a bubble in a pressure gradient is investigated experimentally and numerically with particular emphasis on the behavior at reentrant jet impact and break through the opposite side of the bubble with corresponding energy loss and vorticity generation. High speed photography observations of a bubble generated by electric spark energy deposit in a low ambient pressure tank are coupled with wavelet based Optical Flow Velocimetry (wOFV) and Boundary Element Method (BEM) numerical analysis to examine the flow field resulting from the reentrant jet formation and break through. We study, as an illustration, the effects of the constant pressure gradient due to gravity on the bubble dynamics. Energy losses between the first and second cycle are measured for the bubbles generated under various conditions characterized by a non-dimensional pressure gradient parameter, and the corresponding measured energy loss is used in the numerical simulations. Good correspondence is seen between the image analysis, the wOFV computations, and the BEM results and insight is gained on the involved physics.