Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bubble growth”

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 163 records · Page 9

Analysis of low Reynolds number separation bubbles using semiempirical methods

The formation and growth of transitional separation bubbles can significantly affect boundary-layer development on airfoils operating at low chord Reynolds numbers. Of primary concern is the change in boundary-layer thickness between laminar separation and turbulent reattachment. This can be estimated using semiempirical methods, such as the one devised by Horton (1968), which are based on solutions to the integral forms of the boundary-layer equations. The applicability of these methods at low Reynolds numbers was investigated using hot-wire measurements of bubbles formed on an NACA 66(3)-018 airfoil at chord Reynolds numbers of 50,000-200,000. The momentum thickness growth between separation and transition was found to be similar to that predicted for a laminar half-jet and appears to be influenced by the momentum thickness Reynolds number at separation. This parameter also was found to have a noticeable effect on the Reynolds number based on the length of a bubble's laminar portion.

Schmidt, Gordon S.↗

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↗

A Generalized Eulerian-Lagrangian Analysis, with Application to Liquid Flows with Vapor Bubbles

Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.

Dejong, Frederik J.↗

He bubble-induced phase transformation of W grain boundaries revealed by accelerated molecular dynamics

The growth of He bubbles and the resulting impact on the microstructural evolution of W are of paramount importance to the plasma-facing materials community due to the application of W in Tokamak fusion reactors. Using accelerated molecular dynamics (AMD) techniques, we outline the structural evolution of grain boundaries (GBs) caused by growing He bubbles. It is discovered that when an alternative, low energy, high density GB structure or phase is available, He bubbles can induce a progressive phase transformation of the GB to the higher density phase by the continual nucleation of W Frenkel pairs. Here, we find that the resulting W self-interstitials migrate to sites at the GB which are structurally related to the higher density phase, leading to the transformation. We discuss the implications of this progressive microstructural evolution on the growing He bubble and consider in general how He bubbles will impact the structural evolution of an arbitrary W GB. These findings of GB phase transformation are predicted to impact other damage events in W such as recrystallization, GB migration and defect segregation which must take these findings into account in order to accurately simulate a realistic W microstructure and hence extract experimentally meaningful data.

36 MATERIALS SCIENCE↗

Physical data measurements and mathematical modelling of simple gas bubble experiments in glass melts

In this work consideration is given to the problem of the extraction of physical data information from gas bubble dissolution and growth measurements. The discussion is limited to the analysis of the simplest experimental systems consisting of a single, one component gas bubble in a glassmelt. It is observed that if the glassmelt is highly under- (super-) saturated, then surface tension effects may be ignored, simplifying the task of extracting gas diffusivity values from the measurements. If, in addition, the bubble rise velocity is very small (or very large) the ease of obtaining physical property data is enhanced. Illustrations are given for typical cases.

Weinberg, Michael C.↗

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

The experimental and analytical work that was done to establish justification and feasibility for a shuttle middeck experiment involving mass transfer between a gas bubble and a liquid is described. The experiment involves the observation and measurement of the dissolution of an isolated immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble were successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration was accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model was developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity.

Dewitt, Kenneth J.↗

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

This paper will describe the experimental and analytical work that has been done to establish justification and feasibility for a Shuttle mid-deck experiment involving mass transfer between a gas bubble and a liquid. The experiment involves the observation and measurement of the dissolution of an isolated, immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble have been successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration has been accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model has been developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity. The results will yield information on transport under conditions of pure diffusion.

Dewitt, Kenneth J.↗

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗

Radiation damage effects in beryllium for next generation neutrino beam targetry (Final Technical Report)

Current and future high-power accelerators put severe requirements on materials used for target and beam windows and target facilities have been recognized as a critical challenge in development of future particle accelerators. In accelerators, window and target materials are exposed to extreme conditions, which include bombardment with very high energy protons (1- 100 GeV) and thermomechanical shock waves. Radiation can cause direct damage in the material, and it leads to production of transmutation products (especially helium), both phenomena having a potential adverse effect on the stability and durability of the target/window material. At high enough temperatures, He can aggregate to form gas bubbles, which in turn cause significant dimensional changes (swelling), enable easy crack propagation, and eventually cause failure by fracture. On the other hand, if the temperature is too low, radiation damage accumulates in the form of internal defects (e.g., dislocations), leading to hardening and a decreased ductility of the material. In this project, we will focus on beryllium since it is considered to be one of the candidate materials for beam windows and targets in the next-generation proton accelerators, e.g., the Long Baseline Neutrino Facility (LBNF). Radiation effects in Be have been studied in the context of nuclear fusion reactor applications. However, key differences exist between reactor and accelerator conditions, including neutron vs. proton irradiation, continuous vs. pulsed beam flux, much higher energies of bombarding particles in accelerators, and higher operating temperatures for typical reactors. For example, the impact of beam pulsing on the radiation damage and the He bubble kinetics is largely unknown. While results obtained on Be from fusion research might not be directly transferrable to understanding target materials, there is an opportunity to bring state-of-the-art tools from materials research in nuclear reactors to aid design of target and beam window materials in high-power accelerators. To this end, the overarching goal of this project are to develop an experimentally-validated computational framework capable of predicting radiation damage evolution in beryllium relevant to beam window and target conditions, focusing on He bubble formation and growth as a function of irradiation temperature. Our model will be based on the cluster dynamics formalism, where size distribution of defects and He bubbles is simulated as a function of time, temperature, and radiation dose. Parameters for the model will be taken from published experiments and from high-fidelity atomistic simulations proposed in this project. In addition, we will carry out a series of targeted ex-situ and in-situ dual-beam experiments using low-energy protons to provide critical data for validation of the model on the effects of radiation on He clustering, He bubble distribution, and dislocation loop density/size in proton irradiated Be.

