Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SCALE Code System”

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

Preconceptual Design of Irradiated Fuel Salt Management System

The National Reactor Innovation Center was established by the U.S. Department of Energy to accelerate the demonstration and deployment of advanced nuclear reactors. To meet this mission, the National Reactor Innovation Center has developed the Laboratory for Operations and Testing in the U.S. (LOTUS) to support the first fast-spectrum molten salt reactor demonstration. The goal of this research is to understand the system requirements that may be applied to a molten salt reactor experiment in the LOTUS test bed and to perform systems analyses and preconceptual development of a storage container for the management of the irradiated molten salt fuel. This study proposes a set of irradiated fuel salt management processes for the molten salt reactor, from defueling to storage or disposal. For that, the historic Molten Salt Reactor Experiment (MSRE) at Oak Ridge National Laboratory has served as a crucial source of data for a mock application. This paper proposes a potential geometry for storing irradiated fuel salt based on MSRE data and the Monte Carlo N-Particle radiation transport code. The geometry was analyzed for criticality safety under abnormal situations, such as water ingress and interaction between containers. Based on the determined configuration of the storage system, the decay heat produced from the irradiated fuel salts was calculated to provide a reasonable time scale of the dormancy period. Using the same geometry assumed above, the neutron and gamma dose rates from the irradiated fuel salt at the end of the dormancy period were estimated, which will be useful for understanding the handling and radiation protection requirements for handling the irradiated fuel salt. Lastly, a reactor in the LOTUS test bed could use chloride-based fuel as opposed to the fluoride-based fuel used in MSRE. Thus, a literature study provided data to understand the similarities and differences between the characteristics of fluoride- and chloride-based salts in terms of salt chemistry, material characteristics, postirradiation behavior, neutronics, thermal hydraulics, and heat removal. This study, together with the proposed irradiated fuel salt management processes, preliminary criticality calculations, decay heat, and dose rate estimates, provides a strong basis for future investigation into the management of irradiated fuel salts generated during molten salt reactor technology demonstrations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1998

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗
NASA NTRS2012

Experimenting with the GMAO 4D Data Assimilation

