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

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↗

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↗

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗

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↗

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

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

Buyevich, Yu A.↗