Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase drift”

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

CTF Improved Drag Model and Flow Regime Transition Criteria

The demand for accurate prediction of two-phase flow behavior in a boiling water reactor (BWR) requires a comprehensive understanding of flow regime, void fraction, heat transfer, and pressure drop. The CTF subchannel code, which is used for the Thermal/Hydraulic (T/H) solution in the Consortium for Advanced Simulation of Light Water Reactors (CASL)-developed Virtual Environment for Reactor Application (VERA) core simulator, is being further developed for BWR applications. In support of this goal, the present work highlights some of the two-phase closure model developments towards improving the CTF void fraction prediction, especially for subcooled boiling. The drift-flux approach has been well-developed for upward dispersed two-phase flows and proven to be accurate in predicting void fraction in bubbly and slug flow regimes. In this work, these kinematic constitutive relations for the drift-flux velocity have been implemented into CTF to describe the interfacial drag of bubbly flow as an alternative to the existing model for better void fraction prediction. The success of these constitutive relations also relies on a good flow regime map that accounts for flow conditions and channel geometry. A more reliable flow regime transition criteria that account for the flow condition has also been implemented in this study for modeling the flow regime transition criteria. The newly implemented models are shown to give improved void fraction predictions in comparison to experimental data.

Hizoum, Belgacem↗

Elucidating the reversible exsolution–dissolution behaviour of high-entropy oxides in crystalline and amorphous phases

High-entropy oxides (HEOs), as a subclass of high-entropy materials (HEMs), offer a versatile platform for catalysis by leveraging entropy-stabilized solid solutions with tunable compositions, lattice structures, and electronic properties. While exsolution–dissolution of metal species in crystalline HEOs has emerged as a promising strategy for reversible active sites regeneration, the dynamic behaviour of HEOs possessing amorphous nature remains under-explored, particularly the difference with crystalline counterparts. In this work, we systematically investigate the architecture-dependent exsolution–dissolution behavior of HEOs by comparing a crystalline-phase HEO (c-HEO) and an amorphous-phase HEO (a-HEO), both comprising Ni, Mg, Cu, Zn, and Co as principal metal elements. Using a combination of in situ variable-temperature X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), electron microscopy, and in situ CO diffuse reflectance infrared Fourier transform spectroscopy (CO-DRIFTS), the structural evolution of the two HEO phases under redox conditions was elucidated. Both materials exhibit reversible exsolution of metallic species or alloys in reducing environments, followed by re-incorporation into the host lattice upon oxidation. Remarkably, the a-HEO demonstrates more facile and dynamic self-healing behavior, with alloy exsolution and dissolution occurring under milder conditions because of its enhanced reducibility and structural disorder. This study provides critical insights into the design of next-generation regenerable catalysts based on amorphous HEOs, highlighting the role of phase structure in governing reversible metal-site formation dynamics and catalytic performance.

