A boundary element model of the transport of a semi-infinite bubble through a microvessel bifurcation

Evolution of Thin-Liquid Films Surrounding Bubbles in Microfluidics and Their Impact on the Pressure Drop and Fluid Movement

The evolution of thin-liquid films in a microchannel is one of the most critical and intricate phenomena to understand two-phase movement, evaporation, micromixing, heat transfer, chemical synthesis, biological processes, and efficient energy devices. In this paper, we demonstrate experimentally the effect of a liquid film on the removal of an initially dry and lodged bubble in laser-etched poly(methyl methacrylate) microfluidic networks and discuss the evolution of the liquid film in accordance with the bubble superficial velocity and the effect of liquid properties and branch angle on the evolution of the liquid film and the pressure drop. During the removal of a dry bubble, four stages have been observed in the bubble velocity profile and they directly relate to the evolution of the liquid film. The correlation of maximum bubble velocity has been derived as a function of bubble length, fluid viscosity, surface tension, geometry of the cross-sectional area, and dimensions of the microchannel and agrees with the experimental results. The bubble moving distance required for the full deposition of a continuous and stable thin-liquid film is affected by the liquid viscosity and network branch angle. The liquid with a higher viscosity will increase the pressure drop for removing dry bubbles from microfluidic networks, while this effect will be hampered by increasing the microfluidic network complexity. The deposition of the thin-liquid film surrounding bubbles significantly decreases the pressure drop required to remove bubbles from microfluidics. Compared with deionized water, the glycerol solution is prone to acting as the lubricating liquid due to its strong H-bond interaction with the channel wall and the reduction in interfacial energy of the gas-water interface.

Combined gas embolization and chemotherapy can result in complete tumor regression in a murine hepatocellular carcinoma model

Hepatocellular carcinoma (HCC) is an intractable cancer with a high mortality rate. Transarterial chemoembolization (TACE), a non-curative method, is the first line therapy for intermediate stage patients. This effectively extends patient survival but requires a complicated intraarterial catheterization procedure and is poorly suited to repeated administration. Here, we investigate gas chemoembolization, a less invasive, more easily administered transient occlusion method that circumvents these limitations. We examined the efficacy of repeated embolization combined with systemically administered doxorubicin, the most common chemotherapeutic in TACE, or tirapazamine, a hypoxia-activated cytotoxic agent, in an ectopic xenograft model of HCC. Emboli were generated in situ using acoustic droplet vaporization, the noninvasive focused ultrasound-mediated conversion of intravenously administered perfluorocarbon microdroplets into microbubbles. Gas embolization alone significantly reduced the Ki67 index and tumor viability (11.6 ± 6.71% non-necrotic vs 100% in control; p < 0.01) after 3 treatments, as assessed by histological analysis. Mice treated for three weeks exhibited significant tumor regression compared to control (23.8 ± 5.37% of initial volume vs 427 ± 49.7% in controls, p < 0.01), irrespective of the chosen chemotherapeutic agent. However, an additional three weeks of monitoring post-treatment elucidated a significant difference in the tumor recurrence rate, with combined gas embolization and doxorubicin resulting in the best treatment outcomes (60% complete regression). While doxorubicin administration resulted in significant cardiotoxicity (p < 0.01), it strongly interacted with the droplet shells, reducing the systemic dose by 11.4%. Overall, gas chemoembolization shows promise as a developmental therapy and merits further study in more complex tumor models.

Minimally invasive gas embolization using acoustic droplet vaporization in a rodent model of hepatocellular carcinoma

