Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

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 253 records · Page 14

The WACMOS-ET project – Part 1: Tower-scale evaluation of four remote-sensing-based evapotranspiration algorithms

The WAter Cycle Multi-mission Observation Strategy – EvapoTranspiration (WACMOS-ET) project has compiled a forcing data set covering the period 2005–2007 that aims to maximize the exploitation of European Earth Observations data sets for evapotranspiration (ET) estimation. The data set was used to run four established ET algorithms: the Priestley–Taylor Jet Propulsion Laboratory model (PT-JPL), the Penman–Monteith algorithm from the MODerate resolution Imaging Spectroradiometer (MODIS) evaporation product (PM-MOD), the Surface Energy Balance System (SEBS) and the Global Land Evaporation Amsterdam Model (GLEAM). In addition, in situ meteorological data from 24 FLUXNET towers were used to force the models, with results from both forcing sets compared to tower-based flux observations. Model performance was assessed on several timescales using both sub-daily and daily forcings. The PT-JPL model and GLEAM provide the best performance for both satellite- and tower-based forcing as well as for the considered temporal resolutions. Simulations using the PM-MOD were mostly underestimated, while the SEBS performance was characterized by a systematic overestimation. In general, all four algorithms produce the best results in wet and moderately wet climate regimes. In dry regimes, the correlation and the absolute agreement with the reference tower ET observations were consistently lower. While ET derived with in situ forcing data agrees best with the tower measurements ( R 2 = 0.67), the agreement of the satellite-based ET estimates is only marginally lower ( R 2 = 0.58). Results also show similar model performance at daily and sub-daily (3-hourly) resolutions. Overall, our validation experiments against in situ measurements indicate that there is no single best-performing algorithm across all biome and forcing types. In conclusion, an extension of the evaluation to a larger selection of 85 towers (model inputs resampled to a common grid to facilitate global estimates) confirmed the original findings.

54 ENVIRONMENTAL SCIENCES↗

Theoretical prediction of the band offsets at the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions using the self-consistent ab initio DFT/GGA-1/2 method

The band offsets of the ZnO/anatase TiO{sub 2} and GaN/ZnO heterojunctions are calculated using the density functional theory/generalized gradient approximation (DFT/GGA)-1/2 method, which takes into account the self-energy corrections and can give an approximate description to the quasiparticle characteristics of the electronic structure of semiconductors. We present the results of the ionization potential (IP)-based and interfacial offset-based band alignments. In the interfacial offset-based band alignment, to get the natural band offset, we use the surface calculations to estimate the change of reference level due to the interfacial strain. Based on the interface models and GGA-1/2 calculations, we find that the valence band maximum and conduction band minimum of ZnO, respectively, lie 0.64 eV and 0.57 eV above those of anatase TiO{sub 2}, while lie 0.84 eV and 1.09 eV below those of GaN, which agree well with the experimental data. However, a large discrepancy exists between the IP-based band offset and the calculated natural band offset, the mechanism of which is discussed. Our results clarify band alignment of the ZnO/anatase TiO{sub 2} heterojunction and show good agreement with the GW calculations for the GaN/ZnO heterojunction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Teaching a neural network to attach and detach electrons from molecules

Abstract Interatomic potentials derived with Machine Learning algorithms such as Deep-Neural Networks (DNNs), achieve the accuracy of high-fidelity quantum mechanical (QM) methods in areas traditionally dominated by empirical force fields and allow performing massive simulations. Most DNN potentials were parametrized for neutral molecules or closed-shell ions due to architectural limitations. In this work, we propose an improved machine learning framework for simulating open-shell anions and cations. We introduce the AIMNet-NSE (Neural Spin Equilibration) architecture, which can predict molecular energies for an arbitrary combination of molecular charge and spin multiplicity with errors of about 2–3 kcal/mol and spin-charges with error errors ~0.01e for small and medium-sized organic molecules, compared to the reference QM simulations. The AIMNet-NSE model allows to fully bypass QM calculations and derive the ionization potential, electron affinity, and conceptual Density Functional Theory quantities like electronegativity, hardness, and condensed Fukui functions. We show that these descriptors, along with learned atomic representations, could be used to model chemical reactivity through an example of regioselectivity in electrophilic aromatic substitution reactions.

36 MATERIALS SCIENCE↗

Exploring the interaction between the MW and LMC with a large sample of blue horizontal branch stars from the DESI survey

