Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Atomistic Materials Modeling”

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.

1,810 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

NASA's Advanced Control Law Program for the F-8 Digital Fly-By-Wire Aircraft

This paper briefly describes the NASA F-8 Digital Fly-By-Wire (DFBW) and Langley Research Center's role in investigating and promoting advanced control laws for possible flight experimentation. The primary purpose of the paper is to provide background material for the papers that follow which describe some of the F-8 DFBW adaptive control law studies. The paper provides a brief description of the Phase II DFBW F-8 aircraft and its control system, discusses some of the advanced control law study objectives and guidelines, and provides some mathematical models which are useful in the control analysis problem.

Aerospace Materials↗
DOE OSTI2026

Disorder-induced spin-cluster magnetism in a doped kagome spin liquid candidate

The search for new quantum spin liquid materials relies on systems with strong frustration such as spins on an ideal kagome lattice. However, lattice imperfections can have substantial effects which are as yet not well understood. In recent work, the two-dimensional kagome system YCu 3 ⁢(OH) 6 ⁢[(Cl 𝑥 ⁢Br (1−𝑥) ) 3−𝑦 ⁢(OH) 𝑦 ] has emerged as a leading candidate hosting a Dirac spin liquid which appears to survive at least for 𝑥 < 0.4, associated with alternating-bond-hexagon (ABH) disorder. Here in magnetic samples with 𝑥 = 0.58, 𝑦 = 0.1 we report unusual in-plane ferromagnetic canting (FM) of the in-plane antiferromagnet (AFM), with an unusually wide regime of short-ranged order, and propose theoretical models to explain this behavior. First, we show that Kitaev-type exchanges naturally arise on the kagome lattice to second order in the known Dzyaloshinskii-Moriya exchanges, and that these interactions can produce the unusual in-plane FM canting from antichiral AFM. Second, we propose a phenomenological model of weakly FM-canted spin clusters to describe the short-ranged regime and analyze quantum fluctuations in an ABH toy model to show how ABH disorder can stabilize this regime. Here, the combination of experimental observation and theory suggests that kagome-Kitaev interactions and ABH disorder are necessary for describing the magnetic fluctuations in this family of materials, with potential implications for the proposed proximate spin liquid phase.

