Interface capturing simulations of acoustically driven bubble dynamics in and near tissue

Tissue injury during therapeutic ultrasound or lithotripsy is thought, in cases, to be due to the action of cavitation bubbles. Assessing this and mitigating it is challenging since bubble dynamics in the complex confinement of tissues or in small blood vessels are challenging to predict. Simulations tools require specialized algorithms to simultaneously represent strong acoustic waves and shocks, topologically complex liquid‐vapor phase boundaries, and the complex viscoelastic material dynamics of tissue. We discuss advances in a simulation tool for such situations. A single‐mesh Eulerian solver is used to solve the governing equations. Special sharpening terms maintain the liquid‐vapor interface in face of the finite numerical dissipation included in the scheme to accurately capture shocks. A recent enhancement to this formulation has significantly improved this interface capturing procedure, which is demonstrated for simulation of the Rayleigh collapse of a bubble. The solver also transports elastic stresses and can thus be used to assess the effects of elastic properties on bubble dynamics. A shock‐induced bubble collapse adjacent to a model elastic tissue is used to demonstrate this and draw some conclusions regarding the injury suppressing role that tissue elasticity might play.

Bubble Dynamics In Nucleate Pool Boiling On Thin Wires In Microgravity

A temperature-controlled pool boiling (TCPB) device has been developed to study the bubble behavior and heat transfer in pool boiling phenomenon both in normal gravity and in microgravity. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. The fluid is R113 at 0.1 MPa and subcooled by 26 degrees C nominally for all cases. Three modes of heat transfer, namely single-phase natural convection, nucleate boiling, and two-mode transition boiling, are observed in the experiment both in microgravity aboard the 22nd Chinese recoverable satellite and in normal gravity on the ground before and after the space flight. Dynamic behaviors of vapor bubbles observed in these experiments are reported and analyzed in the present paper. In the regime of fully developed nucleate boiling, the interface oscillation due to coalescence of adjacent tiny bubbles is the primary reason of the departure of bubbles in microgravity. On the contrary, in the discrete bubble regime, it's observed that there exist three critical bubble diameters in microgravity, dividing the whole range of the observed bubbles into four regimes. Firstly, tiny bubbles are continually forming and growing on the heating surface before departing slowly from the wire when their sizes exceed some value of the order of 10(-1) mm. The bigger bubbles with about several millimeters in diameter stay on the wire, oscillate along the wire, and coalesce with adjacent bubbles. The biggest bubble with diameter of the order of 10 mm, which was formed immediately after the onset of boiling, stays continuously

Effect of pulse heating parameters on the microscale bubble dynamics at a microheater surface

We study the effects of pulse heating parameters on the micro bubble behavior of a platinum microheater (100 mu m x 20 mu m) immersed in a methanol pool. The experiment covers the heat fluxes of 10-37 MW/m(2) and pulse frequencies of 25-500 Hz. The boiling incipience is initiated at the superheat limit of methanol, corresponding to the homogeneous nucleation. Three types of micro boiling patterns are identified. The first type is named as the bubble explosion and regrowth, consisting of a violent explosive boiling and shrinking, followed by a slower bubble regrowth and subsequent shrinking, occurring at lower heat fluxes. The second type, named as the bubble breakup and attraction, consists of the violent explosive boiling, bubble breakup and emission, bubble attraction and coalescence process, occurring at higher heat fluxes than those of the first type. The third type, named as the bubble size oscillation and large bubble formation, involves the initial explosive boiling, followed by a short periodic bubble growth and shrinking. Then the bubble continues to increase its size, until a constant bubble size is reached which is larger than the microheater length. 

Bubble Dynamics and Heat Transfer in Microgravity Pool Boiling

Boiling is an extremely complicated and illusive process. Microgravity experiments offer a unique opportunity to study the complex interactions without external forces, such as buoyancy, which can affect the bubble dynamics and the related heat transfer. Furthermore, they can also provide a means to study the actual influence of gravity on the boiling. Two research projects on pool boiling in microgravity have been conducted aboard the Chinese recoverable satellites. Ground-based experiments both in normal gravity and in short-term microgravity in the Drop Tower Beijing and numerical simulations have also been performed. Steady boiling of R113 on thin platinum wires was studied with a temperature-controlled heating method, while quasi-steady boiling of FC-72 on a plane plate was investigated with an exponentially increasing heating voltage. It was found that the bubble dynamics in microgravity has a distinct difference from that in normal gravity, and that the heat transfer characteristic is depended upon the bubble dynamics. Lateral motions of bubbles on the heaters were observed before their departure in microgravity. The surface oscillation of the merged bubbles due to lateral coalescence between adjacent bubbles drove it to detach from the heaters. Slight enhancement of heat transfer on wires is observed in microgravity, while diminution is evident for high heat flux in the plate case.

