Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BUBBLE”

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

Molecular dynamics study of interstitial He clusters in nickel

This study presents a molecular dynamics analysis focusing on the behavior of interstitial helium (He) clusters in nickel (Ni), examining their formation, stability, and migration energetics. Consistent with previous research, we found that the binding energies of interstitial helium within a helium cluster are positive and increase with the cluster size, indicating a preference for helium atoms to cluster together. However, our findings also reveal that while the formation energy increases monotonically with cluster size, the increase in binding energy is non-monotonic. Importantly, small He clusters were observed to be thermally unstable at reactor operational temperatures (approximately 600 K), with the He 2 cluster exhibiting instability even at room temperature. With a binding energy of 0.49 eV for a He 4 cluster, we hypothesize that for helium bubbles to form via homogeneous nucleation (i.e., through trap mutation) at reactor operating temperatures, the helium concentration must be high enough to facilitate the formation of helium clusters of at least size 4 or larger. As expected, interstitial helium and small helium clusters are highly mobile. This mobility was observed not only at room temperature but also at temperatures as low as approximately 200 K. Furthermore, the mean squared displacement method has been utilized to determine the migration barriers and the corresponding prefactors for clusters ranging from He 1 to He 6

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Improved departure from nucleate boiling prediction in rod bundles using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is a regulatory limit for the licensing of pressurized water reactors (PWRs) worldwide. Despite the abundance of predictive tools available to the reactor thermal-hydraulics community, the path for an accurate CHF model remains elusive. This work approaches the prediction of DNB through a physics-informed machine learning-aided framework (PIMLAF) with the objective of achieving superior predictive capabilities for a rod bundle. In view of the limitations in existing macro-scale physics-driven tools, an improved mechanistic model is first proposed, leveraging key concepts in the liquid sublayer dryout and bubble crowding mechanisms. Furthermore, the proposed mechanistic model is able to predict DNB in different heater geometries for a broad range of flow conditions without the need for recalibration. This model is then incorporated as the physics-informed component of the hybrid framework PIMLAF, which takes advantage of established understanding in the field (i.e., domain knowledge [DK]) and uses machine learning (ML) to capture undiscovered information from the mismatch between the actual and DK-predicted output. Two bundle-related case studies using the PWR subchannel and bundle tests (PSBT) database are carried out to illustrate the PIMLAF’s improved performance over traditional approaches for both interpolation and extrapolation purposes. In light of the PIMLAF’s promising potential to reduce prediction error, reactor vendors are encouraged to leverage their in-house experimental efforts and apply the hybrid framework to potentially achieve margin reductions in the minimum DNB ratio (MDNBR) for the designs of interest.

42 ENGINEERING↗
DOE OSTI2024

Thermodynamically consistent Cahn–Hilliard–Navier–Stokes equations using the metriplectic dynamics formalism

Cahn–Hilliard–Navier–Stokes (CHNS) systems describe flows with two-phases, e.g., a liquid with bubbles. Obtaining constitutive relations for general dissipative processes for such systems, which are thermodynamically consistent, can be a challenge. We show how the metriplectic 4-bracket formalism (Morrison and Updike, 2024) achieves this in a straightforward, in fact algorithmic, manner. First, from the noncanonical Hamiltonian formulation for the ideal part of a CHNS system we obtain an appropriate Casimir to serve as the entropy in the metriplectic formalism that describes the dissipation (e.g. viscosity, heat conductivity and diffusion effects). General thermodynamics with the concentration variable and its thermodynamics conjugate, the chemical potential, are included. Having expressions for the Hamiltonian (energy), entropy, and Poisson bracket, we describe a procedure for obtaining a metriplectic 4-bracket that describes thermodynamically consistent dissipative effects. The 4-bracket formalism leads naturally to a general CHNS system that allows for anisotropic surface energy effects. Furthermore, this general CHNS system reduces to cases in the literature, to which we can compare.

Cahn–Hilliard↗
DOE OSTI2021

Initial calculations for source term of Molten Salt Reactors

