Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large-scale transport simulation”

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,256 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Full-core high-burnup BWR LOCA fuel performance analysis and FFRD susceptibility

The susceptibility of the boiling water reactor (BWR) Limerick Unit 1 to fuel fragmentation, relocation, and dispersal during a postulated large-break loss-of-coolant accident (LBLOCA) was calculated using a multiphysics framework. The simulations include full-core, rod-resolved neutronic, thermal hydraulic, and fuel performance models using the VERA, TRACE, and BISON codes. This work focused on the transient BISON simulations, which include both the normal operation and LBLOCA periods in the same simulations. Cladding integrity was assessed using two correlations that are included with BISON. make page break Several new BWR-specific features were recently added to BISON. This work represents the first time these features have been included in a core-scale set of simulations. This study hence evaluates the performance of these new models for an operating reactor with realistic operating conditions. Simulation results showed that cladding integrity was maintained (i.e., no rods burst). Finally, future work to improve BWR and PWR predictions using this framework is suggested.

BISON↗
DOE OSTI2025

Thermalization and criticality on an analogue–digital quantum simulator

Abstract Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators 1,2 . Unlocking the full potential of such systems towards this goal requires flexible initial state preparation, precise time evolution and extensive probes for final state characterization. Here we present a quantum simulator comprising 69 superconducting qubits that supports both universal quantum gates and high-fidelity analogue evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. This hybrid platform features more versatile measurement capabilities compared with analogue-only simulators, which we leverage here to reveal a coarsening-induced breakdown of Kibble–Zurek scaling predictions 3 in theXYmodel, as well as signatures of the classical Kosterlitz–Thouless phase transition 4 . Moreover, the digital gates enable precise energy control, allowing us to study the effects of the eigenstate thermalization hypothesis 5–7 in targeted parts of the eigenspectrum. We also demonstrate digital preparation of pairwise-entangled dimer states, and image the transport of energy and vorticity during subsequent thermalization in analogue evolution. These results establish the efficacy of superconducting analogue–digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Science & Technology - Other Topics↗
DOE OSTI2026

Revealing the Hidden Third Dimension of Point Defects in Two-Dimensional MXenes

Point defects govern many important functional properties of two-dimensional (2D) materials. However, resolving the three-dimensional (3D) arrangement of these defects in multi-layer 2D materials remains a fundamental challenge, hindering rational defect engineering. Here, we overcome this limitation using an artificial intelligence-guided electron microscopy workflow to map the 3D topology and clustering of atomic vacancies in Ti3C2TX MXene. Our approach reconstructs the 3D coordinates of vacancies across hundreds of thousands of lattice sites, generating robust statistical insight into their distribution that can be correlated with specific synthesis pathways. This large-scale data enables us to classify a hierarchy of defect structures-from isolated vacancies to nanopores-revealing their preferred formation and interaction mechanisms, as corroborated by molecular dynamics simulations. This work provides a generalizable framework for understanding and ultimately controlling point defects across large volumes, paving the way for the rational design of defect-engineered functional 2D materials.

2D materials↗
DOE OSTI2026

Unveiling field-transverse spin anisotropy in the spin liquid candidate α-RuCl3 via spin Hall magnetoresistance

α-RuCl3 has emerged as a possible candidate for a quantum spin liquid (QSL) that promises exotic quasiparticles relevant for fault-tolerant quantum computation. Here, we report spin-sensitive transport measurements using a proximal spin Hall metal, platinum (Pt), to probe magnetic moments in the insulator α-RuCl3. We observe spin Hall magnetoresistance (SMR), where both the transverse and longitudinal resistivities exhibit angular oscillations between the in-plane magnetic field and the current, driven by the interplay between the spin Hall effect in Pt and local magnetic moments in α-RuCl3. These oscillations occur from 1.5 to 18 T, covering the zigzag antiferromagnetic, putative QSL, and supposedly partially field-polarized phases. The phase of the SMR oscillations suggests that the local moments, whether static or fluctuating, develop spin anisotropy with a quantization axis in-plane and largely transverse to the magnetic field across all fields from 1.5 to 18 T. Temperature dependence indicates that the spin anisotropy operates at an energy scale similar to that of the reported QSL signatures in α-RuCl3.

Idzuchi, Hiroshi [Tohoku University, Japan]↗
DOE OSTI2026

Predicting Atomistic Transitions with Transformers