Hepatocellular carcinoma is the third leading cause of cancer-related deaths worldwide. Many patients are not eligible for curative therapies, such as surgical resection of the tumor or a liver transplant. Transarterial embolization is one therapy clinically used in these cases; however, this requires a long procedure and careful placement of an intraarterial catheter. Gas embolization has been proposed as a fast, easily administered, more spatially selective, and less invasive alternative. Here, we demonstrate the feasibility and efficacy of using acoustic droplet vaporization to noninvasively generate gas emboli within vasculature. Intravital microscopy experiments were performed using the rat cremaster muscle to visually observe the formation of occlusions. Large gas emboli were produced within the vasculature in the rat cremaster, effectively occluding blood flow. Following these experiments, the therapeutic efficacy of gas embolization was investigated in an ectopic xenograft model of hepatocellular carcinoma in mice. The treatment group exhibited a significantly lower final tumor volume (ANOVA, p = 0.008) and growth rate than control groups - tumor growth was completely halted. Additionally, treated tumors exhibited significant necrosis as determined by histological analysis. To our knowledge, this study is the first to demonstrate the therapeutic efficacy of gas embolotherapy in a tumor model.

Bubble Dislodgment in a Capillary Network with Microscopic Multichannels and Multibifurcation Features

Bubble lodgment in a complex capillary network is a common issue in many industrial and biological processes. Research work reported in the literature only investigated bubble dislodgment in single channels and did not consider the effect of network complexity on the dislodgment. This paper focuses on the pressure required to dislodge single bubbles from a microscopic capillary network and investigates the factors affecting the dislodging pressure to facilitate the precise control of bubble flows in porous media. A capillary network with multibifurcation and a smoothly changed diameter is designed to closely mimic the structure of the physiological vascular networks. Over 600 bubble dislodgment experiments have been conducted to understand the effect of the network structure, channel dimensions, and bubble length on the dislodging pressure. The results indicate that the network structure is a dominant factor affecting the dislodging pressure that increases with the increase in network complexity. The effect of bubble length on the dislodging pressure depends on the bubble length. When the bubble length is less than a certain value, which is around 2 mm in this study, the dislodging pressure increases significantly with the decrease of bubble length. When the bubble length is larger than 2 mm, the dislodging pressure is independent of the bubble length. A model has been proposed to explain the bubble dislodgment in complex capillary networks. The impact of the network structure on the bubble dislodging pressure is characterized by a parameter c _j. The model indicates that the dislodging pressure is the function of bubble length, channel dimension, and network structure. The analysis of model parameters NB _j and MA _j shows that parameter c _j, rather than the channel size, dominates the dislodging pressure for bubbles with a length greater than 2 mm, and the increase rate of the dislodging pressure is significantly affected by both channel size and parameter c _j.

Dynamics of partially obstructed breakup of bubbles in microfluidic Y-junctions

For revealing the dynamics of partially obstructed breakup of bubbles in microfluidic Y-junctions, the combination of dimensionless power-law and geometric model was applied to study the effects of capillary number, bubble length, and channel angle on the bubble rupture process. In the squeezing process, the gas-liquid interface curve follows the parabolic model. For the evolution of the bubble neck during breakup, the increase of the bubble length, the channel angle, and the capillary number leads to the decrease of the focus distance α. The chord m increases with the increase of the capillary number and the decrease of the bubble length, and it would reach the maximum value when the channel angle is 90°. In the fast pinch-off stage during bubble breakup, the bubble's neck curve no longer conforms to the parabolic model so the focus and chord no longer exist. For the evolution of the bubble head during breakup, the value of γ approaches 1 with the increase of the capillary number and the bubble length, and with the close of the channel angle to 90°. It is found that the quadrilateral model can be applied for the partially obstructed rupture of bubbles in the symmetrical microfluidic Y-junction.

Bubbling and foaming assisted clearing of mucin plugs in microfluidic Y-junctions

