Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Multi-scale Mathematical Model of Growth and Coalescence of Bubbles Beneath the Anode in an Aluminum Reduction Cell

Modeling of bubble shapes is challenging because of the wide range of length scales in aluminum reduction cells. A 3D multi-scale mathematical model was developed to understand the nucleation, growth and coalescence of bubbles beneath the anode and to investigate the transition of bubbles from the micro- to macro level. The motion of micro-bubbles is examined using the discrete bubble model (DBM) within a Lagrangian reference. An algorithm for the transition from discrete micro-bubbles to large bubbles, which are fully resolved by the volume of fluid (VOF) approach, is achieved using a user-defined function. The two-way coupling between discrete bubbles and continuous fluids is achieved by inter-phase momentum exchange. The model involves three kinds of bubble coalescence: micro-bubble coalescence is taken into account by the DBM; large bubbles swallowing up micro-bubbles are solved by the discrete-continuum transition model; the coalescence between large bubbles is handled by the VOF method. Numerical results show that the coverage and thickness of bubbles agree well with the experimental data in the literature. Bubble thickness stops increasing when the bubble elongates with the anode bottom, which is 4.0 to 4.5 mm. Meanwhile, two asymmetrical thick heads can be observed under the anode. Bubble release frequency increases with increasing current density.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Single-Bubble Flow Boiling Phenomena - Interface Tracking Simulation

Boiling, as one of the most efficient heat transfer mechanisms, is widely used in various engineering systems. Better understanding and modeling of this process remains a major challenge in multiphase flow research. The distribution of vapor in boiling system affects the heat transfer rate and may cause unfavorable conditions, such as heater burn-out. A number of flow regimes have been observed experimentally. In general, there are three regimes under boiling conditions: nucleate boiling, transition boiling, and film boiling. The nucleate boiling regime is categorized into partial boiling regime and fully developed nucleate boiling, according to behavior of bubble dynamics and heat transfer characteristics. In the first region, as the wall temperature increases, the fraction of the wall surface subject to nucleate boiling increases until bubble formation occupies the entire heated surface. Heat transfer in this region is a complex mixture of single-phase forced convection and nucleate boiling. As bubble density increases, bubbles turn to coalesce and form insulating vapor patches, which impedes liquid flow back to the surface. The heat transfer rate under this condition reduces dramatically. This point is called Critical Heat Flux (CHF), which leads to boiling crisis. The behavior of boiling crisis is dependent on local fluid conditions due to bubble growth and departure. Experiments are the ultimate tool for boiling model development and verification, but they still lack some capabilities needed for better understanding of the local mechanism of this phenomena. It's difficult to capture many important details at the pressure/temperature conditions of interest using experimental techniques, such as full three dimensional measurements of velocity and temperature distribution and estimates of heat flux partitioning during the boiling process. The interface tracking simulation (ITS) is one of the promising approaches to describe heat transfer of boiling phenomena and their underlying mechanisms. The advances of high performance computing (HPC) in recent years made it possible to apply DNS to a wide variety of adiabatic flows. Some work has been done to implement evaporation / condensation models as well. Lee and Nydahl performed a numerical simulation of bubble growth and departure assuming that the bubble has the shape of a truncated sphere and a micro-layer exists during bubble growth. A 2-D cylindrical bubbles simulation has been done by Son with a simplified model to predict the heat flux in a thin liquid micro-layer. A direct simulation of 2D film boiling is conducted by Juric and Tryggvason. The film boiling, below critical pressure with a modified Level Set Method, is studied by Son and Dhir. Ose and Kunugi have developed a boiling and condensation model on subcooled boiling phenomena. Fuchs performed numerical simulations of bubble growth and departure considering wall conduction effects. Fully transient heat and fluid flow is modeled with a free surface of a rising bubble and a periodic calculation of the whole cycling of a growing, detaching, and rising bubble. Mukherjee and Kandikar have done a numerical simulation of an evaporating meniscus on a moving heated surface with level set method. The heat transfer associated with the advancing and receding contact angle has been studied. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate single-bubble growth scenarios. For a single bubble simulation in infinite liquid, the analytical solution exists and can be compared with the numerical results for verification purposes. The presented capabilities will be utilized in the large scale DNS of boiling flows in the future. The finite-element based code, PHASTA is used for the simulation. It is a parallel, hierarchic, higher-order accurate, adaptive, stabilized, transient analysis flow solver, which has been shown to be an effective tool for a wide range of single and two-phase flow problems. It has been shown to be highly scalable on top supercomputers ({approx}85% scaling on up to 3 million mesh partitions using a 92 billion element mesh . PHASTA has been shown to reliably predict the details of adiabatic bubbly flows as well as single-phase flows with heat transfer. To develop the capability of physics-based boiling process, in the presented work the energy and momentum equations are coupled to simulate a single-bubble growth scenario. For a single bubble growth in superheated water, the analytical solution exists and can be compared with the numerical results for verification purposes. The second-order accurate spatial and time discretization is used for the flow solver in the presented work. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2000

