Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational inelasticity”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Open-source library for performance-portable neutrino reaction rates: Application to neutron star mergers

A realistic and detailed description of neutrinos in binary neutron star (BNS) mergers is essential to build reliable models of such systems. To this end, we present bns_nurates, a novel open-source numerical library designed for the efficient on-the-fly computation of neutrino interactions, with particular focus on regimes relevant to BNS mergers. bns_nurates targets a higher level of accuracy and realism in the implementation of commonly employed reactions by accounting for relevant microphysics effects on the interactions, such as weak magnetism and mean field effects. It also includes the contributions of inelastic neutrino scattering off electrons and positrons and (inverse) nucleon decays. Finally, it offers a way to reconstruct the neutrino distribution function in the framework of moment-based transport schemes. As a first application, we compute both energy-dependent and energy-integrated neutrino emissivities and opacities for conditions extracted from a BNS merger simulation with m1 transport scheme. We find some qualitative differences in the results when considering the impact of the additional relevant reactions and of microphysics effects. For example, neutrino-electron/positron scattering reactions are important for the energy exchange of heavy-type neutrinos as they do not undergo semileptonic charged-current processes, when μ± are not accounted for. Moreover, weak magnetism and mean field effects can significantly modify the contribution of β processes for electron-type (anti)neutrinos, increasing at the same time the importance of (inverse) neutron decays. Here, the improved treatment for the reaction rates also modifies the conditions at which neutrinos decouple from matter in the system, potentially affecting their emission spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Accelerating resonant spectroscopy simulations using multishifted biconjugate gradient

Resonant spectroscopies, which involve intermediate states with finite lifetimes, provide important insights into collective excitations in quantum materials that are otherwise inaccessible. However, theoretical understanding in this area is often limited by the numerical challenges of solving Kramers-Heisenberg-type response functions for large-scale systems. To address this, we introduce a multishifted biconjugate gradient algorithm that exploits the shared structure of Krylov subspaces across spectra with varying incident energies, effectively reducing the computational complexity to that of linear spectroscopies. Both mathematical proofs and numerical benchmarks confirm that this algorithm substantially accelerates spectral simulations, achieving constant complexity independent of the number of incident energies, while ensuring accuracy and stability. This development provides a scalable, versatile framework for simulating advanced spectroscopies in quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Development of a MOOSE-based crystal plasticity model with irradiation defect evolution for irradiation creep in 316

Irradiation creep and irradiation swelling are two of the lifetime limiting factors for structural materials in nuclear reactors. These mechanical effects are driven by irradiation defect evolution and the interaction of those defects with dislocations in the microstructure. We present here a coupled cluster dynamics – crystal plasticity approach to model irradiation swelling and creep behavior in 316 SS. The time-dependent evolution of irradiation defects is calculated with a cluster dynamics approach and passed to the crystal plasticity model to compute the dislocation evolution. We show the impact of the irradiation defect evolution on the stress state in the material, which drives inelastic deformation through dislocation motion. The inelastic deformation in the 316 SS is dependent on the dose rate, where the inelastic deformation driven by the early-stage irradiation defect evolution determines the mechanical behavior of the 316 SS.

316 Stainless Steels↗
DOE OSTI2026

Witnessing Spin-Orbital Entanglement Using Resonant Inelastic X-Ray Scattering

Entanglement plays a central role in quantum technologies, yet its characterization and control in materials remain challenging. Recent developments in spectrum-based entanglement witnesses have enabled new strategies for quantifying many-body entanglement in macroscopic materials. Here, in this work, we develop a protocol for detecting spin-orbital entanglement using experiment-accessible resonant inelastic x-ray scattering. Central to our approach is the construction of a Hermitian generator from measurable spectra, which allows us to compute the quantum Fisher information (QFI) available in spin-orbital systems. The resulting QFI provides upper bounds for 𝑘-producible states and thus serves as a robust witness of spin-orbital entanglement. To account for realistic experimental limitations, we further extend our framework to include relaxed QFI bounds applicable to measurements lacking full polarization resolution.

entanglement detection↗
DOE OSTI2023

Design of Monolithic Bi-Layer High-Z PAL-Si Hard X-ray CMOS Image Sensors for Quantum Efficiency Enhancement

