Search NASA⌕ Search

SEARCH · Search NASA

Results for “nodes”

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

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

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

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

PySIDT: Subgraph Isomorphic Decision Trees for Molecular Property Prediction

Accurate molecular property prediction is important across all fields of chemistry. Deep neural networks (DNNs) have become increasingly popular due to their ability to train automatically, avoiding the incredibly tedious process of constructing and extending traditional property estimation schemes. However, DNNs require large amounts of training data, are challenging to interpret, require large amounts of memory to load even during inference, and have severe difficulties incorporating qualitative chemical knowledge, which are often desired for molecular property prediction tasks. Here, in this study, we present PySIDT (https://github.com/zadorlab/PySIDT), a software for training and running inference on Subgraph Isomorphic Decision Trees (SIDTs). SIDTs are graph-based decision trees made of nodes associated with molecular substructures. Inference is done by descending target molecular structures down the decision tree to nodes with matching subgraph isomorphic substructures and making predictions based on the final (most specific) nodes matched. SIDTs scale down well to dataset sizes much smaller than is feasible for DNNs. As trees of molecular substructures, SIDTs are inherently readable and easy to visualize, making them easy to analyze. They are also straightforward to extend and retrain, facilitate uncertainty estimation, and enable easy integration of expert knowledge. We demonstrate the SIDT approach discussing its application to a diverse range of molecular prediction tasks: rate coefficient estimation, diffusion coefficient estimation, thermochemistry estimation, transition state bond stretch prediction, p K a prediction, stability of molecular structures, stability of surface structures, and prediction of surface lateral interaction energetics. Additionally, we demonstrate the power of the SIDT algorithms in two direct learning curve vanilla comparisons with the popular DNN-based software Chemprop and the popular gradient boosted trees-based software XGBoost on enthalpy of formation and rate coefficient prediction tasks. In particular, in the enthalpy of formation case, vanilla PySIDT is able to outperform vanilla Chemprop and XGBoost across the full range of training/validation set sizes out to 11,560 data points.

Johnson, Matthew Sean [Sandia National Laboratorie↗
DOE OSTI2020

Direct Identification of Mixed-Metal Centers in Metal–Organic Frameworks: Cu 3 (BTC) 2 Transmetalated with Rh 2+ Ions

Raman spectroscopy was used to establish direct evidence of heterometallic metal centers in a metal–organic framework (MOF). The Cu 3 (BTC) 2 MOF HKUST-1 (BTC 3– = benzenetricarboxylate) was transmetalated by heating it in a solution of RhCl 3 to substitute Rh 2+ ions for Cu 2+ ions in the dinuclear paddlewheel nodes of the framework. In addition to the Cu–Cu and Rh–Rh stretching modes, Raman spectra of (Cu x Rh 1– x ) 3 (BTC) 2 show the Cu–Rh stretching mode, indicating that mixed-metal Cu–Rh nodes are formed after transmetalation. Density functional theory studies confirmed the assignment of a Raman peak at 285 cm –1 to the Cu–Rh stretching vibration. Electron paramagnetic resonance spectroscopy experiments further supported the conclusion that Rh 2+ ions are substituted into the paddlewheel nodes of Cu 3 (BTC) 2 to form an isostructural heterometallic MOF, and electron microscopy studies showed that Rh and Cu are homogeneously distributed in (Cu x Rh 1– x ) 3 (BTC) 2 on the nanoscale.

02 PETROLEUM↗
DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗
DOE OSTI2020

Structure and Reactivity of Single-Site Vanadium Catalysts Supported on Metal–Organic Frameworks

Single-site catalysts with active species anchored over metal–organic framework (MOF) supports have garnered significant interest in recent years. Catalysts with vanadium oxide (VO x ) species immobilized over Zr-NU-1000 and Hf-MOF-808 (which have nodes containing six Zr (IV) or Hf (IV) ions, respectively) were recently characterized experimentally and were reported to be active for the selective oxidation of benzyl alcohol to benzaldehyde. Here, we report a detailed computational investigation of these VO x -incorporated MOF catalysts (V-MOF) by employing density functional theory. Based on the mode of VO x attachment, various structures are explored, and their relative stabilities and computed IR spectral features are reported. Mechanisms for the selective oxidation reaction are investigated. The analysis of electron flow in the turnover-limiting C–H activation step shows that the hydrogen-atom abstraction process involves a concerted proton-coupled electron-transfer mechanism. Here, the results of this study suggest that in addition to the identity of the metal in the node, the MOF node architecture plays a crucial role in driving the catalytic activities of V-MOFs.

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

Modulating Chemical Environments of Metal–Organic Framework-Supported Molybdenum(VI) Catalysts for Insights into the Structure–Activity Relationship in Cyclohexene Epoxidation

Solid supports are crucial in heterogeneous catalysis due to their profound effects on catalytic activity and selectivity. However, elucidating the specific effects arising from such supports remains challenging. We selected a series of metal–organic frameworks (MOFs) with 8-connected Zr 6 nodes as supports to deposit molybdenum(VI) onto to study the effects of pore environment and topology on the resulting Mo-supported catalysts. As characterized by X-ray absorption spectroscopy (XAS) and single-crystal X-ray diffraction (SCXRD), we modulated the chemical environments of the deposited Mo species. For Mo-NU-1000, the Mo species monodentately bound to the Zr 6 nodes were anchored in the microporous c-pore, but for Mo-NU-1008 they were bound in the mesopore of Mo-NU-1008. Both monodentate and bidentate modes were found in the mesopore of Mo-NU-1200. Cyclohexene epoxidation with H 2 O 2 was probed to evaluate the support effect on catalytic activity and to unveil the resulting structure–activity relationships. SCXRD and XAS studies demonstrated the atomically precise structural differences of the Mo binding motifs over the course of cyclohexene epoxidation. No apparent structural change was observed for Mo-NU-1000, whereas the monodentate mode of Mo species in Mo-NU-1008 and the monodentate and bidentate Mo species in Mo-NU-1200 evolved to a new bidentate mode bound between two adjacent oxygen atoms from the Zr 6 node. This work demonstrates the great advantage of using MOF supports for constructing heterogeneous catalysts with modulated chemical environments of an active species and elucidating structure–activity relationships in the resulting reactions.

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

Proton, Electron, and Hydrogen-Atom Transfer Thermodynamics of the Metal–Organic Framework, Ti-MIL-125, Are Intrinsically Correlated to the Structural Disorder

Interfacial charge transfer reactions involving protons and/or electrons are fundamental to heterogeneous catalysis and many other reactions relevant to energy, chemical, and biological sectors. Metal–organic frameworks (MOFs) with redox-active metal-oxo nodes have emerged as candidate materials to examine these reactions with near-atomic-level precision, given their crystalline nature. Here, we employed a colloidally stable, Ti-based MOF, Ti-MIL-125, with different crystal sizes to examine catalytically relevant charge transfer thermodynamics. The Ti 8 (μ 2 -O) 8 (μ 2 -OH) 4 nodes structurally mimic TiO 2 , which has shown some PCET reactivity toward reactions of H 2 , O 2 , and others. In this report, we have demonstrated that a change in crystal size induces different amounts of structural disorder to the Ti-oxo node, further changing the thermodynamics of proton/electron/hydrogen-atom transfer reactions. Using electrochemical open-circuit potential (E OCP ) measurements, we have determined that all crystallites undergo a 1H + /1e – redox reaction, which, given the stoichiometry, can be considered as a net H atom transfer (HAT) reaction. The thermodynamics of this HAT reaction, the Ti 3+ O–H bond dissociation free energy (BDFE), was dependent on the crystal size of the MOF, as the decrease in crystal size induced more structural disorder. Our computational calculations have indicated that this difference in BDFE is due to a local change in the geometry of Ti cations, rather than the commonly invoked defects, such as the “missing-linker” defect sites. Individual proton/electron transfer (PT/ET) thermodynamics were also highly dependent on the crystal sizes. These were probed using pK a or band gaps (E g ), respectively. These findings suggest that, particularly when MOFs are nanosized with a large amount of structural disorder, they should no longer be considered “true” single-site catalysts; this is an implicit, but widespread assumption within the MOF-based catalysis field. Implications of these findings will be contrasted with structurally similar metal oxides like TiO 2 and other redox-active MOFs.

Bond dissociation free energy↗
DOE OSTI2021

Explore Spatio‐Temporal Learning of Large Sample Hydrology Using Graph Neural Networks

Abstract Streamflow forecasting over gauged and ungauged basins play a vital role in water resources planning, especially under the changing climate. Increased availability of large sample hydrology data sets, together with recent advances in deep learning techniques, has presented new opportunities to explore temporal and spatial patterns in hydrological signatures for improving streamflow forecasting. The purpose of this study is to adapt and benchmark several state‐of‐the‐art graph neural network (GNN) architectures, including ChebNet, Graph Convolutional Network (GCN), and GraphWaveNet, for end‐to‐end graph learning. We explicitly represent river basins as nodes in a graph, learn the spatiotemporal nodal dependencies, and then use the learned relations to predict streamflow simultaneously across all nodes in the graph. The efficacy of the developed GNN models is investigated using the Catchment Attributes and MEteorology for Large‐sample Studies (CAMELS) data set under two settings, fixed graph topology (transductive learning), and variable graph topology (inductive learning), with the latter applicable to prediction in ungauged basins (PUB). Results indicate that GNNs are generally robust and computationally efficient, achieving similar or better performance than a baseline model trained using the long short‐term memory (LSTM) network. Further analyses are conducted to interpret the graph learning process at the edge and node levels and to investigate the effect of different model configurations. We conclude that graph learning constitutes a viable machine learning‐based method for aggregating spatiotemporal information from a multitude of sources for streamflow forecasting

Sun, Alexander Y.↗
DOE OSTI2022

A Comparison of Linear Solvers for Resolving Flow in Three-Dimensional Discrete Fracture Networks

We compare various methods for resolving steady flow within three-dimensional discrete fracture networks, including direct methods, Krylov subspace methods with and without preconditioning, and multi-grid methods. We compared the performance of the methods based on compute times and scaling of the solution as a function of the number of grid nodes and log-variance of the hydraulic aperture. The methods are applied to three test cases: (a) variable density of networks with a truncated power-law distribution of fracture lengths, (b) a fixed network composed of monodisperse fracture sizes but varied permeability/aperture heterogeneity, (c) and a network based on field site in Nevada, US. We chose these cases to allow us to study the impact of the mesh size and flow properties, as well as to demonstrate our conclusions on a large-scale, realistic problem (more than 40 million mesh nodes). A direct solution using Cholesky factorization outperformed other methods for every example but was closely followed in performance by some algebraic multigrid (AMG) preconditioned Krylov subspace methods. Among the Krylov methods, conjugate gradients (CG) with an AMG preconditioner performs the best. Generally, Cholesky factorization is recommended, but CG with an AMG preconditioner may be suitable for very large problems beyond 40 million nodes where the entire linear system cannot reside in memory.

58 GEOSCIENCES↗
DOE OSTI2023

Synoptic Drivers of Atmospheric River Induced Precipitation Near Dronning Maud Land, Antarctica

Atmospheric rivers (ARs) that reach the Antarctic Ice Sheet (AIS) transport anomalous moisture from lower latitudes and can impact the AIS via extreme precipitation and increased downward longwave radiation. ARs contribute significantly to the interannual variability of precipitation over the AIS and thus are likely to play a key role in understanding future changes in the surface mass balance of the AIS. Dronning Maud Land (DML) is one of four maxima in AR frequency over coastal East Antarctica, with AR precipitation explaining 77% of the interannual variability in precipitation for this region. We employ a 16-node self-organizing map (SOM) trained with MERRA-2 sea-level pressure anomalies to identify synoptic-scale environments associated with landfalling ARs in and around DML. Node composites of atmospheric variables reveal common drivers of precipitation associated with ARs reaching DML including anomalous high-low surface pressure couplets, anomalously high integrated water vapor, and coastal barrier jets. Using a quasi-geostrophic framework, we find that upward vertical motion associated with the occlusion process of attendant cyclones dominates atmospheric lift in AR environments. We further identify mechanisms that explain the variability in AR precipitation intensity across nodes, such as the lift associated with the occlusion process of attendant cyclones and the spatial coincidence of ascent induced by the occlusion process and frontogenesis. Here, the latter suggests that ARs making landfall during the mature phase of cyclogenesis result in higher precipitation intensity compared to landfalling ARs that occur during the occluded phase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Examining Infrasound Propagation at High Spatial Resolution Using a Nodal Seismic Array

Infrasound—acoustic waves in the atmosphere below 20 Hz—is a useful monitoring tool. Topography and atmospheric structure strongly control infrasound propagation, and at common source–receiver distances neither of these effects can be ignored when quantitative source constraints are sought. Detailed spatial measurements of the infrasound wavefield would inform propagation models and improve source estimates. However, the “large-N” deployment strategy now well-known in seismology has not yet been realized for infrasound studies. Here, we use the 900-node seismic array from the 2014 Imaging Magma Under St. Helens (iMUSH) experiment as a proxy for a large-N infrasound network, by leveraging acoustic–seismic coupled arrivals. The active-source component of iMUSH consisted of 23 shallowly buried explosions around Mount Saint Helens volcano; these explosions produced epicentral infrasound recorded on the nodes. We find that the bulk presence of ground-coupled infrasound on the nodes is controlled by wind noise and source–receiver distance, with observed arrivals for eight explosions. Explosions with the most extensive coupling produce complex spatial waveform patterns across the array. These patterns are related to both topographic and atmospheric propagation effects, as well as spatially variable site (coupling) effects. We compare our observations to simple topographic diffraction and high-resolution wind advection models, and full-wave numerical simulations. We find strong spatial correlations between (a) coupled arrival strength and modeled topographic obstruction and (b) coupled arrival time and along-path winds. Our seismoacoustic analyses and results are applicable to other existing and future nodal seismic data sets and can expand the utility of such deployments.

58 GEOSCIENCES↗
DOE OSTI2024

Nonlinear optical diode effect in a magnetic Weyl semimetal

Diode effects are of great interest for both fundamental physics and modern technologies. Electrical diode effects (nonreciprocal transport) have been observed in Weyl systems. Optical diode effects arising from the Weyl fermions have been theoretically considered but not probed experimentally. Here, we report the observation of a nonlinear optical diode effect (NODE) in the magnetic Weyl semimetal CeAlSi, where the magnetization introduces a pronounced directionality in the nonlinear optical second-harmonic generation (SHG). We demonstrate a six-fold change of the measured SHG intensity between opposite propagation directions over a bandwidth exceeding 250 meV. Supported by density-functional theory, we establish the linearly dispersive bands emerging from Weyl nodes as the origin of this broadband effect. We further demonstrate current-induced magnetization switching and thus electrical control of the NODE. Our results advance ongoing research to identify novel nonlinear optical/transport phenomena in magnetic topological materials and further opens new pathways for the unidirectional manipulation of light.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Binding and excitations in Si x H y molecular systems using quantum Monte Carlo

In this work, we present high-accuracy correlated calculations of small SixHy molecular systems in both the ground and excited states. We employ quantum Monte Carlo (QMC) together with a variety of many-body wave function approaches based on basis set expansions. The calculations are carried out in a valence-only framework using recently derived correlation consistent effective core potentials. Our primary goal is to understand the fixed-node diffusion QMC errors in both the ground and excited states with single-reference trial wave functions. Using a combination of methods, we demonstrate the very high accuracy of the QMC atomization energies being within ≈0.07 eV or better when compared with essentially exact results. By employing proper choices for trial wave functions, we have found that the fixed-node QMC biases for total energies are remarkably uniform ranging between 1% and 3.5% with absolute values at most ≈0.2 eV across the systems and several types of excitations such as singlets and triplets as well as low-lying and Rydberg-like states. Our results further corroborate that Si systems, and presumably also related main group IV and V elements of the periodic table (Ge, Sn, etc), exhibit some of the lowest fixed-node biases found in valence-only electronic structure QMC calculations.

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

Single photon sources with near unity collection efficiencies by deterministic placement of quantum dots in nanoantennas

Deterministic coupling between photonic nodes in a quantum network is an essential step toward implementing various quantum technologies. The omnidirectionality of free-standing emitters, however, makes this coupling highly inefficient, in particular if the distant nodes are coupled via low numerical aperture (NA) channels such as optical fibers. This limitation requires placing quantum emitters in nanoantennas that can direct the photons into the channels with very high efficiency. Moreover, to be able to scale such technologies to a large number of channels, the placing of the emitters should be deterministic. In this work, we present a method for directly locating single free-standing quantum emitters with high spatial accuracy at the center of highly directional bullseye metal–dielectric nanoantennas. We further employ non-blinking, high quantum yield colloidal quantum dots for on-demand single-photon emission that is uncompromised by instabilities or non-radiative exciton recombination processes. Taken together, this approach results in a record-high collection efficiency of 85% of the single photons into a low NA of 0.5, setting the stage for efficient coupling between on-chip, room temperature nanoantenna-emitter devices and a fiber or a remote free-space node without the need for additional optics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Precise Fermi level engineering in a topological Weyl semimetal via fast ion implantation

The precise controllability of the Fermi level is a critical aspect of quantum materials. For topological Weyl semimetals, there is a pressing need to fine-tune the Fermi level to the Weyl nodes and unlock exotic electronic and optoelectronic effects associated with the divergent Berry curvature. However, in contrast to two-dimensional materials, where the Fermi level can be controlled through various techniques, the situation for bulk crystals beyond laborious chemical doping poses significant challenges. Here, we report the milli-electron-volt (meV) level ultra-fine-tuning of the Fermi level of bulk topological Weyl semimetal tantalum phosphide using accelerator-based high-energy hydrogen implantation and theory-driven planning. By calculating the desired carrier density and controlling the accelerator profiles, the Fermi level can be experimentally fine-tuned from 5 meV below, to 3.8 meV below, to 3.2 meV above the Weyl nodes. High-resolution transmission electron microscopy reveals the crystalline structure is largely maintained under irradiation, while electrical transport indicates that Weyl nodes are preserved and carrier mobility is also largely retained. Our work demonstrates the viability of this generic approach to tune the Fermi level in semimetal systems and could serve to achieve property fine-tuning for other bulk quantum materials with ultrahigh precision.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Wave function analysis with a maximum flow algorithm

An efficient algorithm for computing the maximum-flow path in a network is applied to the identification of the dominant configuration state functions (CSFs) in a graphically contracted function (GCF), configuration interaction, wave function. The flow network is a space of spin-adapted CSFs represented by a Shavitt graph, wherein the nodes correspond to orbital occupations and spin quantum numbers. The graph nodes are connected by arcs, and an arc density is defined as sums of the associated squared CSF coefficients. A max-min approach determines an upper bound to the maximum possible incoming flow for each graph node. A backtracking step generates a candidate walk and is followed by a limited search of alternative branching paths for the dominant CSF. The arc density contributions are removed from the graph, and the algorithm is reapplied to the updated graph. This list of generated walks can be partitioned in order to guarantee that the dominant CSFs have been identified. All of the steps in this algorithm are computationally efficient and do not depend on the potentially large dimension of the underlying linear CSF expansion space. An analysis of low-lying valence states of C-2 illustrates the method.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Coupling of CTF and RELAP5-3D Within an Enhanced Fidelity Nuclear Power Plant Simulator

A robust and accurate multiphysics engineering simulator is being developed to model the core behavior and system response of pressurized water reactors. This simulator relies on the NESTLE and CTF computer codes to model the neutronics and thermal hydraulics (TH), respectively, inside the core on a nodal scale and on the Reactor Excursion and Leak Analysis Program—Three Dimensional (RELAP5-3D) to model the entire nuclear steam supply system. The RELAP5-3D model includes highly detailed nodalization and multidimensional flow modeling throughout the vessel. Previously, pin-resolved data generated via the Virtual Environment for Reactor Analysis core simulator were used to improve the accuracy of the NESTLE core predictions. The engineering simulator being developed as part of this work uses the 3KEYMASTER platform to couple the enhanced NESTLE model to a nodal-fidelity CTF model to balance run time with accuracy; NESTLE provides node-dependent powers to CTF, and CTF provides node-dependent coolant densities and fuel temperatures to NESTLE.An overlapping domain approach is used for the core TH in which RELAP5-3D provides core boundary conditions based on the system response and CTF provides a node-dependent coolant heating rate to the RELAP5-3D core solution. In the preliminary TH demonstration discussed in this paper, CTF and RELAP5-3D provided similar steady-state core predictions, indicating the hydraulic compatibility between the codes, as well as reasonable and expected behavior under hypothetical transient conditions. This provides an initial step in ongoing efforts toward a robust, multiscale TH/neutronics engineering simulator capability.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