Marangoni Effects on Near-Bubble Microscale Transport During Boiling of Binary Fluid Mixtures

In earlier investigations, Marangoni effects were observed to be the dominant mechanism of boiling transport in 2-propanol/water mixtures under reduced gravity conditions. In this investigation we have examined the mechanisms of binary mixture boiling by exploring the transport near a single bubble generated in a binary mixture between a heated surface and cold surface. The temperature field created in the liquid around the bubble produces vaporization over the portion of its interface near the heated surface and condensation over portions of its interface near the cold surface. Experiments were conducted using different mixtures of water and 2-propanol under 1g conditions and under reduced gravity conditions aboard the KC135 aircraft. Since 2-propanol is more volatile than water, there is a lower concentration of 2-propanol near the hot surface and a higher concentration of 2-propanol near the cold plate relative to the bulk quantity. This difference in interface concentration gives rise to strong Marangoni effects that move liquid toward the hot plate in the near bubble region for 2-propanol and water mixtures. In the experiments in this study, the pressure of the test system was maintained at about 5 kPa to achieve the full spectrum of boiling behavior (nucleate boiling, critical heat flux and film boiling) at low temperature and heat flux levels. Heat transfer data and visual documentation of the bubble shape were extracted from the experimental results. In the 1-g experiments at moderate to high heat flux levels, the bubble was observed to grow into a mushroom shape with a larger top portion near the cold plate due to the buoyancy effect. The shape of the bubble was somewhat affected by the cold plate subcooling and the superheat of the heated surface. At low superheat levels for the heated surface, several active nucleation sites were observed, and the vapor stems from them merged to form a larger bubble. The generation rate of vapor is moderate in this regime and the bubble shape is cylindrical in appearance. In some instances, the bubble interface appeared to oscillate. At higher applied heat flux levels, the top of the bubble became larger, apparently to provide more condensing interface area adjacent to the cold plate. Increasing the applied heat flux ultimately led to dry-out of the heated surface, with conditions just prior to dryout corresponding to the maximum heat flux (CHF). A more stable bubble was observed when the system attained the minimum heat flux (for film boiling). In this regime, most of the surface under the bottom of the bubble was dry with nucleate boiling sometimes occuring around the contact perimeter of the bubble at heated surface. Different variations (e.g. gap between two plates, molar concentration of the liquid mixture) of the experiments were examined to determine parametric effects on the boiling process and to determine the best conditions for the KC135 reduced gravity tests. Variation of the gap was found to have a minor impact on the CHF. However, reducing the gap between the hot and cold surface was observed to significantly reduce the minimum heat flux for fixed molar concentration of 2-propanol. In the reduced gravity experiments aboard the KC135 aircraft, the bubble formed in the 6.4 mm gap was generally cylindrical or barrel shaped and it increased its extent laterally as the surface superheat increased. In reduced gravity experiments, dryout of the heated surface under the bubble was observed to occur at a lower superheated temperature than for 1g conditions. Observed features of the boiling process and heat transfer data under reduced gravity will be discussed in detail. The results of the reduced gravity experiments will also be compared to those obtained in comparable 1g experiments. In tandem with the experiments we are also developing a computational model of the transport in the liquid surrounding the bubble during the boiling process. The computational model uses a level set method to model motion of the interface. It will incorporate a macroscale treatment of the transport in the liquid gap between the surfaces and a microscale treatment of transport in the regions between the bubble interface and the solid surfaces. The features of the model will be described in detail. Future research directions suggested by the results to date will also be discussed.

