Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analysis codes”

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 541 records · Page 30

FastrSHWFS Analysis

This code analyzes the reflection geometry of modified Shack-Hartmann Wavefront Sensor (SHWFS) masks, with a focus on the fastrSHWFS designs. It takes a height map of a reflective mask, divides it into sub apertures, removes the global focus term by referencing a central block, and fits local surface planes to each sub aperture. From these fits it computes reflection angles and propagates the corresponding rays downstream to determine spot positions and the overall reflected beam footprint as a function of distance. The implementation is parameterized in sub aperture size, number of sub apertures, focal length, and mask dimensions, so it can be adapted to different mask designs beyond the two current fastrSHWFS masks. The code: Converts mask bitmaps into physical height in microns and reconstructs a focus subtracted mask, divides the mask into a grid of sub apertures and assigns pixel and micron coordinates to each block, excludes user specified unused sub apertures from the analysis, fits a plane to each sub aperture to obtain local surface normals, propagates reflected rays to a range of z positions to compute ray spots and the full reflected beam width, computes the effective mask tilt angle with respect to the incoming beam as a function of propagation distance, and provides helper routines for plotting the grid over the mask and for visualizing the beam geometry. These tools are intended for iterating on fastrSHWFS mask designs and for planning downstream optical layouts, for example choosing lens positions and apertures that capture the reflected beam given known focal plane distances and beam widths.

