Search NASA⌕ Search

SEARCH · Search NASA

Results for “Models, Structural”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Structure of beta-decaying states in the deformed, neutron-rich nucleus 104 Nb

Excited structures in 104 Mo were populated by β decays of the ground and isomeric states in the neutron-rich nucleus 104 Nb. The beams were produced by the CARIBU facility at Argonne National Laboratory, re-accelerated by the ATLAS accelerator and implanted on a moving-tape system in the middle of the GAMMASPHERE array. Separate decay schemes for the two β-decaying states in 104 Nb were constructed for the first time. The structure of the isomers are discussed in the framework of the deformed Nilsson model and systematics of known quasiparticle structures in neighboring nuclei.

Nilsson model↗
DOE OSTI2024

Structure and disorder resulting in ultralow thermal conductivity in the defect chalcopyrite AgInSnSe 4

Quaternary chalcogenides are attractive for a variety of technological fields of interest due to the diverse physical properties these materials possess. Herein, the structural and thermal properties of the disordered chalcopyrite AgInSnSe 4 are reported revealing a distinct relation between its structure and physical properties. Modeling of the experimental temperature-dependent thermal properties revealed an ultralow thermal conductivity (0.47 Wm -1 K -1 at room temperature) due to lattice anharmonicity, a low speed of sound and a low Debye temperature. Moreover, first principles electronic structure calculations revealed that the Ag-Se tetrahedra within the crystal structure possess relatively weak bonds due to occupied antibonding states from p-d orbital hybridization, which suppress the thermal conductivity. Here, in order to further quantify our findings, we extended our analyses to include comparisons with other ternary and quaternary adamantine materials, all of which possess a fourfold tetrahedral coordination of atoms, in revealing the origin of the thermal properties in AgInSnSe 4 .

36 MATERIALS SCIENCE↗
DOE OSTI2026

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗
DOE OSTI2024

Numerical Simulation of Infrasound Resonance in Underground Tunnel Structures

Remote observation of infrasound resonant signals emanating from underground tunnel structures could potentially allow the remote quantification of the geometry of the underground structures in which the signals were generated. However, the sensitivity of these observations to tunnel geometry and changes in that geometry are unknown. In this report we outline a numerical simulation study with the following three objectives: 1. Can we model infrasound resonance using numerical simulations? 2. What is the sensitivity of the modeled observations to changes in tunnel geometry or boundary conditions? And 3. Can we accurately model resonant observations of explosions occurring in Redmond Salt Mine in Central Utah. In this report we outline affirmative answers to the first two objectives, but we were unable to accurately model the Redmond Explosions due to numerical instability in the model of the complex structure of the mine. Recommendations for future work emphasize the need to acquire additional datasets and to explore a more data-based approach in which changes in data signatures are detected as a first step towards developing a method which inverts resonant infrasound signals for tunnel geometry.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

A State-Space Model for Stability Boundary Analysis of Grid-Following Voltage Source Converters Considering Grid Conditions

With the growing significance of renewable energy resources and energy storage systems, the number of grid-connected inverters has been rising at an increasingly rapid pace. Generally, these inverters are directly integrated with the distribution network by synchronizing with the grid voltage at the point of common coupling. However, the low grid strength and varying R/X ratios, as the common characteristics of most distribution networks or weak grids, can lead to dynamic interactions that comprise stability and limit the power transfer capacity of grid-connected inverters. To ensure stable operation of the inverters, researchers must determine the stability boundary, described as the maximum power transfer capacity of grid-connected inverters under the premise of maintaining system small-signal stability. For this purpose, we propose to formulate a state-space model of the system in the synchronously rotating dq-frame of reference and perform eigenvalue analysis to determine the stability boundary. With a detailed model of the control structure and parameters of the grid-connected inverters, the stability boundary is identified as a surface with respect to different grid strengths and R/X ratios. Case study results of proposed eigenvalue analysis are compared with those of admittance model-based stability analysis as well as time-domain simulation using a switching model in Matlab/Simulink, validating the effectiveness and accuracy of the proposed eigenvalue analysis for stability boundary identification.

grid-connected inverters↗
DOE OSTI2024

Quantum stabilization and flat hydrogen-based bands of nitrogen-doped lutetium hydride

