Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code assessment”

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,765 records · Page 98

Compare predictions of transient fission gas release by empirical and mechanistic models to experiments in high burnup UO 2 fuel

Understanding and predicting fuel performance at high burnup require improving our understanding of transient fission gas release. High-burnup operations enable new mechanisms of fission gas release, which affect fuel performance. The Nuclear Regulatory Commission has recently published its interpretation of existing fuel fragmentation, relocation, and dispersal data in a research information letter. There, transient fission gas release was identified as one of the main factors that contributes to fuel fragmentation, relocation, and dispersal, and therefore limits fuel extension to high burnup. However, transient fission gas release is a complex phenomenon that cannot be fully described by simple empirical descriptions. This report summarizes the development of a mechanistic model for high-burnup transient fission gas release in the fuel performance code BISON. This research was supported by the Nuclear Energy Advanced Modeling and Simulation program during fiscal year 2023 to improve our understanding of high-burnup transient fission gas release and ability to predict it as a function of operation history. To support the development of a mechanistic transient fission gas release model, the existing Simple Integrated Fission Gas Release and Swelling (Sifgrs) model in BISON has been completely refactored to make it more modular and extensible. This effort supports the model's application to high-burnup conditions, its extension to other fuel forms, and the continuous improvement of its current features. Once refactoring was completed, models for high-burnup structure formation, fission gas transfer from non-restructured fuel to high-burnup structure, high-burnup structure intragranular and intergranular fission gas behavior, high-burnup structure bubble evolution, fuel pulverization, and the resulting transient fission gas release were tested and implemented in the Simple Integrated Fission Gas Release and Swelling (Sifgrs) model or tightly coupled to it. The new mechanistic model was then compared to an empirical model developed in parallel by a Nuclear Energy University Program project using a Studsvik high-burnup loss-of-coolant-accident assessment case. Finally, the report details the preliminary BISON results for a benchmark activity organized by the Nuclear Energy Agency to evaluate fuel performance codes' predictive capabilities for burst fission gas release. This work represents an important step toward a mechanistic understanding of fission gas release in high-burnup conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Advanced Pellet-Cladding Interaction Modeling using the US DOE CASL Fuel Performance Code: Peregrine

The US DOE’s Consortium for Advanced Simulation of LWRs (CASL) program has undertaken an effort to enhance and develop modeling and simulation tools for a virtual reactor application, including high fidelity neutronics, fluid flow/thermal hydraulics, and fuel and material behavior. The fuel performance analysis efforts aim to provide 3-dimensional capabilities for single and multiple rods to assess safety margins and the impact of plant operation and fuel rod design on the fuel thermo-mechanical-chemical behavior, including Pellet-Cladding Interaction (PCI) failures and CRUD-Induced Localized Corrosion (CILC) failures in PWRs. [1-3] The CASL fuel performance code, Peregrine, is an engineering scale code that is built upon the MOOSE/ELK/FOX computational FEM framework, which is also common to the fuel modeling framework, BISON [4,5]. Peregrine uses both 2-D and 3-D geometric fuel rod representations and contains a materials properties and fuel behavior model library for the UO2 and Zircaloy system common to PWR fuel derived from both open literature sources and the FALCON code [6]. The primary purpose of Peregrine is to accurately calculate the thermal, mechanical, and chemical processes active throughout a single fuel rod during operation in a reactor, for both steady state and off-normal conditions.

CASL, Consortium for Advanced Simulation of LWRs, ↗

Assessment of different coarse graining strategies to simulate polydisperse gas-solids flow

Continuum methods require the additional development of solids stress closures for polydisperse powders based on complex kinetic theories that are non-trivial to develop, code, and numerically converge for the wide range of fluidization regimes from very dilute to dense/frictional flow limit. On the other hand, it is straightforward to model the flow of polydisperse granular materials by treating particles as discrete rigid bodies that are tracked following simple physical laws of motion. The coarsening of these discrete methods by lumping several particles in a parcel alleviates the significant computational cost associated with these discrete methods while introducing some inaccuracies in the numerical results. In this research, we explore two different coarse graining methods that can be applied to polydisperse powders, namely the same statistic weight method (SSW) and the same size parcel method (SSP), and assess their accuracy by comparison with the finest simulation results obtained with a discrete element method (DEM). For Geldart group B powders fluidized at a relative low superficial velocity, the numerical results indicate that the SSW is more accurate than the SSP method. For type A powders fluidized at relatively high velocity, these two methods predict similar results. Interestingly, up to four times increase in the speed of simulation of the SSP method was obtained because the original polydisperse powder is scaled to a mono-disperse system in terms of particle-particle collision. Lastly, these results suggest that the SSP method is more favorable for the simulation of fluidized beds due to its accuracy and efficiency while the SSW method may be used for granular flow and dense fluidized bed systems where capturing the size segregation of particles due to collision is important.

