SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1990

The Importance of the Measurement of the Unsteady Wake of the Aeroassist Flight Experiment Vehicle: A White Paper

The High Frequency Radiometer (HFR) is the only instrument on the Aeroassist Flight Experiment (AFE) with sufficient temporal resolution to discern the frequency of unsteady wake oscillations. Determining both the frequency and amplitude of wake unsteadiness during AFE atmospheric entry is essential for reliably predicting the geometry and motion of the wake of future Aeroassisted Space Transfer Vehicles (ASTV). These parameters directly affect the location and size of the pay- load and the weight of the required afterbody heat protection. The purpose of the AFE is to validate the technologies required for the design of ASTVs, which will be used to exploit Earth-lunar space. This validation will be conducted at a combination of vehicle size, altitude, and velocity not obtainable in ground-based facilities. The AFE will provide the experimental flight data needed to improve our understanding of hypersonic-wake physics and to validate computational predictions of the aerodynamic and heating loads, including afterbody radiative heating loads, on an ASTV. Reliable prediction of ASTV wake flows will ensure that payloads are located within the shear-layer envelope and will determine the amount of thermal protection the payloads require. Specifically, understanding the temporal nature of the wake unsteadiness is important for two reasons. Most importantly, analysis of ground-based experiments suggests that wake unsteadiness results in a variation of as much as +/- 5 deg in the shear-flow turning angle. This angle must be reliably predicted to avoid shear-layer impingement on the vehicle afterbody, which would result in heating rates of about 10 W/sq cm, of the same order as on the forebody stagnation point. Secondly, the energy associated with wake unsteadiness will reduce the static enthalpy of the wake fluid and cause an error of as much as 30% in the amount of predicted wake radiative heating. Therefore, the HFR flight data, which will quantify the frequency and amplitude of the wake unsteadiness, are required for the verification of computational models of ASTV flowfields. Measurement of wake unsteadiness has been endorsed by the Peer Science Steering Group and the AFE Computational Fluid Dynamics (CFD) Working Group. The Peer Science Steering Group has stressed that a reliable measurement of the wake unsteadiness is fundamental because incorporating unsteadiness represents a substantial challenge to the CFD community, and a reliable measurement will raise confidence in the computer simulation. The AFE CFD Working Group has stated that 44 any data which would prove the existence and effects of unsteady flow would be extremely valuable.

Strawa, A. W.↗

DOE OSTI2022

Development, construction and qualification tests of the mechanical structures of the electromagnetic calorimeter of the Mu2e experiment at Fermilab

The “muon-to-electron conversion” (Mu2e) experiment at Fermilab will search for the Charged Lepton Flavour Violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of physics beyond the Standard Model. Mu2e detectors comprise a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. The calorimeter provides excellent electron identification, complementary information to aid pattern recognition and track reconstruction, and a fast calorimetric online trigger. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1340 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named “disks”) positioned at the relative distance of 70 cm downstream the aluminum target along the muon beamline satisfies the Mu2e physics requirements.The hostile Mu2e operational conditions, in terms of radiation levels (total ionizing dose of 12 krad and a neutron fluence of 5x1010 n/cm2 @ 1 MeVeq (Si)/y), magnetic field intensity (1 T) and vacuum level (10$^{-4}$ Torr) have posed tight constraints on the design of the detector mechanical structures and materials choice. The support structure of the two 670 crystal matrices employs two aluminum hollow rings and parts made of open-cell vacuum-compatible carbon fiber. The photosensors and service front-end electronics for each crystal are assembled in a unique mechanical unit inserted in a machined copper holder. The 670 units are supported by a machined plate made of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal contact with the copper holders to constitute a low resistance thermal bridge. The data acquisition electronics is hosted in aluminum custom crates positioned on the external lateral surface of the two disks. The crates also integrate the electronics cooling system as lines running in parallel to the front-end system.The constraints on the calorimeter mechanical structures design, the development from the conceptual design to the specifications of all the structural components, the status of components production and the components quality assurance tests are presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

FENIX: Towards a Fully Integrated Multiphysics Framework for Plasma Facing Component Modeling

