Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code verification”

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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Automatic Testcase Generation for Flight Software

The TacSat3 project is applying Integrated Systems Health Management (ISHM) technologies to an Air Force spacecraft for operational evaluation in space. The experiment will demonstrate the effectiveness and cost of ISHM and vehicle systems management (VSM) technologies through onboard operation for extended periods. We present two approaches to automatic testcase generation for ISHM: 1) A blackbox approach that views the system as a blackbox, and uses a grammar-based specification of the system's inputs to automatically generate *all* inputs that satisfy the specifications (up to prespecified limits); these inputs are then used to exercise the system. 2) A whitebox approach that performs analysis and testcase generation directly on a representation of the internal behaviour of the system under test. The enabling technologies for both these approaches are model checking and symbolic execution, as implemented in the Ames' Java PathFinder (JPF) tool suite. Model checking is an automated technique for software verification. Unlike simulation and testing which check only some of the system executions and therefore may miss errors, model checking exhaustively explores all possible executions. Symbolic execution evaluates programs with symbolic rather than concrete values and represents variable values as symbolic expressions. We are applying the blackbox approach to generating input scripts for the Spacecraft Command Language (SCL) from Interface and Control Systems. SCL is an embedded interpreter for controlling spacecraft systems. TacSat3 will be using SCL as the controller for its ISHM systems. We translated the SCL grammar into a program that outputs scripts conforming to the grammars. Running JPF on this program generates all legal input scripts up to a prespecified size. Script generation can also be targeted to specific parts of the grammar of interest to the developers. These scripts are then fed to the SCL Executive. ICS's in-house coverage tools will be run to measure code coverage. Because the scripts exercise all parts of the grammar, we expect them to provide high code coverage. This blackbox approach is suitable for systems for which we do not have access to the source code. We are applying whitebox test generation to the Spacecraft Health INference Engine (SHINE) that is part of the ISHM system. In TacSat3, SHINE will execute an on-board knowledge base for fault detection and diagnosis. SHINE converts its knowledge base into optimized C code which runs onboard TacSat3. SHINE can translate its rules into an intermediate representation (Java) suitable for analysis with JPF. JPF will analyze SHINE's Java output using symbolic execution, producing testcases that can provide either complete or directed coverage of the code. Automatically generated test suites can provide full code coverage and be quickly regenerated when code changes. Because our tools analyze executable code, they fully cover the delivered code, not just models of the code. This approach also provides a way to generate tests that exercise specific sections of code under specific preconditions. This capability gives us more focused testing of specific sections of code.

Bushnell, David Henry↗
NASA NTRS2005

Compression-Loaded Composite Panels With Elastic Edge Restraints and Initial Prestress

A parametric study of the effects of test-fixture-induced initial prestress and elastic edge restraints on the prebuckling and buckling responses of a compression-loaded, quasi-isotropic curved panel is presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code with high-fidelity models. The results presented show that a wide range of prebuckling and buckling behavior can be obtained by varying parameters that represent circumferential loaded-edge restraint and rotational unloaded-edge restraint provided by a test fixture and that represent the mismatch in specimen and test-fixture radii of curvature. For a certain range of parameters, the panels exhibit substantial nonlinear prebuckling deformations that yield buckling loads nearly twice the corresponding buckling load predicted by a traditional linear bifurcation buckling analysis for shallow curved panels. In contrast, the results show another range of parameters exist for which the nonlinear prebuckling deformations either do not exist or are relatively benign, and the panels exhibit buckling loads that are nearly equal to the corresponding linear bifurcation buckling load. Overall, the results should be of particular interest to scientists, engineers, and designers involved in simulating flight-hardware boundary conditions in structural verification and certification tests, involved in validating structural analysis tools, and interested in tailoring buckling performance.

Hilburger, Mark W.↗
NASA NTRS2004

Effects of Elastic Edge Restraints and Initial Prestress on the Buckling Response of Compression-Loaded Composite Panels