Seth, Arnab [Georgia Institute of Technology, Atla↗
DOE OSTI2026

Electronic glasses from a broken gauge symmetry in disorder-free systems

Glass phases can be stabilized by quenched disorders, as in most spin-glass materials, or self-generated through kinetic freezing in disorder-free systems. A canonical example of the latter is structural glasses, which have been extensively studied for many decades. Yet, how the rugged energy landscape of a glass phase is spontaneously generated in disorder-free systems remains one of the key questions in glass physics. Here, in this work, we present a general electronic mechanism for the emergence of glassy phase using the example of itinerant electrons coupled to XY spins on a lattice. This model can also be viewed as the mean-field theory of a superconducting system with attractive density-density interactions. Intriguingly, the electron gauge symmetry in the strong pairing limit gives rise to a macroscopic degeneracy of XY spins. In the presence of electron hopping that breaks the gauge symmetry, the lifting of the extensive degeneracy leads to a glass phase with disordered pairings. Our findings highlight a scenario in which a glassy state originates from the breaking of quantum gauge symmetry without quenched disorders.

XY model↗
DOE OSTI2026

High-Throughput Uniformity and Defect Monitoring in Low-Temperature Electrolysis Porous Transport Layers Using X-Ray Radiography

Effective quality control (QC) for manufacturing proton exchange membrane water electrolysis (PEMWE) components is critical to enabling widespread adoption of the technology for hydrogen generation. This study investigates X-ray radiography as a novel, high-throughput, potentially in-line QC technique for detecting defects and assessing material property distributions in titanium-based porous transport layers (PTLs) which constitute a crucial component of low temperature PEMWE stacks. We obtain radiographs of a set of fifteen PTLs and model their absorbance of the broadband radiation as a second-order polynomial to account for the non-monoenergetic radiation source used in this study. The resulting model serves as a basis for predicting the areal density and porosity distributions of the PTLs. We find radiography successful in detecting multiple instances of defects, including holes/depressions, cracks, and excess material on the surface or in the pores of the material, demonstrating its potential as a robust in-line QC tool for PTL manufacturing.

08 HYDROGEN↗
DOE OSTI2026

SENTRA: A Modular Computational Graph Framework for Critical Mineral and Materials Supply Chains: Part I: Network Construction Latent-Quantity Estimation, and Temporal Graph Forecasting

Global supply chains for critical minerals and materials are complex, evolving networks of countries, products, production stages, and trade relationships. Existing analytical approaches are limited by fragmented data and static network representations that do not capture the dynamic production dependencies linking raw materials, intermediate products, and final goods across multiple countries. Trade and production statistics provide only a partial view of domestic production, inventories, and material flows, making it difficult to identify indirect sourcing pathways, hidden dependencies, and embedded foreign exposures. This paper introduces the Supply Chain Exposure Network Tracking and Risk Assessment (SENTRA) framework, a modular graph-based computational framework for constructing, analyzing, and forecasting dynamic supply chain networks. As the first paper in a three-part methodological series, it establishes the computational foundation of SENTRA by constructing a temporal attributed multi-relational graph whose nodes represent product–country pairs and whose edges encode observed trade and within-country value-chain relationships. Statistical estimation and constrained optimization recover latent production, final demand, and product input dependency coefficients while enforcing economic accounting constraints. Graph-derived exposure measures quantify direct, transshipment, value-chain, and multi-hop supply chain dependencies independently of the forecasting model. A temporal graph forecasting architecture based on a relational graph neural network then forecasts the evolution of the graph under mass-balance constraints with distribution-free conformal uncertainty quantification. Validation on the global aluminum supply chain shows that the learned graph representations recover economically meaningful supply chain structure, accurately forecast out-of-sample trade relationships, and produce well-calibrated prediction intervals. Subsequent papers apply this computational foundation to exposure assessment, disruption analysis, and scenario-based policy analysis, and extend the framework to multimaterial supply chain modeling and decision support.

36 MATERIALS SCIENCE↗
NASA NTRS2026

Development and Validation of a High-Vacuum Thermal Conductivity Testbed for Aerospace Interface Materials

Thermal Interface Materials (TIMs) are critical components in spacecraft thermal management systems, where thermal performance is strongly influenced by vacuum conditions, interface contact resistance, and layered metallic joint behavior. However, manufacturer-reported thermal conductivity values are often derived under idealized conditions and may not accurately represent performance within operational aerospace applications. To address this limitation, the Testbed for Advanced Interface Materials in Vacuum (TAIMV) was developed as a modular vacuum-compatible thermal conductivity characterization platform capable of evaluating aerospace-relevant TIM configurations under both ambient and high-vacuum environments. The testbed was derived from the ASTM C1044-16 guarded hot plate methodology and incorporates interchangeable layers of stainless steel coupon geometries, independently controlled main and guard heaters, embedded resistance temperature detectors (RTDs), thermocouples, multi-layer insulation (MLI), and a temperature-controlled cold plate to characterize through-thickness thermal gradients across layered interfaces. In the current configuration, interface compression is limited to the nominal contact pressure generated by the experimental stack assembly. Initial experimental campaigns were conducted at ambient pressure and below 1×10-5 torr for vacuum cases using multiple interface materials including Braycote 601EF and Krytox-based greases across a range of thermal operating conditions. In parallel, a coupled numerical Python thermal model was developed to predict temperature distribution throughout the stack while accounting for conduction, radiation, and parasitic heat transfer pathways and effective interface resistance effects. Experimental measurements and numerical predictions showed consistent thermal trends across multiple operating conditions and environmental states. Results also revealed measurable differences between ambient and vacuum thermal behavior, demonstrating the importance of interface resistance, parasitic heat transfer mechanisms, and stack geometry in determining effective thermal performance within layered thermal interfaces. The presented work establishes a foundation for future thermal model correlation efforts and expanded characterization of aerospace thermal interface materials under representative environmental conditions. Future work will focus on the integration of a load cell system to enable controlled pressure-dependent characterization of thermal interface materials under compressive loading. This capability will allow investigation of the influence of contact pressure on effective thermal conductivity, interface resistance, and thermal performance within layered aerospace thermal interfaces under representative operational conditions.

Thermal Development Testing↗
DOE OSTI2024

Collective excitations and low-energy ionization signatures of relativistic particles in silicon detectors

Abstract Solid-state detectors with a low energy threshold have several applications, including searches of non-relativistic halo dark-matter particles with sub-GeV masses. When searching for relativistic, beyond-the-Standard-Model particles with enhanced cross sections for small energy transfers, a small detector with a low energy threshold may have better sensitivity than a larger detector with a higher energy threshold. In this paper, we calculate the low-energy ionization spectrum from high-velocity particles scattering in a dielectric material. We consider the full material response including the excitation of bulk plasmons. We generalize the energy-loss function to relativistic kinematics, and benchmark existing tools used for halo dark-matter scattering against electron energy-loss spectroscopy data. Compared to calculations commonly used in the literature, such as the Photo-Absorption-Ionization model or the free-electron model, including collective effects shifts the recoil ionization spectrum towards higher energies, typically peaking around 4–6 electron-hole pairs. We apply our results to the three benchmark examples: millicharged particles produced in a beam, neutrinos with a magnetic dipole moment produced in a reactor, and upscattered dark-matter particles. Our results show that the proper inclusion of collective effects typically enhances a detector’s sensitivity to these particles, since detector backgrounds, such as dark counts, peak at lower energies.

Physics↗
NASA NTRS1998

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff↗
NASA NTRS2004

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda↗
DOE OSTI2025

Revealing the evolution of order in materials microstructures using multi-modal computer vision

The development of high-performance materials for microelectronics, energy storage, and extreme environments depends on our ability to describe and direct property-defining microstructural order. Our present understanding is typically derived from laborious manual analysis of imaging and spectroscopy data, which is difficult to scale, challenging to reproduce, and lacks the ability to reveal latent associations needed for mechanistic models. Here, we demonstrate a multi-modal machine learning (ML) approach to describe order from electron microscopy analysis of the complex oxide La 1−x Sr x FeO 3 . We construct a hybrid pipeline based on fully and semi-supervised classification, allowing us to evaluate both the characteristics of each data modality and the value each modality adds to the ensemble. We observe distinct differences in the performance of uni- and multi-modal models, from which we draw general lessons in describing crystal order using computer vision.

36 MATERIALS SCIENCE↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2026

BLEECAM™ (Benchmarking Life Cycle Environmental, Economic, and Social Metrics for Critical and Advanced Minerals and Materials) [SWR-25-125]

The National Laboratory of the Rockies' (NLR) Benchmarking Life Cycle Environmental, Economic, and Social Metrics for Critical and Advanced Minerals and Materials (BLEECAM™) is an open-source, integrated decision-support tool for evaluating the impacts, risks, and trade-offs across U.S. and global materials supply chains. Funded by the U.S. Department of Energy, BLEECAM supports supply chain and market analysis. The tool integrates multi-objective supply chain optimization, system dynamics, network design, lifecycle assessment, techno-economic modeling, and social impact assessment methods to evaluate how supply chains evolve over time, geography, and deployment scenarios. BLEECAM also supports analysis related to energy infrastructure, data centers and digital infrastructure, advanced manufacturing, and other sectors that depend on critical materials.

Khalifa, SherifA. [National Laboratory of the Rock↗
DOE OSTI2025

Electronic structure prediction of medium and high entropy alloys across composition space

We propose machine learning (ML) models to predict the electron density — the fundamental unknown of a material’s ground state — across the composition space of concentrated alloys. From this, other physical properties can be inferred, enabling accelerated exploration. A significant challenge is that the number of descriptors and sampled compositions required for accurate prediction grows rapidly with species. To address this, we employ Bayesian Active Learning (AL), which minimizes training data requirements by leveraging uncertainty quantification capabilities of Bayesian Neural Networks. Compared to the strategic tessellation of the composition space, Bayesian-AL reduces the number of training data points by a factor of 2.5 for ternary (SiGeSn) and 1.7 for quaternary (CrFeCoNi) systems. We also introduce easy-to-optimize, body-attached-frame descriptors, which respect physical symmetries while keeping descriptor-vector size nearly constant as alloy complexity increases. Our ML models demonstrate high accuracy and generalizability in predicting both electron density and energy across composition space.

materials science↗
DOE OSTI2026

The Maximal Entanglement Limit in Statistical and High-energy Physics

These lectures advocate the idea that quantum entanglement provides a unifying foundation for both statistical physics and high-energy interactions. I argue that, at sufficiently long times or high energies, most quantum systems approach a Maximal Entanglement Limit (MEL) in which phases of quantum states become unobservable, reduced density matrices acquire a thermal form, and probabilistic descriptions emerge without invoking ergodicity or classical randomness. Within this framework, the emergence of probabilistic parton model, thermalization in the break-up of confining strings and in high-energy collisions, and the universal small-x behavior of structure functions arise as direct consequences of entanglement and geometry of high-dimensional Hilbert space.

36 MATERIALS SCIENCE↗
DOE OSTI2026

X-ray tomography of damage dynamics in advanced materials using a laser wakefield accelerator

Additively manufactured (AM) metals offer the potential for customizable, cost-effective components, but qualification and certification are crucial. Key to this process is understanding pore dynamics under stress, typically analyzed using micro-computed tomography. This study introduces laboratory-scale “betatron” x-rays from laser wakefield acceleration as a high-throughput alternative for x-ray tomography of advanced materials, such as AM AlSi10Mg alloys. Coupled with 3D finite element modeling, this method provides detailed insights into stress-porosity interactions. The approach delivers high-resolution scans, revealing that pore shape and local triaxiality significantly influence fracture dynamics, supporting advanced material characterization. This work also demonstrates the potential and versatility of laser-betatron x-ray μCT for generating large datasets to accelerate our understanding of the stochastic, process-specific nature of pore formation in AM alloys.

Senthilkumaran, Vigneshvar↗
NASA NTRS1995

RTE: A Computer Code for Rocket Thermal Evaluation

The numerical model for a rocket thermal analysis code (RTE) is discussed. RTE is a comprehensive thermal analysis code for thermal analysis of regeneratively cooled rocket engines. The input to the code consists of the composition of fuel/oxidant mixture and flow rates, chamber pressure, coolant temperature and pressure. dimensions of the engine, materials and the number of nodes in different parts of the engine. The code allows for temperature variation in axial, radial and circumferential directions. By implementing an iterative scheme, it provides nodal temperature distribution, rates of heat transfer, hot gas and coolant thermal and transport properties. The fuel/oxidant mixture ratio can be varied along the thrust chamber. This feature allows the user to incorporate a non-equilibrium model or an energy release model for the hot-gas-side. The user has the option of bypassing the hot-gas-side calculations and directly inputting the gas-side fluxes. This feature is used to link RTE to a boundary layer module for the hot-gas-side heat flux calculations.

Mohammad H N Naraghi↗