Computational tools have a crucial role to play in accelerating the deployment of fusion as a clean, reliable, abundant, and sustainable energy source. Multiphysics, high-fidelity simulation capabilities can help model, study, and predict intricate interactions between materials performance, plasma exposure, neutron irradiation, and engineering processes. As such, they can assist in the resolution of scientific and engineering challenges underpinning design, construction, and commission of fusion power plants. To address these needs, ongoing efforts are leveraging the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework and delivering new computational tools for the fusion community. These tools inherit crucial attributes from MOOSE. They are open-source, modular, integrated with nuclear industry-standard software quality assurance processes, and enable multiphysics, multi-fidelity, fully integrated, zero- to three-dimensional, and massively parallel simulations. After a short overview of these capabilities, we will present the development of Fusion ENergy Integrated multiphys-X (FENIX), a MOOSE-based application designed to enable plasma facing component design and performance evaluation. Throughout their lifetime, plasma facing components are exposed to extreme thermal loads, repeated thermal shocks, and irradiation by plasma ions, neutral particles, and high-energy neutrons. Consequently, designing a plasma facing component with acceptable lifetime degradation is extremely challenging. FENIX aims to model the multiphysics environment in which plasma facing components evolve to accelerate their design studies. To that end, FENIX couples existing MOOSE capabilities such as heat transfer, thermomechanics, and thermal hydraulics, with tritium transport via the MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8), with neutronics via the MOOSE-based high-fidelity neutron-photon transport and fluid dynamics code Cardinal, and finally with Particle-in-Cell plasma simulation capabilities being developed in this project. In this study, we present the current FENIX capabilities and preliminary results of its application to model the Tritium Plasma Experiment set up at Idaho National Laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗

DOE OSTI2025

Autonomous sensor suite for evaluating fish-turbine interactions and environmental impacts in marine renewable energy and hydropower

Marine renewable energy (MRE) harnesses ocean-based resources such as waves, tides, currents, and thermal or salinity gradients for sustainable power generation. It has the potential to complement existing renewable resources, support remote communities, and contribute to decarbonization efforts. However, understanding the hydrodynamic forces created by MRE devices and their impacts on marine life is critical for responsible deployment. Here, to address these concerns, advanced sensor devices, including the Marine Sensor Fish (MSF), Sensor Fish Mini (SF Mini), and Flexible Sensor Fish (FSF), were developed to measure interactions between aquatic organisms and MRE systems. This paper details the design, manufacturing, calibration, and field deployment of these sensor suites, highlighting their ability to capture key physical stressors such as shear forces, pressure changes, and collision impacts. The MSF successfully evaluated turbine interactions at a tidal turbine in the Salish Sea, capturing data on turbulence, collision impact, and pressure gradients. The SF Mini validated hydrodynamic conditions in scaled hydraulic models, supporting computational fluid dynamics simulations. The FSF, with its flexible silicone body, measured species-specific impacts in turbulent environments. This research demonstrates the potential of Sensor Fish technology to advance sustainable marine energy systems by reducing biological impacts and informing environmentally sustainable designs.

Ecological impacts↗

DOE OSTI2020

Mercury speciation using microcolumns with a direct mercury analyzer (DMA). Development of analytical method for the Savannah River Site liquid system

Mercury in the Savannah River Site (SRS) Liquid Waste System (LWS) exists in various forms, including: a) ionic inorganic mercury, organomercury (e.g., methyl Hg), and other less abundant species dissolved in LWS fluids, b) mercury solids such as oxides, hydroxides, amalgams, sulfides, and sorbed mercury, c) accumulations of dense liquid elemental mercury, and d) vapor phase elemental and organomercury mercury in tank headspace gas and in evaporators. Strategic and proactive management of the estimated 60,000 Kg mercury in the LWS requires deployment of efficient and effective paradigms for sampling and analysis of the total quantity of mercury in the various physical-chemical forms. Scoping tests of selective gas-phase sorbents for mercury speciation were performed to support LWS objectives. This research specifically focused on developing and testing streamlined methods to obtain high-quality analytical results for gaseous mercury species. The proposed methods are straightforward, using small columns (microcolumns) filled with materials that selectively sorb one or more target forms of mercury. Simple differentiation of mercury species and quantification is then achieved using an efficient thermal-desorption-based total mercury analyzer as the final step. Benefits associated with transitioning from multistep procedures such as EPA Methods 245.1, 1630 and 1631to streamlined microcolumn approaches include reduced analysis time, labor, and waste generation. Such transitioning would maximize the value of the current deployment of direct mercury analysis (EPA Method 7473) at SRS. This work has two target end uses: 1) simplifying the analysis steps and reducing costs for liquid samples (e.g., developing a total organic mercury (TOM) method), and 2) providing options for high-quality onsite analysis and rapid turnaround for gas samples in support of SRS LWS industrial hygiene (IH) objectives.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2024

Design and assembling status of the Mu2e electromagnetic calorimeter mechanical structures