V. Carey↗
NASA NTRS2000

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

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.

McQuillen, J.↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Wall Area of Influence and Growing Wall Heat Transfer due to Sliding Bubbles in Subcooled Boiling Flow

A variety of dynamical features of sliding bubbles and their impact on wall heat transfer were observed at subcooled flow boiling conditions in a vertical square test channel. Among the wide range of parameters observed, we particularly focus in this paper on (i) the sliding bubbles’ effect on wall heat transfer (supplemantry discussion to the authors’ previous work in Yoo et al. (2016a,b)) and (ii) the wall area influenced by sliding bubbles in subcooled boiling flow. At first, this study reveals that the degree of wall heat transfer improvement due to sliding bubbles depended less on the wall superheat condition as the mass flux increased. Also, the sliding bubble trajectory was found to be one of the critical factors in order to properly describe the wall heat transfer associated with sliding bubbles. In particular, the wall area influenced by sliding bubbles depended strongly on both sliding bubble trajectory and sliding bubble size; the sliding bubble trajectory was also observed to be closely related to the sliding bubble size. Importantly, these results indicate the limitation of current approach in CFD analyses especially for the wall area of bubble influence. In addition, the analyses on the temporal fraction of bubbles’ residence (FR) along the heated wall show that the sliding bubbles typically travel through narrow path with high frequency while the opposite was observed downstream. That is, both FR and sliding bubble trajectory depended substantially on the distance from nucleation site, which is expected to be similar for the quenching heat transfer mode induced by sliding bubbles.

42 ENGINEERING↗
DOE OSTI2021

Hydrate Formation on Marine Seep Bubbles and the Implications for Water Column Methane Dissolution

Abstract Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate‐coated bubbles to deliver methane to shallower depths in the water column than hydrate‐free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble‐release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate‐free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas‐water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water‐column observations that most bubbles become hydrate‐coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1999

Bubble Generation in a Continuous Liquid Flow Under Reduced Gravity Conditions

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.

Pais, Salvatore Cezar↗
NASA NTRS2002

The Speed of Axial Propagation of a Cylindrical Bubble Through a Cylindrical Vortex

Inspired by the rapid elongation of air columns injected into vortices by dolphins, we present an exact inviscid solution for the axial speed (assumed steady) of propagation of the tip of a semi-infinite cylindrical bubble along the axis of a cylindrical vortex. The bubble is assumed to be held at constant pressure by being connected to a reservoir, the lungs of the dolphin, say. For a given bubble pressure, there is a modest critical rotation rate above which steadily propagating bubbles exist. For a bubble at ambient pressure, the propagation speed of the bubble (relative to axial velocity within the vortex) varies between 0.5 and 0.6 of the maximum rotational speed of the vortex. Surprisingly, the bubble tip can propagate (almost as rapidly) even when the pressure minimum in the vortex core is greater than the bubble pressure; in this case, solutions exhibit a dimple on the nose of the bubble. A situation important for incipient vortex cavitation, and one which dolphins also demonstrate, is elongation of a free bubble, i.e., one whose internal pressure may vary. Under the assumption that the acceleration term is small (checked a posteriori), the steady solution is applied at each instant during the elongation. Three types of behavior are then possible depending on physical parameters and initial conditions: (A) Unabated elongation with slowly increasing bubble pressure, and nearly constant volume. Volume begins to decrease in the late stages. (B1) Elongation with decreasing bubble pressure. A limit point of the steady solution is encountered at a finite bubble length. (B2) Unabated elongation with decreasing bubble pressure and indefinite creation of volume. This is made possible by the existence of propagating solutions at bubble pressures below the minimum vortex pressure. As the bubble stretches, its radius initially decreases but then becomes constant; this is also observed in experiments on incipient vortex cavitation.

