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

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2026

Surrogate modeling of Monte Carlo radiation transport with convolutional neural networks for shielding optimization

Here, we present a machine learning (ML)-based surrogate model using convolutional neural networks (CNN) designed to emulate the attenuation of neutron fields as they pass through various shielding materials. This model can compute the outgoing neutron flux almost instantaneously and achieves reasonable accuracy compared to traditional Monte Carlo (MC)-based codes, which are computationally intensive. This emulator alleviates the complexity of neutron radiation transport through shielding materials by reducing the dimensionality and enables shielding optimization for a known radiation environment. This optimization process, which would have taken an unrealistic timeline due to several complex radiation transport simulations, can now be achieved in minutes, thus increasing computational capabilities in radiation shielding assessment. We demonstrate the applications of this emulator in computing effective dose rates and optimizing shielding solutions for a heavy-ion accelerator facility, such as the Facility for Rare Isotope Beams, where secondary neutrons produced via beam interactions dominate the radiation environment.

accelerator shielding↗
DOE OSTI2025

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modulation of thermal conductivity of iron-doped ß-Ga2O3 by helium-ion irradiation

This study examines the impact of helium-ion irradiation on the thermal conductivity of ß-Ga2O3. A laser-based spatial domain thermoreflectance technique is used to investigate thermal conductivity map for both un-irradiated and irradiated ß-Ga2O3, which are then validated against simulation results derived from density functional theory-based phonon transport simulations. Since helium bubble evolution was ob- served at the nanoscale using transmission electron microscopy, the simulation study was carried out on eight distinct helium-induced sites in ß-Ga2O3. Our findings indicate a reduction in thermal conductivity for the irradiated samples. Experimental results show a significant reduction in thermal conductivity in irradiated samples, with de- creases of approximately 25% along the [100] direction and 40% along [001] directions. Phonon transport simulations closely replicate these findings, particularly when helium occupying interstitial sites, predicting reductions of ˜53% along [100] and ˜50% along [001] directions. This work underscores the role of irradiation-induced microstructural changes in the heat transport properties of ß-Ga2O3 which is crucial for its application in sensor devices in extreme environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2025

Development of a high fidelity CFD model for solvent evaporation and transport in porous structure during battery electrode drying

An efficient battery manufacturing process is the key to the mass production of Electric Vehicles (EV), in which drying is one of the most energy-intensive steps significantly influencing the battery cell performance. An accurate 3D CFD model for drying is essential for predicting the drying mechanism and optimizing its parameters. By optimizing the drying process, it is possible to reduce energy consumption and cost during battery manufacturing, minimize binder loading and maximize active material loading to achieve superior electrochemical performances and facilitate wider and faster public adoption of EV. This project aims to optimize the drying process during electrode manufacturing by leveraging high-fidelity, porous electrode simulations for solvent evaporation. By optimizing this process, we seek to reduce energy consumption during battery manufacturing, while minimizing binder loading and maximizing active material loading, with the overall goal of enhancing electrical vehicle performance.

