Gas Bubble Dynamics under Mechanical Vibrations

NASA Astrophysics Data System (ADS)

Mohagheghian, Shahrouz; Elbing, Brian

2017-11-01

The scientific community has a limited understanding of the bubble dynamics under mechanical oscillations due to over simplification of Navier-Stockes equation by neglecting the shear stress tensor and not accounting for body forces when calculating the acoustic radiation force. The current work experimental investigates bubble dynamics under mechanical vibration and resulting acoustic field by measuring the bubble size and velocity using high-speed imaging. The experimental setup consists of a custom-designed shaker table, cast acrylic bubble column, compressed air injection manifold and an optical imaging system. The mechanical vibrations resulted in accelerations between 0.25 to 10 times gravitational acceleration corresponding to frequency and amplitude range of 8 - 22Hz and 1 - 10mm respectively. Throughout testing the void fraction was limited to <5%. The bubble size is larger than resonance size and smaller than acoustic wavelength. The amplitude of acoustic pressure wave was estimated using the definition of Bjerknes force in combination with Rayleigh-Plesset equation. Physical behavior of the system was capture and classified. Bubble size, velocity as well as size and spatial distribution will be presented.

Bubble-Induced Color Doppler Feedback for Histotripsy Tissue Fractionation.

PubMed

Miller, Ryan M; Zhang, Xi; Maxwell, Adam D; Cain, Charles A; Xu, Zhen

2016-03-01

Histotripsy therapy produces cavitating bubble clouds to increasingly fractionate and eventually liquefy tissue using high-intensity ultrasound pulses. Following cavitation generated by each pulse, coherent motion of the cavitation residual nuclei can be detected using metrics formed from ultrasound color Doppler acquisitions. In this paper, three experiments were performed to investigate the characteristics of this motion as real-time feedback on histotripsy tissue fractionation. In the first experiment, bubble-induced color Doppler (BCD) and particle image velocimetry (PIV) analysis monitored the residual cavitation nuclei in the treatment region in an agarose tissue phantom treated with two-cycle histotripsy pulses at [Formula: see text] using a 500-kHz transducer. Both BCD and PIV results showed brief chaotic motion of the residual nuclei followed by coherent motion first moving away from the transducer and then rebounding back. Velocity measurements from both PIV and BCD agreed well, showing a monotonic increase in rebound time up to a saturation point for increased therapy dose. In a second experiment, a thin layer of red blood cells (RBC) was added to the phantom to allow quantification of the fractionation of the RBC layer to compare with BCD metrics. A strong linear correlation was observed between the fractionation level and the time to BCD peak rebound velocity over histotripsy treatment. Finally, the correlation between BCD feedback and histotripsy tissue fractionation was validated in ex vivo porcine liver evaluated histologically. BCD metrics showed strong linear correlation with fractionation progression, suggesting that BCD provides useful quantitative real-time feedback on histotripsy treatment progression.

Bubble-induced Color Doppler Feedback for Histotripsy Tissue Fractionation

PubMed Central

Miller, Ryan M.; Zhang, Xi; Maxwell, Adam; Cain, Charles; Xu, Zhen

2016-01-01

Histotripsy therapy produces cavitating bubble clouds to increasingly fractionate and eventually liquefy tissue using high intensity ultrasound pulses. Following cavitation generated by each pulse, coherent motion of the cavitation residual nuclei can be detected using metrics formed from ultrasound color Doppler acquisitions. In this paper, three experiments were performed to investigate the characteristics of this motion as real-time feedback on histotripsy tissue fractionation. In the first experiment, bubble-induced color Doppler (BCD) and particle image velocimetry (PIV) analysis monitored the residual cavitation nuclei in the treatment region in an agarose tissue phantom treated with 2-cycle histotripsy pulses at > 30 MPa using a 500 kHz transducer. Both BCD and PIV results showed brief chaotic motion of the residual nuclei followed by coherent motion first moving away from the transducer and then rebounding back. Velocity measurements from both PIV and BCD agreed well, showing a monotonic increase in rebound time up to a saturation point for increased therapy dose. In a second experiment, a thin layer of red blood cells (RBC) was added to the phantom to allow quantification of the fractionation of the RBC layer to compare with BCD metrics. A strong linear correlation was observed between the fractionation level and the time to BCD peak rebound velocity over histotripsy treatment. Finally, the correlation between BCD feedback and histotripsy tissue fractionation was validated in ex vivo porcine liver evaluated histologically. BCD metrics showed strong linear correlation with fractionation progression, suggesting that BCD provides useful quantitative real-time feedback on histotripsy treatment progression. PMID:26863659

Flow in linearly sheared two-dimensional foams: From bubble to bulk scale.

PubMed

Katgert, Gijs; Latka, Andrzej; Möbius, Matthias E; van Hecke, Martin

2009-06-01

We probe the flow of two-dimensional (2D) foams, consisting of a monolayer of bubbles sandwiched between a liquid bath and glass plate, as a function of driving rate, packing fraction, and degree of disorder. First, we find that bidisperse, disordered foams exhibit strongly rate-dependent and inhomogeneous (shear-banded) velocity profiles, while monodisperse ordered foams are also shear banded but essentially rate independent. Second, we adapt a simple model [E. Janiaud, D. Weaire, and S. Hutzler, Phys. Rev. Lett. 97, 038302 (2006)] based on balancing the averaged drag forces between the bubbles and the top plate F[over ]_{bw} and the averaged bubble-bubble drag forces F[over ]_{bb} by assuming that F[over ]_{bw} approximately v;{2/3} and F[over ]_{bb} approximately ( partial differential_{y}v);{beta} , where v and ( partial differential_{y}v) denote average bubble velocities and gradients. This model captures the observed rate-dependent flows for beta approximately 0.36 , and the rate independent flows for beta approximately 0.67 . Third, we perform independent rheological measurements of F[over ]_{bw} and F[over ]_{bb} , both for ordered and disordered systems, and find these to be fully consistent with the forms assumed in the simple model. Disorder thus leads to a modified effective exponent beta . Fourth, we vary the packing fraction phi of the foam over a substantial range and find that the flow profiles become increasingly shear banded when the foam is made wetter. Surprisingly, the model describes flow profiles and rate dependence over the whole range of packing fractions with the same power-law exponents-only a dimensionless number k that measures the ratio of the prefactors of the viscous drag laws is seen to vary with packing fraction. We find that k approximately (phi-phi_{c});{-1} , where phi_{c} approximately 0.84 corresponds to the 2D jamming density, and suggest that this scaling follows from the geometry of the deformed facets between bubbles in contact. Overall, our work shows that the presence of disorder qualitatively changes the effective bubble-bubble drag forces and suggests a route to rationalize aspects of the ubiquitous Herschel-Bulkley (power-law) rheology observed in a wide range of disordered materials.

Foam flow in a model porous medium: II. The effect of trapped gas.

PubMed

Jones, S A; Getrouw, N; Vincent-Bonnieu, S

2018-05-09

Gas trapping is an important mechanism in both Water or Surfactant Alternating Gas (WAG/SAG) and foam injection processes in porous media. Foams for enhanced oil recovery (EOR) can increase sweep efficiency as they decrease the gas relative permeability, and this is mainly due to gas trapping. However, gas trapping mechanisms are poorly understood. Some studies have been performed during corefloods, but little work has been carried out to describe the bubble trapping behaviour at the pore scale. We have carried out foam flow tests in a micromodel etched with an irregular hexagonal pattern. Image analysis of the foam flow allowed the bubble centres to be tracked and local velocities to be obtained. It was found that the flow in the micromodel is dominated by intermittency and localized zones of trapped gas. The quantity of trapped gas was measured both by considering the fraction of bubbles that were trapped (via velocity thresholding) and by measuring the area fraction containing immobile gas (via image analysis). A decrease in the quantity of trapped gas was observed for both increasing total velocity and increasing foam quality. Calculations of the gas relative permeability were made with the Brooks Corey equation, using the measured trapped gas saturations. The results showed a decrease in gas relative permeabilities, and gas mobility, for increasing fractions of trapped gas. It is suggested that the shear thinning behaviour of foam could be coupled to the saturation of trapped gas.

Bubble Augmented Propulsor Mixture Flow Simulation near Choked Flow Condition

NASA Astrophysics Data System (ADS)

Choi, Jin-Keun; Hsiao, Chao-Tsung; Chahine, Georges

2013-03-01

The concept of waterjet thrust augmentation through bubble injection has been the subject of many patents and publications over the past several decades, and computational and experimental evidences of the augmentation of the jet thrust through bubble growth in the jet stream have been reported. Through our experimental studies, we have demonstrated net thrust augmentation as high as 70%for air volume fractions as high as 50%. However, in order to enable practical designs, an adequately validated modeling tool is required. In our previous numerical studies, we developed and validated a numerical code to simulate and predict the performance of a two-phase flow water jet propulsion system for low void fractions. In the present work, we extend the numerical method to handle higher void fractions to enable simulations for the high thrust augmentation conditions. At high void fractions, the speed of sound in the bubbly mixture decreases substantially and could be as low as 20 m/s, and the mixture velocity can approach the speed of sound in the medium. In this numerical study, we extend our numerical model, which is based on the two-way coupling between the mixture flow field and Lagrangian tracking of a large number of bubbles, to accommodate compressible flow regimes. Numerical methods used and the validation studies for various flow conditions in the bubble augmented propulsor will be presented. This work is supported by Office of Naval Research through contract N00014-11-C-0482 monitored by Dr. Ki-Han Kim.

Bubble Generation in a Continuous Liquid Flow Under Reduced Gravity Conditions

NASA Technical Reports Server (NTRS)

Pais, Salvatore Cezar

1999-01-01

The present work reports a study of bubble generation under reduced gravity conditions for both co-flow and cross-flow configurations. Experiments were performed aboard the DC-9 Reduced Gravity Aircraft at NASA Glenn Research Center, using an air-water system. Three different flow tube diameters were used: 1.27, 1.9, and 2.54 cm. Two different ratios of air injection nozzle to tube diameters were considered: 0.1 and 0.2. Gas and liquid volumetric flow rates were varied from 10 to 200 ml/s. It was experimentally observed that with increasing superficial liquid velocity, the bubbles generated decreased in size. The bubble diameter was shown to increase with increasing air injection nozzle diameters. As the tube diameter was increased, the size of the detached bubbles increased. Likewise, as the superficial liquid velocity was increased, the frequency of bubble formation increased and thus the time to detach forming bubbles decreased. Independent of the flow configuration (for either single nozzle or multiple nozzle gas injection), void fraction and hence flow regime transition can be controlled in a somewhat precise manner by solely varying the gas and liquid volumetric flow rates. On the other hand, it is observed that uniformity of bubble size can be controlled more accurately by using single nozzle gas injection than by using multiple port injection, since this latter system gives rise to unpredictable coalescence of adjacent bubbles. A theoretical model, based on an overall force balance, is employed to study single bubble generation in the dynamic and bubbly flow regime. Under conditions of reduced gravity, the gas momentum flux enhances bubble detachment; however, the surface tension forces at the nozzle tip inhibits bubble detachment. Liquid drag and inertia can act either as attaching or detaching force, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with performed experiments. However, at higher superficial,liquid velocities, the bubble neck length begins to significantly deviate from the value of the air injection nozzle diameter and thus the theory no longer predicts the experiment behavior. Effects of fluid properties, injection geometry and flow conditions on generated bubble size are investigated using the theoretical model. It is shown that bubble diameter is larger in a reduced gravity environment than in a normal gravity environment at similar flow condition and flow geometry.

Observations of bubbles in natural seep flares at MC 118 and GC 600 using in situ quantitative imaging

NASA Astrophysics Data System (ADS)

Wang, Binbin; Socolofsky, Scott A.; Breier, John A.; Seewald, Jeffrey S.

2016-04-01

This paper reports the results of quantitative imaging using a stereoscopic, high-speed camera system at two natural gas seep sites in the northern Gulf of Mexico during the Gulf Integrated Spill Research G07 cruise in July 2014. The cruise was conducted on the E/V Nautilus using the ROV Hercules for in situ observation of the seeps as surrogates for the behavior of hydrocarbon bubbles in subsea blowouts. The seeps originated between 890 and 1190 m depth in Mississippi Canyon block 118 and Green Canyon block 600. The imaging system provided qualitative assessment of bubble behavior (e.g., breakup and coalescence) and verified the formation of clathrate hydrate skins on all bubbles above 1.3 m altitude. Quantitative image analysis yielded the bubble size distributions, rise velocity, total gas flux, and void fraction, with most measurements conducted from the seafloor to an altitude of 200 m. Bubble size distributions fit well to lognormal distributions, with median bubble sizes between 3 and 4.5 mm. Measurements of rise velocity fluctuated between two ranges: fast-rising bubbles following helical-type trajectories and bubbles rising about 40% slower following a zig-zag pattern. Rise speed was uncorrelated with hydrate formation, and bubbles following both speeds were observed at both sites. Ship-mounted multibeam sonar provided the flare rise heights, which corresponded closely with the boundary of the hydrate stability zone for the measured gas compositions. The evolution of bubble size with height agreed well with mass transfer rates predicted by equations for dirty bubbles.

A Study of Bubble and Slug Gas-Liquid Flow in a Microgravity Environment

NASA Technical Reports Server (NTRS)

McQuillen, J.

2000-01-01

The influence of gravity on the two-phase flow dynamics is obvious.As the gravity level is reduced,there is a new balance between inertial and interfacial forces, altering the behavior of the flow. In bubbly flow,the absence of drift velocity leads to spherical-shaped bubbles with a rectilinear trajectory.Slug flow is a succession of long bubbles and liquid slug carrying a few bubbles. There is no flow reversal in the thin liquid film as the long bubble and liquid slug pass over the film. Although the flow structure seems to be simpler than in normal gravity conditions,the models developed for the prediction of flow behavior in normal gravity and extended to reduced gravity flow are unable to predict the flow behavior correctly.An additional benefit of conducting studies in microgravity flows is that these studies aide the development of understanding for normal gravity flow behavior by removing the effects of buoyancy on the shape of the interface and density driven shear flows between the gas and the liquid phases. The proposal calls to study specifically the following: 1) The dynamics of isolated bubbles in microgravity liquid flows will be analyzed: Both the dynamics of spherical isolated bubbles and their dispersion by turbulence, their interaction with the pipe wall,the behavior of the bubbles in accelerated or decelerated flows,and the dynamics of isolated cylindrical bubbles, their deformation in accelerated/decelerated flows (in converging or diverging channels), and bubble/bubble interaction. Experiments will consist of the use of Particle Image Velocimetry (PIV) and Laser Doppler Velocimeters (LDV) to study single spherical bubble and single and two cylindrical bubble behavior with respect to their influence on the turbulence of the surrounding liquid and on the wall 2) The dynamics of bubbly and slug flow in microgravity will be analyzed especially for the role of the coalescence in the transition from bubbly to slug flow (effect of fluid properties and surfactant), to identify clusters that promote coalescence and transition the void fraction distribution in bubbly and slug flow,to measure the wall friction in bubbly flow. These experiments will consist of multiple bubbles type flows and will utilize hot wire and film anemometers to measure liquid velocity and wall shear stress respectively and double fiber optic probes to measure bubble size and velocity as a function of tube radius and axial location.

