Interfacial bubbles formed by plunging thin liquid films in a pool

NASA Astrophysics Data System (ADS)

Salkin, Louis; Schmit, Alexandre; David, Richard; Delvert, Alexandre; Gicquel, Eric; Panizza, Pascal; Courbin, Laurent

2017-06-01

We show that the immersion of a horizontally suspended thin film of liquid in a pool of the same fluid creates an interfacial bubble, that is, a bubble at the liquid-air interface. Varying the fluid properties, the film's size, and its immersion velocity, our experiments unveil two formation regimes characterized by either a visco-capillary or an inertio-capillary mechanism that controls the size of a produced bubble. To rationalize these results, we compare the pressure exerted by the air flow under a plunging film with the Laplace pressure needed to generate film dimpling, which subsequently yields air entrapment and the production of a bubble. This physical model explains the power-law variations of the bubble size with the governing dimensionless number for each regime.

A numerical simulation of the water vapor bubble rising in ferrofluid by volume of fluid model in the presence of a magnetic field

NASA Astrophysics Data System (ADS)

Shafiei Dizaji, A.; Mohammadpourfard, M.; Aminfar, H.

2018-03-01

Multiphase flow is one of the most complicated problems, considering the multiplicity of the related parameters, especially the external factors influences. Thus, despite the recent developments more investigations are still required. The effect of a uniform magnetic field on the hydrodynamics behavior of a two-phase flow with different magnetic permeability is presented in this article. A single water vapor bubble which is rising inside a channel filled with ferrofluid has been simulated numerically. To capture the phases interface, the Volume of Fluid (VOF) model, and to solve the governing equations, the finite volume method has been employed. Contrary to the prior anticipations, while the consisting fluids of the flow are dielectric, uniform magnetic field causes a force acting normal to the interface toward to the inside of the bubble. With respect to the applied magnetic field direction, the bubble deformation due to the magnetic force increases the bubble rising velocity. Moreover, the higher values of applied magnetic field strength and magnetic permeability ratio resulted in the further increase of the bubble rising velocity. Also it is indicated that the flow mixing and the heat transfer rate is increased by a bubble injection and applying a magnetic field. The obtained results have been concluded that the presented phenomenon with applying a magnetic field can be used to control the related characteristics of the multiphase flows. Compared to the previous studies, implementing the applicable cases using the common and actual materials and a significant reduction of the CPU time are the most remarkable advantages of the current study.

μ-PIV measurements of the ensemble flow fields surrounding a migrating semi-infinite bubble.

PubMed

Yamaguchi, Eiichiro; Smith, Bradford J; Gaver, Donald P

2009-08-01

Microscale particle image velocimetry (μ-PIV) measurements of ensemble flow fields surrounding a steadily-migrating semi-infinite bubble through the novel adaptation of a computer controlled linear motor flow control system. The system was programmed to generate a square wave velocity input in order to produce accurate constant bubble propagation repeatedly and effectively through a fused glass capillary tube. We present a novel technique for re-positioning of the coordinate axis to the bubble tip frame of reference in each instantaneous field through the analysis of the sudden change of standard deviation of centerline velocity profiles across the bubble interface. Ensemble averages were then computed in this bubble tip frame of reference. Combined fluid systems of water/air, glycerol/air, and glycerol/Si-oil were used to investigate flows comparable to computational simulations described in Smith and Gaver (2008) and to past experimental observations of interfacial shape. Fluorescent particle images were also analyzed to measure the residual film thickness trailing behind the bubble. The flow fields and film thickness agree very well with the computational simulations as well as existing experimental and analytical results. Particle accumulation and migration associated with the flow patterns near the bubble tip after long experimental durations are discussed as potential sources of error in the experimental method.

Numerical simulation of superheated vapor bubble rising in stagnant liquid

NASA Astrophysics Data System (ADS)

Samkhaniani, N.; Ansari, M. R.

2017-09-01

In present study, the rising of superheated vapor bubble in saturated liquid is simulated using volume of fluid method in OpenFOAM cfd package. The surface tension between vapor-liquid phases is considered using continuous surface force method. In order to reduce spurious current near interface, Lafaurie smoothing filter is applied to improve curvature calculation. Phase change is considered using Tanasawa mass transfer model. The variation of saturation temperature in vapor bubble with local pressure is considered with simplified Clausius-Clapeyron relation. The couple velocity-pressure equation is solved using PISO algorithm. The numerical model is validated with: (1) isothermal bubble rising and (2) one-dimensional horizontal film condensation. Then, the shape and life time history of single superheated vapor bubble are investigated. The present numerical study shows vapor bubble in saturated liquid undergoes boiling and condensation. It indicates bubble life time is nearly linear proportional with bubble size and superheat temperature.

Freeze/Thaw-Induced Embolism: Probability of Critical Bubble Formation Depends on Speed of Ice Formation

DOE PAGES

Sevanto, Sanna; Holbrook, N. Michele; Ball, Marilyn C.

2012-06-06

Bubble formation in the conduits of woody plants sets a challenge for uninterrupted water transportation from the soil up to the canopy. Freezing and thawing of stems has been shown to increase the number of air-filled (embolized) conduits, especially in trees with large conduit diameters. Despite numerous experimental studies, the mechanisms leading to bubble formation during freezing have not been addressed theoretically. We used classical nucleation theory and fluid mechanics to show which mechanisms are most likely to be responsible for bubble formation during freezing and what parameters determine the likelihood of the process. Our results confirm the common assumptionmore » that bubble formation during freezing is most likely due to gas segregation by ice. If xylem conduit walls are not permeable to the salts expelled by ice during the freezing process, osmotic pressures high enough for air seeding could be created. The build-up rate of segregated solutes in front of the ice-water interface depends equally on conduit diameter and freezing velocity. Therefore, bubble formation probability depends on these variables. The dependence of bubble formation probability on freezing velocity means that the experimental results obtained for cavitation threshold conduit diameters during freeze/thaw cycles depend on the experimental setup; namely sample size and cooling rate. The velocity dependence also suggests that to avoid bubble formation during freezing trees should have narrow conduits where freezing is likely to be fast (e.g., branches or outermost layer of the xylem). Avoidance of bubble formation during freezing could thus be one piece of the explanation why xylem conduit size of temperate and boreal zone trees varies quite systematically.« less

Freeze/Thaw-induced embolism: probability of critical bubble formation depends on speed of ice formation.

PubMed

Sevanto, Sanna; Holbrook, N Michele; Ball, Marilyn C

2012-01-01

Bubble formation in the conduits of woody plants sets a challenge for uninterrupted water transportation from the soil up to the canopy. Freezing and thawing of stems has been shown to increase the number of air-filled (embolized) conduits, especially in trees with large conduit diameters. Despite numerous experimental studies, the mechanisms leading to bubble formation during freezing have not been addressed theoretically. We used classical nucleation theory and fluid mechanics to show which mechanisms are most likely to be responsible for bubble formation during freezing and what parameters determine the likelihood of the process. Our results confirm the common assumption that bubble formation during freezing is most likely due to gas segregation by ice. If xylem conduit walls are not permeable to the salts expelled by ice during the freezing process, osmotic pressures high enough for air seeding could be created. The build-up rate of segregated solutes in front of the ice-water interface depends equally on conduit diameter and freezing velocity. Therefore, bubble formation probability depends on these variables. The dependence of bubble formation probability on freezing velocity means that the experimental results obtained for cavitation threshold conduit diameters during freeze/thaw cycles depend on the experimental setup; namely sample size and cooling rate. The velocity dependence also suggests that to avoid bubble formation during freezing trees should have narrow conduits where freezing is likely to be fast (e.g., branches or outermost layer of the xylem). Avoidance of bubble formation during freezing could thus be one piece of the explanation why xylem conduit size of temperate and boreal zone trees varies quite systematically.

Diffuse interface simulation of bubble rising process: a comparison of adaptive mesh refinement and arbitrary lagrange-euler methods

NASA Astrophysics Data System (ADS)

Wang, Ye; Cai, Jiejin; Li, Qiong; Yin, Huaqiang; Yang, Xingtuan

2018-06-01

Gas-liquid two phase flow exists in several industrial processes and light-water reactors (LWRs). A diffuse interface based finite element method with two different mesh generation methods namely, the Adaptive Mesh Refinement (AMR) and the Arbitrary Lagrange Euler (ALE) methods is used to model the shape and velocity changes in a rising bubble. Moreover, the calculating speed and mesh generation strategies of AMR and ALE are contrasted. The simulation results agree with the Bhagat's experiments, indicating that both mesh generation methods can simulate the characteristics of bubble accurately. We concluded that: the small bubble rises as elliptical with oscillation, whereas a larger bubble (11 mm > d > 7 mm) rises with a morphology between the elliptical and cap type with a larger oscillation. When the bubble is large (d > 11 mm), it rises up as a cap type, and the amplitude becomes smaller. Moreover, it takes longer to achieve the stable shape from the ellipsoid to the spherical cap type with the increase of the bubble diameter. The results also show that for smaller diameter case, the ALE method uses fewer grids and has a faster calculation speed, but the AMR method can solve the case of a large geometry deformation efficiently.

Buoyancy Driven Shear Flows of Bubble Suspensions

NASA Technical Reports Server (NTRS)

Koch, D. L.; Hill, R. J.; Chellppannair, T.; Zenit, R.; Zenit, R.; Spelt, P. D. M.

1999-01-01

