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 685 records · Page 38

Simulation Planning Using Component Based Cost Model

Successful simulations for scientific discovery on high-performance computing platforms require careful planning, including verification of specific application configuration and runtime parameters, estimation of resource requirements, and steering and monitoring of the simulation. However, simulation planning is an aspect of scientific computing that is extremely sparse in available literature or training. In this paper we focus on the resource management aspect of such planning through formulation of a component-based cost model. We illustrate the methodology and formulation through FLASH, a highly configurable simulation code used in multiple scientific domains.

Dubey, Anshu↗

Photonuclear Physics in SCALE

Radiation shielding modeling involves a wide range of applications, such as power reactors, critical experiments, portable sources, linear accelerators, and fusion systems. Particle energy ranges vary from thermal energies to fast energy spectra. MAVRIC, as the radiation shielding sequence of the SCALE code system [1] developed and maintained by Oak Ridge National Laboratory, is capable of modeling both neutron and photon interactions in matter and is widely used by the US Nuclear Regulatory Commission, Department of Energy, the National Nuclear Security Administration, and radiation shielding practitioners all around the world. Although most applications don’t involve highly energetic photons that can cause photonuclear reactions and create coupled chains of high energy, ionized electrons and bremsstrahlung radiation within a material medium, when present, these reactions can have a significant impact upon the downstream particle fluxes from incident photons above 20 MeV. For example, linear accelerators used for medical isotope production are examples of applications in which high-energetic photons are generated and used for driving or controlling other components of a complex facility. Radiation shielding analyses of such a facility would require accounting for any and all significant sources of radiation. Unfortunately, MAVRIC included no data and methods to model either photonuclear reactions or the bremsstrahlung photons until this work, which enabled accounting for both. Photonuclear reaction cross sections were imported from the Evaluated Nuclear Data File libraries, and corresponding methods were implemented in the SCALE code system; see Section 2 for how the nuclear data are imported and used. Bremsstrahlung photons were modeled by assuming thick-target approximation whereby all generated electrons and positrons due to photon interactions are assumed to be absorbed within the same material in which they were born, yielding a quick particle transport simulation method by using fast lookup tables for generating subsequent gammas; see Section 4 for details. A suite of verification models was used to test implemented methods and associated data. Tests yielded satisfactory results for the new features. Enhancements will be released as part of the next major SCALE release (v7.0), which is anticipated to be deployed in 2025, with beta releases available in the interim.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SCALE Modeling of the Fast Spectrum Heat Pipe Reactor

As part of the severe accident analysis collaboration with Sandia National Laboratories (SNL) and the Nuclear Regulatory Commission (NRC), SCALE models were developed for a fast-spectrum heat pipe reactor. These models were based on the Idaho National Laboratory (INL) Design A concept, which is an alternative design to the Los Alamos National Laboratory (LANL) Special Purpose Reactor (SPR), also known as the Megapower reactor. The model contains 1,134 heat pipes, surrounded by hexagonal fuel elements, with a potassium working fluid; the fuel is UO 2 with 19.75 wt% 235 U enrichment. The model contains axial beryllium oxide (BeO) reflectors above and below the active fuel region along with a radial alumina reflector containing 12 B 4 C control drums. The center of the core is left unfueled to make room for two shutdown control rods, one annular and one solid. The active region of the core was discretized into twenty axial and five radial zones to analyze spatial variations in power and burnup. Infinite lattice unit cell sensitivity studies were used to perform verification between the SCALE and INL models. The eigenvalue results agreed well with the reported results to within roughly 50 percent mille (pcm). Full-core model verification was performed by analyzing system eigenvalues with differing configurations of control drum and shutdown rod positions. These full core results all had eigenvalue differences less than 310 pcm. Control drum and shutdown rod worths were also compared, with differences of 3.2% or less. Using the verified model, the isotopic inventory and decay heat, as well as temperature feedback coefficients, were calculated and provided to SNL as input to the MELCOR severe accident code to analyze potential releases from this class of reactor. The results of the MELCOR analysis are provided in a different report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A software environment for effective reliability management for pulsed power design

