SEARCH · Search NASA

Results for “parameterization”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Multiscale study of helium diffusion in Ni-Cr alloys: Short-range trapping versus long-range channeling

Ni-Cr alloys are widely employed as structural materials in fast nuclear reactors but are vulnerable to high-temperature helium (He) embrittlement (HTHE) under fast neutron irradiation. A comprehensive understanding of He diffusion in Ni-Cr alloys, which governs the kinetics of HTHE, is therefore essential for developing resilient materials and preventing failure. In this work, we reveal the underlying mechanisms of He diffusion in pure Ni and Ni-Cr alloys by integrating density functional theory (DFT) with atomic kinetic Monte Carlo (AKMC) simulations. Our findings uncover a non-monotonic dependence of He diffusivity on Cr concentration, contradicting the monotonic trends predicted by DFT-parameterized theories. At low Cr concentrations, He diffusion is dominated by short-range trapping, characterized by multiple trapping sites and a distinct mechanism within the first nearest neighbor of Cr, differing from that in pure Ni. At high Cr concentrations, these local traps become interconnected, forming long-range fast diffusion channels that enhance He mobility. The competition between localized trapping and extended channeling results in a diffusivity that first decreases, then increases with rising Cr content. These atomic-scale insights offer critical guidance for the design of radiation-tolerant Ni-based alloys. Moreover, the combined DFT-AKMC methodology and the concept of random walker diffusion through interconnected energy basins present a broadly applicable framework for studying transport phenomena in disordered systems.

36 - MATERIALS SCIENCE↗

DOE OSTI2026

Chiral anomaly: From vacuum to Columbia plot

Here, we use a low-energy effective approach, the extended linear sigma model, to study realizations of the U(1) A anomaly with different operators, linear and quadratic in the ’t Hooft determinant. After discussing the parameterization in agreement with vacuum’s phenomenology, we investigate the influence of these different anomaly terms on the Columbia plot: the square of the ’t Hooft determinant favors a cross-over for small quark masses. Finally, we also discuss the extension of the ’t Hooft determinant to cases in which different mesonic multiplets interact with each other. Novel chiral anomalous interaction terms involving excited (pseudo)scalar states, pseudovector, and pseudotensor mesons are expressed via a mathematical extension of the determinant, denoted as a polydeterminant.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

A theoretical study of solid solution strengthening in the refractory medium entropy alloy Nb 45 Ta 25 Ti 15 Hf 15

The refractory medium-entropy alloy (RMEA) Nb 45 Ta 25 Ti 15 Hf 15 exhibits exceptional tensile ductility and fracture toughness at ambient temperature, but its engineering applications are limited by a lack of high temperature strength. Using a machine-learning interatomic potential (MLIP) with near-density functional theory (DFT) accuracy, we conducted molecular dynamics (MD) and statics simulations of the behavior of dislocations with both screw and edge characters. We also analyze experimentally measured yield strengths using the Rao-Suzuki model and the Maresca-Curtin model modified to include a temperature-dependent shear modulus and a bulk modulus-dependent misfit volume, thereby uncovering the mechanisms underlying the yielding of this RMEA. Compared with the published experimental yield strength, the models parameterized by the MLIP effectively reproduce the experimental results over a wide temperature range. The models and MD simulations indicate that yielding is governed by screw dislocations, with dipole dragging as the dominant mechanism. In MD simulations, we observed a potential softening mechanism not considered by the Rao-Suzuki screw model: slow migration of interstitial jogs along the dislocation core, which could lead to the annihilation of vacancy and interstitial jog pairs by their combination.

BCC complex concentrated alloys↗

DOE OSTI2025

Development and assessment of hierarchical multi-reward reinforcement learning based potential for silicene with state-of-the-art models