Horner, Jeffrey Scott [Sandia National Laboratorie↗
NASA NTRS2003

Analysis of Surface and Bulk Behavior in Ni-Pd Alloys

The most salient features of the surface structure and bulk behavior of Ni-Pd alloys have been studied using the BFS method for alloys. Large-scale atomistic simulations were performed to investigate surface segregation profiles as a function of temperature, crystal face, and composition. Pd enrichment of the first layer was observed in (111) and (100) surfaces, and enrichment of the top two layers occurred for (110) surfaces. In all cases, the segregation profile shows alternate planes enriched and depleted in Pd. In addition, the phase structure of bulk Ni-Pd alloys as a function of temperature and composition was studied. A weak ordering tendency was observed at low temperatures, which helps explain the compositional oscillations in the segregation profiles. Finally, based on atom-by-atom static energy calculations, a comprehensive explanation for the observed surface and bulk features will be presented in terms of competing chemical and strain energy effects.

Ordering↗
DOE OSTI2026

Characterization of Flashback and Flame-Holding in a Jet-in-Crossflow Mixing Configuration With Methane–Hydrogen Fuel Blends

An approach that combines experimental and numerical analyses has been implemented to characterize the fundamentals of flashback events and flame-holding phenomena during high-hydrogen combustion in a jet-in-crossflow (JICF) configuration. Such flame dynamics are visualized experimentally using nanosecond (ns)-based hydroxyl planar laser-induced fluorescence (OH-PLIF) and chemiluminescence diagnostics techniques. The JICF burner has an optically accessible premixing tube allowing the optical diagnostics. The testing was conducted for varied premixer velocities (V) and equivalence ratio (ϕ) for 90–100% (H 2 , by mole) H 2 /CH 4 reactant mixtures at atmospheric temperature and pressure conditions. Two distinct flashback events were identified – conventional rich flashback and lean flashback, recorded while increasing ϕ and decreasing ϕ, respectively. The cause of lean flashback was attributed to lower momentum flux ratio which bends the jet sharply, closer toward the injection plane. The mean OH-PLIF images characterized the flame-holding behavior, where the flame was found to be stabilized on the leeward side only or on both windward and leeward sides as a lifted flame near the fuel port. A large eddy simulation (LES) with detailed chemistry combustion modeling approach was implemented along with an OH* submechanism, and it showed qualitative agreement with the integrated line-of-sight chemiluminescence results as well as the planar OH-PLIF measurements.

chemiluminescence diagnostic↗
NASA NTRSDate not supplied

High-Resolution Global Ocean Simulation in the SWOT ERA: Introducing ECCO Llc4320v2

The Surface Water and Ocean Topography (SWOT) mission is opening up new opportunities to evaluate and drive improvement in km-scale global ocean and atmosphere-ocean models via its simultaneous swath observations of sea-surface height, wind speed, and significant wave height down to 2 km scales globally. Motivated by this, we have initiated a collaborative effort at NASA Ames to run a new global ocean-sea-ice simulation in MITgcm at a nominal 1/48-deg resolution of the year 2023 that enables direct comparisons with SWOT observations. New features relative to a prior simulation at a similar resolution of the year 2012 include much more realistic tidal forcing, thermodynamic coupling in ice-shelf cavities, high-frequency hourly atmospheric forcing, much improved river runoff, double the vertical resolution in the deep ocean (173 vs 90 levels), updated bathymetry including improvements in finescale bathymetry derived from SWOT, and updated vertical mixing scheme parameters that improve the equatorial Pacific circulation.

ECCO↗
DOE OSTI2023

Fabrication of Large-Area Metal-on-Carbon Catalytic Condensers for Programmable Catalysis

Catalytic condensers stabilize charge on either side of a high-k dielectric film to modulate the electronic states of a catalytic layer for the electronic control of surface reactions. Here, carbon sputtering provided for fast, large-scale fabrication of metal–carbon catalytic condensers required for industrial application. Carbon films were sputtered on HfO 2 dielectric/p-type Si with different thicknesses (1, 3, 6, and 10 nm), and the enhancement of conductance and capacitance of carbon films was observed upon increasing the carbon thickness following thermal treatment at 400 °C. After Pt deposition on the carbon films, the Pt catalytic condenser exhibited a high capacitance of ∼210 nF/cm 2 that was maintained at a frequency ∼1000 Hz, satisfying the requirement for a dynamic catalyst to implement catalytic resonance. Temperature-programmed desorption of carbon monoxide yielded CO desorption peaks that shifted in temperature with the varying potential applied to the condenser (−6 or +6 V), indicating a shift in the binding energy of carbon monoxide on the Pt condenser surface. A substantial increase in capacitance (∼2000 nF/cm 2 ) of the Pt-on-carbon devices was observed at elevated temperatures of 400 °C that can modulate ∼10% of charge per metal atom when 10 V potential was applied. A large catalytic condenser of 42 cm 2 area Pt/C/HfO 2 /Si exhibited a high capacitance of 9393 nF with a low leakage current/capacitive current ratio (<0.1), demonstrating the practicality and versatility of the facile, large-scale fabrication method for metal–carbon catalytic condensers.

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

Biocatalyst discovery and design for plastics deconstruction: A multi‐scale perspective

Plastic waste accumulation poses significant environmental challenges due to a lack of economical solutions for the molecular deconstruction of diverse synthetic polymers. Biological‐based degradation offers promise but is hindered by the crystallinity, hydrophobicity, and additive complexity of plastics, which restrict biocatalyst access and activity. To address these problems, we propose a multi‐scale framework that combines detailed materials characterization, optimization of plastic‐biomolecular interfacial interactions, and enhancement of biocatalytic kinetics to develop effective plastic‐deconstructing enzymes. This approach leverages principles from reaction kinetics, transport and interfacial phenomena, and enzyme engineering to systematically address barriers across diverse plastic types. Our framework aims to accelerate the discovery and optimization of biocatalysts capable of scalable, selective, and efficient deconstruction of plastic waste. These advances hold potential to enable sustainable biological recycling and upcycling pathways, contributing to global efforts in mitigating plastic pollution and promoting circular material economies.

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

A Probabilistic Approach to Load Modeling for Central HVAC Systems in Large Commercial Buildings for Retrofit Decisions Under Uncertainty

Retrofitting central HVAC systems in large commercial buildings with advanced technologies like heat recovery chillers (HRCs) offers a significant opportunity to enhance energy efficiency. However, analyzing these retrofits is challenging with traditional whole-building simulation tools, which require intensive calibration and struggle to model innovative system configurations and controls. To overcome these limitations, this study proposes a load profilebased retrofit analysis framework that provides better decisions under uncertainty. The main focus of this paper is the development of a probabilistic load profile model that can be used in the framework by using exploratory data analysis (EDA) of measured building data to properly quantify its inherent variability. A non-parametric Gaussian Process (GP) model was employed to capture the time- and weather-dependent characteristics of the heating load while explicitly modeling its uncertainty. The model's effectiveness is demonstrated through strong predictive performance on unseen data and physically interpretable insights into load behavior. This data-driven, probabilistic load profile serves as a robust and flexible input for subsequent system simulations, enabling a more confident and statistically sound analysis of retrofit potential.

Ham, S W↗
DOE OSTI2023

Dense nuclear matter equation of state from heavy-ion collisions

The nuclear equation of state (EOS) is at the center of numerous theoretical and experimental efforts in nuclear physics. With advances in microscopic theories for nuclear interactions, the availability of experiments probing nuclear matter under conditions not reached before, endeavors to develop sophisticated and reliable transport simulations to interpret these experiments, and the advent of multi-messenger astronomy, the next decade will bring new opportunities for determining the nuclear matter EOS, elucidating its dependence on density, temperature, and isospin asymmetry. Among controlled terrestrial experiments, collisions of heavy nuclei at intermediate beam energies (from a few tens of MeV/nucleon to about 25 GeV/nucleon in the fixed-target frame) probe the widest ranges of baryon density and temperature, enabling studies of nuclear matter from a few tenths to about 5 times the nuclear saturation density and for temperatures from a few to well above a hundred MeV, respectively. Collisions of neutron-rich isotopes further bring the opportunity to probe effects due to the isospin asymmetry. However, capitalizing on the enormous scientific effort aimed at uncovering the dense nuclear matter EOS, both at RHIC and at FRIB as well as at other international facilities, depends on the continued development of state-of-the-art hadronic transport simulations. Furthermore, this white paper highlights the essential role that heavy-ion collision experiments and hadronic transport simulations play in understanding strong interactions in dense nuclear matter, with an emphasis on how these efforts can be used together with microscopic approaches and neutron star studies to uncover the nuclear EOS.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Exact block encoding of imaginary time evolution with universal quantum neural networks

We develop a constructive approach to generate quantum neural networks capable of representing the exact thermal states of all many-body qubit Hamiltonians. The Trotter expansion of the imaginary time propagator is implemented through an exact block encoding by means of a unitary, restricted Boltzmann machine architecture. Marginalization over the hidden-layer neurons (auxiliary qubits) creates the nonunitary action on the visible layer. Then, we introduce a unitary deep Boltzmann machine architecture in which the hidden-layer qubits are allowed to couple laterally to other hidden qubits. We prove that this wave-function is closed under the action of the imaginary time propagator and, more generally, can represent the action of a universal set of quantum gate operations. We provide analytic expressions for the coefficients for both architectures, thus enabling exact network representations of thermal states without stochastic optimization of the network parameters. In the limit of large imaginary time, the yields the ground state of the system. The number of qubits grows linearly with the number of interactions and total imaginary time for a fixed interaction order. Both networks can be readily implemented on quantum hardware via midcircuit measurements of auxiliary qubits. If only one auxiliary qubit is measured and reset, the circuit depth scales linearly with imaginary time and number of interactions, while the width is constant. Alternatively, one can employ a number of auxiliary qubits linearly proportional to the number of interactions, and circuit depth grows linearly with imaginary time only. Every midcircuit measurement has a postselection success probability, and the overall success probability is equal to the product of the probabilities of the midcircuit measurements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗
DOE OSTI2026

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)↗
DOE OSTI2026