Accurate knowledge of the atomistic transition pathways in materials and material surfaces is crucial for many material science problems. However, conventional simulation techniques used to find these transitions are extremely computationally intensive. Even with large-scale, accelerated material simulations, the computational cost constrains the applicable domain in practice. Machine learning models, with the potential to learn the complex emergent behaviors governing atomistic transitions as a fast surrogate model, have great promise to predict transitions with a vastly reduced computational cost. Here, we demonstrate how transformers can be trained to predict atomistic transitions in nano-clusters. We show how we evaluate physical validity of the predictions and how a multitude of additional, different microstates can be generated by slightly varying the data provided to the model.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Symmetry‐Breaking Strategy Yields Dopant‐Free Small Molecule Hole Transport Materials for Inorganic Perovskite Solar Cells with 20.58% Efficiency and Outstanding Stability

Abstract Inorganic perovskites are known for their excellent photothermal stability; however, the photothermal stability of all‐inorganic n‐i‐p perovskite solar cells (PSCs) is compromised due to ion diffusion and free radical‐induced degradation caused by the use of doped spiro‐OMeTAD hole transport materials (HTMs). In this study, two isomeric donor–acceptor–donor (D–A–D) type small molecules, namely HBT and HiBT, were developed and used as dopant‐free HTMs, using 2,1,3‐benzothiadiazole or benzo[d][1,2,3]thiadiazole as acceptor moieties. The HiBT molecule, with its symmetry‐breaking features, exhibits a large dipole moment, enhanced coordination‐active sites, and a well‐aligned energy level structure, all of which contribute to passivating perovskite surface defects and improving free charge separation. As a result, inorganic CsPbI 3 PSCs with HiBT HTM achieved an impressive power conversion efficiency (PCE) of 20.58%, the highest reported for dopant‐free HTM‐based inorganic PSCs. Moreover, the enhanced hydrophobic properties of HiBT molecules, coupled with their ability to passivate perovskite surface defects, contribute to significantly improved device stability. The unencapsulated devices based on HiBT HTM retained over 83% and 80% of their initial efficiency after being stored at 85 °C for 50 days and undergoing maximum power point (MPP) tracking at 85 °C for 1100 h, respectively. These results highlight that the symmetry‐breaking strategy is an exceptionally effective approach for designing efficient, dopant‐free small molecule HTMs, significantly contributing to both the high efficiency and enhanced stability of all‐inorganic PSCs.

Chemistry↗
DOE OSTI2025

Exploiting ν-dependence of projected energy correlators in HICs

We extend the recently derived factorization formula for energy-energy correlators to study the analytic structure of general ν-point projected energy correlators in heavy ion collisions. The ν-point projected energy correlators (or, ν-correlators) are an analytically continued family of the integer N-point projected energy correlators, which probe correlations between N final-state particles. By tracking the largest separation (χ) between the N particles, in vacuum, their structure is closely related to the DGLAP splitting functions and exhibits a classical scaling behavior ∼ 1/χ which is modified by resummation through the anomalous dimensions. We show that, in a thermal medium, the ν-correlators display non-trivial angular scaling already at the leading order in perturbation theory. We find that for non-integer values, particularly ν < 1, medium-induced jet function is enhanced compared to ν > 1. This is particularly manifested in the ratios of ν-correlators with respect to the two-point energy correlator, which encodes intrinsic angular information for ν < 1 when compared to large ν values. Moreover, for small-ν values, the ν-correlators appear to saturate at ν = 0.01 . We further confirm our leading-order numerical computations against simulated events from JEWEL for the parton-level production cross-section. Finally, we qualitatively discuss the effect of BFKL resummation for various values of ν.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 2

Topics covered include: Ringwoodite-olivine assemblages in Dhofar L6 melt veins; Amorphization of forsterite grains due to high energy heavy ion irradiation: Implications for grain processing in ISM; Validation of AUTODYN in replicating large-scale planetary impact events; A network of geophysical observatories for mars; Modelling catastrophic floods on the surface of mars; Impact into coarse grained spheres; The diderot meteorite: The second chassignite; Galileo global color mosaics of Io; Ganymede's sulci on global and regional scales; and The cold traps near the south pole of the moon.

Source record↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2025

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2024

Learning Nonlinear Reduced Models from Data with Operator Inference

This review discusses Operator Inference, a nonintrusive reduced modeling approach that incorporates physical governing equations by defining a structured polynomial form for the reduced model, and then learns the corresponding reduced operators from simulated training data. The polynomial model form of Operator Inference is sufficiently expressive to cover a wide range of nonlinear dynamics found in fluid mechanics and other fields of science and engineering, while still providing efficient reduced model computations. The learning steps of Operator Inference are rooted in classical projection-based model reduction; thus, some of the rich theory of model reduction can be applied to models learned with Operator Inference. This connection to projection-based model reduction theory offers a pathway toward deriving error estimates and gaining insights to improve predictions. Furthermore, through formulations of Operator Inference that preserve Hamiltonian and other structures, important physical properties such as energy conservation can be guaranteed in the predictions of the reduced model beyond the training horizon. This review illustrates key computational steps of Operator Inference through a large-scale combustion example.

