SEARCH · Search NASA

Results for “normalizing flows”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Thermal oxidation of nuclear graphite and pyrolytic carbon coatings

The oxidation of pyrolytic carbon (PyC) deposited via fluidized bed chemical vapor deposition was characterized and compared with that of standard nuclear-grade graphite. The materials were heated at 700 to 1000 °C in a thermogravimetric analysis system under 20% v/v O 2 flow, allowing for direct comparison of dynamic oxidative mass change in each material. Further, three different PyC samples fabricated under different conditions exhibited variation in total mass loss and mass loss rate, varying by as much as 709 mg/cm 2 in total mass loss and 14.2 (mg/cm 2 )/min in mass loss rate at a single temperature. These variations highlight the correlation between PyC microstructure/defect density and oxidation susceptibility. Additionally, changes in the microstructure and composition between PyC and graphite were characterized via scanning electron microscopy and correlated to the mass loss results. The results of this work have implications toward the safety of tristructural isotropic (TRISO) and other coated particle fuels, especially under off-normal conditions, given the limited information that exists about the oxidation behavior of PyC.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multi-frame x-ray radiography and image tracking for quantification of expansion in laser-driven tin ejecta microjets

One regime of experimental particle-laden flow study involves ejecta microjets—often defined as a stream of micrometer-scale particles generated through shock interaction with a non-uniform surface and generally travel above 1 km/s. In order to capture the change in characteristics as a function of propagation time, we apply a multi-frame x-ray radiography platform to observe and track the jet transport dynamics. A synchrotron x-ray source allows us to perform quantitative analyses and comparisons between the eight images captured by the imaging system. Observation of a single jet through time allows the use of a cross correlation algorithm to independently track various regions within the jet and quantify the jet expansion over time using normalized area and normalized areal density values. Through a comparison with the calculated values of ballistic transport, these findings show less expansion than expected for ballistically transporting particles. In conclusion, this work combines multi-frame synchrotron radiography with image tracking to establish a foundation for future studies on jet transport and particle interaction dynamics.

