Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical integration”

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

Expert‐in‐the‐loop design of integral nuclear data experiments

Abstract Nuclear data are fundamental inputs to radiation transport codes used for reactor design and criticality safety. The design of experiments to reduce nuclear data uncertainty has been a challenge for many years, but advances in the sensitivity calculations of radiation transport codes within the last two decades have made optimal experimental design possible. The design of integral nuclear experiments poses numerous challenges not emphasized in classical optimal design, in particular, constrained design spaces (in both a statistical and engineering sense), severely under‐determined systems, and optimality uncertainty. We present a design pipeline to optimize critical experiments that uses constrained Bayesian optimization within an iterative expert‐in‐the‐loop framework. We show a successfully completed experiment campaign designed with this framework that involved two critical configurations and multiple measurements that targeted compensating errors in 239 Pu nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Fast Fourier transform evaluation of the Fresnel integral for gravitational-wave lensing

Gravitational waves (GWs) exhibit wave-optics effects when their wavelength is comparable to the scale of the gravitational lens. This may occur in lensing from galactic subhalos in GWs emitted by binary black-hole mergers and is gaining interest as a novel probe of dark matter. Predictions for observables in these cases ultimately rely on evaluating a Fresnel integral that quantifies the effect of lensing on the amplitude of a GW at a given frequency. However, numerical evaluation of this Fresnel integral is tricky, and several algorithms and publicly available codes that implement it have been developed. Here, we show that the dependence of this integral on the lens position can be written as a two-dimensional Fourier transform. Modern FFT techniques then enable rapid evaluation at all-sky positions simultaneously for general lenses without symmetry. Vectorization of FFT routines allows for derivatives with respect to model parameters to be obtained with only incremental additional computational cost. If the lens is axisymmetric, further speedups can be achieved with recently developed techniques for nonuniform fast Hankel transforms. To demonstrate, we make available Fresnel Integral Optimization with Nonuniform Transforms (fiona), an efficient and accurate code that is significantly faster than current methods for dense source grids, reaching 2 orders of magnitude speedups for ∼10 6 GW-emitting points. As part of FIONA , we developed code that provides vectorized nonuniform fast Hankel transforms that may have other uses (e.g., calculation of cosmological two-point correlation functions) beyond those considered here.

dark matter↗
DOE OSTI2021

Feasibility of full-core pin resolved CFD simulations of small modular reactor with momentum sources

Complex flow structure interactions and heat transfer processes take place in nuclear reactor cores. Given the extreme pressure/temperature and radioactive conditions inside the core, numerical simulations offer an attractive and sometimes more feasible approach to study the related flow and heat transfer phenomena in addition to the experiments. Under the Exascale Computing Project, the full-core simulation of a small modular reactor (SMR) has been pursued coupling Computational Fluid Dynamics (CFD) and neutronics. A key aspect of the modeling of SMR fuel assemblies is the presence of spacer grids and the mixing promoted by mixing vanes or the equivalent. A reduced order methodology is adopted based on momentum sources to mimic the mixing of the vanes. The momentum sources have been carefully calibrated with detailed Large Eddy Simulations (LES) of spacer grids performed with Nek5000. Modeling the spacer grid and mixing vanes (SGMV) effect without body-fitted computational grid avoids the excessive costs in resolving the local geometric details, and thus supports the simulation to be scaled up to the full core. Besides the progress on momentum source modeling, this paper also features the first full-core pin resolved CFD simulation ever performed to the authors' knowledge. This represents a significant advancement in capability for the CFD of nuclear reactors, which will hopefully serve as an inspiration for further integrating high-fidelity numerical simulations in actual engineering designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Geologic Carbon Storage of Anthropogenic CO2 under the Colorado Plateau in Emery County, Utah