This paper provides an overview of the current MSR design space and lists unique features of the various designs under consideration. Some general considerations for source terms calculation for Molten Salt Reactors (MSRs) are explained. Applicability and limitations of terminology currently defined for legacy light water reactor (LWR) systems are discussed in the view of MSRs and the need for updated terminology is discussed. Calculations carried out for the Molten Salt Reactor Experiment (MSRE) are discussed with a qualitative comparison to the designs presented. The nature of the fission products (FPs) and actinides for Low enriched uranium, thorium and fast U/Pu fuel cycles employed in representative molten salt reactor systems are discussed. Computational results are obtained from a code (Serpent 2) with online reprocessing. Divergence in source terms when fission product bubbling is demonstrated. The source release for each molten salt reactor during postulated accidents is also presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Effect of feed composition on the production of off-gases during vitrification of simulated low-activity nuclear waste

During the vitrification of nuclear waste, hazardous and radioactive emissions are generated from the feed-to-glass conversion reactions, in addition to discharges from forced air bubbling and air inleakage. Although the major gaseous emissions are water vapor, nitrogen, and carbon dioxide, various monitored environmental pollutants are also released, such as nitrogen oxides or sulfur dioxide. In addition, reactions between organics and nitrates in the feed may also form products of incomplete combustion such as carbon monoxide and acetonitrile. Although off-gas emissions are commonly measured during both laboratory- and pilot-scale melter testing, no predictive tool is currently available to a priori estimate the composition of gaseous emissions during nuclear waste vitrification. This work forms a basis for the development of such predictive tool by measuring gas evolution from a broad range of simulated low-activity waste melter feeds using evolved gas analysis data and developing correlations between the feed and off-gas compositions. Using reaction stoichiometry and regression analysis, we demonstrate that next to the content of nitrogen and organic carbon in the feed, the gaseous emissions are affected by the feed reduction-oxidation conditions – the more the feed is reduced, the less nitrogen monoxide, and more carbon monoxide and acetonitrile evolves. In conclusion, the results presented in this work provide a first step towards reducing the amount of expensive physical melter testing and the regression analysis provides a simple tool for rapid optimization of feed composition with respect to off-gas composition.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Synchrotron X-ray based particle image velocimetry to measure multiphase streamflow and densitometry

Synchrotron-based radiation techniques are increasingly used in the field of flow condition measurements. As a complicated multiphase flow, cavitation has been investigated in a wide variety of industrial fields, to prevent damages on spillways at large dams or intensify the bacteria eradication rate during water treatments. To control the detrimental effects of cavitating flows, the dynamics of and mechanisms affecting cavitation development shall be visualized and identified in real-time. Due to the large reflection and scattering effects of multiphase flow, standard visible-light imaging technologies cannot penetrate vapor bubbles and provide extremely limited information. In this study, the synchrotron source emitted high-flux X-ray pulses are used to capture flow motion and visualize internal structures with sufficiently high spatial and temporal resolutions. The image processing procedures on raw x-ray images are developed using 2D Fourier transform (FT) and wavelet transform (WT) to implement the contrast enhancement and de-noising. Through the edge detection algorithm, tracking either seeded particles or phase interfaces inside the opaque multiphase flow can be achieved to perform the particle image velocimetry. Further, with decomposing particle-only and particle-eliminated subplots from the X-ray phase-contrast image, the time-resolved velocity and void fraction fields are obtained simultaneously, paving a way for further flow condition analysis such as densitometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Temporal characterization of fundamental plasma parameters in pulsed liquid electrode plasma (LEP) optical emission spectrometry

