Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Physical quantities”

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 595 records · Page 33

Machine-Learned Manifold-Based Models for Large Eddy Simulation of Turbulent Combustion

Reduced-order manifold combustion models are commonly used to lower the cost of reacting Large Eddy Simulations (LES) and may be derived either from physical principles as in flamelet models or from data-driven methods like Principal Component Analysis (PCA). In either case, neural networks are increasingly used as part of these models to provide a nonlinear mapping between a small set of pre-defined variables that parameterize the manifold and outputs of interest, such as reaction rates. In this work, we propose a new manifold-based modeling approach that combines the definition of the manifold-parameterizing variables (linear combinations of species), the nonlinear mapping to the outputs, and closure of filtered quantities for LES into the structure of a single neural network. This allows the process used to train the neural network to simultaneously optimize both the functional form of the model and the identities of the inputs to the model. The new approach can flexibly incorporate thermochemical data from any combustion system; if trained on data from 1D flames it can be interpreted as an optimized flamelet model, but it can also be used to learn models from data from more complex configurations. This work presents a priori evaluations of the new approach in both contexts. Evaluation using data from 1D premixed flames demonstrates the physical interpretability of the manifold variables generated by the new approach. Evaluation using data from direct numerical simulations of turbulent flames shows improved predictions relative to either flamelet or PCA-based models in a more complex configuration.

47 OTHER INSTRUMENTATION↗

Model-Free Measurement of Local Entropy Production and Extractable Work in Active Matter

Time-reversal symmetry breaking and entropy production are universal features of nonequilibrium phenomena. Despite its importance in the physics of active and living systems, the entropy production of systems with many degrees of freedom has remained of little practical significance because the high dimensionality of their state space makes it difficult to measure. Here, in this work, we introduce a local measure of entropy production and a numerical protocol to estimate it. We establish a connection between the entropy production and extractability of work in a given region of the system and show how this quantity depends crucially on the degrees of freedom being tracked. We validate our approach in theory, simulation, and experiments by considering systems of active Brownian particles undergoing motility-induced phase separation, as well as active Brownian particles and E.coli in a rectifying device in which the time-reversal asymmetry of the particle dynamics couples to spatial asymmetry to reveal its effects on a macroscopic scale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Informing nuclear physics via machine learning methods with differential and integral experiments

Information from differential nuclear-physics experiments and theory is often too uncertain to accurately define nuclear-physics observables such as cross sections or energy spectra. Integral experimental data, representing the applications of these observables, are often more precise but depend simultaneously on too many of them to unambiguously identify issues in the observable with human expert analysis alone. Here, we explore how we can leverage physics knowledge gained from differential experimental data, nuclear theory, integral experiments, and neutron-transport calculations to better understand nuclear-physics observables in the context of the application area represented by integral experiments. We support this task with machine-learning methods to discern trends in a large amount of convoluted data. Differential and integral information was used in an analysis augmented by the random forest and the Shapley additive explanations metric. We chose as an application area one that is represented by criticality measurements and pulsed-sphere neutron-leakage spectra. We show one representative example ( 241 Pu fission observables) where the combination of differential and integral information allowed to resolve issues in data representing these observables. As a starting point, the machine learning (ML) algorithms highlighted several observables as leading potentially to bias in simulating integral experiments. Differential information, paired with sensitivity to integral quantities, allowed us then to pinpoint one specific observable ( 241 Pu fission cross section) as the main driver of bias. The comparison to integral experiments, on the other hand, allowed us to indicate a likely reliable experiment among several discrepant ones for this observables. In other cases (e.g., 239 Pu observables), we were not able to resolve the confounding introduced by integral experiments but instead highlighted the need for targeted new experiments and theory developments to better constrain the nuclear-physics space for the application area represented by integral experiments. We were able to combine information from differential experimental data, nuclear-physics theory, integral experiments, and neutron-transport simulations of the latter experiments with the help of the random forest algorithm and expert judgment. This combination of knowledge allows to improve our description of nuclear-physics observables as applied to a particular application area.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SCALE Lattice Physics Code Assessments of Accident Tolerant Fuel

