Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bubble dynamics”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Compositionally-driven surface nanostructuring on refractory compositionally complex alloys under low energy helium bombardment

Additive manufacturing enables user-defined control of the compositional complexity, opening new design spaces for complex concentrated alloys (CCAs). Refractory CCAs may offer enhanced performance in the divertor region of a fusion reactor environment where plasma-facing materials will be subject to high temperatures, low energy He and D particles, and 14 MeV neutrons. In this work, specimens with the nominal composition of NbTaMoTi, NbTaMo, NbTaTi, and NbTa were fabricated via directed energy deposition (DED) then bombarded with 40 eV He ions to a fluence of 2x10 26 m −2 at ∼ 1000 K. Post irradiation, the surface morphology and composition were examined with electron microscopy and x-ray photoelectron spectroscopy to offer information on the spatial distribution of surface nano-structuring. Sub-surface He bubbles driving the nano-structuring were examined with electron microscopy techniques. Analysis indicates the local composition directly influences the He bubble and tendril size, while the region with the highest complexity showed the shortest nano-structuring. Molecular Dynamics simulations complement the experimental results, showing comparable helium bubble growth and migration as a function of compositional complexity. This work demonstrates the compositional dependence of surface nanostructure formation of compositionally complex alloys, important for future design of complex plasma facing materials in fusion reactors.

Complex concentrated alloys↗
DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

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↗
DOE OSTI2020

Suppression of Helium Bubble Nucleation in Beryllium Exposed Tungsten Surfaces.

One of the most severe obstacles to increasing the longevity of tungsten-based plasma facing components, such as divertor tiles, is the surface deterioration driven by sub-surface helium bubble formation and rupture. Supported by experimental observations at PISCES, this work uses molecular dynamics simulations to identify the microscopic mechanisms underlying suppression of helium bubble formation by the introduction of plasma-borne beryllium. Simulations of the initial surface material (crystalline W), early-time Be exposure (amorphous W-Be) and final WBe 2 intermetallic surfaces were used to highlight the effect of Be. Significant differences in He retention, depth distribution and cluster size were observed in the cases with beryllium present. Helium resided much closer to the surface in the Be cases with nearly 80% of the total helium inventory located within the first 2 nm. Moreover, coarsening of the He depth profile due to bubble formation is suppressed due to a one-hundred fold decrease in He mobility in WBe2, relative to crystalline W. This is further evidenced by the drastic reduction in He cluster sizes even when it was observed that both the amorphous W-Be and WBe 2 intermetallic phases retain nearly twice as much He during cumulative implantation studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

An Assessment of Nonhydrostatic and Hydrostatic Dynamical Cores at Seasonal Time Scales in the Energy Exascale Earth System Model (E3SM)

Abstract In global atmospheric modeling, the differences between nonhydrostatic (NH) and hydrostatic (H) dynamical cores are negligible in dry simulations when grid spacing is larger than 10 km. However, recent studies suggest that those differences can be significant at far coarser resolution when moisture is included. To better understand how NH and H differences manifest in global fields, we perform and analyze an ensemble of 28 and 13 km seasonal simulations with the NH and H dynamical cores in the Energy Exascale Earth System Model global atmosphere model, where the differences between H and NH configurations are minimized. A set of idealized rising bubble experiments is also conducted to further investigate the differences. Although NH and H differences are not significant in global statistics and zonal averages, significant differences in precipitation amount and patterns are observed in parts of the tropics. The most prominent differences emerge near India and the Western Pacific in the boreal summer, and the central‐southern Indian Ocean and Pacific in the boreal winter. Tropical differences influence surrounding regions through modification of the regional circulation and can propagate to the extratropics, leading to significant temperature and geopotential differences over the middle to high latitudes. While the dry bubble experiments show negligible deviation between H and NH dynamics until grid spacing is below 6.25 km, precipitation amount and vertical velocity are different in the moist case even at 25 km resolution.

58 GEOSCIENCES↗
DOE OSTI2025

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗
DOE OSTI2024

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Condensation of collective polar vortex modes