A parametric study of the effects of test-fixture-induced initial prestress and elastic edge restraints on the prebuckling and buckling responses of a compression-loaded, quasi-isotropic curved panel is presented. The numerical results were obtained by using a geometrically nonlinear finite element analysis code with high-fidelity models. The results presented show that a wide range of prebuckling and buckling behavior can be obtained by varying parameters that represent circumferential loaded-edge restraint and rotational unloaded-edge restraint provided by a test fixture and that represent the mismatch in specimen and test-fixture radii of curvature. For a certain range of parameters, the panels exhibit substantial nonlinear prebuckling deformations that yield buckling loads nearly twice the corresponding buckling load predicted by a traditional linear bifurcation buckling analysis for shallow curved panels. In contrast, the results show another range of parameters exist for which the nonlinear prebuckling deformations either do not exist or are relatively benign, and the panels exhibit buckling loads that are nearly equal to the corresponding linear bifurcation buckling load. Overall, the results should also be of particular interest to scientists, engineers, and designers involved in simulating flight-hardware boundary conditions in structural verification and certification tests, involved in validating structural analysis tools, and interested in tailoring buckling performance.

Hilburger, Mark W.↗
DOE OSTI2024

Structural complexity of snapshots of two-dimensional Fermi-Hubbard systems

The development of quantum gas microscopy for two-dimensional optical lattices has provided an unparalleled tool to study the Fermi-Hubbard model (FHM) with ultracold atoms. Spin-resolved projective measurements, or snapshots, have played a significant role in quantifying correlation functions, theory verification, and thus the uncovering of underlying physical phenomena such as antiferromagnetism at commensurate filling on bipartite lattices and other charge and spin correlations, as well as dynamical properties at various densities. Here we employ a recent concept, the multiscale structural complexity, and show that when computed for the snapshots (of single spin species, local moments, or total density) it can provide a theory-free property, immediately accessible to experiments. Specifically, after benchmarking results for Ising and $XY$ models, we study the structural complexity of snapshots of the repulsive FHM in the two-dimensional square lattice as a function of doping and temperature. We generate projective measurements using determinant quantum Monte Carlo and compare their complexities against those from the experiment. We demonstrate that these complexities are linked to relevant physical observables such as the entropy and double occupancy. Their behaviors capture the development of correlations and relevant length scales in the system. Furthermore, we provide an open-source code in python which can be implemented into data analysis routines in experimental settings for the square lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Center for Extended Magnetohydrodrynamic Modeling

This work used the NIMROD code to solve two tokamak regimes found on the DIII-D tokamak. The first regime studied is the Quiescent H-mode regime that has the advantage of enabling tokamak high performance without the deleterious effects of Edge Localized Modes (ELMs). ELMs are a possible impediment to the success of tokamaks as a burning plasma reactor and the QH mode regime provides a mechanism for avoiding it. The QH mode achieves this state due to the appearance of long wavelength instability that saturates instead of reaching an ELM state. NIMROD simulations were able to describe many of the features of a particular discharge, including the importance of shear flow in reaching the state, the role of the parallel heat flux to decorrelate the temperature and velocity fluctations, and in the correlation of density and velocity fluctations to provide the density pumpout. This work resulted in two direct publications and an invited talk. Progress towards these publications resulted in two more publications (verification and solving the Grad-Shafranov equation with scrape-off-layer currents). The second regime studied was the Giant Sawtooth regime whereby energetic particles affect the stability of the 1/1 mode. This work resulted in performance improvements of NIMROD's energetic particle algorithms, and improvements to the drift-kinetic equation infrastructure used by both continuum and energetic particle closure methods. The studies proved more computationally demanding than expected because as the 1/1 is stabilized, the near marginal conditions cause long simulations times. Publication of this work is in progress under other funding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Synthesizing Safety Conditions for Code Certification Using Meta-Level Programming

In code certification the code consumer publishes a safety policy and the code producer generates a proof that the produced code is in compliance with the published safety policy. In this paper, a novel viewpoint approach towards an implementational re-use oriented framework for code certification is taken. It adopts ingredients from Necula's approach for proof-carrying code, but in this work safety properties can be analyzed on a higher code level than assembly language instructions. It consists of three parts: (1) The specification language is extended to include generic pre-conditions that shall ensure safety at all states that can be reached during program execution. Actual safety requirements can be expressed by providing domain-specific definitions for the generic predicates which act as interface to the environment. (2) The Floyd-Hoare inductive assertion method is refined to obtain proof rules that allow the derivation of the proof obligations in terms of the generic safety predicates. (3) A meta-interpreter is designed and experimentally implemented that enables automatic synthesis of proof obligations for submitted programs by applying the modified Floyd-Hoare rules. The proof obligations have two separate conjuncts, one for functional correctness and another for the generic safety obligations. Proof of the generic obligations, having provided the actual safety definitions as context, ensures domain-specific safety of program execution in a particular environment and is simpler than full program verification.