This report highlights accident-tolerant fuel (ATF) code assessment activities performed under Project NRCHQ- 60-17-T-0017, Lattice Physics Enhancements and Assessment. ATF covers a broad range of advanced fuel and clad designs for light water reactors (LWRs) to enhance performance under several accident conditions. Several ATF concepts are anticipated to be deployed as lead test rods (LTRs) or lead test assemblies (LTAs) within the next five years. The purpose of this work is to assess the predictive capabilities of NRC neutronics codes that underpin various licensing calculations. ATF designs use different fuel and clad materials compared to standard UO 2 fuel and zirconium-alloy claddings (hereafter UO 2 -Zry). These new materials and geometrical designs need to be assessed to quantify the impact of nuclear data uncertainties on quantities of interest (QOIs) in licensing calculations and the impact of modeling approximations which may be valid for UO 2 -Zry but not for ATF. This report outlines a systematic approach for ATF neutronics code assessment which includes sensitivity and uncertainty (S/U) analysis of nuclear data, identification of experimental benchmark and gaps for code validation, investigation of modeling approximations, and code-to-code comparisons of calculated QOIs against high-fidelity reference continuous energy (CE) Monte Carlo (MC) calculations. This report focuses on the assessment of the SCALE/Polaris lattice physics code for reactor safety analysis. Polaris lattice physics calculations generate few group (FG) cross sections for PARCS full-core calculations. (Full-core analysis, spent fuel analysis, and severe accident analysis will be performed in future work.) The selected ATF concepts for this report include Cr 2 O 3 and Al 2 O 3 -Cr 2 O 3 -doped UO 2 fuel, U 3 Si 2 fuel, FeCrAl cladding, SiC cladding, and Cr-coated cladding.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Select Proliferation Studies on TRISO-fueled, Heat-Pipe-cooled Microreactors

Nuclear microreactors carry the potential to open up new markets for the nuclear industry, as their expected cost competitiveness in non-traditional market segments (e.g., mines, military bases, extraterrestrial surfaces, and remote areas), and their inherent safety features make them deployable when other power sources are unavailable or difficult to exploit. For countries that have not traditionally participated in nuclear power, microreactors represent a clean energy solution [1]. However, their use, especially in non-weapons states, may entail challenges in terms of maintaining international nuclear safeguards [2]. Furthermore, the deployment locations where microreactors may prove most cost competitive would be difficult to access by state and International Atomic Energy Agency (IAEA) inspectors [3]. In addition to the isolated nature of potential deployment sites, the low-power characteristic of microreactors suggests that numerous microreactors would need to be deployed to meet energy demands. That, coupled with the unique physics of many current microreactor designs, opens up a new area of research with respect to nonproliferation and safeguards concerns [2]. Whereas traditional facilities are inspected as isolated cases when looking for signs of diversion or misuse; microreactors may need to be assessed in the context of the entire fleet to which they belong. International safeguards necessitate timely detection of any significant quantities (SQs) of material that are being diverted (e.g., 1 SQ of special nuclear material diverted over the course of a 1-year period) [4]. For low-enriched uranium, the IAEA defines 1 SQ as corresponding to 75 kg of 235U. The purpose of the present paper is to explore the detectability threshold for material diversion in microreactors by relying on critical control drum angles, excess reactivity, and the reactor lifetime as the selected operational parameters. For this assessment, a heat-pipe-cooled microreactor was regarded as the base design. While the conclusions reached in this paper are not readily extendable to the design of actual microreactors, the analysis herein enables conclusions to be drawn regarding the level of accuracy needed for reference calculations in order to detect diversion scenarios by utilizing the selected operational parameters (i.e., mainly control drum angles, critical insertion angle, and the reactor lifetime).

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct Monte Carlo Simulations of Chemical Equilibrium Composition of Molecular Fluid Mixtures Under Shock Conditions