In this work the gas volume fraction and the root-mean-squared fluid velocity are measured in buoyancy driven shear flows of bubble suspensions in a tall, inclined, rectangular channel. The experiments are performed under conditions where We << 1a nd Re >> 1, for which comparisons are made with kinetic theory and numerical simulations. Here Re = gamma(a(exp 2)/nu is the Reynolds number and We = rho(gamma(exp 2))a(exp 3)/sigma is the Weber number; gamma is the shear rate, a is the bubble radius, nu is the kinematic viscosity of the liquid, rho is the density of the liquid, and sigma is the surface tension of the gas/liquid interface. Kang et al. calculated the bubble phase pressure and velocity variance of sheared bubble suspensions under conditions where the bubbles are spherical and the liquid phase velocity field can be approximated using potential flow theory, i.e. We= 0 and Re >> 1. Such conditions can be achieved in an experiment using gas bubbles, with a radius of O(0.5mm), in water. The theory requires that there be no average relative motion of the gas and liquid phases, hence the motivation for an experimental program in microgravity. The necessity of performing preliminary, Earth based experiments, however, requires performing experiments where the gas phase rises in the liquid, which significantly complicates the comparison of experiments with theory. Rather than comparing experimental results with theory for a uniform, homogeneous shear flow, experiments can be compared directly with solutions of the averaged equations of motion for bubble suspensions. This requires accounting for the significant lift force acting on the gas phase when the bubbles rise parallel to the average velocity of the sheared suspension. Shear flows can be produced in which the bubble phase pressure gradient, arising from shear induced collisions amongst the bubbles, balances a body force (centrifugal or gravitational) on the gas phase. A steady, non-uniform gas volume fraction can be measured, from which the bubble phase pressure gradient can be obtained and compared to theory and numerical simulations. The presence of bounding walls further complicates the experiments, since the detailed interactions of the bubbles with bounding walls is not well understood, especially in the presence of gravity, where the momentum and energy exchange depends on the inclination of the wall.

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.

Simulations of viscous and compressible gas-gas flows using high-order finite difference schemes

NASA Astrophysics Data System (ADS)

Capuano, M.; Bogey, C.; Spelt, P. D. M.

2018-05-01

A computational method for the simulation of viscous and compressible gas-gas flows is presented. It consists in solving the Navier-Stokes equations associated with a convection equation governing the motion of the interface between two gases using high-order finite-difference schemes. A discontinuity-capturing methodology based on sensors and a spatial filter enables capturing shock waves and deformable interfaces. One-dimensional test cases are performed as validation and to justify choices in the numerical method. The results compare well with analytical solutions. Shock waves and interfaces are accurately propagated, and remain sharp. Subsequently, two-dimensional flows are considered including viscosity and thermal conductivity. In Richtmyer-Meshkov instability, generated on an air-SF6 interface, the influence of the mesh refinement on the instability shape is studied, and the temporal variations of the instability amplitude is compared with experimental data. Finally, for a plane shock wave propagating in air and impacting a cylindrical bubble filled with helium or R22, numerical Schlieren pictures obtained using different grid refinements are found to compare well with experimental shadow-photographs. The mass conservation is verified from the temporal variations of the mass of the bubble. The mean velocities of pressure waves and bubble interface are similar to those obtained experimentally.

The buoyancy-driven motion of a single skirted bubble or drop rising through a viscous liquid

NASA Astrophysics Data System (ADS)

Ohta, Mitsuhiro; Sussman, Mark

2012-11-01

The buoyancy-driven motion of a single skirted bubble or drop rising through a viscous liquid is computationally explored by way of 3d-axisymmetric computations. The Navier-Stokes equations for incompressible two-fluid flow are solved numerically in which the coupled level-set and volume-of-fluid method is used to simulate the deforming bubble/drop boundary and the interface jump conditions on the deforming boundary are enforced through a sharp interface numerical treatment. Dynamic, block structured adaptive grid refinement is employed in order to sufficiently resolve the thin skirts. Results on the sensitivity of the thickness of trailing bubble/drop skirts to the density ratio and viscosity ratio are reported. It is shown that both the density ratio (not the density difference) and the viscosity ratio effect the skirt thickness. Previous theory for predicting skirt thickness can be refined as a result of our calculations. It is also discovered that the formation of thin skirts for bubbles and drops have little effect on the rise velocity. In other words, the measured Re number for cases without skirt formation have almost the same values for Re as cases with a thin skirt.

Experimental study of the flow pattern around a bubble confined in a microfluidic Hele-Shaw cell

NASA Astrophysics Data System (ADS)

Tsoumpas, Yannis; Fajolles, Christophe; Malloggi, Florent

2017-11-01

The flow field around a bubble moving with respect to a surrounding liquid in a Hele-Shaw cell can usually be characterized by a recirculating flow, which is typically attributed to a Marangoni effect due to surface tension gradients generated by a non-uniform distribution of surfactants (or temperature) along the liquid-gas interface. In the present study, we try to visualize such a flow employing 3D micro-particle tracking velocimetry. We perform experiments on an immobile flattened air bubble that is surrounded by a flow of aqueous solution of surfactant (SDS), in a microfluidic chamber described in the work of Sungyon Lee et al.. The suspending fluid is seeded with spherical micro-particles, with those captured by the recirculating flow orbiting in a three-dimensional trajectory in the vicinity of the liquid-air interface. We address the effect of velocity of the surrounding fluid, surfactant concentration and bubble radius on the recirculating flow pattern. The case of a liquid-liquid interface, with a hexadecane drop as the dispersed phase, is also discussed. The authors would like to acknowledge the financial support of Enhanced Eurotalents program (an FP7 Marie Skłodowska-Curie COFUND program) & ANR (ANR-13-BS09-0011).

Direct numerical simulation of incompressible multiphase flow with phase change

NASA Astrophysics Data System (ADS)

Lee, Moon Soo; Riaz, Amir; Aute, Vikrant

2017-09-01

Simulation of multiphase flow with phase change is challenging because of the potential for unphysical pressure oscillations, spurious velocity fields and mass flux errors across the interface. The resulting numerical errors may become critical when large density contrasts are present. To address these issues, we present a new approach for multiphase flow with phase change that features, (i) a smooth distribution of sharp velocity jumps and mass flux within a narrow region surrounding the interface, (ii) improved mass flux projection from the implicit interface onto the uniform Cartesian grid and (iii) post-advection velocity correction step to ensure accurate velocity divergence in interfacial cells. These new features are implemented in combination with a sharp treatment of the jumps in pressure and temperature gradient. A series of 1-D, 2-D, axisymmetric and 3-D problems are solved to verify the improvements afforded by the new approach. Axisymmetric film boiling results are also presented, which show good qualitative agreement with heat transfer correlations as well as experimental observations of bubble shapes.

Effect of surface mobility on the particle sliding along a bubble or a solid sphere.

PubMed

Wang, Weixing; Zhou, Zhiang; Nandakumar, K; Xu, Zhenghe; Masliyah, Jacob H

2003-03-01

The sliding velocity of glass beads on a spherical surface, made either of an air bubble or of a glass sphere held stationary, is measured to investigate the effect of surface mobility on the particle sliding velocity. The sliding process is recorded with a digital camera and analyzed frame by frame. The sliding glass bead was found to accelerate with increasing angular position on the collector's surface. It reaches a maximum velocity at an angular position of about 100 degrees and then, under certain conditions, the glass bead leaves the surface of the collector. The sliding velocity of the glass bead depends strongly on the surface mobility of a bubble, decreasing with decreasing surface mobility. By a mobile surface we mean one which cannot set up resistive forces to an applied stress on the surface. The sliding velocity on a rigid surface, such as a glass sphere, is much lower than that on a mobile bubble surface. The sliding velocity can be described through a modified Stokes equation. A numerical factor in the modified Stokes equation is determined by fitting the experimental data and is found to increase with decreasing surface mobility. Hydrophobic glass beads sliding on a hydrophobic glass sphere were found to stick at the point of impact without sliding if the initial angular position of the impact is less than some specific angle, which is defined as the critical sticking angle. The sticking of the glass beads can be attributed to the capillary contracting force created by the formation of a cavity due to spontaneous receding of the nonwetting liquid from the contact zone. The relationship between the critical sticking angle and the particle size is established based on the Yushchenko [J. Colloid Interface Sci. 96 (1983) 307] analysis.

Limits of the potential flow approach to the single-mode Rayleigh-Taylor problem

NASA Astrophysics Data System (ADS)

Ramaprabhu, P.; Dimonte, Guy; Young, Yuan-Nan; Calder, A. C.; Fryxell, B.

2006-12-01

We report on the behavior of a single-wavelength Rayleigh-Taylor flow at late times. The calculations were performed in a long square duct (λ×λ×8λ) , using four different numerical simulations. In contradiction with potential flow theories that predict a constant terminal velocity, the single-wavelength Rayleigh-Taylor problem exhibits late-time acceleration. The onset of acceleration occurs as the bubble penetration depth exceeds the diameter of bubbles, and is observed for low and moderate density differences. Based on our simulations, we provide a phenomenological description of the observed acceleration, and ascribe this behavior to the formation of Kelvin-Helmholtz vortices on the bubble-spike interface that diminish the friction drag, while the associated induced flow propels the bubbles forward. For large density ratios, the formation of secondary instabilities is suppressed, and the bubbles remain terminal consistent with potential flow models.

Numerical simulations of non-spherical bubble collapse.

PubMed

Johnsen, Eric; Colonius, Tim

2009-06-01

A high-order accurate shock- and interface-capturing scheme is used to simulate the collapse of a gas bubble in water. In order to better understand the damage caused by collapsing bubbles, the dynamics of the shock-induced and Rayleigh collapse of a bubble near a planar rigid surface and in a free field are analysed. Collapse times, bubble displacements, interfacial velocities and surface pressures are quantified as a function of the pressure ratio driving the collapse and of the initial bubble stand-off distance from the wall; these quantities are compared to the available theory and experiments and show good agreement with the data for both the bubble dynamics and the propagation of the shock emitted upon the collapse. Non-spherical collapse involves the formation of a re-entrant jet directed towards the wall or in the direction of propagation of the incoming shock. In shock-induced collapse, very high jet velocities can be achieved, and the finite time for shock propagation through the bubble may be non-negligible compared to the collapse time for the pressure ratios of interest. Several types of shock waves are generated during the collapse, including precursor and water-hammer shocks that arise from the re-entrant jet formation and its impact upon the distal side of the bubble, respectively. The water-hammer shock can generate very high pressures on the wall, far exceeding those from the incident shock. The potential damage to the neighbouring surface is quantified by measuring the wall pressure. The range of stand-off distances and the surface area for which amplification of the incident shock due to bubble collapse occurs is determined.

Numerical simulations of non-spherical bubble collapse

PubMed Central

JOHNSEN, ERIC; COLONIUS, TIM

2009-01-01

A high-order accurate shock- and interface-capturing scheme is used to simulate the collapse of a gas bubble in water. In order to better understand the damage caused by collapsing bubbles, the dynamics of the shock-induced and Rayleigh collapse of a bubble near a planar rigid surface and in a free field are analysed. Collapse times, bubble displacements, interfacial velocities and surface pressures are quantified as a function of the pressure ratio driving the collapse and of the initial bubble stand-off distance from the wall; these quantities are compared to the available theory and experiments and show good agreement with the data for both the bubble dynamics and the propagation of the shock emitted upon the collapse. Non-spherical collapse involves the formation of a re-entrant jet directed towards the wall or in the direction of propagation of the incoming shock. In shock-induced collapse, very high jet velocities can be achieved, and the finite time for shock propagation through the bubble may be non-negligible compared to the collapse time for the pressure ratios of interest. Several types of shock waves are generated during the collapse, including precursor and water-hammer shocks that arise from the re-entrant jet formation and its impact upon the distal side of the bubble, respectively. The water-hammer shock can generate very high pressures on the wall, far exceeding those from the incident shock. The potential damage to the neighbouring surface is quantified by measuring the wall pressure. The range of stand-off distances and the surface area for which amplification of the incident shock due to bubble collapse occurs is determined. PMID:19756233

Evaluation of Interfacial Forces and Bubble-Induced Turbulence Using Direct Numerical Simulation

NASA Astrophysics Data System (ADS)

Feng, Jinyong

High fidelity prediction of multiphase flows is important in a wide range of engineering applications. While some multiphase flow scenarios can be successfully modeled, many questions remain unanswered regarding the interaction between the bubbles and the turbulence, and present significant challenges in the development of closure laws for the multiphase computational fluid dynamics (M-CFD) models. To address these challenges, we propose to evaluate the interfacial forces and bubble-induced turbulence in both laminar and turbulent flow field with direct numerical simulation (DNS) approach. Advanced finite-element based flow solver (PHASTA) with level-set interface tracking method is utilized for these studies. The proportional-integral-derivative (PID) controller is adopted to ensure the statistically steady state bubble position and perform the detailed study of the turbulent field around the bubble. Selected numerical capabilities and post-processing codes are developed to achieve the research goals. The interface tracking approach is verified and validated by comparing the interfacial forces with the experiment-based data and correlations. The sign change of transverse lift force is observed as the bubble becomes more deformable. A new correlation is proposed to predict the behavior of the drag coefficient over the wide range of conditions. The wall effect on the interfacial forces are also investigated. In homogeneous turbulent flow, the effect of bubble deformability, turbulent intensity and relative velocity on the bubble-induced turbulence are analyzed. The presented method and novel results will complement the experimental database, provide insight to the bubbleinduced turbulence mechanism and help the development of M-CFD closure models.

The Effect of the Density Ratio on the Nonlinear Dynamics of the Unstable Fluid Interface

NASA Technical Reports Server (NTRS)

Abarzhi, S. I.

2003-01-01

Here we report multiple harmonic theoretical solutions for a complete system of conservation laws, which describe the large-scale coherent dynamics in RTI and RMI for fluids with a finite density ratio in the general three-dimensional case. The analysis yields new properties of the bubble front dynamics. In either RTI or RMI, the obtained dependencies of the bubble velocity and curvature on the density ratio differ qualitatively and quantitatively from those suggested by the models of Sharp (1984), Oron et al. (2001), and Goncharov (2002). We show explicitly that these models violate the conservation laws. For the first time, our theory reveals an important qualitative distinction between the dynamics of the RT and RM bubbles.

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.

Investigation of the interaction dynamics of a pair of laser-induced bubbles generated at the same time through double-exposure strobe method and numerical simulations

NASA Astrophysics Data System (ADS)

Han, Bing; Liu, Liu; Ni, Xiao-Wu

2017-08-01

In order to understand the interaction dynamics of a pair of laser-induced bubbles, a double-exposure strobe photography experimental setup is build up to study the temporal evolution of the bubble pairs and to measure the transient bubble-interface moving speed. The interaction mechanisms of the bubble pairs are discussed together with the numerical results obtained through OpenFOAM. It is shown that the direction and the velocity of the jetting could be controlled by the relative size and the relative initiation distance of the bubble pair, when the bubbles are generated at the same time, i.e., in-phase. The liquid jet is considered to be a penetrating jet. The jet is originated from the smaller bubble and clearly protruding outside of the bigger bubble. The parameter space of the relative size and the initiation distance of the bubble pair allowing the formation of the penetrating jet are very narrow. It is concluded that the liquid jet induced by the bubble interactions resulted from the collapse and the rebound of the smaller bubble nearby the bigger bubble. This is defined as the "catapult effect." Such a directional liquid transportation is a promising tool as a micro-injector or a micro-pump. The investigation results could be also supplementary to the understandings of the bubble dynamics.

The Seismic Velocity In Gas-charged Magma

NASA Astrophysics Data System (ADS)

Sturton, S.; Neuberg, J. W.

2001-12-01

Long-period and hybrid events, seen at the Soufrière Hills Volcano, Montserrat, show dominant low frequency content suggesting the seismic wavefield is formed as a result of interface waves at the boundary between a fluid and a solid medium. This wavefield will depend on the impedance contrast between the two media and therefore the difference in seismic velocity. For a gas-charged magma, increasing pressure with depth reduces the volume of gas exsolved, increasing the seismic velocity with depth in the conduit. The seismic radiation pattern along the conduit can then be modelled. Where single events merge into tremor, gliding lines can sometimes be seen in the spectra and indicate either changes in the seismic parameters with time or varying triggering rates of single events.The differential equation describing the time dependence of bubble growth by diffusion is solved numerically for a stationary magma column undergoing a decompression event. The volume of gas is depth dependent and increases with time as the bubbles grow and expand. It is used to calculate the depth and time dependence of the density, pressure and seismic velocity. The effect of different viscosities associated with different magma types and concentration of water in the melt on the rate of bubble growth is explored. Crystal growth, which increases the concentration of water in the melt, affects the amount of gas that can be exsolved.

Dynamics of Rear Stagnant Cap formation at the surface of spherical bubbles rising in surfactant solutions at large Reynolds numbers under conditions of small Marangoni number and slow sorption kinetics.

PubMed

Dukhin, S S; Kovalchuk, V I; Gochev, G G; Lotfi, M; Krzan, M; Malysa, K; Miller, R

2015-08-01

On the surface of bubbles rising in a surfactant solution the adsorption process proceeds and leads to the formation of a so called Rear Stagnant Cap (RSC). The larger this RSC is the stronger is the retardation of the rising velocity. The theory of a steady RSC and steady retarded rising velocity, which sets in after a transient stage, has been generally accepted. However, a non-steady process of bubble rising starting from the initial zero velocity represents an important portion of the trajectory of rising, characterized by a local velocity profile (LVP). As there is no theory of RSC growth for large Reynolds numbers Re » 1 so far, the interpretation of LVPs measured in this regime was impossible. It turned out, that an analytical theory for a quasi-steady growth of RSC is possible for small Marangoni numbers Ma « 1, i.e. when the RSC is almost completely compressed, which means a uniform surface concentration Γ(θ)=Γ(∞) within the RSC. Hence, the RSC angle ψ(t) is obtained as a function of the adsorption isotherm parameters and time t. From the steady velocity v(st)(ψ), the dependence of non-steady velocity on time is obtained by employing v(st)[ψ(t)] via a quasi-steady approximation. The measurement of LVP creates a promising new opportunity for investigation of the RSC dynamics and adsorption kinetics. While adsorption and desorption happen at the same localization in the classical methods, in rising bubble experiments desorption occurs mainly within RSC while adsorption on the mobile part of the bubble surface. The desorption flux from RSC is proportional to αΓ(∞), while it is usually αΓ. The adsorption flux at the mobile surface above RSC can be assumed proportional to βC0, while it is usually βC0(1-Γ/Γ(∞)). These simplifications may become favorable in investigations of the adsorption kinetics for larger molecules, in particular for globular proteins, which essentially stay at an interface once adsorbed. Copyright © 2014 Elsevier B.V. All rights reserved.

Chemical reactions at aqueous interfaces

NASA Astrophysics Data System (ADS)

Vecitis, Chad David

2009-12-01

Interfaces or phase boundaries are a unique chemical environment relative to individual gas, liquid, or solid phases. Interfacial reaction mechanisms and kinetics are often at variance with homogeneous chemistry due to mass transfer, molecular orientation, and catalytic effects. Aqueous interfaces are a common subject of environmental science and engineering research, and three environmentally relevant aqueous interfaces are investigated in this thesis: 1) fluorochemical sonochemistry (bubble-water), 2) aqueous aerosol ozonation (gas-water droplet), and 3) electrolytic hydrogen production and simultaneous organic oxidation (water-metal/semiconductor). Direct interfacial analysis under environmentally relevant conditions is difficult, since most surface-specific techniques require relatively `extreme' conditions. Thus, the experimental investigations here focus on the development of chemical reactors and analytical techniques for the completion of time/concentration-dependent measurements of reactants and their products. Kinetic modeling, estimations, and/or correlations were used to extract information on interfacially relevant processes. We found that interfacial chemistry was determined to be the rate-limiting step to a subsequent series of relatively fast homogeneous reactions, for example: 1) Pyrolytic cleavage of the ionic headgroup of perfluorooctanesulfonate (PFOS) and perfluorooctanoate (PFOA) adsorbed to cavitating bubble-water interfaces during sonolysis was the rate-determining step in transformation to their inorganic constituents carbon monoxide, carbon dioxide, and fluoride; 2) ozone oxidation of aqueous iodide to hypoiodous acid at the aerosol-gas interface is the rate-determining step in the oxidation of bromide and chloride to dihalogens; 3) Electrolytic oxidation of anodic titanol surface groups is rate-limiting for the overall oxidation of organics by the dichloride radical. We also found chemistry unique to the interface, for example: 1) Adsorption of dilute PFOS(aq) and PFOA(aq) to acoustically cavitating bubble interfaces was greater than equilibrium expectations due to high-velocity bubble radial oscillations; 2) Relative ozone oxidation kinetics of aqueous iodide, sulfite, and thiosulfate were at variance with previously reported bulk aqueous kinetics; 3) Organics that directly chelated with the anode surface were oxidized by direct electron transfer, resulting in immediate carbon dioxide production but slower overall oxidation kinetics. Chemical reactions at aqueous interfaces can be the rate-limiting step of a reaction network and often display novel mechanisms and kinetics as compared to homogeneous chemistry.

Bubble velocity, diameter, and void fraction measurements in a multiphase flow using fiber optic reflectometer

NASA Astrophysics Data System (ADS)

Lim, Ho-Joon; Chang, Kuang-An; Su, Chin B.; Chen, Chi-Yueh

2008-12-01

A fiber optic reflectometer (FOR) technique featuring a single fiber probe is investigated for its feasibility of measuring the bubble velocity, diameter, and void fraction in a multiphase flow. The method is based on the interference of the scattered signal from the bubble surface with the Fresnel reflection signal from the tip of the optical fiber. Void fraction is obtained with a high accuracy if an appropriate correction is applied to compensate the underestimated measurement value. Velocity information is accurately obtained from the reflected signals before the fiber tip touches the bubble surface so that several factors affecting the traditional dual-tip probes such as blinding, crawling, and drifting effects due to the interaction between the probe and bubbles can be prevented. The coherent signals reflected from both the front and rear ends of a bubble can provide velocity information. Deceleration of rising bubbles and particles due to the presence of the fiber probe is observed when they are very close to the fiber tip. With the residence time obtained, the bubble chord length can be determined by analyzing the coherent signal for velocity determination before the deceleration starts. The bubble diameters are directly obtained from analyzing the signals of the bubbles that contain velocity information. The chord lengths of these bubbles measured by FOR represent the bubble diameters when the bubble shape is spherical or represent the minor axes when the bubble shape is ellipsoidal. The velocity and size of bubbles obtained from the FOR measurements are compared with those obtained simultaneously using a high speed camera.

Microjet formation in a capillary by laser-induced cavitation

NASA Astrophysics Data System (ADS)

Peters, Ivo R.; Tagawa, Yoshiyuki; van der Meer, Devaraj; Prosperetti, Andrea; Sun, Chao; Lohse, Detlef

2010-11-01

A vapor bubble is created by focusing a laser pulse inside a capillary that is partially filled with water. Upon creation of the bubble, a shock wave travels through the capillary. When this shock wave meets the meniscus of the air-water interface, a thin jet is created that travels at very high speeds. A crucial ingredient for the creation of the jet is the shape of the meniscus, which is responsible for focusing the energy provided by the shock wave. We examine the formation of this jet numerically using a boundary integral method, where we prepare an initial interface at rest inside a tube with a diameter ranging from 50 to 500 μm. To simulate the effect of the bubble we then apply a short, strong pressure pulse, after which the jet forms. We investigate the influence of the shape of the meniscus, and pressure amplitude and duration on the jet formation. The jet shape and velocity obtained by the simulation compare well with experimental data, and provides good insight in the origin of the jet.

Taylor bubbles at high viscosity ratios: experiments and numerical simulations

NASA Astrophysics Data System (ADS)

Hewakandamby, Buddhika; Hasan, Abbas; Azzopardi, Barry; Xie, Zhihua; Pain, Chris; Matar, Omar

2015-11-01

The Taylor bubble is a single long bubble which nearly fills the entire cross section of a liquid-filled circular tube, often occurring in gas-liquid slug flows in many industrial applications, particularly oil and gas production. The objective of this study is to investigate the fluid dynamics of three-dimensional Taylor bubble rising in highly viscous silicone oil in a vertical pipe. An adaptive unstructured mesh modelling framework is adopted here which can modify and adapt anisotropic unstructured meshes to better represent the underlying physics of bubble rising and reduce computational effort without sacrificing accuracy. The numerical framework consists of a mixed control volume and finite element formulation, a `volume of fluid'-type method for the interface-capturing based on a compressive control volume advection method, and a force-balanced algorithm for the surface tension implementation. Experimental results for the Taylor bubble shape and rise velocity are presented, together with numerical results for the dynamics of the bubbles. A comparison of the simulation predictions with experimental data available in the literature is also presented to demonstrate the capabilities of our numerical method. EPSRC Programme Grant, MEMPHIS, EP/K0039761/1.