This article experimentally investigates the inception of an innovative hard X-ray photon energy attenuation layer (PAL) to advance high-energy X-ray detection (20–50 keV). A bi-layer design with a thin film high-Z PAL on the top and Si image sensor on the bottom has previously demon-strated quantum yield enhancement via computational methods by the principle of photon energy down conversion (PEDC), where high-energy X-ray photon energies are attenuated via inelastic scattering down to ≤10 keV, which is suitable for efficient photoelectric absorption by Si. Quantum yield enhancement has been experimentally confirmed via a preliminary demonstration using PAL-integrated Si-based CMOS image sensors (Si CIS). Furthermore, substituting the high-Z PAL with a lower-Z material—Sn—and alternatively coupling it with a conventional scintillator ma-terial—Lutetium-yttrium oxyorthosilicate (LYSO)—have been compared to demonstrate the most prominent efficacy of monolithic integration of high-Z PAL on Si CIS to detect hard X-rays, paving the way for next-generation high-energy X-ray detection methods.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Neutron Radiography and Computed Tomography of Biological Systems at the Oak Ridge National Laboratory's High Flux Isotope Reactor

Neutrons have historically been used for a broad range of biological applications employing techniques such as small-angle neutron scattering, neutron spin echo, diffraction, and inelastic scattering. Unlike neutron scattering techniques that obtain information in reciprocal space, attenuation-based neutron imaging measures a signal in real space that is resolved on the order of tens of micrometers. The principle of neutron imaging follows the Beer-Lambert law and is based on the measurement of the bulk neutron attenuation through a sample. Greater attenuation is exhibited by some light elements (most notably, hydrogen), which are major components of biological samples. Contrast agents such as deuterium, gadolinium, or lithium compounds can be used to enhance contrast in a similar fashion as it is done in medical imaging, including techniques such as optical imaging, magnetic resonance imaging, X-ray, and positron emission tomography. For biological systems, neutron radiography and computed tomography have increasingly been used to investigate the complexity of the underground plant root network, its interaction with soils, and the dynamics of water flux in situ. Moreover, efforts to understand contrast details in animal samples, such as soft tissues and bones, have been explored. This manuscript focuses on the advances in neutron bioimaging such as sample preparation, instrumentation, data acquisition strategy, and data analysis using the High Flux Isotope Reactor CG-1D neutron imaging beamline. The aforementioned capabilities will be illustrated using a selection of examples in plant physiology (herbaceous plant/root/soil system) and biomedical applications (rat femur and mouse lung).

Bilheux, Hassina↗
DOE OSTI2026

Semi-inclusive deep-inelastic scattering on a polarized spin-1 target. I. Cross section and spin observables

We develop the theoretical framework for semi-inclusive deep-inelastic scattering on a polarized spin-1 target and apply it to scattering on the polarized deuteron with spectator nucleon tagging. In Part I (this article) we present the general form of the semi-inclusive cross section and polarization observables for the spin-1 target. A relativistically covariant formulation in terms of four-vectors and invariant polarization parameters is employed. The target polarization is described by a spin density matrix with vector and tensor polarization. The spin and azimuthal angle dependence of the semi-inclusive cross section is derived and parametrized in terms of invariant structure functions. To validate the result, the structure functions are expressed as photon-target helicity amplitudes with known symmetry properties. The expressions presented here are kinematic (no assumptions about particle production dynamics) and valid in all regions of the deep-inelastic final state (current and target fragmentation regions). In Part II (following article), we consider deep-inelastic scattering on the polarized deuteron with spectator nucleon tagging as a special case of target fragmentation. The semi-inclusive structure functions are computed by separating nuclear and hadronic structure, and the polarization observables are explored as functions of the tagged nucleon momentum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Duality and form factors in the thermally deformed two-dimensional tricritical Ising model

The thermal deformation of the critical point action of the 2D tricritical Ising model gives rise to an exact scattering theory with seven massive excitations based on the exceptional E_7 E 7 Lie algebra. The high and low temperature phases of this model are related by duality. This duality guarantees that the leading and sub-leading magnetisation operators, \sigma(x) σ ( x ) and \sigma'(x) σ ′ ( x ) , in either phase are accompanied by associated disorder operators, \mu(x) μ ( x ) and \mu'(x) μ ′ ( x ) . Working specifically in the high temperature phase, we write down the sets of bootstrap equations for these four operators. For \sigma(x) σ ( x ) and \sigma'(x) σ ′ ( x ) , the equations are identical in form and are parameterised by the values of the one-particle form factors of the two lightest \mathbb{Z}_2 ℤ 2 odd particles. Similarly, the equations for \mu(x) μ ( x ) and \mu'(x) μ ′ ( x ) have identical form and are parameterised by two elementary form factors. Using the clustering property, we show that these four sets of solutions are eventually not independent; instead, the parameters of the solutions for \sigma(x)/\sigma'(x) σ ( x ) / σ ′ ( x ) are fixed in terms of those for \mu(x)/\mu'(x) μ ( x ) / μ ′ ( x ) . We use the truncated conformal space approach to confirm numerically the derived expressions of the matrix elements as well as the validity of the \Delta Δ -sum rule as applied to the off-critical correlators. We employ the derived form factors of the order and disorder operators to compute the exact dynamical structure factors of the theory, a set of quantities with a rich spectroscopy which may be directly tested in future inelastic neutron or Raman scattering experiments.