A new Monte Carlo simulation method is implemented which holds the number and kind of atoms constant rather than the number and kind of molecules for a molecular fluid mixture. Chemical reactions are allowed as a correlated move of atoms such as an interchange of atoms between molecules. Equilibrium chemical composition is determined as an average over the simulation along with equation of state quantities. Results are presented for N 2 + O 2 → 2NO at 30GPa and 3000K. The equilibrium composition is shown to be sensitive to the potentials between unlike species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic evidence for scission neutrons

We discuss the neck rupture stage of fission and the emission of particles afterwards, two extremely rapid and highly non-equilibrium processes. Currently, the neck rupture cannot be directly probed by experiment, highlighting the importance of reliable theoretical predictions for this stage of fission. Time-dependent density functional theory (TDDFT) is used to simulate the spontaneous fission of 252 Cf. In conjunction with statistical models, inputs from microscopic calculations can be used to make predictions for prompt neutron and gamma emission spectra, quantities which can be experimentally measured. Here we characterize the scission mechanism within TDDFT and estimate the number of scission neutrons and their kinetic energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rational surfaces, flows and radial structure in the TJ-II stellarator

Abstract In this work, we report on the results obtained by measuring several turbulent quantities well inside the plasma edge by means of a Langmuir probe during dynamical rotational transform scans in the TJ-II stellarator, while applying a radial electric field to the edge plasma using a biasing probe. By calculating the intermittence parameter from floating potential measurements, we are able to identify a major low order rational surface and hence relate the probe measurements to the local value of the rotational transform. Based on the former, we are able to show that the poloidal plasma velocity (and hence radial electric field) has a significant radial structure that is clearly related to the rotational transform profile and in particular the lowest order rational surfaces in the range studied. The poloidal velocity is also affected by the edge biasing. The particle flux Γ was also found to exhibit a radial pattern, as did the flow shear suppression term ω E × B , but the relation of the former to the low-order rational surfaces was less clear. We surmise that this lack of direct correspondence is due to an unknown term in the turbulence evolution equation: the instability growth rate, γ . We make use of a reduced Magnetohydrodynamic turbulence model to interpret the results. Overall, a picture is obtained in which the plasma self-organizes towards a state with a clear radial pattern of the radial electric field, in line with expectations from some numerical studies describing the spontaneous formation of an ‘ E × B staircase’, consisting of alternating layers with fast and slow radial transport. In this state, the radial profiles of various quantities (density, temperature, pressure) will not be smooth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Silicon-Equivalent Neutron Dose Estimations

Radiation transport calculations are often used to estimate dose or exposure to components and personnel surrounding a radiation source. The sources for these calculations are decaying radionuclides within the various nuclear materials. Under long-term radiation exposure, some components may show signs of material damage such as hardening, cracking, discoloration, etc. resulting from the dose received from nearby radioactive sources. Radiation dose quantities are a useful means to correlate material damage or degradation to radiation interactions. The degradation is often dependent on how incident radiation interacts and deposits energy into the material. This dependence implies that radiation dose should be properly defined for the material effect of interest. Generic dose quantities, such as total energy deposited, may not necessarily correlate well for degradation of any given material.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Deep convolutional neural networks for multi-scale time-series classification and application to disruption prediction in fusion devices

The multi-scale, mutli-physics nature of fusion plasmas makes predicting plasma events challenging. Recent advances in deep convolutional neural network architectures (CNN) utilizing dilated convolutions enable accurate predictions on sequences which have long-range, multi-scale characteristics, such as the time-series generated by diagnostic instruments observing fusion plasmas. Here we apply this neural network architecture to the popular problem of disruption prediction in fusion tokamaks, utilizing raw data from a single diagnostic, the Electron Cyclotron Emission imaging (ECEi) diagnostic from the DIII-D tokamak. ECEi measures a fundamental plasma quantity (electron temperature) with high temporal resolution over the entire plasma discharge, making it sensitive to a number of potential pre-disruptions markers with different temporal and spatial scales. Promising, initial disruption prediction results are obtained training a deep CNN with large receptive field ({$\sim$}30k), achieving an $F_1$-score of {$\sim$}91\% on individual time-slices using only the ECEi data.

