Helium diffusion and bubble evolution in tungsten nanotendrils
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
The interaction of monoclonal antibodies (mAbs) with air/water interfaces plays a crucial role in their overall stability in solution. We aim to understand this behavior using pendant bubble measurements to track the dynamic tension reduction and x-ray reflectivity to obtain the electron density profiles (EDPs) at the surface. Native immunoglobulin G mAb is a rigid molecule with a flat, “Y” shape, and simulated EDPs are obtained by rotating a homology construct at the surface. Comparing simulations with experimental EDPs, we obtain surface orientation probability maps showing mAbs transition from flat-on Y-shape configurations to side-on or end-on configurations with increasing concentration. The modeling also shows the presence of β sheets at the surface. Overall, the experiments and the homology modeling elucidate the orientational phase space during different stages of adsorption of mAbs at the air/water interface. These finding will help define new strategies for the manufacture and storage of antibody-based therapeutics.
Morphology is a fundamental attribute when investigating voids and bubbles in UO 2 . Here we use molecular dynamics and Monte Carlo simulations to predict the lowest energy shapes for voids and bubbles in UO 2 . The energies of the {100}, {110} and {111} surfaces have been calculated and used to predict the equilibrium void shape from Wulff construction. This equilibrium shape is compared to low energy faceted voids exhibiting different relative proportions of each family of terminating surfaces. It is found that the equilibrium Wulff shape does not represent the lowest energy morphology for nm void sizes at temperatures between 300 K and 1200 K. Furthermore, the lowest energy faceted voids are slightly more energetically favourable than spherical voids, and as Xe is added, and bubble pressure increases, the faceted morphology becomes even more favourable than the spherical shape.
Numerical simulations of a freely bubbling cylindrical fluidized bed are carried out using a coupled computational fluid dynamics and discrete element method (CFD-DEM) model and compared to recent experimental data. The experiments were conducted using high-resolution and high-frequency magnetic resonance imaging providing high-fidelity data of the bubbling within a central 10 mm slice of the bed. Qualitatively, we nd more regular (i.e., less chaotic) structures in the simulated beds than observed experimentally. Quantitatively, however, the bubble diameter and number of bubbles as a function of height within the bed is predicted well by the base model. Unfortunately, the regularity in the simulations manifests as a considerable discrepancy in the speed of the (dense) emulsion phase. The simulated velocity probability distribution functions show an accumulation of low-speed regions and deficiency of high-speed regions. A simple parametric study of the base model is also carried out considering many of the most common CFD-DEM modeling parameters. It is found that the fluid grid size, geometry resolution, transfer kernel and drag law did not have a significant effect on bubble or particle dynamics.
The mechanisms responsible for the collapse of helium-filled bubbles during the passage of shock waves in monocrystalline copper are revealed. Both internal pressure (caused by pre-existing helium atoms) and bubble size are varied in molecular dynamics simulations to understand the atomistic scale deformation as they are subjected to shock compression at pressures of 48, 123, and 170 GPa, corresponding to particle velocities of 1.0, 2.0, and 2.5 km/s. Both empty and helium filled bubbles serve as dislocation sources, generating intense, localized plastic regions. There are distinct differences in the collapse of empty voids compared to He-filled bubbles, the former requiring less stress and generating a greater density of dislocations for a given shock strength. Additionally, a generalized model for dislocation emission is proposed, where the inclusion of shear stress generated by the helium bubble increases the critical stress to generate dislocations at the defect surface, demonstrating the change in plastic deformation.
The dynamic regime of gas injection through a circular plate orifice into an ideally wetting liquid is considered, when successively detached bubbles may be regarded as separate identities. In normal gravity and at relatively low gas flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole time of its evolution. If the flow rate exceeds a certain threshold value, another stage of the detachment process takes place in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical stem that lengthens as the bubble rises above the plate. The bubble shape resembles then that of a mushroom and the upper envelope continues to grow until the gas supply through the stem is completely cut off. Such a stage is always present under conditions of sufficiently low gravity, irrespective of the flow rate. Two major reasons make for bubble detachment: the buoyancy force and the force due to the momentum inflow into the bubble with the injected gas. The former force dominates the process at normal gravity whereas the second one plays a key role under negligible gravity conditions. It is precisely this fundamental factor that conditions the drastic influence on bubble growth and detachment that changes in gravity are able to cause. The frequency of bubble formation is proportional to and the volume of detached bubbles is independent of the gas flow rate in sufficiently low gravity, while at normal and moderately reduced gravity conditions the first variable slightly decreases and the second one almost linearly increases as the flow rate grows. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.
This report summarizes lower length scale computational research conducted to improve the pulverization criterion for high-burnup UO 2 fuel in the BISON fuel performance code. This research was sponsored by the NEAMS program during FY22. Efforts to improve the model primarily focus on calculating the current pressure of bubbles in the high-burnup structure (HBS) region, as well as the critical bubble pressure at which pulverization occurs in the HBS. The phase-field model for predicting initial bubble pressure in the HBS was improved by implementing a more realistic model for defect production, and by coupling the phase-field model for inter-granular bubble evolution with the spatially-resolved cluster dynamics code Xolotl to simulate intra-granular fission gas evolution. An evolution equation for current bubble pressure was added to BISON, an improvement over the previous model that assumed that bubble pressure was static following HBS formation. To improve calculations of critical bubble pressure for pulverization, 3-D phase-field fracture simulations were performed, and a function for critical bubble pressure was fit—based on the simulation results—to replace the previous function, which had been determined using 2-D simulations. The impacts of these modifications to the existing pulverization assessment cases will be reported in the forthcoming engineering scale milestone on high-burnup UO 2 pulverization. A poromechanics-based approach was used to include the effect of bubble over pressurization on the stress state of the pellet at the engineering scale, and the initial strategy for integration of the pulverization criterion with pellet-scale mechanical degradation using a smeared cracking model was developed. This report also describes the initial research that was conducted to inform a forthcoming transient fission gas release (tFGR) model in BISON. Focusing on the fission gas release caused by pulverization of the outer pellet rim, a model for the amount of fission gas release was developed as a function of fuel porosity, bubble size, and fragment size. Implementation of this model in BISON and its impact of engineering-scale predictions will also be described in the aforementioned engineering-scale milestone on high-burnup UO 2 pulverization.
A model of the dynamic regime of gas injection through a submerged plate orifice into an ideally wetting liquid is developed in the circumstance when successively detached bubbles may be regarded as independent objects. Two major factors favor bubble detachment: the buoyancy force and a force due to the momentum inflow into the bubble with injected gas. In normal and moderately reduced gravity, the first factor dominates. At relatively low flow rates, a growing bubble is modeled as a spherical segment touching the orifice perimeter during the whole period of its evolution till detachment. If the flow rate exceeds a critical value, another stage of bubble evolution occurs in which an almost spherical gas envelope is connected with the orifice by a nearly cylindrical gaseous stem that lengthens as the bubble rises above the plate. The bubble continues to grow until the gas supply through the stem is completely cut off, after which back flow of gas into the stem from the bubble becomes possible. In microgravity, the second factor prevails, and the latter stage is always present irrespective of the flow rate. However, the gas envelope rises and the stem lengthens very slowly. This difference in the underlying physical mechanisms provides for key properties of bubble growth and detachment being drastically different in appreciable and sufficiently reduced gravity. The frequency of bubble formation slightly decreases with and the detachment volume is almost proportional to the gas flow rate in the first case, in accordance with familiar relations. In the second case, the first variable is proportional to the flow rate whereas the second one is independent of it. Effects of other parameters, such as the orifice radius, gas and liquid densities, and surface tension are discussed.
Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic ray proton (CRp), cosmic ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics, and heating, and model their expected cavity-power versus radio-luminosity relation (P cav –L R ). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within ∼30 Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the P cav –L R plane, but the values are broadly consistent with observed ranges for FR-I sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.
In collaboration with Argonne National Laboratory (ANL), Los Alamos National Laboratory (LANL) is assisting in the design and development of portions of the second ANL Bubble Experiment to be performed in late 2020 at ANL. The ANL Bubble Experiment, as called in this report, is a series of direct electron irradiations of a uranyl sulfate solution to produce radiolysis-induced gas bubbles of hydrogen and oxygen. The gas bubbles formed in the solution enhance mixing and heat transfer. The study of the dynamics, shape, and size of radiolysis-induced gas bubbles is of great importance to understand and characterize the thermal and fluid behavior of the solution, especially for the solution based, neutron-induced fission production technique for Mo-99. During the first experiment performed in 2014, ANL’s 35 MeV electron linear accelerator provided average powers of 6, 12, and 15 kW using a rastered beam to homogeneously heat the 15 x 15 x 80-cm uranyl-sulfate solution. Gas bubble size, shape, velocity, solution temperature, and hydrogen and oxygen concentrations were recorded during the irradiations. The details of the experimental setup and results for the first experiment are described in the ANL reports, Design and Construction of Experiment for Direct Electron Irradiation of a Uranyl Sulfate Solution: Bubble Formation and Thermal Hydraulics Studies and Experimental Results for Direct Electron Irradiation of a Uranyl Sulfate Solution: Bubble Formation and Thermal Hydraulics Studies. The purpose of this report is to describe the issues and lessons learned associated with the first experiment and elaborate on designs to improve and obtain more accurate experimental results for the upcoming experiment.
Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.
Not provided.
In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.
SHINE Medical Technologies is planning to use neutron-induced fission in a subcritical low-enriched-uranium uranyl sulfate target solution for production of 99 Mo. During this operation, the solution will undergo self-heating due to fissioning of the uranium, radiolytic decomposition of the water, and circulation due to thermal gradients generated in the solution, and will be cooled by cooling tubes running through the annulus and from cooling outside the annulus. Because the formation of the radiolysis-induced bubbles (H 2 and O 2 ) and their size and dynamics will impact the operational parameters of the liquid target, an understanding of bubble behavior is critical for the ability to predict the behavior of the target solution during this operation. It is also important to be able to predict the thermal gradients and the circulation in the vessel. Researchers at Argonne National Laboratory have designed an experimental setup to study radiolytic gas formation in uranyl sulfate under direct electron beam irradiation and have conducted initial experiments. Results of those experiments provided invaluable information on thermal hydraulic behavior of the solution and some information on bubble formation and behavior, but those initial experiments fell short in the measurements of the gas generation rates and bubble behavior. To address the shortcomings of the original experiment, the irradiation setup was modified to improve our abilities to measure gas generation rates and measure the temperature distribution in the solution with better precision. Modifications to the experimental setup and preparation for the irradiations are described below.
Experimental results and analysis of the effects of boundary layer tripping on dynamic stall of a transiently pitching airfoil are presented. At low Mach numbers, the tripped airfoil exhibits qualitative similarity with the behavior of the untripped airfoil. However, the local supersonic flow at Mach numbers greater than 0.3 is significantly modified by the trip leading to vastly different shock/boundary layer interactions, dynamic stall onset and vortex formation angles. The formation of the laminar separation bubble is found to have a favorable influence in delaying dynamic stall on the untripped airfoil flow. In both Mach number regimes, the tripped flow actually stalls at slightly lower angles of attack. Further experimentation with three trips on an oscillating airfoil showed that the dynamic stall process is very sensitive to the state of the turbulence in the boundary layer. This sensitivity points to a need for finer turbulence modeling techniques for use in dynamic stall computations.
The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.
Nanostructured Ferritic Alloys (NFAs) containing oxide nanoparticles are prospective advanced structural materials in future fusion reactors. Lacking fusion neutron sources with an adequate flux, modeling is critical to understand radiation damage in the NFAs, including helium effects. Machine learning interatomic potentials (MLPs) have been demonstrated to be capable of describing atomic interactions in chemically complex systems with comparable accuracy to density functional theory (DFT) but with much less computational costs. Hence, MLPs are chosen to study He bubble accumulation in chemically complex NFAs. As a preliminary step, we develop a Fe–He potential, based on ~10,000 atomic configurations. The developed MLP accurately predicts the bulk properties of BCC-Fe compared to DFT, including the lattice constant (2.832 Å), elastic constants (c 11 = 271 GPa, c 12 = 141, c 44 = 93), and phonon frequencies (maximum error < 3.5%). The MLP predicts He binding energies in He n V bubbles and He n clusters with a mean absolute error (MAE) of only ~70 meV, which is ~3 to ~5 times smaller than the MAE of binding energies estimated using empirical potentials. The MLP also accurately predicts the migration barrier of interstitial He (64 meV). Furthermore, the MLP correctly predicts the maximum number of He atoms (critical size) at which a He n V bubble (He n cluster) can kick-out a self-interstitial atom (SIA), forming a bigger He n V 2 bubble (He n V bubble). Subsequently, we use the MLP to perform molecular dynamics simulations to investigate the effect of temperature on the critical size of He n V bubbles and He n clusters. We find that the critical size decreases with increasing temperature — indicating that smaller bubbles and clusters can kick-out an SIA at higher temperatures, thereby escalating the He bubble growth. We report the current Fe–He MLP can be further developed to include all chemical interactions in NFAs.
Consideration is given to the time-dependent evolutions of the free surface profile (bubble shapes) of a cylindrical container, partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry in low and microgravity environments. The dynamics of the bubble shapes are calculated for four cases: linear time-dependent functions of spin-up and spin-down in low and microgravity, linear time-dependent functions of increasing and decreasing gravity at high and low rotating cylinder speeds, time-dependent step functions of spin-up and spin-down in low gravity, and sinusoidal function oscillation of the gravity environment in high and low rotating cylinder speeds. It is shown that the computer algorithms developed by Hung et al. (1988) may be used to simulate the profile of time-dependent bubble shapes under variations of centrifugal, capillary, and gravity forces.