Eusterbrock, Jutta↗
NASA NTRS2009

The Integrated Hazard Analysis Integrator

Hazard analysis addresses hazards that arise in the design, development, manufacturing, construction, facilities, transportation, operations and disposal activities associated with hardware, software, maintenance, operations and environments. An integrated hazard is an event or condition that is caused by or controlled by multiple systems, elements, or subsystems. Integrated hazard analysis (IHA) is especially daunting and ambitious for large, complex systems such as NASA s Constellation program which incorporates program, systems and element components that impact others (International Space Station, public, International Partners, etc.). An appropriate IHA should identify all hazards, causes, controls and verifications used to mitigate the risk of catastrophic loss of crew, vehicle and/or mission. Unfortunately, in the current age of increased technology dependence, there is the tendency to sometimes overlook the necessary and sufficient qualifications of the integrator, that is, the person/team that identifies the parts, analyzes the architectural structure, aligns the analysis with the program plan and then communicates/coordinates with large and small components, each contributing necessary hardware, software and/or information to prevent catastrophic loss. As viewed from both Challenger and Columbia accidents, lack of appropriate communication, management errors and lack of resources dedicated to safety were cited as major contributors to these fatalities. From the accident reports, it would appear that the organizational impact of managers, integrators and safety personnel contributes more significantly to mission success and mission failure than purely technological components. If this is so, then organizations who sincerely desire mission success must put as much effort in selecting managers and integrators as they do when designing the hardware, writing the software code and analyzing competitive proposals. This paper will discuss the necessary and sufficient requirements of one of the significant contributors to mission success, the IHA integrator. Discussions will be provided to describe both the mindset required as well as deleterious assumptions/behaviors to avoid when integrating within a large scale system.

Morris, A. Terry↗
NASA NTRS2006

Partial Automation of Requirements Tracing

Requirements Tracing on Target (RETRO) is software for after-the-fact tracing of textual requirements to support independent verification and validation of software. RETRO applies one of three user-selectable information-retrieval techniques: (1) term frequency/inverse document frequency (TF/IDF) vector retrieval, (2) TF/IDF vector retrieval with simple thesaurus, or (3) keyword extraction. One component of RETRO is the graphical user interface (GUI) for use in initiating a requirements-tracing project (a pair of artifacts to be traced to each other, such as a requirements spec and a design spec). Once the artifacts have been specified and the IR technique chosen, another component constructs a representation of the artifact elements and stores it on disk. Next, the IR technique is used to produce a first list of candidate links (potential matches between the two artifact levels). This list, encoded in Extensible Markup Language (XML), is optionally processed by a filtering component designed to make the list somewhat smaller without sacrificing accuracy. Through the GUI, the user examines a number of links and returns decisions (yes, these are links; no, these are not links). Coded in XML, these decisions are provided to a "feedback processor" component that prepares the data for the next application of the IR technique. The feedback reduces the incidence of erroneous candidate links. Unlike related prior software, RETRO does not require the user to assign keywords, and automatically builds a document index.

Hayes, Jane↗
NASA NTRS2018

NOS3: NASA Operational Simulator for Small Satellites