Geologic Carbon Storage (GCS) is a promising technology for storing large volumes of anthropogenic CO2 effectively and permanently. Numerical simulations are an integral part of site selection and characterization for any potential GCS site. As part of the DOE-funded CarbonSAFE Rocky Mountains Phase I project, a regional GCS analysis was undertaken to understand the efficacy of storing CO2 emissions from the power generation and heavy industry in central Utah’s favorable geology. In this study, the injection of CO2 for geologic storage was simulated in the Navajo Sandstone Formation in Emery County, Utah. Carbon dioxide was sourced from regional power generation stations and heavy industries throughout Utah, with an emphasis on emissions reduction at the Hunter Power Plant near Castle Dale, Utah. A simulation grid was extracted from the project’s geological model encompassing an area around Price, Huntington, and Castle Dale in central Utah. The Navajo Sandstone Member of the Glen Canyon Group was the target of CO2 injection with the overlying Carmel formation providing the primary seal. A suite of simulations was performed assessing the viability of this area for permanent CO2 storage. Results indicate that the area can not only store 46 million metric tons of anthropogenic CO2, meeting the project goals, but this area has the capacity to securely store at least 1.3 billion tons of CO2, suggesting the injection site and surrounding geology are suitable locations for commercial-scale GCS.

58 GEOSCIENCES↗
DOE OSTI2020

Efficient Fourier transforms for transverse momentum dependent distributions

Hadron production at low transverse momenta in semi-inclusive deep inelastic scattering can be described by transverse momentum dependent (TMD) factorization. This formalism has also been widely used to study the Drell-Yan process and back-to-back hadron pair production in $e^+e^-$ collisions. These processes are the main ones for extractions of TMD parton distribution functions and TMD fragmentation functions, which encode important information about nucleon structure and hadronization. One of the most widely used TMD factorization formalism in phenomenology formulates TMD observables in coordinate $b_\perp$-space, the conjugate space of the transverse momentum. The Fourier transform from $b_\perp$-space back into transverse momentum space is sufficiently complicated due to oscillatory integrands that it requires a careful and computationally intensive numerical treatment in order to avoid potentially large numerical errors. Within the TMD formalism, the azimuthal angular dependence is analytically integrated and the two-dimensional $b_\perp$ integration reduces to a one-dimensional integration over the magnitude $b_\perp$. In this paper we develop a fast numerical Hankel transform algorithm for such a $b_\perp$-integration that improves the numerical accuracy of TMD calculations in all standard processes. Libraries for this algorithm are implemented in Python 2.7 and 3, C++, as well as FORTRAN77. All packages are made available open source.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

304L stainless steel↗
DOE OSTI2020

Towards improved speed and accuracy of laser powder bed fusion simulations via representation of multiple time scales

Here, we report the result of investigations into alternative representations of the temporal evolution of the laser powder bed fusion (LPBF) process. In order to provide the reader with sufficient context for the following discussion, the modeling challenge and its underlying cause, as well as current attempts to provide a satisfactory solution, will be discussed. Next, the author will describe two numerical methods (multirate time integration and parallel-in-time) to better represent the temporal scales of the problem. Finally, the results of applying these methods to a two-dimensional version of the heat transfer portion of the LPBF problem will be presented, reporting on both convergence and performance behavior, which indicate the possibility of significant speedup of solution calculations (40–100 times), along with well-characterized effects on solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quasi-optical beam tracing module development for millimeter-wave high-wavenumber collective scattering on the NSTX-U and EAST tokamaks

A Python3-based beam tracing code utilizing Quasi-Optics has been developed to track both incident and receiving beams in high-k collective millimeter wave scattering systems within magnetic fusion plasmas. In contrast to existing ray tracing codes that solely consider refraction, this beam tracing code incorporates diffraction phenomena, providing a more comprehensive calculation. Here, this enhanced capability allows for a more accurate calculation of the scattering volume and spatial resolution in high-k collective scattering systems, crucial for evaluating system performance and facilitating data analysis. Unlike Geometrical Optics, Quasi-Optics employs the complex eikonal method, representing a Gaussian beam as a collection of coupled rays to accurately preserve diffraction characteristics. The developed code is intended for application in NSTX-Upgrade and EAST high-k beam tracing analyses, targeting frequencies of 693 GHz and 270 GHz, respectively. The high-k system's primary objective is the observation of electron-scale instabilities. Employing a symplectic integrator, the code ensures numerical accuracy, assessed through the conservation of the Hamiltonian. With its precision and efficiency, the code facilitates rapid inter-shot analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A machine-learning approach to measure 3D sample properties from 2D Transmission Electron Microscopy images