The fundamental characteristics of liquid electrode plasma (LEP), a pulsed plasma source for optical emission spectrometry, were investigated. Two distinct phases were observed during the process of pulsed plasma generation, namely bubble generation and active plasma discharge. Additionally the ionization efficiency of the LEP, with Mg as a representative analyte, was gauged from the ratio of Mg II 279.553 nm to Mg I 285.213 nm emission and was found to increase from about 5% to 20% in a close-to-linear fashion with the discharge voltage from 800 V to 1200 V. The Mg II/Mg I ratio of the LEP was 2.5 to 3 orders of magnitude less than that typically offered by an inductively coupled plasma (ICP) but was comparable to other solution-based glow discharges. It was found that an off-time interval of more than 150 ms between successive discharge pulses was required to obtain a stable pulse-to-pulse discharge current. Temporally resolved emissions of Mg II 279.6 nm, Mg I 285.2 nm, Fe I 373.5 nm, OH band head at 306 nm, and Hα line at 656.3 nm showed that the background species (OH band and Hα line) reached their maximum emission intensities at around 0.5 ms to 0.7 ms with respect to the onset of the discharge pulse whereas the maximum emissions were observed between 0.7 ms to 0.9 ms for analyte species (Mg and Fe lines). The electron density observed in the present work was in the range from 5.7 × 10 15 cm -3 to 8.2 × 10 15 cm -3 , which was similar to those found in an analytical ICP. The temporal averaged OH rotational temperature was 3300 K, which was comparable to the values of an analytical ICP and solution-based glow discharges. By contrast, the temporally averaged Fe I excitation temperature was around 8900 K, which was even higher than that of an analytical ICP and roughly triple the values obtained by techniques based on glow discharge of liquid samples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Carbon dioxide capture with aqueous amino acids: Mechanistic study of amino acid regeneration by guanidine crystallization and process intensification

CO 2 capture from powerplant-generated flue gas via a phase-changing process involving absorption with aqueous amino acids (e.g., glycine or sarcosine) and bicarbonate crystallization with bis-iminoguanidines (e.g., glyoxal-bis-iminoguanidine or GBIG) is investigated in this paper. This process is of high interest due to its potential to decrease the energy penalty for CO 2 capture by significantly reducing the solvent regeneration energy typically associated with aqueous amine solvents. A critical step in the proposed CO 2 capture mechanism is the regeneration of the amino acid by removal of protons and bicarbonate ions from solution through crystallization of GBIGH 2 2+ bicarbonate salt. Here, we investigated the thermodynamics and kinetics of glycine regeneration by crystallization of GBIGH 2 2+ (HCO 3 – ) 2 (H 2 O) 2 . A theoretical model was developed and compared to experimental data to simulate and predict the glycine regeneration and determine its reaction mechanism. This combined experimental and theoretical study led to the conclusion that, while the GBIGH 2 2+ bicarbonate crystallization step provides most of the thermodynamic driving force for the glycine regeneration, the rate-limiting step is the protonation of GBIG prior to crystallization. The CO 2 loading and amino acid regeneration steps were combined into a single, intensified process using a bubble column reactor. The CO 2 loading capacity of GBIG was experimentally determined to be roughly 1.36 mol CO 2 per mol GBIG. These results provide the fundamental basis for developing an effective carbon capture technology with phase-changing amino acid/guanidine absorbents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Pulsed contact glow discharge electrification for online material quantification

Real-time elemental analysis of liquids in extreme environments remains a significant challenge for sensor development, particularly in applications such as advanced nuclear reactors. Plasma-based optical emission spectroscopy offers a promising solution. In this study, we investigate pulsed contact glow discharge electrolysis (CGDE) as an alternative plasma generation method for emission-based material quantification. Using aqueous NaCl solutions, we examine the role of discharge polarity on plasma behavior, emission characteristics, and spectral output. Our results demonstrate that cathodic discharges generate more energetic plasmas, enabling the detection of ionic species such as cerium, while anodic discharges favor atomic emission. These findings reveal a potential dual-mode sensing capability controlled by circuit configuration rather than hardware changes. Through synchronized measurements of emission, current, bubble, and spectral dynamics, we present a model in which electron emission and ohmic heating drive plasma formation. This model enhances our understanding of pulsed-CGDE and supports its potential as a versatile, adaptable platform for future deployment in high-temperature liquids such as molten salts.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
DOE OSTI2023

Dynamics of oil–water interface at the beginning of the ultrasonic emulsification process