Two-phase flow characteristics of liquid nitrogen in vertically upward 0.5 and 1.0 mm micro-tubes: Visualization studies

NASA Astrophysics Data System (ADS)

Zhang, P.; Fu, X.

2009-10-01

Application of liquid nitrogen to cooling is widely employed in many fields, such as cooling of the high temperature superconducting devices, cryosurgery and so on, in which liquid nitrogen is generally forced to flow inside very small passages to maintain good thermal performance and stability. In order to have a full understanding of the flow and heat transfer characteristics of liquid nitrogen in micro-tube, high-speed digital photography was employed to acquire the typical two-phase flow patterns of liquid nitrogen in vertically upward micro-tubes of 0.531 and 1.042 mm inner diameters. It was found from the experimental results that the flow patterns were mainly bubbly flow, slug flow, churn flow and annular flow. And the confined bubble flow, mist flow, bubble condensation and flow oscillation were also observed. These flow patterns were characterized in different types of flow regime maps. The surface tension force and the size of the diameter were revealed to be the major factors affecting the flow pattern transitions. It was found that the transition boundaries of the slug/churn flow and churn/annular flow of the present experiment shifted to lower superficial vapor velocity; while the transition boundary of the bubbly/slug flow shifted to higher superficial vapor velocity compared to the results of the room-temperature fluids in the tubes with the similar hydraulic diameters. The corresponding transition boundaries moved to lower superficial velocity when reducing the inner diameter of the micro-tubes. Time-averaged void fraction and heat transfer characteristics for individual flow patterns were presented and special attention was paid to the effect of the diameter on the variation of void fraction.

Flowfield measurements in a separated and reattached flat plate turbulent boundary layer

NASA Technical Reports Server (NTRS)

Patrick, William P.

1987-01-01

The separation and reattachment of a large-scale, two-dimensional turbulent boundary layer at low subsonic speed on a flat plate has been studied experimentally. The separation bubble was 55 cm long and had a maximum bubble thickness, measured to the height of the mean dividing streamline, of 17 cm, which was twice the thickness of the inlet boundary layer. A combination of laser velocimetry, hot-wire anemometry, pneumatic probing techniques, and flow visualization were used as diagnostics. Principal findings were that an outer inviscid rotational flow was defined which essentially convected over the blockage associated with the inner, viscously dominated bubble recirculation region. A strong backflow region in which the flow moved upstream 100 percent of the time was measured near the test surface over the central 35 percent of the bubble. A laminar backflow boundary layer having pseudo-turbulent characteristics including a log-linear velocity profile was generated under the highly turbulent backflow. Velocity profile shapes in the reversed flow region matched a previously developed universal backflow profile at the upstream edge of the separation region but not in the steady backflow region downstream. A smoke flow visualization movie and hot-film measurements revealed low frequency nonperiodic flapping at reattachment. However, forward flow fraction data at reattachment and mean velocity profiles in the redeveloping boundary layer downstream of reattachment correlated with backward-facing step data when the axial dimension was scaled by the distance from the maximum bubble thickness to reattachment.

Factors That Modulate Properties of Primary Marine Aerosol Generated From Ambient Seawater on Ships at Sea

NASA Astrophysics Data System (ADS)

Keene, William C.; Long, Michael S.; Reid, Jeffrey S.; Frossard, Amanda A.; Kieber, David J.; Maben, John R.; Russell, Lynn M.; Kinsey, Joanna D.; Quinn, Patricia K.; Bates, Timothy S.

2017-11-01

Model primary marine aerosol (mPMA) was produced by bubbling clean air through flowing natural seawater in a high-capacity generator deployed on ships in the eastern North Pacific and western North Atlantic Oceans. Physicochemical properties of seawater and mPMA were quantified to characterize factors that modulated production. Differences in surfactant organic matter (OM) and associated properties including surface tension sustained plumes with smaller bubble sizes, slower rise velocities, larger void fractions, and older surface ages in biologically productive relative to oligotrophic seawater. Production efficiencies for mPMA number (PEnum) and mass (PEmass) per unit air detrained from biologically productive seawater during daytime were greater and mass median diameters smaller than those in the same seawater at night and in oligotrophic seawater during day and night. PEmass decreased with increasing air detrainment rate suggesting that surface bubble rafts suppressed emission of jet droplets and associated mPMA mass. Relative to bubbles emitted at 60 cm depth, PEnum for bubbles emitted from 100 cm depth was approximately 2 times greater. mPMA OM enrichment factors (EFs) and mass fractions based on a coarse frit, fine frits, and a seawater jet exhibited similar size-dependent variability over a wide range in chlorophyll a concentrations. Results indicate that the physical production of PMA number and mass from the ocean surface varies systematically as interrelated functions of seawater type and, in biologically productive waters, time of day; bubble injection rate, depth, size, and surface age; and physical characteristics of the air-water interface whereas size-resolved OM EFs and mass fractions are relatively invariant.

The sudden coalescene model of the boiling crisis

DOE Office of Scientific and Technical Information (OSTI.GOV)

Carrica, P.M.; Clausse, A.

1995-09-01

A local two-phase flow integral model of nucleate boiling and crisis is presented. The model is based on average balances on a control volume, yielding to a set of three nonlinear differential equations for the local void fraction, bubble number density and velocity. Boiling crisis as critical heat flux is interpreted as a dynamic transition caused by the coalescence of bubbles near the heater. The theoretical dynamic model is compared with experimental results obtained for linear power ramps in a horizontal plate heater in R-113, showing an excellent qualitative agreement.

Possible high sonic velocity due to the inclusion of gas bubbles in water

NASA Astrophysics Data System (ADS)

Banno, T.; Mikada, H.; Goto, T.; Takekawa, J.

2010-12-01

If formation water becomes multi-phase by inclusion of gas bubbles, sonic velocities would be strongly influenced. In general, sonic velocities are knocked down due to low bulk moduli of the gas bubbles. However, sonic velocities may increase depending on the size of gas bubbles, when the bubbles in water or other media oscillate due to incoming sonic waves. Sonic waves are scattered by the bubbles and the superposition of the incoming and the scattered waves result in resonant-frequency-dependent behavior. The phase velocity of sonic waves propagating in fluids containing bubbles, therefore, probably depends on their frequencies. This is a typical phenomenon called “wave dispersion.” So far we have studied about the bubble impact on sonic velocity in bubbly media, such as the formation that contains gas bubbles. As a result, it is shown that the bubble resonance effect is a key to analyze the sonic phase velocity increase. Therefore to evaluate the resonance frequency of bubbles is important to solve the frequency response of sonic velocity in formations having bubbly fluids. There are several analytical solutions of the resonance frequency of bubbles in water. Takahira et al. (1994) derived a equation that gives us the resonance frequency considering bubble - bubble interactions. We have used this theory to calculate resonance frequency of bubbles at the previous work. However, the analytical solution of the Takahira’s equation is based on several assumptions. Therefore we used a numerical approach to calculate the bubble resonance effect more precisely in the present study. We used the boundary element method (BEM) to reproduce a bubble oscillation in incompressible liquid. There are several reasons to apply the BEM. Firstly, it arrows us to model arbitrarily sets and shapes of bubbles. Secondly, it is easy to use the BEM to reproduce a boundary-surface between liquid and gas. The velocity potential of liquid surrounding a bubble satisfies the Laplace equation when the liquid is supposed to be incompressible. We got the boundary integral equation from the Laplace equation and solved the boundary integral equation by the BEM. Then, we got the gradient of the velocity potential from the BEM. We used this gradient to get time derivative of the velocity potential from the Bernouii’s equation. And we used the second order Adams-Bashforth method to execute time integration of the velocity potential. We conducted this scheme iteratively to calculate a bubble oscillation. At each time step, we input a pressure change as a sinusoidal wave. As a result, we observed a bubble oscillation following the pressure frequency. We also evaluated the resonance frequency of a bubble by changing the pressure frequency. It showed a good agreement with the analytical solution described above. Our future work is to extend the calculation into plural bubbles condition. We expect that interaction between bubbles becomes strong and resonance frequency of bubbles becomes small when distance between bubbles becomes small.

Slug Flow Analysis in Vertical Large Diameter Pipes

NASA Astrophysics Data System (ADS)

Roullier, David

The existence of slug flow in vertical co-current two-phase flow is studied experimentally and theoretically. The existence of slug flow in vertical direction implies the presence of Taylor bubbles separated by hydraulically sealed liquid slugs. Previous experimental studies such as Ombere-Ayari and Azzopardi (2007) showed the evidence of the non-existence of Taylor bubbles for extensive experimental conditions. Models developed to predict experimental behavior [Kocamustafaogullari et al. (1984), Jayanti and Hewitt. (1990) and Kjoolas et al. (2017)] suggest that Taylor bubbles may disappear at large diameters and high velocities. A 73-ft tall and 101.6-mm internal diameter test facility was used to conduct the experiments allowing holdup and pressure drop measurements at large L/D. Superficial liquid and gas velocities varied from 0.05-m/s to 0.2 m/s and 0.07 m/s to 7.5 m/s, respectively. Test section pressure varied from 38 psia to 84 psia. Gas compressibility effect was greatly reduced at 84 psia. The experimental program allowed to observe the flow patterns for flowing conditions near critical conditions predicted by previous models (air-water, 1016 mm ID, low mixture velocities). Flow patterns were observed in detail using wire-mesh sensor measurements. Slug-flow was observed for a narrow range of experimental conditions at low velocities. Churn-slug and churn-annular flows were observed for most of the experimental data-points. Cap-bubble flow was observed instead of bubbly flow at low vSg. Wire-mesh measurements showed that the liquid has a tendency to remain near to the walls. The standard deviation of radial holdup profile correlates to the flow pattern observed. For churn-slug flow, the profile is convex with a single maximum near the pipe center while it exhibits a concave shape with two symmetric maxima close to the wall for churn-annular flow. The translational velocity was measured by two consecutive wire-mesh sensor crosscorrelation. The results show linear trends at low mixture velocities and non-linear behaviors at high mixture velocities. The translational velocity trends seem to be related to the flow-pattern observed, namely to the ability of the gas to flow through the liquid structures. A simplified Taylor bubble stability model is proposed. The model allows to estimate under which conditions Taylor bubbles disappear, properly accounting for the diameter effect and velocity effect observed experimentally. In addition, annular flow distribution coefficient relating true holdup to centerline holdup in vertical flow is proposed. The proposed coefficient defines the tendency of the liquid to remain near the walls. This coefficient increases linearly with the void fraction.

Do volcanic gases represent equilibrium volatile concentrations? Some insights from a model of diffusive fractionation during rapid bubble growth

NASA Astrophysics Data System (ADS)

Baker, D. R.

2012-12-01

Measurements of volcanic gas compositions are often presumed to be directly related to equilibrium compositions of fluids exsolved at depth in magmatic systems that rapidly escape into the atmosphere. In particular, changes in the ratios of volatile species concentrations in volcanic gases have been interpreted to reflect influx of new magma batches or changes in the degassing depth. However, other mechanisms can also yield changes in volcanic gas compositions. One such mechanism is diffusive fractionation during rapid bubble growth. Such fractionation can occur because radial growth rates of bubbles in magmas are estimated to be in the range of 10-6 to 10-3 m s-1 and diffusion coefficients of minor volatiles (e.g., Cl, F, S, CO2) are orders of magnitude slower, 10-12 to 10-9 m2 s-1. Thus a bubble that rapidly grows and subsequently loses its volatiles to the surface may contribute a fluid sample whose concentration is affected by the interplay between the kinetics of bubble growth and volatile diffusion in the melt. A finite difference code was developed to calculate the effects of rapid bubble growth on the concentration of minor elements in the bubble for a spherical growth geometry. The bubble is modeled with a fixed growth rate and a constant equilibrium fluid-melt partition coefficient, KD. Bubbles were modeled to grow to a radius of 50 μm, the size at which the dominant bubble growth mechanism appears to change from diffusion to coalescence. The critical variables that control the departure from equilibrium behavior are the K D and the ratio of the growth velocity, V, to the diffusivity, D. Modeling bubble growth in a magma chamber at 100 MPa demonstrates that when KD is in the range of 10 to 1000 at low V/D values (e.g., 103 m-1) the composition of the fluid is at, or near, equilibrium with the melt. However, as V/D increases the bubble composition deviates increasingly from equilibrium. For V/D ratios of 105 and equilibrium KD's of either 50 or 100 (similar to estimates for S), a bubble with a 50 μm radius will contain a fluid whose concentration was apparently determined by a KD of less than 10. These models also demonstrate that the combination of rapid bubble growth with slow diffusion can deplete the melt in the volatile species only within the immediate neighborhood, on the order of 100 μm. If bubbles are spaced further apart the melts may retain significant concentrations of dissolved volatiles, which could lead to secondary and tertiary nucleation events. These models for diffusive fractionation during rapid bubble growth suggest that changes in the ratios of minor elements in volcanic gases may be influenced by bubble growth rate changes. Volatiles with lower diffusivities and volatiles with very high or very low partition coefficients will be more influenced by this process. Diffusive fractionation may be responsible for the drop in the CO2/SO2 ratios sometimes observed prior to large eruptions of Stromboli volcano.

Application of Electromagnetic Induction Technique to Measure the Void Fraction in Oil/Gas Two Phase Flow

NASA Astrophysics Data System (ADS)

Wahhab, H. A. Abdul; Aziz, A. R. A.; Al-Kayiem, H. H.; Nasif, M. S.; Reda, M. N.

2018-03-01

In this work, electromagnetic induction technique of measuring void fraction in liquid/gas fuel flow was utilized. In order to improve the electric properties of liquid fuel, an iron oxide Fe3O4 nanoparticles at 3% was blended to enhance the liquid fuel magnetization. Experiments have been conducted for a wide range of liquid and gas superficial velocities. From the experimental results, it was realized that there is an existing linear relationship between the void fraction and the measured electromotive force, when induction coils were connected in series for excitation coils, regardless of increase or decrease CNG bubbles distribution in liquid fuel flow. Therefore, it was revealed that the utilized method yielded quite reasonable account for measuring the void fraction, showing good agreement with the other available measurement techniques in the two-phase flow, and also with the published literature of the bubbly flow pattern. From the results of the present investigation, it has been proven that the electromagnetic induction is a feasible technique for the actual measurement of void fraction in a Diesel/CNG fuel flow.