The muon-to-electron conversion (Mu2e) experiment at Fermilab will search for the Charged Lepton Flavour Violating neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the unambiguous evidence of physics beyond the Standard Model. The detector has been designed as a state-of-the-art crystal calorimeter and employs 1348 pure Cesium Iodide (CsI) crystals readout by UV-extended silicon photosensors and fast front-end and digitization electronics. A design consisting of two identical annular matrices (named disks ) positioned at the relative distance of 70 cm, downstream the aluminum target along the muon beamline, satisfies the Mu2e physics requirements. The hostile Mu2e operational conditions, in terms of radiation levels (total ionizing dose of 12 krad and a neutron fluence of 5 10^{10} n/cm2 @ 1 MeVeq (Si)/y), magnetic field intensity (1 T) and vacuum level (10−4 Torr) have posed tight constraints on the design of the detector mechanical structures and materials choice. The support structure of the two 674 crystal matrices employs two aluminum hollow rings and parts made of open- cell vacuum-compatible carbon fiber. The photosensors and service front-end electronics for each crystal are assembled in a unique mechanical unit inserted in a machined copper holder. The 674 units are supported by a machined plate made of vacuum-compatible plastic material. The plate also integrates the cooling system made of a network of copper lines flowing a low temperature radiation-hard fluid and placed in thermal con- tact with the copper holders. The data acquisition electronics is hosted in aluminum custom crates positioned on the external lateral surface of the two disks. The crates also integrate the electronics cooling system. In this poster we will review the constraints on the calorimeter mechanical structures, the technological choices, and the status of assembling at Fermilab.

Pasciuto, D.↗

DOE OSTI2022

Observing flow of He II with unsupervised machine learning

Abstract Time dependent observations of point-to-point correlations of the velocity vector field (structure functions) are necessary to model and understand fluid flow around complex objects. Using thermal gradients, we observed fluid flow by recording fluorescence of $${\text{He}}_{2}^{*}$$ He 2 ∗ excimers produced by neutron capture throughout a ~ cm 3 volume. Because the photon emitted by an excited excimer is unlikely to be recorded by the camera, the techniques of particle tracking (PTV) and particle imaging (PIV) velocimetry cannot be applied to extract information from the fluorescence of individual excimers. Therefore, we applied an unsupervised machine learning algorithm to identify light from ensembles of excimers (clusters) and then tracked the centroids of the clusters using a particle displacement determination algorithm developed for PTV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1982

Thresholds for the onset of fluid and magnetofluid turbulence

Linear stability calculations conducted in plasma physics are based on the theory of hydrodynamic stability of neutral fluids. The present investigation is concerned with the validity of procedures based on linear stability analysis in six much-studied situations. All have the property that the fluid goes from laminar to turbulent at critical values of some dimensionless number, such as the Reynolds number or the Rayleigh number. The question is, at what threshold values of the dimensionless number do the unstable motions set in, and can these thresholds be predicted by a linear analysis of the stability of the laminar state as the threshold is approached from the stable side. In three cases linear stability analysis clearly seems to fail. These cases include Plane Poiseuille Flow, Plane Couette Flow, and Cylindrical Pipe Flow (Hagen-Poiseuille Flow). Situations in which predictions provided by linear stability analysis are correct are related to Rotating Couette Flow, Thermally-Driven Convection, and Instability of Laminar Boundary Layers.

Montgomery, D.↗

DOE OSTI2025

Dielectron production in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02TeV

The first measurement of the 𝑒 + ⁢𝑒 − pair production at midrapidity and low invariant mass in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider is presented. The yield of 𝑒 + ⁢𝑒 − pairs is compared with a cocktail of expected hadronic decay contributions in the invariant mass (𝑚𝑒⁢𝑒) and pair transverse momentum (𝑝 T,𝑒⁢𝑒 ) ranges 𝑚 𝑒⁢𝑒 < 3.5 GeV/𝑐 2 and 𝑝 T,𝑒⁢𝑒 < 8 GeV/𝑐. For 0.18 < 𝑚 𝑒⁢𝑒 < 0.5 GeV/𝑐 2 the ratio of data to the cocktail of hadronic contributions amounts to 1.40 ± 0.11⁢(stat.) ± 0.23 ⁢(syst.) ± 0.16⁢ (cocktail) and 1.42 ± 0.11⁢ (stat.) ± 0.23 (syst.)$^{+0.24}_{−0.29}$ ⁢(cocktail), including or not including medium effects in the estimation of the heavy-flavor background, respectively. It is consistent with predictions from two different models for an additional contribution of thermal 𝑒+⁢𝑒− pairs from the hadronic and partonic phases. In the intermediate-mass range (1.2 < 𝑚 𝑒⁢𝑒 < 2.6 GeV/𝑐 2 ), the pair transverse impact parameter of the 𝑒 + ⁢𝑒 − pairs (DCA𝑒⁢𝑒, where “DCA” denotes “distance of closest approach”) is used for the first time in Pb-Pb collisions to separate displaced dielectrons from heavy-flavor hadron decays from a possible (thermal) contribution produced at the interaction point. The data are consistent with a suppression of 𝑒 + ⁢𝑒 − pairs from $𝑐\bar{𝑐}$ and an additional prompt component. Finally, the first direct-photon measurement in the 10% most central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV is reported via the study of virtual direct photons in the transverse momentum range 1 < 𝑝 T < 5 GeV/𝑐. A model including prompt photons, as well as photons from the preequilibrium and fluid-dynamic phases, can reproduce the result, while being at the upper edge of the data uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1994

