Control of bubble ζ-potentials to improve the kinetics of bubble-particle interactions
Not provided.
SEARCH · Search NASA
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.
Not provided.
Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.
Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. Finally, the size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.
The heat transfer effect of sliding bubbles is discussed in two aspects based on new experimental findings: (i) the area of bubble influence and (ii) bubble sliding characteristics influencing the degree of wall heat transfer. One notable finding was that the bubble sliding behavior, characterized by sliding trajectories, played a crucial role in determining the area of bubble influence and the bubble influence factor (K). Also, such sliding characteristic was observed strongly dependent on the sliding bubble size. An empirical model has been proposed based on this finding to evaluate the bubble influence factor (K) for the sliding bubbles. Another experimental finding was that the bubble sliding characteristic, represented by bubble spreading factor (S), had a significant impact on the degree of wall heat transfer induced by sliding bubbles. This is due to the fact that the bubbles spread across the heater width as they slid downstream, which reduced the effective frequency of sliding bubbles per unit wall area. These new findings not only improve our understanding of the heat transfer associated with sliding bubbles, but also are expected to contribute to improving the insight in the existing CFD boiling heat transfer models.
To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized this is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory's Marmot application for phase-field simulation of nuclear materials. Simulations of bubble growth in the high burnup region during steady-state conditions showed that initially overpressurized bubbles decreased in pressure during growth, but still remained above equilibrium pressure. During a LOCA transient, simulations of bubbles in the high burnup region showed that bubble size did not change significantly, and the pressure increase due to the transient was calculated and passed to a phase-field model of fracture.
Micro-bubbles dispersed in bubble column reactors have received great interest in recent years, due to their small size, stability, high gas-liquid interfacial area concentrations and longer residence times. The high gas-liquid interfacial area concentrations lead to high mass transfer rates compared to conventional bubble column reactors. In the present work, experiments have been performed in a down-flow bubble column reactor with micro-bubbles generated and dispersed by a novel mechanism to determine the gas-liquid interfacial area concentrations by measuring the void fraction and bubble size distributions. Gamma-ray densitometry has been employed to determine the axial and radial distributions of void fraction and a high speed camera equipped with a borescope is used to measure the axial and radial variations of bubble sizes. Also, the effects of superficial gas and liquid velocities on the two-phase flow characteristics have been investigated. Further, reconstruction techniques of the radial void fraction profiles from the gamma densitometry's chordal measurements are discussed and compared for a bubble column reactor with dispersed micro-bubbles. The results demonstrate that the new bubble generation technique offers high interfacial area concentrations (1,000 to 4,500 m 2 /m 3 ) with sub-millimeter bubbles (500 to 900 µm) and high overall void fractions (10% – 60%) in comparison with previous bubble column reactor designs. The void fraction data was analyzed using slip velocity model and empirical correlation has been proposed to predict the Sauter mean bubble diameter.
Not Available
A monolithic fuel design based on a U-Mo alloy has been selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). An issue with U-Mo monolithic fuel is the large amount of swelling that takes place during operation. The accurate prediction of fuel evolution under irradiation requires implementation of correct thermodynamic properties into mesoscale and continuum level fuel performance modeling codes. However, the thermodynamic properties of the fission gas bubbles (such as the relationship among bubble size, equilibrium Xe concentration, and bubble pressure) are not well known. This work studies Xe bubbles in ?U-Mo from a diameter of 3 nm up to 8.5 nm and from 400 K up to 700 K. The energetic relationship of Xe bubbles with regard to voids and Xe substitutional atoms is described. The transition is also determined for when a bubble becomes over-pressurized. Finally, an equation of state is fit to the pressure as a function of molar volume and temperature.
Due to the large formation energy of vacancies and noble gas atoms in the form of interstitials or substitutional atoms in nuclear fuel (UO2), the thermodynamic equilibrium concentrations of these species are very low in the nuclear fuel matrix even at very high temperatures, which imposes difficulties upon the quantitative study of bubble evolution via the phase-field method. In this study, a quantitative phase-field model is proposed to deal with this problem. The system’s free energy density is derived according to the principles of thermodynamics, with consideration of the elastic effect and with the use of real material parameters from experiments. The model is useful for the study of the kinetics of gas bubble growth with very dilute concentrations of vacancy and gas atoms in the matrix. This model is applied to study single bubble growth and multiple bubble growth under various concentrations of vacancy and gas atoms and at various temperatures. The elastic effect and the effects of the generation rate of vacancies and gas atoms on bubble growth are analyzed.
Not Available
Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.
Irradiation examination shows that gas bubble swelling kinetics is much faster after irradiation-induced recrystallization than that prior recrystallization in UMo fuels. It implies that gas bubbles in coarse grains and small recrystallized grains have different growth behavior. In this work, a phase-field model of gas bubble evolution integrating microstructure dependent cluster dynamics has been developed, for the first time, to study the gas bubble swelling behavior in the recrystallization zone of UMo fuels. Generation, diffusion, reaction, sink, emission and clustering of vacancies and interstitials are described by the cluster dynamics model while a phase-field model is used to describe the evolution of non-equilibrium gas bubbles including nucleation and growth. With the coupled model, the effect of defect generation rate, clustering rate, interstitial emission and sink rates on grain boundaries on the gas bubble evolution are systematically simulated. A set of model parameters (defect generation rate, clustering rate, interstitial emission and sink rates) is determined by comparing measured and simulated gas bubble swelling kinetics. The results demonstrate that interstitial clustering is one of the important physical mechanisms which results in a fast gas bubble swelling kinetics in the recrystallization zone. The developed model can also be extended to study the associated growth of defect and second phase precipitates often observed in irradiated materials.
Not Available
Vertical vibration is known to cause bubble breakup, clustering and retardation in gas-liquid systems. In a bubble column, vibration increases the mass transfer ratio by increasing the residence time and phase interfacial area through introducing kinetic buoyancy force (Bjerknes effect) and bubble breakup. Previous studies have explored the effect of vibration frequency (f), but minimal effort has focused on the effect of amplitude (A) on mass transfer intensification. Thus, the current work experimentally examines bubble size, void fraction, and mass transfer in a bubble column under relatively high amplitude vibration (1.5 mm < A <9.5 mm) over a frequency range of 7.5–22.5 Hz. Results of the present work were compared with past studies. The maximum stable bubble size under vibration was scaled using Hinze theory for breakage. Results of this work indicate that vibration frequency exhibits local maxima in both mass transfer and void fraction. Moreover, an optimum amplitude that is independent of vibration frequency was found for mass transfer enhancements. Finally, this work suggests physics-based models to predict void fraction and mass transfer in a vibrating bubble column.
Abstract not provided
Not provided.
The relationship between noble metal phase particles and fission gas bubble production in used nuclear fuel is described. The majority of Te atoms within noble metal phase undergo radioactive decay to form stable Xe within a few hours after particle formation. This results in the production of clusters of Xe atoms contained within the solid metal matrix exhibiting an equivalent gas bubble pressure approaching 1 GPa. These high pressure bubbles are stabilized by the UO2 within the bulk of the fuel. However, when these bubbles form near the fuel / cladding interface, in combination with local and temporal damage caused by fission recoil, they are capable of overcoming the fracture strength of the UO2 and rupturing catastrophically. The force of the resulting bubble rupture is sufficient to eject noble metal phase particles several microns into the cladding. This proposed mechanism explains the observance of noble metal phase in cladding and is consistent with a host of morphological features found near the fuel / cladding interface.