The reliable design of magnetically insulated transmission lines (MITLs) for very high current pulsed power machines must be accomplished in the future by utilizing a variety of sophisticated modeling tools. The complexity of the models required is high and the number of sub-models and approximations large. The potential for significant analyst error using a single tool is large, with possible reliability issues associated with the plasma modeling tools themselves or the chosen approach by the analyst to solve a given problem. Here, we report on a software infrastructure design that provides a workable framework for building self-consistent models and constraining feedback to limit analyst error. The framework and associated tools aid the development of physical intuition, the development of increasingly sophisticated models, and the comparison of performance results. The work lays the computational foundation for designing state-of-the-art pulsed-power experiments. The design and useful features of this environment are described. We discuss the utility of the Git source code management system and a GitLab interface for use in project management that extends beyond software development tasks.

42 ENGINEERING↗

Verification and application of resonance broadened quasi-linear (RBQ) model with multiple Alfvénic instabilities

The resonance broadened quasilinear (RBQ) model for the problem of relaxing the hot ion distribution function in constant-of-motion 3D space [Gorelenkov et al., Nucl. Fusion 58, 082016 (2018)] is presented with the self-consistent evolution of multiple Alfvén eigenmode amplitudes. The RBQ model represents the generalization of the earlier published model [Berk et al., Nucl. Fusion 35, 1661 (1995)] by carefully examining the wave particle interaction in the presence of realistic Alfvén eigenmode (AE) structures and pitch angle scattering with the help of the guiding center code ORBIT. One aspect of the generalization is that the RBQ model goes beyond the local perturbative-pendulumlike approximation for the wave particle dynamics near the resonance. An iterative procedure is introduced to account for eigenstructures varying within the resonances. It is found that a radially localized mode structure implies a saturation level 2–3 times smaller than that predicted by an earlier bump-on-tail quasilinear model that employed uniform mode structures. We apply the RBQ code to a DIII-D plasma with an elevated q-profile where the beam ion profiles exhibit stiff transport properties [Collins et al., Phys. Rev. Lett. 116, 095001 (2016)]. Finally, the properties of AE driven fast ion distribution relaxation are studied for validations of the applied RBQ model in DIII-D discharges. Initial results show that the model is robust, is numerically efficient, and can predict fast ion relaxation in present and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of R-matrix calculations for charged-particle reactions in the resolved resonance region for the 7 Be system

R-matrix theory is used to present nuclear reactions in the resolved resonance region. It uses information on bound states and low energy resonances to accurately parametrize cross sections on the resonances as well as the non-resonant background. Since the seminal work of Lane and Thomas (1958), the method has been widely used to analyze experimental cross-section data in a broad range of fields spanning nuclear reaction dynamics, nuclear astrophysics, ion beam analysis and their applications. Different R-matrix codes have been developed and used in these different applications with very little communication among the developers or practitioners on the capabilities, achievements or limitations of the codes. A limited comparison among three R-matrix codes on neutron-induced reactions was performed by the International Atomic Energy Agency (IAEA) International Evaluation of Neutron Cross Section Standards project (2007). Since then, significant progress has been made in their implementation of the R-matrix algorithms, and R-matrix codes have enhanced capabilities. In this paper we introduce, for the first time, the results of a comprehensive effort to verify the most widely used R-matrix codes in the various fields of nuclear science and applications: AMUR, AZURE2, CONRAD, EDA, FRESCO, GECCCOS, and SAMMY. In addition to the description of the capabilities of the codes and their specifications, we discuss the results of a joint exercise which was coordinated by the International Atomic Energy Agency. The aim of the exercise was to compare calculations of charged-particle reaction cross sections for the light composite system 7Be. The calculations were performed by the codes using identical input R-matrix parameters and other specifications and were limited to charged-particle channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Vision and Development of a Design, Implementation, and Verification Automation (DIVA) Software Platform for DNA Construction

