SEARCH · Search NASA

Results for “computer simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

464 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗

DOE OSTI2026

SimH 2 : an integrated techno-economic modeling framework for hydrogen pipeline infrastructure and network optimization

Large-scale hydrogen (H 2 ) pipeline transport design and network optimization have seldom been reported due to the lack of a cost model accounting for the relationship between transport cost and hydrogen mass flow rate. Here, this work introduced a system-level cost model for hydrogen pipeline transport at supercritical state and integrated it with an existing CO 2 pipeline network tool, SimCCS, for hydrogen-specific pipeline design and optimization. The Intermountain West (I-West) region of the U.S., historically dependent on fossil fuel-based economies, is chosen to demonstrate the capabilities of our H 2 pipeline cost model and transport network optimization platform called SimH 2 . Two scenarios are examined: one where the pipeline is not allowed to pass through disadvantaged communities and the other where it is permitted. The results highlight that incorporating disadvantaged-community constraints lead to longer pipeline routes and increased transport costs, reflecting the trade-offs involved in equitable infrastructure development. It is demonstrated that the newly developed SimH 2 tool not only enables the efficient design of H 2 transportation pipelines but also optimizes the network by accounting for local terrain and the presence of disadvantaged areas.

08 HYDROGEN↗

DOE OSTI2026

ClassNMSW- a real-time classification approach for non-recycled municipal solid waste using hyperspectral imaging

Real-time classification of non-recycled municipal solid waste (NMSW) is essential for efficient valorization. This study introduces ClassNMSW, a comprehensive framework for classifying 22 NMSW subclasses under industrial constraints by using hyperspectral imaging (HSI). A primary innovation of this work is the development of a variance-controlled spectral extraction algorithm. Unlike traditional methods that rely on simple averaging, this approach systematically investigates the extent of pixel extraction to minimize the loss of critical chemical information while maximizing data reduction thus ensuring high spectral fidelity with low computational cost. The approach developed in this work integrates automated, computer-vision-based background removal, eliminating the need for the manual thresholding common in current literature. To resolve ambiguities among chemically similar subclasses, a tiered classification and multi-camera fusion strategy (NIR17 and NIR22) is implemented. Results demonstrate that ClassNMSW achieves an object-wise weighted accuracy of 98.70% for single-sensor configurations and 100% under sensor fusion. A novel rolling-window strategy satisfies desired end-to-end latency of <2 s, satisfying the strict deterministic requirements of high-speed industrial sorting environments. The ClassNMSW framework provides a scalable foundation for advancing circularity and resource recovery in large-scale waste valorization operations.

99 - GENERAL AND MISCELLANEOUS↗

DOE OSTI2026

First Wall Design of a Tokamak Pilot Plant Using a Monte Carlo Model for 3-D Heat Flux Deposition

We present a method for calculating the heat fluxes deposited on nonaxisymmetric tokamak first wall components, allowing for a first-of-its-kind model for power handling in the tokamak far scrape-off layer (SOL). The DIV3D Monte Carlo model features strict global power conservation and can calculate the finite cross-field plasma transport into magnetically-shadowed regions, which is significant when dealing with meter-scale shadows introduced by components such as poloidal limiters or antennas. As a case study, we apply the DIV3D model to inform the distribution of first wall poloidal limiters in an ARC-class reactor device. We demonstrate that discrete protection limiters can efficiently reduce peak heat fluxes on recessed breeder wall components in the presence of significant far-SOL plasma fluxes. By varying the toroidal periodicity and radial standoff depth of the limiters, we demonstrate one of the tradeoffs that must be considered in first wall design: more limiters provide greater protection, but at the cost of reduced breeding performance. We also present the impact that radial misalignments between limiters would have on first wall power loading.

Monte Carlo methods↗

DOE OSTI2026

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗

DOE OSTI2026

Integrating AI Data Centers with the Power Grid