Microfluidic Y-junctions were used to study mechanical mechanisms involved in pig gastric mucin (PGM) plug removal from within one of two bifurcation branches with 2-phase air and liquid flow. Water control experiments showed moderate plug removal due to shear from vortex formation in the blockage branch and suggest a PGM yield stress of 35Pa, as determined by computational fluid dynamics. Addition of hexadecyltrimethylammonium bromide (CTAB) surfactant improved clearing effectiveness due to bubbling in 1mm diameter channels and foaming in 500μm diameter channels. Plug removal mechanisms have been identified as vortex shear, bubble scouring, and then foam scouring as air flow rate is increased with constant liquid flow. The onset of bubbling and foaming is attributed to a flow regime transition from slug to slug-annular. Flow rates explored for 1mm channels are typically experienced by bronchioles in generations 8 and 9 of lungs. Results have implications on treatment of cystic fibrosis and other lung diseases.

Small-bubble transport and splitting dynamics in a symmetric bifurcation

Simulations of small bubbles traveling through symmetric bifurcations are conducted to garner information pertinent to gas embolotherapy, a potential cancer treatment. Gas embolotherapy procedures use intra-arterial bubbles to occlude tumor blood supply. As bubbles pass through bifurcations in the blood stream nonhomogeneous splitting and undesirable bioeffects may occur. To aid development of gas embolotherapy techniques, a volume of fluid method is used to model the splitting process of gas bubbles passing through artery and arteriole bifurcations. The model reproduces the variety of splitting behaviors observed experimentally, including the bubble reversal phenomenon. Splitting homogeneity and maximum shear stress along the vessel walls is predicted over a variety of physical parameters. Small bubbles, having initial length less than twice the vessel diameter, were found unlikely to split in the presence of gravitational asymmetry. Maximum shear stresses were found to decrease exponentially with increasing Reynolds number. Vortex-induced shearing near the bifurcation is identified as a possible mechanism for endothelial cell damage.

Computational analysis of microbubble flows in bifurcating airways: role of gravity, inertia, and surface tension

Although mechanical ventilation is a life-saving therapy for patients with severe lung disorders, the microbubble flows generated during ventilation generate hydrodynamic stresses, including pressure and shear stress gradients, which damage the pulmonary epithelium. In this study, we used computational fluid dynamics to investigate how gravity, inertia, and surface tension influence both microbubble flow patterns in bifurcating airways and the magnitude/distribution of hydrodynamic stresses on the airway wall. Direct interface tracking and finite element techniques were used to simulate bubble propagation in a two-dimensional (2D) liquid-filled bifurcating airway. Computational solutions of the full incompressible Navier-Stokes equation were used to investigate how inertia, gravity, and surface tension forces as characterized by the Reynolds (Re), Bond (Bo), and Capillary (Ca) numbers influence pressure and shear stress gradients at the airway wall. Gravity had a significant impact on flow patterns and hydrodynamic stress magnitudes where Bo > 1 led to dramatic changes in bubble shape and increased pressure and shear stress gradients in the upper daughter airway. Interestingly, increased pressure gradients near the bifurcation point (i.e., carina) were only elevated during asymmetric bubble splitting. Although changes in pressure gradient magnitudes were generally more sensitive to Ca, under large Re conditions, both Re and Ca significantly altered the pressure gradient magnitude. We conclude that inertia, gravity, and surface tension can all have a significant impact on microbubble flow patterns and hydrodynamic stresses in bifurcating airways.

Phase change events of volatile liquid perfluorocarbon contrast agents produce unique acoustic signatures

Phase-change contrast agents (PCCAs) provide a dynamic platform to approach problems in medical ultrasound (US). Upon US-mediated activation, the liquid core vaporizes and expands to produce a gas bubble ideal for US imaging and therapy. In this study, we demonstrate through high-speed video microscopy and US interrogation that PCCAs composed of highly volatile perfluorocarbons (PFCs) exhibit unique acoustic behavior that can be detected and differentiated from standard microbubble contrast agents. Experimental results show that when activated with short pulses PCCAs will over-expand and undergo unforced radial oscillation while settling to a final bubble diameter. The size-dependent oscillation phenomenon generates a unique acoustic signal that can be passively detected in both time and frequency domain using confocal piston transducers with an 'activate high' (8 MHz, 2 cycles), 'listen low' (1 MHz) scheme. Results show that the magnitude of the acoustic 'signature' increases as PFC boiling point decreases. By using a band-limited spectral processing technique, the droplet signals can be isolated from controls and used to build experimental relationships between concentration and vaporization pressure. The techniques shown here may be useful for physical studies as well as development of droplet-specific imaging techniques.