The dynamics of extended objects such as domain walls, domain bubbles, vortex structures, etc., can be described by their equations of motion associated with their effective mass and spring constant. Here we analytically derive the equations of motion for the polarization dynamics and elastodynamics for the structural responses of ferroelectric polar vortices, and theoretically extract their effective mass, spring constant, and mode frequencies. We demonstrate two subterahertz phonon modes and predicted their frequencies, both consistent with our recent experimental measurements and phase-field simulations. We show that elastic modulation of the energy function and spring constants leads to a condensation of a collective mode upon a second-order structural transition from symmetric to asymmetric vortices at a critical strain, analogous to the ferroelectric soft phonon mode at a ferroelectric transition. Furthermore, the present work offers a theoretical framework for predicting and manipulating the ultrafast collective dynamics of polar nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Impact of grain boundary and surface diffusion on predicted fission gas bubble behavior and release in UO 2 fuel

In this work, we quantify the impact of grain boundary (GB) and surface diffusion on fission gas bubble evolution and fission gas release in UO 2 nuclear fuel using simulations with a hybrid phase field/cluster dynamics model. Here, we begin with a comprehensive literature review of uranium vacancy and xenon atom diffusivity in UO 2 through the bulk, along GBs, and along surfaces. In our model we represent fast GB and surface diffusion using a heterogeneous diffusivity that is a function of the order parameters that represent bubbles and grains. We find that the GB diffusivity directly impacts the rate of gas release via GB transport, and that the GB diffusivity is likely below 104 times the lower value from Olander and van Uffelen. We also find that the surface diffusivity impacts bubble coalescence and mobility, and that the bubble surface diffusivity is likely below 10 -4 times the value from Zhou and Olander.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Dynamic deuterium retention properties of tungsten measured using laser-induced breakdown spectroscopy

Dynamic deuterium (D) retention properties of tungsten (W) are investigated under steady-state D plasma exposure in the PISCES-A linear plasma device. In contrast to static retention, dynamic retention is quickly released from the material at the termination of the incident plasma flux, and thus in-operando laser-induced breakdown spectroscopy (LIBS) measurements have been conducted during steady-state D plasma exposure. A procedure is, first, established to extract the dynamic retention component from the in-operando LIBS D I 656.1 nm line intensity, which can contain contributions from dynamic and static retention as well as background D/D 2 gas, excited by both steady-state and laser-induced plasmas. Using the developed procedure, the dynamic D retention in W is systematically examined while scanning the following plasma exposure parameters: incident ion energy, E i , sample temperature, T s , incident ion flux, Γ i . No clear E i dependence is seen in the range of E i ~ 45–175 eV, as expected from a small variation in the average implantation depth (~3-6 nm) of D in W. The dynamic retention is found to monotonously decrease with increasing T s from 348 to 573 K, while the total static retention is reported to peak at T s ~ 500–600 K. It is revealed that the dynamic retention linearly increases with increasing Γ i , and then saturates at Γ i ≥ 0.75 × 10 21 m –2 s –1 . Possible physical mechanisms for the observed dependence of the dynamic retention on T s and Γ i are discussed. Furthermore, sequential pure He, followed by pure D, plasma exposures show that the dynamic D retention is not strongly affected by He bubbles in the near-surface region.

super-saturated surface layer↗
DOE OSTI2023

Gaps labeling theorem for the bubble-diamond self-similar graphs

Abstract Motivated by the appearance of fractals in several areas of physics, especially in solid state physics and the physics of aperiodic order, and in other sciences, including the quantum information theory, we present a detailed spectral analysis for a new class of fractal-type diamond graphs, referred to as bubble-diamond graphs, and provide a gap-labeling theorem in the sense of Bellissard for the corresponding probabilistic graph Laplacians using the technique of spectral decimation. Labeling the gaps in the Cantor set by the normalized eigenvalue counting function, also known as the integrated density of states, we describe the gap labels as orbits of a second dynamical system that reflects the branching parameter of the bubble construction and the decimation structure. The spectrum of the natural Laplacian on limit graphs is shown generically to be pure point supported on a Cantor set, though one particular graph has a mixture of pure point and singularly continuous components.

Physics↗
DOE OSTI2024

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↗
DOE OSTI2020

Computationally efficient CFD prediction of bubbly flow using physics-guided deep learning