Wang, Qingju [Univ. of Tennessee, Knoxville, TN (U↗

Implementation of a Drift Flux Model into SAM with Development of a Verification and Validation Test Suite for Modeling of Noncondensable Gas Mixtures

The advanced thermal-hydraulic system code, System Analysis Module (SAM), was originally developed for the modeling of single-phase flow in advanced reactors. It has since been expanded to include a four-equation drift flux model for the modeling of two-phase flows containing a noncondensable gas. The model was expanded to support the modeling of molten salt reactor (MSR) designs in which the fuel is directly dissolved in the circulating coolant. These designs have shown that circulating gas bubbles can play an important role in the management of fission products and the operational behavior of the reactor. A drift flux model was implemented to more accurately capture the localized behavior of the void in the core and its impact on the mass transfer of fission products. A thorough assessment of the new model was performed by developing a verification and validation test suite. Verification problems were designed to test all major terms in the new governing equations. The new model converged to the correct solution at the expected order of accuracy for all verification cases. The validation cases included a wide range of flow and void conditions in different pipe geometries. Although higher void experiments show a slight underprediction of void by the drift flux model, experiments that aim to reproduce Molten Salt Reactor Experiment (MSRE) experimental conditions show good agreement with the model. The gas transport model was activated for a SAM model of the MSRE to demonstrate that it can be used in a more complex model. Finally, this gas transport model will be used along with an interfacial area transport equation being implemented in SAM for the prediction of mass transport behavior in MSR conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Volume IV. The DUNE far detector single-phase technology

The preponderance of matter over antimatter in the early universe, the dynamics of the supernovae that produced the heavy elements necessary for life, and whether protons eventually decay -- these mysteries at the forefront of particle physics and astrophysics are key to understanding the early evolution of our universe, its current state, and its eventual fate. DUNE is an international world-class experiment dedicated to addressing these questions as it searches for leptonic charge-parity symmetry violation, stands ready to capture supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. Central to achieving DUNE's physics program is a far detector that combines the many tens-of-kiloton fiducial mass necessary for rare event searches with sub-centimeter spatial resolution in its ability to image those events, allowing identification of the physics signatures among the numerous backgrounds. In the single-phase liquid argon time-projection chamber (LArTPC) technology, ionization charges drift horizontally in the liquid argon under the influence of an electric field towards a vertical anode, where they are read out with fine granularity. A photon detection system supplements the TPC, directly enhancing physics capabilities for all three DUNE physics drivers and opening up prospects for further physics explorations. The DUNE far detector technical design report (TDR) describes the DUNE physics program and the technical designs of the single- and dual-phase DUNE liquid argon TPC far detector modules. Volume IV presents an overview of the basic operating principles of a single-phase LArTPC, followed by a description of the DUNE implementation. Each of the subsystems is described in detail, connecting the high-level design requirements and decisions to the overriding physics goals of DUNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Phase-space Energization of Ions in Oblique Shocks

Abstract Examining energization of kinetic plasmas in phase space is a growing topic of interest, owing to the wealth of data in phase space compared to traditional bulk energization diagnostics. Via the field-particle correlation (FPC) technique and using multiple means of numerically integrating the plasma kinetic equation, we have studied the energization of ions in phase space within oblique collisionless shocks. The perspective afforded to us with this analysis in phase space allows us to characterize distinct populations of energized ions. In particular, we focus on ions that reflect multiple times off the shock front through shock-drift acceleration, and how to distinguish these different reflected populations in phase space using the FPC technique. We further extend our analysis to simulations of three-dimensional shocks undergoing more complicated dynamics, such as shock ripple, to demonstrate the ability to recover the phase-space signatures of this energization process in a more general system. This work thus extends previous applications of the FPC technique to more realistic collisionless shock environments, providing stronger evidence of the technique’s utility for simulation, laboratory, and spacecraft analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Cooling with Thick Wedges for a Muon Collider

A muon collider is a particle accelerator that collides muons rather than protons or electrons. An important component of such a collider is cooling, in which the emittance of the muon beam is reduced to make it suitable for use in a collider. In the last stage of this process, called final 4D cooling, emittance in the transverse axes is reduced while emittance in the longitudinal axis is allowed to grow. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on single thick wedges. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We characterized the effects of various parameters on the performance of this system, including the initial emittances and Twiss parameters, the length and angle of the wedges, the length of the drift channel, and the length, frequency, gradient, and phase of the RF cavities. We optimized the design parameters of the cooling channel and produced two conceptual designs (corresponding to two possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ≈ 3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ultralow-temperature cryogenic transmission electron microscopy using a new helium flow cryostat stage

Advances in cryogenic electron microscopy have opened new avenues for probing quantum phenomena in correlated materials. This study reports the installation and performance of a new side-entry condenZero cryogenic cooling system for JEOL (Scanning) Transmission Electron Microscopes (S/TEM), utilizing compressed liquid helium (LHe) and designed for imaging and spectroscopy at ultra-low temperatures. The system includes an external dewar mounted on a vibration-damping stage and a pressurized, low-noise helium transfer line with a remotely controllable needle valve, ensuring stable and efficient LHe flow with minimal thermal and mechanical noise. Performance evaluation demonstrates a stable base temperature of 4.37 K measured using a Cernox bare chip sensor on the holder with temperature fluctuations within ±0.004 K. Complementary in-situ electron energy-loss spectroscopy (EELS) via aluminum bulk plasmon analysis was used to measure the local specimen temperature and validate cryogenic operation during experiments. The integration of cryogenic cooling with other microscopy techniques, including electron diffraction and Lorentz TEM, was demonstrated by resolving charge density wave (CDW) transitions in NbSe2 using electron diffraction, and imaging nanometric magnetic skyrmions in MnSi via Lorentz TEM. In conclusion, this platform provides reliable cryogenic operation below 7 K, establishing a low-drift route for direct visualization of electronic and magnetic phase transformations in quantum materials.

Charge density wave↗

Stochastic fluctuation and transport of tokamak edge plasmas with the resonant magnetic perturbation field

Here we present that a statistical method known as the complexity–entropy analysis is useful to characterize a state of plasma turbulence and flux in the resonant magnetic perturbation (RMP) edge localized mode (ELM) control experiment. The stochastic pedestal top temperature fluctuation in the RMP ELM suppression phase is distinguished from the chaotic fluctuation in the natural ELM-free phase. It is discussed that the stochastic temperature fluctuation can be originated from the narrow layer of the field penetration on the pedestal top. The forced magnetic island can emit the resonant drift wave of comparable sizes (relatively low-k) in the RMP ELM suppression phase, and it can result in the generation of stochastic higher wavenumber fluctuations coupled to tangled fields around the island. The analysis of the ion saturation current measurement around the major outer striking point on the divertor shows that it also becomes more stochastic as the stronger plasma response to the RMP field is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

High-fidelity kinetic modeling of instabilities and gyromotion physics in nonuniform low-beta plasmas

A fourth-order accurate continuum kinetic Vlasov solver and a systematic method for constructing customizable kinetic equilibria are demonstrated to be powerful tools for the study of nonuniform collisionless low-beta plasmas. The noise-free methodology is applied to investigate two gradient-driven instabilities in 4D (x,y,vx,vy) phase space: the Kelvin–Helmholtz instability and the lower hybrid drift instability. Nonuniform two-species configurations where ion gyroradii are comparable to gradient scale lengths are explored. The approach sheds light on the evolution of the pressure tensor in Kelvin–Helmholtz instabilities and demonstrates that the associated stress tensor deviates significantly from the gyroviscous stress tensor. Even at high magnetization, first-order approximations to finite-gyromotion physics are shown to be inadequate for the Kelvin–Helmholtz instability, as shear scales evolve to become on par with gyromotion scales. The methodology facilitates exploring transport and energy partitioning properties associated with lower hybrid drift instabilities in low-beta plasma configurations. Distribution function features are captured in detail, including the formation of local extrema in the vicinity of particle-wave resonances. Finally, the approach enables detailed targeted investigations and advances kinetic simulation capability for plasmas in which gyromotion plays an important role.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Implementation of a Drift-Flux Model in SAM for Modeling of Passively Transported GAS in Molten Salt Reactors

The Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is further developing the Multiphysics Object Oriented Simulation Environment (MOOSE) system thermal hydraulics (T/H) code, System Analysis Module (SAM), to include modeling of molten salt reactor (MSR) designs. MSR designs in which the fuel is dissolved in the salt coolant poses unique modeling and simulation (M&S) challenges because fission products—some of which will be noncondensable gasses—travel with the coolant, thus allowing for deposition and interaction throughout the entire system flow loop. This behavior must be modeled to successfully understand MSR system behavior during normal and off-normal operation events. In past work, SAM has been extended to include a passive species transport capability for modeling of MSR fission products that travel with the salt coolant. This work extends the model to include a drift-flux model to capture the tendency of the gas phase to move at a different velocity than the bulk liquid. Additional constitutive models have also been added for calculation of gas volume fraction, bubble radius, and interfacial area, which will be important for calculation of species phase migration in future coupling activities. Testing has been performed to ensure that the new model agrees with expected values for several simple models, but additional future work will be required to test the model for more complex geometries that will be needed for full-scale MSR simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

ANTS

The ANTS code (Alternate Non-Linear Two-phase Solver) is based on a novel non-linear solution algorithm for the solution of the two-phase, subchannel fluid equations. It achieves its performance through decoupling of the two-phase momentum equations (axial and transverse) from the axial phasic mass and energy equations which allows for a nested non-linear iteration scheme. This enables a plane-by-plane solution which the inner iteration focuses on a reduced non-linear equation set for the primitives in phasic mass flow rate, enthalpy and void for each node edge. Single node edges are coupled as part of the outer iteration via surface mass fluxes which appear as source terms in the inner iteration scheme. The outer iteration readily accommodates two-phase flow phenomena closure relationships for subchannel mixing and void drift. A primary feature is the use of a non-staggered mesh computational mesh and steady-state iterative solver in contrast to all existing subchannel codes.

Kropaczek, David J↗

Understanding Selectivity in CO2 Hydrogenation to Methanol for MoP Nanoparticle Catalysts Using In Situ Techniques

Molybdenum phosphide (MoP) catalyzes the hydrogenation of CO, CO2, and their mixtures to methanol, and it is investigated as a high-activity catalyst that overcomes deactivation issues (e.g., formate poisoning) faced by conventional transition metal catalysts. MoP as a new catalyst for hydrogenating CO2 to methanol is particularly appealing for the use of CO2 as chemical feedstock. Herein, we use a colloidal synthesis technique that connects the presence of MoP to the formation of methanol from CO2, regardless of the support being used. By conducting a systematic support study, we see that zirconia (ZrO2) has the striking ability to shift the selectivity towards methanol by increasing the rate of methanol conversion by two orders of magnitude compared to other supports, at a CO2 conversion of 1.4% and methanol selectivity of 55.4%. In situ X-ray Absorption Spectroscopy (XAS) and in situ X-ray Diffraction (XRD) indicate that under reaction conditions the catalyst is pure MoP in a partially crystalline phase. Results from Diffuse Reflectance Infrared Fourier Transform Spectroscopy coupled with Temperature Programmed Surface Reaction (DRIFTS-TPSR) point towards a highly reactive monodentate formate intermediate stabilized by the strong interaction of MoP and ZrO2. This study definitively shows that the presence of a MoP phase leads to methanol formation from CO2, regardless of support and that the formate intermediate on MoP governs methanol formation rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and Validation of Two-Phase Flow Models in MOOSE and Application to Molten Salt Reactors

Two-phase flow in Molten Salt Reactors (MSRs) is important as it impacts reactivity evolution, reactor transient response, and the removal of species dissolved in the molten salt through gas phase transfer. Therefore, accurately predicting the gas distribution and the associated liquid-gas interface area in MSRs is essential for their design and operation. Recently, we integrated a new two-phase model into Idaho National Laboratory (INL)’s Multiphysics Object-Oriented Simulation Environment (MOOSE): a multi-D generalization of a mixture drift-flux model. It provides greater computational efficiency, which is typically preferred for modeling reactor transients. However, the mixture model’s accuracy in capturing void distribution and interfacial area in MSRs still needs to be assessed. This article begins with a description of the mathematical framework for the two-phase model implemented in MOOSE. It then presents validation of these models against relevant experimental data. Finally, the model is applied to the Molten Salt Reactor Experiment case study, analyzing various operational conditions such as different rates of fission product volatilization and diverse cover gas entrainment scenarios at the reactor pump. The article concludes by assessing the suitability of the mixture drift-flux model for capturing the two-phase flow dynamics critical to MSR operations

42 - ENGINEERING↗

Observations and Simulations of Dropout Events and Flux Decays in October 2013: Comparing MEO Equatorial With LEO Polar Orbit

We compare ESA PROBA-V observations of electron flux at LEO with those from the NASA Van Allen Probes mostly at MEO for October 2013. Dropouts are visible at all energy during four storms from both satellites. Equatorially trapped electron fluxes are higher than at LEO by 10 2 (<1 MeV) to 10 5 (>2.5 MeV). We observe a quite isotropic structure of the outer belt during quiet times, contrary to the inner belt, and pitch angle dependence of high energy injection. We find a very good overlap of the outer belt at MEO and LEO at ~0.5 MeV. We use test-particle simulations of the energetic electrons trapped in the terrestrial magnetic field to study the outer radiation belt electron flux changes during geomagnetic storms. We show that the Dst (Disturbance storm time) effect during the main phase of a geomagnetic storm results in a betatron mechanism causing outward radial drift and a deceleration of the electrons. This outward drift motion is energy independent, pitch angle-dependent, and represents a significant distance (~1 L-shell at L = 5 for moderate storms). At fixed L-shell, this causes a decay of the LEO precipitating flux (adiabatic outward motion), followed by a return to the normal state (adiabatic inward motion) during main and recovery phases. Dst effect, associated with magnetopause shadowing and radial diffusion can explain the main characteristics of outer radiation belt electron dropouts in October 2013. We also use Fokker-Planck simulations with event-driven diffusion coefficients at high temporal resolution, to distinguish instantaneous loss from the gradual scattering that depopulates the slot region and the outer belt after storms. Simulations reproduce the slot formation and the gradual loss in the outer belt. The typical energy dependence of these losses leads to the absence of scattering for relativistic and ultra-relativistic electrons in the outer belt, oppositely to dropouts.

79 ASTRONOMY AND ASTROPHYSICS↗

Velocity-space signatures of shock-drift acceleration at quasi-perpendicular collisionless shocks

The shock-drift acceleration of ions at quasi-perpendicular shocks is a well-known kinetic mechanism for the acceleration of a small fraction of incoming ions to high energy. Here, we use a suite of sixteen hybrid simulations of quasi-perpendicular collisionless shocks over the range of Alfvén Mach number 4.3 ≤ M A ≤ 15.8 (corresponding to a range of fast magnetosonic Mach numbers 2.6 ≤ M f ≤ 9.4) and shock-normal angle 45°≤ θ Bn ≤ 90° to identify the velocity-space signature of shock-drift acceleration using the field-particle correlation technique. We show that the features of the ion velocity distribution in the shock foot and ramp regions can be clearly interpreted by analysis of the single-particle trajectory of a reflected ion through the full 3D-3V phase space. The characteristic features of the velocity-space signature of shock-drift acceleration remain qualitatively robust over the full parameter range of our simulations, providing a potential means for its identification using single-point spacecraft measurements. At higher Alfvén Mach numbers M A ≳8 (M f ≳5), kinetic instabilities generate fluctuations of the electromagnetic fields within the shock transition region, leading to fluctuations in and smearing out of the resulting velocity-space signatures, but the signature remains generally robust and identifiable. The results on the shock-drift acceleration of ions presented here represent a novel means to determine more completely the partitioning of upstream bulk flow kinetic energy into plasma heating, particle acceleration, and electromagnetic fields in collisionless shocks.

Howes, Gregory G. [Univ. of Iowa, Iowa City, IA (U↗

Pro‐ L * ‐ A Probabilistic L * Mapping Tool for Ground Observations

Abstract Both ground and space observations are used extensively in the modeling of space weather processes within the Earth’s magnetosphere. In radiation belt physics modeling, one of the key phase‐space coordinates is L *, which indicates the location of the drift paths of energetic electrons. Global magnetic field models allow a subset of locations on the ground (mainly subauroral) to be mapped along field lines to a location in space and transformed into L *, provided that the initial ground location maps to a closed drift path. This allows observations from ground, or low‐altitude space‐based platforms to be mapped into space in order to inform radiation belt modeling. Many data‐based magnetic field models exist; however, these models can significantly disagree on mapped L * values for a single point on the ground, during both quiet times and storms. We present a state of the art probabilistic L * mapping tool, Pro‐ L *, which produces probability distributions for L * corresponding to a given ground location. Pro‐ L * has been calculated for a high resolution magnetic latitude by magnetic local time grid in the Earth’s Northern Hemisphere. We have developed the probabilistic model using 11 years of L * calculations for seven widely used magnetic field models. Usage of the tool is highlighted for both event studies and statistical models, and we demonstrate a number of potential applications.

79 ASTRONOMY AND ASTROPHYSICS↗