Mode properties of low-frequency waves: Kinetic theory versus Hall-MHD

In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.

Krauss-Varban, D.↗

NASA NTRS2012

The Spatial Distribution of Organic Matter and Mineralogical Relationships in Carbonaceous Chondrites

Organic matter present within primitive carbonaceous meteorites represents the complex conglomeration of species formed in a variety of physically and temporally distinct environments including circumstellar space, the interstellar medium, the Solar Nebula & Jovian sub-nebulae and asteroids. In each case, multiple chemical pathways would have been available for the synthesis of organic molecules. Consequently these meteorites constitute a unique record of organic chemical evolution in the Universe and one of the biggest challenges in organic cosmochemistry has been in deciphering this record. While bulk chemical analysis has provided a detailed description of the range and diversity of organic species present in carbonaceous chondrites, there is virtually no hard experimental data as to how these species are spatially distributed and their relationship to the host mineral matrix, (with one exception). The distribution of organic phases is nevertheless critical to understanding parent body processes. The CM and CI chondrites all display evidence of low temperature (< 350K) interaction with aqueous fluids, which based on O isotope data, flowed along thermal gradients within the respective parent bodies. This pervasive aqueous alteration may have led to aqueous geochromatographic separation of organics and synthesis of new organics coupled to aqueous mineral alteration. To address such issues we have applied the technique of microprobe two-step laser desorption / photoionization mass spectrometry (L2MS) to map in situ the spatial distribution of a broad range of organic species at the micron scale in the freshly exposed matrices of the Bells, Tagish Lake and Murchison (CM2) carbonaceous chondrites.

Clemett, S. J.↗

NASA NTRS1999

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗

DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Symmetric Set of Transport Coefficients for Collisional Magnetized Plasma

Braginskii extended magnetohydrodynamics is used to model transport in collisional astrophysical and high energy density plasmas. We show that commonly used approximations to the α ⊥ and β ⊥ transport coefficients [e.g., Epperlein and Haines, Phys. Fluids 29 , 1029 (1986) ] have a subtle inaccuracy that causes significant artificial magnetic dissipation and discontinuities. This is because magnetic transport actually relies on β ∥ - β ⊥ and α ⊥ - α ∥ , rather than α ⊥ and β ⊥ themselves. We provide fit functions that rectify this problem and thus resolve the discrepancies with kinetic simulations in the literature. Furthermore, then implemented in the g orgon code, they reduce the predicted density asymmetry amplitude at laser ablation fronts. Recognizing the importance of α ⊥ - α ∥ and β ∥ - β ⊥ , we recast the set of coefficients. This makes explicit the symmetry of the magnetic and thermal transport, as well as the symmetry of the coefficients themselves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗

DOE OSTI2024

Minimizing Toxicity and Optimizing Lubricity of Ionic Liquids for Eco-Friendly Lubrication

As much as 60 million liters of lubricating fluids end up in the environment annually, and thus, the environmental impact of lubricants is increasingly recognized in addition to meeting the rheological and tribological requirements. Although US Environmental Protection Agency-approved environmentally acceptable base fluids are available, there is a lack of additives that are both nontoxic and effective in wear protection. Here, this study reports the successful development of a new class of ionic liquids (ILs) with demonstrated significantly lower aquatic toxicity and superior friction reduction and wear protection capabilities compared with a commercial lubricant additive and some ILs reported in the literature. Specifically, ammonium phosphate and phosphonium phosphate ILs with four-carbon alkyls have been identified with balanced oil solubility, thermal stability, toxicity, and lubricity, which provide fundamental insights for future development of eco-friendly ILs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