The NASA Operational Simulator for Small Satellites (NOS3) is a suite of open-source software tools to aid in areas such as software development, integration & test (I&T), mission operations/training, verification and validation (V&V), and software systems check-out. NOS3 provides a software development environment, a multi-target build system, operational interface/ground software, dynamics and environment simulations, and software-based hardware models. NOS3 has just recently been open-sourced by NASA and is available for immediate use. It enables the development of flight software (FSW) early in the project life cycle when hardware availability is limited. Small satellite development suffers from extensive lead times on many of the commercial-off-the-shelf (COTS) components as well as limited funding for engineering test units (ETUs). To alleviate the need to provide a hardware test-bed for each developer/tester, NOS3 hardware models are based upon characteristic data or manufacturer's data sheets for each individual component. The NOS3 hardware models' fidelity is such that FSW executes unaware that physical hardware is not present. This allows FSW binaries to be compiled for both the simulation environment and the flight computer without changing the FSW source code. For hardware models that provide data which is dependent upon the environment and spacecraft dynamics, such as a GPS receiver or magnetometer, an open-source tool from NASA GSFC (42 Spacecraft Simulator) is used to provide the necessary data. The underlying infrastructure used to transfer messages between FSW and the hardware models can also be used to monitor, intercept, and inject messages, which has proven to be beneficial for V&V of larger missions such as James Webb Space Telescope (JWST). As hardware is selected and becomes available, drivers can be added to the NOS3 environment to enable hardware-in-the-loop (HWIL) testing. When strict time synchronization is not vital, any number of combinations of hardware components and software-based models can be tested. NOS3 was actively used for FSW development and component testing of the Simulation-to-Flight 1 (STF-1) CubeSat and the Lunar IceCube CubeSat. As NOS3 matures, hardware models have been added for common small satellite components such as GPS receivers, electrical power systems and batteries, and antenna systems.

Suder, Mark↗
NASA NTRS1993

Three-dimensional Navier-Stokes analysis and redesign of an imbedded bellmouth nozzle in a turbine cascade inlet section

Verification of proposed turbopump blading performance will involve evaluation of candidate blades in cascade test facilities. It is necessary to be able to predict the flow fields within these cascades for the results to be applicable to actual engine environments. This work presents the results of a study to predict the flow field for the NASA Lewis Transonic Turbine Blade Cascade Facility, which is similar to those used to evaluate rocket propulsion turbines. A pitchwise nonuniform total pressure distribution was observed at the blade row leading edge plane. A CFD analysis was used to show that the cause of the flow nonuniformity was a pair of vortices that originated in an embedded bellmouth inlet. Further CFD analysis was used to verify that a redesigned inlet section resulted in a flow with acceptable uniformity. A computational analysis was chosen because physical accessibility to the inlet section was limited, and because a computational approach also allows one to examine design changes cheaper and more quickly than an experimental approach would. The PARC code, a general purpose, three-dimensional, Navier-Stokes code with multiblock solution capability, was chosen for the present study. Results are presented detailing the computational requirements needed to accurately predict flows of this nature. Calculations of the original geometry showed total pressure loss regions consistent in strength and in location to experimental measurements. An examination of the results shows that the distortions are caused by a pair of vortices that originate as a result of the interaction of the flow with the imbedded bellmouth. Computations were performed for an inlet geometry which eliminated the imbedded bellmouth by bridging the region between it and the upstream wall. This analysis indicated that eliminating the imbedded bellmouth eliminates the troublesome pair of vortices, resulting in a flow with much greater pitchwise uniformity.

Giel, P. W.↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Calibration of TOMS Radiances From Ground Observations

Verification of a stratospheric ozone recovery remains a high priority for environmental research and policy definition. Models predict an ozone recovery at a much lower rate than the measured depletion rate observed to date. Therefore improved precision of the satellite and ground ozone observing systems are required over the long term to verify its recovery. We show that validation of radiances from the ground can be a very effective means for correcting long term drifts of backscatter type satellite measurements and can be used to cross calibrate all BUV instruments in orbit (TOMS, SBUV/2, GOME, SCIAMACHY, OMI, GOME-2, OMPS). This method bypasses the retrieval algorithms used to derive ozone products from both satellite and ground based measurements that are normally used to validate the satellite data. Radiance comparisons employ forward models, but they are inherently more accurate than the retrieval This method employs very accurate comparisons between ground based zenith sicy radiances and satellite nadir radiances and employs two well established capabilities at the Goddard Space Flight Center, 1) the SSBUV calibration facilities and 2) the radiative transfer codes used for the TOMS and SBUV/2 algorithms and their subsequent refinements. The zenith sky observations are made by the SSBUV where its calibration is maintained to a high degree of accuracy and precision. Radiative transfer calculations show that ground based zenith sky and satellite nadir backscatter ultraviolet comparisons can be made very accurately under certain viewing conditions. Initial ground observations taken from Goddard Space Flight Center compared with radiative transfer calculations has indicated the feasibility of this method. The effect of aerosols and varying ozone amounts are considered in the model simulations and the theoretical comparisons. The radiative transfer simulations show that the ground and satellite radiance comparisons can be made with an uncertainty of less than l\% without the knowledge of the amount ozone viewed by either instrument on ground or in space. algorithms.