A lot of effort has been dedicated in recent years towards understanding the basics of cavitation induced emulsification, mainly in the form of single cavitation bubbles. Regarding bulk acoustic emulsification, a lot less research has been done. In our here presented work we utilize advanced high-speed observation techniques in visible light and X-Rays to build upon that knowledge and advance the understanding of bulk emulsion preparation. During research we discovered that emulsion formation has an acute impact on the behavior of the interface and more importantly on its position relative to the horn, hence their interdependence must be carefully studied. We did this by observing bulk emulsification with 2 cameras simultaneously and corroborating these measurements with observation under X-Rays. Since the ultrasonic horns location also influences interface behavior, we shifted its initial position to different locations nearer to and further away from the oil–water interface in both phases. We found that a few millimeters distance between the horn and interface is not enough for fine emulsion formation, but that they must be completely adjacent to each other, with the horn being located inside the oil–water interface. We also observed some previously undiscovered phenomena, such as the splitting of the interface to preserve continuous emulsion formation, climbing of the interface up the horn and circular interface protrusions towards the horn forming vertical emulsion streams. Interestingly, no visible W/O emulsion was ever formed during our experiments, only O/W regardless of initial horn position.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Characterization of the acoustic cavitation in ionic liquids in a horn-type ultrasound reactor

Most ultrasound-based processes root in empirical approaches. Because nearly all advances have been conducted in aqueous systems, there exists a paucity of information on sonoprocessing in other solvents, particularly ionic liquids (ILs). In this work, we modelled an ultrasonic horn-type sonoreactor and investigated the effects of ultrasound power, sonotrode immersion depth, and solvent’s thermodynamic properties on acoustic cavitation in nine imidazolium-based and three pyrrolidinium-based ILs. The model accounts for bubbles, acoustic impedance mismatch at interfaces, and treats the ILs as incompressible, Newtonian, and saturated with argon. Following a statistical analysis of the simulation results, we determined that viscosity and ultrasound input power are the most significant variables affecting the intensity of the acoustic pressure field (P), the volume of cavitation zones (V), and the magnitude of the maximum acoustic streaming surface velocity (u). V and u increase with the increase of ultrasound input power and the decrease in viscosity, whereas the magnitude of negative P decreases as ultrasound power and viscosity increase. Probe immersion depth positively correlates with V, but its impact on P and u is insignificant. 1-alkyl-3-methylimidazolium-based ILs yielded the largest V and the fastest acoustic jets – 0.77cm 3 and 24.4ms -1 for 1-ethyl-3-methylimidazolium chloride at 60W. 1-methyl-3-(3-sulfopropyl)-imidazolium-based ILs generated the smallest V and lowest u – 0.17cm 3 and 1.7ms -1 for 1-methyl-3-(3-sulfopropyl)-imidazolium p-toluene sulfonate at 20W. Sonochemiluminescence experiments validated the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydrodynamic cavitation

The paper investigates the oil–water emulsification process inside a micro-venturi channel. More specifically, the possible influence of Kelvin-Helmholtz instability on the emulsification process. High-speed visualizations were conducted inside a square venturi constriction with throat dimensions of 450 µm by 450 µm, both under visible light and X-Rays. We show that cavity shedding caused by the instability results in the formation of several cavity vortices. Their rotation causes the deformation of the oil stream into a distinct wave-like shape, combined with fragmentation into larger drops due to cavitation bubble collapse. Later on, the cavity collapse further disperses the larger drops into a finer emulsion. Thus, it turns out that the Kelvin-Helmholtz instability is similarly characteristic for hydrodynamic cavitation emulsification inside a microchannel as is the Rayleigh-Taylor instability for acoustically driven emulsion formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

A vorticity-based criterion to characterise leading edge dynamic stall onset