The Large Magellanic Cloud (LMC) is a Milky Way (MW) satellite that is massive enough to gravitationally attract the MW disc and inner halo, causing significant motion of the inner MW with respect to the outer halo. In this work, we probe this interaction by constructing a sample of 9866 blue horizontal branch (BHB) stars with radial velocities from the DESI spectroscopic survey out to 120 kpc from the Galactic centre. This is the largest spectroscopic set of BHB stars in the literature to date, and it contains four times more stars with Galactocentric distances beyond 50 kpc than previous BHB catalogues. Using the DESI BHB sample combined with SDSS BHBs, we measure the bulk radial velocity of stars in the outer halo and observe that the velocity in the Southern Galactic hemisphere is different by 3.7σ from the North. Modelling the projected velocity field shows that its dipole component is directed at a point 22 deg away from the LMC along its orbit, which we interpret as the travel direction of the inner MW. The velocity field includes a monopole term that is –24 km s –1 ⁠, which we refer to as compression velocity. This velocity is significantly larger than predicted by the current models of the MW and LMC interaction. This work uses DESI data from its first 2 yr of observations, but we expect that with upcoming DESI data releases, the sample of BHB stars will increase and our ability to measure the MW–LMC interaction will improve significantly.

Galaxy: evolution↗

Updated SAM Model for the Molten Salt Reactor Experiment (MSRE)

The development of reference standard problems based on prototypical reactor designs is of particular importance to verify the adequacy of computer codes and evaluation models for specific reactor types. To support the multiphysics coupled simulation of molten-salt-fueled reactor (MSR) using SAM and Griffin computer codes for safety and licensing analysis, much efforts have been put into enhancing code capabilities and developing reference models for the MSR primary loop in SAM. In this work, a previously developed Molten Salt Reactor Experiment (MSRE) primary loop model was updated to include a two-dimensional (2-D) core region and external core components in one-dimension (1-D) or zero-dimension (0-D). To ensure accurate feedback calculation in multi-physics simulations, the delayed neutron precursor tracking model and solid graphite model were added in the SAM model. In addition, the 2-D and 1-D domains are tightly coupled using the recently developed single-solve approach in SAM. The updated model has been tested under both steady-state and transient scenarios to demonstrate its potential for the multi-physics simulation of MSRE with coupled SAM and Griffin.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High Efficiency 2-Bladed Utility Wind Turbine Enhanced by CoFlow Jet Active Flow Control

The aerodynamic performance and power output of a high efficiency Co-Flow Jet (CFJ) 2-bladed Wind Turbine is studied in this paper. CFJ is a Zero-Net-Mass-Flux (ZNMF) active flow control method that dramatically increases airfoil lift coefficient and suppresses flow separation at a low energy expenditure. The 3D Reynolds Averaged Navier-Stokes (RANS) equations with one-equation Spalart-Allmaras (SA) turbulence model are solved to simulate the 3D flows of the wind turbines. The CFJ-Wind Turbine in this paper is modified from the NREL 3-bladed 3.4 MW reference wind turbine with the same RPM, blade length, distribution of chord, twist, and bending. The CFJ 2-bladed wind turbine utilizes the CFJ-NACA 64 series airfoils at most of the blade span from 31.88\% to the tip. The CFJ injection and suction slots are implemented along the full blade span to enhance the effectiveness and power output. The study shows that the 2-bladed CFJ turbine significantly increases power output compared with the 3-bladed baseline turbine at all wind speeds by 15.8\% at the wind speed of 9.8 m/s and by 25.9\% at 4 m/s respectively. The efficiency enhancement is substantially more at the low wind speed than at high wind speed. This would be very beneficial to increase the annual power output or capacity factor as wind turbines work at the speeds lower than the rated speed for most of the time. The results show that the flow field around the blade surface presents a significant enhanced suction peak near the blade leading edge, which results in a very high tangential force and thus increased wind turbine net power coefficient. Consistent with the previous study for the small NREL Phase VI research wind turbine, this paper demonstrates that CFJ 2-bladed wind turbine is an effective technology to extract more power from the wind at all speeds. This study is just an initial effort to demonstrate the advantage of the CFJ blade with no design optimization. It is believed that there is more potential to further enhance the efficiency improvement.

active flow control↗

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

Simple Approximation for the Ideal Reference State of Gases Adsorbed on Solid-State Surfaces