Bojkov, B. R.↗
DOE OSTI2023

Earthquake-induced impact of base-isolated buildings: theory, numerical modeling, and design solutions

Earthquake shaking more intense than that used to size the horizontal clearance between a base–isolated building and near–rigid perimeter moat wall will result in hard impact, producing high–frequency, high–amplitude acceleration response in the structure and supported equipment. Herein this paper provides a design solution for the damaging effects of hard impact by installing a compliant engineered element in the load path between the base–isolated building and the moat wall, resulting in soft impact and a much smaller acceleration response. The engineered element assumed herein is a commercial–off–the–shelf marine fender with mechanical properties determined by physical testing. The attachment of a flexible engineered element, with well–defined stiffness and damping, to a near–rigid moat wall, simplifies the numerical modeling of the building–moat wall system and eliminates the need to bound the lateral stiffness of the wall for impact calculations. The simple model of the engineered element can be implemented in commercial finite element codes. Theory is developed for two–sided impact of a single–degree–of–freedom oscillator. Analytical solutions are derived for the shifted first–mode frequency of the impacted oscillator and for its free–vibration response. The shifted first–mode frequency is a function of the composite lateral stiffness of the isolator–engineered element assembly and its earthquake–induced displacement. Local peaks in the spectral response of the impacted oscillator form at odd integer multiples of the shifted first–mode frequency. The analytical solutions can be used to verify, in part, the numerical model used for impact analysis.

42 ENGINEERING↗
DOE OSTI2017

Galaxy bias from galaxy–galaxy lensing in the DES science verification data

We present a measurement of galaxy-galaxy lensing around a magnitude-limited ($i_{AB} < 22.5$) sample of galaxies from the Dark Energy Survey Science Verification (DES-SV) data. We split these lenses into three photometric-redshift bins from 0.2 to 0.8, and determine the product of the galaxy bias $b$ and cross-correlation coefficient between the galaxy and dark matter overdensity fields $r$ in each bin, using scales above 4 Mpc/$h$ comoving, where we find the linear bias model to be valid given our current uncertainties. We compare our galaxy bias results from galaxy-galaxy lensing with those obtained from galaxy clustering (Crocce et al. 2016) and CMB lensing (Giannantonio et al. 2016) for the same sample of galaxies, and find our measurements to be in good agreement with those in Crocce et al. (2016), while, in the lowest redshift bin ($z\sim0.3$), they show some tension with the findings in Giannantonio et al. (2016). We measure $b\cdot r$ to be $0.87\pm 0.11$, $1.12 \pm 0.16$ and $1.24\pm 0.23$, respectively for the three redshift bins of width $\Delta z = 0.2$ in the range $0.2<z <0.8$, defined with the photometric-redshift algorithm BPZ. Using a different code to split the lens sample, TPZ, leads to changes in the measured biases at the 10-20\% level, but it does not alter the main conclusion of this work: when comparing with Crocce et al. (2016) we do not find strong evidence for a cross-correlation parameter significantly below one in this galaxy sample, except possibly at the lowest redshift bin ($z\sim 0.3$), where we find $r = 0.71 \pm 0.11$ when using TPZ, and $0.83 \pm 0.12$ with BPZ.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Extension and Demonstration of NEAMS Multiphysics Tools to Lead-Cooled, Sodium-Cooled, and Molten Salt Fast Reactor Applications