Mechanics↗
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

Electrode Erosion and Prefire Studies Towards Fusion Scale Pulsed Power

This study presents a comprehensive investigation of electrode erosion and discharge behavior in spark gap switches over long switching cycle lifetimes. Brass, copper–tungsten (CuW), and stainless steel electrodes are tested under controlled conditions to quantify material degradation, debris accumulation, and changes in breakdown voltage. High-resolution imaging and statistical analysis of spark channel locations and gap breakdown voltages reveal how surface evolution influences long-term performance and reliability. These results provide essential data for lifetime modeling and inform design strategies for pulsed power systems in emerging applications such as private sector fusion energy and large-scale facilities like Sandia’s Z Machine and proposed ZX upgrades, where high repetition reliability and predictable behavior are critical.

electrical breakdown↗
DOE OSTI2026

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗
NASA NTRS1990

Discharge Transient Coupling in Large Space Power Systems

Experiments have shown that plasma environments can induce discharges in solar arrays. These plasmas simulate the environments found in low earth orbits where current plans call for operation of very large power systems. The discharges could be large enough to couple into the power system and possibly disrupt operations. Here, the general concepts of the discharge mechanism and the techniques of coupling are discussed. Data from both ground and flight experiments are reviewed to obtain an expected basis for the interactions. These concepts were applied to the Space Station solar array and distribution system as an example of the large space power system. The effect of discharges was found to be a function of the discharge site. For most sites in the array discharges would not seriously impact performance. One location at the negative end of the array was identified as a position where discharges could couple to charge stored in system capacitors. This latter case could impact performance.

N John Stevens↗
DOE OSTI2024

Tailoring Ion Transport in Li 3‐3y Ho 1+y Cl 6‐x Br x via Transition‐Metal Free Structural Planes and Charge Carrier Distribution

Abstract Localized atomistic disorder in halide‐based solid electrolytes (SEs) can be leveraged to boost Li + mobility. In this study, Li + transport in structurally modified Li 3 HoCl 6 , via Br − introduction and Li + deficiency, is explored. The optimized Li 3‐3 y Ho 1+ y Cl 6‐ x Br x achieves an ionic conductivity of 3.8 mS cm −1 at 25 °C, the highest reported for holmium halide materials. 6,7 Li nuclear magnetic resonance and relaxometry investigations unveil enhanced ion dynamics with bromination, attaining a Li + motional rate neighboring 116 MHz. X‐ray diffraction analyses reveal mixed‐anion‐induced phase transitions with disproportionate octahedral expansions and distortions, creating Ho‐free planes with favorable energetics for Li + migration. Bond valence site energy analysis highlights preferred Li + transport pathways, particularly in structural planes devoid of Ho 3+ blocking effects. Molecular dynamics simulations corroborate enhanced Li + diffusion with Br − introduction into Li 3 HoCl 6 . Li‐Ho electrostatic repulsions in the (001) plane presumably drive Li + diffusion into the Ho‐free (002) layer, enabling rapid intraplanar Li + motion and exchange between the 2d and 4h sites. Li 3‐3 y Ho 1+ y Cl 6‐ x Br x also demonstrates good battery cycling stability. These findings offer valuable insights into the intricate correlations between structure and ion transport and will help guide the design of high‐performance fast ion conductors for all‐solid‐state batteries.

Chemistry↗
DOE OSTI2025

Charm and bottom hadrons in hot hadronic matter

Heavy quarks, and the hadrons containing them, are excellent probes of the QCD medium formed in high-energy heavy-ion collisions, as they provide essential information on the transport properties of the medium and how quarks color-neutralize into hadrons. Large theoretical and phenomenological efforts have been dedicated thus far to assess the diffusion of charm and bottom quarks in the quark–gluon plasma and their subsequent hadronization into heavy-flavor (HF) hadrons. However, the fireball formed in heavy-ion collisions also features an extended hadronic phase, and therefore any quantitative analysis of experimental observables needs to account for the rescattering of charm and bottom hadrons. This is further reinforced by the presence of a QCD cross-over transition and the notion that the interaction strength is maximal in the vicinity of the pseudo-critical temperature. We review existing approaches for evaluating the interactions of open HF hadrons in a hadronic heat bath and the pertinent results for scattering amplitudes, spectral functions and transport coefficients. While most of the work to date has focused on D -mesons, we also discuss excited states as well as HF baryons and the bottom sector. Both the HF hadro-chemistry and bottom observables will play a key role in future experimental measurements. We also conduct a survey of transport calculations in heavy-ion collisions that have included effects of hadronic HF diffusion and assess its impact on various observables.

effective field theories↗