Reference states are useful as models for facilitating calculations of equilibrium constants, and they may also serve as standard states that are convenient for organizing and tabulating thermodynamic data; however, standard state conventions and appropriate reference states for adsorbed species have received less attention than those for pure substances and solutes. Here, we compare seven choices of reference states for calculations of equilibrium constants and transition state theory rate constants for flat surfaces, in particular (1) an ideal 2D harmonic oscillator, (2) an ideal rigid-molecule harmonic oscillator, (3) an ideal 2D harmonic oscillator with separable surface modes, (4) a 2D ideal gas, (5) an ideal 2D hindered translator, (6) an ideal 2D hindered translator with lowest-order barriers, and (7) a simple ideal 2D hindered translator proposed in this work. The advantage of models 5–7 is that they can treat both mobile and localized adsorbates in a consistent way, whereas models 1–3 are only appropriate for localized adsorbates, and model 4 is only appropriate for a freely translating adsorbate. Furthermore, models 6 and 7 reduce the computational cost without the user having to calculate barrier heights for diffusion. An advantage of the simple ideal 2D hindered translator is that it has a physical high-temperature limit. We also propose a reference state for nonflat surfaces. Here, the user is encouraged to choose a reference state based on the appropriateness of the model and the practicality of the calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

INTERSTELLAR PICKUP ION PRODUCTION IN THE GLOBAL HELIOSPHERE AND HELIOSHEATH

Interstellar pickup ions (PUIs) play a significant part in mediating the solar wind (SW) interaction with the interstellar medium. In this paper, we examine the details of spatial variation of the PUI velocity distribution function (VDF) in the SW by solving the PUI transport equation. We assume the PUI distribution is isotropic resulting from strong pitch-angle scattering by wave–particle interaction. A three-dimensional model combining the MHD treatment of the background SW and neutrals with a kinetic treatment of PUIs throughout the heliosphere and the surrounding local interstellar medium has been developed. The model generates PUI power-law tails via second-order Fermi process. We analyze how PUIs transform across the heliospheric termination shock and obtain the PUI phase space distribution in the inner heliosheath including continuing velocity diffusion. Our simulated PUI spectra are compared with observations made by New Horizons , Ulysses , Voyager 1, 2 , and Cassini , and a satisfactory agreement is demonstrated. Some specific features in the observations, for example, a cutoff of PUI VDF at v = V {sub SW} and a f ∝ v {sup -5} tail in the reference frame of the SW, are well represented by the model.

79 ASTRONOMY AND ASTROPHYSICS↗

Blade-Resolved, Single-Turbine Simulations Under Atmospheric Flow

The ultimate goal of the ExaWind project is to enable scientific discovery through predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The objective of this milestone is to demonstrate a 'production run' of a megawatt-scale turbine to exercise and demonstrate the turbine simulation capabilities that have been developed during the previous two years of the ExaWind project. Demonstrating this simulation capability is the critical first step in enabling Nalu-Wind to be used as a predictive modeling tool that can advance the scientific understanding of wind turbine performance. In Section 2, the milestone description is provided. Next, in Section 3, we present the NREL 5-megawatt (MW) reference turbine that was used to demonstrate Nalu-Wind's capabilities and describe the four computational meshes that were generated. In Section 4, we present results from simulations of the turbine operating under uniform inflow and homogeneous isotropic turbulence (HIT) conditions. Finally, in Section 5, notable achievements, conclusions, and a path forward are described.

17 WIND ENERGY↗

Final Techno-Economic Analysis of 550 MWe Supercritical PC Power Plant CO 2 Capture with Linde-BASF Advanced PCC Technology