Sun, Yuchen [Lawrence Livermore National Laborator↗

DOE OSTI2024

Macroscopic modeling of gas permeability in hierarchical micro/nanoporous media: A unified characterization of rarefaction using Klinkenberg theory and equivalent diameter

Estimating gas transport through a hierarchical micro/nanoporous system is challenging due to non-equilibrium gas dynamics. The primary difficulty lies in determining the rarefaction level, because identifying a representative flow dimension in a complex porous system with multiple pore scales is not straightforward. Our study performed a pore-level analysis for gas permeability in dual-scale porous media with varying porosity, throat size, and secondary pore size under different rarefaction conditions. We found that secondary porosity negatively affects permeability due to increased friction forces, with this influence growing as the secondary pore size and porosity increase until the secondary pore becomes comparable to the throat. However, rarefaction reduces the effects of secondary pores due to boundary slip. Traditional Knudsen number (Kn) calculations based on Darcy-defined height failed to accurately describe the rarefaction effects on gas permeability. Instead, we introduced an equivalent diameter to calculate the Kn, which provided an accurate normalization of apparent gas permeability independent of pore geometry. Furthermore, the extended Kozeny–Carman–Klinkenberg model developed in our previous study successfully yielded a macroscopic model for apparent gas permeability in hierarchical micro/nanoporous systems as a function of the traditional Darcy height and porosity.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2026

Thermal-Hydraulic-Mechanical Modeling and Simulation of Sodium-Potassium–Cooled MARVEL Microreactor Core

The U.S. Department of Energy's Microreactor Program, with Idaho National Laboratory's development of a nuclear microreactor applications test bed named MARVEL, aims to support R&D for the deployment of small, transportable reactors across civilian, industrial, and defense sectors. The MARVEL microreactor, an 85-kWth thermal fission reactor, incorporates TRIGA nuclear fuel and a sodium-potassium eutectic as its primary coolant, designed for safety and efficiency, with natural circulation eliminating the risk of critical heat flux conditions. The reliance on natural circulation for primary cooling means the reactor avoids using fuel spacers to minimize core pressure drop, which could disrupt the primary coolant's natural flow. However, the reactor core’s tight P/D ratio of 1.05, in the absence of fuel spacers, could pose a risk of fuel rod contact and increased peak cladding temperatures. To ensure the reactor safety, this study conducted computational modeling and simulations to investigate the reactor's thermal-hydraulic-mechanical characteristics, including the reactor core heat transfer coefficients, the potential for rod-to-rod contact, and assessed its impact on peak cladding temperature and overall reactor safety. The computational analyses of the MARVEL microreactor core revealed that the thermal deformation of fuel rods under worst-case scenario may lead to the fuel rod contact, but the peak cladding temperatures will remain significantly lower than the safety criteria, ensuring the safety operation reactor without fuel spacers under normal operating conditions.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2026

Monitored Fluctuating Hydrodynamics

We introduce a hydrodynamic framework for describing monitored classical stochastic processes. We study the conditional ensembles for these monitored processes—i.e., we compute spacetime correlation functions conditioned on a fixed, typical measurement record. In the presence of global symmetries we show that these conditional ensembles can undergo measurement-induced “sharpening” phase transitions as a function of the monitoring rate; moreover, even weak monitoring can give rise to novel critical phases, derived entirely from a classical perspective. We give a simple hydrodynamic derivation of the known “charge-fuzzy phase” for weakly monitored diffusive many-body quantum systems. We show that although the unmonitored symmetric and asymmetric exclusion processes are in different universality classes of transport, the fluctuations in their conditional ensembles flow to the same fixed point with emergent relativistic invariance under monitoring. On the other hand, weakly monitored systems with non-Abelian symmetries enter a novel strongly coupled fixed point with nontrivial dynamical exponent, which we characterize. Our formalism naturally accounts for monitoring general observables, such as currents or density gradients, and allows for a direct calculation of information-theoretic diagnostics of sharpening transitions, including the Shannon entropy of the measurement record.

classical statistical mechanics↗

DOE OSTI2024

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗

DOE OSTI2024

Alternating conduction and convection drying of paper – an experimental analysis with a continuous data acquisition approach

In conventional multi-cylinder drying of paper and board, both conductive drying from steam-heated dryer cylinders and convective drying by flowing air over the paper surface in the pockets are used. Conductive drying from steam-heated drying cylinders is a critical component in providing the necessary thermal energy to paper and board as they dry. Steam temperature and internal and external resistances at the contacting surface are critical process parameters influencing the conductive drying process. An experimental setup was developed to study the alternating conductive and convective drying of paper and board. Paper sheet moisture, temperature, and temperature distribution within the heated platen and the instantaneous heat flux as the sheet was being dried were measured. The instantaneous heat flux, contact heat transfer coefficient, and drying rates were determined as drying proceeds. Experimental results, as well as comparisons to literature and commercial data, are presented. The conductive heat transfer coefficients determined were compared to traditional correlations normally used in the modeling of paper drying. Similarly, the convective heat and mass transfer coefficients are also determined and compared to literature data. In addition to the evaluation of alternating conductive and convective drying characteristics of paper and board, the potential inclusion of auxiliary energy components will also be included. Experimental results from the conduction and convection drying system are presented. Furthermore, this data will be useful in process development, intensification of manufacturing processes, and modeling and simulation of paper drying processes.

42 ENGINEERING↗

DOE OSTI1995

Nuclear Safety [Vol. 36, No. 1, January-June 1995]

Nuclear Safety is a journal that covers significant issues in the field of nuclear safety. Its primary scope is safety in the design, construction, operation, and decommissioning of nuclear power reactors worldwide and the research and analysis activities that promote this goal, but it also encompasses the safety aspects of the entire nuclear fuel cycle, including fuel fabrication, spent-fuel processing and handling, and nuclear waste disposal, the handling of fissionable materials and radioisotopes, and the environmental effects of all these activities. Table of Contents for this issue follows. THE CHORNOBYL ACCIDENT: 1 The Chornobyl Accident Revisited, Part II: The State of the Nuclear Fuel Located Within the Chornobyl Sarcophagus, A A. Borovoi and A. R. Sich; GENERAL SAFETY CONSIDERATIONS: 33 Nuclear Power Safety in Central and Eastern Europe, R. Wilson; 46 Safety of Nuclear Power Reactors in the Former Eastern European Countries, S. Chakraborty; 53 Technical Note: On the Definition of Common-Cause Failures, H. Paula; ACCIDENT ANALYSIS: 58 Modeling and Analysis of Core-Debris Recriticality During Hypothetical Severe Accidents in the Advanced Neutron Source Reactor, S.-H. Kim, V. Georgevich, D. B. Simpson, C. O. Slater, and R. P. Taleyarkhan; 68 Ignitability of Hydrogen/Oxygen/Diluent Mixtures in the Presence of Hot Surfaces, R. K. Kumar and G. W. Koroll; 94 Coupled RELAP5 and CONTAIN Accident Analysis Using PVM, K. A. Smith, A. J. Baratta, and G. E. Robinson; CONTROL AND INSTRUMENTATION: 109 Application of Fuzzy Logic in Nuclear Reactor Control Part I: An Assessment of State-of-the-Art, A. S. Heger, N. K. Alang-Rashid, and M. Jamshidi; DESIGN FEATURES: 122 Twenty-Third DOE/NRC Nuclear Air-Cleaning and Treatment Conference, R. R. Bellamy, J. J. Hayes, and M. W. First; ENVIRONMENTAL EFFECTS: 135 Atmospheric Dispersion and the Radiological Consequences of Normal Airborne Effluents from a Nuclear Power Plant, D. Fang, C. Z. Sun, and L. Yang; 142 Calculation of Distribution Coefficients for Radionuclides in Soils and Sediments, I. Puigdomenech and U. Bergstrom: OPERATING EXPERIENCES: 155 Reactor Shutdown Experience, Compiled by J. W. Cletcher; U.S. NUCLEAR REGULATORY COMMISSION INFORMATION AND ANALYSES: 158 Operating Experience Feedback Report—Reliability of Safety-Related Steam Turbine-Driven Standby Pumps Used in U.S. Commercial Nuclear Power Plants, J. R. Boardman; 166 Turbine Building Hazards, H. L Ornstein; RECENT DEVELOPMENTS: 169 Reports, Standards, and Safety Guides, D. S. Queener; 175 Proposed Rule Changes as of Dec. 31,1994; ANNOUNCEMENTS: 32 Harvard School of Public Health In-Place Filter Testing Workshop; 134 International Conference on Advances in the Operational Safety of Nuclear Power Plants; 193 30th Tennessee Industries Week; 193 DOE Technical Standards Program 1995 Workshop; 194 Multiphase Flow Experiments and Instrumentation; 180 The Authors; 185 Indexes to Nuclear Safety, Volumes 34 and 35.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

X-ray fluorescence mapping of brain tissue reveals the profound extent of trace element dysregulation in stroke pathophysiology

Abstract The brain is a privileged organ with regard to its trace element composition and maintains a robust barrier system to sequester this specialized environment from the rest of the body and the vascular system. Stroke is caused by loss of adequate blood flow to a region of the brain. Without adequate blood flow ischaemic changes begin almost immediately, triggering an ischaemic cascade, characterized by ion dysregulation, loss of function, oxidative damage, cellular degradation, and breakdown of the barrier that helps maintain this environment. Ion dysregulation is a hallmark of stroke pathophysiology and we observe that most elements in the brain are dysregulated after stroke. X-ray fluorescence-based detection of physiological changes in the neurometallome after stroke reveals profound ion dysregulation within the lesion and surrounding tissue. Not only are most elements significantly dysregulated after stroke, but the level of dysregulation cannot be predicted from a cell-level description of dysregulation. X-ray fluorescence imaging reveals that the stroke lesion retains <25% of essential K+ after stroke, but this element is not concomitantly elevated elsewhere in the organ. Moreover, elements like Na+, Ca2+, and Cl− are vastly elevated above levels available in normal brain tissue (>400%, >200%, and >150%, respectively). We hypothesize that weakening of the blood–brain barrier after stroke allows elements to freely diffuse down their concentration gradient so that the stroke lesion is in equilibrium with blood (and the compartments containing brain interstitial fluid and cerebrospinal fluid). The change observed for the neurometallome likely has consequences for the potential to rescue infarcted tissue, but also presents specific targets for treatment.

Biochemistry & Molecular Biology↗

DOE OSTI2024

Stability and transport of gyrokinetic critical pedestals

Abstract A gyrokinetic threshold model for pedestal width–height scaling prediction is applied to multiple devices. A shaping and aspect ratio scan is performed on National Spherical Torus Experiment (NSTX) equilibria, finding Δ ped = 0.92 A 1.04 κ − 1.24 0.38 δ β θ , ped 1.05 for the wide-pedestal branch with pedestal width Δ ped , aspect ratio A , elongation κ , triangularity δ , and normalized pedestal height β θ , ped . The width–transport scaling is found to vary significantly if the pedestal height is varied either with a fixed density or fixed temperature, showing how fueling and heating sources affect the pedestal density and temperature profiles for the kinetic-ballooning-mode (KBM) limited profiles. For an NSTX equilibrium, at fixed density, the wide branch is Δ ped = 0.028 ( q e / Γ e − 1.7 ) 1.5 ∼ η e 1.5 and at fixed temperature Δ ped = 0.31 ( q e / Γ e − 4.7 ) 0.85 ∼ η e 0.85 , where q e and Γ e are turbulent electron heat and particle fluxes and η e = ∇ ln ⁡ T e / ∇ ln ⁡ n e for an electron temperature T e and density n e . Pedestals close to the KBM limit are shown to have modified turbulent transport coefficients compared to the strongly driven KBMs. The role of flow shear is studied as a width–height scaling constraint and pedestal saturation mechanism for a standard and lithiated wide pedestal discharge. Finally, the stability, transport, and flow shear constraints are combined and examined for an NSTX experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Scanning Electrochemical Microscopy for Kinetic Investigations in Viscous Deep Eutectic Solvents: Identifying Practical Approach Curves and Deviations from Electron Transfer Models

Determining heterogeneous electrochemical electron transfer (ET) kinetics in electrolytes with a wide range of physical properties is of great interest for achieving high-performance redox flow batteries. Among such electrolytes, concentrated hydrogen-bonded electrolytes (CoHBEs), including deep eutectic solvents (DESs), have recently garnered significant attention. Unfortunately, traditional Tafel analysis using macroelectrodes often encounters issues with mass transfer limitations in CoHBEs with high viscosities, thereby restricting kinetic analysis to a narrow potential window. Here, in this work, we introduce a methodology for evaluating ET kinetics in viscous DES using the scanning electrochemical microscopy (SECM). We first determined practical solutions to SECM tip positioning in ethaline DES, which yield pseudopositive feedback responses. Lattice Boltzmann method (LBM) simulations helped us rationalize the impact of the fluid and concentration fields, as well as tip geometry, tip approach velocity v, and the solvent viscosity ηs, on the shape of the approach curves. In addition to successfully recreating approach curves over a variety of conditions, we found that approaching a conductor ensured a practical point where the normalized tip response (Ni T = 2) converged at L = 0.7 within ∼10% error regardless of tip velocity. With positioning capabilities at hand, we investigated the kinetics of Fe 3+ /Fe 2+ redox couple in aqueous and the ethaline media. The experimental kinetic results were interpreted using the Butler–Volmer (BV) and Marcus–Hush–Chidsey (MHC) models. For ethaline, a nonideal kinetic behavior was observed, potentially attributed to solvent dynamics within DESs or to the interplay of chloride anions in the charge transfer process.

electrodes↗

DOE OSTI2025

Investigating the relationship between simulation parameters and flow variables in simulating atmospheric gravity waves for wind energy applications

Wind farms, particularly offshore clusters, are becoming larger than ever before. Besides influencing the surface wind flow and the inflow for downstream wind farms, large wind farms can trigger atmospheric gravity waves in the inversion layer and the free atmosphere aloft. Wind-farm-induced gravity waves can cause adverse pressure gradients upstream of the wind farm, which contribute to the global blockage effect, and can induce favorable pressure gradients above and downstream of the wind farm that enhance wake recovery. Numerical modeling is a powerful means of studying these wind-farm-induced atmospheric gravity waves, but it comes with the challenge of handling spurious reflections of these waves from domain boundaries. Typically, approaches which employ radiation boundary conditions and forcing zones are used to avoid these reflections. However, the simulation setup of these approaches relies heavily on ad hoc processes. For instance, the widely used Rayleigh damping method requires ad hoc tuning to produce a setup that may only produce satisfactory results for a particular case. To provide more systematic guidance on setting up realistic simulations of atmospheric gravity waves, we conduct a large-eddy simulation (LES) study of flow over a 2D hill and through a wind farm canopy that explores the optimum domain size and damping layer setup depending on the fundamental parameters which determine the flow characteristics. In this work, we only consider linearly stratified conditions (i.e., no inversion layer), thereby focusing on internal gravity waves in the free atmosphere and their reflections from the domain boundaries. This type of flow is governed by a single Froude number, which dictates most of the internal wave properties, such as wavelength, amplitude, and direction. This, in turn, will dictate the optimum domain size and Rayleigh damping layer setup. We find the effective horizontal and vertical wavelengths (the representative wavelengths of the entire wave spectrum) to be the appropriate length scales to size the domain and damping layer thickness, and the optimal Rayleigh damping coefficient scales with the Brunt–Väisälä frequency. Considering Froude numbers seen in wind farm applications, we propose recommendations to limit the reflections to less than 10 % of the total upward-propagating wave energy. Typically, damping is done at the top boundary, but given the non-periodic lateral boundary conditions of practical wind farm simulation domains, we find that damping the inflow–outflow boundaries is of equal importance to damping the top boundary. The Brunt–Väisälä frequency-normalized damping coefficient should be between 1 and 10. The damping layer thickness should be at least one effective vertical wavelength; damping layers exceeding 1.5 times the vertical wavelength are found to be unnecessary. The domain length and height should accommodate at least one effective horizontal and vertical wavelength, respectively. Moreover, Rayleigh damping does not damp the waves completely, and the non-damped energy might accumulate over the simulation time.

17 WIND ENERGY↗

DOE OSTI2025

The Turbulent Pressure Spectrum Within the Roughness Sublayer of a Subarctic Forest Canopy

The turbulent static pressure spectrum E pp (k x ) as a function of longitudinal wavenumber k x in the roughness sublayer of forested canopies is of interest to a plethora of problems such as pressure transport in the turbulent kinetic energy budget, pressure pumping from snow or forest floor, and coupling between flow within and above canopies. Long term static pressure measurements above a sub-arctic forested canopy for near-neutral conditions during the winter and spring were collected and analyzed for three snow cover conditions: trees and ground covered with snow, trees are snow free but the ground is covered with snow, and snow free cover. In all three cases, it is shown that E pp (k x ) obeys the attached eddy hypothesis at low wavenumbers (1/δ < k x < 1/z) —with E pp (k x ) ∝ u$^4_*$k$^{-1}_x$ and Kolmogorov scaling in the inertial subrange at higher wavenumbers—with E pp (k x ) ∝ ϵ 4/3 k x -7/3 , where u * is the friction velocity at the canopy top, is the mean turbulent kinetic energy dissipation rate, z is the distance from the snow top, and is the boundary layer depth. The implications of these two scaling laws to the normalized root-mean squared pressure C p = o p /u$^2_*$ and its newly proposed logarithmic scaling with normalized wall-normal distance z/δ are discussed for snow covered and snow free vegetation conditions. The work here also shows that k$^{-1}_x$ the in the E pp (k x ) appears more extensive and robust than its longitudinal velocity counterpart.

Aslan, Toprak [Finnish Meteorological Inst., Helsi↗

DOE OSTI2024

Off-design operation and performance of pumped thermal energy storage

In this article, we describe off-design models and control strategies for a Pumped Thermal Energy Storage (PTES) system that uses liquid thermal energy storage: specifically molten salt for hot storage and methanol for cold storage. Off-design conditions arise when load-following, or due to variations in storage tank temperatures or ambient temperatures. We propose a control strategy that uses inventory control to manage the mass flow rate in the thermodynamic cycles, which facilitates load following. We also propose a control strategy for the storage fluid mass flow rates, which are varied to ensure the molten salt is maintained at its design temperature. This maximizes efficiency and minimizes problems with salt freezing or degradation. The cold storage fluid mass flow rate is varied so that the cold tanks have the same state-of-charge as the hot tanks. This leads to variations in cold fluid temperature, but these variations are shown to be acceptably small (e.g. 7.5% increase), and this control method is shown to be simpler and more efficient than an alternative strategy where tanks become unbalanced. The ambient temperature and storage tank temperatures are moved ±50 °C from the design values and the impact on power, duration, and tank temperatures is quantified. Results demonstrate that the proposed control strategy is stable and self-correcting - that is, storage temperatures converge on stable values after two-to-three charge-discharge cycles. When inputs return to design values, the system returns to its design point after two charge-discharge cycles. We also demonstrate that inventory control enables delivery of the target power output even when off-design conditions exist that would normally reduce the power output.

25 ENERGY STORAGE↗

DOE OSTI2025

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Impacts of Irradiation Structural Behavior on Thermal Hydraulics Safety Analysis to Support MURR LEU Conversion

The University of Missouri Research Reactor (MURR) located in Columbia, Missouri is one of six U.S. High Performance Research Reactors (USHPRR), including one critical facility, that is actively collaborating with U. S. Department of Energy (DOE) National Nuclear Security Administration (NNSA) Material Management and Minimization (M3) Office of Reactor Conversion and Uranium Supply to convert from the use of highly enriched uranium (HEU; ≥20 wt% U-235) to low-enriched uranium (LEU; <20 wt% U-235) fuel. A new type of very high-density LEU fuel based on an alloy of uranium and 10 wt% molybdenum (U-10Mo) is expected to allow the conversion to LEU of MURR, as well as four other USHPRR. MURR has been working with the Reactor Conversion Pillar at Argonne to perform fuel element design and fuel cycle performance analyses, steady-state thermal hydraulics safety analysis, and accident safety analyses in preparation for the conversion of MURR and to support a preliminary safety analysis report for conversion to LEU fuels. Subsequent analyses have also been performed, including transition cycles where all-fresh LEU fuel elements are introduced upon conversion and progressing through reactor operations the core is brought to equilibrium. Thermal hydraulics safety analyses performed as part of the above have employed an assumption on channel gap reduction due to burnup-related phenomena including fuel swelling, irradiation creep, and oxide layer buildup. Recently, a series of structural analyses have been performed on the MURR LEU fuel plates and an element due to significant differences between the plate and element designs of the MURR HEU and LEU fuels. In addition, NUREG-1537 indicates that structural phenomena are to be evaluated. Two separate types of structural analyses were performed for the MURR LEU fuel element: fluid-structure interaction (FSI) and irradiation thermo-mechanical. The FSI analysis evaluated the effects of hydraulic forces on the MURR LEU fuel element to quantify the flow-induced plate deflection, and a minimal impact to the channel gap thickness was predicted under prototypic and bounding conditions. The irradiation thermo-mechanical analysis evaluated the effects of fuel swelling, irradiation creep, and thermal expansion for the MURR LEU plates and the element for prototypic thermal and irradiation conditions based on a high-fidelity approach multiphysics approach. Overall, this thermo-mechanical analysis predicts smaller gap thickness changes in previously limiting regions. Larger changes are predicted in the middle of channels, and for end channels where power density is not typically a maximum. An additional thermo-mechanical analysis was performed for the outermost HEU fuel plate, which showed a similar magnitude of deflection as the outermost LEU plate. Due to substantial differences between the channel gap reductions assumed for the previous safety analyses and those predicted by the irradiation thermo-mechanical analysis, a need to evaluate their impact on the thermal hydraulics safety analyses arose. This report presents the results from the steady-state safety analyses for normal operation as well as the accident analyses for the two most limiting accident scenarios.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

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↗