We explore the electronic and structural properties of Lu-H-N structures with specific N,H ordering due to their remarkably narrow hydrogen-based bands at the Fermi level. Although exhibits a classical instability persisting up to 17 GPa, it is anharmonically stable near ambient pressure when accounting for quantum nuclear effects. The presence of flat bands near is understood to arise from destructive quantum interference between and surrounding orbitals, with certain types of defects leaving the flat bands unaffected. The results suggest that if this system is superconducting, there is an optimal pressure near ambient where the superconducting is maximized by anharmonically stabilized low-frequency and nonadiabatically coupled high-frequency hydrogen modes. Despite the metastability of this structure, its dynamical stability when modeled beyond a classical harmonic approach suggests the narrow bands near do not cause structural distortions in certain Lu-H-N stoichiometries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Integration of the fundamental knowledge on solvent-packing interactions into the multiscale framework for column scale design and optimization

The interfacial area, also known as the effective mass transfer area, is a key factor for determining the mass transfer for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Most of the widely used empirical and semi-empirical models for interfacial area were derived indirectly through absorption mass transfer with simplifications based on fast chemical kinetics. This report presents the comprehensive unique multiscale approach to develop a surrogate model for effective mass transfer area in structured packed columns that accounts local hydrodynamics as well as variation in physical properties, and changes in solid surface characteristics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Unconventional low temperature decomposition of a saturated hydrocarbon over atomically-dispersed titanium-aluminum-boron catalyst

Sonochemically-synthesized atomically-dispersed titanium-aluminum-boron nanopowder (TiAlB NP) exhibits a remarkable low-temperature catalytic activation of aliphatic C-H bonds at 750 K followed by C-C bond activation thus emerging as a potent low-cost alternative to expensive platinum group metals. Here, the model saturated hydrocarbon, exo-tetrahydrodicyclopentadiene (C10H16), undergoes catalytic decomposition on TiAlB NPs in a chemical microreactor to produce 1,3-cyclopentadiene (c-C5H6), cyclopentene (c-C5H8), and molecular hydrogen (H2) as detected in situ via isomer-selective, single-photon ionization time-of-flight mass spectrometry. Extensive electronic structure theory calculations on model clusters of the catalyst decode a unique synergy among the atomic constituents of the catalyst and chemical bonding in this stepwise, retro Diels Alder reaction: Ti, although insensitive to C-H activation in its metallic state, initiates the catalysis via chemisorption of the hydrocarbon, adjacent B centers readily abstract hydrogen atoms and store them during the catalytic cycle, while Al stabilizes the catalyst structure yet providing space for critical docking sites for the departing hydrocarbons.

Biswas, Souvick↗
DOE OSTI2025

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2024

BEAST DB: Grand-Canonical Database of Electrocatalyst Properties

We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.

database↗
DOE OSTI2025

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Random Forest Prediction of Crystal Structure from Electron Diffraction Patterns

Transmission electron microscopy (TEM) diffraction patterns are regularly used to determine the structure of crystalline materials. Electron diffraction is the most common method to solve for unknown or partially known crystal structures, as it provides direct and interpretable feedback on the orientation of crystal grains under the beam [1]. However, it remains a challenge to determine the crystal structure of a new material or even a new phase of an existing material. Analysis of such materials commonly requires manual exploration and comparison with simulated diffraction patterns. This is often a time consuming process with no obvious start point when many similar structures are possible, and this method cannot be used to determine crystal structure or orientation from structures not included in the diffraction libraries. Therefore, we have developed a machine learning model to determine the crystal structure of a material from its electron diffraction pattern.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparing computational times for simulations when using PBPK model template and stand-alone implementations of PBPK models

Introduction We previously developed a PBPK model template that consists of a single model “superstructure” with equations and logic found in many physiologically based pharmacokinetic (PBPK) models. Using the template, one can implement PBPK models with different combinations of structures and features. Methods To identify factors that influence computational time required for PBPK model simulations, we conducted timing experiments using various implementations of PBPK models for dichloromethane and chloroform, including template and stand-alone implementations, and simulating four different exposure scenarios. For each experiment, we measured the required computational time and evaluated the impacts of including various model features (e.g., number of output variables calculated) and incorporating various design choices (e.g., different methods for estimating blood concentrations). Results We observed that model implementations that treat body weight and dependent quantities as constant (fixed) parameters can result in a 30% time savings compared with options that treat body weight and dependent quantities as time-varying. We also observed that decreasing the number of state variables by 36% in our PBPK model template led to a decrease of 20–35% in computational time. Other factors, such as the number of output variables, the method for implementing conditional statements, and the method for estimating blood concentrations, did not have large impacts on simulation time. In general, simulations with PBPK model template implementations of models required more time than simulations with stand-alone implementations, but the flexibility and (human) time savings in preparing and reviewing a model implemented using the PBPK model template may justify the increases in computational time requirements. Conclusion Our findings concerning how PBPK model design and implementation decisions impact computational speed can benefit anyone seeking to develop, improve, or apply a PBPK model, with or without the PBPK model template.