Numerical Simulation of Bubble Transport in a Bifurcating Microchannel: A Preliminary Study

In this paper, we present the computational fluid dynamics (CFD) simulations of bubble transport in a first generation bifurcating microchannel. In the present study, the human arteriole is modeled as a two-dimensional (2D) rectangular bifurcating microchannel. The microchannel is filled with blood and a single perfluorocarbon (PFC) bubble is introduced in the parent channel. The simulations are carried out to identify the lodging and dislodging pressures for two nondimensional bubble sizes, Ld (ratio of the dimensional bubble length to the parent tube diameter), that is for Ld = 1 and Ld = 2. Subsequently, the bubble transport and splitting behavior due to the presence of symmetry and asymmetry in the daughter channels of the microchannel is studied for these bubble sizes. The splitting behavior of the bubble under the effect of gravity is also assessed and reported here. For the symmetric bifurcation model, the splitting ratio (SR) (ratio of bubble volume in bottom daughter channel to bubble volume in top daughter channel), of the bubble was found to be 1. For the asymmetric model, the splitting ratio was found to be less than 1. The loss in the bubble volume in the asymmetric model was attributed to surface tension effects and the resistance offered by the flow, which led to the bubble sticking and sliding along the walls of the channel. With the increase in roll angle, Φ (angle which the plane makes with the horizontal to study the effects of gravity), there was a decline in the splitting ratio.

Phase-change contrast agents for imaging and therapy

Phase-change contrast agents (PCCAs) for ultrasound-based applications have resulted in novel ways of approaching diagnostic and therapeutic techniques beyond what is possible with microbubble contrast agents and liquid emulsions. When subjected to sufficient pressures delivered by an ultrasound transducer, stabilized droplets undergo a phase-transition to the gaseous state and a volumetric expansion occurs. This phenomenon, termed acoustic droplet vaporization, has been proposed as a means to address a number of in vivo applications at the microscale and nanoscale. In this review, the history of PCCAs, physical mechanisms involved, and proposed applications are discussed with a summary of studies demonstrated in vivo. Factors that influence the design of PCCAs are discussed, as well as the need for future studies to characterize potential bioeffects for administration in humans and optimization of ultrasound parameters.

Microbubble transport through a bifurcating vessel network with pulsatile flow

Motivated by two-phase microfluidics and by the clinical applications of air embolism and a developmental gas embolotherapy technique, experimental and theoretical models of microbubble transport in pulsatile flow are presented. The one-dimensional time-dependent theoretical model is developed from an unsteady Bernoulli equation that has been modified to include viscous and unsteady effects. Results of both experiments and theory show that roll angle (the angle the plane of the bifurcating network makes with the horizontal) is an important contributor to bubble splitting ratio at each bifurcation within the bifurcating network. When compared to corresponding constant flow, pulsatile flow was shown to produce insignificant changes to the overall splitting ratio of the bubble despite the order one Womersley numbers, suggesting that bubble splitting through the vasculature could be modeled adequately with a more modest constant flow model. However, bubble lodging was affected by the flow pulsatility, and the effects of pulsatile flow were evident in the dependence of splitting ratio of bubble length. The ability of bubbles to remain lodged after reaching a steady state in the bifurcations is promising for the effectiveness of gas embolotherapy to occlude blood flow to tumors, and indicates the importance of understanding where lodging will occur in air embolism. The ability to accurately predict the bubble dynamics in unsteady flow within a bifurcating network is demonstrated and suggests the potential for bubbles in microfluidics devices to encode information in both steady and unsteady aspects of their dynamics.