36 MATERIALS SCIENCE↗

Calibration of the Diffusivity Predictions of Centipede Using Approximate Bayesian Computation and Applications in Nyx (Engineering Scale) and Xolotl-MARMOT (Meso-Scale) Simulations

Fission gas evolution and release in UO 2 nuclear fuel are important fuel performance metrics and occur in several distinct stages: 1) nucleation, growth and resolution of intra-granular bubbles, 2) diffusion to grain boundaries and 3) nucleation and growth of bubbles at grain boundaries, which eventually form a connected network (percolation) enabling release of gas from grain boundaries through connections to triple junctions, grain edges or free surfaces. The NE-SciDAC project is developing several computational tools to model this problem, which are connected in a hierarchical multi-scale framework. The information transfer in the multi-scale framework is a critical step that, in addition to best-estimates, should include uncertainty quantification. Despite taking a first-principles multi-scale approach, there is a need to perform parameter calibration to ensure consistency with available experimental data. In the present study, uncertainty quantification (UQ) and parameter calibration is demonstrated for one of the lower length scale codes in the multi-scale framework (Centipede) and then the results, including instances of the propagated uncertainties, are used in other codes within the framework, specifically Nyx and Xolotl-MARMOT. We calibrated the model parameters in Centipede, a computer code used to predict diffusivities of uranium (U) and xenon (Xe) in the context of the simulation of fission gas in uranium oxide (UO 2 ) nuclear fuel. The Centipede code depends on 183 parameters, all of which are subject to uncertainty. The three data sets used in our calibration effort are taken from the literature. This data is available as a set of measurements, including measurement errors. Our goal is to calibrate a statistical model that predicts both the value of the measurement and the uncertainty associated with the measurement. We perform a Bayesian calibration of the model parameters using a dedicated approximate Bayesian computation (ABC) likelihood function. To avoid excessive computational costs, we replace the expensive Centipede simulation code by a higher-order surrogate model, constructed using only the 9 most important parameters. These important parameters are identified by a preliminary global sensitivity analysis (GSA) study. Among the important parameters are T0 (the temperature at which UO 2 is perfectly stoichiometric) and Hf_pO2 (the temperature dependence of the oxygen (O) partial pressure) that should be considered as operating conditions to be estimated along with the other parameters. We consider two different cases: one where we define one set of these operating conditions for all data sets, and one where we define distinct operating condition parameters for each data set. The Xe diffusivities predicted by the latter case show distinct features that could not be observed in the former. Next, we use the diffusivity predictions by Centipede as input to Nyx, a reduced order fuel performance code focused on gas behavior alone, in order to estimate quantities associated with inter-granular bubble formation at conditions specified by the experiments. Finally, the diffusivities obtained from the calibrated Centipede runs were used in coupled Xolotl-MARMOT simulations of intra- and inter-granular gas evolution. The results are compared to simulations using the baseline diffusivities from Turnbull et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical Modeling of Three-Dimensional Fluid Flow with Phase Change

We present a numerical method to compute phase change dynamics of three-dimensional deformable bubbles. The full Navier-Stokes and energy equations are solved for both phases by a front tracking/finite difference technique. The fluid boundary is explicitly tracked by discrete points that are connected by triangular elements to form a front that is used to keep the stratification of material properties sharp and to calculate the interfacial source terms. Two simulations are presented to show robustness of the method in handling complex phase boundaries. In the first case, growth of a vapor bubble in zero gravity is studied where large volume increase of the bubble is managed by adaptively increasing the front resolution. In the second case, growth of a bubble under high gravity is studied where indentation at the rear of the bubble results in a region of large curvature which challenges the front tracking in three dimensions.

Esmaeeli, Asghar↗

Barophysiology and Biophysics