Abstract DNA construction, while a prerequisite to many biological endeavors, is often a time-consuming distraction from an individual’s primary research objectives. We envisioned that with the right software infrastructure and cultural mindset, a single person could execute in parallel the batched DNA construction tasks of an entire research institute, at scales realizing efficiency gains through process and laboratory automation. In pursuit of this vision, we developed the Design, Implementation, and Verification Automation (DIVA) software platform. DIVA’s web interface enables researchers to design DNA constructs (using visual biological computer-aided design tools and biological parts repositories), submit designs for construction to dedicated staff, and track DNA construction as it progresses. DIVA supports the dedicated staff through the DNA construction process and records both successful and unsuccessful attempts toward improving the overall process. The platform is publicly available at public-diva.jbei.org and its open-source code through github.com/JBEI/DIVA.

Plahar, Hector [DOE Agile BioFoundry , , ,; DOE Jo↗

Validation of electron backscattering in the MCNP6.3 code

A new electron backscattering validation suite has been developed for the MCNP® code. The calculations in this suite cover both condensed history and single-event transport methods, seven different elemental materials ranging from beryllium to uranium, and incident electron energies from 200 eV to 14.1 MeV. In general, the accuracy of the MCNP code for electron backscattering calculations is energy-dependent. For incident electron energies below 2 keV, the condensed history method is not applied while the single-event method shows poor agreement with experimental measurements due to poor accuracy of the underlying atomic data. For incident energies between 2 keV and 256 keV, the single-event method typically shows better agreement with experiments, while for incident energies above 256 keV the condensed history method gives better accuracy. Additionally, agreement with experimental data is generally worse for low- Z materials such as beryllium and carbon, which may reflect deficiencies in the data or physics underlying electron transport through light elements in the MCNP code. The new electron backscattering validation suite, along with the previously introduced electron stopping power and Lockwood energy deposition validation suites, provide MCNP users with a quantified understanding of the uncertainties in electron transport calculations. Future improvements to electron transport in the MCNP code are necessary for low energies (≤ 2 keV) and for low- Z materials.

36 MATERIALS SCIENCE↗

North Dakota Integrated Carbon Storage Complex Feasibility Study. Final report

In spring 2017, the Energy & Environmental Research Center (EERC) initiated an effort to determine the feasibility of developing a commercial-scale CO 2 geologic storage complex able to store 50+ million tonnes (Mt) of CO 2 in central North Dakota safely, permanently, and economically. The objective was to fulfill the goals of the U.S. Department of Energy (DOE) Carbon Storage Assurance Facility Enterprise (CarbonSAFE) Initiative and address technical and nontechnical challenges specific to commercial-scale deployment of a CO2 storage project. The findings clearly show that the concept of capturing CO 2 from a lignite-fired electrical generation facility in central North Dakota and safely and permanently storing the CO 2 in the deep subsurface is indeed technically, economically, and socially feasible. This project evaluated two study areas and their respective geologic storage complexes located adjacent to separate coal-fired facilities in North Dakota: the Basin Electric Power Cooperative (BEPC)-owned Great Plains Synfuels Plant (GPSP) and the Minnkota Power Cooperative (Minnkota)-owned Milton R. Young Station (MRYS). These locations, one with CO 2 capture in place and an existing CO 2 pipeline, are bolstered by progressive North Dakota pore space ownership and long-term liability laws. These elements and a motivated team created an ideal synergistic scenario for ensuring success of the CarbonSAFE Initiative and promoting North Dakota’s statewide vision for carbon management. The project included drilling two new geologic characterization wells, integrating an existing 3-D seismic survey, creating a geologic model subsequently used for injection simulation, a risk assessment, public outreach, and generating a site development plan based on results. In addition, the performance of select National Risk Assessment Partnership tools was evaluated. The geologic characterization wells were drilled ~5600 feet deep to the Broom Creek Formation; ~350 feet of core was retrieved from each well. The core included the Broom Creek (targeted injection zone) and a portion of the overlying Opeche Shale (seal). The Flemmer-1 well, west of Beulah, North Dakota, yielded 169 feet of sandstone. The BNI-1 well located south of Center, North Dakota, yielded 124 feet of sandstone. In each case, laboratory analysis of the sandstone showed permeability in the 300–1000-mD range, with porosity of 20%–30%. The Flemmer-1 well was sited within the boundaries of an existing 3-D seismic survey. Colocating the well with the seismic survey maximized the relationship between new and legacy data and developed a first-of-its-kind interpretation of the geologic fabric of the Broom Creek. Geologic characterization data were integrated into a 5544-mi 2 geocellular model that encompassed both new wells and stratigraphy from the surface to the Amsden Formation (underlying the Broom Creek). The model was later expanded vertically to include the deeper Black Island–Deadwood interval to examine its potential viability as a storage target. The geocellular model provided the foundation for dynamic simulation of CO 2 into the Broom Creek. Results of the simulation suggest that the Broom Creek could accept the DOE target rate of 2 Mt/yr of CO 2 into as few as two wells. To bracket the expected capture from MRYS, simulations were also investigated for a 4-Mt/yr rate near MRYS. Although more wells are needed (two additional), the Broom Creek still has the storage resource to accept the CO 2 at the increased rate. A risk assessment exercise was conducted to identify and assess technical and nontechnical risks that could prevent potential candidate storage complexes within the study area from serving as commercial storage sites. The assessment identified and evaluated six technical risk categories: 1) CO 2 injectivity, 2) storage capacity, 3) lateral migration of CO 2 , 4) lateral pressure propagation, 5) vertical migration of CO 2 or formation brine, and 6) induced seismicity. Following two rounds of analysis and scoring, no risks were determined to preclude continued efforts to develop carbon capture, utilization, and storage (CCUS) in central North Dakota. The risk assessment results will be used to guide future site characterization, modeling and simulation, and monitoring activities. A specific nontechnical strategic risk based on challenges that may be realized in amalgamation of pore space resulted in vertically expanding the geologic model to incorporate the potential for stacked storage in multiple saline reservoirs. By including the Black Island–Deadwood interval (the basal sedimentary reservoir in this region), the amalgamated areal extent could be reduced by as much as 45%. Working with a smaller geographic area reduces risks and costs associated with monitoring and pore space leasing. An economic evaluation incorporating capture; transport (<5 mi); Class VI wells; permitting; and monitoring, verification, and accounting suggests implementing commercial-scale CCUS is economically attractive if the federal tax benefits of 45Q are included. This is validated by Minnkota’s continued pursuit of CO 2 capture and geologic storage at MRYS through its Project Tundra initiative, indicating that there is a business case for CCUS in central North Dakota. Currently, North Dakota is the only state with underground injection control (UIC) Class VI primacy. Built into the North Dakota Century Code is a series of regulatory requirements that guide the process to obtain a Class VI CO 2 storage facility permit. As part of this project, a site development plan was compiled to assure compliance with North Dakota’s requirements to permit a commercial-scale CO 2 storage operation and includes a prospective time line encompassing a general breakdown of activities. In total, an estimated 30 months would be needed to execute the necessary steps to attain a North Dakota CO 2 storage facility permit. Outreach was an integral part of the project and encompassed any project-related activity that had contact or exposure beyond the project team. The goals of outreach were to foster an environment from which stakeholders could make informed decisions about the project and gauge community receptiveness to a CCUS project. A consistent set of messages and outreach products were developed in conjunction with an outreach advisory board that integrated project partners and team members. To gauge public acceptability of geologically storing CO 2 , 5611 households in the project area were invited to participate in an online survey. The survey results indicate that the public attitude regarding CCUS is neutral to positive, with strong sentiment that CO 2 capture and storage may be an approach to maintain the economic vitality of the region. To achieve project objectives, critical support in the form of financial backing, engineering evaluations, site access, outreach collaboration, operations data, risk assessment/evaluation, and software access was provided by BEPC, the North Dakota Industrial Commission Lignite Research Council, ALLETE Clean Energy, BNI Energy, North American Coal Corporation, Minnkota, Prairie Public Broadcasting, Computer Modelling Group Ltd., and Schlumberger.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