42 ENGINEERING↗

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗

Industrial Geospatial Analysis Tool for Energy Evaluation (IGATE-E)

IGATE-E is an energy analysis tool for industrial energy evaluation. The tool applies statistical modeling to multiple publicly available datasets and provides information at the geospatial resolution of zip code using bottom up approaches. Within each zip code, the current version of the tool estimates electrical energy consumption of manufacturing industries based on each type of industries using DOE s Industrial Assessment Center database (IAC-DB) and DOE s Energy Information Administration Manufacturing Energy Consumption Survey database (EIA-MECS DB), in addition to other commercially available databases such as the Manufacturing News database (MNI, Inc.). Ongoing and future work include adding modules for the predictions of fuel energy consumption streams, manufacturing process steps energy consumption, major energy intensive processes (EIPs) within each industry type among other metrics of interest. The tool provides validation against DOE s EIA-MECS state level energy estimations and permits several statistical examinations. IGATE-E is intended to be a decision support and planning tool to a wide spectrum of energy analysts, researchers, government organizations, private consultants, industry partners, and alike.

Energy Analysis Tool↗

SCALE Continuous-Energy Eigenvalue Sensitivity Coefficient Calculations

Sensitivity coefficients describe the fractional change in a system response that is induced by changes to system parameters and nuclear data. The Tools for Sensitivity and UNcertainty Analysis Methodology Implementation (TSUNAMI) code within the SCALE code system makes use of eigenvalue sensitivity coefficients for an extensive number of criticality safety applications, including quantifying the data-induced uncertainty in the eigenvalue of critical systems, assessing the neutronic similarity between different critical systems, and guiding nuclear data adjustment studies. The need to model geometrically complex systems with improved fidelity and the desire to extend TSUNAMI analysis to advanced applications has motivated the development of a methodology for calculating sensitivity coefficients in continuous-energy (CE) Monte Carlo applications. The Contributon-Linked eigenvalue sensitivity/Uncertainty estimation via Tracklength importance CHaracterization (CLUTCH) and Iterated Fission Probability (IFP) eigenvalue sensitivity methods were recently implemented in the CE-KENO framework of the SCALE code system to enable TSUNAMI-3D to perform eigenvalue sensitivity calculations using continuous-energy Monte Carlo methods. This work provides a detailed description of the theory behind the CLUTCH method and describes in detail its implementation. This work explores the improvements in eigenvalue sensitivity coefficient accuracy that can be gained through the use of continuous-energy sensitivity methods and also compares several sensitivity methods in terms of computational efficiency and memory requirements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Eucalyptus – An Analysis Suite for Fault Trees with Uncertainty Quantification

Eucalyptus is a novel code developed at Lawrence Livermore National Laboratory to incorporate uncertainty quantification into Fault Tree Analysis (FTA). This tool addresses the challenge of imperfect knowledge in “grey-box” systems by allowing analysts to incorporate and propagate uncertainty from component-level assessments to system-level effects. Eucalyptus facilitates a consistent evaluation of the impact of subject matter expert judgment and knowledge gaps on overall system response by Monte Carlo generation of possible system fault trees, sampling probabilities of the existence of subsystems and components. Here, the code supports the specification of fault trees through text and allows export to various formats, including auto-generated images, easing analysis and reducing errors. It has undergone extensive verification testing, demonstrating its reliability and readiness for deployment, and leverages on-node parallelism for rapid analysis. Example analyses are shown that include the identification of system failure paths and quantification of the value of further information about system components.

Fault Tree Analysis↗

Assessing plasma face component thermal response to rotating 3D magnetic fields for SPARC tokamak

