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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupled Experimental/Computational Investigation of the Dynamics of Interacting Magnetized Plasmas

The interaction, or interpenetration, of magnetized plasmas of different density and/or pressure occurs in a wide variety of natural and man-made systems. Such systems include extragalactic jets propagating into the intergalactic medium, solar coronal mass ejections into background solar wind, compact toroid (CT) fueling of magnetic fusion plasmas, and jets of capsule shell impurities into DT fusion fuel, which can lead to enhanced impurity mix in inertial fusion implosions. These plasmas may take the form of jets, with open, helical magnetic structure, or plasma “bubbles” with closed magnetic fields (B-fields), such as spheromaks or CT’s. Such structures, both open and closed B-field cases, can transport heat, particles and magnetic flux or magnetic helicity into background plasma regions. For example, the origin of extragalactic magnetic fields may be due, at least in part, to transport by astrophysical jets. The goal of this proposed work was to elucidate the detailed plasma and magnetic field dynamics of high-density plasma jets (open B-field) and bubbles (closed B-field) propagating into lower density background magnetized plasma through controlled laboratory experiments and closely coupled nonlinear MHD modeling. These experiments were conducted in the HelCat (Helicon-Cathode) linear plasma device at the University of New Mexico (UNM). Plasma jets and bubbles were launched via an existing compact coaxial plasma gun, mounted on the HelCat device. This gun produced plasmas tens of cm in scale and lasting tens of microseconds, thereby allowing detailed multipoint, space- and time-resolved measurements to be made routinely. The experiments were directly modeled using the extended magnetohydrodynamic (XMHD) PERSEUS code, developed at Cornell University [23,24]. Both experimental and numerical modeling work are ongoing. The main results to date are reported here. Additional supplemental funding for one year (8/1/2019 – 7/31/2020) supported numerical investigation of photoionization processes important in many low temperature plasmas, including the HelCat device. Initial results of this modeling work is also reported.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

In-situ ion irradiation of a spent UO2 fuel: evolution of fission products and nanograins with radiation dose

This study investigates the behavior of fission gas bubbles, five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) and nanograins in spent uranium dioxide (UO2) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles, 5MPs and nanograins in UO2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. The study also reveals nanograin coarsening under irradiation. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. This study provides critical insights into the dynamic behavior of fission products and nanograins under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

In-situ ion irradiation of fission products in a spent UO 2 fuel

This study investigates the behavior of fission gas bubbles and five metal precipitates (5MPs) (Mo, Ru, Rh, Tc, Pd) in spent uranium dioxide (UO 2 ) fuel under various ion irradiation doses, temperatures, and flux conditions. Utilizing in-situ ion irradiation and advanced transmission electron microscopy, we analyzed the evolution of fission gas bubbles and 5MPs in UO 2 samples from the Belgium Reactor 3 (BR-3). Our findings reveal significant shrinkage of fission gas bubbles and 5MPs with increasing irradiation dose, accompanied by a decrease in pair density. We demonstrate that ion irradiation induces a homogeneous re-solution process where individual atoms are ejected from bubbles and precipitates, leading to their dissolution and subsequent re-precipitation in the matrix. In conclusion, this study provides critical insights into the dynamic behavior of fission products under irradiation, facilitating the development of predictive models and contributing to the optimization of nuclear fuel performance and safety.

Fission products↗
DOE OSTI2025

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Active microrheology and dynamic phases for pattern-forming systems with competing interactions

We consider the driven dynamics of a probe particle moving through an assembly of particles with competing long-range repulsive and short-range attractive interactions, which form crystal, stripe, labyrinth, and bubble states as the ratio of attraction to repulsion is varied. We show that the probe particle exhibits a depinning-like threshold from an elastic regime, where the probe particle is trapped by interactions with the other particles, to a plastic flow regime, where the probe particle can break bonds in the surrounding medium. For a fixed particle density, the depinning threshold and sliding velocity of the probe particle vary nonmonotonically as the attraction term is increased. A velocity minimum appears near the crystal to stripe crossover, and there is a significant increase in the depinning threshold in the bubble regime when the probe particle is strongly confined inside the bubbles. For fixed attractive interaction but increasing particle density, the behavior is also nonmonotonic, and there are jumps and drops in the velocity and depinning threshold corresponding to points at which the system transitions between different structures. There are also several distinct flow states that can be characterized by the amount of plastic deformation induced by the probe particle in the surrounding medium. Each flow state generates a different amount of effective drag on the probe particle, and there can be jumps in the velocity-force curve at transitions between the states. In conclusion, we also find that when oriented stripes are present, the probe particle can move along the stripe in an edge transport state that has a finite Hall angle.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

risingMicrobubbleLattice

This code is a simulation case to be run with OpenFOAM, an open-source computational fluid dynamics software. A gmsh mesh file is also included. Specifically, this simulation demonstrates the transport of a single microbubble rising through an ordered lattice due to an applied flow. The bubble deforms as it squeezes through the pores of the lattice. Is