CLimate Impact: Determining Etiology thRough pAthways (CLDERA)

Climate impacts have broad economic, health, political, and national security ramifications. Societally relevant impacts are typically farther downstream, are the product of multiple interacting processes, and can arise over small regions and timeframes because their sources are short-term and localized. Short-term forcings (as can be seen in volcanic eruptions, climatic tipping points (e.g., the collapse of rainforests or the disappearance of sea ice), or in increasingly plausible climate interventions) fundamentally possess low signal-to-noise and could benefit from accounting for the multiple conditional processes through which a downstream impact arises. Under the Grand Challenge LDRD CLDERA (CLimate impacts: Discovering Etiology thRough pAthways), we have developed tools to enable downstream impact attribution from geographically and temporally localized source forcings in the climate. CLDERA developed methods that can distinguish how a localized source drives the climate system to respond with particular impacts. The how is embodied in pathways – the spatio-temporally evolving chain of physical processes that connects a source to a series of increasingly distant impacts. Novel analytic methods in pursuit of downstream impact attribution were developed and demonstrated on simulations and observations of the 1991 eruption of Mt. Pinatubo in the Philippines. As described within this report we have • developed stratospheric expertise and aerosol modeling capabilities in E3SM, • created original methods to detect and model pathways from source-to-impact, and • advanced climate attribution through novel methods, cases, and approaches. Further, CLDERA developed a tiered verification process consisting of controlled datasets to prototype, verify, and refine the original method development. CLDERA increased Sandia’s footprint in the climate analytics community and developed new climate collaborations whilst also creating a cadre of climate analysts at Sandia. The products from CLDERA have been extensive with a total of 9 journal articles published, 12 articles submitted and under review, and an additional 8 articles in preparation. We have produced 1750 simulated years and developed 9 code-bases. This report details these accomplishments and serves as a summary of the work completed during the CLDERA Grand Challenge.