Transmission Electron Microscopy (TEM) is a powerful tool for the characterization of materials at the nanoscale; however, its inherent two-dimensional (2D) nature poses significant challenges to accurately measure three-dimensional (3D) properties. We introduce a supervised machine-learning model that predicts 3D structural information, such as sample thickness and curvature, from a series of conventional 2D TEM images. The model, a U-Net convolutional neural network, is trained on a large synthetic dataset generated from dynamical diffraction simulations that model TEM’s complex, nonlinear image formation, accounting for sample thickness and curvature. This physically realistic framework enables exploration of a broad parameter space impractical to sample experimentally. We demonstrate that the trained model has accurate predictions for experimental single-crystal silicon samples, achieving performance comparable to established measurement techniques. This work highlights the critical role of robust, simulation-based training in overcoming the limitations of real-world imaging artifacts and inconsistent sample geometries. By integrating machine learning with numerical simulations, we offer an efficient and scalable framework for quantitative TEM analysis, paving the way for more sophisticated 3D characterization of complex materials.

Dynamical diffraction↗
DOE OSTI2025

A review on the state of thermal hydraulics research on air ingress scenarios in High-Temperature Gas-cooled Reactors following a D-LOFC

With the expectation of near-immediate carbon neutrality, widespread implementation of proven High-Temperature Gas-cooled Reactors (HTGRs) embodies a viable solution pathway given their inherent, passive safety features and high thermal efficiency. This study provides an overview of the current state of research involving the thermal hydraulics associated with air ingress from a depressurized loss of forced cooling (D-LOFC) in HTGRs. Accurately characterizing and predicting the physical phenomena underlying air ingress is of paramount concern, as the integrity of the fuel and core graphite support structures are threatened by the presence of oxygen. Broadly speaking, the air ingress scenario can be delineated into three main stages: (1) Depressurization, (2) Density-Driven Flow, and (3) Natural Convection. In tandem with the underlying fundamental theory, this review collates and synthesizes the existing body of contemporary research concerning the air ingress scenario following a D-LOFC. As evinced by this review, our current understanding and predictive abilities have benefited from extensive research, predominantly concentrated on the rate of air ingestion into the core. Here, additional research is necessary to holistically capture the phenomenology of an air ingress scenario following a D-LOFC by considering an additional variable: the oxygen content of the ingressing air. The latter variable requires investigation into the complex interactions of the fully integrated system. Additionally, while numerical tools are evolving domestically through the Nuclear Energy Advanced Modeling and Simulation program, a sufficiently validated code remains absent.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Harnessing the Second-Order Metal−Insulator Transition for Neuromorphic Computing

Vanadium oxides are widely studied phase change materials for brain-inspired computing architectures. Systems like VO 2 and V 2 O 3 exhibit first-order metal−insulator transitions (MITs) with hysteresis and percolative switching, increasing stochasticity and device variability. Here, we focus on the less explored Magnéli phase V 4 O 7 , which undergoes a continuous, non-hysteretic, second-order MIT. This surprisingly enables highly reproducible volatile resistive switching in spiking-neuron-type devices. We synthesize V 4 O 7 films, characterize their structural and transport properties, and demonstrate voltage and current-driven threshold switching with electrothermal feedback. In a Pearson–Anson oscillator, V 4 O 7 devices produce stable, tunable spiking across 20–200 kHz, with consistent operation among multiple devices. We introduce a numerical analog leaky-integrate-and-fire (aLIF) model that captures waveform shapes and their dependence on load resistance, temperature, and voltage. Furthermore, these findings suggest that second-order MIT materials like V 4 O 7 are promising for deterministic, scalable spiking neuron arrays for neuromorphic computing.

V4O7↗
DOE OSTI2023

