Search NASA⌕ Search

SEARCH · Search NASA

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 19 records

Large-scale atomistic simulations of helium-3 bubble growth in complex palladium alloys

Palladium is an attractive material for hydrogen and hydrogen-isotope storage applications due to its properties of large storage density and high diffusion of lattice hydrogen. When considering tritium storage, the material’s structural and mechanical integrity is threatened by both the embrittlement effect of hydrogen and the creation and evolution of additional crystal defects (e.g., dislocations, stacking faults) caused by the formation and growth of helium-3 bubbles. Using recently developed inter-atomic potentials for the palladium-silver-hydrogen system, we perform large-scale atomistic simulations to examine the defect-mediated mechanisms that govern helium bubble growth. Our simulations show the evolution of a distribution of material defects, and we compare the material behavior displayed with expectations from experiment and theory. We also present density functional theory calculations to characterize ideal tensile and shear strengths for these materials, which enable the understanding of how and why our developed potentials either meet or confound these expectations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a mechanistic model for sliding bubbles growth prediction in subcooled boiling flow

In this paper, a new mechanistic model for predicting the growth of sliding bubbles in subcooled boiling flow has been developed. The basis of this model is the assumption that the fundamental heat transfer mechanisms acting on a stagnant bubble growing at a nucleation site can also be applied to the growth of sliding bubbles. The physical assumptions of the existing bubble growth models were modified to reflect the intrinsic physics associated with the sliding bubble growth. The heat transfer mechanisms considered are microlayer evaporation, evaporation of superheated liquid, and condensation. The present model has been successfully validated against a broad set of subcooled flow boiling data, generated by three separate research groups using different working fluids and different heater surfaces. The model successfully reproduced the sliding bubble growth dependencies on mass flux, wall superheat, and liquid subcooling level as observed in the experiments. Finally, on the other hand, in a few cases where the bubble growth inertia was substantially involved in the sliding bubble growth, the predictive ability was limited.

42 ENGINEERING↗

Effective parameterization of phase-field models of fission gas bubble growth

Fission gas bubbles are one of the most important microstructural features of ceramic nuclear fuels. As gas bubbles grow and interconnect, they allow release of gases, with important consequences for fuel performance. Phase-field modeling has been increasingly used to simulate the evolution of fission gas bubble microstructural because of its capability to capture complex microstructural features. However, computational performance limitations have made it difficult to simulate all the defects present in fuels during operation. For this reason, phase-field models have often simulated only vacancies and used multiple approaches to include the effect of vacancy-interstitial recombination and sinks in a simplified way. Here, we compare some of the most prevalent approaches, including source-only and source/sink. The kinetics of bubble growth using these approaches are analyzed analytically, and simulations with these approaches are compared to a full vacancy-interstitial model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a Mechanistic Model for Predicting Sliding Vapor Bubble Growth

A new mechanistic model for predicting the growth of sliding vapor bubble in subcooled boiling flows has been developed. This model employs a physical base that the fundamental mechanisms applied to the bubble growth at a nucleation site can also be applied to the growth of sliding bubble. The heat transfer mechanisms considered are microlayer evaporation, evaporation of superheated liquid, and condensation. Modifications were made to the physical assumptions/reasoning applied to the existing growth models to reflect the intrinsic physics associated with the growth of sliding bubble. The present model has been successfully validated against a broad set of subcooled flow boiling data generated by three separate research groups using different working fluids and heater surfaces. The physical dependency of bubble growth during sliding on the relevant flow boiling parameters (i.e., mass flux, wall superheat, and subcooling level) observed in these experiments was well reproduced by the present model.

42 ENGINEERING↗

Dissociated prismatic loop punching by bubble growth in FCC metals

Abstract Materials performance can be significantly degraded due to bubble generation. In this work, the bubble growth process is elaborated in Cu by atomistic modeling to bridge the gap of experimental observations. Upon continuous He implantation, bubble growth is accommodated first by nucleation of dislocation network from bubble surface, then formation of dissociated prismatic dislocation loop (DPDL), and final DPDL emission in $$\langle 110\rangle$$ ⟨ 110 ⟩ directions. As the DPDL is found capable of collecting He atoms, this process is likely to assist the formation of self-organized bubble superlattice, which has been reported from experiments. Moreover, the pressurized bubble in solid state manifests the shape of an imperfect octahedron, built by Cu $$\{111\}$$ { 111 } surfaces, consistent with experiments. These atomistic details integrating experimental work fill the gap of mechanistic understanding of athermal bubble growth in Cu. Importantly, by associating with nanoindentation testings, DPDL punching by bubble growth arguably applies to various FCC (face-centered cubic) metals such as Au, Ag, Ni, and Al.