Studies of Two-Phase Gas-Liquid Flow in Microgravity. Ph.D. Thesis, Dec. 1994

NASA Technical Reports Server (NTRS)

Bousman, William Scott

1995-01-01

Two-phase gas-liquid flows are expected to occur in many future space operations. Due to a lack of buoyancy in the microgravity environment, two-phase flows are known to behave differently than those in earth gravity. Despite these concerns, little research has been conducted on microgravity two-phase flow and the current understanding is poor. This dissertation describes an experimental and modeling study of the characteristics of two-phase flows in microgravity. An experiment was operated onboard NASA aircraft capable of producing short periods of microgravity. In addition to high speed photographs of the flows, electronic measurements of void fraction, liquid film thickness, bubble and wave velocity, pressure drop and wall shear stress were made for a wide range of liquid and gas flow rates. The effects of liquid viscosity, surface tension and tube diameter on the behavior of these flows were also assessed. From the data collected, maps showing the occurrence of various flow patterns as a function of gas and liquid flow rates were constructed. Earth gravity two-phase flow models were compared to the results of the microgravity experiments and in some cases modified. Models were developed to predict the transitions on the flow pattern maps. Three flow patterns, bubble, slug and annular flow, were observed in microgravity. These patterns were found to occur in distinct regions of the gas-liquid flow rate parameter space. The effect of liquid viscosity, surface tension and tube diameter on the location of the boundaries of these regions was small. Void fraction and Weber number transition criteria both produced reasonable transition models. Void fraction and bubble velocity for bubble and slug flows were found to be well described by the Drift-Flux model used to describe such flows in earth gravity. Pressure drop modeling by the homogeneous flow model was inconclusive for bubble and slug flows. Annular flows were found to be complex systems of ring-like waves and a substrate film. Pressure drop was best fitted with the Lockhart- Martinelli model. Force balances suggest that droplet entrainment may be a large component of the total pressure drop.

MOBI: Microgravity Observations of Bubble Interactions

NASA Technical Reports Server (NTRS)

Koch, Donald L.; Sangani, Ashok

2004-01-01

One of the greatest uncertainties affecting the design of multiphase flow technologies for space exploration is the spatial distribution of phases that will arise in microgravity or reduced gravity. On Earth, buoyancy-driven motion predominates whereas the shearing of the bubble suspension controls its behavior in microgravity. We are conducting a series of ground-based experiments and a flight experiment spanning the full range of ratios of buoyancy to shear. These include: (1) bubbles rising in a quiescent liquid in a vertical channel; (2) weak shear flow induced by slightly inclining the channel; (3) moderate shear flow in a terrestrial vertical pipe flow; and (4) shearing of a bubble suspension in a cylindrical Couette cell in microgravity. We consider nearly monodisperse suspensions of 1 to 1.8 mm diameter bubbles in aqueous electrolyte solutions. The liquid velocity disturbance produced by bubbles in this size range can often be described using an inviscid analysis. Electrolytic solutions lead to hydrophilic repulsion forces that stabilize the bubble suspension without causing Marangoni stresses. We will discuss the mechanisms that control the flow behavior and phase distribution in the ground-based experiments and speculate on the factors that may influence the suspension flow and bubble volume fraction distribution in the flight experiment.

Experimental study on the void fraction of air-water two-phase flow in a horizontal circular minichannel

NASA Astrophysics Data System (ADS)

Sudarja, Indarto, Deendarlianto, Haq, Aqli

2016-06-01

Void fraction is an important parameter in two-phase flow. In the present work, the adiabatic two-phase air-water flow void fraction in a horizontal minichannel has been studied experimentally. A transparent circular channel with 1.6 mm inner diameter was employed as the test section. Superficial gas and liquid velocities were varied in the range of 1.25 - 66.3 m/s and 0.033 - 4.935 m/s, respectively. Void fraction data were obtained by analyzing the flow images being captured by using a high-speed camera. Here, the homogeneous (β) and the measured void fractions (ɛ), respectively, were compared to the existing correlations. It was found that: (1) for the bubbly and slug flows, the void fractions increases with the increase of JG, (2) for churn, slug-annular, and annular flow patterns, there is no specific correlation between JG and void fraction was observed due to effect of the slip between gas and liquid, and (3) whilst for bubbly and slug flows the void fractions are close to homogeneous line, for churn, annular, and slug-annular flows are far below the homogeneous line. It indicates that the slip ratios for the second group of flow patterns are higher than unity.

A Study of Nucleate Boiling with Forced Convection in Microgravity

NASA Technical Reports Server (NTRS)

Merte, Herman, Jr.

1999-01-01

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Blast wave attenuation in liquid foams: role of gas and evidence of an optimal bubble size.

PubMed

Monloubou, Martin; Bruning, Myrthe A; Saint-Jalmes, Arnaud; Dollet, Benjamin; Cantat, Isabelle

2016-09-28

Liquid foams are excellent systems to mitigate pressure waves such as acoustic or blast waves. The understanding of the underlying dissipation mechanisms however still remains an active matter of debate. In this paper, we investigate the attenuation of a weak blast wave by a liquid foam. The wave is produced with a shock tube and impacts a foam, with a cylindrical geometry. We measure the wave attenuation and velocity in the foam as a function of bubble size, liquid fraction, and the nature of the gas. We show that the attenuation depends on the nature of the gas and we experimentally evidence a maximum of dissipation for a given bubble size. All features are qualitatively captured by a model based on thermal dissipation in the gas.

Hindrance Velocity Model for Phase Segregation in Suspensions of Poly-dispersed Randomly Oriented Spheroids

NASA Astrophysics Data System (ADS)

Faroughi, S. A.; Huber, C.

2015-12-01

Crystal settling and bubbles migration in magmas have significant effects on the physical and chemical evolution of magmas. The rate of phase segregation is controlled by the force balance that governs the migration of particles suspended in the melt. The relative velocity of a single particle or bubble in a quiescent infinite fluid (melt) is well characterized; however, the interplay between particles or bubbles in suspensions and emulsions and its effect on their settling/rising velocity remains poorly quantified. We propose a theoretical model for the hindered velocity of non-Brownian emulsions of nondeformable droplets, and suspensions of spherical solid particles in the creeping flow regime. The model is based on three sets of hydrodynamic corrections: two on the drag coefficient experienced by each particle to account for both return flow and Smoluchowski effects and a correction on the mixture rheology to account for nonlocal interactions between particles. The model is then extended for mono-disperse non-spherical solid particles that are randomly oriented. The non-spherical particles are idealized as spheroids and characterized by their aspect ratio. The poly-disperse nature of natural suspensions is then taken into consideration by introducing an effective volume fraction of particles for each class of mono-disperse particles sizes. Our model is tested against new and published experimental data over a wide range of particle volume fraction and viscosity ratios between the constituents of dispersions. We find an excellent agreement between our model and experiments. We also show two significant applications for our model: (1) We demonstrate that hindered settling can increase mineral residence time by up to an order of magnitude in convecting magma chambers. (2) We provide a model to correct for particle interactions in the conventional hydrometer test to estimate the particle size distribution in soils. Our model offers a greatly improved agreement with the results obtained with direct measurement methods such as laser diffraction.

On shapes and motion of an elongated bubble in downward liquid pipe flow

NASA Astrophysics Data System (ADS)

Fershtman, A.; Babin, V.; Barnea, D.; Shemer, L.

2017-11-01

In stagnant liquid, or in a steady upward liquid pipe flow, an elongated (Taylor) bubble has a symmetric shape. The translational velocity of the bubble is determined by buoyancy and the liquid velocity profile ahead of it. In downward flow, however, the symmetry of the bubble nose can be lost. Taylor bubble motion in downward flow is important in numerous applications such as chemical plants and cooling systems that often contain countercurrent gas-liquid flow. In the present study, the relation between the Taylor bubble shape and its translational velocity is investigated experimentally in a vertical pipe for various downward liquid flow rates. At higher downward velocities, the bubble may be forced by the background flow to propagate downward against buoyancy. In order to include those cases as well in our experimental analysis, the bubbles were initially injected into stagnant liquid, whereas the downward flow was initiated at a later stage. This experimental procedure allowed us to identify three distinct modes of translational velocities for a given downward background liquid flow; each velocity corresponds to a different bubble shape. Hydrodynamic mechanisms that govern the transition between the modes observed in the present study are discussed.

A new method for ultrasound detection of interfacial position in gas-liquid two-phase flow.

PubMed

Coutinho, Fábio Rizental; Ofuchi, César Yutaka; de Arruda, Lúcia Valéria Ramos; Neves, Flávio; Morales, Rigoberto E M

2014-05-22

Ultrasonic measurement techniques for velocity estimation are currently widely used in fluid flow studies and applications. An accurate determination of interfacial position in gas-liquid two-phase flows is still an open problem. The quality of this information directly reflects on the accuracy of void fraction measurement, and it provides a means of discriminating velocity information of both phases. The algorithm known as Velocity Matched Spectrum (VM Spectrum) is a velocity estimator that stands out from other methods by returning a spectrum of velocities for each interrogated volume sample. Interface detection of free-rising bubbles in quiescent liquid presents some difficulties for interface detection due to abrupt changes in interface inclination. In this work a method based on velocity spectrum curve shape is used to generate a spatial-temporal mapping, which, after spatial filtering, yields an accurate contour of the air-water interface. It is shown that the proposed technique yields a RMS error between 1.71 and 3.39 and a probability of detection failure and false detection between 0.89% and 11.9% in determining the spatial-temporal gas-liquid interface position in the flow of free rising bubbles in stagnant liquid. This result is valid for both free path and with transducer emitting through a metallic plate or a Plexiglas pipe.

A New Method for Ultrasound Detection of Interfacial Position in Gas-Liquid Two-Phase Flow

PubMed Central

Coutinho, Fábio Rizental; Ofuchi, César Yutaka; de Arruda, Lúcia Valéria Ramos; Jr., Flávio Neves; Morales, Rigoberto E. M.

2014-01-01

Ultrasonic measurement techniques for velocity estimation are currently widely used in fluid flow studies and applications. An accurate determination of interfacial position in gas-liquid two-phase flows is still an open problem. The quality of this information directly reflects on the accuracy of void fraction measurement, and it provides a means of discriminating velocity information of both phases. The algorithm known as Velocity Matched Spectrum (VM Spectrum) is a velocity estimator that stands out from other methods by returning a spectrum of velocities for each interrogated volume sample. Interface detection of free-rising bubbles in quiescent liquid presents some difficulties for interface detection due to abrupt changes in interface inclination. In this work a method based on velocity spectrum curve shape is used to generate a spatial-temporal mapping, which, after spatial filtering, yields an accurate contour of the air-water interface. It is shown that the proposed technique yields a RMS error between 1.71 and 3.39 and a probability of detection failure and false detection between 0.89% and 11.9% in determining the spatial-temporal gas-liquid interface position in the flow of free rising bubbles in stagnant liquid. This result is valid for both free path and with transducer emitting through a metallic plate or a Plexiglas pipe. PMID:24858961

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. I - Velocity selection. II - Stability

NASA Technical Reports Server (NTRS)

Tanveer, S.

1989-01-01

An asymptotic theory is presented for the determination of velocity and linear stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. First the bubble velocity relative to the fluid velocity at infinity is determined for small surface tension by means of a transcendentally small correction to the asymptotic series solution. In addition, a linear stability analysis shows that only the solution branch corresponding to the largest possible bubble velocity for given surface tension is stable, while all the others are unstable.

Motion-sensitized SPRITE measurements of hydrodynamic cavitation in fast pipe flow.

PubMed

Adair, Alexander; Mastikhin, Igor V; Newling, Benedict

2018-06-01

The pressure variations experienced by a liquid flowing through a pipe constriction can, in some cases, result in the formation of a bubble cloud (i.e., hydrodynamic cavitation). Due to the nature of the bubble cloud, it is ideally measured through the use of non-optical and non-invasive techniques; therefore, it is well-suited for study by magnetic resonance imaging. This paper demonstrates the use of Conical SPRITE (a 3D, centric-scan, pure phase-encoding pulse sequence) to acquire time-averaged void fraction and velocity information about hydrodynamic cavitation for water flowing through a pipe constriction. Copyright © 2018 Elsevier Inc. All rights reserved.

Geomagnetically conjugate observation of plasma bubbles and thermospheric neutral winds at low latitudes

NASA Astrophysics Data System (ADS)

Fukushima, D.; Shiokawa, K.; Otsuka, Y.; Nishioka, M.; Kubota, M.; Tsugawa, T.; Nagatsuma, T.; Komonjinda, S.; Yatini, C. Y.

2015-03-01

This is the first paper that reports simultaneous observations of zonal drift of plasma bubbles and the thermospheric neutral winds at geomagnetically conjugate points in both hemispheres. The plasma bubbles were observed in the 630 nm nighttime airglow images taken by using highly sensitive all-sky airglow imagers at Kototabang, Indonesia (geomagnetic latitude (MLAT): 10.0°S), and Chiang Mai, Thailand (MLAT: 8.9°N), which are nearly geomagnetically conjugate stations, for 7 h from 13 to 20 UT (from 20 to 03 LT) on 5 April 2011. The bubbles continuously propagated eastward with velocities of 100-125 m/s. The 630 nm images at Chiang Mai and those mapped to the conjugate point of Kototabang fit very well, which indicates that the observed plasma bubbles were geomagnetically connected. The eastward thermospheric neutral winds measured by two Fabry-Perot interferometers were 70-130 m/s at Kototabang and 50-90 m/s at Chiang Mai. We compared the observed plasma bubble drift velocity with the velocity calculated from the observed neutral winds and the model conductivity, to investigate the F region dynamo contribution to the bubble drift velocity. The estimated drift velocities were 60-90% of the observed velocities of the plasma bubbles, suggesting that most of the plasma bubble velocity can be explained by the F region dynamo effect.

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. Part 1: Velocity selection

NASA Technical Reports Server (NTRS)

Tanveer, Saleh

1989-01-01

An asymptotic theory is presented for the determination of velocity and linear stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. In the first part, the bubble velocity U relative to the fluid velocity at infinity is determined for small surface tension T by determining transcendentally small correction to the asymptotic series solution. It is found that for any relative bubble velocity U in the interval (U(c),2), solutions exist at a countably infinite set of values of T (which has zero as its limit point) corresponding to the different branches of bubble solutions. U(c) decreases monotonically from 2 to 1 as the bubble area increases from 0 to infinity. However, for a bubble of arbitrarily given size, as T approaches 0, solution exists on any given branch with relative bubble velocity U satisfying the relation 2-U = cT to the 2/3 power, where c depends on the branch but is independent of the bubble area. The analytical evidence further suggests that there are no solutions for U greater than 2. These results are in agreement with earlier analytical results for a finger. In Part 2, an analytic theory is presented for the determination of the linear stability of the bubble in the limit of zero surface tension. Only the solution branch corresponding to the largest possible U for given surface tension is found to be stable, while all the others are unstable, in accordance with earlier numerical results.