The rapid expansion of artificial intelligence (AI) has triggered an unprecedented surge in electricity demand, with US data center energy use projected to double or triple 2023 levels by 2028. This exponential growth places strain on grid infrastructure, which can hinder timely construction of desired computing capacity. To bridge this supply-demand gap, utilities and AI developers are increasingly turning to demand flexibility, a strategy that incentivizes shifting or reducing power use during peak periods of grid stress. Data centers are uniquely equipped for flexible operations due to their digital workloads, built-in redundancy, and onsite energy assets. This article outlines four primary mechanisms to enable data center flexibility: computational load flexibility (shifting tasks temporally or geographically), flexible use of core facility infrastructure adjustments, energy storage utilization, and onsite electricity generation. To encourage adoption, utilities are deploying new tariff designs, including voluntary interruptible service riders, mandated flexibility requirements, and streamlined interconnection processes for flexible loads. For the highly capitalized and rapidly growing AI industry, the primary motivators for embracing these strategies are expediting facility interconnection, satisfying emerging regulatory mandates, and mitigating community resistance. While demand flexibility cannot substitute the long-term need for new bulk power generation, it serves as an essential, immediate solution for enabling near-term deployment. By transforming data centers from grid stressors into stabilizing assets, flexible operations can ensure reliable grid integration, ease market pressures, and support a resilient power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Synthesis‐Related Nanoscale Defects in Mo‐Based Janus Monolayers Revealed by Cross‐Correlated AFM and TERS Imaging

2D Janus transition metal dichalcogenides (TMDs) are promising candidates for various applications including non-linear optics, energy harvesting, and catalysis. These materials are usually synthesized via chemical conversion of pristine TMDs. Nanometer-scale characterization of the obtained Janus materials’ morphology and local composition is crucial for both the synthesis optimization and the future device applications. In this work, we present the results of cross-correlated atomic force microscopy (AFM) and tip-enhanced Raman spectroscopy (TERS) study of Janus monolayers synthesized by the hydrogen plasma-assisted chemical conversion of MoSe 2 and MoS 2 . We demonstrate that the choice of both the growth substrate and the starting TMD influences the residual strain, thereby shaping the nanoscale morphology of the resulting Janus material. Furthermore, by employing TERS imaging, we show the presence of nanoscale islands (≈20 nm across) of MoSe 2 - ${\mathrm{Mo}}_{{\mathrm{Se}}}^{\mathrm{S}}$ (MoS 2 -${\mathrm{Mo}}_{\mathrm{S}}^{{\mathrm{Se}}}$) vertical heterostructures originating from the bilayer nanoislands in the precursor monolayer crystals. The understanding of the origins of nanoscale defects in Janus TMDs revealed in this study can help with further optimization of the Janus conversion process towards uniform and wrinkle-/crack-free Janus materials. Moreover, this work shows that cross-correlated AFM and TERS imaging is a powerful and accessible method for studying nanoscale composition and defects in Janus TMD monolayers.

2D materials↗

DOE OSTI2026

Enter the AHU (36th Chamber of ASHRAE): A Multi-site Field Study of ASHRAE G36

Despite being recognized as the best practice for advanced building controls, ASHRAE Guideline 36 (G36) has seen slow adoption in retrofit cases. Decisionmakers lack credible field evidence to justify the time and person-power investment. Most prior analyses have relied on software simulations, which overlook implementation challenges and fail to persuade owners to move from models to real-world deployment. This paper presents a multi-site field study of G36 performance, drawing on measured results from 17 projects across diverse building types and climate zones. The analysis disaggregates outcomes by the most widely adopted air handling unit (AHU) based G36 strategies, including trim-and-respond approaches to supply air temperature (SAT) and duct static pressure (DSP) reset, and economizer controls. Results for controls re programming implementations are encouraging, with HVAC savings ranging from 2% - 49%, with a median of 18%. The range aligns with simulation study findings showing 1% - 46% savings through these three strategies. These findings show that even without capital investment in control infrastructure upgrades, existing building owners are reaping significant benefits from updating HVAC sequences of operation to industry best-practice solutions. The results give practitioners and decisionmakers a reference point on what to expect, which strategies deliver, and how field performance may compare to simulation.