The Global Modeling and Assimilation Office (GMAO) has been working to promote its prototype four-dimensional variational (4DVAR) system to a version that can be exercised at operationally desirable configurations. Beyond a general circulation model (GeM) and an analysis system, traditional 4DV AR requires availability of tangent linear (TL) and adjoint (AD) models of the corresponding GeM. The GMAO prototype 4DVAR uses the finite-volume-based GEOS GeM and the Grid-point Statistical Interpolation (GSI) system for the first two, and TL and AD models derived ITom an early version of the finite-volume hydrodynamics that is scientifically equivalent to the present GEOS nonlinear GeM but computationally rather outdated. Specifically, the TL and AD models hydrodynamics uses a simple (I-dimensional) latitudinal MPI domain decomposition, which has consequent low scalability and prevents the prototype 4DV AR ITom being used in realistic applications. In the near future, GMAO will be upgrading its operational GEOS GCM (and assimilation system) to use a cubed-sphere-based hydrodynamics. This versions of the dynamics scales to thousands of processes and has led to a decision to re-derive the TL and AD models for this more modern dynamics, thus taking advantage of a two-dimensional MPI decomposition and improved scalability properties. With the aid of the Transformation of Algorithms in FORTRAN (l'AF) automatic adjoint generation tool and some hand-coding, a version of the cubed-sphere-based TL and AD models, with a simplified vertical diffusion scheme, is now available, enabling multiple configurations of standard implementations of 4DV AR in GEOS. Concurrent to this development, collaboration with the National Centers for Environmental Prediction (NCEP) and the Earth System Research Laboratory (ESRL) has allowed GMAO to implement a hybrid-ensemble capability within the GEOS data assimilation system. Both 3Dand 4D-ensemble capabilities are presently available thus allowing GMAO to now evaluate the performance and benefit of various ensemble and variational assimilation strategies. This presentation will cover the most recent developments taking place at GMAO and show results from various comparisons from traditional techniques to more recent ensemble-based ones.

Todling, R.↗
DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The EGS Collab Project: Learnings from Experiment 1

The primary objective of the EGS Collab Project sponsored by DOE is to increase the understanding needed to efficiently implement enhanced geothermal systems (EGS). One goal of the EGS Collab project is to create a collaborative research environment in which to study stimulation of crystalline rock at the 10 meter scale. Key to this effort is the collection of high quality data to allow comparison to numerical coupled process models in an effort to build confidence in the codes and modeling techniques used. In response to this, the EGS Collab team has created an underground test bed at the Sanford Underground Research Facility (SURF) in Lead SD at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently designing a second test bed aimed at investigating shear stimulation (Experiment 2). At the Experiment 1 location, we have characterized our host rock using laboratory testing and numerous field-based geophysical and geological techniques, and created a well-instrumented test bed to allow us to carefully monitor stimulation events and flow tests. In addition to the installed geophysical sensors, we have used tracer tests, differences in the ambient microbial communities at flow collection locations, and cold water injection to inform us about dynamic flow pathways. In Experiment 1, we have hydraulically stimulated the host rock in a number times at several locations in one well, creating new fractures that connect to existing fractures between the injection and production boreholes. We have performed long-term ambient and chilled water injection tests as an analog to EGS, and have monitored system changes resulting from these water injections through geophysical monitoring, flow and pressure measurements, tracer tests, and microbiology. Here, we summarize the tests performed, issues identified including poroelastic and thermoelastic effects, Joule-Thomson effects, restarting effects, indications of flow channeling, and the primary learnings from Experiment 1.

Enhanced Geothermal Systems, EGS Collab, stimulati↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
NASA NTRS1973

Technical consulting services for earth resources program

Evaluations of several sets of test data which show a normal amount of variability in gamma and Dmax on KODAK Aerial Duplicating Film (ESTAR Base) SO-467 are reported. Kodak is now furnishing to engineers at GSFC a full set of test data for each coating of SO-467 film. Glass-enclosed sensitometer step tablets which are being made for the Herrnfeld instrument at GSFC are discussed. These tablets include a gray scale and uniformity check areas as well as tribar resolution patterns imaged in a durable chromium layer. It is recommended that the Sioux Falls EDC procure four new enlarging printers and study the value and economics of a system for printer automation based on edge coding.

Source record↗
NASA NTRS1992

Modal Test/Analysis Correlation of Space Station Structures Using Nonlinear Sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.↗
NASA NTRS1992

Modal test/analysis correlation of Space Station structures using nonlinear sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlations. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.↗
NASA NTRS1995

Introduction to Penelope

A formal program verification is a (mathematical) proof that a program executed according to its intended model meets some specification. This proves that the algorithm defined by the program is correct in the precise technical sense of being consistent with a particular specification. A program correct in this sense is free from a large and important class of errors, even though its behavior may still produce unintended results--either because the implementation of the programming language itself does not match the model of execution, or because the specification does not correctly express the user's intentions. Penelope is a prototype system for interactively developing and verifying programs that are written in a rich subset of sequential Ada. Penelope can be used to develop a program and its correctness proof incrementally, and in concert with one another. Incrementality is used in a number of ways to help make verification more tractable and more productive. For example, if an already-verified program is modified, one can attempt to prove the modified version by replaying and modifying the original verification. Penelope's specification language, Larch/Ada, belongs to the family of Larch interface languages. Larch/Ada scales up properly, in the sense that it is demonstrably sound to decompose a system hierarchically and reason locally about the implementation of each piece. Penelope has been applied in various demonstration projects--for specification (guidance control, distributed operating systems), verification (of off-the-shelf code), and formal development (by non-expert as well as expert users). Some features of Penelope have been embodied in Ada Wise, a lint-like non-interactive tool that warns of the potential for certain dynamic semantic errors in Ada programs.

Guaspari, David↗
NASA NTRS2003

X-33 Rev-F Turbulent Aeroheating Results From Test 6817 in NASA Langley 20-Inch Mach 6 Air Tunnel and Comparisons With Computations

Measurements and predictions of the X-33 turbulent aeroheating environment have been performed at Mach 6, perfect-gas air conditions. The purpose of this investigation was to compare measured turbulent aeroheating levels on smooth models, models with discrete trips, and models with arrays of bowed panels (which simulate bowed thermal protections system tiles) with each other and with predictions from two Navier-Stokes codes, LAURA and GASP. The wind tunnel testing was conducted at free stream Reynolds numbers based on length of 1.8 x 10(exp 6) to 6.1 x 10(exp 6) on 0.0132 scale X-33 models at a = 40-deg. Turbulent flow was produced by the discrete trips and by the bowed panels at ill but the lowest Reynolds number, but turbulent flow on the smooth model was produced only at the highest Reynolds number. Turbulent aeroheating levels on each of the three model types were measured using global phosphor thermography and were found to agree to within .he estimated uncertainty (plus or minus 15%) of the experiment. Computations were performed at the wind tunnel free stream conditions using both codes. Turbulent aeroheating levels predicted using the LAURA code were generally 5%-10% lower than those from GASP, although both sets of predictions fell within the experimental accuracy of the wind tunnel data.

Hollis, Brian R.↗
NASA NTRS2004

Full Scale Rotor Aeroacoustic Predictions and the Link to Model Scale Rotor Data

The NASA Aeroacoustic Prediction System (NAPS) is used to establish a link between model-scale and full-scale rotor predictions and is partially validated against measured wind tunnel and flight aeroacoustic data. The prediction approach of NAPS couples a comprehensive rotorcraft analysis with acoustic source noise and propagation codes. The comprehensive analysis selected for this study is CAMRAD-II, which provides the performance/trim/wake solution for a given rotor or flight condition. The post-trim capabilities of CAMRAD-II are used to compute high-resolution sectional airloads for the acoustic tone noise analysis, WOPMOD. The tone noise is propagated to observers on the ground with the propagation code, RNM (Rotor Noise Model). Aeroacoustic predictions are made with NAPS for an isolated rotor and compared to results of the second Harmonic Aeroacoustic Rotor Test (HART-II) program, which tested a 40% dynamically and Mach-scaled BO-105 main rotor at the DNW. The NAPS is validated with comparisons for three rotor conditions: a baseline condition and two Higher Harmonic Control (HHC) conditions. To establish a link between model and full-scale rotor predictions, a full-scale BO-105 main rotor input deck for NAPS is created from the 40% scale rotor input deck. The full-scale isolated rotor predictions are then compared to the model predictions. The comparisons include aerodynamic loading, acoustic levels, and acoustic pressure time histories for each of the three conditions. With this link established, full-scale predictions are made for a range of descent flight conditions and compared with measured trends from the recent Rotorcraft Operational Noise Abatement Procedures (RONAP) flight test conducted by DLR and ONERA. Additionally, the effectiveness of two HHC conditions from the HART-II program is demonstrated for the full-scale rotor in flight.

Boyd, D. Douglas, Jr.↗
NASA NTRS2011

Vertical Distribution of Aersols and Water Vapor Using CRISM Limb Observations

Near-infrared spectra taken in a limb-viewing geometry by the Compact Reconnaissance Imaging Spectrometer for Mars (CRISM) on-board the Mars Reconnaissance Orbiter (MRO) provide a useful tool for probing atmospheric structure. Specifically, the observed radiance as a function of wavelength and height above the limb allows the vertical distribution of both dust and ice aerosols to be retrieved. These data serve as an important supplement to the aerosol profiling provided by the MRO/MCS instrument allowing independent validation and giving additional information on particle physical and scattering properties through multi-wavelength studies. A total of at least ten CRISM limb observations have been taken so far covering a full Martian year. Each set of limb observations nominally contains about four dozen scans across the limb giving pole-to-pole coverage for two orbits at roughly 100 and 290 W longitude over the Tharsis and Syrtis/Hellas regions, respectively. At each longitude, limb scans are spaced roughly 10 degrees apart in latitude, with a vertical spatial resolution on the limb of roughly 800 m. Radiative transfer modeling is used to model the observations. We compute synthetic CRISM limb spectra using a discrete-ordinates radiative transfer code that accounts for multiple scattering from aerosols and accounts for spherical geometry of the limb observations by integrating the source functions along curved paths in that coordinate system. Retrieved are 14-point vertical profiles for dust and water ice aerosols with resolution of 0.4 scale heights between one and six scale heights above the surface. After the aerosol retrieval is completed, the abundances of C02 (or surface pressure) and H20 gas are retrieved by matching the depth of absorption bands at 2000 nm for carbon dioxide and at 2600 run for water vapor. In addition to the column abundance of water vapor, limited information on its vertical structure can also be retrieved depending on the signal available from aerosol scattering.

Smith, Michael D.↗
NASA NTRS2000

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗
NASA NTRS1992

Observations on computational methodologies for use in large-scale, gradient-based, multidisciplinary design

Various computational methodologies relevant to large-scale multidisciplinary gradient-based optimization for engineering systems design problems are examined with emphasis on the situation where one or more discipline responses required by the optimized design procedure involve the solution of a system of nonlinear partial differential equations. Such situations occur when advanced CFD codes are applied in a multidisciplinary procedure for optimizing an aerospace vehicle design. A technique for satisfying the multidisciplinary design requirements for gradient information is presented. The technique is shown to permit some leeway in the CFD algorithms which can be used, an expansion to 3D problems, and straightforward use of other computational methodologies.

Newman, P. A.↗
DOE OSTI2020

TeraChem: A graphical processing unit-accelerated electronic structure package for large-scale ab initio molecular dynamics

TeraChem was born in 2008 with the goal of providing fast on-the-fly electronic structure calculations to facilitate ab initio molecular dynamics studies of large biochemical systems such as photoswitchable proteins and multichromophoric antenna complexes. Originally developed for videogaming applications, graphics processing units (GPUs) offered a low-cost parallel computer architecture that became more accessible for general-purpose GPU computing with the release of CUDA in 2007. The evaluation of the electron repulsion integrals (ERIs) is a major bottleneck in electronic structure codes and provides an attractive target for acceleration on GPUs. Thus, highly efficient routines for evaluation of and contractions between the ERIs and density matrices were implemented in TeraChem. Here, electronic structure methods were developed and implemented to leverage these integral contraction routines, resulting in the first quantum chemistry package designed from the ground up for GPUs. This GPU acceleration makes TeraChem capable of performing large-scale ground and excited state calculations in the gas and condensed phase. Today, TeraChem's speed forms the basis for a suite of quantum chemistry applications, including optimization and dynamics of proteins, automated and interactive chemical discovery tools, and large-scale nonadiabatic dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗
DOE OSTI2025

Maps of growing season gross primary production and net ecosystem exchange for Council Road Mile Marker 71, Seward Peninsula, Alaska, [2017-2023]

This data archive is in support of the Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) publication "Integrating Characteristic Arctic Vegetation in a Land Surface Model Improves Representation of Carbon Dynamics Across a Tundra Landscape", by Murphy et al. (2025a). Murphy et al. (2025a) evaluated whether incorporating observed Arctic vegetation heterogeneity into ELM, the land model of the Department of Energy’s Energy Exascale Earth System Model (E3SM), improved simulations of tundra carbon cycling. The associated model archive can be found at Murphy et al. (2025b). The study focused on the spatial patterns and net landscape-level growing season productivity and carbon uptake. As part of this evaluation, observationally derived maps of average growing season (June–August) net ecosystem exchange (NEE) and gross primary production (GPP) were developed for the same domain. These maps, which form the dataset described here, integrate eddy covariance flux tower, remote sensing, and vegetation community data to provide spatially explicit benchmarks for model evaluation. The maps provide spatially explicit estimates of average growing season NEE and GPP across 13 tundra vegetation communities within the study domain. By combining flux tower observations with Airborne Visible-Infrared Imaging Spectrometer-Next Generation (AVIRIS-NG) hyperspectral imagery and drone-based normalized difference vegetation index (NDVI), these maps capture the heterogeneity of carbon fluxes associated with different Arctic vegetation types. While they represent average seasonal conditions rather than interannual variability, the maps provide a unique dataset for evaluating model performance, comparing vegetation community contributions to landscape-scale carbon cycling, and supporting regional analyses of Arctic carbon dynamics. This data archive contains 5 m resolution maps of vegetation communities, vegetation community average growing season GPP, and vegetation community average growing season NEE (three *.tif files), a User’s Guide (*pdf file), and Table 1 of the User’s Guide displaying vegetation community coverage and average growing season NEE and GPP values (*.csv file).

Murphy, Bailey [ORNL] (ORCID:0000000203995221)↗