We develop a new interatomic force field for Silicene, a 2D material with a buckled hexagonal lattice structure with high polymorphism. We introduce new parameterizations of a Tersoff model using a hierarchical multi-reward reinforcement learning (RL) methodology coupled with a continuous Monte Carlo Tree Search optimization. Our model significantly outperforms existing methods by enhancing the accuracy of predictions for the structural and thermodynamic properties of seven silicene polymorphs-including structure, energy, equation of state, elasticity, and phonon dispersion-when compared to established models. We further make a comprehensive comparison of the various models in predicting the mechanical and thermal properties of silicene. We trace the origin of the improved performance to the description of the angular dependence in the bond-order term, suggesting that modifying the angular terms in short-range models is essential to capture the structural diversity in low dimensional systems.

2D materials↗

DOE OSTI2026

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗

DOE OSTI2025

SympGNNs: Symplectic Graph Neural Networks for identifying high-dimensional Hamiltonian systems and node classification

Existing neural network models to learn Hamiltonian systems, such as SympNets, although accurate in low-dimensions, struggle to learn the correct dynamics for high-dimensional many-body systems. Herein, we introduce Symplectic Graph Neural Networks (SympGNNs) that can effectively handle system identification in high-dimensional Hamiltonian systems, as well as node classification. SympGNNs combine symplectic maps with permutation equivariance, a property of graph neural networks. Specifically, we propose two variants of SympGNNs: (i) G-SympGNN and (ii) LA-SympGNN, arising from different parameterizations of the kinetic and potential energy. We demonstrate the capabilities of SympGNN on two physical examples: a 40-particle coupled Harmonic oscillator, and a 2000-particle molecular dynamics simulation in a two-dimensional Lennard-Jones potential. Furthermore, we demonstrate the performance of SympGNN in the node classification task, achieving accuracy comparable to the state-of-the-art. Finally, we also empirically show that SympGNN can overcome the oversmoothing and heterophily problems, two key challenges in the field of graph neural networks.

Deep learning↗

DOE OSTI2025

Direct local parametrization of nuclear state densities using the back-shifted Bethe formula

Level densities are often parametrized using the back-shifted Bethe formula (BBF) for nuclei that possess experimental data for s-wave neutron resonance average spacings and a complete discrete level sequence at low excitation energies. However, these parametrizations require the additional modeling of the dependence of the spin-cutoff parameter on excitation energy. Here, in this work, we avoid the need to model the spin distribution of level densities by using the experimental data to parametrize directly the state densities, for which the BBF does not depend on the spin-cutoff parameter. This approach allows for a local parameterization of state densities that is independent of the spin-cutoff parameter. We provide these parameters in a tabulated form for applications in nuclear reaction calculations and for testing microscopic approaches to state densities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗

DOE OSTI2024

Modeling inclusive electron-nucleus scattering with Bayesian artificial neural networks

We introduce a Bayesian protocol based on artificial neural networks that is suitable for modeling inclusive electron-nucleus scattering on a variety of nuclear targets with quantified uncertainties. Unlike previous applications in the field, which directly parameterize the cross sections, our approach employs artificial neural networks to represent the longitudinal and transverse response functions. In contrast to cross sections, which depend on the incoming energy, scattering angle, and energy transfer, the response functions are determined solely by the energy and momentum transfer to the system, allowing the angular component to be treated analytically. We assess the accuracy and predictive power of our framework against the extensive data in the quasielastic inclusive electron-scattering database. Additionally, we present novel extractions of the longitudinal and transverse response functions and compare them with previous experimental analysis and nuclear ab-initio calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Improved nuclear-structure corrections to the hyperfine splitting of electronic and muonic deuterium

We calculate the nuclear-structure correction to the hyperfine splitting in both electronic and muonic deuterium using interactions from chiral effective field theory. We explore the sensitivity to different parameterizations of the nucleon-nucleon force, study the convergence pattern in the order-by-order chiral expansion, and estimate remaining uncertainties. Our results are consistent with earlier calculations from pionless effective field theory, offering new insights for a robust uncertainty quantification. Thanks to the order-of-magnitude reduction in uncertainty achieved with chiral effective field theory, the two-photon exchange contribution in electronic deuterium agrees with experimental extractions within 0.5σ, in contrast to the 2.6σ discrepancy observed in muonic deuterium. This study lays the groundwork for extending TPE calculations to HFS in heavier atomic systems.