convolutional neural networks↗

Dark matter capture in celestial objects: treatment across kinematic and interaction regimes

Signatures of dark matter in celestial objects have become of increasing interest due to their powerful detection prospects. To test any of these signatures, the fundamental quantity needed is the rate in which dark matter is captured by celestial objects. Depending on whether dark matter is light, heavy, or comparable in mass to the celestial-body scattering targets, there are different considerations when calculating the capture rate. Furthermore, if dark matter has strong or weak interactions, the physical behaviour important for capture varies. Using both analytic approximations and simulations, here we demonstrate how to treat dark matter capture in a range of celestial objects for arbitrary dark matter mass and interaction strength. We release our calculation framework as a public package available in both Python and Mathematica versions, called Asteria.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalized Hydrodynamics: A Perspective

Conventional hydrodynamics describes systems with few long-lived excitations. In one dimension, however, many experimentally relevant systems feature a large number of long-lived excitations and conserved quantities even at high temperature, because they are proximate to integrable limits. Such models cannot be treated using conventional hydrodynamics. The framework of generalized hydrodynamics (GHD) was recently developed to treat the dynamics of one-dimensional models: It combines ideas from integrability, hydrodynamics, and kinetic theory to come up with a quantitative theory of transport. GHD has successfully settled several long-standing questions about one-dimensional transport; it has also been leveraged to study dynamical questions beyond the transport of conserved quantities and to systems that are not integrable. In this article, we introduce the main ideas and predictions of GHD, survey some of the most recent theoretical extensions and experimental tests of the GHD framework, and discuss some open questions in transport that the GHD perspective might elucidate. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗

Thermochemical Data Fusion Using Graph Representation Learning

Large databases are required for “Big Data” applications in catalysis and materials science. Thermochemical databases can be created by combining data from various sources and by correcting low-fidelity datasets to higher accuracy with minimal computation. To achieve this “data fusion”, thermochemical quantities of interest, calculated at various levels of density functional theory (DFT), need to be mapped to the same, high levels of theory. In this work, a graph theoretical, statistical framework is proposed for such tasks. Subgraph frequencies are shown to provide a natural representation for learning these fusion maps. The maps are linear and are learnt with automated descriptor selection. Using a dataset of as few as ~1% from the QM9 database of 133,885 molecules, these models can predict multiple thermochemical quantities at a higher level of theory with an accuracy of 1 kcal/mol. Here, the method is explainable, generalizable, and provides a diagnostic tool for outlier identification

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Preparation and Characterization of Model Homotopic Catalysts: Rh Adatoms, Nanoparticles, and Mixed Oxide Surfaces on Fe 3 O 4 (001)

The atomic-level characterization of active sites is essential to understanding the mechanisms behind catalytic reactions. In this study, using scanning tunneling microscopy (STM) and X-ray photoelectron spectroscopy (XPS), we follow the morphological changes of a model Rh catalysts supported on Fe 3 O 4 (001) as a function of temperature and Rh coverage. Here, we identify the preparation conditions leading to model catalysts containing homotopic or nearly homotopic Rh species bound as adatoms, small clusters, substitutional within the Fe 3 O 4 (001), and nanoparticles. Adsorbates such as CO and CO 2 are further subsequently used to characterize the properties of different Rh sites. Using temperature programmed desorption (TPD), we demonstrate that adatoms, clusters, and nanoparticles exhibit high-temperature CO desorption (250-600 K). Strong binding on such sites further allows for CO oxidation to CO 2 via the Mars-van-Krevelen mechanism. In contrast, CO 2 was found to interact weakly with all Rh sites. Differences in desorption temperature enable the use of CO and CO 2 as titration methods for nanoparticles and Fe 3 O 4 (001), respectively. A small quantity of CO 2 was found to be reduced to CO on Rh adatoms and clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Beam correction for multi-pass arcs in FFA@CEBAF: status update