Deshpande, Reva↗

DOE OSTI2026

Neutron detector response modeling in NOvA

Neutrons can present a significant challenge for neutrino experiments in which energy reconstruction is critical. With the ability to escape detection completely and with a weak correlation between their kinetic energy and any eventual energy deposition, it is difficult to fully account for neutrons produced in neutrino interactions. This in turn leads to significant model dependence when evaluating neutron-related systematic uncertainties. The NOvA experiment is a long-baseline neutrino oscillation experiment with a high-statistics sample of antineutrino data collected by its near detector. We report an excess relative to data of simulated neutron candidates with low energy depositions when using standard Geant4 physics lists. The simulation excess is traced to an overabundance of secondary photons produced from interactions of neutrons with kinetic energy greater than \SI{20}{\mega\eV}. Improved agreement with data is obtained by applying the data-driven neutron-on-carbon \menate model for neutrons between \SI{20}{\mega\eV} and ${\sim}$\SI{100}{\mega\eV}. With \menate, the residual oversimulation is more uniform across the calorimetric neutron energy spectrum, suggesting possible overproduction of primary neutrons by the GENIE neutrino interaction generator. These results motivate the adoption of \menate-supplemented Geant4 simulation as the nominal simulation in the production of future \nova simulation.

Abubakar, S.↗

DOE OSTI2026

Realization of two-dimensional discrete time crystals with anisotropic Heisenberg coupling

A discrete time crystal (DTC) is an out-of-equilibrium phase of matter that spontaneously breaks discrete time-translation symmetry. Previous studies have been limited to a set of models with Ising-like couplings - and mostly only in one dimension - thus precluding our understanding of the existence (or not) of DTCs in models with more realistic interactions. In this work, by combining the latest generation of IBM quantum processors with state-of-the-art tensor network methods, we demonstrate the existence of a DTC in a two-dimensional system governed by anisotropic Heisenberg interactions. We uncover a rich phase diagram encompassing spin-glass, ergodic, and time-crystalline phases, and identify the interplay of initialization, interaction anisotropy, and driving protocols in stabilizing the DTC phase. By extending the study of Floquet matter beyond simplified models, we lay the groundwork for exploring how driven systems bridge the gap between quantum coherence and emergent non-equilibrium thermodynamics.

Phase transitions and critical phenomena↗

DOE OSTI2024

Copacabana: a probabilistic membership assignment method for galaxy clusters

Cosmological analyses using galaxy clusters in optical/near-infrared photometric surveys require robust characterization of their galaxy content. Precisely determining which galaxies belong to a cluster is crucial. In this paper, we present the COlor Probabilistic Assignment of Clusters And BAyesiaN Analysis (Copacabana) algorithm. Copacabana computes membership probabilities for all galaxies within an aperture centred on the cluster using photometric redshifts, colours, and projected radial probability density functions. We use simulations to validate Copacabana and we show that it achieves up to 89 per cent membership accuracy with a mild dependence on photometric redshift uncertainties and choice of aperture size. We find that the precision of the photometric redshifts has the largest impact on the determination of the membership probabilities followed by the choice of the cluster aperture size. We also quantify how much these uncertainties in the membership probabilities affect the stellar mass–cluster mass scaling relation, a relation that directly impacts cosmology. Using the sum of the stellar masses weighted by membership probabilities (⁠μ * ⁠) as the observable, we find that Copacabana can reach an accuracy of 0.06 dex in the measurement of the scaling relation at low redshift for a Legacy Survey of Space and Time type survey. These results indicate the potential of Copacabana and μ * to be used in cosmological analyses of optically selected clusters in the future.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