54 ENVIRONMENTAL SCIENCES↗

Monte Carlo calculation and verification of the geometrical factors for the NPDGamma experiment

The NPDGamma experiment measures the parity-violating asymmetry in gamma-ray emission in the capture of polarized neutrons on liquid parahydrogen. The sensitivity to the asymmetry for each detector in the array is used as a parameter in the extraction of the physics asymmetry from the measured data. The detector array is approximately cylindrically symmetric around the target and a step-wise sinusoidal function has been used for the sensitivity in the previous iteration of the NPDGamma experiment, but deviations from cylindrical symmetry necessitate the use of a Monte Carlo model to determine corrections to the geometrical factors. For the calculations, source code modifications to MCNPX were done in order to calculate the sensitivity of each cesium iodide detector to the physics asymmetry. We describe the MCNPX model and results from calculations and how the results are validated through measurement of the parity violating asymmetry of gamma-rays from neutron capture on chlorine.

Grammer, K. B.↗

The Source Physics Experiment (SPE) Science Plan (Version 2.4c)

The Source Physics Experiment (SPE) series is a long-term NNSA research and development effort designed to improve U.S. arms control and nuclear nonproliferation verification and monitoring capabilities. The findings from the SPE will advance the United States’ nuclear explosion monitoring capabilities, particularly with respect to detection, discrimination and determination of yields associated with small nuclear explosions that can be lost amid the noisy seismo-acoustic background from other sources. The data generated from the SPE, a series of well-designed and recorded chemical explosions, will contribute to the development and validation of first-principles explosive source generated seismo-acoustic modeling codes. These codes will then facilitate the update of semi-empirical methods, currently based on historic test site data, such that key explosion observables can be reproduced, thus improving confidence in nuclear test monitoring in new areas and/or under novel emplacement conditions. The overall SPE project is comprised of both the development of the new explosion simulation codes and the chemical explosion test series. The chemical explosion test series will generate the empirical data required to both develop and validate the new simulation codes.

58 GEOSCIENCES↗

Data Evaluation of Actinide Cross Sections: 230 Th and 231 Th