The SHARP toolkit is a high-fidelity reactor simulation tool developed under the U.S. Department of Energy, Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. SHARP toolkit is comprised of the neutronics module PROTEUS thermal hydraulics module Nek5000, and structural mechanics module Diablo. During FY17 and FY18, the PROTEUS and Nek5000 components of SHARP were applied to solve challenging sodium-cooled fast reactor (SFR) problems. In particular, selected hot channel factors (HCF) for a prototype metal-fueled SFR design (the AFR-100) were analyzed in high fidelity, and the “SHARP zooming capability” for SFRs was developed and demonstrated to reduce computational expense for full core problems in cases where detailed data is needed in selected fuel assemblies. After the previous success applying SHARP to challenging SFR problems, the focus in FY19 expanded to additional fast reactor applications: lead cooled fast reactors (LFR), sodium cooled fast reactors (SFR), and fast molten salt reactors (MSR). The specific technical tasks were (1) assessment of hot channel factors for LFR, for which no data currently exists, (2) implementation of a gamma transport capability in the PROTEUS solvers to extend the zooming capability to non-fueled assemblies, (3) demonstration of zooming capability in assemblies of the Versatile Test Reactor (VTR), and (4) development of a new coupled capability for simulating fast MSRs, specifically modeling precursor transport inside and outside the core. The progress on tasks (1) and (3) are discussed primarily in this document. Tasks (2) and (4) are briefly mentioned in this document, and the interested reader is referred to separate reports for detailed information. First-of-a-kind hot channel factor (HCF) estimation for LFR with high fidelity codes (PROTEUS/Nek5000) was successfully demonstrated in this study. Selected HCF were computed and compared with SFR data (AFR-100, EBR-II). The findings confirm that different reactor types, design parameters and uncertainties lead to different HCFs. Careful estimation of HCF for a specific design is necessary to obtain appropriate HCFs. In addition to improvement in HCF accuracy, high fidelity tools generate data to help the designer better understand the mechanism of the impact from these uncertainties. For example, the impact of cladding thickness manufacturing tolerance resulted in non-intuitive effects in the corner pins of the LFR assembly. This procedure of computing HCF using high fidelity models shows promise and flexibility for being repeated for any arbitrary reactor of choice, since these tools are targeted to solve many reactor types and geometries. Progress was made towards extending the previously demonstrated SHARP zooming capability to non-fueled SFR assemblies. In particular, a gamma transport capability was implemented in both high fidelity PROTEUS solvers in order to accurately account for heat deposition caused by gamma particles, which accounts for ~10% of total core power. Neutronics verification cases were carried out for a candidate Versatile Test Reactor (VTR) design using the new gamma transport capability in PROTEUS. Comparisons were made with continuous energy MCNP calculations and shown to agree well. The models for the full core design with heterogeneous control and fuel assemblies is in progress and will be completed in FY20. In the fast molten salt reactor (MSR) domain, a new capability was developed to model the flow of neutron precursors both inside and outside of the core, which impacts the effective delayed neutron fraction of the system. The PROTEUS-SN and PROTEUS-NODAL solvers were coupled to Nek5000 in order to simulate precursor flow. The workflows were demonstrated on test problems. The fast MSR modeling capability was supported under two other work packages, and comprehensive details are included in a separate report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Performance Results on CPU/GPU Exascale Architectures for OMEGA: The Ocean Model for E3SM Global Applications

The US Department of Energy (DOE) conducts climate simulations on some of the world’s largest supercomputers. These exascale machines use heterogeneous architectures with both CPUs and GPUs, and scientific codes must adapt to make full use of this computing power. Los Alamos National Lab is developing Omega: The Ocean Model for E3SM Global Applications, which is specifically designed for modern exascale computers. It uses external libraries that have been optimized for a variety of architectures to run on different supercomputers. Omega is an unstructured-mesh ocean model based on TRiSK numerical methods. It will be the new ocean component of the DOE’s Energy Exascale Earth System Model (E3SM). The algorithms in Omega follow those of the current ocean component, MPAS-Ocean, but it will be written in C++ rather than Fortran to take advantage of the Kokkos performance portability library. Omega spatial operators are written as Kokkos kernels to run efficiently on both CPUs and GPUs. Work on Omega began in 2023 with a new C++ framework for unstructured mesh partitioning, halo exchanges, parallel IO, and Kokkos interfaces. The current version, Omega-0, is being developed to solve the shallow water equations and at present includes all of the tendency terms but not time stepping. Here we share the results of Omega-0 verification and performance testing. Verification includes unit tests implemented with CTest as well as convergence tests in Polaris, an in-house python package with a large suite of test problems. Performance tests compare simulations conducted on CPUs versus GPUs and across different architectures: tests are run on Frontier, which has AMD “Optimized 3rd Gen EPYC” CPUs and AMD MI250X GPUs, as well as Perlmutter, which is composed of AMD EPYC 7763 CPUs and NVIDIA A100 GPUs.

58 GEOSCIENCES↗