Decompression is an important aspect of extravehicular activity (EVA). Errors can result in decompression sickness (DCS) if the protective measures are too liberal, while valuable on-orbit time is dissipated in prophylactic methodologies that are excessively conservative. Nucleation is an important consideration in many natural events, and its control is very important in many industrial procedures. The amount of Extravehicular Activity (EVA) that will be required during the construction of the International Space Station exceeds all of the other activity combined. The requirements in astronaut time and consumables (breathing oxygen and air) will be considerable. In an attempt to mitigate these requirements, Project ARGO was investigated in 1990 to investigate the effects of gravitational forces on the musculoskeletal system. This work has led to the present plans for the reduction of prebreathe duration. Over the past decade, research has been directed towards an understanding of the biophysical basis of the formation and growth of the decompression gas phase with the goal of improving the efficiency of the EVA process. In the past, we have direct work towards a more complete understanding of gas bubble formation and growth and exercise-enhanced washout during oxygen prebreathe.

Powell, Michael R.↗

Deciphering the influence of high-flux helium plasma exposure on tungsten nanostructure growth using real-time, in situ spectroscopic ellipsometry

In situ analysis of surfaces during high-flux plasma exposure represents a long-standing challenge in the study of plasma–material interactions. While post-mortem microscopy can provide a detailed picture of structural and compositional changes, in situ techniques can capture the dynamic evolution of the surface. In this study, we demonstrate how spectroscopic ellipsometry can be applied to the real-time characterization of W nanostructure (also known as “fuzz”) growth during exposure to low temperature, high-flux He plasmas. Strikingly, over a wide range of sample temperatures and helium fluences, the measured ellipsometric parameters (Ψ,Δ) collapse onto a single curve that can be directly correlated with surface morphologies characterized by ex situ helium ion microscopy. The initial variation in the (Ψ,Δ) parameters appears to be governed by small changes in surface roughness (<50 nm) produced by helium bubble nucleation and growth, followed by the emergence of 50 nm diameter W tendrils. This basic behavior appears to be reproducible over a wide parameter space, indicating that the spectroscopic ellipsometry may be of general practical use as a diagnostic to study surface morphologies produced by high-flux He implantation in refractory metals. An advantage of the methods outlined here is that they are applicable at low incident ion energies, even below the sputtering threshold. As an example of this application, we apply in situ ellipsometry to examine how W fuzz growth is affected both by varying ion energy and the temperature of the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Dissolution of multicomponent bubbles

The behavior of an isolated, stationary, multicomponent gas bubble in a glassmelt containing several dissolved gases is considered. The relevant mass-transport equations are formulated and calculations are performed for the case of two diffusing gases using a quasi-stationary model and a numerical solution of the exact mass-transfer equations. The results obtained from these two approaches are compared. The factors which govern the dissolution or growth of a bubble are thermodynamic and kinetic in origin. The tendency of a bubble to grow or shrink at long times is controlled by departure from overall equilibrium, whereas the short-time bubble dynamics may be dominated by kinetic effects. As a result of the existence of these dual influences, maxima and/or minima occur in the functional dependence of the bubble radius on time.

Weinberg, M. C.↗

Generalized rate theory for void and bubble swelling and its application to δ -plutonium

A rate theory for void and bubble swelling is derived that allows both vacancies and self-interstitial atoms to be generated by thermal activation at all sinks. In addition, they can also be produced by displacement damage from external and internal radiation. This generalized rate theory (GRT) is applied to swelling of gallium-stabilized δ-plutonium in which α-decay causes the displacement damage. Since the helium atoms produced also become trapped in vacancies, a distinction is made between empty and occupied vacancies. The growth of helium bubbles observed by transmission electron microscopy in weapons-grade and in material enriched with Pu238 is analyzed, using different values for the formation energy of self-interstitial atoms (SIA) and two different sets of relaxation volumes for the vacancy and for the SIA. One set allows preferential capture of SIA at dislocations, while the other set gives equal preference to both vacancy and SIA. It is found that the helium bubble diameters observed are in better agreement with GRT predictions if no preferential capture occurs at dislocations. Therefore, helium bubbles in δ-plutonium will not evolve into voids. Furthermore, the helium density within the bubbles remains sufficiently high to cause thermal emission of SIA. Based on a helium density between two to three helium atoms per vacant site, the sum of formation and migration energies must be around 2.0 eV for SIA in δ-plutonium.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗

Gas Diffusion in Fluids Containing Bubbles

Mathematical model describes movement of gases in fluid containing many bubbles. Model makes it possible to predict growth and shrink age of bubbles as function of time. New model overcomes complexities involved in analysis of varying conditions by making two simplifying assumptions. It treats bubbles as point sources, and it employs approximate expression for gas concentration gradient at liquid/bubble interface. In particular, it is expected to help in developing processes for production of high-quality optical glasses in space.

Zak, M.↗