Cortes Cubero, Axel↗
DOE OSTI2026

Altered morphology and diffusivity of water confined in MXenes: Machine learning–accelerated computations combined with experiments

Nanoconfined water exhibits unique properties compared to bulk water due to limited quantities, frustrated hydrogen bonding, and surface interactions, which are fundamental for energy storage and transport applications. We integrate machine learning–accelerated ab initio molecular dynamics with x-ray diffraction (XRD) and inelastic neutron scattering (INS) to systematically analyze the thermodynamic and dynamic behavior of water confined between functionalized (-F, -O, and -OH) two-dimensional (2D) Ti 3 C 2 T x MXene layers. As water intercalates between layers, the interlayer spacing exhibits layer-dependent staging characteristics. The water polarization can be flipped by the count and morphology of intercalated molecules interacting with MXene surface groups, resulting in varying electrostatic potential profiles. On the basis of interfacial electrostatic potential, hydrogen bond lifetime, and molecular orientation, we establish a linear combination of exponential model describing water diffusivity. These computational insights align well with experimental x-ray and neutron measurements, suggesting strategies for tuning water morphology and transport by tailoring MXene surface chemistry and water content for electrochemical energy storage and nanofluidic applications.

Tang, Jiawei [Southeast Univ., Nanjing (China)] (O↗
DOE OSTI2020

Scalar and tensor neutrino interactions

We constrain general Dirac neutrino interactions based on the Standard Model Effective Field Theory framework extended with right-handed neutrinos N (SMNEFT) using deep inelastic and coherent elastic neutrino scattering, nuclear beta decay, and meson decay data, and high energy electron-proton and proton-proton collider data. We compute the one-loop anomalous dimensions of the low-energy effective field theory (LEFT) below the electroweak scale and of SMNEFT above the electroweak scale. The tree-level matching between LEFT and SMNEFT is performed at the electroweak scale. Currently, the most stringent limits on scalar and tensor interactions arise from pseudoscalar meson decays and the LHC measurements at the per mille level. In the future, the upcoming High-Luminosity LHC (HL-LHC) has the potential to reach the 10 -4 level and LHeC can play an important role under certain theoretical assumptions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The mechanics of brittle granular materials with coevolving grain size and shape

The influence of particle shape on the mechanics of sand is widely recognized, especially in mineral processing and geomechanics. However, most existing continuum theories for engineering applications do not encompass the morphology of the grains and its evolution during comminution. Similarly, the relatively few engineering models accounting for grain-scale processes tend to idealize particles as spheres, with their diameters considered as the primary and sole geometric descriptor. This paper inspires a new generation of constitutive laws for crushable granular continua with arbitrary, yet evolving, particle morphology. We explore the idea of introducing multiple grain shape descriptors into Continuum Breakage Mechanics (CBM), a theory originally designed to track changes in particle size distributions during confined comminution. We incorporate the influence of these descriptors on the elastic strain energy potential and treat them as dissipative state variables. In analogy with the original CBM, and in light of evidence from extreme fragmentation in nature, the evolution of the additional shape descriptors is postulated to converge towards an attractor. Comparisons with laboratory experiments, discrete element analyses and particle-scale fracture models illustrate the encouraging performance of the theory. The theory provides insights into the feedback among particle shape, compressive yielding and inelastic deformation in crushable granular continua. These results inspire new questions that should guide future research into crushable granular systems using particle-scale imaging and computations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Two-particle correlation functions in cluster perturbation theory: Hubbard spin susceptibilities

Cluster perturbation theory (CPT) is a computationally economic method commonly used to estimate the momentum- and energy-resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. Finally, we benchmark the method with the one-dimensional Fermi-Hubbard model by comparing with known results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Nonlinear soliton confinement in weakly coupled antiferromagnetic spin chains

In this work, we analyze the low-energy dynamics of quasi-one-dimensional, large-$\textit{S}$ quantum antiferromagnets with easy-axis anisotropy, using a semiclassical nonlinear sigma model. The saddle point approximation leads to a sine-Gordon equation which supports soliton solutions. These correspond to the movement of spatially extended domain walls. Long-range magnetic order is a consequence of a weak interchain coupling. Below the ordering temperature, the coupling to nearby chains leads to an energy cost associated with the separation of two domain walls. From the kink-antikink two-soliton solution, we compute the effective confinement potential. At distances large compared to the size of the solitons the potential is linear, as expected for pointlike domain walls. At small distances the gradual annihilation of the solitons weakens the effective attraction and renders the potential quadratic. From numerically solving the effective one-dimensional Schrödinger equation with this nonlinear confinement potential we compute the soliton bound state spectrum. We apply the theory to CaFe 2 O 4 , an anisotropic $\textit{S}$ = 5/2 magnet based upon antiferromagnetic zigzag chains. Using inelastic neutron scattering, we are able to resolve seven discrete energy levels for spectra recorded slightly below the Néel temperature $T_N$ ≈ 200 K. These modes are well described by our nonlinear confinement model in the regime of large spatially extended solitons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Phonon Thermal Transport in UO 2 via Self-Consistent Perturbation Theory

Computing thermal transport from first-principles in UO 2 is complicated due to the challenges associated with Mott physics. Here, we use irreducible derivative approaches to compute the cubic and quartic phonon interactions in UO 2 from first principles, and we perform enhanced thermal transport computations by evaluating the phonon Green’s function via self-consistent diagrammatic perturbation theory. Our predicted phonon lifetimes at T = 600 K agree well with our inelastic neutron scattering measurements across the entire Brillouin zone, and our thermal conductivity predictions agree well with previous measurements. Finally, both the changes due to thermal expansion and self-consistent contributions are nontrivial at high temperatures, though the effects tend to cancel, and interband transitions yield a substantial contribution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct prediction of inelastic neutron scattering spectra from the crystal structure*

Abstract Inelastic neutron scattering (INS) is a powerful technique to study vibrational dynamics of materials with several unique advantages. However, analysis and interpretation of INS spectra often require advanced modeling that needs specialized computing resources and relevant expertise. This difficulty is compounded by the limited experimental resources available to perform INS measurements. In this work, we develop a machine-learning based predictive framework which is capable of directly predicting both one-dimensional INS spectra and two-dimensional INS spectra with additional momentum resolution. By integrating symmetry-aware neural networks with autoencoders, and using a large scale synthetic INS database, high-dimensional spectral data are compressed into a latent-space representation, and a high-quality spectra prediction is achieved by using only atomic coordinates as input. Our work offers an efficient approach to predict complex multi-dimensional neutron spectra directly from simple input; it allows for improved efficiency in using the limited INS measurement resources, and sheds light on building structure-property relationships in a variety of on-the-fly experimental data analysis scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Reactivity Coefficient Measurements to Aid in Reducing Compensating Errors in Plutonium Nuclear Data

Compensating errors between several nuclear data observables in a nuclear data library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors between fast (0.1–5 MeV) 239Pu nuclear data for prompt fission neutron spectra (PFNS), average prompt fission neutron multiplicities, and neutron induced fission, capture, elastic, and inelastic cross sections. This work will focus on the design and execution of void reactivity coefficient measurements in the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). Two different base configurations were designed and measured, one with high neutron leakage, and one with low neutron leakage. Both were primarily made up of plutonium metal (Zero Power Physics Reactor plates) without interstitial moderators and reflected by half-inch aluminum. Design optimization showed that void reactivity coefficient measurements in three locations per configuration was most impactful to reduce nuclear data uncertainties due to the varying impacts from elastic and inelastic scattering, as well as fission and capture. The locations for measurements were chosen based on preliminary studies which balanced measurement uncertainty and measurement practicality. The measurements were also selected to have sensitivities maximally complementary to previous arrangements. Comparisons across nuclear data libraries highlight the potential impact.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

One Dimensional Penstock Flow Models for Hydropower Digital Twin

This report details the inelastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The flow model is coupled to a mechanistic turbine model to compute the rotation rate of the turbine, which is dependent on the flow rate, gravitational head, and various parameters of the machine, such as inlet vane angle. Results are given for an example hydropower facility in which the turbine is subjected to a specified time-varying load, and the results on the flow rate are shown. Future applications for the flow model are also briefly reviewed.

13 HYDRO ENERGY↗