36 MATERIALS SCIENCE↗

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The importance of long-timescale simulations for driven systems: An example of He bubble growth at a W GB

Accelerated Molecular Dynamics (AMD) is used to study complex systems under realistic conditions by extending the timescales accessible by Molecular Dynamics. However, some studies rely instead on driving atomic systems harder with higher temperature, faster growth, etc. Here, we study He bubble growth at a W grain boundary as an illustration of harnessing AMD methods to avoid consequences of over-driving the system. The growth mechanisms observed for a He bubble grown under realistic conditions are compared to bubbles-grown orders of magnitude faster, at rates typical of conventional molecular dynamics simulations. We find that progressive growth mechanisms and bubble structures depend on the rate at which the bubble is grown providing further evidence that care must be taken when simulating the dynamics of driven systems such as this one.

36 MATERIALS SCIENCE↗

A thermally-limited bubble growth model for the relaxation time of superheated fuels

Here, we propose a novel approach to evaluate the relaxation time of vapor bubble growth in the context of the flash boiling of a superheated liquid. In alternative to the empirical correlation derived from superheated water experiments almost fifty years ago, the new model describes the thermally-dominated growth of vapor bubbles in terms that are dependent on the local Jakob number (the ratio of sensible heat to latent heat during phase change) and the number density of vapor bubbles. The model is tested by plugging the resulting relaxation time into the Homogenous Relaxation Model (HRM). Flash-boiling simulations carried out with HRM are compared with n-pentane (C 5 H 12 ) injection and boil-off experiments conducted with a real-size, axial-hole, transparent gasoline injector discharging into a constant-pressure vessel. The long-distance microscopy images from the experiments, processed to derive the projected liquid volume (PLV) of the spray, provide a unique set of time-resolved validation data for direct fuel injection simulations. At conditions ranging from flaring to mild and minimal flash boiling, we show that switching to the new relaxation time improves the agreement with the measured PLV profiles with respect to the standard empirical model. Particularly at flaring conditions, the predicted increase in gas cooling caused by rapid vapor production is shown to be more consistent with the observed boil-off.

42 ENGINEERING↗

New helium bubble growth mode at a symmetric grain-boundary in tungsten: accelerated molecular dynamics study

This work, with an emphasis on helium irradiation rates appropriate for fusion-plasma conditions, advances the understanding of helium evolution at grain boundaries in W, an important consideration in the understanding of W as a plasma-facing component. Using accelerated molecular dynamics, helium bubble nucleation and growth at a symmetric Σ5[100](310) tilt grain-boundary in W is studied. The simulations reveal that the growth mode associated with bubble growth at the grain-boundary leads to a suppression of the helium supply to the bubble and hence to arrested growth. Such an unconventional bubble growth mode may dominate in materials with a high density of sinks.

36 MATERIALS SCIENCE↗

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗

Material migration in W and Mo during bubble growth and fuzz formation

Growth of helium (He) induced bubbles and fuzz in tungsten (W) and molybdenum (Mo) is investigated using samples of W films on Mo substrates and Mo films on W substrates exposed to He-containing plasma in the temperature range of 340 to 1075 K, fluence range of 1.0–14 × 10 25 He·m -2 , and incident ion energy of <50 eV. No fuzz (only up to 2 nm diameter bubbles) and no material transport occur in W films at ≤750 K, while precursors-of or fully-developed fuzz and material mixing occur in W and Mo films at ≥800 K. This suggests that fuzz forms in multi-material systems as long as one material meets the conditions for fuzz formation, namely T s / T m ~ 0.27–0.5 where T s and T m are the sample exposure and material melting temperatures, respectively. Larger He bubbles, more material mixing, and further-developed fuzz occur at higher temperature due to increased mobility of He atoms and small He clusters. Accumulation of substrate material at the surface of fuzzy W and Mo thin-film (<80 nm) samples suggests fuzz growth by material transport from the bubble layer in the bulk up to the fiber tip, likely by a two-step process: (i) diffusion of punched dislocation loops in the bulk toward the fuzz base and (ii) diffusion of adatoms along the fuzz base and fiber surface (with effective transport of adatoms upwards due to trapping of adatoms at curved surfaces of fiber tips and/or due to the continuous generation of adatoms at the fuzz base). While the bubble size and fuzz thickness increase with reduced W concentration in Mo thin-film samples at 838 K likely due to an increase in trap mutation and dislocation loop punching in Mo compared to W, the fuzz thickness decreases with reduced W concentration at 1075 K despite an increase in the bubble size likely due to slower diffusion of interstitial loops in Mo.