We propose a more conservative, physically-intuitive criterion, namely, the boundary enstrophy flux ( BEF), to characterise leading-edge-type dynamic stall onset in incompressible flows. Our results are based on wall-resolved large-eddy simulations of pitching aerofoils, with fine spatial and temporal resolution around stall onset. We observe that |BEF| reaches a maximum within the stall onset regime identified. By decomposing the contribution to BEF from the flow field, we find that the dominant contribution arises from the laminar leading edge region, due to the combined effect of large clockwise vorticity and favourable pressure gradient. A relatively small contribution originates from the transitional/turbulent laminar separation bubble (LSB) region, due to LSB-induced counter-clockwise vorticity and adverse pressure gradient. This results in BEF being nearly independent of the integration length as long as the region very close to the leading edge is included. This characteristic of BEF yields a major advantage in that the effect of partial or complete inclusion of the noisy LSB region can be filtered out, without changing the BEF peak location in time significantly. Next, we analytically relate BEF to the net wall shear and show that its critical value (=max(|BEF|) ) corresponds to the instant of maximum net shear prevailing at the wall. Finally, we have also compared BEF with the leading edge suction parameter ( LESP) and find that the former reaches its maximum value between 0.3° and 0.8° of rotation earlier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Scale interactions and anisotropy in Rayleigh–Taylor turbulence

We study energy scale transfer in Rayleigh–Taylor (RT) flows by coarse graining in physical space without Fourier transforms, allowing scale analysis along the vertical direction. Two processes are responsible for kinetic energy flux across scales: baropycnal work Λ, due to large-scale pressure gradients acting on small scales of density and velocity; and deformation work Π, due to multiscale velocity. Our coarse-graining analysis shows how these fluxes exhibit self-similar evolution that is quadratic-in-time, similar to the RT mixing layer. We find that Λ is a conduit for potential energy, transferring energy non-locally from the largest scales to smaller scales in the inertial range where Π takes over. In three dimensions, Π continues a persistent cascade to smaller scales, whereas in two dimensions Π rechannels the energy back to larger scales despite the lack of vorticity conservation in two-dimensional (2-D) variable density flows. This gives rise to a positive feedback loop in 2-D RT (absent in three dimensions) in which mixing layer growth and the associated potential energy release are enhanced relative to 3-D RT, explaining the oft-observed larger α values in 2-D simulations. Despite higher bulk kinetic energy levels in two dimensions, small inertial scales are weaker than in three dimensions. Moreover, the net upscale cascade in two dimensions tends to isotropize the large-scale flow, in stark contrast to three dimensions. Furthermore, our findings indicate the absence of net upscale energy transfer in three-dimensional RT as is often claimed; growth of large-scale bubbles and spikes is not due to ‘mergers’ but solely due to baropycnal work Λ.

42 ENGINEERING↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Speed of Sound Measurements of Binary Mixtures of trans -1,2-Difluoroethylene (R-1132(E)) with Difluoromethane (R-32) or 2,3,3,3-Tetrafluoropropene (R-1234yf)

The speed of sound of two binary mixtures containing 0.336 mole fraction trans-1,2- difluoroethylene (R-1132(E)) with difluoromethane (R-32) and 0.435 mole fraction R-1132(E) with 2,3,3,3-tetrafluoropropene (R-1234yf) were measured with a dual-path pulse-echo instrument. The speed of sound was measured along several pseudo-isochores for each blend at temperatures ranging from 230 K to 325 K for blends with R-32 and to 342 K for blends with R- 1234yf. Measurements started at pressures just above each mixtures bubble point pressure and were limited to 8 MPa to avoid potential disproportionation reactions of the R-1132(E). The data were compared to multi-fluid models incorporating Helmholtz-energy-explicit equations of state (EOS) for each pure fluid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Theoretical Investigations of Anharmonic Effects and Phonon Transport in the Cubic Phase of Crystalline Perovskite CsPbCl 3

The role of anharmonic effects on lattice dynamics and thermal transport has been investigated in the cubic phase of the CsPbCl 3 perovskite. Limitations of the harmonic approximation, which lead to phonon instabilities in the Brillouin zone, were addressed based on self-consistent phonon theory in combination with first-principles calculations and by incorporation of frequency renormalization effects from bubble and loop diagrams. This theoretical approach demonstrates significant improvements compared to results obtained using a self-consistent phonon dynamic matrix. Both the cubic-to-tetragonal phase transition temperature as well as the predicted lattice thermal conductivity values show good agreement with selected sets of available experimental data. In conclusion, the accuracy of the predicted thermal conductivity data is also tested against systems significantly larger than those allowed by quantum calculations, by using molecular dynamic simulations with machine-learning force fields.

Lattices↗