INVESTIGATION OF BUBBLE DYNAMICS AND HEATING DURING FOCUSED ULTRASOUND INSONATION IN TISSUE-MIMICKING MATERIALS

The deposition of ultrasonic energy in tissue can cause tissue damage due to local heating. For pressures above a critical threshold, cavitation will occur in tissue and bubbles will be created. These oscillating bubbles can induce a much larger thermal energy deposition in the local region. Traditionally, clinicians and researchers have not exploited this bubble-enhanced heating since cavitation behavior is erratic and very difficult to control. The present work is an attempt to control and utilize this bubble-enhanced heating. First, by applying appropriate bubble dynamic models, limits on the asymptotic bubble size distribution are obtained for different driving pressures at 1 MHz. The size distributions are bounded by two thresholds: the bubble shape instability threshold and the rectified diffusion threshold. The growth rate of bubbles in this region is also given, and the resulting time evolution of the heating in a given insonation scenario is modeled. In addition, some experimental results have been obtained to investigate the bubble-enhanced heating in an agar and graphite based tissue- mimicking material. Heating as a function of dissolved gas concentrations in the tissue phantom is investigated. Bubble-based contrast agents are introduced to investigate the effect on the bubble-enhanced heating, and to control the initial bubble size distribution. The mechanisms of cavitation-related bubble heating are investigated, and a heating model is established using our understanding of the bubble dynamics. By fitting appropriate bubble densities in the ultrasound field, the peak temperature changes are simulated. The results for required bubble density are given. Finally, a simple bubbly liquid model is presented to estimate the shielding effects which may be important even for low void fraction during high intensity focused ultrasound (HIFU) treatment.

EXPERIMENTAL INVESTIGATION OF LINEAR BUBBLE DYNAMICS IN A VISCOELASTIC XANTHAN GEL

Oceanic bubble plumes caused by ship wakes or breaking waves disrupt sonar communi- cation because of the dramatic change in sound speed and attenuation in the bubbly fluid. Experiments in bubbly fluids have suffered from the inability to quantitatively characterize the fluid because of continuous air bubble motion. Conversely, single bubble experiments, where the bubble is trapped by a pressure field or stabilizing object, are limited in usable frequency range, apparatus complexity, or the invasive nature of the stabilizing object (wire, plate, etc.). Suspension of a bubble in a viscoelastic Xanthan gel allows acoustically forced oscilla- tions with negligible translation over a broad frequency band. Assuming only linear, radial motion, laser scattering from a bubble oscillating below, through, and above its resonance is measured. As the bubble dissolves in the gel, different bubble sizes are measured in the range 240 – 470 μm radius, corresponding to the frequency range 6 – 14 kHz. Equalization of the cell response in the raw data isolates the frequency response of the bubble. Compari- son to theory for a bubble in water shows good agreement between the predicted resonance frequency and damping, such that the bubble behaves as if it were oscillating in water.

Benthic megafauna and CO2 bubble dynamics observed by underwater photography during a controlled sub-seabed release of CO2

In 2012, a controlled sub-seabed release of carbon dioxide (CO2) was conducted in Ardmucknish Bay, a shallow (12 m) coastal bay on the west coast of Scotland. During the experiment, CO2 gas was released 12 m below the seabed for 37 days, causing significant disruption to sediment and water carbonate chemistry as the gas passed up through the sediment and into the overlying water. One of the aims of the study was to investigate how the impacts caused by leakage from geological CO2 Capture and Storage (CCS) could be detected and quantified in the context of natural heterogeneity and dynamics. To do this underwater photography was used to analyze (i) the benthic megafaunal response to the CO2 release and (ii) the dynamics of the CO2 bubble streams, emerging from the seabed into the overlying water column. The frequently observed megafauna species in the study area were Virgularia mirabilis (Cnidaria), Turritella communis (Mollusca), Asterias rubens (Echinodermata), Pagurus bernhardus (Crustacea), Liocarcinus depurator (Crustacea), and Gadus morhua (Osteichthyes). No discernable abnormal behavior was observed for these megafauna, in any of the zones investigated, during or after the CO2 release. Time-lapse photography revealed that the intensity and presence of the CO2 bubble plume was affected by the tides, with the most active bubbling seen at low tides and the larger hydrostatic pressure at high tide suppressing CO2 bubbling from the seabed.

Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics

The interaction of a single bubble with a single vortex ring in water has been studied experimentally. Measurements of both the bubble dynamics and vorticity dynamics have been done to help understand the two-way coupled problem. The circulation strength of the vortex ring (Gamma) has been systematically varied, while keeping the bubble diameter (D-b) constant, with the bubble volume to vortex core volume ratio (V-R) also kept fixed at about 0.1. The other important parameter in the problem is a Weber number based on the vortex ring strength. (We = 0.87 rho(Gamma/2 pi a)(2)/(sigma/D-b); a = vortex core radius, sigma = surface tension), which is varied over a large range, We = 3-406. The interaction between the bubble and ring for each of the We cases broadly falls into four stages. Stage I is before capture of the bubble by the ring where the bubble is drawn into the low-pressure vortex core, while in stage II the bubble is stretched in the azimuthal direction within the ring and gradually broken up into a number of smaller bubbles. Following this, in stage III the bubble break-up is complete and the resulting smaller bubbles slowly move around the core, and finally in stage IV the bubbles escape. Apart from the effect of the ring on the bubble, the bubble is also shown to significantly affect the vortex ring, especially at low We (We similar to 3). In these low-We cases, the convection speed drops significantly compared to the base case without a bubble, while the core appears to fragment with a resultant large decrease in enstrophy by about 50 %. In the higher-We cases (We > 100), there are some differences in convection speed and enstrophy, but the effects are relatively small. The most dramatic effects of the bubble on the ring are found for thicker core rings at low We (We similar to 3) with the vortex ring almost stopping after interacting with the bubble, and the core fragmenting into two parts. The present idealized experiments exhibit many phenomena also seen in bubbly turbulent flows such as reduction in enstrophy, suppression of structures, enhancement of energy at small scales and reduction in energy at large scales. These similarities suggest that results from the present experiments can be helpful in better understanding interactions of bubbles with eddies in turbulent flows.

Dynamics of discrete bubble in nucleate pool boiling on thin wires in microgravity

A space experiment on bubble behavior and heat transfer in subcooled pool boiling phenomenon has been performed utilizing the temperature-controlled pool boiling (TCPB) device both in normal gravity in the laboratory and in microgravity aboard the 22(nd) Chinese recoverable satellite. The fluid is degassed R113 at 0.1 MPa and subcooled by 26 degrees C nominally. A thin platinum wire of 60 mu m in diameter and 30 mm in length is simultaneously used as heater and thermometer. Only the dynamics of the vapor bubbles, particularly the lateral motion and the departure of discrete vapor bubbles in nucleate pool boiling are reported and analyzed in the present paper. It's found that these distinct behaviors can be explained by the Marangoni convection in the liquid surrounding vapor bubbles. The origin of the Marangoni effect is also discussed.

Suppression of air bubble in microfluidic channel using electro-acoustic excitation

A simple method to study the air bubble dynamics and to burst the air bubbles formed on the electrode– electrolyte interface in a parallel gate electrode fluidic channel is demonstrated. Upon application of a voltage across the electrodes,volume of water contained between them begins to electrolyzing depending on the conductivity, as well as it boils due to heating effect. This results in bubble formation within. These bubbles grow in radius with higher potential difference applied across the electrodes. As an approach towards removing these bubbles, an alternating current is applied at low potential difference of a 5 volts and high frequency at few megahertz. The alternating electric field had a heating effect on the bubbles where the energy input due to current heats up water and bursts the bubble. The bubbles of size up to 480μm were burst at 2500 V/m using this approach.

Model Order Reduction for Determining Bubble Parameters to Attain a Desired Fluid Surface Shape

In this paper, a new methodology for predicting fluid free surface shape using Model Order Reduction (MOR) is presented. Proper Orthogonal Decomposition combined with a linear interpolation procedure for its coefficient is applied to a problem involving bubble dynamics near to a free surface. A model is developed to accurately and efficiently capture the variation of the free surface shape with different bubble parameters. In addition, a systematic approach is developed within the MOR framework to find the best initial locations and pressures for a set of bubbles beneath the quiescent free surface such that the resultant free surface attained is close to a desired shape. Predictions of the free surface in two-dimensions and three-dimensions are presented.