Shariff, Karim↗
DOE OSTI2020

Force Balance Model Assessment for Mechanistic Prediction of Sliding Bubble Velocity in Vertical Subcooled Boiling Flow

The bubble sliding after departing from nucleation site is frequently observed in flow boiling systems and, the crucial impact on wall heat transfer has been evidenced through many experiments. As a result, the heat transfer modeling associated with sliding bubble has become one of the subjects of great attention in CFD boiling heat transfer community. The modeling efforts are primarily aimed at improving the existing Heat Flux Partitioning (HFP) model via the implementation of sliding bubble-induced heat transfer. The performance of HFP model depends inherently on the fidelity of sub-models used to predict the fundamental bubble parameters (e.g., bubble departure/lift-off diameter). In the same context, the accurate prediction of sliding bubble parameters (e.g., sliding bubble growth, sliding bubble velocity) is essential to achieving the successful heat transfer modeling associated with sliding bubble. Of many sliding bubble parameters, this paper deals with the sliding bubble velocity. Specifically, the force balance model was assessed in view of the sliding bubble velocity prediction. The parametric effect of key sub-models (e.g., drag force, bubble growth models) used in the force balance equation was investigated. The experimental data from Maity (2000) and Yoo et al. (2016) were used for demonstrating the model predictive performance. It was found that the force balance model proposed in this study was able to predict well the bubble sliding velocity based on the accurate prediction of bubble growth during sliding. The predictive performance was proven with the experimental data measured under various subcooled boiling conditions of both water and refrigerant (NOVEC-7000) flowing upward in vertical channels.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Understanding the effect of minor alloying elements on helium bubble formation in ferritic-martensitic steels

Ferritic-martensitic steels are promising structural materials for advanced nuclear reactors. To minimize long-term radioactivity, reduced-activation ferritic-martensitic steels have been developed by substituting high-activation elements like Ni and Mo with low-activation elements such as W. However, the impact of these alloying modifications on helium bubble formation, which plays a key role in material swelling, remains unclear. Here, in this study, we compared helium bubble formation in ferritic-martensitic steel T91 and reduced-activation ferritic-martensitic steel F82H. Both materials were irradiated with sequential 100 keV, 150 keV, and 200 keV helium ions to a dose of 0.5 dpa and a helium concentration of 9,000 appm at 500°C. The helium bubbles in F82H exhibited a larger average size and a lower density than those in T91, suggesting differences in minor alloying elements may influence the bubble growth. Here, to investigate the effects of these alloying elements, we characterized radiation-induced segregation near bubbles and grain boundaries. Prominent Ni-Mn-Si enriched clusters were found near bubbles in T91, while only Mn-Si enriched clusters were found near bubbles in F82H. In addition, the obvious Cr enrichment near grain boundaries was absent around bubbles in both steels. The different segregation trends among elements revealed the variations in element diffusion mechanisms and the different sink biases between bubbles and grain boundaries. Cr enrichment near grain boundaries is mostly driven by interstitial-mediated diffusion. However, since bubble growth relies on net vacancy flux, vacancy-mediated diffusion plays a dominant role in controlling element segregation near bubbles. Therefore, Cr enrichment was not found near bubbles. Because of preferential vacancy-drag diffusion for Ni, Si and Mn, these elements were enriched near bubbles. Due to the strong binding energies of vacancies with these solute atoms, the vacancy diffusivity can be reduced near these solutes. Therefore, the more prominent Ni-Si-Mn clustered near helium bubbles in T91 lead to stronger suppression of helium bubble growth compared to F82H.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Cavitation inception by the backscattering of pressure waves from a bubble interface