Conformational Isomerization of Imide Anions Governs Solvation and Transport in Water-in-Salt Electrolytes

The behavior of highly concentrated electrolytes departs radically from the dilute-solution theory, yet the molecular origin of this transformation remains unresolved. Here, we identify the conformational isomerization of molecular ions as a decisive, previously unrecognized control parameter governing structure and transport in crowded aqueous electrolytes. Across a series of fluorosulfonimide anions, we show that increasing concentration drives a collective shift from extended transoid to compact cisoid conformers, revealed by small-angle X-ray scattering, vibrational spectroscopy, pulsed-field gradient NMR, and molecular dynamics simulations. This conformational transition triggers a collapse of the hydrogen-bonded water network and the emergence of densely packed ionic domains with confined water, producing a qualitative change in Li+ transport from solvent-mediated diffusion to network-confined hopping. Anion size and asymmetry systematically tune the onset of this transition, demonstrating that molecular geometry dictates mesoscale organization and dynamics in the ion-rich regime. Our results establish ion conformation, not merely composition or coordination, as a fundamental thermodynamic variable in concentrated solutions, providing a chemical framework that unifies solvation structure and transport in water-in-salt electrolytes and suggesting new principles for designing dense ionic media.

Nguyen, Huong TD↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗
DOE OSTI2026

Fingerprinting Superconductors by Disentangling Andreev and Quasiparticle Currents Across Tunable Tunnel Junctions

Tunneling Andreev reflection (TAR) spectroscopy provides a new approach to identify superconducting pairing symmetry at the atomic scale. Using atomistic superconducting transport simulations, we reveal the mechanism by which TAR can distinguish between pairing symmetries, which is complementary to traditional conductance-based techniques. In particular, owing to the additivity of the excess tunneling decay rate, the TAR spectrum is a weighted average of the contributions from quasiparticle currents, Andreev reflection, and higher-order scattering processes, and their relative weights depend on both the superconducting order parameter and the coupling strength. Within the local tunneling model, TAR dominates mid-gap conductance for s-wave superconductors, is suppressed for d-wave, and coexists with quasiparticle tunneling in sign-changing symmetries if the expectation value for the superconducting gap remains finite. Meanwhile, higher-order processes generally enhance the TAR signal when GN exceeds approximately 0.1G0. As a result, TAR provides a rich spectral fingerprint of the underlying pairing symmetry and electronic structure, enabling atomically resolved identification of unconventional superconducting states.

Maksymovych, Petro [Clemson University]↗