Guo, Jack [Lawrence Livermore National Laboratory ↗
DOE OSTI2025

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗
DOE OSTI2024

Finite Bubble Statistics Constrain Late Cosmological Phase Transitions

We consider first order cosmological phase transitions (PTs) happening at late times below standard model temperatures T PT ≲ GeV. The inherently stochastic nature of bubble nucleation and the finite number of bubbles associated with a late-time PT lead to superhorizon fluctuations in the PT completion time. We compute how such fluctuations eventually source curvature fluctuations with universal properties, independent of the microphysics of the PT dynamics. Using cosmic microwave background (CMB) and large scale structure measurements, we constrain the energy released in a dark-sector PT. For 0.1 eV ≲ T PT ≲ keV this constraint is stronger than both the current bound from additional neutrino species Δ⁢N eff , and in some cases, even CMB-S4 projections. Future measurements of CMB spectral distortions and pulsar timing arrays will also provide competitive sensitivity for keV ≲ T PT ≲ GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗
DOE OSTI2025

Dynamics of Downdrafts Around a Growing Convective Cloud: A Numerical Study

We examine the dynamics of cloud-edge downdrafts over the growth phase of isolated cumuli, combining Eulerian and Lagrangian analyses. As in previous studies, our results show that growing cumuli are surrounded by downdrafts linked to cloud-scale quasi-toroidal circulations at all times at middle and upper cloud levels consistent with the thermal chain description of convective clouds. These toroidal circulations are responsible for the most intense cloud-edge downdrafts in our simulations. In the upper cloud half, roughly 30%–50% of the upward mass flux is typically compensated within a radius of about twice the updraft radius in quasi-laminar simulations forced by a warm bubble in an initially quiescent flow. In a turbulent cloud forced by surface fluxes, this compensation fraction is around 10%–30% over the same region. In contrast to the buoyancy-centered view of subsiding shells, Eulerian and Lagrangian vertical momentum budget analyses show that the most intense cloud-edge downdrafts in the turbulent setup, and after spin-up of the toroidal circulation in the quasi-laminar experiments, are predominantly mechanically forced (i.e., driven by dynamic pressure accelerations). This is consistent throughout the entire growth phase of the cumulus clouds and across tests with varying assumptions, including drier and moister environments. Despite dynamic pressure perturbations being the main driver of toroidal downdrafts, the downdraft speed (relative to the corresponding updraft velocity) exceeds the prediction of the non-buoyant Hill's spherical vortex—a simple model frequently used for cloud-scale circulations—by more than 30%.

Pardo, Lianet Hernández [Goethe Univ., Frankfurt (↗
DOE OSTI2025

Influence of Pore Water on the Fracture of Silica Sand at Particle Scale

The fracture of particles has a significant influence on the engineering behavior of granular materials. There are no reported conclusive answers in the literature to explain the influence of pore water on the fracture of sand, gravel, railroad ballast rock, aggregates, rockfill, and so on. Here, in this paper, two novel techniques were adopted to investigate the influence of pore water on the fracture of natural silica sand at the particle scale. The three-dimensional (3D) synchrotron microcomputed tomography (SMT) imaging technique was used to acquire and analyze multiple 3D images of specimens composed of wetted silica sand that were subjected to confined one-dimensional (1D) compression loading up to the fracture stage. A few representative particles were identified and digitally separated from the liquid water and gas phases. The 3D SMT images offer clear experimental evidence of the phase transition of water from liquid to gas within the opening cracks. In addition, a computational fluid dynamics numerical simulation framework for water flow into an opening crack was adopted and the results show that the fast opening of the crack induced cavitation, water phase transition, and water hammer effects. The results of the numerical simulations agree with the confined 1D compression experiments where the onset of the bubbles (gas phase) within the crack is mainly generated during the secondary cavitation stage, and the generated bubble near the crack mouth opening resulted from the primary cavitation, which was directly caused by the crack expansion. The results reported in this paper present a new cavitation phenomenon that occurs during particle fracture, which offers a physics-based explanation of why water-wetted sand fractures at smaller stresses than dry and can pave the way for more in-depth future studies on the fracture of water-wetted sand.

cavitation↗
DOE OSTI2026

Time-resolved unsteady shock motion and spectral analysis of a hypersonic shock-wave/boundary-layer interaction over a double-cone

Direct experimental measurements of double-cone shock-wave/boundary-layer interaction dynamics at Mach 7.2 were collected using a schlieren-based shock-detection algorithm at the University of Texas at San Antonio. The points of flow separation and reattachment were spatially and spectrally evaluated at 300 kHz for time-resolved characterization of the unsteadiness of this flow field. The separation bubble was found to rapidly expand and contract about a large mean length. Spectral analysis of the velocity fluctuations of this motion revealed this behavior to be a turbulent phenomenon, with spectral energy trends consistent with the Kolmogorov theory for the turbulent energy cascade.

Glasby, Ryan [ORNL]↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