Thermal response simulations of plasma-facing components (PFCs) in the SPARC tokamak, performed with the HEAT code, show that three-dimensional (3D) heat loads resulting from stationary n=1 perturbations require highly radiative scenarios, with up to 95% of the power crossing the separatrix (P SOL ) being radiated, to maintain PFC temperatures within acceptable operational limits, whereas the application of slowly rotating 3D fields substantially reduces the thermal loads. The HEAT module, developed to predict heat loads from non-axisymmetric plasmas, is extended to model time-dependent heat flux patterns generated by rotating 3D fields, and a comprehensive thermal analysis is performed on PFCs subjected to both the maximum and minimum power loads, as well as to rotating heat flux distributions, to evaluate the temperature evolution for varying perturbation amplitudes and rotation frequencies. The extension of this analysis to 3D fields with toroidal mode number n=2 shows that this configuration leads to weaker localized heat flux peaks relative to the n=1 case, enabling safe operation with less than 80% of the power radiated when static 3D fields of low amplitude are applied, while using slowly rotating fields at higher amplitudes. These results indicate that n=2 perturbations are generally less detrimental to divertor power exhaust, emphasizing the strong dependence of divertor power exhaust on the characteristics of the applied 3D fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Justification for Use of the HELP Model to Estimate Infiltration Rates for the E-Area Low-Level Waste Facility Performance Assessment

Whiteside et al. conducted an in-depth literature review and performed basic infiltration model simulations to evaluate the Hydrologic Evaluation of Landfill Performance (HELP) model against five alternative computational codes utilizing the Richards equation for variably-saturated flow. Of the five alternatives considered (LEACHM, UNSAT-H, SVFlux, HYDRUS-2D3D, and VADOSE/W), only HYDRUS-2D3D and VADOSE/W were selected for further evaluation, which consisted of a side-by-side comparison versus HELP of calculated infiltration rates through a one-layer soil column with two percent slope.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Code Standard Evaluation for ASME AG-1-2019, Section FL

ASME AG-1 Section FL Deep Bed Sand Filters includes a requirement to use Special Test Sections to provide a representative assessment of the filter efficiency. The SRPPF project design basis is to test the entire filter bed (versus a representative sample) following the same basis methodology to test the existing filters at the Savannah River Site. Using a combination of: 1) the filter testing experience at SRS since the early 1970’s, 2) the methodology found in ASME AG-1 for testing high capacity HEPA filter housings and 3) requirements for stack or duct sampling systems published by the Environmental Protection Agency (EPA), the Special Test Sections methodology for Deep Bed Sand (DBS) filters in ASME AG-1 Section FL has been replaced with a methodology to test the entire DBS filter system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Development of remote sensing technology in New Zealand, part 1. Mapping land use and environmental studies in New Zealand, part 2. Indigenous forest assessment, part 3. Seismotectonic, structural, volcanologic and geomorphic study of New Zealand, part 4

The author has identified the following significant results. As part of the tape reformatting process, a simple coded picture output program was developed. This represents Pixel's radiance level by one of a 47 character set on a nonoverprinting line printer. It not only has aided in locating areas for the reformatting process, but has also formed the foundation for a supervised clustering package. This in turn has led to a simplistic but effective thematic mapping package.

Probine, M. C.↗

Assessment of the operating characteristics of the SSME LOX turbopump pump-end bearing

A bearing/shaft model of the SSME LOX turbopump was developed using the SHABERTH bearing/shaft math modeling computer code. A previously developed bearing/shaft thermal model of the SSME LOX turbopump turbine and bearing was used in conjunction with SHABERTH to evaluate the thermomechanical operating characteristics of the LOX turbopump end bearings. Results show that for the two unmounted diametrical clearances evaluated (4.0 mils and 6.3 mils), the inboard pump end bearing supports about 81% of the isolator load for the small clearance and 77% of the isolator load for the larger clearance. Bearing clearance changes due to thermal effects were 40% for the 4.0 mil diametrical clearance case and 19% for the 6.3 mil clearance case evaluated. The thermal analysis included evaluation of bearing temperatures for a subcooled case and a saturated case. Results indicate that no drastic temperature change occurred between the two cases. Since the rolling element and race surfaces of the subcooled case were at temperatures sufficiently high enough to be vapor blanketed, exceeding saturation temperature at the bearing inlet did not increase surface temperatures greatly.

New, L. S.↗

Fast estimation of dense gas dispersion from multiple continuous CO 2 surface leakage sources for risk assessment