The secondary cavitation that occurs by the backscattering of focused ultrasound from a primary cavitation bubble caused by the negative pressure part of the ultrasound (Maxwell, et al., 2011) might be useful for the energy exchange due to bubble oscillations in High Intensity Focused Ultrasound (HIFU). The present study is concerned with the cavitation inception by the backscattering of ultrasound from a bubble. In the present experiment, a laser-induced bubble which is generated by a pulsed focused laser beam with high intensity is utilized as a primary cavitation bubble. After generating the bubble, focused ultrasound is emitted to the bubble. The acoustic field and the bubble motion are observed with a high-speed video camera. It is confirmed that the secondary cavitation bubble clouds are generated by the backscattering from the laser-induced bubble. The growth of cavitation bubble clouds is analyzed with the image processing method. The experimental results show that the height and width of the bubble clouds grow in stepwise during their evolution. The direct numerical simulations are also conducted for the backscattering of incident pressure waves from a bubble in order to evaluate a pressure field near the bubble. It is shown that the ratio of a bubble collapse time t{sub 0} to a characteristic time of wave propagation t{sub S}, η = t{sub 0}/t{sub s}, is an important determinant for generating negative pressure region by backscattering. The minimum pressure location by the backscattering in simulations is in good agreement with the experiment.

42 ENGINEERING↗
DOE OSTI2022

Plasma breakdown in bubbles passing between two pin electrodes

The ignition of plasmas in liquids has applications from medical instrumentation to manipulation of liquid chemistry. Formation of plasmas directly in a liquid often requires prohibitively large voltages to initiate breakdown. Producing plasma streamers in bubbles submerged in a liquid with higher permittivity can significantly lower the voltage needed to initiate a discharge by reducing the electric field required to produce breakdown. The proximity of the bubble to the electrodes and the shape of the bubbles play critical roles in the manner in which the plasma is produced in, and propagates through, the bubble. In this paper, we discuss results from a three-dimensional direct numerical simulation (DNS) used to investigate the shapes of bubbles formed by injection of air into water. Comparisons are made to results from a companion experiment. A two-dimensional plasma hydrodynamics model was then used to capture the plasma streamer propagation in the bubble using a static bubble geometry generated by the DNS The simulations showed two different modes for streamer formation depending on the bubble shape. In an elliptical bubble, a short electron avalanche triggered a surface ionization wave (SIWs) resulting in plasma propagating along the surface of the bubble. In a circular bubble, an electron avalanche first traveled through the middle of the bubble before two SIWs began to propagate from the point closest to the grounded electrode where a volumetric streamer intersected the surface. In an elliptical bubble approaching a powered electrode in a pin-to-pin configuration, we experimentally observed streamer behavior that qualitatively corresponds with computational results. Optical emission captured over the lifetime of the streamer curve along the path of deformed bubbles, suggesting propagation of the streamer along the liquid/gas boundary interface. Plasma generation supported by the local field enhancement of the deformed bubble surface boundaries is a mechanism that is likely responsible for initiating streamer formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quantitatively Monitoring Bubble-Flow at a Seep Site Offshore Oregon: Field Trials and Methodological Advances for Parallel Optical and Hydroacoustical Measurements

Two lander-based devices, the Bubble-Box and GasQuant-II, were used to investigate the spatial and temporal variability and total gas flow rates of a seep area offshore Oregon, United States. The Bubble-Box is a stereo camera–equipped lander that records bubbles inside a rising corridor with 80 Hz, allowing for automated image analyses of bubble size distributions and rising speeds. GasQuant is a hydroacoustic lander using a horizontally oriented multibeam swath to record the backscatter intensity of bubble streams passing the swath plain. The experimental set up at the Astoria Canyon site at a water depth of about 500 m aimed at calibrating the hydroacoustic GasQuant data with the visual Bubble-Box data for a spatial and temporal flow rate quantification of the site. For about 90 h in total, both systems were deployed simultaneously and pressure and temperature data were recorded using a CTD as well. Detailed image analyses show a Gaussian-like bubble size distribution of bubbles with a radius of 0.6–6 mm (mean 2.5 mm, std. dev. 0.25 mm); this is very similar to other measurements reported in the literature. Rising speeds ranged from 15 to 37 cm/s between 1- and 5-mm bubble sizes and are thus, in parts, slightly faster than reported elsewhere. Bubble sizes and calculated flow rates are rather constant over time at the two monitored bubble streams. Flow rates of these individual bubble streams are in the range of 544–1,278 mm 3 /s. One Bubble-Box data set was used to calibrate the acoustic backscatter response of the GasQuant data, enabling us to calculate a flow rate of the ensonified seep area (~1,700 m 2 ) that ranged from 4.98 to 8.33 L/min (5.38 × 10 6 to 9.01 × 10 6 CH 4 mol/year). Such flow rates are common for seep areas of similar size, and as such, this location is classified as a normally active seep area. For deriving these acoustically based flow rates, the detailed data pre-processing considered echogram gridding methods of the swath data and bubble responses at the respective water depth. The described method uses the inverse gas flow quantification approach and gives an in-depth example of the benefits of using acoustic and optical methods in tandem.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2012