This topical report presents the techno-economic evaluation of a 550 MWe supercritical pulverized coal (PC) power plant utilizing Illinois No. 6 coal as fuel, integrated with 1) a previously presented (for a subcritical PC plant) Linde-BASF post-combustion CO 2 capture (PCC) plant incorporating BASF’s OASE® blue aqueous amine-based solvent (LB1) [Ref. 6] and 2) a new Linde-BASF PCC plant incorporating the same BASF OASE® blue solvent that features an advanced stripper interstage heater design (SIH) to optimize heat recovery in the PCC process. The process simulation and modeling for this report is performed using Aspen Plus V8.8. Technical information from the PCC plant is determined using BASF’s proprietary thermodynamic and process simulation models. The simulations developed and resulting cost estimates are first validated by reproducing the results of DOE/NETL Case 12 representing a 550 MWe supercritical PC-fired power plant with PCC incorporating a monoethanolamine (MEA) solvent as used in the DOE/NETL Case 12 reference [Ref. 2]. The results of the techno-economic assessment are shown comparing two specific options utilizing the BASF OASE® blue solvent technology (LB1 and SIH) to the DOE/NETL Case 12 reference. The results are shown comparing the energy demand for PCC, the incremental fuel requirement, and the net higher heating value (HHV) efficiency of the PC power plant integrated with the PCC plant. A comparison of the capital costs for each PCC plant configuration corresponding to a net 550 MWe power generation is also presented. Lastly, a cost of electricity (COE) and cost of CO 2 captured assessment is shown illustrating the substantial cost reductions achieved with the Linde-BASF PCC plant utilizing the advanced SIH configuration in combination with BASF’s OASE® blue solvent technology as compared to the DOE/NETL Case 12 reference. The key factors contributing to the reduction of COE and the cost of CO 2 captured, along with quantification of the magnitude of the reductions achieved by each of these factors, are also discussed. Additionally, a high-level techno-economic analysis of one more highly advanced Linde-BASF PCC configuration case (LB1-CREB) is also presented to demonstrate the significant impact of innovative PCC plant process design improvements on further reducing COE and cost of CO 2 captured for overall plant cost and performance comparison purposes. Overall, the net efficiency of the integrated 550 MWe supercritical PC power plant with CO 2 capture is increased from 28.4% with the DOE/NETL Case 12 reference to 30.9% with the Linde-BASF PCC plant previously presented utilizing the BASF OASE® blue solvent [Ref. 6], and is further increased to 31.4% using Linde-BASF PCC plant with BASF OASE® blue solvent and an advanced SIH configuration. The Linde-BASF PCC plant incorporating the BASF OASE® blue solvent also results in significantly lower overall capital costs, thereby reducing the COE and cost of CO 2 captured from $147.25/MWh and $56.49/MT CO 2 , respectively, for the reference DOE/NETL Case 12 plant, to $128.49/MWh and $41.85/MT CO 2 for process case LB1, respectively, and $126.65/MWh and $40.66/MT CO 2 for process case SIH, respectively. With additional innovative Linde-BASF PCC process configuration improvements, the COE and cost of CO2 captured can be further reduced to $125.51/MWh and $39.90/MT CO 2 for LB1-CREB. Most notably, the Linde-BASF process options presented here have already demonstrated the potential to lower the cost of CO2 captured below the DOE target of $40/MT CO 2 at the 550 MWe scale for second generation PCC technologies.

20 FOSSIL-FUELED POWER PLANTS↗

Characterization of the Uranium-Bearing Products of the Ammonia Injection Remediation Method - 15707

Commissioned in 1943, the Hanford Site is the home of the world's first full-scale nuclear production facility. After more than 40 years of processing, improper discharge and failed storage tanks have resulted in the contamination of the Hanford vadose zone with radiologically contaminated waste. Remediation methods considered for the area include sequestration of the mobile uranium by injection of ammonia (NH{sub 3}) gas. Injection of this reactive gas increases the pore water pH, promoting the dissolution of soil minerals; a subsequent re-establishment of natural conditions is believed to result in the recrystallization of those soil minerals and the co-precipitation of uranium phases. Characterization of the precipitates formed when the ammonia (NH{sub 3}) injection method, proposed for the Hanford vadose zone, is applied to synthetic pore water on a laboratory scale. This involves: Identification of the uranium-bearing phases; A study of the effect of pore water constitution on the phases produced. Synthetic pore water solutions were prepared to mimic selected constituents of the pore water from the Hanford 200 Area: Varying constituent concentrations; Observations led to modifications to increase relative uranium content: 200 ppm → 500 ppm U; Substituted sodium silicate for silicic acid. Dried samples were analyzed by SEM-EDS and powder X-ray diffraction (XRD). A suspension in ethanol was prepared and analyzed by Transmission Electron Microscope (TEM) with Selective Area Electron Diffraction (SAED). Studied and applied geochemical modeling software (Visual Minteq and Geo-chemist's Workbench) for speciation and prediction of potential species formed. Discussion/Conclusions: - SEM-EDS Analysis: Scanning electron microscopy revealed areas of high average atomic weight which would be confirmed to be uranium-rich by energy dispersive spectroscopy analysis. The samples prepared with a more than two-fold increase in uranium concentration did not show the crystal-like forms spotted in the lower [U] samples. The samples prepared with silicic acid substituted in showed little to no uranium precipitation [not shown]. - Diffraction Analysis: Powder XRD analysis of the low [U] samples confirmed the presence of repeating crystalline patterns in the precipitate samples. Nitratine (NaNO{sub 3}) showed a massive presence, potentially obscuring significant peaks in the low [U] samples. Though it was not predicted, there was a tentative match for cejkaite (Na{sub 4}(UO{sub 2})(CO{sub 3}){sub 3}). Elevated uranium samples showed none of the needle-like crystalline shapes [diffraction data pending]. TEM - SAED confirmed the presence of a polycrystalline phase. Data analysis is complex but a comparison of average d-spacings with those of reference minerals could further support a match with cejkaite. - Geochemical Modeling: The system conditions were input into geochemical modeling software to predict the forms present throughout the reaction process. The model was used to demonstrate that the pH of the system could be increased using [NH{sub 4}{sup +}] in solution to better represent the system. Using the [NH{sub 4}{sup +}] increase, the changes in significant species with increasing pH were observed. Future Work: Repeat samples will be prepared and analyzed to bolster any potential identification. Promising samples will be submitted to the PNNL's EMSL facility, which is better equipped to perform the required analyses. An expanded study on the impact of carbonate (CO{sub 3}{sup 2-}) on the phases produced.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Time-lapse seismic inversion for CO 2 saturation with SeisCO2Net: An application to Frio-II site