A new evaluation of the ENDL cross section set for Thorium (Z=90) is developed using the TALYS statistical model cross section code. The primary goal of this effort is to produce an evaluation that attempts to match as closely as possible fission cross sections developed through surrogate reaction techniques on actinide targets 230 Th and 231 Th. This evaluation effort and the processing needed to render its results into data libraries is a necessary step in making the efforts of nuclear experimentalists useful to the broad community of researchers engaged in simulations of nuclear fusion for basic and applied science. Another aspect, verification and validation against various AGEX experiments, is also presented. The end-product is an updated library that includes the latest measurements of fission cross sections for comparison against those measured via traditional techniques. All the steps in the evaluation, processing, validation and verification, and library release are described in the following sections. For completeness, the appendix contains all the parameters used in the TALYS cross section evaluation for neutrons incident on 230 Th and 231 Th.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Source Physics Experiment (SPE) Science Plan

The Source Physics Experiment (SPE) series is a long-term NNSA research and development effort designed to improve U.S. arms control and nuclear nonproliferation verification and monitoring capabilities. The findings from the SPE will advance the United States’ nuclear explosion monitoring capabilities, particularly with respect to detection, discrimination and determination of yields associated with small nuclear explosions that can be lost amid the noisy seismo-acoustic background from other sources. The data generated from the SPE, a series of well-designed and recorded chemical explosions, will contribute to the development and validation of first-principles explosive source generated seismo-acoustic modeling codes. These codes will then facilitate the update of semi-empirical methods, currently based on historic test site data, such that key explosion observables can be reproduced, thus improving confidence in nuclear test monitoring in new areas and/or under novel emplacement conditions. The overall SPE project is comprised of both the development of the new explosion simulation codes and the chemical explosion test series. The chemical explosion test series will generate the empirical data required to both develop and validate the new simulation codes.

58 GEOSCIENCES↗

Standards, dissemination, and best practices in systems biology

In this study, the reproducibility of scientific research is crucial to the success of the scientific method. Here, we review the current best practices when publishing mechanistic models in systems biology. We recommend, where possible, to use software engineering strategies such as testing, verification, validation, documentation, versioning, iterative development, and continuous integration. In addition, adhering to the Findable, Accessible, Interoperable, and Reusable modeling principles allows other scientists to collaborate and build off of each other’s work. Existing standards such as Systems Biology Markup Language, CellML, or Simulation Experiment Description Markup Language can greatly improve the likelihood that a published model is reproducible, especially if such models are deposited in well-established model repositories. Where models are published in executable programming languages, the source code and their data should be published as open-source in public code repositories together with any documentation and testing code. For complex models, we recommend container-based solutions where any software dependencies and the run-time context can be easily replicated.

59 BASIC BIOLOGICAL SCIENCES↗

Fission Product Transport in TRISO Particles and Pebbles

This document demonstrates completion of the goals described in the technical narrative of the FOA project titled: ”Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing” regarding fission product transport in the Kairos-proposed fuel pebble by INL and Kairos Power. Showcased in this report are code developments and simulations in BISON that extend the state of the art in computation and understanding of fission product transport in a TRISO fuel particle and pebble. These enhancements lay the foundation for making predictions of fission product transport that can be used as input in the fuel licensing process. This was achieved by installing existing fuel material models originally used in PARFUME, developing a new failure probability method that is efficient and multi-dimensional, employing material homogenization, and expanding verification and validation simulations to demonstrate the efficacy of the work. All this work is leveraged to spotlight the main deliverable; a three-dimensional model and corresponding demonstration simulation of a pebble, which will serve as the starting point for models used to predict fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Python Script to Read MCNP6.3 Surface-Source Files

This report provides a Python script to read an MCNP ® surface source file created with the SSW card with SYM = 0 (i.e., the default symmetry treatment). For background: the general format of an MCNP surface-source file is described in; however, that document did not provide coding and/or a tool to interrogate such files. The current format will not be given in this document other than through the record-read statements necessary for the script to function. The reader capability in this report is augmented with the ability to directly write a couple demonstrative outputs: 1. A comma-separated value (CSV) file containing particle phase-space state information and 2. A Matplotlib histogram of the energy distribution of the particles. This report also describes accompanying verification work that shows the script performing as required with MCNP6.2, MCNP6.3, and (expected) MCNP6.4 surface-source files. However, users of the enclosed script must still verify that the script is behaving correctly for their own work.

97 MATHEMATICS AND COMPUTING↗