Surface leakage of CO 2 , and associated potential impacts on health, safety, and the environment (HSE) are considered hazards of geologic carbon sequestration (GCS). There are two challenges associated with impact assessment of CO 2 surface dispersion. First, the fact that CO 2 is a dense gas makes its dispersion in air a complex process. Rigorous numerical solutions for modeling concentration distributions are relatively time-consuming. Second, impact assessment requires consideration of uncertainty, e.g., quantification of how much uncertainty is propagated through input parameters to model outputs by carrying out large numbers of model runs. In order to assess the potential consequences of surface leakage of CO 2 , it is useful to have a model that executes very quickly for repeated model calculations (e.g., in Monte Carlo mode) of the atmospheric dispersion of CO 2 (concentrations as a function of space and time). In addition, the model should be able to handle multiple surface leakage sources. In this work, we have extended the nomograph approach of Britter and McQuaid (1988) for estimating dense gas plume length from single leakage source to multiple leakage sources. The method is very fast and therefore amenable to general system-level GCS risk assessment including uncertainty quantification within the framework of the National Risk Assessment Partnership (NRAP) Integrated Assessment Model (IAM). The method is conservative in that it assumes the wind could be from any direction, and it handles multiple sources by a simple superposition approach. The method produces results in reasonable agreement with a sophisticated computational fluid dynamics (CFD) code, but runs in a small fraction of the time.

58 GEOSCIENCES↗

Assessment of crack growth in a space shuttle main engine first-stage, high-pressure fuel turbopump blade

A two-dimensional finite element fracture mechanics analysis of a space shuttle main engine (SSME) turbine blade firtree was performed using the MARC finite element code. The analysis was conducted under combined effects of thermal and mechanical loads at steady-state conditions. Data from a typical engine stand cycle of the SSME were used to run a heat transfer analysis and, subsequently, a thermal structural fracture mechanics analysis. Temperature and stress contours for the firtree under these operating conditions were generated. High stresses were found at the firtree lobes where crack initiation was triggered. A life assessment of the firtree was done by assuming an initial and a final crack size.

Abdul-Aziz, Ali↗

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Oak Ridge OpenSHMEM Benchmark Suite

The assessment of application performance is a fundamental task in high-performance computing (HPC). The OpenSHMEM Benchmark (OSB) suite is a collection of micro-benchmarks and mini-applications/compute kernels that have been ported to use OpenSHMEM. Some, like the NPB OpenSHMEM benchmarks, have been published before while most others have been used for evaluations but never formally introduced or discussed. This suite puts them together and is useful for assessing the performance of different use cases of OpenSHMEM. This offers system implementers a useful means of measuring performance and assessing the effects of new features as well as implementation strategies. The suite is also useful for application developers to assess the performance of the growing number of OpenSHMEM implementations that are emerging. In this paper, we describe the current set of codes available within the OSB suite, how they are intended to be used, and, where possible, a snapshot of their behavior on one of the OpenSHMEM implementations available to us. We also include detailed descriptions of every benchmark and kernel, focusing on how OpenSHMEM was used. This includes details on the enhancements we made to the benchmarks to support multithreaded variants. We encourage the OpenSHMEM community to use, review, and provide feedback on the benchmarks.

Naughton III, Thomas↗

Analysis of OpenMP 4.5 Offloading in Implementations: Correctness and Overhead

The OpenMP language features have been evolving to meet the rapid development in hardware platforms. This journal focuses on evaluating implementations of OpenMP 4.5 target offload features in compilers such as Clang, XL and GCC that are an integral part of the software harness on supercomputers and clusters. We use Summit (Top supercomputer in the world as of November 2018) as one of our experimental setup. Such an effort is particularly critical on such supercomputers as that is being widely used by application developers to run their scientific codes at scale. Our tests not only evaluate the OpenMP implementations but also expose ambiguities within the OpenMP 4.5 specification. We also assess the overhead of the different OpenMP runtimes in relationship to the different directives and clauses. This helps in assessing the interaction of different OpenMP directives independent of other application artifacts. We are aware that the implementations are constantly evolving and Summit is advertised as having only partial OpenMP 4.x support. This is a synergistic effort to help identify and fix bugs in features’ implementations that are required by applications and prevent deployment delays later on. Going forward, we also plan to interact with standard benchmarking organizations like SPEC/HPG to donate our tests and mini-apps/kernels for potential inclusion in the next release versions of SPEC benchmark suite.

97 MATHEMATICS AND COMPUTING↗