Chiral effective field theory↗

DOE OSTI2026

Benchmarking state-of-the-art theory and empirical models of pionless neutrino-argon scattering in GENIE

Upcoming experiments need improved simulations of neutrino scattering. This work uses the popular GENIE event generator to test a variety of neutrino interaction models against recent MicroBooNE measurements of pionless charged-current interactions. The GENIE code can easily interchange model components, including nucleon form factor parameterizations, quasielastic cross-section calculations, treatments of the nuclear ground state and hadronic final-state interactions. Leveraging this software capability in comparisons with MicroBooNE data, the performance of some of GENIE’s most theoretically sophisticated model components is evaluated and contrasted with more empirically-driven alternatives.

Neutrino cross section↗

DOE OSTI2025

Assessment of flamelet/progress variable methods for supersonic combustion

Tabulated chemistry models, including the flamelet/progress variable approach, have been successfully used for a variety of turbulent flame simulations. The progress variable describes the progress of reactions in a system and parameterizes a lookup table of thermochemical variables. This approach reduces the cost of simulations, transporting only one scalar (progress variable) instead of the many species mass fractions required for detailed chemistry. Originally developed for low Mach number flame simulations, recent works have focused on extensions of this approach to compressible flames, supersonic combustion, and detonations, with applications such as scramjet combustors and rotating detonation engines. Unlike low Mach simulations, compressible flow simulations require solving the energy transport equation, which is coupled to the equation of state. This leads to additional modeling challenges regarding the thermodynamics and its impact on the chemistry. The validity of modeling assumptions, for example the relationship between energy and temperature, also varies with the combustion regime. The present work provides a detailed assessment of the existing strategies for chemistry tabulation for compressible/supersonic combustion, including detonations. A priori analysis indicates that approximations which are reasonable for weakly compressible flames may break down for shock-induced combustion. Furthermore, the analysis identifies specific assumptions and approximations that do not hold for detonations, emphasizing that care must be taken when applying tabulated chemistry models outside their intended combustion regimes.

Detonations↗

DOE OSTI2025

Characterizing leaf-scale fluorescence with spectral invariants

Sun-induced chlorophyll fluorescence (SIF) is increasingly recognized as a non-destructive probe for tracking terrestrial photosynthesis. Emerging developments in spectral invariants theory provide an innovative and efficient approach for representing SIF radiative transfer processes at the canopy scale. However, modeling leaf-scale fluorescence based on the spectral invariants properties (SIP) remains underexplored. In this study, the spectral invariants theory is employed for the first time to model the leaf-scale total, backward and forward fluorescence (leaf-SIP SIF). The leaf-SIP SIF model separates the leaf-scale radiative transfer process into two distinct components: the wavelength-dependent one associated with leaf biochemical properties, and the wavelength-independent component linked to leaf structural characteristics. The leaf structure-related effects are characterized by two spectrally invariant parameters: the photon recollision probability (p) and the scattering asymmetry parameter (q), which are parameterized using the directly measurable leaf dry matter. Evaluation against field measurements shows that the proposed leaf-SIP SIF model has a good performance, with coefficient of determination (R 2 ) of 0.89, 0.89, 0.90 and root mean squared errors (RMSE) of 1.28, 0.69, 0.74 Wm -2 µm -1 sr -1 , respectively for the total, backward, and forward fluorescence (660–800 nm). The leaf-SIP SIF model with a more concise formulation demonstrates comparable performance with the widely used Fluspect model. Further, the leaf-SIP SIF model provides a simple and efficient approach for simulating leaf-scale fluorescence, with the potential to be integrated into a unified SIP-based model framework for simulating the radiative transfer processes across the soil-leaf-canopy-atmosphere continuum.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗

DOE OSTI2026

Uncertain quantum computing futures and potential energy and physical resource impacts at scale