Gerard, BenjaminL [Lawrence Livermore National Lab↗

ALFITeX: A code introducing improvements in the analysis of alpha spectra

A new code for the deconvolution of complex alpha-particle spectra, called ALFITeX, has been developed. The procedure is based on the Levenberg-Marquardt algorithm, with the curve fitting being the mathematical function formed by the convolution of a Gaussian with two left-handed exponential curves in the low-energy-tail region. ALFITeX has been checked with its application to the calculation of the alpha-particle emission probabilities of several nuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

EFFECT OF DISCRETIZATION ERROR ON THE QUANTIFICATION OF MODEL UNCERTAINTY IN TWO-PHASE FLOW SIMULATIONS

Verification, validation and uncertainty quantification (VVUQ) have become a common practice in thermal-hydraulics analysis. In general, these activities deal with propagations of uncertainties in computer code simulations, e.g., through system analysis codes. However, most existing such activities in thermal-hydraulics analysis have been primarily focused on input and model uncertainties, while numerical errors were largely overlooked. Numerical errors can appear in many different forms, e.g., round-off error, statistical sampling error, and discretization error. In thermal-hydraulics analysis, especially two-phase flow simulations commonly encountered in reactor safety analysis, the lack of consideration of discretization error is mainly due to the difficulty in estimating them. Accurate estimations of discretization error require continuous mesh refinement and/or implementation of high-order numerical schemes in system analysis code, both of which are difficult to achieve in existing codes. In this work, we will build a computer code that incorporates both first-order and second-order numerical methods to solve the two-phase flow problems. The first-order method resembles the one used in many existing system analysis codes; and the second-order method works as the reference to estimate numerical errors. Numerical verification of spatial discretization schemes will be presented in the form of mesh convergence study. It will also be demonstrated via case studies that, in practical scenarios, discretization errors can be as large as, or even larger than, model uncertainties.

42 ENGINEERING↗

Modeling of a Generic Pebble Bed High-temperature Gas-cooled Reactor (PB-HTGR) with SAM

This report presents the modeling of the core of a generic pebble-bed reactor (PBR) at the system level using the System Analysis Module (SAM) code. This work is an extension of a previous work by the authors (Ooi et al. (2021)) that used the so-called 2-D ring model approach to model the PBMR-400. With the new approach, the pebble bed of the reactor is modeled with multiple PBCoreChannel components with spherical heat structures which allows the code to calculate thermal fluid parameters with built-in closure relations. The new core-channel approach is an improvement to the 2-D ring model approach as it does not introduce geometric distortions to the model and thus reduces the uncertainties of the predictions. In addition to thermal fluid simulations, point kinetics (PKE) are included to the model. Simulations are performed under a steady-state normal operation condition and a load-following transient scenario. This particular transient scenario is chosen as it tests both the thermal fluid and neutronics aspects of the model. The predicted results from both the steady-state and transient scenarios are compared with the results by Stew- art et al. (2021) who performed similar simulations with a Griffin-Pronghorn coupled tool. Despite the differences between the codes, with SAM being a system-analysis code and Pronghorn being a porous-medium code, both sets of results compare favorably. The over- all profiles and trends of the predicted temperatures and reactivities from the SAM and Griffin/Pronghorn simulations are similar, with some differences in their predicted values. The first part of the report covers the significance of a relatively fast-running approach that is capable of modeling the pebble bed reactor at the system-level while simultaneously capturing the radial thermal behavior of the core. Then, the modeling approach used in this work is discussed in details. Lastly, the results and comparisons with the Griffin/Pronghorn simulation are presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Antiproton bounds on dark matter annihilation from a combined analysis using the DRAGON2 code

Abstract Early studies of the AMS-02 antiproton ratio identified a possible excess over the expected astrophysical background that could be fit by the annihilation of a weakly interacting massive particle (WIMP). However, recent efforts have shown that uncertainties in cosmic-ray propagation, the antiproton production cross-section, and correlated systematic uncertainties in the AMS-02 data, may combine to decrease or eliminate the significance of this feature. We produce an advanced analysis using the DRAGON2 code which, for the first time, simultaneously fits the antiproton ratio along with multiple secondary cosmic-ray flux measurements to constrain astrophysical and nuclear uncertainties. Compared to previous work, our analysis benefits from a combination of: (1) recently released AMS-02 antiproton data, (2) updated nuclear fragmentation cross-section fits, (3) a rigorous Bayesian parameter space scan that constrains cosmic-ray propagation parameters.We find no statistically significant preference for a dark matter signal and set strong constraints on WIMP annihilation tobb̅, ruling out annihilation at the thermal cross-section for dark matter masses below ∼ 200 GeV. We do find a positive residual that is consistent with previous work, and can be explained by a ∼ 70 GeV WIMP annihilating below the thermal cross-section. However, our default analysis finds this excess to have a local significance of only 2.8σ, which is decreased to 1.8σwhen the look-elsewhere effect is taken into account.

Astronomy & Astrophysics↗

Implementation of Additional Models into the MACCS Code for Nearfield Consequence Analysis

The NRC’s Non-Light Water Reactor Vision and Strategy report discusses the MACCS code readiness for nearfield analyses. To increase the nearfield capabilities of MACCS, the plume meander model from Ramsdell and Fosmire was integrated into MACCS and the MACCS plume meander model based on U.S. NRC Regulatory Guide 1.145 was updated. Test cases were determined to verify the plume meander model implementation into MACCS 4.1. The results using the implemented MACCS plume meander models match the comparisons with other codes and analytical calculations. This verifies that the additional MACCS plume meander models have been successfully implemented into MACCS 4.1. This report documents the verification of these model implementations into MACCS and a comparison of the results using

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Analysis of the ASME Code Rules for Subsection III-5-HHB (Composite Materials) for Current HTR Design Requirements

This document includes the critical analysis review of the American Society of Mechanical Engineers (ASME) Section III Division 5 Subsection HH Subpart B (HHB), including Mandatory Appendices, that was published in 2023. In the context of this document, reference to “the code” is specific to this subsection unless otherwise specified. A specific composites task group within the ASME Nonmetallic Design and Materials Working Group, with the support of external experts, was established to perform a gap analysis review. The significant findings are summarized here. The committee response with suggested action items are detailed in the body of the report.

36 MATERIALS SCIENCE↗

Gamma Spectrometry Code Rodeo for Uranium Enrichment—FY 2021 Report

In FY22, the participating Department of Energy (DOE) laboratories continued the work on the project “Gamma Spectrometry Code Rodeo for Uranium Enrichment”. Oak Ridge National Laboratory (ORNL) and Lawrence Livermore National Laboratory (LLNL) continued their efforts to acquire spectra using uranium sources. A collimator (with side and back shields) was custom designed and constructed at ORNL. Spectra were acquired at ORNL using the uranium enrichment standards supplied by New Brunswick Laboratory (NBL) and the collimated M400 detector. Spectra were also acquired using steel absorbers of different thickness placed between the source and the collimated M400 detector. The plan was to ship the collimator to LLNL so that spectra could be acquired using the LLNL detector in collimated geometries. This activity was suspended since problems were noticed with the performance of the LLNL M400 detector. Also, it was brought to light to the project team that some of the aspects of the M400 detector had been updated by the vendor H3D in response to requests by the IAEA. In February 2021, the ORNL and LLNL M400 detectors were returned to the H3D factory in order to resolve the problems with the LLNL detector and to update ORNL and LLNL detectors to the same configuration as the M400 supplied to the IAEA. H3D repaired the LLNL M400 and performed the IAEA updates on both the ORNL and LLNL detectors. ORNL and LLNL received the detectors in April 2021. Data acquisition was continued using the upgraded detectors. ORNL re-collected the spectra using point sources and the NBL standards in uncollimated and collimated geometries. The collimator was shipped from ORNL to LLNL. Spectra were acquired at LLNL through Q4 of FY22. When the ORNL spectra were examined closely, it became apparent that there were spurious artifacts present in some of the spectra. This rendered suspect the ORNL spectra collected in June/July 2021 time frame. ORNL contacted H3D and per H3D’s advice, the M400 detector was returned to H3D for repairs in mid-August 2021. H3D diagnosed the problem and established that the root cause was the extra tight packaging that created a strain on the CZT crystals, causing one of them to break down. The repaired M400 was returned to ORNL in early September 2021. ORNL re-started data acquisition for the third time. Data was continued to be acquired at ORNL and LLNL through September 2021. The ORNL and LLNL data will be shared with the principal investigators from the analyzing laboratories. The GADRAS code and FRAM have been developed to have the capability to analyze the M400 spectra. Good quality spectra from FY21 have been analyzed using modified versions of GADRAS and FRAM. The spectra have also been analyzed using the GEM code. The analysis results from GADRAS, FRAM and GEM codes are presented and discussed in this report. Some preliminary results from FRAM and GEM analysis were presented in February 2021 at the U-Pu Isotopics Workshop sponsored by the IAEA. CZTU code is being modified to analyze M400 spectra. Data collection, code development and analysis will be continued in FY22 as per the FY22 Project Work Plan (PWP).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification and validation of the PLTEMP/ANL code for thermal hydraulic analysis of experimental and test reactors

The document compiles in a single volume several verification and validation works done for the PLTEMP/ANL code during the years of its development and improvement. Some works that are available in the open literature are simply referenced at the outset, and are not included in the document. PLTEMP has been used in conversion safety analysis reports of several US and foreign research reactors that have been licensed and converted. A list of such reactors is given. Each chapter of the document deals with the verification or validation of a specific model. The model verification is usually done by comparing the code with hand calculation, Microsoft spreadsheet calculation, or Mathematica calculation. The model validation is done by comparing the code with experimental data or a more validated code like the RELAP5 code.

Kalimullah, M.↗

Demonstration of LOTUS multiphysics BEPU analysis framework for LB-LOCA simulations

The LOCA Toolkit for U.S. light water reactors (LOTUS) currently under development at Idaho National Labs (INL) is a plug and play, multiphysics environment to be used in support of system level plant transient analysis. New proposed rule changes in LOCA safety regulations (10 CFR 50.46c) require the inclusion of cladding hydrogen content in the evaluation of equivalent cladding reacted (ECR) and peak cladding temperature (PCT). The potential for a future implementation of said rule change motivates the current development of a unique union of core design, system analysis, and fuel performance codes within a framework allowing for uncertainty quantification (UQ). A demonstration of LOTUS capabilities to address this potential need is performed with an integration of the core design code PHISICS, the fuel performance code FRAPCON, and the system analysis code RELAP5-3D. UQ is performed via Wilks’ method to compute the one sided 95/95 confidence values of the aforementioned safety metrics. Results demonstrate the benefits of the Multiphysics Best Estimate Plus Uncertainty (MP-BEPU) methodology and provide useful visualization of the limiting cases.

10CFR50.46c↗

TRANSP integrated modeling code for interpretive and predictive analysis of tokamak plasmas

This paper provides a comprehensive review of the TRANSP code, a sophisticated tool for interpretive and predictive analysis of tokamak plasmas, detailing its major capabilities and features. It describes the equations for particle, power, and momentum balance analysis, as well as the poloidal field diffusion equations. The paper outlines the spatial and time grids used in TRANSP and details the equilibrium assumptions and solvers. Various models for heating and current drive and radiation, including updates to the NUBEAM model, are discussed. The handling of large-scale events such as sawtooth crashes and pellet injections is examined, along with the predictive capabilities for advancing plasma profiles. The integration of TRANSP with the ITER Integrated Modeling and Analysis Suite (IMAS) is highlighted, demonstrating enhanced data access and analysis capabilities. Additionally, the paper discusses best practices and continuous integration techniques to enhance TRANSP's robustness. The suite of TRANSP tools, designed for efficient data analysis and simulation, further supports the optimization of tokamak operations and coupling with other tokamak codes. Continuous development and support ensure that TRANSP remains a major code for the analysis of experimental data for controlled thermonuclear fusion, maintaining its critical role in supporting the optimization of tokamak operations and advancing fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

LANL CSSE L2: Case Study of In Situ Data Analysis in ASC Integrated Codes [Slides]

The massive, complex datasets generated by scientific simulations on next-generation computer architectures at extreme scale pose new challenges for data management and scientific understanding. A critical element in meeting these new challenges is in situ data analysis where information is reduced and analyzed as close to the point of generation as possible, for example, at run time or as the data stream to longer term storage. This milestone will demonstrate in situ data analysis in xRAGE, an ASC Eulerian Radiation Hydrodynamics code. We will work with the Eulerian Applications Project to identify target problems that will exercise physics of interest to the program, and we will work closely with users to maximize the achieved usability and productivity improvements. We will investigate enhancements to the in situ process that enable and improve techniques of data reduction and data triage in addition to improving in situ performance by using portable acceleration technology. We will report on their viability and performance in xRAGE or another ASC code. These products will provide ASC IC developers with reference codes and tools that will benefit DSW by advancing the technology for predictive science at extreme scales.

97 MATHEMATICS AND COMPUTING↗

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Utilizing many-core accelerators for halo and center finding within a cosmology simulation

Efficiently finding and computing statistics about “halos” (regions of high density) are essential analysis steps for N-body cosmology simulations. However, in state-of-the-art simulation codes, these analysis operators do not currently take advantage of the shared-memory data-parallelism available on multi-core and many-core architectures. The Hybrid / Hardware Accelerated Cosmology Code (HACC) is designed as an MPI+X code, but the analysis operators are parallelized only among MPI ranks, because of the difficulty in porting different X implementations (e.g., OpenMP, CUDA) across all architectures on which it is run. In this paper, we present portable data-parallel algorithms for several variations of halo finding and halo center finding algorithms. These are implemented with the PISTON component of the VTK-m framework, which uses Nvidia's Thrust library to construct data-parallel algorithms that allow a single implementation to be compiled to multiple backends to target a variety of multi-core and many-core architectures. Finally, we compare the performance of our halo and center finding algorithms against the original HACC implementations on the Moonlight, Stampede, and Titan supercomputers. The portability of Thrust allowed the same code to run efficiently on each of these architectures. On Titan, the performance improvements using our code have enabled halo analysis to be performed on a very large data set (81923 particles across 16,384 nodes of Titan) for which analysis using only the existing CPU algorithms was not feasible.

Sewell, Christopher↗

Cross-Code Verification of Neutronics Analysis Tools at INL Applied for 238 Pu Production in the Advanced Test Reactor

Here, analyses are completed for experiments prior to experiment irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Various codes are used to qualify all experiments planned for insertion in the reactor, thereby ensuring that all safety and programmatic requirements are satisfied preirradiation. Among the common experiment analysis tools at INL are MCNP5 coupled to ORIGEN2 (MOPY) and MC21. MOPY uses MCNP5 for transport calculations along with calculations for fluxes and select reaction rates, and then ORIGEN2 handles the step-by-step and postirradiation depletion. MC21 handles all in-reactor transport and step-by-step, during-irradiation, depletion calculations, and then ORIGEN (SCALE 6.2.3) is used for decay and dose calculations postirradiation. The MOPY results, along with those obtained via two variations of the MC21 model, were compared in terms of 238 Pu production in the ATR’s H10 position. For the MOPY model, the MC21 model utilizing the HELIOS-based fission product (FP) library, and the MC21 model utilizing the expanded 1300 FP library, the during-cycle irradiation in-core heating results were sufficiently equivalent; however, the MOPY model and the MC21 model with the HELIOS library showed some differences relating to the respective FP libraries. Ultimately, the MC21 model with a 1300 FP library produced the most consistent results throughout the cycle, whereas the MC21 model that utilized the (smaller) HELIOS library was able to handle during-irradiation analysis but lacked certain short-lived FPs that significantly contributed to the total decay heat at shutdown. MOPY, on the other hand, was found to overpredict fission gas production, as a result of limitations in the ORIGEN2 code.

ATR↗

SE2-ANL Subchannel Analysis with New Heating Calculation Based on VARIANT Transport Code

The multi-assembly, steady-state sub-channel analysis code SE2-ANL has been widely used in thermal-hydraulic analysis of sodium cooled fast reactors (SFR). It is used to compute the core-wide temperature profiles for a given flow distribution and to determine the assembly flow rates yielding a desired (assembly-wise) peak temperature distribution. It is a modified version of SUPERENERGY-2 coupled with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In the current version of SE2-ANL, the neutron and gamma fluxes are calculated using the finite-difference diffusion theory option of the DIF3D code with six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. In order to determine the pin power distributions more accurately by eliminating these limitations, a new coupled heating calculation procedure has recently developed based on the neutron and gamma transport calculations with the VARIANT code. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at the center of each fuel pin and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The SE2-ANL code has been modified to interface with the new, coupled neutron and gamma heating calculation scheme. In addition, several improvements have been made on the computational methods and models of SE2-ANL. To eliminate the limit on the axial mesh size due to the numerical instability problem of the current explicit scheme, the numerical scheme for discretizing the energy equation in the axial direction has been changed to the {theta}-method of the SLTHEN code. An automated flow allocation scheme has also been implemented in SE2-ANL to determine assembly flow rates such that the peak cladding mid-wall temperatures of individual orifice zones are equalized over the burn cycle. The impacts of the new heating calculation scheme based on VARIANT transport solutions on the cladding and fuel temperatures have been tested against the existing method of SE2-ANL. Using the existing and new heating calculation schemes, full-core subchannel analyses were performed for the metal core design of the 1000 MWt Advanced Burner Reactor (ABR-1000), and the resulting coolant, cladding and fuel temperatures were compared. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The Secret Life of Quarks, Final Report for the University of North Carolina at Chapel Hill

This final report summarizes activities and results at the University of North Carolina as part of the the SciDAC-2 Project The Secret Life of Quarks: National Computational Infrastructure for Lattice Quantum Chromodynamics. The overall objective of the project is to construct the software needed to study quantum chromo- dynamics (QCD), the theory of the strong interactions of subatomic physics, and similar strongly coupled gauge theories anticipated to be of importance in the LHC era. It built upon the successful efforts of the SciDAC-1 project National Computational Infrastructure for Lattice Gauge Theory, in which a QCD Applications Programming Interface (QCD API) was developed that enables lat- tice gauge theorists to make effective use of a wide variety of massively parallel computers. In the SciDAC-2 project, optimized versions of the QCD API were being created for the IBM Blue- Gene/L (BG/L) and BlueGene/P (BG/P), the Cray XT3/XT4 and its successors, and clusters based on multi-core processors and Infiniband communications networks. The QCD API is being used to enhance the performance of the major QCD community codes and to create new applications. Software libraries of physics tools have been expanded to contain sharable building blocks for inclusion in application codes, performance analysis and visualization tools, and software for au- tomation of physics work flow. New software tools were designed for managing the large data sets generated in lattice QCD simulations, and for sharing them through the International Lattice Data Grid consortium. As part of the overall project, researchers at UNC were funded through ASCR to work in three general areas. The main thrust has been performance instrumentation and analysis in support of the SciDAC QCD code base as it evolved and as it moved to new computation platforms. In support of the performance activities, performance data was to be collected in a database for the purpose of broader analysis. Third, the UNC work was done at RENCI (Renaissance Computing Institute), which has extensive expertise and facilities for scientific data visualization, so we acted in an ongoing consulting and support role in that area.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Code Verification Capabilities and Assessments in Support of ASC V&V Level 2 Milestone #6035

This document provides a summary of the code verification activities supporting the FY17 Level 2 V&V milestone entitled “Deliver a Capability for V&V Assessments of Code Implementations of Physics Models and Numerical Algorithms in Support of Future Predictive Capability Framework Pegposts.” The physics validation activities supporting this milestone are documented separately. The objectives of this portion of the milestone are: 1) Develop software tools to support code verification analysis; 2) Document standard definitions of code verification test problems; and 3) Perform code verification assessments (focusing on error behavior of algorithms). This report and a set of additional standalone documents serve as the compilation of results demonstrating accomplishment of these objectives.

97 MATHEMATICS AND COMPUTING↗