Seismic monitoring of geological CO 2 storage (GCS) involves highly nonlinear seismic inversion and petrophysical inversion, making it challenging to estimate CO 2 volume efficiently and detect possible early CO 2 leakages. Deep learning (DL) using convolutional neural networks (CNNs) has shown promise in solving highly nonlinear seismic inversion problems. However, direct estimation of CO 2 plume extent/saturation from time-lapse seismic gathers using DL is still underexplored, with no reported field applications to date. The investigation of field data is primarily hindered by scarcity of field data for neural network training. Other obstacles include highly nonlinear seismic-petrophysics inverse relationship, and presence of noise in field seismic data. We introduce SeisCO2Net, a deep CNN that predicts CO 2 saturation maps directly from time-lapse full waveform shot gathers. For training, we use site-specific geological information, fluid flow physics, rock physics, and seismic modeling to generate synthetic datasets that closely resemble the CO 2 storage site. Synthetic tests show promising results, inspiring us to apply SeisCO2Net's trained weights on field data collected at Frio-II GCS site by leveraging transfer learning principles. As reference, we compare SeisCO2Net's predicted CO 2 saturation maps with results obtained from physics-based inversion. Our analyses show both methods display similar CO 2 plume shapes, reasonable CO 2 plume characteristics, and comparable saturation values. Our results suggest pre-training CNNs on physics-informed synthetic datasets and then applying the learned weights to field data is a viable approach to estimating field CO 2 saturation. This method effectively addresses the scarcity of field training data, thus encouraging the feasibility of long-term GCS monitoring.

58 GEOSCIENCES↗

Appalachian Basin Play Fairway Analysis: Thermal Quality Analysis in Low-Temperature Geothermal Play Fairway Analysis (GPFA-AB)

This collection of files are part of a larger dataset uploaded in support of Low Temperature Geothermal Play Fairway Analysis for the Appalachian Basin (GPFA-AB). Phase 1 of the GPFA-AB project identified potential Geothermal Play Fairways within the Appalachian basin of Pennsylvania, West Virginia and New York. This was accomplished through analysis of 4 key criteria: thermal quality, natural reservoir productivity, risk of seismicity, and heat utilization. Each of these analyses represent a distinct project task, with the fifth task encompassing combination of the 4 risks factors. Supporting data for all five tasks has been uploaded into the Geothermal Data Repository node of the National Geothermal Data System (NGDS). This submission comprises the data for Thermal Quality Analysis (project task 1) and includes all of the necessary shapefiles, rasters, datasets, code, and references to code repositories that were used to create the thermal resource and risk factor maps as part of the GPFA-AB project. The identified Geothermal Play Fairways are also provided with the larger dataset. Figures (.png) are provided as examples of the shapefiles and rasters. The regional standardized 1 square km grid used in the project is also provided as points (cell centers), polygons, and as a raster. Two ArcGIS toolboxes are available: 1) RegionalGridModels.tbx for creating resource and risk factor maps on the standardized grid, and 2) ThermalRiskFactorModels.tbx for use in making the thermal resource maps and cross sections. These toolboxes contain item description documentation for each model within the toolbox, and for the toolbox itself. This submission also contains three R scripts: 1) AddNewSeisFields.R to add seismic risk data to attribute tables of seismic risk, 2) StratifiedKrigingInterpolation.R for the interpolations used in the thermal resource analysis, and 3) LeaveOneOutCrossValidation.R for the cross validations used in the thermal interpolations. Some file descriptions make reference to various 'memos'. These are contained within the final report submitted October 16, 2015. Each zipped file in the submission contains an 'about' document describing the full Thermal Quality Analysis content available, along with key sources, authors, citation, use guidelines, and assumptions, with the specific file(s) contained within the .zip file highlighted. UPDATE: Newer version of the Thermal Quality Analysis has been added here: https://gdr.openei.org/submissions/879 (Also linked below) Newer version of the Combined Risk Factor Analysis has been added here: https://gdr.openei.org/submissions/880 (Also linked below)