Considerable attention has recently focused on the vast energy and water demands of supercomputing, namely large-scale data centers that underpin artificial intelligence (AI), one of the great disruptors of contemporary society. Looking ahead some years from now, quantum computing is poised to disrupt established computing paradigms once again. Scientists and engineers are now working intensely to bring this century-old dream of physicists to fruition. Yet, as quantum computers begin to be integrated with classical supercomputing architectures, the implications for energy and physical resource use also need to be understood, especially how they compare to today’s AI data centers. These impacts have not yet been quantified by the research community – a notable gap in the literature, even if commercial-scale deployment of Quantum-Accelerated Computing Infrastructure (QuACI) is not expected for a few more years. This study is the first to conduct such an assessment. Using publicly available information from academic sources and private industry, we characterize multiple configurations of superconducting qubit-based, fault-tolerant quantum computers (FTQC) that could plausibly be deployed at scale in the 2030s and into the 2040s. By parameterizing these FTQC systems at a process level, we conduct a prospective scenario analysis to quantify their energy and physical resource needs. While these estimates are uncertain, given the current state of quantum technologies and their unknown future trajectories, important insights can already be drawn. One key finding is that while the electricity needs for a fleet of FTQCs are within the bounds of previous modeling studies that have explored high electricity demand futures, the needs for certain physical resources, namely water and helium-3, could pose bottlenecks to QuACI scale-up.

Computing↗

DOE OSTI2024

A simple model for short-range ordering kinetics in multi-principal element alloys

Short-range ordering (SRO) in multi-principal element alloys influences material properties such as strength and corrosion. While some degree of SRO is expected at equilibrium, predicting the kinetics of its formation is challenging. We present a simplified isothermal concentration-wave (CW) model to estimate an effective relaxation time of SRO formation. Estimates from the CW model agree to within a factor of five with relaxation times obtained from kinetic Monte Carlo (kMC) simulations when above the highest ordering instability temperature. Further, the advantage of the CW model is that it only requires mobility and thermodynamic parameters, which are readily obtained from alloy mobility databases and Metropolis Monte Carlo simulations, respectively. The simple parameterization of the CW model and its analytical nature makes it an attractive tool for the design of processing conditions to promote or suppress SRO in multicomponent alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Building a predictive model for polycyclic aromatic hydrocarbon dosimetry in organotypically cultured human bronchial epithelial cells using benzo[ a ]pyrene

The airway epithelium is a primary route of exposure for inhaled toxicants, and organotypic culture models represent an important advancement for toxicity testing compared to simple in vitro models that may lack metabolic capability and multicellular structure/communication associated with the bronchial epithelium in vivo. A quantitative understanding of chemical dosimetry is key for interpreting and extrapolating study results; however, dosimetry is understudied in organotypic models limiting ability to predict toxicity. We developed a dosimetry model for primary human bronchial epithelial cells (HBECs) cultured at the air-liquid interface (ALI) using benzo[a]pyrene (BAP), a representative polycyclic aromatic hydrocarbon. Dose and time course evaluation of metabolite formation and enzyme activity and expression were utilized to parameterize a cellular dosimetry model to improve the utility of ALI-HBECs for assessing chemical risk. Dosimetry analysis demonstrated absorption of BAP into cells and an increase in Phase 1 and 2 metabolites over time that correlated with regulation of metabolizing enzymes. BAP was cleared from cells by 48 h after exposure, and the primary metabolites generated in ALI-HBECs were BAP-3-phenol, BAP-4,5-dihydrodiol, BAP-7,8-dihydrodiol, BAP-9,10-dihydrodiol, BAP-7,8,9,10-tetrol, BAP-3-phenol-glucuronide, BAP-4,5-dihydrodiol-glucuronide, and BAP-9,10-dihydrodiol-glucuronide. The resulting dosimetry model described BAP and 7,8-dihydrodiol toxicokinetics in ALI-HBECs and suggested active excretion of 7,8-dihydrodiol. Overall, this study demonstrates metabolic competency of ALI-HBECs for BAP metabolism, demonstrates the usefulness of complex in vitro systems for human-relevant toxicity data, and exhibits how in silico models can be utilized for understanding the dosimetry of test compounds to aid in in vitro to human extrapolation of toxicity data for risk assessments.

Benzo[a]pyrene↗