Physics↗

Phase-field simulations of fission gas bubble growth and interconnection in U-(Pu)-Zr nuclear fuel

Abstract The growth and interconnection of fission gas bubbles in the hotter central regions of U-(Pu)-Zr nuclear fuel has been simulated with a phase-field model. The Cahn-Hilliard equation was used to represent the two-phase microstructure, with a single defect species. The volume fraction of the bubble phase and surface area of the bubble-matrix interface were determined during growth and interconnection. Surface area increased rapidly during the initial stages of growth, then slowed and finally decreased as bubble interconnection began and coarsening acted to reduce surface area. The fraction of the bubbles vented to a simulation domain boundary, f V , was quantified as a measure of the microstructure’s interconnectivity and plotted as a function of porosity p . The defect species diffusivity was varied; although changes in diffusivity significantly affected the microstructure, the plots of f V vs. p did not change significantly. The percolation threshold p c was calculated to be approximately 0.26, depending on the assumed diffusivity and using an initial bubble number density based on experimental observations. This is slightly smaller than the percolation threshold for continuum percolation of overlapping 3D spheres. The simulation results were used to parameterize two different engineering-scale swelling models for U-(Pu)-Zr in the nuclear fuel performance code BISON.

Aagesen, Larry K. (ORCID:000000034936676X)↗

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

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.

36 MATERIALS SCIENCE↗

Theoretical Model of Helium Bubble Growth and Density in Plasma-Facing Metals

We present a theoretically-motivated model of helium bubble density as a function of volume for high-pressure helium bubbles in plasma-facing tungsten. The model is a good match to the empirical correlation we published previously [Hammond et al., Acta Mater. 144, 561–578 (2018)] for small bubbles, but the current model uses no adjustable parameters. The model is likely applicable to significantly larger bubbles than the ones examined here, and its assumptions can be extended trivially to other metals and gases. We expect the model to be broadly applicable and useful in coarse-grained models of gas transport in metals.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Incorporating bubble growth volume feedback to improve simulation of the response of a structure containing liquid and gas to sudden energy input

The SNS target module is a stainless-steel structure that contains and directs mercury. The mercury is struck with short, intense proton pulses to create neutrons. The pulses deposit energy and cause loads on the target structure and mercury cavitation. Helium bubbles are introduced into the mercury to reduce loads and cavitation. This gas complicates the prediction of the target module’s physical response. The current method used to simulate the structural response to the proton pulse incorporates a simple but effective approach to approximate mercury cavitation effects. This method cannot account for the changes from injected non-condensable gas. The ability to model the target vessel response is a prerequisite for estimating its life. Known approaches for simulating a bubbly fluid mixture’s response are computationally expensive and impractical for detailed engineering models. A constitutive model is proposed for target mercury with injected gas bubbles, with the intent of providing more accurate mercury–vessel response behavior without the computing costs required to track individual bubbles. The model’s assumptions and a computational model embodiment are described. The bubbly mercury computational model is subjected to test problems to assess its potential utility. Here, the model was able to better predict the effect of gas on wave propagation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations

We investigate helium flux effects on helium transport and surface evolution in plasma-facing tungsten using molecular dynamics. The simulations span two orders of magnitude, from ITER-relevant levels to those more typical of simulations published to date. Simulation times of up to 2.5 µs (corresponding to actual fluences of 1.4 x 10 21 m -2 ) are achieved, revealing concerted bubble-bursting events that are responsible for significant and very sudden changes in surface morphology. The depth distribution of helium depends very strongly on helium flux: helium self-trapping becomes more probable near the surface at high flux, and a layer of near-surface bubbles forms. Helium retention prior to the onset of bubble bursting is also substantially lower at low flux than it is at high flux. Surface features at low fluence are correlated with the positions of bubbles, but at high fluence, bubbles tend to coalesce, venting to the surface at one or more locations and leaving large interconnected cavities below the surface. Ruptured bubbles may serve as pathways deeper into the material, allowing helium to bypass the layer of near-surface bubbles and fill deeper, potentially much larger, bubbles that can produce more substantial surface features. Deeper bubbles also emit prismatic dislocation loops that can fill in cavities closer to the surface. Our results suggest that nearly all molecular dynamics simulations published to date are hampered by finite-size effects, and that helium flux is a very important parameter in determining the behavior of helium in plasma-facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