THE FERMI BUBBLES. II. THE POTENTIAL ROLES OF VISCOSITY AND COSMIC-RAY DIFFUSION IN JET MODELS

The origin of the Fermi bubbles recently detected by the Fermi Gamma-ray Space Telescope in the inner Galaxy is mysterious. In the companion paper Guo and Mathews (Paper I), we use hydrodynamic simulations to show that they could be produced by a recent powerful active galactic nucleus (AGN) jet event. Here, we further explore this scenario to study the potential roles of shear viscosity and cosmic-ray (CR) diffusion on the morphology and CR distribution of the bubbles. We show that even a relatively low level of viscosity ({mu}{sub visc} {approx}> 3 g cm{sup -1} s{sup -1}, or {approx}0.1%-1% of Braginskii viscosity in this context) could effectively suppress the development of Kelvin-Helmholtz instabilities at the bubble surface, resulting in smooth bubble edges as observed. Furthermore, viscosity reduces circulating motions within the bubbles, which would otherwise mix the CR-carrying jet backflow near bubble edges with the bubble interior. Thus viscosity naturally produces an edge-favored CR distribution, an important ingredient to produce the observed flat gamma-ray surface brightness distribution. Generically, such a CR distribution often produces a limb-brightened gamma-ray intensity distribution. However, we show that by incorporating CR diffusion that is strongly suppressed across the bubble surface (as inferred from sharp bubble edges) but is close to canonical values in the bubble interior, we obtain a reasonably flat gamma-ray intensity profile. The similarity of the resulting CR bubble with the observed Fermi bubbles strengthens our previous result in Paper I that the Fermi bubbles were produced by a recent AGN jet event. Studies of the nearby Fermi bubbles may provide a unique opportunity to study the potential roles of plasma viscosity and CR diffusion on the evolution of AGN jets and bubbles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Helium Bubble Growth and Retention Characteristics in Aging Palladium Tritide

The helium-related characteristics of aged palladium tritide, including tritium loss, the formation of He bubbles from tritium decay-induced He atoms, and He thermal desorption behavior are computed using a continuum model of bubble evolution that incorporates the log-normal bubble spacing distribution deduced from 3 He nuclear magnetic resonance (NMR) measurements. This spacing distribution produces significant differences between mean (expectation) values of the bubble size and pressure and the average (median) values of these aging characteristics obtained by a previous calculation. The new calculations find much of the He is retained in a large quantity of smaller bubbles with higher bubble pressures, improving the overall stability of the bubbles and extending the age of He retention. By contrast, the integrated tritide material characteristics of swelling and tritium pressure are found relatively insensitive to details of the bubble spacing distribution. Inter-bubble fracture is predicted to begin with small, closely-spaced bubbles and progress to include larger bubbles near the onset of rapid He release. The critical concentration for this onset increases almost linearly with material tensile strength and decreases with shear modulus. Release of He from the aging solid is modeled as a growing network of fractured inter-bubble ligaments. For small particle or small grain material, the physical size of the fracture cluster, relative to the material dimensions, becomes important. In contrast with normal aging, inter-bubble fracture occurring during thermal desorption appears to begin with large, widely-spaced bubbles, indicating elevated temperature techniques may be of limited use for evaluating development of the bubble fracture network.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2000