This work examines the multi-pass steering of six electron beams in an FFA arc ranging from approximately 10.5 GeV to 22 GeV. Shown here is an algorithm based on singular value decomposition (SVD) to successfully steer all six beams through the arc given precise knowledge of all beam positions at each of one hundred and one diagnostic locations with one hundred individual corrector magnets: that is successive application of SVD to different 100 × 101 response matrices—one for each beam energy. Further, a machine learning scheme is developed which only requires knowledge of the energy-averaged beam position at each location to provide equivalent steering. Extension of this scheme to other beam optics quantities as well as transverse and longitudinal coupling is explored.

Accelerator Physics↗

R-matrix Resolved Resonance Region Evaluation of 140,142 Ce

Oak Ridge National Laboratory completed the resolved resonance region (RRR) evaluation of the two most abundant cerium isotopes, 140 Ce (88.45%) and 142 Ce (11.11%), as requested by the US Nuclear Criticality Safety Program. These evaluations are based on recent high-resolution transmission and capture measurements performed on nat Ce and highly enriched 142 Ce samples at the JRC-Geel Linear Accelerator facility, as well as measured thermal constants available from the EXFOR database. Starting from the resonance parameters of the ENDF/B-VIII.0 library followed by a preliminary R-matrix analysis, an updated set of resonance parameters and corresponding covariance information were derived by fitting these measured data using the Reich–Moore approximation of the R -matrix theory, as implemented in the SAMMY code system. The 140 Ce RRR upper energy limit was kept at 200 keV, whereas the 142 Ce resonance region was extended from 13 to 26 keV. Updated statistical properties were obtained for the new evaluations and compared to those derived from the ENDF/B-VIII.0 nuclear data library. The new evaluation work improved some of the discrepancies found in previous work, such as the capture resonance integral and stellar Maxwellian-averaged cross sections. These integral quantities were mainly derived from the fit of the latest measured data, especially the neutron capture yield data for 142 Ce isotope.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

LDRD Abbreviated report: High-Order General-Discrete-Ordinates Method Enabling Efficient Deterministic Transport in Hydrodynamic Simulations

Deterministic transport simulations for national-security and energy applications often operate in high-dimensional phase-space, where accuracy and cost both become major challenges. A common numerical artifact in such problems is the “ray-effect,” which appears as unphysical streaks. Beyond misinterpretation, these artifacts can contaminate tightly coupled physics, such as fluid dynamics, radiation-hydrodynamics, and laser-plasma interactions, eroding the predictive capability of entire multiphysics workflows. Our objective was to make high-dimension studies practical on modern hardware while mitigating the ray-effect without relying on prohibitively expensive sampling approaches such as Monte Carlo methods. We developed the Generic Discretization Library (GenDiL), a Graphics Processing Unit (GPU)-first framework that uses high-order Discontinuous Galerkin (DG) methods and matrix-free algorithms to reduce memory usage and improve computational efficiency, critical for phase-space simulations. GenDiL supports phase-space adaptivity in both mesh size and polynomial order (hp-adaptivity) to place resolution only where it is needed. A central capability is Local Dimensional Refinement (LDR), which couples lower-dimension continuum models to higher-dimension kinetic models through stable and conservative interfaces, so that high-fidelity physics is applied only in regions where it is essential. Building on the GenDiL framework, we developed the General SN (GSN) family of algorithms as a true generalization of the polar SN approach (discrete ordinates, often denoted SN). Rather than tying discrete ordinates to a specific polar change of coordinates, GSN formulates transport on an arbitrary change of coordinates chosen to reduce ray-effect. We studied two complementary variants: an analytic variant, where the coordinate map is prescribed in advance by a closed-form function; and a data-driven variant, where a quantity of interest, such as the net flux, guides the coordinate system. GenDiL provides the library infrastructure for efficient GPU execution, but the GSN concept is algorithmic and independent of any one library. Across representative high-dimension tests, including non-symmetric solutions, both variants delivered strong ray-effect mitigation at practical cost, moving four- to six-dimensional analysis toward repeatable, routine studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