Numerical Analysis of the Influence of Low Frequency Vibration on Bubble Growth

PubMed Central

Han, D.; Kedzierski, Mark A.

2017-01-01

Numerical simulation of bubble growth during pool boiling under the influence of low frequency vibration was performed to understand the influence of common vibrations such as those induced by wind, highway transportation, and nearby mechanical devices on the performance of thermal systems that rely on boiling. The simulations were done for saturated R123 boiling at 277.6 K with a 15 K wall superheat. The numerical volume-of-fluid method (fixed grid) was used to define the liquid-vapor interface. The basic bubble growth characteristics including the bubble departure diameter and the bubble departure time were determined as a function of the bubble contact angle (20°–80°), the vibration displacement (10 µm–50 µm), the vibration frequency (5 Hz–25 Hz), and the initial vibration direction (positive or negative). The bubble parameters were shown to be strongly dependent on the bubble contact angle at the surface. For example, both the bubble departure diameter and the bubble departure time increased with the contact angle. At the same vibration frequency and the initial vibration direction, the bubble departure diameter and the bubble departure time both decreased with increasing vibration displacement. In addition, the vibration frequency had a greater effect on the bubble growth characteristics than did the vibration displacement. The vibration frequency effect was strongly influenced by the initial vibration direction. The pressure contour, the volume fraction of vapor phase, the temperature profile, and the velocity vector were investigated to understand these dynamic bubble behaviors. The limitation of the computational fluid dynamics approach was also described. PMID:28747812

Stream-wise distribution of skin-friction drag reduction on a flat plate with bubble injection

NASA Astrophysics Data System (ADS)

Qin, Shijie; Chu, Ning; Yao, Yan; Liu, Jingting; Huang, Bin; Wu, Dazhuan

2017-03-01

To investigate the stream-wise distribution of skin-friction drag reduction on a flat plate with bubble injection, both experiments and simulations of bubble drag reduction (BDR) have been conducted in this paper. Drag reductions at various flow speeds and air injection rates have been tested in cavitation tunnel experiments. Visualization of bubble flow pattern is implemented synchronously. The computational fluid dynamics (CFD) method, in the framework of Eulerian-Eulerian two fluid modeling, coupled with population balance model (PBM) is used to simulate the bubbly flow along the flat plate. A wide range of bubble sizes considering bubble breakup and coalescence is modeled based on experimental bubble distribution images. Drag and lift forces are fully modeled based on applicable closure models. Both predicted drag reductions and bubble distributions are in reasonable concordance with experimental results. Stream-wise distribution of BDR is revealed based on CFD-PBM numerical results. In particular, four distinct regions with different BDR characteristics are first identified and discussed in this study. Thresholds between regions are extracted and discussed. And it is highly necessary to fully understand the stream-wise distribution of BDR in order to establish a universal scaling law. Moreover, mechanism of stream-wise distribution of BDR is analysed based on the near-wall flow parameters. The local drag reduction is a direct result of near-wall max void fraction. And the near-wall velocity gradient modified by the presence of bubbles is considered as another important factor for bubble drag reduction.

Three-dimensional single-mode nonlinear ablative Rayleigh-Taylor instability

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yan, R.; Aluie, H.; Laboratory for Laser Energetics, University of Rochester, Rochester, New York 14627

The nonlinear evolution of the single-mode ablative Rayleigh-Taylor instability is studied in three dimensions. As the mode wavelength approaches the cutoff of the linear spectrum (short-wavelength modes), it is found that the three-dimensional (3D) terminal bubble velocity greatly exceeds both the two-dimensional (2D) value and the classical 3D bubble velocity. Unlike in 2D, the 3D short-wavelength bubble velocity does not saturate. The growing 3D bubble acceleration is driven by the unbounded accumulation of vorticity inside the bubble. The vorticity is transferred by mass ablation from the Rayleigh-Taylor spikes to the ablated plasma filling the bubble volume.

Effect of disjoining pressure on terminal velocity of a bubble sliding along an inclined wall.

PubMed

Del Castillo, Lorena A; Ohnishi, Satomi; White, Lee R; Carnie, Steven L; Horn, Roger G

2011-12-15

The influence of salt concentration on the terminal velocities of gravity-driven single bubbles sliding along an inclined glass wall has been investigated, in an effort to establish whether surface forces acting between the wall and the bubble influence the latter's mobility. A simple sliding bubble apparatus was employed to measure the terminal velocities of air bubbles with radii ranging from 0.3 to 1.5 mm sliding along the interior wall of an inclined Pyrex glass cylinder with inclination angles between 0.6 and 40.1°. Experiments were performed in pure water, 10 mM and 100 mM KCl solutions. We compared our experimental results with a theory by Hodges et al. which considers hydrodynamic forces only, and with a theory developed by two of us which considers surface forces to play a significant role. Our experimental results demonstrate that the terminal velocity of the bubble not only varies with the angle of inclination and the bubble size but also with the salt concentration, particularly at low inclination angles of ∼1-5°, indicating that double-layer forces between the bubble and the wall influence the sliding behavior. This is the first demonstration that terminal velocities of sliding bubbles are affected by disjoining pressure. Copyright © 2011 Elsevier Inc. All rights reserved.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, B.; The Peac Institute of Multiscale Sciences, Chengdu, Sichuan 610207; Wang, L.

With large-scale molecular dynamics simulations, we investigate shock response of He nanobubbles in single crystal Cu. For sufficient bubble size or internal pressure, a prismatic dislocation loop may form around a bubble in unshocked Cu. The internal He pressure helps to stabilize the bubble against plastic deformation. However, the prismatic dislocation loops may partially heal but facilitate nucleation of new shear and prismatic dislocation loops. For strong shocks, the internal pressure also impedes internal jetting, while a bubble assists local melting; a high speed jet breaks a He bubble into pieces dispersed among Cu. Near-surface He bubbles may burst andmore » form high velocity ejecta containing atoms and small fragments, while the ejecta velocities do not follow the three-dimensional Maxwell-Boltzmann distributions expected for thermal equilibrium. The biggest fragment size deceases with increasing shock strength. With a decrease in ligament thickness or an increase in He bubble size, the critical shock strength required for bubble bursting decreases, while the velocity range, space extension and average velocity component along the shock direction, increase. Small bubbles are more efficient in mass ejecting. Compared to voids and perfect single crystal Cu, He bubbles have pronounced effects on shock response including bubble/void collapse, Hugoniot elastic limit (HEL), deformation mechanisms, and surface jetting. HEL is the highest for perfect single crystal Cu with the same orientations, followed by He bubbles without pre-existing prismatic dislocation loops, and then voids. Complete void collapse and shear dislocations occur for embedded voids, as opposed to partial collapse, and shear and possibly prismatic dislocations for He bubbles. He bubbles lower the threshhold shock strength for ejecta formation, and increase ejecta velocity and ejected mass.« less

Trends in long-period seismicity related to magmatic fluid compositions

USGS Publications Warehouse

Morrissey, M.M.; Chouet, B.A.

2001-01-01