Bernstein, Amanda S.↗
DOE OSTI2025

Modeling neutral defects in III-V ternary alloys with a special quasirandom structure: Analysis of As- and III-site point defects in InGaAs

While first-principles density functional theory modeling has become a vital tool to investigate defect properties in semiconductors, the lack of crystalline periodicity in pseudobinary random composition alloys, such as In 1−𝑥 ⁢Ga 𝑥 ⁢As, complicates such analyses. We present a simulation strategy to systematically take into account the variability in the local defect environment in order to predict statistical properties of neutral intrinsic defects in In 1−𝑥⁢ Ga 𝑥 ⁢As. We use a comprehensive sampling from a modest-sized 64-atom special quasirandom structure (SQS) to define a statistically representative set of defects, and use a 512-atom hypercell, a 2 × 2 × 2 supercell of SQS supercells, to achieve cell-size convergence. We articulate an equivalent site principle and describe how it constrains atomic chemical reference energies in computation of defect formation energies in pseudobinary alloys. A simple protocol for estimating reference energies for the Ga and In atoms sharing the III site succeeds in obtaining the equivalence of defects at Ga-sites and In sites in the SQS supercell, (<30 meV differences in average formation energies). For III-site defects, such as the As antisite As III , the statistical variability in formation energies is modest, ≈ 0.1–0.2 eV. The variability in formation energy at As-site defects, such as the As vacancy 𝑣 As , can be much larger, >1 eV. The As antisite is shown to be a low-energy defect and the most likely to be present in as-grown materials, just as in GaAs. All other defects are higher-energy defects unlikely to be important in native material, but potentially important in radiation-damaged material. With a strong variability in defect energies, especially on the As-site, explicit consideration of statistical variability due to compositional randomness will be imperative for meaningful and quantitative comparisons to experiment.

Density functional theory↗
DOE OSTI2025

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗
DOE OSTI2024

A Perspective on Scalable AI on High-Performance Computing and Leadership Class Supercomputing Facilities [Industrial and Governmental Activities]

Many scientific applications that support the mission of the US Department of Energy (US-DoE) require modeling complex engineering and/or physical systems. Here, examples of such complex systems arise from: (a) materials science to develop new compounds with exceptional mechanical and thermodynamical properties (e.g., resistance to mechanical stresses and high temperatures), (b) structural and nuclear engineering to model the temporal evolution of the structural damage of concrete shields exposed to continuous neutron and gamma radiations emitted by the nuclear reactor core, (c) urban sciences (e.g., transportation and smart buildings), and (d) power grid systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Atomistic Insights into Lithium–Glyme Solvate Ionic Liquids: Effects of Chain Length and Anion Coordination

For this study, mixtures of lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) in diglyme (G2), triglyme (G3) and tetraglyme (G4) at solvate ionic liquid (SIL) concentrations were investigated using classical molecular dynamics (cMD) simulations with a physically motivated force-field specifically developed for modeling these systems. The structural and dynamical properties of the mixtures were computed and analyzed. Lithium solvation shells, radial distribution functions, and X-ray structure factors were studied across the different SIL systems. Translational diffusion and rotational relaxation times were also evaluated, exhibiting similar trends with increasing glyme chain length. The results are consistent with experimental data and in good agreement with previous computational studies on G3 and G4. These findings validate the accuracy of the force field in modeling glyme systems and its use for describing the [Li(G2) 4/3 ][TFSI] mixture. Additionally, the thermal and electrochemical stability of these electrolytes were systematically examined. The thermal stability appears to be governed by cooperative interactions among glyme molecules, while the electrochemical stability is primarily influenced by Li + -anion interactions, which vary significantly with glyme chain length. Overall, the study sheds light on the crucial role of the anion in these glyme-based SILs and offers valuable insights into Li + -glyme systems at SIL concentrations, highlighting their promise as potential Li-ion battery electrolytes.

anions↗