To realize efficient computational fluid dynamics (CFD) prediction of two-phase flow, a multi-scale framework was proposed in this paper by applying a physics-guided data-driven approach. Instrumental to this framework, Feature Similarity Measurement (FSM) technique was developed for error estimation in two-phase flow simulation using coarse-mesh CFD, to achieve a comparable accuracy as fine-mesh simulations with fast-running feature. In this work, by defining physics-guided parameters and variable gradients as physical features, FSM has the capability to capture the underlying local patterns in the coarse-mesh CFD simulation. Massive low-fidelity data and respective high-fidelity data are used to explore the underlying information relevant to the main simulation errors and the effects of phenomenological scaling. By learning from previous simulation data, a surrogate model using deep feedforward neural network (DFNN) can be developed and trained to estimate the simulation error of coarse-mesh CFD. In a demonstration case of two-phase bubbly flow, the DFNN model well captured and corrected the unphysical “peaks” in the velocity and void fraction profiles near the wall in the coarse-mesh configuration, even for extrapolative predictions. The research documented supports the feasibility of the physics-guided deep learning methods for coarse mesh CFD simulations which has a potential for the efficient industrial design.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

The interaction of edge dislocations with hydrogen-helium bubbles in tungsten

In fusion reactors, plasmas-facing components undergo degradation due to displacement damage and injection of gas impurities. Significant amounts of hydrogen (H) and helium (He) impurities can be introduced into materials through plasma exposure and nuclear transmutation, and their synergistic effects lead to the formation of mixed He/H bubbles under irradiation, changing their hardening mechanisms. In this study, molecular dynamics (MD) models were developed to study the interaction between ½<111> edge dislocations with mixed He/H-bubbles in W. The critical resolved shear stress (CRSS) that is required for dislocation breakaway increases with He/H-to-vacancy ratio at 600 K, indicating strong pinning effects of mixed He/H bubbles. However, as the temperature increases up to 1400 K, the mixed He/H bubbles become unstable as the H atoms are increasingly emitted from the mixed He/H bubbles and migrate into W matrix, which lowers bubble pressure and CRSS. Overall, these results highlight the synergistic effects of He and H on dislocation behavior in irradiated W during deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Sizing and Selection of Pressure Relief Valves for High-Pressure Thermal–Hydraulic Systems

This study covers the critical concerns related to the sizing, selection, installation, maintenance, and testing of pressure safety valves (PSVs). The aim is to ensure the safety of pressurized systems, hydrostatic transmission systems, and hydraulic plants, including process plants, thermal power plants, and nuclear reactor systems. PSVs are devices that ensure the safety and reliability of pressurized vessels, lines, and systems during overpressure events. The task of selecting which PSV features are of greatest value for a specific purpose is complex—especially in the design of a high-pressure experimental thermal–hydraulic facility for hydrostatic and transient testing of the reactor system—when the systems are in the design and development phases and require qualification and demonstration to prove that they have reached a given level of technological readiness. The present study highlights the required steps for users to follow the associated rules, guidelines, and recommendations. As a part of this research, case studies are presented to help readers better understand the applicable strategy and standards. A discussion and a review of PSV performance degradation and failure are summarized to provide a better understanding of varied process applications and conditions, including fluid flow dynamics, boundary-layer formation and pressure drops, gas bubble formation and collapse, geometric configurations, inlet/outlet piping, abrupt pressure fluctuations, and acoustic resonance. Moreover, this study discusses the servicing and testing of PSVs in a multiphase pressurized system. Overall, it provides a basic overview of how PSVs ensure the safety of pressurized systems, supported by case studies and industrial practices.

02 PETROLEUM↗
DOE OSTI2024

Electron‐Assisted Generation and Straight Movement of Skyrmion Bubble in Kagome TbMn 6 Sn 6

Abstract Topological magnetic textures are promising candidates as binary data units for the next‐generation memory device. The precise generation and convenient control of nontrivial spin topology at zero field near room temperature endows the critical advantages in skyrmionic devices but is not simultaneously integrated into one material. Here, in the Kagome plane of quantum TbMn 6 Sn 6 , the expedient generation of the skyrmion bubbles in versatile forms of lattice, chain, and isolated one by converging the electron beam, where the electron intensity gradient contributes to the dynamic generation from local anisotropy variation near spin reorientation transition (SRT) is reported. Encouragingly, by utilizing the dynamic shift of the SRT domain interface, the straight movement is actualized with the skyrmion bubble slave to the SRT domain interface forming an elastic composite object, avoiding the usual deflection from the skyrmion Hall effect. The critical contribution of the SRT domain interface via conveniently electron‐assisted heating is further theoretically validated in micromagnetic simulation, highlighting the compatible application possibility in advanced devices.

Li, Zhuolin↗