Sound speeds and densities are calculated for three different types of fluids: gas-gas mixture; ash-gas mixture; and bubbly liquid. These fluid properties are used to calculate the impedance contrast (Z) and crack stiffness (C) in the fluid-driven crack model (Chouet: J. Geophys. Res., 91 (1986) 13,967; 101 (1988) 4375; A seismic model for the source of long-period events and harmonic tremor. In: Gasparini, P., Scarpa, R., Aki, K. (Eds.), Volcanic Seismology, IAVCEI Proceedings in Volcanology, Springer, Berlin, 3133). The fluid-driven crack model describes the far-field spectra of long-period (LP) events as modes of resonance of the crack. Results from our calculations demonstrate that ash-laden gas mixtures have fluid to solid density ratios comparable to, and fluid to solid velocity ratios lower than bubbly liquids (gas-volume fractions 20% gas-volume fraction yields values of Q-1r similar to those for a rectangular crack. As with gas-gas and ash-gas mixtures, an increase in mass fraction narrows the bandwidth of the dominant mode and shifts the spectra to lower frequencies. Including energy losses due to dissipative processes in a bubbly liquid increases attenuation. Attenuation may also be higher in ash-gas mixtures and foams if the effects of momentum and mass transfer between the phases were considered in the calculations. ?? 2001 Elsevier Science B. V. All rights reserved.

Three-dimensional single-mode nonlinear ablative Rayleigh-Taylor instability

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yan, R.; Betti, R.; Sanz, J.

The nonlinear evolution of the single-mode ablative Rayleigh-Taylor instability is studied in three dimensions. As the mode wavelength approaches the cutoff of the linear spectrum (short-wavelength modes), it is found that the three-dimensional (3D) terminal bubble velocity greatly exceeds both the two-dimensional (2D) value and the classical 3D bubble velocity. Unlike in 2D, the 3D short-wavelength bubble velocity does not saturate. The growing 3D bubble acceleration is driven by the unbounded accumulation of vorticity inside the bubble. As a result, the vorticity is transferred by mass ablation from the Rayleigh-Taylor spikes to the ablated plasma filling the bubble volume.

Effect of disjoining pressure on terminal velocity of a bubble sliding along an inclined wall

PubMed Central

Del Castillo, Lorena A.; Ohnishi, Satomi; White, Lee R.; Carnie, Steven L.; Horn, Roger G.

2011-01-01

The influence of salt concentration on the terminal velocities of gravity-driven single bubbles sliding along an inclined glass wall has been investigated, in an effort to establish whether surface forces acting between the wall and the bubble influence the latter’s mobility. A simple sliding bubble apparatus was employed to measure the terminal velocities of air bubbles with radii ranging from 0.3 to 1.5 mm sliding along the interior wall of an inclined Pyrex glass cylinder with inclination angles between 0.6 and 40.1°. Experiments were performed in pure water, 10 mM and 100 mM KCl solutions. We compared our experimental results with a theory by Hodges et al. [1] which considers hydrodynamic forces only, and with a theory developed by two of us [2] which considers surface forces to play a significant role. Our experimental results demonstrate that the terminal velocity of the bubble not only varies with the angle of inclination and the bubble size but also with the salt concentration, particularly at low inclination angles of ∼1–5°, indicating that double-layer forces between the bubble and the wall influence the sliding behavior. This is the first demonstration that terminal velocities of sliding bubbles are affected by disjoining pressure. PMID:21924429

A compact x-ray system for two-phase flow measurement

NASA Astrophysics Data System (ADS)

Song, Kyle; Liu, Yang

2018-02-01

In this paper, a compact x-ray densitometry system consisting of a 50 kV, 1 mA x-ray tube and several linear detector arrays is developed for two-phase flow measurement. The system is capable of measuring void fraction and velocity distributions with a spatial resolution of 0.4 mm per pixel and a frequency of 1000 Hz. A novel measurement model has been established for the system which takes account of the energy spectrum of x-ray photons and the beam hardening effect. An improved measurement accuracy has been achieved with this model compared with the conventional log model that has been widely used in the literature. Using this system, void fraction and velocity distributions are measured for a bubbly and a slug flow in a 25.4 mm I.D. air-water two-phase flow test loop. The measured superficial gas velocities show an error within  ±4% when compared with the gas flowmeter for both conditions.

Close entrainment of massive molecular gas flows by radio bubbles in the central galaxy of Abell 1795

NASA Astrophysics Data System (ADS)

Russell, H. R.; McNamara, B. R.; Fabian, A. C.; Nulsen, P. E. J.; Combes, F.; Edge, A. C.; Hogan, M. T.; McDonald, M.; Salomé, P.; Tremblay, G.; Vantyghem, A. N.

2017-12-01

We present new ALMA observations tracing the morphology and velocity structure of the molecular gas in the central galaxy of the cluster Abell 1795. The molecular gas lies in two filaments that extend 5-7 kpc to the N and S from the nucleus and project exclusively around the outer edges of two inner radio bubbles. Radio jets launched by the central active galactic nucleus have inflated bubbles filled with relativistic plasma into the hot atmosphere surrounding the central galaxy. The N filament has a smoothly increasing velocity gradient along its length from the central galaxy's systemic velocity at the nucleus to -370 km s^{-1}, the average velocity of the surrounding galaxies, at the furthest extent. The S filament has a similarly smooth but shallower velocity gradient and appears to have partially collapsed in a burst of star formation. The close spatial association with the radio lobes, together with the ordered velocity gradients and narrow velocity dispersions, shows that the molecular filaments are gas flows entrained by the expanding radio bubbles. Assuming a Galactic XCO factor, the total molecular gas mass is 3.2 ± 0.2 × 109 M⊙. More than half lies above the N radio bubble. Lifting the molecular clouds appears to require an infeasibly efficient coupling between the molecular gas and the radio bubble. The energy required also exceeds the mechanical power of the N radio bubble by a factor of 2. Stimulated feedback, where the radio bubbles lift low-entropy X-ray gas that becomes thermally unstable and rapidly cools in situ, provides a plausible model. Multiple generations of radio bubbles are required to lift this substantial gas mass. The close morphological association then indicates that the cold gas either moulds the newly expanding bubbles or is itself pushed aside and shaped as they inflate.

Marangoni bubble motion in zero gravity. [Lewis zero gravity drop tower

NASA Technical Reports Server (NTRS)

Thompson, R. L.; Dewitt, K. J.

1979-01-01

It was shown experimentally that the Marangoni phenomenon is a primary mechanism for the movement of a gas bubble in a nonisothermal liquid in a low gravity environment. A mathematical model consisting of the Navier-Stokes and thermal energy equations, together with the appropriate boundary conditions for both media, is presented. Parameter perturbation theory is used to solve this boundary value problem; the expansion parameter is the Marangoni number. The zeroth, first, and second order approximations for the velocity, temperature and pressure distributions in the liquid and in the bubble, and the deformation and terminal velocity of the bubble are determined. Experimental zero gravity data for a nitrogen bubble in ethylene glycol, ethanol, and silicone oil subjected to a linear temperature gradient were obtained using the NASA Lewis zero gravity drop tower. Comparison of the zeroth order analytical results for the bubble terminal velocity showed good agreement with the experimental measurements. The first and second order solutions for the bubble deformation and bubble terminal velocity are valid for liquids having Prandtl numbers on the order of one, but there is a lack of appropriate data to test the theory fully.

Time-resolved fast-neutron radiography of air-water two-phase flows in a rectangular channel by an improved detection system

NASA Astrophysics Data System (ADS)

Zboray, Robert; Dangendorf, Volker; Mor, Ilan; Bromberger, Benjamin; Tittelmeier, Kai

2015-07-01

In a previous work, we have demonstrated the feasibility of high-frame-rate, fast-neutron radiography of generic air-water two-phase flows in a 1.5 cm thick, rectangular flow channel. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany, using an multi-frame, time-resolved detector developed for fast neutron resonance radiography. The results were however not fully optimal and therefore we have decided to modify the detector and optimize it for the given application, which is described in the present work. Furthermore, we managed to improve the image post-processing methodology and the noise suppression. Using the tailored detector and the improved post-processing, significant increase in the image quality and an order of magnitude lower exposure times, down to 3.33 ms, have been achieved with minimized motion artifacts. Similar to the previous study, different two-phase flow regimes such as bubbly slug and churn flows have been examined. The enhanced imaging quality enables an improved prediction of two-phase flow parameters like the instantaneous volumetric gas fraction, bubble size, and bubble velocities. Instantaneous velocity fields around the gas enclosures can also be more robustly predicted using optical flow methods as previously.

Pore Size Control in Aluminium Foam by Standardizing Bubble Rise Velocity and Melt Viscosity

NASA Astrophysics Data System (ADS)

Avinash, G.; Harika, V.; Sandeepika, Ch; Gupta, N.

2018-03-01

In recent years, aluminium foams have found use in a wide range of applications. The properties of these foams, as good structural strength with light weight have made them as a promising structural material for aerospace industry. Foaming techniques (direct and indirect) are used to produce these foams. Direct foaming involves blowing of gas to create gas bubbles in the melt whereas indirect foaming technique uses blowing agents as metallic hydrides, which create hydrogen bubbles. Porosity and its distribution in foams directly affect its properties. This demands for more theoretical studies, to control such cellular structure and hence properties. In present work, we have studied the effect of gas bubble rise velocity and melt viscosity, on pore size and its distribution in aluminium foam. A 15 PPI aluminium foam, prepared using indirect foaming technique having porosity ~86 % was used for study. In order to obtain metal foam, the bubble must not escape from the melt and should get entrapped during solidification. Our calculations suggest that bubble rise velocity and melt viscosity are responsible for vertical displacement of bubble in the melt. It is observed that melt viscosity opposes bubble rise velocity and help the bubbles to stay in the melt, resulting in porous structure.

Proof of concept of the CaO/Ca(OH)2 reaction in a continuous heat-exchanger BFB reactor for thermochemical heat storage in CSP plants

NASA Astrophysics Data System (ADS)

Rougé, Sylvie; Criado, Yolanda A.; Huille, Arthur; Abanades, J. Carlos

2017-06-01

The CaO/Ca(OH)2 hydration/dehydration reaction has long been identified as a attractive method for storing CSP heat. However, the technology applications are still at laboratory scale (TG or small fixed beds). The objective of this work is to investigate the hydration and dehydration reactions performance in a bubbling fluidized bed (BFB) which offers a good potential with regards to heat and mass transfers and upscaling at industrial level. The reactions are first investigated in a 5.5 kW batch BFB, the main conditions are the bed temperature (400-500°C), the molar fraction of steam in the fluidizing gas (0-0.8), the fluidizing gas velocity (0.2-0.7 m/s) and the mass of lime in the batch (1.5-3.5 kg). To assist in the interpretation of the experimental results, a standard 1D bubbling reactor model is formulated and fitted to the experimental results. The results indicate that the hydration reaction is mainly controlled by the slow kinetics of the CaO material tested while significant emulsion-bubble mass-transfer resistances are identified during dehydration due to the much faster dehydration kinetics. In the continuity of these preliminary investigations, a continuous 15.5 kW BFB set-up has been designed, manufactured and started with the objective to operate the hydration and dehydration reactions in steady state during a few hours, and to investigate conditions of faster reactivity such as higher steam molar fractions (up to 1), temperatures (up to 600°C) and velocities (up to 1.5 m/s).

Cavitation Bubble Streaming in Ultrasonic-Standing-Wave Field

NASA Astrophysics Data System (ADS)

Nomura, Shinfuku; Mukasa, Shinobu; Kuroiwa, Masaya; Okada, Yasuyuki; Murakami, Koichi

2005-05-01

The mechanism of cavitation bubble streaming by ultrasonic vibration in a water tank was experimentally investigated. A standard ultrasonic cleaner unit with a resonant frequency of 40 kHz was used as an ultrasonic generator. The behavior of the streaming was visualized by the schlieren method and sonochemical luminescence, and the velocity of the streaming was measured by laser Doppler velocity measurement equipment (LDV). The cavitation bubble streaming has two structures. A cavitation cloud, which consists of many cavitation bubbles, is shaped like a facing pair of bowls with a diameter of approximately 1/3 the wavelength of the standing wave, and moves inside the standing-wave field with a velocity of 30 to 60 mm/s. The cavitation bubbles move intensely in the cloud with a velocity of 5 m/s at an ultrasonic output power of 75 W. The streaming is completely different from conventional acoustic streaming. Also the cavitation bubble is generated neither at the pressure node nor at the antinode.

Analysis of the three-dimensional structure of a bubble wake using PIV and Galilean decomposition

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hassan, Y.A.; Schmidl, W.D.; Ortiz-Villafuerte, J.

1999-07-01

Bubbly flow plays a key role in a variety of natural and industrial processes. An accurate and complete description of the phase interactions in two-phase bubbly flow is not available at this time. These phase interactions are, in general, always three-dimensional and unsteady. Therefore, measurement techniques utilized to obtain qualitative and quantitative data from two-phase flow should be able to acquire transient and three-dimensional data, in order to provide information to test theoretical models and numerical simulations. Even for dilute bubble flows, in which bubble interaction is at a minimum, the turbulent motion of the liquid generated by the bubblemore » is yet to be completely understood. For many years, the design of systems with bubbly flows was based primarily on empiricism. Dilute bubbly flows are an extension of single bubble dynamics, and therefore improvements in the description and modeling of single bubble motion, the flow field around the bubble, and the dynamical interactions between the bubble and the flow will consequently improve bubbly flow modeling. The improved understanding of the physical phenomena will have far-reaching benefits in upgrading the operation and efficiency of current processes and in supporting the development of new and innovative approaches. A stereoscopic particle image velocimetry measurement of the flow generated by the passage of a single air-bubble rising in stagnant water, in a circular pipe is presented. Three-dimensional velocity fields within the measurement zone were obtained. Ensemble-averaged instantaneous velocities for a specific bubble path were calculated and interpolated to obtain mean three-dimensional velocity fields. A Galilean velocity decomposition is used to study the vorticity generated in the flow.« less

Monitoring and analyzing waste glass compositions

DOEpatents

Schumacher, R.F.

1994-03-01

A device and method are described for determining the viscosity of a fluid, preferably molten glass. The apparatus and method use the velocity of rising bubbles, preferably helium bubbles, within the molten glass to determine the viscosity of the molten glass. The bubbles are released from a tube positioned below the surface of the molten glass so that the bubbles pass successively between two sets of electrodes, one above the other, that are continuously monitoring the conductivity of the molten glass. The measured conductivity will change as a bubble passes between the electrodes enabling an accurate determination of when a bubble has passed between the electrodes. The velocity of rising bubbles can be determined from the time interval between a change in conductivity of the first electrode pair and the second, upper electrode pair. The velocity of the rise of the bubbles in the glass melt is used in conjunction with other physical characteristics, obtained by known methods, to determine the viscosity of the glass melt fluid and, hence, glass quality. 2 figures.

Monitoring and analyzing waste glass compositions

DOEpatents

Schumacher, Ray F.

1994-01-01

A device and method for determining the viscosity of a fluid, preferably molten glass. The apparatus and method uses the velocity of rising bubbles, preferably helium bubbles, within the molten glass to determine the viscosity of the molten glass. The bubbles are released from a tube positioned below the surface of the molten glass so that the bubbles pass successively between two sets of electrodes, one above the other, that are continuously monitoring the conductivity of the molten glass. The measured conductivity will change as a bubble passes between the electrodes enabling an accurate determination of when a bubble has passed between the electrodes. The velocity of rising bubbles can be determined from the time interval between a change in conductivity of the first electrode pair and the second, upper electrode pair. The velocity of the rise of the bubbles in the glass melt is used in conjunction with other physical characteristics, obtained by known methods, to determine the viscosity of the glass melt fluid and, hence, glass quality.

Experimental Study on the Flow Regimes and Pressure Gradients of Air-Oil-Water Three-Phase Flow in Horizontal Pipes

PubMed Central

Al-Hadhrami, Luai M.; Shaahid, S. M.; Tunde, Lukman O.; Al-Sarkhi, A.

2014-01-01

An experimental investigation has been carried out to study the flow regimes and pressure gradients of air-oil-water three-phase flows in 2.25 ID horizontal pipe at different flow conditions. The effects of water cuts, liquid and gas velocities on flow patterns and pressure gradients have been studied. The experiments have been conducted at 20°C using low viscosity Safrasol D80 oil, tap water and air. Superficial water and oil velocities were varied from 0.3 m/s to 3 m/s and air velocity varied from 0.29 m/s to 52.5 m/s to cover wide range of flow patterns. The experiments were performed for 10% to 90% water cuts. The flow patterns were observed and recorded using high speed video camera while the pressure drops were measured using pressure transducers and U-tube manometers. The flow patterns show strong dependence on water fraction, gas velocities, and liquid velocities. The observed flow patterns are stratified (smooth and wavy), elongated bubble, slug, dispersed bubble, and annular flow patterns. The pressure gradients have been found to increase with the increase in gas flow rates. Also, for a given superficial gas velocity, the pressure gradients increased with the increase in the superficial liquid velocity. The pressure gradient first increases and then decreases with increasing water cut. In general, phase inversion was observed with increase in the water cut. The experimental results have been compared with the existing unified Model and a good agreement has been noticed. PMID:24523645

Experimental study on the flow regimes and pressure gradients of air-oil-water three-phase flow in horizontal pipes.

PubMed

Al-Hadhrami, Luai M; Shaahid, S M; Tunde, Lukman O; Al-Sarkhi, A

2014-01-01

An experimental investigation has been carried out to study the flow regimes and pressure gradients of air-oil-water three-phase flows in 2.25 ID horizontal pipe at different flow conditions. The effects of water cuts, liquid and gas velocities on flow patterns and pressure gradients have been studied. The experiments have been conducted at 20 °C using low viscosity Safrasol D80 oil, tap water and air. Superficial water and oil velocities were varied from 0.3 m/s to 3 m/s and air velocity varied from 0.29 m/s to 52.5 m/s to cover wide range of flow patterns. The experiments were performed for 10% to 90% water cuts. The flow patterns were observed and recorded using high speed video camera while the pressure drops were measured using pressure transducers and U-tube manometers. The flow patterns show strong dependence on water fraction, gas velocities, and liquid velocities. The observed flow patterns are stratified (smooth and wavy), elongated bubble, slug, dispersed bubble, and annular flow patterns. The pressure gradients have been found to increase with the increase in gas flow rates. Also, for a given superficial gas velocity, the pressure gradients increased with the increase in the superficial liquid velocity. The pressure gradient first increases and then decreases with increasing water cut. In general, phase inversion was observed with increase in the water cut. The experimental results have been compared with the existing unified Model and a good agreement has been noticed.

Multiple bubbles in a Hele-Shaw cell

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vasconcelos, G.L.

A new class of exact solutions is reported for an infinite stream of identical groups of bubbles moving with a constant velocity [ital U] in a Hele-Shaw cell when surface tension is neglected. It is suggested that the existence of these solutions might explain some of the complex behavior observed in recent experiments on rising bubbles in a Hele-Shaw cell. Solutions for a finite number of bubbles in a channel are also obtained. In this case, it is shown that solutions with an arbitrary bubble velocity [ital U][gt][ital V], where [ital V] is the fluid velocity at infinity, can inmore » general be obtained from a simple transformation of the solutions for [ital U]=2[ital V].« less

Effect of the Trendelenburg position on the distribution of arterial air emboli in dogs

NASA Technical Reports Server (NTRS)

Butler, Bruce D.; Laine, Glen A.; Leiman, Basil C.; Warters, Dave; Kurusz, Mark

1988-01-01

The effect of Trendelenburg position (TP) on the distribution of arterial air emboli in dogs was examined in a two-part investigation. In the first part, the effects of the bubble size and the vessel angle on the bubble velocity and the direction of flow were investigated in vitro, using a simulated carotid artery preparation. It was found that larger bubbles increased in velocity in the same direction as the blood flow at 0-, 10-, and 30-deg vessel angles, and decreased when the vessel was positioned at 90 deg. Smaller bubbles did not change velocity from 0 to 30 deg, but acted to increase the velocity, in the same direction as the flood flow, at 90 deg. The second series of experiments examined the effect of 0 to 30 deg TP on carotid-artery distribution of gas bubbles injected into the left ventricle or ascending aorta of anesthetized dogs. It was found that, regardless of the degree of the TP, the bubbles passed into the carotid artery simultaneously with the passage into the abdominal aorta. It is concluded that the TP does not prevent arterial bubbles from reaching the brain.

Time-resolved Fast Neutron Radiography of Air-water Two-phase Flows

NASA Astrophysics Data System (ADS)

Zboray, Robert; Dangendorf, Volker; Mor, Ilan; Tittelmeier, Kai; Bromberger, Benjamin; Prasser, Horst-Michael

Neutron imaging, in general, is a useful technique for visualizing low-Z materials (such as water or plastics) obscured by high-Z materials. However, when significant amounts of both materials are present and full-bodied samples have to be examined, cold and thermal neutrons rapidly reach their applicability limit as the samples become opaque. In such cases one can benefit from the high penetrating power of fast neutrons. In this work we demonstrate the feasibility of time-resolved, fast neutron radiography of generic air-water two-phase flows in a 1.5 cm thick flow channel with Aluminum walls and rectangular cross section. The experiments have been carried out at the high-intensity, white-beam facility of the Physikalisch-Technische Bundesanstalt, Germany. Exposure times down to 3.33 ms have been achieved at reasonable image quality and acceptable motion artifacts. Different two-phase flow regimes such as bubbly slug and churn flows have been examined. Two-phase flow parameters like the volumetric gas fraction, bubble size and bubble velocities have been measured.

The relation between pre-eruptive bubble size distribution, ash particle morphology, and their internal density: Implications to volcanic ash transport and dispersion models

NASA Astrophysics Data System (ADS)

Proussevitch, Alexander

2014-05-01

Parameterization of volcanic ash transport and dispersion (VATD) models strongly depends on particle morphology and their internal properties. Shape of ash particles affects terminal fall velocities (TFV) and, mostly, dispersion. Internal density combined with particle size has a very strong impact on TFV and ultimately on the rate of ash cloud thinning and particle sedimentation on the ground. Unlike other parameters, internal particle density cannot be measured directly because of the micron scale sizes of fine ash particles, but we demonstrate that it varies greatly depending on the particle size. Small simple type ash particles (fragments of bubble walls, 5-20 micron size) do not contain whole large magmatic bubbles inside and their internal density is almost the same as that of volcanic glass matrix. On the other side, the larger compound type ash particles (>40 microns for silicic fine ashes) always contain some bubbles or the whole spectra of bubble size distribution (BSD), i.e. bubbles of all sizes, bringing their internal density down as compared to simple ash. So, density of the larger ash particles is a function of the void fraction inside them (magmatic bubbles) which, in turn, is controlled by BSD. Volcanic ash is a product of the fragmentation of magmatic foam formed by pre-eruptive bubble population and characterized by BSD. The latter can now be measured from bubble imprints on ash particle surfaces using stereo-scanning electron microscopy (SSEM) and BubbleMaker software developed at UNH, or using traditional high-resolution X-Ray tomography. In this work we present the mathematical and statistical formulation for this problem connecting internal ash density with particle size and BSD, and demonstrate how the TFV of the ash population is affected by variation of particle density.

Methodology Development of a Gas-Liquid Dynamic Flow Regime Transition Model

NASA Astrophysics Data System (ADS)

Doup, Benjamin Casey

Current reactor safety analysis codes, such as RELAP5, TRACE, and CATHARE, use flow regime maps or flow regime transition criteria that were developed for static fully-developed two-phase flows to choose interfacial transfer models that are necessary to solve the two-fluid model. The flow regime is therefore difficult to identify near the flow regime transitions, in developing two-phase flows, and in transient two-phase flows. Interfacial area transport equations were developed to more accurately predict the dynamic nature of two-phase flows. However, other model coefficients are still flow regime dependent. Therefore, an accurate prediction of the flow regime is still important. In the current work, the methodology for the development of a dynamic flow regime transition model that uses the void fraction and interfacial area concentration obtained by solving three-field the two-fluid model and two-group interfacial area transport equation is investigated. To develop this model, detailed local experimental data are obtained, the two-group interfacial area transport equations are revised, and a dynamic flow regime transition model is evaluated using a computational fluid dynamics model. Local experimental data is acquired for 63 different flow conditions in bubbly, cap-bubbly, slug, and churn-turbulent flow regimes. The measured parameters are the group-1 and group-2 bubble number frequency, void fraction, interfacial area concentration, and interfacial bubble velocities. The measurements are benchmarked by comparing the prediction of the superficial gas velocities, determined using the local measurements with those determined from volumetric flow rate measurements and the agreement is generally within +/-20%. The repeatability four-sensor probe construction process is within +/-10%. The repeatability of the measurement process is within +/-7%. The symmetry of the test section is examined and the average agreement is within +/-5.3% at z/D = 10 and +/-3.4% at z/D = 32. Revised source/sink terms for the two-group interfacial area transport equations are derived and fit to area-averaged experimental data to determine new model coefficients. The average agreement between this model and the experiment data for the void fraction and interfacial area concentration is 10.6% and 15.7%, respectively. This revised two-group interfacial area transport equation and the three-field two-fluid model are used to solve for the group-1 and group-2 interfacial area concentration and void fraction. These values and a dynamic flow regime transition model are used to classify the flow regimes. The flow regimes determined using this model are compared with the flow regimes based on the experimental data and on a flow regime map using Mishima and Ishii's (1984) transition criteria. The dynamic flow regime transition model is shown to predict the flow regimes dynamically and has improved the prediction of the flow regime over that using a flow regime map. Safety codes often employ the one-dimensional two-fluid model to model two-phase flows. The area-averaged relative velocity correlation necessary to close this model is derived from the drift flux model. The effects of the necessary assumptions used to derive this correlation are investigated using local measurements and these effects are found to have a limited impact on the prediction of the area-averaged relative velocity.

Gas transfer in a bubbly wake flow

NASA Astrophysics Data System (ADS)

Karn, A.; Gulliver, J. S.; Monson, G. M.; Ellis, C.; Arndt, R. E. A.; Hong, J.

2016-05-01

The present work reports simultaneous bubble size and gas transfer measurements in a bubbly wake flow of a hydrofoil, designed to be similar to a hydroturbine blade. Bubble size was measured by a shadow imaging technique and found to have a Sauter mean diameter of 0.9 mm for a reference case. A lower gas flow rate, greater liquid velocities, and a larger angle of attack all resulted in an increased number of small size bubbles and a reduced weighted mean bubble size. Bubble-water gas transfer is measured by the disturbed equilibrium technique. The gas transfer model of Azbel (1981) is utilized to characterize the liquid film coefficient for gas transfer, with one scaling coefficient to reflect the fact that characteristic turbulent velocity is replaced by cross-sectional mean velocity. The coefficient was found to stay constant at a particular hydrofoil configuration while it varied within a narrow range of 0.52-0.60 for different gas/water flow conditions.

Numerical simulation of bubble plumes and an analysis of their seismic attributes

NASA Astrophysics Data System (ADS)

Li, Canping; Gou, Limin; You, Jiachun

2017-04-01

To study the bubble plume's seismic response characteristics, the model of a plume water body has been built in this article using the bubble-contained medium acoustic velocity model and the stochastic medium theory based on an analysis of both the acoustic characteristics of a bubble-contained water body and the actual features of a plume. The finite difference method is used for forward modelling, and the single-shot seismic record exhibits the characteristics of a scattered wave field generated by a plume. A meaningful conclusion is obtained by extracting seismic attributes from the pre-stack shot gather record of a plume. The values of the amplitude-related seismic attributes increase greatly as the bubble content goes up, and changes in bubble radius will not cause seismic attributes to change, which is primarily observed because the bubble content has a strong impact on the plume's acoustic velocity, while the bubble radius has a weak impact on the acoustic velocity. The above conclusion provides a theoretical reference for identifying hydrate plumes using seismic methods and contributes to further study on hydrate decomposition and migration, as well as on distribution of the methane bubble in seawater.

Isotopic fractionation of volatile species during bubble growth in magmas

NASA Astrophysics Data System (ADS)

Watson, E. B.

2016-12-01

Bubbles grow in decompressing magmas by simple expansion and also by diffusive supply of volatiles to the bubble/melt interface. The latter phenomenon is of significant geochemical interest because diffusion can fractionate isotopes, raising the possibility that the isotopic character of volatile components in bubbles may not reflect that of volatiles dissolved in the host melt over the lifetime of a bubble—even in the complete absence of equilibrium vapor/melt isotopic fractionation. None of the foregoing is conceptually new, but recent experimental studies have established the existence of isotope mass effects on diffusion in silicate melts for several elements (Li, Mg, Ca, Fe), and this finding has now been extended to the volatile (anionic) element chlorine (Fortin et al. 2016; this meeting). Knowledge of isotope mass effects on diffusion of volatile species opens the way for quantitative models of diffusive fractionation during bubble growth. Significantly different effects are anticipated for "passive" volatiles (e.g., noble gases and Cl) that are partitioned into existing bubbles but play little role in nucleation and growth, as opposed to "active" volatiles whose limited solubilities lead to bubble nucleation during magma decompression. Numerical solution of the appropriate diffusion/mass-conservation equations reveals that the isotope effect on passive volatiles partitioned into bubbles growing at a constant rate in a static system depends (predictably) upon R/D, Kd and D1/D2 (R = growth rate; D = diffusivity; Kd = bubble/melt partition coefficient; D1/D2 = diffusivity ratio of the isotopes of interest). Constant R is unrealistic, but other scenarios can be explored by including the solubility and EOS of an "active" volatile (e.g., CO2) in numerical simulations of bubble growth. For plausible decompression paths, R increases exponentially with time—leading, potentially, to larger isotopic fractionation of species partitioned into the growing bubble.

Simulations of Bubble Motion in an Oscillating Liquid

NASA Astrophysics Data System (ADS)

Kraynik, A. M.; Romero, L. A.; Torczynski, J. R.

2010-11-01

Finite-element simulations are used to investigate the motion of a gas bubble in a liquid undergoing vertical vibration. The effect of bubble compressibility is studied by comparing "compressible" bubbles that obey the ideal gas law with "incompressible" bubbles that are taken to have constant volume. Compressible bubbles exhibit a net downward motion away from the free surface that does not exist for incompressible bubbles. Net (rectified) velocities are extracted from the simulations and compared with theoretical predictions. The dependence of the rectified velocity on ambient gas pressure, bubble diameter, and bubble depth are in agreement with the theory. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

An analytical approach to the rise velocity of periodic bubble trains in non-Newtonian fluids.

PubMed

Frank, X; Li, H Z; Funfschilling, D

2005-01-01

The present study aims at providing insight into the acceleration mechanism of a bubble chain rising in shear-thinning viscoelastic fluids. The experimental investigation by the Particle Image Velocimetry (PIV), birefringence visualisation and rheological simulation shows that two aspects are central to bubble interactions in such media: the stress creation by the passage of bubbles, and their relaxation due to the fluid's memory forming an evanescent corridor of reduced viscosity. Interactions between bubbles were taken into account mainly through a linear superposition of the stress evolution behind each bubble. An analytical approach together with the rheological consideration was developed to compute the rise velocity of a bubble chain in function of the injection period and bubble volume. The model predictions compare satisfactorily with the experimental investigation.

Geomagnetic conjugate observations of plasma bubbles and thermospheric neutral winds at equatorial latitudes

NASA Astrophysics Data System (ADS)

Fukushima, D.; Shiokawa, K.; Otsuka, Y.; Nishioka, M.; Kubota, M.; Tsugawa, T.; Nagatsuma, T.

2012-12-01

Plasma bubbles are plasma-density depletion which is developed by the Rayleigh-Taylor instability on the sunset terminator at equatorial latitudes. They usually propagate eastward after the sunset. The eastward propagation of the plasma bubbles is considered to be controlled by background eastward neutral winds in the thermosphere through the F-region dynamo effect. However, it is not clear how the F-region dynamo effect contributes to the propagation of the plasma bubbles, because plasma bubbles and background neutral winds have not been simultaneously observed at geomagnetic conjugate points in the northern and southern hemispheres. In this study, geomagnetic conjugate observations of the plasma bubbles at low latitudes with thermospheric neutral winds were reported. The plasma bubbles were observed at Kototabang (0.2S, 100.3E, geomagnetic latitude (MLAT): 10.0S), Indonesia and at Chiang Mai (18.8N, 98.9E, MLAT: 8.9N), Thailand, which are geomagnetic conjugate stations, on 5 April, 2011 from 13 to 22 UT (from 20 to 05 LT). These plasma bubbles were observed in the 630-nm airglow images taken by using highly-sensitive all-sky airglow imagers at both stations. They propagated eastward with horizontal velocities of about 100-125 m/s. Background thermospheric neutral winds were also observed at both stations by using two Fabry-Perot interferometers (FPIs). The eastward wind velocities were about 70-130 m/s at Kototabang, and about 50-90 m/s at Chiang Mai. We estimated ion drift velocities by using these neutral winds observed by FPIs and conductivities calculated from the IRI and MSIS models. The estimated velocities were about 60-90 % of the drift velocities of plasma bubbles. This result shows that most of the plasma bubble drift can be explained by the F-region dynamo effect, and additional electric field effect may come in to play.

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. Part 2: Stability

NASA Technical Reports Server (NTRS)

Tanveer, Saleh

1989-01-01

The analysis is extended to determine the linear stability of a bubble in a Hele-Shaw cell analytically. Only the solution branch corresponding to largest possible bubble velocity U for given surface tension is found to be stable, while all the others are unstable, in accordance with earlier numerical results.

Field-driven chiral bubble dynamics analysed by a semi-analytical approach

NASA Astrophysics Data System (ADS)

Vandermeulen, J.; Leliaert, J.; Dupré, L.; Van Waeyenberge, B.

2017-12-01

Nowadays, field-driven chiral bubble dynamics in the presence of the Dzyaloshinskii-Moriya interaction are a topic of thorough investigation. In this paper, a semi-analytical approach is used to derive equations of motion that express the bubble wall (BW) velocity and the change in in-plane magnetization angle as function of the micromagnetic parameters of the involved interactions, thereby taking into account the two-dimensional nature of the bubble wall. It is demonstrated that the equations of motion enable an accurate description of the expanding and shrinking convex bubble dynamics and an expression for the transition field between shrinkage and expansion is derived. In addition, these equations of motion show that the BW velocity is not only dependent on the driving force, but also on the BW curvature. The absolute BW velocity increases for both a shrinking and an expanding bubble, but for different reasons: for expanding bubbles, it is due to the increasing importance of the driving force, while for shrinking bubbles, it is due to the increasing importance of contributions related to the BW curvature. Finally, using this approach we show how the recently proposed magnetic bubblecade memory can operate in the flow regime in the presence of a tilted sinusoidal magnetic field and at greatly reduced bubble sizes compared to the original device prototype.

Bubble coalescence in a Newtonian fluid

NASA Astrophysics Data System (ADS)

Garg, Vishrut; Basaran, Osman

2017-11-01

Bubble coalescence plays a central role in the hydrodynamics of gas-liquid systems such as bubble column reactors, spargers, and foams. Two bubbles approaching each other at velocity V coalesce when the thin film between them ruptures, which is often the rate-limiting step. Experimental studies of this system are difficult, and recent works provide conflicting results on the effect of V on coalescence times. We simulate the head-on approach of two bubbles of equal radii R in an incompressible Newtonian fluid (density ρ, viscosity μ, and surface tension σ) by solving numerically the free boundary problem comprised of the Navier Stokes and continuity equations. Simulations are made challenging by the existence of highly disparate lengthscales, i.e. film thickness and drop radii, which are resolved by using the method of elliptic mesh generation. For a given liquid, the bubbles are shown to coalesce for all velocities below a critical value. The effects of Ohnesorge number Oh = μ /√{ ρσR } on coalescence time and critical velocity are also investigated.

Zonal drift velocities of the ionospheric plasma bubbles over brazilian region using oi630nm airglow digital images

NASA Astrophysics Data System (ADS)

Arruda, D. C. S.; Sobral, J. H. A.; Abdu, M. A.; Castilho, V. M.; Takahashi, H.

The zonal drift velocities of the ionospheric plasma bubbles over the Brazilian region are analyzed in this study that is based on OI630nm airglow digital images. These digital images were obtained by an all-sky imager system between October 1998 and August 2000, at Cachoeira Paulista (22.5°S, 45°W), a low latitude region. In this period, 138 nights of OI 630 nm airglow experiments were carried out of which 30 nights detected the ionospheric plasma bubbles. These 30 nights correspond to magnetically quiet days (ΣK_P<24+) and were grouped according approximately to their season. KEY WORDS: Imager System, Ionospheric Plasma Bubbles, Zonal drift velocities, OI630nm.

The elasticity of soap bubbles containing wormlike micelles.

PubMed

Sabadini, Edvaldo; Ungarato, Rafael F S; Miranda, Paulo B

2014-01-28

Slow-motion imaging of the rupture of soap bubbles generally shows the edges of liquid films retracting at a constant speed (known as the Taylor-Culick velocity). Here we investigate soap bubbles formed from simple solutions of a cationic surfactant (cetyltrimethylammonium bromide - CTAB) and sodium salicylate. The interaction of salicylate ions with CTAB leads to the formation of wormlike micelles (WLM), which yield a viscoelastic behavior to the liquid film of the bubble. We demonstrate that these elastic bubbles collapse at a velocity up to 30 times higher than the Taylor-Culick limit, which has never been surpassed. This is because during the bubble inflation, the entangled WLM chains stretch, storing elastic energy. This extra energy is then released during the rupture of the bubble, yielding an additional driving force for film retraction (besides surface tension). This new mechanism for the bursting of elastic bubbles may have important implications to the breakup of viscoelastic sprays in industrial applications.

A fractional Fourier transform analysis of a bubble excited by an ultrasonic chirp.

PubMed

Barlow, Euan; Mulholland, Anthony J

2011-11-01

The fractional Fourier transform is proposed here as a model based, signal processing technique for determining the size of a bubble in a fluid. The bubble is insonified with an ultrasonic chirp and the radiated pressure field is recorded. This experimental bubble response is then compared with a series of theoretical model responses to identify the most accurate match between experiment and theory which allows the correct bubble size to be identified. The fractional Fourier transform is used to produce a more detailed description of each response, and two-dimensional cross correlation is then employed to identify the similarities between the experimental response and each theoretical response. In this paper the experimental bubble response is simulated by adding various levels of noise to the theoretical model output. The method is compared to the standard technique of using time-domain cross correlation. The proposed method is shown to be far more robust at correctly sizing the bubble and can cope with much lower signal to noise ratios.

Close relationship between a dry-wet transition and a bubble rearrangement in two-dimensional foam

PubMed Central

Furuta, Yujiro; Oikawa, Noriko; Kurita, Rei

2016-01-01

Liquid foams are classified into a dry foam and a wet foam, empirically judging from the liquid fraction or the shape of the gas bubbles. It is known that physical properties such as elasticity and diffusion are different between the dry foam and the wet foam. Nevertheless, definitions of those states have been vague and the dry-wet transition of foams has not been clarified yet. Here we show that the dry-wet transition is closely related to rearrangement of the gas bubbles, by simultaneously analysing the shape change of the bubbles and that of the entire foam in two dimensional foam. In addition, we also find a new state in quite low liquid fraction, which is named “superdry foam”. Whereas the shape change of the bubbles strongly depends on the change of the liquid fraction in the superdry foam, the shape of the bubbles does not change with changing the liquid fraction in the dry foam. Our results elucidate the relationship between the transitions and the macroscopic mechanical properties. PMID:27874060

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, W. H.; HEDPS and CAPT, Peking University, Beijing 100871; Wang, L. F.

In this research, the temporal evolution of the bubble tip velocity in Rayleigh-Taylor instability (RTI) at arbitrary Atwood numbers and different initial perturbation velocities with a discontinuous profile in irrotational, incompressible, and inviscid fluids (i.e., classical RTI) is investigated. Potential models from Layzer [Astrophys. J. 122, 1 (1955)] and perturbation velocity potentials from Goncharov [Phys. Rev. Lett. 88, 134502 (2002)] are introduced. It is found that the temporal evolution of bubble tip velocity [u(t)] depends essentially on the initial perturbation velocity [u(0)]. First, when the u(0)

Thermocapillary Bubble Migration: Thermal Boundary Layers for Large Marangoni Numbers

NASA Technical Reports Server (NTRS)

Balasubramaniam, R.; Subramanian, R. S.

1996-01-01

The migration of an isolated gas bubble in an immiscible liquid possessing a temperature gradient is analyzed in the absence of gravity. The driving force for the bubble motion is the shear stress at the interface which is a consequence of the temperature dependence of the surface tension. The analysis is performed under conditions for which the Marangoni number is large, i.e. energy is transferred predominantly by convection. Velocity fields in the limit of both small and large Reynolds numbers are used. The thermal problem is treated by standard boundary layer theory. The outer temperature field is obtained in the vicinity of the bubble. A similarity solution is obtained for the inner temperature field. For both small and large Reynolds numbers, the asymptotic values of the scaled migration velocity of the bubble in the limit of large Marangoni numbers are calculated. The results show that the migration velocity has the same scaling for both low and large Reynolds numbers, but with a different coefficient. Higher order thermal boundary layers are analyzed for the large Reynolds number flow field and the higher order corrections to the migration velocity are obtained. Results are also presented for the momentum boundary layer and the thermal wake behind the bubble, for large Reynolds number conditions.

Effects of tissue stiffness, ultrasound frequency, and pressure on histotripsy-induced cavitation bubble behavior.

PubMed

Vlaisavljevich, Eli; Lin, Kuang-Wei; Warnez, Matthew T; Singh, Rahul; Mancia, Lauren; Putnam, Andrew J; Johnsen, Eric; Cain, Charles; Xu, Zhen

2015-03-21

Histotripsy is an ultrasound ablation method that controls cavitation to fractionate soft tissue. In order to effectively fractionate tissue, histotripsy requires cavitation bubbles to rapidly expand from nanometer-sized initial nuclei into bubbles often larger than 50 µm. Using a negative pressure high enough to initiate a bubble cloud and expand bubbles to a sufficient size, histotripsy has been shown capable of completely fractionating soft tissue into acelluar debris resulting in effective tissue removal. Previous work has shown that the histotripsy process is affected by tissue mechanical properties with stiffer tissues showing increased resistance to histotripsy fractionation, which we hypothesize to be caused by impeded bubble expansion in stiffer tissues. In this study, the hypothesis that increases in tissue stiffness cause a reduction in bubble expansion was investigated both theoretically and experimentally. High speed optical imaging was used to capture a series of time delayed images of bubbles produced inside mechanically tunable agarose tissue phantoms using histotripsy pulses produced by 345 kHz, 500 kHz, 1.5 MHz, and 3 MHz histotripsy transducers. The results demonstrated a significant decrease in maximum bubble radius (Rmax) and collapse time (tc) with both increasing Young's modulus and increasing frequency. Furthermore, results showed that Rmax was not increased by raising the pressure above the intrinsic threshold. Finally, this work demonstrated the potential of using a dual-frequency strategy to modulate the expansion of histotripsy bubbles. Overall, the results of this study improve our understanding of how tissue stiffness and ultrasound parameters affect histotripsy-induced bubble behavior and provide a rational basis to tailor acoustic parameters for treatment of the specific tissues of interest.

Effects of tissue stiffness, ultrasound frequency, and pressure on histotripsy-induced cavitation bubble behavior

NASA Astrophysics Data System (ADS)

Vlaisavljevich, Eli; Lin, Kuang-Wei; Warnez, Matthew T.; Singh, Rahul; Mancia, Lauren; Putnam, Andrew J.; Johnsen, Eric; Cain, Charles; Xu, Zhen

2015-03-01

Histotripsy is an ultrasound ablation method that controls cavitation to fractionate soft tissue. In order to effectively fractionate tissue, histotripsy requires cavitation bubbles to rapidly expand from nanometer-sized initial nuclei into bubbles often larger than 50 µm. Using a negative pressure high enough to initiate a bubble cloud and expand bubbles to a sufficient size, histotripsy has been shown capable of completely fractionating soft tissue into acelluar debris resulting in effective tissue removal. Previous work has shown that the histotripsy process is affected by tissue mechanical properties with stiffer tissues showing increased resistance to histotripsy fractionation, which we hypothesize to be caused by impeded bubble expansion in stiffer tissues. In this study, the hypothesis that increases in tissue stiffness cause a reduction in bubble expansion was investigated both theoretically and experimentally. High speed optical imaging was used to capture a series of time delayed images of bubbles produced inside mechanically tunable agarose tissue phantoms using histotripsy pulses produced by 345 kHz, 500 kHz, 1.5 MHz, and 3 MHz histotripsy transducers. The results demonstrated a significant decrease in maximum bubble radius (Rmax) and collapse time (tc) with both increasing Young’s modulus and increasing frequency. Furthermore, results showed that Rmax was not increased by raising the pressure above the intrinsic threshold. Finally, this work demonstrated the potential of using a dual-frequency strategy to modulate the expansion of histotripsy bubbles. Overall, the results of this study improve our understanding of how tissue stiffness and ultrasound parameters affect histotripsy-induced bubble behavior and provide a rational basis to tailor acoustic parameters for treatment of the specific tissues of interest.

Convection in a volcanic conduit recorded by bubbles

USGS Publications Warehouse

Carey, Rebecca J.; Manga, Michael; Degruyter, Wim; Gonnermann, Helge M.; Swanson, Donald; Houghton, Bruce F.; Orr, Tim R.; Patrick, Matthew R.

2013-01-01

Microtextures of juvenile pyroclasts from Kīlauea’s (Hawai‘i) early A.D. 2008 explosive activity record the velocity and depth of convection within the basaltic magma-filled conduit. We use X-ray microtomography (μXRT) to document the spatial distribution of bubbles. We find small bubbles (radii from 5 μm to 70 μm) in a halo surrounding larger millimeter-size bubbles. This suggests that dissolved water was enriched around the larger bubbles—the opposite of what is expected if bubbles grow as water diffuses into the bubble. Such volatile enrichment implies that the volatiles within the large bubbles were redissolving into the melt as they descended into the conduit by the downward motion of convecting magma within the lava lake. The thickness of the small bubble halo is ∼100–150 μm, consistent with water diffusing into the melt on time scales on the order of 103 s. Eruptions, triggered by rockfall, rapidly exposed this magma to lower pressures, and the haloes of melt with re-dissolved water became sufficiently supersaturated to cause nucleation of the population of smaller bubbles. The required supersaturation pressures are consistent with a depth of a few hundred meters and convection velocities of the order of 0.1 m s−1, similar to the circulation velocity observed on the surface of the Halema‘uma‘u lava lake.

Research on the Conductivity-Based Detection Principles of Bubbles in Two-Phase Flows and the Design of a Bubble Sensor for CBM Wells.

PubMed

Wu, Chuan; Wen, Guojun; Han, Lei; Wu, Xiaoming

2016-09-17

The parameters of gas-liquid two-phase flow bubbles in field coalbed methane (CBM) wells are of great significance for analyzing coalbed methane output, judging faults in CBM wells, and developing gas drainage and extraction processes, which stimulates an urgent need for detecting bubble parameters for CBM wells in the field. However, existing bubble detectors cannot meet the requirements of the working environments of CBM wells. Therefore, this paper reports findings on the principles of measuring the flow pattern, velocity, and volume of two-phase flow bubbles based on conductivity, from which a new bubble sensor was designed. The structural parameters and other parameters of the sensor were then computed, the "water film phenomenon" produced by the sensor was analyzed, and the appropriate materials for making the sensor were tested and selected. After the sensor was successfully devised, laboratory tests and field tests were performed, and the test results indicated that the sensor was highly reliable and could detect the flow patterns of two-phase flows, as well as the quantities, velocities, and volumes of bubbles. With a velocity measurement error of ±5% and a volume measurement error of ±7%, the sensor can meet the requirements of field use. Finally, the characteristics and deficiencies of the bubble sensor are summarized based on an analysis of the measurement errors and a comparison of existing bubble-measuring devices and the designed sensor.

Research on the Conductivity-Based Detection Principles of Bubbles in Two-Phase Flows and the Design of a Bubble Sensor for CBM Wells

PubMed Central

Wu, Chuan; Wen, Guojun; Han, Lei; Wu, Xiaoming

2016-01-01

The parameters of gas-liquid two-phase flow bubbles in field coalbed methane (CBM) wells are of great significance for analyzing coalbed methane output, judging faults in CBM wells, and developing gas drainage and extraction processes, which stimulates an urgent need for detecting bubble parameters for CBM wells in the field. However, existing bubble detectors cannot meet the requirements of the working environments of CBM wells. Therefore, this paper reports findings on the principles of measuring the flow pattern, velocity, and volume of two-phase flow bubbles based on conductivity, from which a new bubble sensor was designed. The structural parameters and other parameters of the sensor were then computed, the “water film phenomenon” produced by the sensor was analyzed, and the appropriate materials for making the sensor were tested and selected. After the sensor was successfully devised, laboratory tests and field tests were performed, and the test results indicated that the sensor was highly reliable and could detect the flow patterns of two-phase flows, as well as the quantities, velocities, and volumes of bubbles. With a velocity measurement error of ±5% and a volume measurement error of ±7%, the sensor can meet the requirements of field use. Finally, the characteristics and deficiencies of the bubble sensor are summarized based on an analysis of the measurement errors and a comparison of existing bubble-measuring devices and the designed sensor. PMID:27649206

Quantification of Methane Gas Flux and Bubble Fate on the Eastern Siberian Arctic Shelf Utilizing Calibrated Split-beam Echosounder Data.

NASA Astrophysics Data System (ADS)

Weidner, E. F.; Mayer, L. A.; Weber, T. C.; Jerram, K.; Jakobsson, M.; Chernykh, D.; Ananiev, R.; Mohammad, R.; Semiletov, I. P.

2016-12-01

On the Eastern Siberian Arctic Shelf (ESAS) subsea permafrost, shallow gas hydrates, and trapped free gas hold an estimated 1400 Gt of methane. Recent observations of methane bubble plumes and high concentrations of dissolved methane in the water column indicate methane release via ebullition. Methane gas released from the shallow ESAS (<50 m average depth) has high potential to be transported to the atmosphere. To directly and quantitatively address the magnitude of methane flux and the fate of rising bubbles in the ESAS, methane seeps were mapped with a broadband split-beam echosounder as part of the Swedish-Russian-US Arctic Ocean Investigation of Climate-Cryosphere-Carbon Interactions program (SWERUS-C3). Acoustic measurements were made over a broad range of frequencies (16 to 29 kHz). The broad bandwidth provided excellent discrimination of individual targets in the water column, allowing for the identification of single bubbles. Absolute bubble target strength values were determined by compensating apparent target strength measurements for beam pattern effects via standard calibration techniques. The bubble size distribution of seeps with individual bubble signatures was determined by exploiting bubble target strength models over the broad range of frequencies. For denser seeps, with potential higher methane flux, bubble size distribution was determined via extrapolation from seeps in similar geomorphological settings. By coupling bubble size distributions with rise velocity measurements, which are made possible by split-beam target tracking, methane gas flux can be estimated. Of the 56 identified seeps in the SWERUS data set, individual bubbles scatterers were identified in more than half (31) of the seeps. Preliminary bubble size distribution results indicate bubble radii range from 0.75 to 3.0 mm, with relatively constant bubble size distribution throughout the water column. Initial rise velocity observations indicate bubble rise velocity increases with decreasing depth, seemingly independent of bubble radius.

Bubble and Slug Flow at Microgravity Conditions: State of Knowledge and Open Questions

NASA Technical Reports Server (NTRS)

Colin, C.; Fabre, J.; McQuillen, J.

1996-01-01

Based on the experiments carried out over the past decade at microgravity conditions, an overview of our current knowledge of bubbly and slug flows is presented. The transition from bubble to slug flow, the void fraction and the pressure drop are discussed from the data collected in the literature. The transition from bubble to slug flow may be predicted by introducing a critical void fraction that depends on the fluid properties and the pipe diameter; however, the role of coalescence which controls this transition is not clearly understood. The void fraction may be accurately calculated using a drift-flux model. It is shown from local measurements that the drift of the gas with respect to the mixture is due to non-uniform radial distribution of void fraction. The pressure drop happens to be controlled by the liquid flow for bubbly flow whereas for slug flow the experimental results show that pressure drops is larger than expected. From this study, the guidelines for future research in microgravity are given.

Electroweak bubble wall speed limit

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bödeker, Dietrich; Moore, Guy D., E-mail: bodeker@physik.uni-bielefeld.de, E-mail: guymoore@ikp.physik.tu-darmstadt.de

In extensions of the Standard Model with extra scalars, the electroweak phase transition can be very strong, and the bubble walls can be highly relativistic. We revisit our previous argument that electroweak bubble walls can 'run away,' that is, achieve extreme ultrarelativistic velocities γ ∼ 10{sup 14}. We show that, when particles cross the bubble wall, they can emit transition radiation. Wall-frame soft processes, though suppressed by a power of the coupling α, have a significance enhanced by the γ-factor of the wall, limiting wall velocities to γ ∼ 1/α. Though the bubble walls can move at almost the speedmore » of light, they carry an infinitesimal share of the plasma's energy.« less

An iterative fullwave simulation approach to multiple scattering in media with randomly distributed microbubbles

NASA Astrophysics Data System (ADS)

Joshi, Aditya; Lindsey, Brooks D.; Dayton, Paul A.; Pinton, Gianmarco; Muller, Marie

2017-05-01

Ultrasound contrast agents (UCA), such as microbubbles, enhance the scattering properties of blood, which is otherwise hypoechoic. The multiple scattering interactions of the acoustic field with UCA are poorly understood due to the complexity of the multiple scattering theories and the nonlinear microbubble response. The majority of bubble models describe the behavior of UCA as single, isolated microbubbles suspended in infinite medium. Multiple scattering models such as the independent scattering approximation can approximate phase velocity and attenuation for low scatterer volume fractions. However, all current models and simulation approaches only describe multiple scattering and nonlinear bubble dynamics separately. Here we present an approach that combines two existing models: (1) a full-wave model that describes nonlinear propagation and scattering interactions in a heterogeneous attenuating medium and (2) a Paul-Sarkar model that describes the nonlinear interactions between an acoustic field and microbubbles. These two models were solved numerically and combined with an iterative approach. The convergence of this combined model was explored in silico for 0.5 × 106 microbubbles ml-1, 1% and 2% bubble concentration by volume. The backscattering predicted by our modeling approach was verified experimentally with water tank measurements performed with a 128-element linear array transducer. An excellent agreement in terms of the fundamental and harmonic acoustic fields is shown. Additionally, our model correctly predicts the phase velocity and attenuation measured using through transmission and predicted by the independent scattering approximation.

Numerical study of the influence of geometrical characteristics of a vertical helical coil on a bubbly flow

NASA Astrophysics Data System (ADS)

Saffari, H.; Moosavi, R.

2014-11-01

In this article, turbulent single-phase and two-phase (air-water) bubbly fluid flows in a vertical helical coil are analyzed by using computational fluid dynamics (CFD). The effects of the pipe diameter, coil diameter, coil pitch, Reynolds number, and void fraction on the pressure loss, friction coefficient, and flow characteristics are investigated. The Eulerian-Eulerian model is used in this work to simulate the two-phase fluid flow. Three-dimensional governing equations of continuity, momentum, and energy are solved by using the finite volume method. The k- ɛ turbulence model is used to calculate turbulence fluctuations. The SIMPLE algorithm is employed to solve the velocity and pressure fields. Due to the effect of a secondary force in helical pipes, the friction coefficient is found to be higher in helical pipes than in straight pipes. The friction coefficient increases with an increase in the curvature, pipe diameter, and coil pitch and decreases with an increase in the coil diameter and void fraction. The close correlation between the numerical results obtained in this study and the numerical and empirical results of other researchers confirm the accuracy of the applied method. For void fractions up to 0.1, the numerical results indicate that the friction coefficient increases with increasing the pipe diameter and keeping the coil pitch and diameter constant and decreases with increasing the coil diameter. Finally, with an increase in the Reynolds number, the friction coefficient decreases, while the void fraction increases.

The Collapse of Vapor Bubbles in a Spatially Non-Uniform Flow

NASA Technical Reports Server (NTRS)

Hao, Y.; Prosperetti, A.

2000-01-01

Pressure gradients act differently on liquid particles and suspended bubbles and are, therefore, capable of inducing a relative motion between the phases even when no relative velocity initially exists. As a consequence of the enhanced heat transfer in the presence of convection, this fact may have a major impact on the evolution of a vapor bubble. The effect is particularly strong in the case of a collapsing bubble for which, due to the conservation of the system's impulse, the induced relative velocity tends to be magnified when the bubble volume shrinks. A practical application could be, for instance, the enhancement of the condensation rate of bubbles downstream of a heated region, thereby reducing the quality of a flowing liquid-vapor mixture. A simple model of the process, in which the bubble is assumed to be spherical and the flow potential, is developed in the paper.

Formation of soap bubbles by gas jet

NASA Astrophysics Data System (ADS)

Zhou, Maolei; Li, Min; Chen, Zhiyuan; Han, Jifeng; Liu, Dong

2017-12-01

Soap bubbles can be easily generated by various methods, while their formation process is complicated and still worth studying. A model about the bubble formation process was proposed in the study by Salkin et al. [Phys. Rev. Lett. 116, 077801 (2016)] recently, and it was reported that the bubbles were formed when the gas blowing velocity was above one threshold. However, after a detailed study of these experiments, we found that the bubbles could be generated in two velocity ranges which corresponded to the laminar and turbulent gas jet, respectively, and the predicted threshold was only effective for turbulent gas flow. The study revealed that the bubble formation was greatly influenced by the aerodynamics of the gas jet blowing to the film, and these results will help to further understand the formation mechanism of the soap bubble as well as the interaction between the gas jet and the thin liquid film.

Acoustic wave propagation in bubbly flow with gas, vapor or their mixtures.

PubMed

Zhang, Yuning; Guo, Zhongyu; Gao, Yuhang; Du, Xiaoze

2018-01-01

Presence of bubbles in liquids could significantly alter the acoustic waves in terms of wave speed and attenuation. In the present paper, acoustic wave propagation in bubbly flows with gas, vapor and gas/vapor mixtures is theoretically investigated in a wide range of parameters (including frequency, bubble radius, void fraction, and vapor mass fraction). Our finding reveals two types of wave propagation behavior depending on the vapor mass fraction. Furthermore, the minimum wave speed (required for the closure of cavitation modelling in the sonochemical reactor design) is analyzed and the influences of paramount parameters on it are quantitatively discussed. Copyright © 2017 Elsevier B.V. All rights reserved.

The effects of a decompression on seismic parameter profiles in a gas-charged magma

NASA Astrophysics Data System (ADS)

Sturton, Susan; Neuberg, Jürgen

2003-11-01

Seismic velocities in a gas-charged magma vary with depth and time. Relationships between pressure, density, exsolved gas content, and seismic velocity are derived and used in conjunction with expressions describing diffusive bubble growth to find a series of velocity profiles which depend on time. An equilibrium solution is obtained by considering a column of magma in which the gas distribution corresponds to the magmastatic pressure profile with depth. Decompression events of various sizes are simulated, and the resulting disequilibrium between the gas pressure and magmastatic pressure leads to bubble growth and therefore to a change of seismic velocity and density with time. Bubble growth stops when the system reaches a new equilibrium. The corresponding volume increase is accommodated by accelerating the magma column upwards and an extrusion of lava. A timescale for the system to return to equilibrium can be obtained. The effect of changes in magma viscosity and bubble number density is examined.

One-dimensional bubble model of pulsed discharge in water

NASA Astrophysics Data System (ADS)

Lu, XinPei

2007-09-01

In this paper, a one-dimensional bubble model of pulsed discharge in water is presented. With a total input energy of 0.63J, the simulation results show that when the bubble collapses at the center of the bubble, the plasma pressure oscillates strongly. It oscillates between 800 and 1150atm with an oscillation frequency of about 6.9MHz, while at r =R/2 (R: bubble radius), the gas velocity oscillates intensely at the same frequency. It oscillates between -235 and 229m/s when the bubble radius reaches its minimum. But it does not oscillate at r =R because of the inertia of the surrounding water. The bubble collapses and reexpands with almost the same speed as that of the zero-dimensional (0D) model. This further confirms why the shock wave pressure from the 0D mode has a good agreement with the experimental results since the shock wave pressure is only determined by the bubble wall velocity v(R ).

Measurements of Shear Lift Force on a Bubble in Channel Flow in Microgravity

NASA Technical Reports Server (NTRS)

Nahra, Henry K.; Motil, Brian J.; Skor, Mark

2003-01-01

Under microgravity conditions, the shear lift force acting on bubbles, droplets or solid particles in multiphase flows becomes important because under normal gravity, this hydrodynamic force is masked by buoyancy. This force plays an important role in furnishing the detachment process of bubbles in a setting where a bubble suspension is needed in microgravity. In this work, measurements of the shear lift force acting on a bubble in channel flow are performed. The shear lift force is deduced from the bubble kinematics using scaling and then compared with predictions from models in literature that address different asymptotic and numerical solutions. Basic trajectory calculations are then performed and the results are compared with experimental data of position of the bubble in the channel. A direct comparison of the lateral velocity of the bubbles is also made with the lateral velocity prediction from investigators, whose work addressed the shear lift on a sphere in different two-dimensional shear flows including Poiseuille flow.

DOE Office of Scientific and Technical Information (OSTI.GOV)

White, D.B.

This paper reports on experiments to examine gas migration rates in drilling muds that were performed in a 15-m-long, 200-mm-ID inclinable flow loop where air injection simulates gas entry during a kick. These tests were conducted using a xanthum gum (a common polymer used in drilling fluids) solution to simulate drilling muds as the liquid phase and air as the gas phase. This work represents a significant extension of existing correlations for gas/liquid flows in large pipe diameters with non- Newtonian fluids. Bubbles rise faster in drilling muds than in water despite the increased viscosity. This surprising result is causedmore » by the change in the flow regime, with large slug-type bubbles forming at lower void fractions. The gas velocity is independent of void fraction, thus simplifying flow modeling. Results show that a gas influx will rise faster in a well than previously believed. This has major implications for kick simulation, with gas arriving at the surface earlier than would be expected and the gas outflow rate being higher than would have been predicted. A model of the two-phase gas flow in drilling mud, including the results of this work, has been incorporated into the joint Schlumberger Cambridge Research (SCR)/BP Intl. kick model.« less

Lyα-emitting galaxies as a probe of reionization: large-scale bubble morphology and small-scale absorbers

NASA Astrophysics Data System (ADS)

Kakiichi, Koki; Dijkstra, Mark; Ciardi, Benedetta; Graziani, Luca

2016-12-01

The visibility of Lyα-emitting galaxies during the Epoch of Reionization is controlled by both diffuse H I patches in large-scale bubble morphology and small-scale absorbers. To investigate their impacts on Lyα transfer, we apply a novel combination of analytic modelling and cosmological hydrodynamical, radiative transfer simulations to three reionization models: (I) the `bubble' model, where only diffuse H I outside ionized bubbles is present; (II) the `web' model, where H I exists only in overdense self-shielded gas; and (III) the hybrid `web-bubble' model. The three models can explain the observed Lyα luminosity function equally well, but with very different H I fractions. This confirms a degeneracy between the ionization topology of the intergalactic medium (IGM) and the H I fraction inferred from Lyα surveys. We highlight the importance of the clustering of small-scale absorbers around galaxies. A combined analysis of the Lyα luminosity function and the Lyα fraction can break this degeneracy and provide constraints on the reionization history and its topology. Constraints can be improved by analysing the full MUV-dependent redshift evolution of the Lyα fraction of Lyman break galaxies. We find that the IGM-transmission probability distribution function is unimodal for bubble models and bimodal in web models. Comparing our models to observations, we infer that the neutral fraction at z ˜ 7 is likely to be of the order of tens of per cent when interpreted with bubble or web-bubble models, with a conservative lower limit ˜1 per cent when interpreted with web models.

Vertical rise velocity of equatorial plasma bubbles estimated from Equatorial Atmosphere Radar (EAR) observations and HIRB model simulations

NASA Astrophysics Data System (ADS)

Tulasi Ram, S.; Ajith, K. K.; Yokoyama, T.; Yamamoto, M.; Niranjan, K.

2017-06-01

The vertical rise velocity (Vr) and maximum altitude (Hm) of equatorial plasma bubbles (EPBs) were estimated using the two-dimensional fan sector maps of 47 MHz Equatorial Atmosphere Radar (EAR), Kototabang, during May 2010 to April 2013. A total of 86 EPBs were observed out of which 68 were postsunset EPBs and remaining 18 EPBs were observed around midnight hours. The vertical rise velocities of the EPBs observed around the midnight hours are significantly smaller ( 26-128 m/s) compared to those observed in postsunset hours ( 45-265 m/s). Further, the vertical growth of the EPBs around midnight hours ceases at relatively lower altitudes, whereas the majority of EPBs at postsunset hours found to have grown beyond the maximum detectable altitude of the EAR. The three-dimensional numerical high-resolution bubble (HIRB) model with varying background conditions are employed to investigate the possible factors that control the vertical rise velocity and maximum attainable altitudes of EPBs. The estimated rise velocities from EAR observations at both postsunset and midnight hours are, in general, consistent with the nonlinear evolution of EPBs from the HIRB model. The smaller vertical rise velocities (Vr) and lower maximum altitudes (Hm) of EPBs during midnight hours are discussed in terms of weak polarization electric fields within the bubble due to weaker background electric fields and reduced background ion density levels.