The Investigation of the Effects of Gravity on Single Bubble Sonoluminescence

In single bubble following it's rapid collapse each cycle of oscillation of an ultrasonic field. Since widely varying length and time scales affect the bubble dynamics and optical emission processes, it is difficult to anticipate the importance of the effects of gravity present for observations on earth. Our bubble is driven in an acoustically resonating cavity at it's first harmonic mode. The acoustical radiation pressure (Bjerknes force) will then keep it suspended in the center near the pressure antinode. When driven in a region where the diffusive processes balance the bubble it acts in a nonlinear but regular way, emitting a short (approx. 200ps) burst of light each acoustic cycle. Balancing the Bjerknes force with buoyancy, as in, we can see that the bubble should be displaced from the velocity node approximately 20m at normal gravity. Therefore, water flows past the bubble at the time of collapse. Gravitation also changes the ambient pressure at the bubble's location, as Delta.P = rho.g.h this gives a change of approximately -0.5% in our experiment when going from 1.8g to 0g. Studies of ambient pressure changes were also done in order to assess these effects. Inside a pressure sealed chamber a spherical glass cell is filled with distilled water which has been degassed to 120mmHg. A bubble is then trapped in the center and driven by a piezoelectric transducer at 32.2kHz attached to the side of the cell. An optical system is then set up to take strobbed video images along and light emission data simultaneously. Temperature, pressure, drive voltage, and listener voltage are also monitored. PMT output in Volts The radii of the bubbles for both experiment s are fit using the Rayleigh-Plesset equation and the acoustic drive amplitude and the ambient bubble radius are found. There is little change in the acoustic drive amplitude as we expect, since we are not varying the drive voltage. However. the ambient bubble radius goes up considerably. These changes (increased light output, increased maximum bubble radius, and increased ambient bubble radius) are also observed when the ambient pressure is varied in the laboratory by an amount similar to that due to gravitation. The changes in the ambient bubble radius and light output with a change in ambient pressure are predicted by the "dissociation hypothesis" and have been observed by other groups in the laboratory. It seems clear that buoyancy's effect on light output and bubble radius, are at best on the same order as the effects of ambient pressure.

Dzikowicz, Ben↗
NASA NTRS1995

The Isolated Bubble Regime in Pool Nucleate Boiling

We consider an isolated bubble boiling regime in which vapour bubbles are intermittently produced at a prearranged set of nucleation site on an upward facing overheated wall plane. In this boiling regime, the bubbles depart from the wall and move as separate entities. Except in the matter of rise velocity, the bubbles do not interfere and are independent of one another. However, the rise velocity is dependent on bubble volume concentration in the bulk. Heat transfer properties specific to this regime cannot be described without bubble detachment size, and we apply our previously developed dynamic theory of vapour bubble growth and detachment to determine this size. Bubble growth is presumed to be thermally controlled. Two limiting cases of bubble evolution are considered: the one in which buoyancy prevails in promoting bubble detachment and the one in which surface tension prevails. We prove termination of the isolated regime of pool nucleate boiling to result from one of the four possible causes, depending on relevant parameters values. The first cause consists in the fact that the upward flow of rising bubbles hampers the downward liquid flow, and under certain conditions, prevents the liquid from coming to the wall in an amount that would be sufficient to compensate for vapour removal from the wall. The second cause is due to the lateral coalescence of growing bubbles that are attached to their corresponding nucleation sites, with ensuing generation of larger bubbles and extended vapour patches near the wall. The other two causes involve longitudinal coalescence either 1) immediately in the wall vicinity, accompanied by the establishment of the multiple bubble boiling regime, or 2) in the bulk, with the formation of vapour columns. The longitudinal coalescence in the bulk is shown to be the most important cause. The critical wall temperature and the heat flux density associated with isolated bubble regime termination are found to be functions of the physical and operating parameters and are discussed in detail.

Buyevich, Y. A.↗