15 GEOTHERMAL ENERGY↗

OrbNet Denali: A machine learning potential for biological and organic chemistry with semi-empirical cost and DFT accuracy

We present OrbNet Denali, a machine learning model for an electronic structure that is designed as a drop-in replacement for groundstate density functional theory (DFT) energy calculations. The model is a message-passing graph neural network that uses symmetryadapted atomic orbital features from a low-cost quantum calculation to predict the energy of a molecule. OrbNet Denali is trained on a vast dataset of 2.3 × 10 6 DFT calculations on molecules and geometries. This dataset covers the most common elements in biochemistry and organic chemistry (H, Li, B, C, N, O, F, Na, Mg, Si, P, S, Cl, K, Ca, Br, and I) and charged molecules. OrbNet Denali is demonstrated on several well-established benchmark datasets, and we find that it provides accuracy that is on par with modern DFT methods while offering a speedup of up to three orders of magnitude. For the GMTKN55 benchmark set, OrbNet Denali achieves WTMAD-1 and WTMAD-2 scores of 7.19 and 9.84, on par with modern DFT functionals. For several GMTKN55 subsets, which contain chemical problems that are not present in the training set, OrbNet Denali produces a mean absolute error comparable to those of DFT methods. For the Hutchison conformer benchmark set, OrbNet Denali has a median correlation coefficient of R 2 = 0.90 compared to the reference DLPNO-CCSD(T) calculation and R 2 = 0.97 compared to the method used to generate the training data (ωB97X-D3/def2- TZVP), exceeding the performance of any other method with a similar cost. Similarly, the model reaches chemical accuracy for noncovalent interactions in the S66x10 dataset. For torsional profiles, OrbNet Denali reproduces the torsion profiles of ωB97X-D3/def2-TZVP with an average mean absolute error of 0.12 kcal/mol for the potential energy surfaces of the diverse fragments in the TorsionNet500 dataset.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Heat deposition analysis for the High Flux Isotope Reactor’s HEU and LEU core models

The High Flux Isotope Reactor at Oak Ridge National Laboratory is an 85 MW th pressurized light-water-cooled and -moderated flux-trap type research reactor. The reactor is used to conduct numerous experiments, advancing various scientific and engineering disciplines. As part of an ongoing program sponsored by the US Department of Energy National Nuclear Security Administration Office of Material Management and Minimization, studies are being performed to assess the feasibility of converting the reactor’s highly enriched uranium fuel to low-enriched uranium fuel. To support this conversion project, reference models with representative experiment target loading and explicit fuel plate representation were developed and benchmarked for both fuels to (1) allow for consistent comparison between designs for both fuel types and (2) assess the potential impact of low-enriched uranium conversion. These high-fidelity models were used to conduct heat deposition analyses at the beginning and end of the reactor cycle and are presented herein. This article (1) discusses the High Flux Isotope Reactor models developed to facilitate detailed heat deposition analyses of the reactor’s highly enriched and low-enriched uranium cores, (2) examines the computational approach for performing heat deposition analysis, which includes a discussion on the methodology for calculating the amount of energy released per fission, heating rates, power and volumetric heating rates, and (3) provides results calculated throughout various regions of the highly enriched and low-enriched uranium core at the beginning and end of the reactor cycle. These are the first detailed high-fidelity heat deposition analyses for the High Flux Isotope Reactor’s highly enriched and low-enriched core models with explicit fuel plate representation. Lastly, these analyses are used to compare heat distributions obtained for both fuel designs at the beginning and end of the reactor cycle, and they are essential for enabling comprehensive thermal hydraulics and safety analyses that require detailed estimates of the heat source within all of the reactor’s fuel element regions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