Evolution of neuronal cell classes and types in the vertebrate retina

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs1. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates. By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Scaling patch analysis of turbulent planar plume

Proper scaling in turbulent planar plumes is investigated in this study using a scaling patch approach. Based on the scaled boundary conditions, a proper velocity scale for the mean axial flow is the plume centerline velocity U ref = U ctr , and a proper temperature scale for the temperature excess is Θ ref = T ctr – T ∞ , where T ctr is the plume centerline temperature and T ∞ is the ambient fluid temperature. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity, mean momentum, and mean energy equations, respectively, the following is found: (1) a proper scale for the mean transverse flow is V ref = ( d δ / d x ) U ctr , where d δ / d x is the growth rate of the plume width. (2) A proper scale for the Reynolds shear stress is R vu , ref = U ctr V ref = ( d δ / d x ) U ctr 2 , a mix of the scales for the mean axial and transverse flows. (3) A proper scale for the turbulent heat flux is R v θ , ref = V ref Θ ctr , a mix of the scales for the mean transverse flow and mean temperature excess. The mean transverse flow thus plays a critical role in the scaling of turbulent planar plumes. Approximate functions are developed for the scaled mean transverse flow, Reynolds shear stress, and turbulent temperature flux, and are found to agree favorably with experimental and numerical simulation data. The integral analysis of the mean momentum equation yields a Richardson number Ri, which remains invariant in the axial direction. The Richardson number is defined as Ri = def g β Θ ctr δ t / ( U ctr V ref ) ≈ 1 / 2 , where g is the gravitational acceleration, β is the thermal expansion coefficient, and δ t is the plume half-width based on the mean temperature profile. This Richardson number arises directly from the scaling patch analysis of the mean momentum equation, including both the streamwise and transverse velocity scales.

42 ENGINEERING↗
DOE OSTI2023

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Meta-GGA performance in solids at almost GGA cost

A recent modification, r 2 SCAN, of the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation (GGA) exchange-correlation functional mostly eliminates numerical instabilities and attendant integration grid sensitivities exhibited by SCAN. Here, we show that the successful deorbitalization of SCAN to SCAN-L (SCAN with density Laplacian dependence) carries over directly to yield r 2 SCAN-L. A major benefit is that the high iteration counts that hindered the use of SCAN-L are eliminated in r 2 SCAN-L. It therefore is a computationally much faster meta-GGA than its orbital-dependent antecedent. In this work, validation data for molecular heats of formation, bond lengths, and vibration frequencies (G3/99X, T96-R, and T82-F test sets, respectively) and on lattice constants, and cohesive energies (for 55 solids) and bulk moduli (for 40 solids) are provided. In addition, we show that the overmagnetization of bcc Fe, hcp Co, and fcc Ni persists in r 2 SCAN but does not appear in r 2 SCAN-L. Distinct from SCAN, both r 2 SCAN and r 2 SCAN-L give the correct nonmagnetic ground state for bcc V.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Evidence for Dissociation and Ionization in Shock Compressed Nitrogen to 800 GPa

Triple bonding in the nitrogen molecule (N 2 ) is among the strongest chemical bonds with a dissociation enthalpy of 9.8 eV/molecule. Nitrogen is therefore an excellent test bed for theoretical and numerical methods aimed at understanding how bonding evolves under the influence of the extreme pressures and temperatures of the warm dense matter regime. Here, we report laser-driven shock experiments on fluid molecular nitrogen up to 800 GPa and 4.0 g/cm 3 . Line-imaging velocimetry measurements and impedance matching method with a quartz reference yield shock equation of state data of initially precompressed nitrogen. Comparison with numerical simulations using path integral Monte Carlo and density functional theory molecular dynamics reveals clear signatures of chemical dissociation and the onset of L-shell ionization. Furthermore, combining data along multiple shock Hugoniot curves starting from densities between 0.76 and 1.29 g/cm 3 , our study documents how pressure and density affect these changes in chemical bonding and provides benchmarks for future theoretical developments in this regime, with applications for planetary interior modeling, high energy density science, and inertial confinement fusion research.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