The interaction of bubbles with solidification interfaces. [during coasting phase of sounding rocket flight

NASA Technical Reports Server (NTRS)

Papazian, J. M.; Wilcox, W. R.

1977-01-01

The behavior of bubbles at a dendritic solidification interface was studied during the coasting phase of a sounding rocket flight. Sequential photographs of the gradient freeze experiment showed nucleation, growth and coalescence of bubbles at the moving interface during both the low-gravity and one-gravity tests. In the one-gravity test the bubbles were observed to detach from the interface and float to the top of the melt. However, in the low-gravity tests no bubble detachment from the interface or steady state bubble motion occurred and large voids were grown into the crystal. These observations are discussed in terms of the current theory of thermal migration of bubbles and in terms of their implications on the space processing of metals.

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.

Weakly nonlinear incompressible Rayleigh-Taylor instability growth at cylindrically convergent interfaces

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

Wang, L. F.; He, X. T.; HEDPS, Center for Applied Physics and Technology, Peking University, Beijing 100871

2013-04-15

A weakly nonlinear (WN) model has been developed for the incompressible Rayleigh-Taylor instability (RTI) in cylindrical geometry. The transition from linear to nonlinear growth is analytically investigated via a third-order solutions for the cylindrical RTI initiated by a single-mode velocity perturbation. The third-order solutions can depict the early stage of the interface asymmetry due to the bubble-spike formation, as well as the saturation of the linear (exponential) growth of the fundamental mode. The WN results in planar RTI [Wang et al., Phys. Plasmas 19, 112706 (2012)] are recovered in the limit of high-mode number perturbations. The difference between the WNmore » growth of the RTI in cylindrical geometry and in planar geometry is discussed. It is found that the interface of the inward (outward) development spike/bubble is extruded (stretched) by the additional inertial force in cylindrical geometry compared with that in planar geometry. For interfaces with small density ratios, the inward growth bubble can grow fast than the outward growth spike in cylindrical RTI. Moreover, a reduced formula is proposed to describe the WN growth of the RTI in cylindrical geometry with an acceptable precision, especially for small-amplitude perturbations. Using the reduced formula, the nonlinear saturation amplitude of the fundamental mode and the phases of the Fourier harmonics are studied. Thus, it should be included in applications where converging geometry effects play an important role, such as the supernova explosions and inertial confinement fusion implosions.« less

How are soap bubbles blown? Fluid dynamics of soap bubble blowing

NASA Astrophysics Data System (ADS)

Davidson, John; Lambert, Lori; Sherman, Erica; Wei, Timothy; Ryu, Sangjin

2013-11-01

Soap bubbles are a common interfacial fluid dynamics phenomenon having a long history of delighting not only children and artists but also scientists. In contrast to the dynamics of liquid droplets in gas and gas bubbles in liquid, the dynamics of soap bubbles has not been well documented. This is possibly because studying soap bubbles is more challenging due to there existing two gas-liquid interfaces. Having the thin-film interface seems to alter the characteristics of the bubble/drop creation process since the interface has limiting factors such as thickness. Thus, the main objective of this study is to determine how the thin-film interface differentiates soap bubbles from gas bubbles and liquid drops. To investigate the creation process of soap bubbles, we constructed an experimental model consisting of air jet flow and a soap film, which consistently replicates the conditions that a human produces when blowing soap bubbles, and examined the interaction between the jet and the soap film using the high-speed videography and the particle image velocimetry.

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.

Expanding Taylor bubble under constant heat flux

NASA Astrophysics Data System (ADS)

Voirand, Antoine; Benselama, Adel M.; Ayel, Vincent; Bertin, Yves

2016-09-01

Modelization of non-isothermal bubbles expanding in a capillary, as a contribution to the understanding of the physical phenomena taking place in Pulsating Heat Pipes (PHPs), is the scope of this paper. The liquid film problem is simplified and solved, while the thermal problem takes into account a constant heat flux density applied at the capillary tube wall, exchanging with the liquid film surrounding the bubble and also with the capillary tube outside medium. The liquid slug dynamics is solved using the Lucas-Washburn equation. Mass and energy balance on the vapor phase allow governing equations of bubble expansion to be written. The liquid and vapor phases are coupled only through the saturation temperature associated with the vapor pressure, assumed to be uniform throughout the bubble. Results show an over-heating of the vapor phase, although the particular thermal boundary condition used here always ensures an evaporative mass flux at the liquid-vapor interface. Global heat exchange is also investigated, showing a strong decreasing of the PHP performance to convey heat by phase change means for large meniscus velocities.

Jet and Vortex Projectile Flows in Shock/bubble-on-wall Configuration

NASA Astrophysics Data System (ADS)

Peng, Gaozhu; Zabusky, Norman

2001-11-01

We observe intense coaxial upstream and radial flow structures from a shock in air interacting with a SF6 half-bubble placed against an ideally reflecting wall. Our axisymmetric numerical simulations were done with PPM and models a spherical bubble struck symmetrically by two identical approaching shocks . A "dual" vorticity deposition arises at early time and a coaxial upstream moving primary jet and radial vortex ring flow appears. A coherent vortex ring or vortex projectile (VP), with entrained shocklets originates from the vortex layer produced at the Mach stem (which arises from the primary reflected shock). This VP moves ahead of the jet. The original transmitted wave and other trapped waves in the expanding axial jet causes a collapsing and expanding cavity and other instabilities on the complex bubble interface. We present and analyze our results with different diagnostics: vorticity, density, divergence of velocity, and numerical shadowgraph patterns; global quantification of circulation, enstrophy and r-integrated vorticity; etc. We also discuss data projection and filtering for quantifying and validating complex flows.

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.

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

Spontaneous and Directional Bubble Transport on Porous Copper Wires with Complex Shapes in Aqueous Media.

PubMed

Li, Wenjing; Zhang, Jingjing; Xue, Zhongxin; Wang, Jingming; Jiang, Lei

2018-01-24

Manipulation of gas bubble behaviors is crucial for gas bubble-related applications. Generally, the manipulation of gas bubble behaviors generally takes advantage of their buoyancy force. It is very difficult to control the transportation of gas bubbles in a specific direction. Several approaches have been developed to collect and transport bubbles in aqueous media; however, most reliable and effective manipulation of gas bubbles in aqueous media occurs on the interfaces with simple shapes (i.e., cylinder and cone shapes). Reliable strategies for spontaneous and directional transport of gas bubbles on interfaces with complex shapes remain enormously challenging. Herein, a type of 3D gradient porous network was constructed on copper wire interfaces, with rectangle, wave, and helix shapes. The superhydrophobic copper wires were immersed in water, and continuous and stable gas films then formed on the interfaces. With the assistance of the Laplace pressure gradient between two bubbles, gas bubbles (including microscopic gas bubbles) in the aqueous media were subsequently transported, continuously and directionally, on the copper wires with complex shapes. The small gas bubbles always moved to the larger ones.

Average properties of bidisperse bubbly flows

NASA Astrophysics Data System (ADS)

Serrano-García, J. C.; Mendez-Díaz, S.; Zenit, R.

2018-03-01

Experiments were performed in a vertical channel to study the properties of a bubbly flow composed of two distinct bubble size species. Bubbles were produced using a capillary bank with tubes with two distinct inner diameters; the flow through each capillary size was controlled such that the amount of large or small bubbles could be controlled. Using water and water-glycerin mixtures, a wide range of Reynolds and Weber number ranges were investigated. The gas volume fraction ranged between 0.5% and 6%. The measurements of the mean bubble velocity of each species and the liquid velocity variance were obtained and contrasted with the monodisperse flows with equivalent gas volume fractions. We found that the bidispersity can induce a reduction of the mean bubble velocity of the large species; for the small size species, the bubble velocity can be increased, decreased, or remain unaffected depending of the flow conditions. The liquid velocity variance of the bidisperse flows is, in general, bound by the values of the small and large monodisperse values; interestingly, in some cases, the liquid velocity fluctuations can be larger than either monodisperse case. A simple model for the liquid agitation for bidisperse flows is proposed, with good agreement with the experimental measurements.

Perturbation theory and numerical modelling of weakly and moderately nonlinear incompressible Richtmyer-Meshkov instability

NASA Astrophysics Data System (ADS)

Herrmann, M.; Velikovich, A. L.; Abarzhi, S. I.

2014-10-01

A study of incompressible two-dimensional Richtmyer-Meshkov instability by means of high-order Eulerian perturbation theory and numerical simulations is reported. Nonlinear corrections to Richtmyer's impulsive formula for the bubble and spike growth rates have been calculated analytically for arbitrary Atwood number and an explicit formula has been obtained for it in the Boussinesq limit. Conditions for early-time acceleration and deceleration of the bubble and the spike have been derived. In our simulations we have solved 2D unsteady Navier-Stokes equations for immiscible incompressible fluids using the finite volume fractional step flow solver NGA developed by, coupled to the level set based interface solver LIT,. The impact of small amounts of viscosity and surface tension on the RMI flow dynamics is studied numerically. Simulation results are compared to the theory to demonstrate successful code verification and highlight the influence of the theory's ideal inviscid flow assumption. Theoretical time histories of the interface curvature at the bubble and spike tip and the profiles of vertical and horizontal velocities have been favorably compared to simulation results, which converge to the theoretical predictions as the Reynolds and Weber numbers are increased. Work supported by the US DOE/NNSA.

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.

Measurement of interactions between solid particles, liquid droplets, and/or gas bubbles in a liquid using an integrated thin film drainage apparatus.

PubMed

Wang, Louxiang; Sharp, David; Masliyah, Jacob; Xu, Zhenghe

2013-03-19

A novel device was designed to measure drainage dynamics of thin liquid films confined between a solid particle, an immiscible liquid droplet, and/or gas bubble. Equipped with a bimorph force sensor, a computer-interfaced video capture, and a data acquisition system, the newly designed integrated thin film drainage apparatus (ITFDA) allows for the direct and simultaneous measurements of force barrier, true film drainage time, and bubble/droplet deformation under a well-controlled external force, receding and advancing contact angles, capillary force, and adhesion (detachment) force between an air bubble or oil droplet and a solid, a liquid, or an air bubble in an immiscible liquid. Using the diaphragm of a high-frequency speaker as the drive mechanism for the air bubble or oil droplet attached to a capillary tube, this newly designed device is capable of measuring forces over a wide range of hydrodynamic conditions, including bubble approach and retract velocities up to 50 mm/s and displacement range up to 1 mm. The results showed that the ITFDA was capable of measuring hydrodynamic resistance, film drainage time, and other important physical parameters between air bubbles and solid particles in aqueous solutions. As an example of illustrating the versatility, the ITFDA was also applied to other important systems such as interactions between air bubble and oil droplet, two air bubbles, and two oil droplets in an aqueous solution.

Effect of dissolved gases in water on acoustic cavitation and bubble growth rate in 0.83 MHz megasonic of interest to wafer cleaning.

PubMed

Kang, Bong-Kyun; Kim, Min-Su; Park, Jin-Goo

2014-07-01

Changes in the cavitation intensity of gases dissolved in water, including H2, N2, and Ar, have been established in studies of acoustic bubble growth rates under ultrasonic fields. Variations in the acoustic properties of dissolved gases in water affect the cavitation intensity at a high frequency (0.83 MHz) due to changes in the rectified diffusion and bubble coalescence rate. It has been proposed that acoustic bubble growth rates rapidly increase when water contains a gas, such as hydrogen faster single bubble growth due to rectified diffusion, and a higher rate of coalescence under Bjerknes forces. The change of acoustic bubble growth rate in rectified diffusion has an effect on the damping constant and diffusivity of gas at the acoustic bubble and liquid interface. It has been suggested that the coalescence reaction of bubbles under Bjerknes forces is a reaction determined by the compressibility and density of dissolved gas in water associated with sound velocity and density in acoustic bubbles. High acoustic bubble growth rates also contribute to enhanced cavitation effects in terms of dissolved gas in water. On the other hand, when Ar gas dissolves into water under ultrasound field, cavitation behavior was reduced remarkably due to its lower acoustic bubble growth rate. It is shown that change of cavitation intensity in various dissolved gases were verified through cleaning experiments in the single type of cleaning tool such as particle removal and pattern damage based on numerically calculated acoustic bubble growth rates. Copyright © 2014 Elsevier B.V. All rights reserved.

Numerical Simulations of Inclusion Behavior in Gas-Stirred Ladles

NASA Astrophysics Data System (ADS)

Lou, Wentao; Zhu, Miaoyong

2013-06-01

A computation fluid dynamics-population balance model (CFD-PBM) coupled model has been proposed to investigate the bubbly plume flow and inclusion behavior including growth, size distribution, and removal in gas-stirred ladles, and some new and important phenomena and mechanisms were presented. For the bubbly plume flow, a modified k- ɛ model with extra source terms to account for the bubble-induced turbulence was adopted to model the turbulence, and the bubble turbulent dispersion force was taken into account to predict gas volume fraction distribution in the turbulent gas-stirred system. For inclusion behavior, the phenomena of inclusions turbulent random motion, bubbles wake, and slag eye forming on the molten steel surface were considered. In addition, the multiple mechanisms both that promote inclusion growth due to inclusion-inclusion collision caused by turbulent random motion, shear rate in turbulent eddy, and difference inclusion Stokes velocities, and the mechanisms that promote inclusion removal due to bubble-inclusion turbulence random collision, bubble-inclusion turbulent shear collision, bubble-inclusion buoyancy collision, inclusion own floatation near slag-metal interface, bubble wake capture, and wall adhesion were investigated. The importance of different mechanisms and total inclusion removal ratio under different conditions, and the distribution of inclusion number densities in ladle, were discussed and clarified. The results show that at a low gas flow rate, the inclusion growth is mainly attributed to both turbulent shear collision and Stokes collision, which is notably affected by the Stokes collision efficiency, and the inclusion removal is mainly attributed to the bubble-inclusion buoyancy collision and inclusion own floatation near slag-metal interface. At a higher gas flow rate, the inclusions appear as turbulence random motion in bubbly plume zone, and both the inclusion-inclusion and inclusion-bubble turbulent random collisions become important for inclusion growth and removal. With the increase of the gas flow rate, the total removal ratio increases, but when the gas flow rate exceeds 200 NL/min in 150-ton ladle, the total removal ration almost does not change. For the larger size inclusions, the number density in bubbly plume zone is less than that in the sidewall recirculation zones, but for the small size inclusions, the distribution of number density shows the opposite trend.

Time-resolved particle image velocimetry measurements of the 3D single-mode Richtmyer-Meshkov instability

NASA Astrophysics Data System (ADS)

Xu, Qian

The Richtmyer-Meshkov Instability (RMI) (Commun. Pure Appl. Math 23, 297-319, 1960; Izv. Akad. Nauk. SSSR Maekh. Zhidk. Gaza. 4, 151-157, 1969) occurs due to an impulsive acceleration acting on a perturbed interface between two fluids of different densities. In the experiments presented in this thesis, single mode 3D RMI experiments are performed. An oscillating speaker generates a single mode sinusoidal initial perturbation at an interface of two gases, air and SF6. A Mach 1.19 shock wave accelerates the interface and generates the Richtmyer-Meshkov Instability. Both gases are seeded with propylene glycol particles which are illuminated by an Nd: YLF pulsed laser. Three high-speed video cameras record image sequences of the experiment. Particle Image Velocimetry (PIV) is applied to measure the velocity field. Measurements of the amplitude for both spike and bubble are obtained, from which the growth rate is measured. For both spike and bubble experiments, amplitude and growth rate match the linear stability theory at early time, but fall into a non-linear region with amplitude measurements lying between the modified 3D Sadot et al. model ( Phys. Rev. Lett. 80, 1654-1657, 1998) and the Zhang & Sohn model (Phys. Fluids 9. 1106-1124, 1997; Z. Angew. Math Phys 50. 1-46, 1990) at late time. Amplitude and growth rate curves are found to lie above the modified 3D Sadot et al. model and below Zhang & Sohn model for the spike experiments. Conversely, for the bubble experiments, both amplitude and growth rate curves lie above the Zhang & Sohn model, and below the modified 3D Sadot et al. model. Circulation is also calculated using the vorticity and velocity fields from the PIV measurements. The calculated circulation are approximately equal and found to grow with time, a result that differs from the modified Jacobs and Sheeley's circulation model (Phys. Fluids 8, 405-415, 1996).

Interfacial gauge methods for incompressible fluid dynamics

PubMed Central

Saye, Robert

2016-01-01

Designing numerical methods for incompressible fluid flow involving moving interfaces, for example, in the computational modeling of bubble dynamics, swimming organisms, or surface waves, presents challenges due to the coupling of interfacial forces with incompressibility constraints. A class of methods, denoted interfacial gauge methods, is introduced for computing solutions to the corresponding incompressible Navier-Stokes equations. These methods use a type of “gauge freedom” to reduce the numerical coupling between fluid velocity, pressure, and interface position, allowing high-order accurate numerical methods to be developed more easily. Making use of an implicit mesh discontinuous Galerkin framework, developed in tandem with this work, high-order results are demonstrated, including surface tension dynamics in which fluid velocity, pressure, and interface geometry are computed with fourth-order spatial accuracy in the maximum norm. Applications are demonstrated with two-phase fluid flow displaying fine-scaled capillary wave dynamics, rigid body fluid-structure interaction, and a fluid-jet free surface flow problem exhibiting vortex shedding induced by a type of Plateau-Rayleigh instability. The developed methods can be generalized to other types of interfacial flow and facilitate precise computation of complex fluid interface phenomena. PMID:27386567

A Hele-Shaw-Cahn-Hilliard Model for Incompressible Two-Phase Flows with Different Densities

NASA Astrophysics Data System (ADS)

Dedè, Luca; Garcke, Harald; Lam, Kei Fong

2017-07-01

Topology changes in multi-phase fluid flows are difficult to model within a traditional sharp interface theory. Diffuse interface models turn out to be an attractive alternative to model two-phase flows. Based on a Cahn-Hilliard-Navier-Stokes model introduced by Abels et al. (Math Models Methods Appl Sci 22(3):1150013, 2012), which uses a volume-averaged velocity, we derive a diffuse interface model in a Hele-Shaw geometry, which in the case of non-matched densities, simplifies an earlier model of Lee et al. (Phys Fluids 14(2):514-545, 2002). We recover the classical Hele-Shaw model as a sharp interface limit of the diffuse interface model. Furthermore, we show the existence of weak solutions and present several numerical computations including situations with rising bubbles and fingering instabilities.

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

High-Speed Synchrotron X-ray Imaging Studies of the Ultrasound Shockwave and Enhanced Flow during Metal Solidification Processes

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

Tan, Dongyue; Lee, Tung Lik; Khong, Jia Chuan

2015-03-31

The highly dynamic behavior of ultrasonic bubble implosion in liquid metal, the multiphase liquid metal flow containing bubbles and particles, and the interaction between ultrasonic waves and semisolid phases during solidification of metal were studied in situ using the complementary ultrafast and high-speed synchrotron X-ray imaging facilities housed, respectively, at the Advanced Photon Source, Argonne National Laboratory, US, and Diamond Light Source, UK. Real-time ultrafast X-ray imaging of 135,780 frames per second revealed that ultrasonic bubble implosion in a liquid Bi-8 wt pctZn alloy can occur in a single wave period (30 kHz), and the effective region affected by themore » shockwave at implosion was 3.5 times the original bubble diameter. Furthermore, ultrasound bubbles in liquid metal move faster than the primary particles, and the velocity of bubbles is 70 similar to 100 pct higher than that of the primary particles present in the same locations close to the sonotrode. Ultrasound waves can very effectively create a strong swirling flow in a semisolid melt in less than one second. The energetic flow can detach solid particles from the liquid-solid interface and redistribute them back into the bulk liquid very effectively.« less

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

Perfluorinated Surfactant Chain-Length Effects on Sonochemical Kinetics

NASA Astrophysics Data System (ADS)

Campbell, Tammy Y.; Vecitis, Chad D.; Mader, Brian T.; Hoffmann, Michael R.

2009-08-01

The sonochemical degradation kinetics of the aqueous perfluorochemicals (PFCs) perfluorobutanoate (PFBA), perfluorobutanesulfonate (PFBS), perfluorohexanoate (PFHA), and perfluorohexanesulfonate (PFHS) have been investigated. Surface tension measurements were used to evaluate chain-length effects on equilibrium air-water interface partitioning. The PFC air-water interface partitioning coefficients, KeqPF, and maximum surface concentrations, ΓmaxPF, were determined from the surface pressure equation of state for PFBA, PFBS, PFHA, and PFHS. Relative KeqPF values were dependent upon chain length KeqPFHS ≅ 2.1KeqPFHA ≅ 3.9KeqPFBS ≅ 5.0KeqPFBA, whereas relative ΓmaxPF values had minimal chain length dependence ΓmaxPFHS ≅ ΓmaxPFHA ≅ ΓmaxPFBS ≅ 2.2ΓmaxPFBA. The rates of sonolytic degradation were determined over a range of frequencies from 202 to 1060 kHz at dilute (<1 μM) initial PFC concentrations and are compared to previously reported results for their C8 analogs: perfluorooctanesulfonate (PFOS) and perfluorooctanoate (PFOA). Under all conditions, the time-dependent PFC sonolytic degradation was observed to follow pseudo-first-order kinetics, i.e., below kinetic saturation, suggesting bubble-water interface populations were significantly below the adsorption maximum. The PFHX (where X = A or S) sonolysis rate constant was observed to peak at an ultrasonic frequency of 358 kHz, similar to that for PFOX. In contrast, the PFBX degradation rate constants had an apparent maximum at 610 kHz. Degradation rates observed for PFHX are similar to previously determined PFOX rates, kapp,358PFOX ≅ kapp,358PFHX. PFOX is sonolytically pyrolyzed at the transiently cavitating bubble-water interface, suggesting that rates should be proportional to equilibrium interfacial partitioning. However, relative equilibrium air-water interfacial partitioning predicts that KeqPFOX ≅ 5KeqPFHX. This suggests that at dilute PFC concentrations, adsorption to the bubble-water interface is ultrasonically enhanced due to high-velocity radial bubble oscillations. PFC sonochemical kinetics are slower for PFBS and further diminished for PFBA as compared to longer analogs, suggesting that PFBX surface films are of lower stability due to their greater water solubility.

Liquid phase stabilization versus bubble formation at a nanoscale curved interface

NASA Astrophysics Data System (ADS)

Schiffbauer, Jarrod; Luo, Tengfei

2018-03-01

We investigate the nature of vapor bubble formation near a nanoscale-curved convex liquid-solid interface using two models: an equilibrium Gibbs model for homogenous nucleation, and a nonequilibrium dynamic van der Waals-diffuse-interface model for phase change in an initially cool liquid. Vapor bubble formation is shown to occur for sufficiently large radius of curvature and is suppressed for smaller radii. Solid-fluid interactions are accounted for and it is shown that liquid-vapor interfacial energy, and hence Laplace pressure, has limited influence over bubble formation. The dominant factor is the energetic cost of creating the solid-vapor interface from the existing solid-liquid interface, as demonstrated via both equilibrium and nonequilibrium arguments.

Bubble bursting at an interface

NASA Astrophysics Data System (ADS)

Kulkarni, Varun; Sajjad, Kumayl; Anand, Sushant; Fezzaa, Kamel

2017-11-01

Bubble bursting is crucial to understanding the life span of bubbles at an interface and more importantly the nature of interaction between the bulk liquid and the outside environment from the point of view of chemical and biological material transport. The dynamics of the bubble as it rises from inside the liquid bulk to its disappearance on the interface after bursting is an intriguing process, many aspects of which are still being explored. In our study, we make detailed high speed imaging measurements to examine carefully the hole initiation and growth in bursting bubbles that unearth some interesting features of the process. Previous analyses available in literature are revisited based on our novel experimental visualizations. Using a combination of experiments and theory we investigate the role of various forces during the rupturing process. This work aims to further our current knowledge of bubble dynamics at an interface with an aim of predicting better the bubble evolution from its growth to its eventual integration with the liquid bulk.

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.

Propagation of misfit dislocations from buffer/Si interface into Si

DOEpatents

Liliental-Weber, Zuzanna [El Sobrante, CA; Maltez, Rogerio Luis [Porto Alegre, BR; Morkoc, Hadis [Richmond, VA; Xie, Jinqiao [Raleigh, VA

2011-08-30

Misfit dislocations are redirected from the buffer/Si interface and propagated to the Si substrate due to the formation of bubbles in the substrate. The buffer layer growth process is generally a thermal process that also accomplishes annealing of the Si substrate so that bubbles of the implanted ion species are formed in the Si at an appropriate distance from the buffer/Si interface so that the bubbles will not migrate to the Si surface during annealing, but are close enough to the interface so that a strain field around the bubbles will be sensed by dislocations at the buffer/Si interface and dislocations are attracted by the strain field caused by the bubbles and move into the Si substrate instead of into the buffer epi-layer. Fabrication of improved integrated devices based on GaN and Si, such as continuous wave (CW) lasers and light emitting diodes, at reduced cost is thereby enabled.

Investigation of Nucleate Boiling Mechanisms Under Microgravity Conditions

NASA Technical Reports Server (NTRS)

Dhir, V. K.; Qiu, D. M.; Ramanujapu, N.; Hasan, M. M.

1999-01-01

The present work is aimed at the experimental studies and numerical modeling of the bubble growth mechanisms of a single bubble attached to a heating surface and of a bubble sliding along an inclined heated plate. Single artificial cavity of 10 microns in diameter was made on the polished Silicon wafer which was electrically heated at the back side in order to control the surface nucleation superheat. Experiments with a sliding bubble were conducted at different inclination angles of the downward facing heated surface for the purpose of studying the effect of magnitude of components of gravity acting parallel to and normal to the heat transfer surface. Information on the bubble shape and size, the bubble induced liquid velocities as well as the surface temperature were obtained using the high speed imaging and hydrogen bubble techniques. Analytical/numerical models were developed to describe the heat transfer through the micro-macro layer underneath and around a bubble formed at a nucleation site. In the micro layer model the capillary and disjoining pressures were included. Evolution of the bubble-liquid interface along with induced liquid motion was modeled. As a follow-up to the studies at normal gravity, experiments are being conducted in the KC-135 aircraft to understand the bubble growth/detachment under low gravity conditions. Experiments have been defined to be performed under long duration of microgravity conditions in the space shuttle. The experiment in the space shuttle will provide bubble growth and detachment data at microgravity and will lead to validation of the nucleate boiling heat transfer model developed from the preceding studies conducted at normal and low gravity (KC-135) conditions.

Forced convection flow boiling and two-phase flow phenomena in a microchannel

NASA Astrophysics Data System (ADS)

Na, Yun Whan

2008-07-01

The present study was performed to numerically analyze the evaporation phenomena through the liquid-vapor interface and to investigate bubble dynamics and heat transfer behavior during forced convective flow boiling in a microchannel. Flow instabilities of two-phase flow boiling in a microchannel were studied as well. The main objective of this research is to investigate the fundamental mechanisms of two-phase flow boiling in a microchannel and provide predictive tools to design thermal management systems, for example, microchannel heat sinks. The numerical results obtained from this study were qualitatively and quantitatively compared with experimental results in the open literature. Physical and mathematical models, accounting for evaporating phenomena through the liquid-vapor interface in a microchannel at constant heat flux and constant wall temperature, have been developed, respectively. The heat transfer mechanism is affected by the dominant heat conduction through the thin liquid film and vaporization at the liquid-vapor interface. The thickness of the liquid film and the pressure of the liquid and vapor phases were simultaneously solved by the governing differential equations. The developed semi-analytical evaporation model that takes into account of the interfacial phenomena and surface tension effects was used to obtain solutions numerically using the fourth-order Runge-Kutta method. The effects of heat flux 19 and wall temperature on the liquid film were evaluated. The obtained pressure drops in a microchannel were qualitatively consistent with the experimental results of Qu and Mudawar (2004). Forced convective flow boiling in a single microchannel with different channel heights was studied through a numerical simulation to investigate bubble dynamics, flow patterns, and heat transfer. The momentum and energy equations were solved using the finite volume method while the liquid-vapor interface of a bubble is captured using the VOF (Volume of Fluid) technique. The effects of different constant heat fluxes and different channel heights on the boiling mechanisms were investigated. The effects of liquid velocity on the bubble departure diameter were analyzed. The obtained results showed that the wall superheats at the position of nucleate boiling are relatively independent of the mass flow rates at the same channel height. The obtained results, however, showed that the heat flux at the onset of nucleate boiling strongly depends on the channel height. With a decrease of the channel height and an increase of the liquid velocity at the channel inlet, the departure diameter of a bubble was smaller. The periodic flow patterns, such as the bubbly flow, elongated slug flow, and churn flow were observed in the microchannel. Flow instabilities of two-phase flow boiling in a trapezoidal microchannel using a three-dimensional model were investigated. Fluctuation behaviors of flow boiling parameters such as wall temperature and inlet pressure caused by periodic flow patterns were studied at different heat fluxes and mass fluxes. The numerical results showed large amplitude and short period oscillations for wall temperature and inlet pressure fluctuations. Stable and unstable flow boiling regime with short period oscillations were investigated. Those flow boiling regimes were not listed in stable and unstable boiling regime map proposed by Wang et al. (2007).

Bubble dynamics in microchannels: inertial and capillary migration forces

NASA Astrophysics Data System (ADS)

Rivero-Rodriguez, Javier; Scheid, Benoit

2018-05-01

This work focuses on the dynamics of a train of unconfined bubbles flowing in microchan- nels. We investigate the transverse position of a train of bubbles, its velocity and the associated pressure drop when flowing in a microchannel depending on the internal forces due to viscosity, inertia and capillarity. Despite the small scales of the system, inertia, referred to as inertial migration force, play a crucial role in determining the transverse equilibrium position of the bubbles. Beside inertia and viscosity, other effects may also affect the transverse migration of bubbles such as the Marangoni surface stresses and the surface deformability. We look at the influence of surfactants in the limit of infinite Marangoni effect which yields rigid bubble interface. The resulting migration force may balance external body forces if present such as buoyancy, Dean or magnetic ones. This balance not only determines the transverse position of the bubbles but, consequently, the surrounding flow structure, which can be determinant for any mass/heat transfer process involved. Finally, we look at the influence of the bubble deformation on the equilibrium position and compare it to the inertial migration force at the centred position, explaining the stable or unstable character of this position accordingly. A systematic study of the influence of the parameters - such as the bubble size, uniform body force, Reynolds and capillary numbers - has been carried out using numerical simulations based on the Finite Element Method, solving the full steady Navier-Stokes equations and its asymptotic counterpart for the limits of small Reynolds and/or capillary numbers.

Effervescence in champagne and sparkling wines: From bubble bursting to droplet evaporation

NASA Astrophysics Data System (ADS)

Séon, T.; Liger-Belair, G.

2017-01-01

When a bubble reaches an air-liquid interface, it ruptures, projecting a multitude of tiny droplets in the air. Across the oceans, an estimated 1018 to 1020 bubbles burst every second, and form the so called sea spray, a major player in earth's climate system. At a smaller scale, in a glass of champagne about a million bubbles nucleate on the wall, rise towards the surface and burst, giving birth to a particular aerosol that holds a concentrate of wine aromas. Based on the model experiment of a single bubble bursting in simple liquids, we depict each step of this effervescence, from bubble bursting to drop evaporation. In particular, we propose simple scaling laws for the jet velocity and the top drop size. We unravel experimentally the intricate roles of bubble shape, capillary waves, gravity, and liquid properties in the jet dynamics and the drop detachment. We demonstrate how damping action of viscosity produces faster and smaller droplets and more generally how liquid properties enable to control the bubble bursting aerosol characteristics. In this context, the particular case of Champagne wine aerosol is studied in details and the key features of this aerosol are identified. We demonstrate that compared to a still wine, champagne fizz drastically enhances the transfer of liquid into the atmosphere. Conditions on bubble radius and wine viscosity that optimize aerosol evaporation are provided. These results pave the way towards the fine tuning of aerosol characteristics and flavor release during sparkling wine tasting, a major issue of the sparkling wine industry.

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.

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.

Interfacial gauge methods for incompressible fluid dynamics

DOE PAGES

Saye, R.

2016-06-10

Designing numerical methods for incompressible fluid flow involving moving interfaces, for example, in the computational modeling of bubble dynamics, swimming organisms, or surface waves, presents challenges due to the coupling of interfacial forces with incompressibility constraints. A class of methods, denoted interfacial gauge methods, is introduced for computing solutions to the corresponding incompressible Navier-Stokes equations. These methods use a type of "gauge freedom" to reduce the numerical coupling between fluid velocity, pressure, and interface position, allowing high-order accurate numerical methods to be developed more easily. Making use of an implicit mesh discontinuous Galerkin framework, developed in tandem with this work,more » high-order results are demonstrated, including surface tension dynamics in which fluid velocity, pressure, and interface geometry are computed with fourth-order spatial accuracy in the maximum norm. Applications are demonstrated with two-phase fluid flow displaying fine-scaled capillary wave dynamics, rigid body fluid-structure interaction, and a fluid-jet free surface flow problem exhibiting vortex shedding induced by a type of Plateau-Rayleigh instability. The developed methods can be generalized to other types of interfacial flow and facilitate precise computation of complex fluid interface phenomena.« less

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.

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.

Bubble Induced Disruption of a Planar Solid-Liquid Interface During Controlled Directional Solidification in a Microgravity Environment

NASA Technical Reports Server (NTRS)

Grugel, Richard N.; Brush, Lucien N.; Anilkumar, Amrutur V.

2013-01-01

Pore Formation and Mobility Investigation (PFMI) experiments were conducted in the microgravity environment aboard the International Space Station with the intent of better understanding the role entrained porosity/bubbles play during controlled directional solidification. The planar interface in a slowing growing succinonitrile - 0.24 wt% water alloy was being observed when a nitrogen bubble traversed the mushy zone and remained at the solid-liquid interface. Breakdown of the interface to shallow cells subsequently occurred, and was later evaluated using down-linked data from a nearby thermocouple. These results and other detrimental effects due to the presence of bubbles during solidification processing in a microgravity environment are presented and discussed.

Bubble clustering in a glass of stout beer

NASA Astrophysics Data System (ADS)

Iwatsubo, Fumiya; Watamura, Tomoaki; Sugiyama, Kazuyasu

2017-11-01

To clarify why the texture in stout beer poured into a pint glass descends, we investigated local time development of the void fraction and velocity of bubbles. The propagation of the number density distribution, i.e. the texture, appearing near the inclined wall is observed. We visualized individual advective bubbles near the inclined wall by microscope and measured the local void fraction using brightness of images while the velocity of bubbles by means of Particle Tracking Velocimetry. As the result of measurements, we found the local void fraction and the bubbles advection velocity increase and decrease repeatedly with a time delay. We conclude the texture pattern is composed of fluid blobs which contain less bubbles; extruding and suction flows respectively toward and from the interior of the container form respectively in front and back of the blobs.

Investigation of Mechanisms Associated with Nucleate Boiling Under Microgravity Conditions

NASA Technical Reports Server (NTRS)

Dhir, Vijay K.

1996-01-01

The focus of the present work is to experimentally study and to analytically/numerically model the mechanisms of growth of bubbles attached to, and sliding along, a heated surface. To control the location of the active cavities, the number, the spacing, and the nucleation superheat, artificial cavities will be formed on silicon wafers. In order to study the effect of magnitude of components of gravitational acceleration acting parallel to, and normal to the surface, experiments will be conducted on surfaces inclined at different angles including a downward facing surface. Information on the temperature field around bubbles, bubble shape and size, and bubble induced liquid velocities will be obtained through the use of holography, video/high speed photography and hydrogen bubble techniques, respectively. Analytical/numerical models will be developed to describe the heat transfer including that through the micro-macro layer underneath and around a bubble. In the micro layer model capillary and disjoining pressures will be included. Evolution of the interface along with induced liquid motion will be modelled. Subsequent to the world at normal gravity, experiments will be conducted in the KC-135 or the Lear jet especially to learn about bubble growth/detachment under low gravity conditions. Finally, an experiment will be defined to be conducted under long duration of microgravity conditions in the space shuttle. The experiment in the space shuttle will provide microgravity data on bubble growth and detachment and will lead to a validation of the nucleate boiling heat transfer model developed from the preceding studies performed at normal and low gravity (KC-135 or Lear jet) conditions.

Cavitation and bubble dynamics: the Kelvin impulse and its applications

PubMed Central

Blake, John R.; Leppinen, David M.; Wang, Qianxi

2015-01-01

Cavitation and bubble dynamics have a wide range of practical applications in a range of disciplines, including hydraulic, mechanical and naval engineering, oil exploration, clinical medicine and sonochemistry. However, this paper focuses on how a fundamental concept, the Kelvin impulse, can provide practical insights into engineering and industrial design problems. The pathway is provided through physical insight, idealized experiments and enhancing the accuracy and interpretation of the computation. In 1966, Benjamin and Ellis made a number of important statements relating to the use of the Kelvin impulse in cavitation and bubble dynamics, one of these being ‘One should always reason in terms of the Kelvin impulse, not in terms of the fluid momentum…’. We revisit part of this paper, developing the Kelvin impulse from first principles, using it, not only as a check on advanced computations (for which it was first used!), but also to provide greater physical insights into cavitation bubble dynamics near boundaries (rigid, potential free surface, two-fluid interface, flexible surface and axisymmetric stagnation point flow) and to provide predictions on different types of bubble collapse behaviour, later compared against experiments. The paper concludes with two recent studies involving (i) the direction of the jet formation in a cavitation bubble close to a rigid boundary in the presence of high-intensity ultrasound propagated parallel to the surface and (ii) the study of a ‘paradigm bubble model’ for the collapse of a translating spherical bubble, sometimes leading to a constant velocity high-speed jet, known as the Longuet-Higgins jet. PMID:26442141

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.

Molecular dynamics simulations of ejecta production from sinusoidal tin surfaces under supported and unsupported shocks

NASA Astrophysics Data System (ADS)

Wu, Bao; Wu, FengChao; Zhu, YinBo; Wang, Pei; He, AnMin; Wu, HengAn

2018-04-01

Micro-ejecta, an instability growth process, occurs at metal/vacuum or metal/gas interface when compressed shock wave releases from the free surface that contains surface defects. We present molecular dynamics (MD) simulations to investigate the ejecta production from tin surface shocked by supported and unsupported waves with pressures ranging from 8.5 to 60.8 GPa. It is found that the loading waveforms have little effect on spike velocity while remarkably affect the bubble velocity. The bubble velocity of unsupported shock loading remains nonzero constant value at late time as observed in experiments. Besides, the time evolution of ejected mass in the simulations is compared with the recently developed ejecta source model, indicating the suppressed ejection of unmelted or partial melted materials. Moreover, different reference positions are chosen to characterize the amount of ejecta under different loading waveforms. Compared with supported shock case, the ejected mass of unsupported shock case saturates at lower pressure. Through the analysis on unloading path, we find that the temperature of tin sample increases quickly from tensile stress state to zero pressure state, resulting in the melting of bulk tin under decaying shock. Thus, the unsupported wave loading exhibits a lower threshold pressure causing the solid-liquid phase transition on shock release than the supported shock loading.

Influences of Scavenging and Removal of Surfactants by Bubble Processing on Primary Marine Aerosol Production from North Atlantic Seawater

NASA Astrophysics Data System (ADS)

Duplessis, P.; Chang, R.; Frossard, A. A.; Keene, W. C.; Maben, J. R.; Long, M. S.; Beaupre, S. R.; Kieber, D. J.; Kinsey, J. D.; Zhu, Y.; Lu, X.; Bisgrove, J.

2017-12-01

Primary marine aerosol particles (PMA) are produced by bursting bubbles from breaking waves at the air-sea interface and significantly modulate atmospheric chemical transformations and cloud properties. Surfactants in bulk seawater rapidly (seconds) adsorb onto fresh bubble surfaces forming organic films that influence size, rise velocity, bursting behavior, and associated PMA emissions. During a cruise on the R/V Endeavor in September and October 2016, PMA production from biologically productive and oligotrophic seawater was investigated at four stations in the western North Atlantic Ocean. PMA were produced in a high-capacity generator via turbulent mixing of seawater and clean air in a Venturi nozzle. When the flow of fresh seawater through the generator was turned off, surfactant depletion via bubble processing resulted in greater PMA mass production efficiencies per unit air detrained but had no consistent influence on number production efficiencies. The greater (factor of 3) production efficiencies of organic matter associated with PMA generated with the Venturi relative to those generated with frits during previous campaigns contributed to a faster depletion of surfactants from the seawater reservoir and corresponding divergence in response.

Bursting Bubbles from Combustion of Thermoplastic Materials in Microgravity

NASA Technical Reports Server (NTRS)

Butler, K. B.

1999-01-01

Many thermoplastic materials in common use for a wide range of applications, including spacecraft, develop bubbles internally as they burn due to chemical reactions taking place within the bulk. These bubbles grow and migrate until they burst at the surface, forceably ejecting volatile gases and, occasionally, molten fuel. In experiments in normal gravity, Kashiwagi and Ohlemiller observed vapor jets extending a few centimeters from the surface of a radiatively heated polymethylmethacrylate (PMMA) sample, with some molten material ejected into the gas phase. These physical phenomena complicated the combustion process considerably. In addition to the non-steady release of volatiles, the depth of the surface layer affected by oxygen was increased, attributed to the roughening of the surface by bursting events. The ejection of burning droplets in random directions presents a potential fire hazard unique to microgravity. In microgravity combustion experiments on nylon Velcro fasteners and on polyethylene wire insulation, the presence of bursting fuel vapor bubbles was associated with the ejection of small particles of molten fuel as well as pulsations of the flame. For the nylon fasteners, particle velocities were higher than 30 cm/sec. The droplets burned robustly until all fuel was consumed, demonstrating the potential for the spread of fire in random directions over an extended distance. The sequence of events for a bursting bubble has been photographed by Newitt et al.. As the bubble reaches the fluid surface, the outer surface forms a dome while the internal bubble pressure maintains a depression at the inner interface. Liquid drains from the dome until it breaks into a cloud of droplets on the order of a few microns in size. The bubble gases are released rapidly, generating vortices in the quiescent surroundings and transporting the tiny droplets. The depression left by the escaping gases collapses into a central jet, which rises with a high velocity and may break up, releasing one or more relatively large drops (on the order of a millimeter in these experiments). A better understanding of bubble development and bursting processes, the effects of bursting behavior on burning rate of the bulk material, and the circumstances under which large droplets are expelled, as well as their trajectories, sizes, and burning rates, is sought through computer modeling compared with experiment.

Role of entrapped vapor bubbles during microdroplet evaporation

NASA Astrophysics Data System (ADS)

Putnam, Shawn A.; Byrd, Larry W.; Briones, Alejandro M.; Hanchak, Michael S.; Ervin, Jamie S.; Jones, John G.

2012-08-01

On superheated surfaces, the air bubble trapped during impingement grows into a larger vapor bubble and oscillates at the frequency predicted for thermally induced capillary waves. In some cases, the entrapped vapor bubble penetrates the droplet interface, leaving a micron-sized coffee-ring pattern of pure fluid. Vapor bubble entrapment, however, does not influence the evaporation rate. This is also true on laser heated surfaces, where a laser can thermally excite capillary waves and induce bubble oscillations over a broad range of frequencies, suggesting that exciting perturbations in a pinned droplets interface is not an effective avenue for enhancing evaporative heat transfer.

A Finite Element Method for Simulation of Compressible Cavitating Flows

NASA Astrophysics Data System (ADS)

Shams, Ehsan; Yang, Fan; Zhang, Yu; Sahni, Onkar; Shephard, Mark; Oberai, Assad

2016-11-01

This work focuses on a novel approach for finite element simulations of multi-phase flows which involve evolving interface with phase change. Modeling problems, such as cavitation, requires addressing multiple challenges, including compressibility of the vapor phase, interface physics caused by mass, momentum and energy fluxes. We have developed a mathematically consistent and robust computational approach to address these problems. We use stabilized finite element methods on unstructured meshes to solve for the compressible Navier-Stokes equations. Arbitrary Lagrangian-Eulerian formulation is used to handle the interface motions. Our method uses a mesh adaptation strategy to preserve the quality of the volumetric mesh, while the interface mesh moves along with the interface. The interface jump conditions are accurately represented using a discontinuous Galerkin method on the conservation laws. Condensation and evaporation rates at the interface are thermodynamically modeled to determine the interface velocity. We will present initial results on bubble cavitation the behavior of an attached cavitation zone in a separated boundary layer. We acknowledge the support from Army Research Office (ARO) under ARO Grant W911NF-14-1-0301.

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.

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.

Carbon dioxide induced bubble formation in a CH4-CO2-H2O ternary system: a molecular dynamics simulation study.

PubMed

Sujith, K S; Ramachandran, C N

2016-02-07

The extraction of methane from its hydrates using carbon dioxide involves the decomposition of the hydrate resulting in a CH4-CO2-H2O ternary solution. Using classical molecular dynamics simulations, we investigate the evolution of dissolved gas molecules in the ternary system at different concentrations of CO2. Various compositions considered in the present study resemble the solution formed during the decomposition of methane hydrates at the initial stages of the extraction process. We find that the presence of CO2 aids the formation of CH4 bubbles by causing its early nucleation. Elucidation of the composition of the bubble revealed that in ternary solutions with high concentration of CO2, mixed gas bubbles composed of CO2 and CH4 are formed. To understand the role of CO2 in the nucleation of CH4 bubbles, the structure of the bubble formed was analyzed, which revealed that there is an accumulation of CO2 at the interface of the bubble and the surrounding water. The aggregation of CO2 at the bubble-water interface occurs predominantly when the concentration of CO2 is high. Radial distribution function for the CH4-CO2 pair indicates that there is an increasingly favorable direct contact between dissolved CH4 and CO2 molecules in the bubble-water interface. It is also observed that the presence of CO2 at the interface results in the decrease in surface tension. Thus, CO2 leads to greater stability of the bubble-water interface thereby bringing down the critical size of the bubble nuclei. The results suggest that a rise in concentration of CO2 helps in the removal of dissolved CH4 thereby preventing the accumulation of methane in the liquid phase. Thus, the presence of CO2 is predicted to assist the decomposition of methane hydrates in the initial stages of the replacement process.

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.

Two-phase damping and interface surface area in tubes with vertical internal flow

NASA Astrophysics Data System (ADS)

Béguin, C.; Anscutter, F.; Ross, A.; Pettigrew, M. J.; Mureithi, N. W.

2009-01-01

Two-phase flow is common in the nuclear industry. It is a potential source of vibration in piping systems. In this paper, two-phase damping in the bubbly flow regime is related to the interface surface area and, therefore, to flow configuration. Experiments were performed with a vertical tube clamped at both ends. First, gas bubbles of controlled geometry were simulated with glass spheres let to settle in stagnant water. Second, air was injected in stagnant alcohol to generate a uniform and measurable bubble flow. In both cases, the two-phase damping ratio is correlated to the number of bubbles (or spheres). Two-phase damping is directly related to the interface surface area, based on a spherical bubble model. Further experiments were carried out on tubes with internal two-phase air-water flows. A strong dependence of two-phase damping on flow parameters in the bubbly flow regime is observed. A series of photographs attests to the fact that two-phase damping in bubbly flow increases for a larger number of bubbles, and for smaller bubbles. It is highest immediately prior to the transition from bubbly flow to slug or churn flow regimes. Beyond the transition, damping decreases. It is also shown that two-phase damping increases with the tube diameter.

Modeling of Vapor Bubble Growth Under Nucleate Boiling Conditions in Reduced Gravity

NASA Technical Reports Server (NTRS)

Buyevich, Yu A.; Webbon, Bruce W.

1995-01-01

A dynamic model is developed to describe the evolution of a vapor bubble growing at a nucleation site on a superheated surface under arbitrary gravity. The bubble is separated from the surface by a thin microlayer and grows due to the evaporation from the microlayer interface. The average thickness of the microlayer increases as the bubble expands along the surface if the evaporation rate is lower than some critical value. The corresponding threshold value of the surface temperature has to be associated with the burn-out crisis. Two main reasons make for bubble separation, which are the buoyancy force and a force caused by the vapor momentum that comes to the bubble with vapor molecules. The latter force is somewhat diminished if condensation takes place at the upper bubble surface in subcooled liquids. The action of the said forces is opposed by inertia of the additional mass of liquid as the bubble center rises above the surface and by inertia of liquid being expelled by the growing bubble in radial directions. An extra pressure force arises due to the liquid inflow into the microlayer with a finite velocity. The last force helps in holding the bubble close to the surface during an initial stage of bubble evolution. Two limiting regimes with distinctly different properties can be singled out, depending on which of the forces that favor bubble detachment dominates. Under conditions of moderately reduced gravity, the situation is much the same as in normal gravity, although the bubble detachment volume increases as gravity diminishes. In microgravity, the buoyancy force is negligible. Then the bubble is capable of staying near the surface for a long time, with intensive evaporation from the microlayer. It suggests a drastic change in the physical mechanism of heat removal as gravity falls below a certain sufficiently low level. Inferences of the model and conclusions pertaining to effects caused on heat transfer processes by changes in bubble hydrodynamics induced by gravity are discussed in connection with experimental evidence, both available in current and in as yet unpublished literature.

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 breakup phenomena in a venturi tube

NASA Astrophysics Data System (ADS)

Fujiwara, Akiko

2005-11-01

Microbubble has distinguished characteristics of large surface area to unit volume and small buoyancy, and it has advantages in many engineering fields. Recently microbubble generators with low energy and high performance are required to wide applications. In the present study, we propose one new effective technique to generate tiny bubbles with less than 200 μm diameter utilizing venturi tube under high void fraction condition. The objective of the present study is to elucidate the mechanism of bubble breakup phenomena in the venturi tube and to clarify the effects of parameters which are necessary to realize an optimum system experimentally. Experiment was conducted with void fraction of 4% and variation of liquid velocity from 9 to 26 m/s at the throat. Under low velocity condition, bubbles which were observed with a high speed camera parted gradually in a wide region. On the contrary under high velocity condition, bubbles expanded after passing through the throat and shrank rapidly. Since the speed of sound in gas-liquid system is extremely lower than that of single-phase flow, the bubble breakup phenomenon in the venturi tube is explained as the supersonic flow in a Laval nozzle. By rapid pressure recovery in diverging area, expanding bubbles collapse violently. The tiny bubbles are generated due to the surface instability of shrinking bubbles.

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.

Computational Fluid Dynamics-Population Balance Model Simulation of Effects of Cell Design and Operating Parameters on Gas-Liquid Two-Phase Flows and Bubble Distribution Characteristics in Aluminum Electrolysis Cells

NASA Astrophysics Data System (ADS)

Zhan, Shuiqing; Wang, Junfeng; Wang, Zhentao; Yang, Jianhong

2018-02-01

The effects of different cell design and operating parameters on the gas-liquid two-phase flows and bubble distribution characteristics under the anode bottom regions in aluminum electrolysis cells were analyzed using a three-dimensional computational fluid dynamics-population balance model. These parameters include inter-anode channel width, anode-cathode distance (ACD), anode width and length, current density, and electrolyte depth. The simulations results show that the inter-anode channel width has no significant effect on the gas volume fraction, electrolyte velocity, and bubble size. With increasing ACD, the above values decrease and more uniform bubbles can be obtained. Different effects of the anode width and length can be concluded in different cell regions. With increasing current density, the gas volume fraction and electrolyte velocity increase, but the bubble size keeps nearly the same. Increasing electrolyte depth decreased the gas volume fraction and bubble size in particular areas and the electrolyte velocity increased.

Three-dimensional features on oscillating microbubbles streaming flows

NASA Astrophysics Data System (ADS)

Rossi, Massimiliano; Marin, Alvaro G.; Wang, Cheng; Hilgenfeldt, Sascha; Kähler, Christian J.

2013-11-01

Ultrasound-driven oscillating micro-bubbles have been used as active actuators in microfluidic devices to perform manifold tasks such as mixing, sorting and manipulation of microparticles. A common configuration consists in side-bubbles, created by trapping air pockets in blind channels perpendicular to the main channel direction. This configuration results in bubbles with a semi-cylindrical shape that creates a streaming flow generally considered quasi two-dimensional. However, recent experiments performed with three-dimensional velocimetry methods have shown how microparticles can present significant three-dimensional trajectories, especially in regions close to the bubble interface. Several reasons will be discussed such as boundary effects of the bottom/top wall, deformation of the bubble interface leading to more complex vibrational modes, or bubble-particle interactions. In the present investigation, precise measurements of particle trajectories close to the bubble interface will be performed by means of 3D Astigmatic Particle Tracking Velocimetry. The results will allow us to characterize quantitatively the three-dimensional features of the streaming flow and to estimate its implications in practical applications as particle trapping, sorting or mixing.

Modeling of single film bubble and numerical study of the plateau structure in foam system

NASA Astrophysics Data System (ADS)

Sun, Zhong-guo; Ni, Ni; Sun, Yi-jie; Xi, Guang

2018-02-01

The single-film bubble has a special geometry with a certain amount of gas shrouded by a thin layer of liquid film under the surface tension force both on the inside and outside surfaces of the bubble. Based on the mesh-less moving particle semi-implicit (MPS) method, a single-film double-gas-liquid-interface surface tension (SDST) model is established for the single-film bubble, which characteristically has totally two gas-liquid interfaces on both sides of the film. Within this framework, the conventional surface free energy surface tension model is improved by using a higher order potential energy equation between particles, and the modification results in higher accuracy and better symmetry properties. The complex interface movement in the oscillation process of the single-film bubble is numerically captured, as well as typical flow phenomena and deformation characteristics of the liquid film. In addition, the basic behaviors of the coalescence and connection process between two and even three single-film bubbles are studied, and the cases with bubbles of different sizes are also included. Furthermore, the classic plateau structure in the foam system is reproduced and numerically proved to be in the steady state for multi-bubble connections.

CFD analysis of hydrodynamic studies of a bubbling fluidized bed

NASA Astrophysics Data System (ADS)

Rao, B. J. M.; Rao, K. V. N. S.; Ranga Janardhana, G.

2018-03-01

Fluidization velocity is one of the most important parameter to characterize the hydrodynamic studies of fluidized bed asit determines different flow regimes. Computational Fluid Dynamics simulations are carriedfor a cylindrical bubbling fluidized bed with a static bed height 1m with 0.150m diameter of gasification chamber. The parameter investigated is fluidization velocity in range of 0.05m/s to 0.7m/s. Sand with density 2600kg/m3 and with a constant particle diameter of sand 385μm is employed for all the simulations. Simulations are conducted using the commercial Computational Fluid Dynamics software, ANSYS-FLUENT.The bubbling flow regime is appeared above the air inlet velocity of 0.2m/s. Bubbling character is increased with increase in inlet air velocities indicated by asymmetrical fluctuations of volume fractions in radial directions at different bed heights

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.

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.

The calculation of weakly non-spherical cavitation bubble impact on a solid

NASA Astrophysics Data System (ADS)

Aganin, A. A.; Guseva, T. S.; Kosolapova, L. A.; Khismatullina, N. A.

2016-11-01

The effect of small spheroidal non-sphericity of a cavitation bubble touching a solid at the beginning of its collapse on its impact on the solid of a copper-nickel alloy is investigated. The impact on the solid is realized by means of a high-speed liquid jet arising at collapse on the bubble surface. The shape of the jet, its velocity and pressure are calculated by the boundary element method. The spatial and temporal characteristics of the pressure pulses on the solid surface are determined by the CIP-CUP method on dynamically adaptive grids without explicitly separating the gas-liquid interface. The solid surface layer dynamics is evaluated by the Godunov method. The results are analyzed in dimensionless variables obtained with using the water hammer pressure, the time moment and the jet-solid contact area radius at which the jet begins to spread on the solid surface. It is shown that in those dimensionless variables, the dependence of the spatial and temporal characteristics of the solid surface pressure pulses on the initial bubble shape non-sphericity is relatively small. The nonsphericity also slightly influences the main qualitative features of the dynamic processes inside the solid, whereas its effect on their quantitative characteristics can be significant.

Advanced subgrid-scale modeling for convection-dominated species transport at fluid interfaces with application to mass transfer from rising bubbles

NASA Astrophysics Data System (ADS)

Weiner, Andre; Bothe, Dieter

2017-10-01

This paper presents a novel subgrid scale (SGS) model for simulating convection-dominated species transport at deformable fluid interfaces. One possible application is the Direct Numerical Simulation (DNS) of mass transfer from rising bubbles. The transport of a dissolving gas along the bubble-liquid interface is determined by two transport phenomena: convection in streamwise direction and diffusion in interface normal direction. The convective transport for technical bubble sizes is several orders of magnitude higher, leading to a thin concentration boundary layer around the bubble. A true DNS, fully resolving hydrodynamic and mass transfer length scales results in infeasible computational costs. Our approach is therefore a DNS of the flow field combined with a SGS model to compute the mass transfer between bubble and liquid. An appropriate model-function is used to compute the numerical fluxes on all cell faces of an interface cell. This allows to predict the mass transfer correctly even if the concentration boundary layer is fully contained in a single cell layer around the interface. We show that the SGS-model reduces the resolution requirements at the interface by a factor of ten and more. The integral flux correction is also applicable to other thin boundary layer problems. Two flow regimes are investigated to validate the model. A semi-analytical solution for creeping flow is used to assess local and global mass transfer quantities. For higher Reynolds numbers ranging from Re = 100 to Re = 460 and Péclet numbers between Pe =104 and Pe = 4 ṡ106 we compare the global Sherwood number against correlations from literature. In terms of accuracy, the predicted mass transfer never deviates more than 4% from the reference values.

Numerical analysis of bubble-cluster formation in an ultrasonic field

NASA Astrophysics Data System (ADS)

Kim, Donghyun; Son, Gihun

2016-11-01

Bubble-cluster formation in an ultrasonic field is investigated numerically solving the conservation equations of mass, momentum and energy. The liquid-gas interface is calculated using the volume-of-fluid method with variable gas density to consider the bubble compressibility. The effect of liquid-gas phase change is also included as the interface source terms of the mass and energy equations. The numerical approach is tested through the simulation of the expansion and contraction motion of a compressed bubble adjacent to a wall. When the bubble is placed in an ultrasonic field, it oscillates radially and then collapses violently. Numerical simulation is also performed for bubble-cluster formation induced by an ultrasonic generator, where the generated bubbles are merged into a macrostructure along the acoustic flow field. The effects of ultrasonic power and frequency, liquid properties and pool temperature on the bubble-cluster formation are investigated. This work was supported by the Korea Institute of Energy Research.

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)

Non-equilibrium phase stabilization versus bubble nucleation at a nanoscale-curved Interface

NASA Astrophysics Data System (ADS)

Schiffbauer, Jarrod; Luo, Tengfei

Using continuum dynamic van der Waals theory in a radial 1D geometry with a Lennard-Jones fluid model, we investigate the nature of vapor bubble nucleation near a heated, nanoscale-curved convex interface. Vapor bubble nucleation and growth are observed for interfaces with sufficiently large radius of curvature while phase stabilization of a superheated fluid layer occurs at interfaces with smaller radius. The hypothesis that the high Laplace pressure required for stable equilibrium of very small bubbles is responsible for phase stability is tested by effectively varying the parameter which controls liquid-vapor surface tension. In doing so, the liquid-vapor surface tension- hence Laplace pressure-is shown to have limited effect on phase stabilization vs. bubble nucleation. However, the strong dependence of nucleation on leading-order momentum transport, i.e. viscous dissipation, near the heated inner surface is demonstrated. We gratefully acknowledge ND Energy for support through the ND Energy Postdoctoral Fellowship program and the Army Research Office, Grant No. W911NF-16-1-0267, managed by Dr. Chakrapani Venanasi.

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

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.

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.

Ionospheric Research with Miniaturized Plasma Sensors Aboard FalconSAT-3

NASA Astrophysics Data System (ADS)

Habash Krause, L.; Herrero, F. A.; Chun, F. K.; McHarg, M. G.

2003-12-01

Investigations into a novel technique to measure ionosphere-thermosphere parameters have culminated in the Flat Plasma Spectrometer (FLAPS) experiment, presently under development through a collaboration between NASA Goddard Space Flight Center (GSFC) and the U. S. Air Force Academy (USAFA). FLAPS is capable of providing measurements of the full neutral wind vector, full ion-drift velocity vector, neutral and ion temperatures, and deviations from thermalization. In addition, coarse mass spectroscopy is possible using an energy analysis technique. The suite of instruments is comprised of a set of 16 individual neutral and ion analyzers, each of which is designed to perform a specific function. Advances in miniaturization technology have enabled a design in which the 16-sensor suite resides on a circular microchannel plate with an effective area of 25 cm2. The FLAPS electronics package, consisting of low voltage and high voltage power supplies, a microprocessor, and Application Specific Integrated Circuit (ASIC) amplifiers, requires a volume of 290 cm3, power of 1.5 W, and a mass of 500 g. The suite requires a +5V regulated power line from the spacecraft, and the telemetry interface is a 5.0 V TTL-compatible serial connection. Data collection rates vary from 1 to 1000 (192 Byte) spectra per second. The motivation for the FLAPS experiment is driven by objectives that fall into both basic science and technology demonstration categories. Scientifically, there is strong interest in the effects of ionosphere-thermosphere coupling and non-thermalized plasma on the processes associated with equatorial F-region ionospheric plasma bubbles. These bubbles have been known to scintillate transionospheric propagation of radio waves, often resulting in disruptions of space-based communication and navigation systems. FLAPS investigations will assist in quantifying the impact of various processes on the instigation or suppression of plasma bubbles; certain outstanding questions include 1) What is the relevance of meridional winds in suppression of plasma bubble growth? 2) What role does a velocity space instability driven by non-thermalized plasma play in the generation of small scale (<1 km) bubbles? 3) What process is responsible for turbulence in plasma beyond the edges of a bubble structure? Technologically, the need for small yet capable instruments arises from the desire to make multipoint in situ measurements of "microscopic" plasma parameters to provide insight into "macroscopic" phenomena. Examples include coherency of spatial boundaries of large-scale ( ˜100 km) plasma bubbles, three dimensional structure of the equatorial wind and temperature anomaly, and vertical velocity gradients in the low latitude ionosphere. This paper provides an overview of the experiment motivation and instrument design of the FLAPS experiment.

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

Recognition and measurement gas-liquid two-phase flow in a vertical concentric annulus at high pressures

NASA Astrophysics Data System (ADS)

Li, Hao; Sun, Baojiang; Guo, Yanli; Gao, Yonghai; Zhao, Xinxin

2018-02-01

The air-water flow characteristics under pressure in the range of 1-6 MPa in a vertical annulus were evaluated in this report. Time-resolved bubble rising velocity and void fraction were also measured using an electrical void fraction meter. The results showed that the pressure has remarkable effect on the density, bubble size and rise velocity of the gas. Four flow patterns (bubble, cap-bubble, cap-slug, and churn) were also observed instead of Taylor bubble at high pressure. Additionally, the transition process from bubble to cap-bubble was investigated at atmospheric and high pressures, respectively. The results revealed that the flow regime transition criteria for atmospheric pressure do not work at high pressure, hence a new flow regime transition model for annular flow channel geometry was developed to predict the flow regime transition, which thereafter exhibited high accuracy at high pressure condition.

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.

Topics in viscous potential flow of two-phase systems

NASA Astrophysics Data System (ADS)

Padrino Inciarte, Juan Carlos

Two-phase flows are ubiquitous, from natural and domestic environments to industrial settings. However, due to their complexity, modeling these fluid systems remains a challenge from both the perspective of fundamental questions on the dynamics of an individual, smooth interface, and the perspective of integral analyses, which involve averaging of the conservation laws over large domains, thereby missing local details of the flow. In this work, we consider a set of five problems concerning the linear and non-linear dynamics of an interface or free surface and the study of cavitation inception. Analyses are carried out by assuming the fluid motion to be irrotational, that is, with zero vorticity, and the fluids to be viscous, although results from rotational analyses are presented for the purpose of comparison. The problems considered here are the following: First, we analyze the non-linear deformation and break-up of a bubble or drop immersed in a uniaxial extensional flow of an incompressible viscous fluid. The method of viscous potential flow, in which the flow field is irrotational and viscosity enters through the balance of normal stresses at the interface, is used in the analysis. The governing equations are solved numerically to track the motion of the interface by coupling a boundary element method with a time-integration routine. When break-up occurs, the break-up time computed here is compared with results obtained elsewhere from numerical simulations of the Navier.Stokes equations, which thus keeps vorticity in the analysis, for several combinations of the relevant dimensionless parameters of the problem. For the bubble, for Weber numbers 3 ≤ We ≤ 6, predictions from viscous potential flow shows good agreement with the results from the Navier.Stokes equations for the bubble break-up time, whereas for larger We, the former underpredicts the results given by the latter. Including viscosity increases the break-up time with respect to the inviscid case. For the drop, as expected, increasing the viscous effects of the irrotational motion produces large, elongated drops that take longer to break up in comparison with results for inviscid fluids. In the second problem, we compute the force acting on a spherical bubble of variable radius moving within a liquid with an outer spherical boundary. Viscous potential flow and the dissipation method, which is another purely irrotational approach stemming from the mechanical energy equation, are both systematically implemented. This exposes the role of the choice of the outer boundary condition for the stress on the drag, an issue not explained in the literature known to us. By means of the well-known "cell-model" analysis, the results for the drag are then applied to the case of a swarm of rising bubbles having a certain void fraction. Computations from the dissipation method for the drag coefficient and rise velocity for a bubble swarm agree with numerical solutions; evaluation against experimental data for high Reynolds and low Weber numbers shows that all the models considered, including those given in the literature, overpredict the bubble swarm rise velocity. In the next two problems, we apply the analysis of viscous potential flow and the dissipation method to study the linear dynamics of waves of "small" amplitude acting either on a plane or on a spherical interface separating a liquid from a dynamically inactive fluid. It is shown that the viscous irrational theories exhibit the features of the wave dynamics by comparing with the exact solution. The range of parameters for which good agreement with the exact solution exists is presented. The general trend shows that for long waves the dissipation method results in the best approximation, whereas for short waves, even for very viscous liquids, viscous potential flow demonstrates better agreement. Finally, the problem of cavitation inception for the flow of a viscous liquid past a stationary sphere is studied by means of the theory of stress-induced cavitation. The flow field for a single phase needed in the analysis is found from three different methods, namely, the numerical solution of the Navier--Stokes equations, the irrotational motion of a viscous fluid, and, in the limit of no inertia, the Stokes flow formulation. The new predictions are then compared with those obtained from the classical pressure criterion. The main finding is that at a fixed cavitation number more viscous liquids are at greater risk to cavitation.

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

Barber, Jacqueline; Aix-Marseille Universite; Brutin, David

Boiling in microchannels remains elusive due to the lack of full understanding of the mechanisms involved. A powerful tool in achieving better comprehension of the mechanisms is detailed imaging and analysis of the two-phase flow at a fundamental level. Boiling is induced in a single microchannel geometry (hydraulic diameter 727 {mu}m), using a refrigerant FC-72, to investigate the effect of channel confinement on bubble growth. A transparent, metallic, conductive deposit has been developed on the exterior of the rectangular microchannel, allowing simultaneous uniform heating and visualisation to be achieved. The data presented in this paper is for a particular casemore » with a uniform heat flux applied to the microchannel and inlet liquid mass flowrate held constant. In conjunction with obtaining high-speed images and videos, sensitive pressure sensors are used to record the pressure drop across the microchannel over time. Bubble nucleation and growth, as well as periodic slug flow, are observed in the microchannel test section. The periodic pressure fluctuations evidenced across the microchannel are caused by the bubble dynamics and instances of vapour blockage during confined bubble growth in the channel. The variation of the aspect ratio and the interface velocities of the growing vapour slug over time, are all observed and analysed. We follow visually the nucleation and subsequent both 'free' and 'confined' growth of a vapour bubble during flow boiling of FC-72 in a microchannel, from analysis of our results, images and video sequences with the corresponding pressure data obtained. (author)« less

On the origin of the driving force in the Marangoni propelled gas bubble trapping mechanism.

PubMed

Miniewicz, A; Quintard, C; Orlikowska, H; Bartkiewicz, S

2017-07-19

Gas bubbles can be trapped and then manipulated with laser light. In this report, we propose the detailed optical trapping mechanism of gas bubbles confined inside a thin light-absorbing liquid layer between two glass plates. The necessary condition of bubble trapping in this case is the direct absorption of light by the solution containing a dye. Due to heat release, fluid whirls propelled by the surface Marangoni effect at the liquid/gas interface emerge and extend to large distances. We report the experimental microscopic observation of the origin of whirls at an initially flat liquid/air interface as well as at the curved interface of a liquid/gas bubble and support this finding with advanced numerical simulations using the finite element method within the COMSOL Multiphysics platform. The simulation results were in good agreement with the observations, which allowed us to propose a simple physical model for this particular trapping mechanism, to establish the origin of forces attracting bubbles toward a laser beam and to predict other phenomena related to this effect.

Sharp Interface Tracking in Rotating Microflows of Solvent Extraction

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

Glimm, James; Almeida, Valmor de; Jiao, Xiangmin

2013-01-08

The objective of this project is to develop a specialized sharp interface tracking simulation capability for predicting interaction of micron-sized drops and bubbles in rotating flows relevant to optimized design of contactor devices used in solvent extraction processes of spent nuclear fuel reprocessing. The primary outcomes of this project include the capability to resolve drops and bubbles micro-hydrodynamics in solvent extraction contactors, determining from first principles continuum fluid mechanics how micro-drops and bubbles interact with each other and the surrounding shearing fluid for realistic flows. In the near term, this effort will play a central role in providing parameters andmore » insight into the flow dynamics of models that average over coarser scales, say at the millimeter unit length. In the longer term, it will prove to be the platform to conduct full-device, detailed simulations as parallel computing power reaches the exaflop level. The team will develop an accurate simulation tool for flows containing interacting droplets and bubbles with sharp interfaces under conditions that mimic those found in realistic contactor operations. The main objective is to create an off-line simulation capability to model drop and bubble interactions in a domain representative of the averaged length scale. The technical approach is to combine robust interface tracking software, subgrid modeling, validation quality experiments, powerful computational hardware, and a team with simulation modeling, physical modeling and technology integration experience. Simulations will then fully resolve the microflow of drops and bubbles at the microsecond time scale. This approach is computationally intensive but very accurate in treating important coupled physical phenomena in the vicinity of interfaces. The method makes it possible to resolve spatial scales smaller than the typical distance between bubbles and to model some non-equilibrium thermodynamic features such as finite critical tension in cavitating liquids« less

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.

Effects of Force Fields on Interface Dynamics, in view of Two-Phase Heat Transfer Enhancement and Phase Management for Space Applications

NASA Astrophysics Data System (ADS)

Di Marco, P.; Saccone, G.

2017-11-01

On earth, gravity barely influences the dynamics of interfaces. For what concerns bubbles, buoyancy governs the dynamics of boiling mechanism and thus affects boiling heat transfer capacity. While, for droplets, the coupled effects of wettability and gravity affects interface exchanges. In space, in the lack of gravity, rules are changed and new phenomena come into play. The present work is aimed to study the effects of electric field on the shape and behaviour of bubbles and droplets in order to understand how to handle microgravity applications; in particular, the replacement of gravity with electric field and their coupled effects are evaluated. The experiments spread over different setups, gravity conditions, working fluids, interface conditions. Droplets and bubbles have been analysed with and without electric field, with and without (adiabatic) heat and mass transfer across the interface. Furthermore, the results of the 4 ESA Parabolic Flight Campaigns (PFC 58, 60, 64 & 66), for adiabatic bubbles, adiabatic droplets and evaporating droplets, will be summarized, discussed, and compared with the ground tests.

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.

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.