EvoDiffMol: evolutionary diffusion framework for 3D molecular design with optimized properties

Designing molecules with specific target properties remains a fundamental challenge in computational chemistry. While existing approaches show promise, most rely on simplified representations like SMILES strings or 2D graphs that lack essential three-dimensional geometric information. We present EvoDiffMol, a computational framework that integrates evolutionary algorithms with three-dimensional diffusion models for property-driven molecular generation. The method operates through adaptive evolutionary optimization, where population-based selection guides the generation process toward desired property landscapes. EvoDiffMol supports both unconstrained molecular design and scaffold-constrained generation that preserves fixed substructures while optimizing complementary regions. Comprehensive evaluation demonstrates exceptional performance, achieving the highest drug-likeness score (0.94) among all compared state-of-the-art methods while maintaining excellent validity, uniqueness, and novelty. Beyond single property optimization, the framework demonstrates flexible multi-property optimization capabilities, simultaneously controlling multiple molecular descriptors including synthetic accessibility, lipophilicity, topological polar surface area, and clinically relevant ADMET properties such as cardiotoxicity (hERG) and intestinal permeability (Caco-2). This adaptability spans from simple descriptors to practical pharmaceutical endpoints without requiring complete model retraining. The framework achieves precise control over target property values, generating molecules with properties closely matching specified targets for both single and multiple descriptors. Scaffold-constrained experiments preserve fixed molecular cores while maintaining effective property optimization. The three-dimensional representation offers advantages in maintaining structural validity during iterative optimization, with potential for geometry-aware applications in materials science and drug discovery.

3D molecular generation↗

DOE OSTI2026

Machine Learning for Multipactor Susceptibility Prediction in Planar RF Gaps

Multipactor discharge is a nonlinear electron avalanche that limits the performance of high-power radio-frequency (RF) and vacuum electronic devices. Predicting multipactor susceptibility traditionally relies on Monte Carlo or particle-in-cell (PIC) simulations, which become computationally expensive for large parametric studies. In this work, we present a supervised machine-learning (ML) framework for prediction of multipactor susceptibility in a two-surface planar geometry. The models are trained using high-fidelity PIC simulation generated susceptibility data and learn the relationship between operational parameters, geometry, and material-dependent secondary electron emission properties. The proposed approach enables rapid reconstruction of susceptibility charts while preserving the physical structure of multipactor growth regions.

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Stabilizer Scars

Quantum many-body scars are eigenstates in nonintegrable isolated quantum systems that defy typical thermalization paradigms, violating the eigenstate thermalization hypothesis and quantum ergodicity. Here, we identify exact analytic scar solutions in a 2+1 dimensional lattice gauge theory in a quasi-1D limit as zero-magic resource stabilizer states. Our results also highlight the importance of magic resources for gauge theory thermalization, revealing a connection between computational complexity and quantum ergodicity.

eigenstate thermalization↗

DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗

DOE OSTI2026

Distributed Tomographic Reconstruction with Quantization

Conventional tomographic reconstruction typically depends on centralized servers for both data storage and computation, leading to concerns about memory limitations and data privacy. Distributed reconstruction algorithms mitigate these issues by partitioning data across multiple nodes, reducing server load and enhancing privacy. However, these algorithms often encounter challenges related to memory constraints and communication overhead between nodes. In this paper, we introduce a decentralized Alternating Directions Method of Multipliers (ADMM) with configurable quantization. By distributing local objectives across nodes, our approach is highly scalable and can efficiently reconstruct images while adapting to available resources. To overcome communication bottlenecks, we propose two quantization techniques based on K-means clustering and JPEG compression. Numerical experiments with benchmark images illustrate the tradeoffs between communication efficiency, memory use, and reconstruction accuracy.

Miao, Runxuan↗

DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗

DOE OSTI2026

Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions

Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.

Walker, Christopher C. [Oak Ridge National Laborat↗