Search NASA⌕ Search

SEARCH · Search NASA

Results for “SCALE Code System”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

About the Impact of the Unresolved Resonance Region in Monte Carlo Simulations of Sodium Fast Reactors

In the last few years, and within the framework of different European projects, KENO-VI code from SCALE system has been employed to perform detailed continuous-energy Mon-te Carlo transport calculations for advanced fast reactors. The core characterization of both the sodium-cooled ASTRID and the lead-cooled ALFRED reactors was performed during the FP7 cross-cutting ESNII+ project; more recently, core calculations for the sodium-cooled Superphénix reactor and the improved European Sodium Fast Reactor design were performed within the HORIZON2020 ESFR-SMART project. In all cases, the effective multiplication factor predicted by KENO-VI was systematically higher (around 400-500 pcm) than the values computed by MCNP and Serpent Monte Carlo codes, using the same nu-clear data library.In order to provide insight into the origin of the observed discrepancies, a simplified 2D MOX-fueled SFR pin-cell benchmark has been launched. The multiplication factor, as well as 1-group and VITAMINJ 175-group cross-sections computed by KENO-VI, Serpent and MCNP codes employing ENDF/B-VII.1 data library, have been compared.Significant differences between KENO-VI and the other codes have been found in the un-resolved resonance regions of 239Pu and 241Pu capture and production cross sections, while negligible differences appeared outside those energy ranges. On the other hand, calculations without using probability tables have shown very good agreement. Quantita-tive comparison is presented and analyzed, along with a discussion of the impact of the probability-table treatment in the three codes for MOX-fueled systems with typical SFR spectrum.

Jimenez carrascosa, A↗

Extending the Nuclide Inventory Validation Basis for High-Burnup Fuel with New Radiochemical Assay Data

Efforts are underway at Oak Ridge National Laboratory to improve the nuclide inventory validation basis for spent nuclear fuel at high burnups. Recently conducted radiochemical assay experiments provided new measurement data for nine samples of fuel irradiated in a pressurized water reactor, with estimated sample burnups in the 30 to 70 GWd/t range. This type of destructive assay data is essential for validating computational methods, tools, and nuclear data applied in nuclear safety analyses and for improving our understanding of the bias and uncertainty in code predictions. The measurement data include key actinides and fission products that span a gamut of needs and interests for nuclear science and engineering applications in criticality safety, reactor physics, nuclide inventory, decay heat, and radiation shielding. The SCALE 6.3 code system with ENDF/B-VII.1 cross-section libraries was used to simulate the irradiation histories of the measured fuel samples. The calculated nuclide concentrations are compared to corresponding measurement data. The significance of the comparisons is discussed, emphasizing how the addition of the new measurement data fills gaps in the validation basis at high burnups and contributes to the decrease in bias and uncertainty for predicted nuclide concentrations. The discussion addresses the effect of the sample burnup used in the simulation—which is based on reactor operator records or on calibration to measured data for burnup indicator fission products—on the validation results.

Nuclide inventory↗

Kernel Density Estimators for SAMPLER Analyses [Slides]

This presentation discusses Sampler, specifically in terms of uncertainty quantification and states that for many simulation tasks nominal input parameters have complicated influence on the uncertainty in the final result. Sampler provides a convenient means to perform stochastic uncertainty quantification, brute force sampling parameter space and tallying changes in outputs, and it can be coupled with any SCALE modules. Perturbed responses can be collected and analyzed with standard statistical measures; however, two problems have been identified. The first problem is that for Monte Carlo simulations the individual response uncertainties are ignored, and the second problem is that features of distributions are easily hidden in histograms. They found that the solution is to add Kernel Density Estimator (KDE) analysis to the Sampler output. In conclusion, Sampler is a unique tool for stochastic uncertainty propagation, and it can be coupled with any module in SCALE. Current functionality is limited to histograms, and it may not adequately account for sample uncertainties and/or not adequately display underlying distribution features. Work is in progress to implement KDE in the Sampler module and no additional flagging is needed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear Data Libraries Sensitivity Studies for ORSA Using SCALE

Subcritical assemblies offer valuable training capabilities in nuclear criticality safety (NCS) for individuals handling fissile material. At Oak Ridge National Laboratory(ORNL), the Oak Ridge Subcritical Assembly (ORSA), a new experimental facility, is being established to provide hands-on training for the Nuclear Criticality Safety Program (NCSP).It is essential to accurately determine the neutron multiplication factor (keff) to ensure that ORSA remains subcritical and safe during operations. This study investigated the sensitivity of keff to variations across Evaluated Nuclear Data File (ENDF/B) libraries, consisting of ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1. The analysis was conducted using the CSAS6 sequence in the SCALE-6.3 code system. Individual isotopes in the ORSA model were replaced one at a time with ENDF/B-VII.1 as the base library and changing to ENDF/B-VIII.0or ENDF/B-VIII.1. The results demonstrated that the changes in keffof the nuclides associated with the ORSA model were mostly within the uncertainty of the base model (~24 pcm), except for primary nuclides like Uranium-235and H-poly with few other nuclides. The relative delta keff values of Uranium-235 and H-poly, expressed in pcm, were +280 and -212 in ENDF/B-VIII.0 and +332 and -275 in ENDF/B-VIII.1, respectively, which were notable changes in reactivity. These results demonstrated that ORSA was largely insensitive to variations across these nuclear data libraries.

Hong, Evan [North Carolina State University]↗

ADETS User Manual

ADETS (Automated DEpletion Transport System) calculates coupled neutronic/isotopic results for nuclear systems and produces a large number of criticality and burnup results based on various material feed/removal specifications, power(s), and time intervals. ADETS is a fully automated tool that links the LANL MCNP Monte Carlo transport code with the SCALE system, specifically ORIGEN-ARP (radioactive decay and burnup code). Additionally, ADETS can compute the dose at various distances of the reactor using a point gamma source, whose characteristics are computed by ORIGEN-ARP. In addition, ADETS is fully coupled with the Uncertanty Quantification, Probabilistic Risk Assessment, Parameter Optimization and Data Analysis code RAVEN (developed at the Idaho National Laboratory as well). RAVEN is released with ADETS but MCNP and SCALE (ORIGEN-ARP) are not included.

42 ENGINEERING↗

Expanded Validation of Uranium Systems with the KENO Monte Carlo Codes and SCALE 6.2.4 [Slides]

There is likely a trend in bias for uranium-fueled systems as a function of spectrum, but probably not as a function of enrichment. Some enrichment trends appear, but are potentially the result of outlier experimental results; further investigation is needed where possible. There are good examples in this data of markedly different biases from different systems. Future work includes review and vetting for inclusion in the VALID library and analysis with TSUNAMI.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Expanded Validation of Uranium Systems with the KENO Monte Carlo Codes and SCALE 6.2.4

The development of accident-tolerant fuel and advanced reactors fueled with high-assay low-enriched uranium has created an increased focus on systems with uranium enrichments between 5 and 20 wt % 235 U. SCALE criticality safety validation reports generated at Oak Ridge National Laboratory contain very few benchmark experiments in this range. The work presented herein is an initial investigation of expanded validation of the entire range of enrichments for uranium-fueled systems, including consideration of over 300 newly developed KENO models. The vast majority of the new models are between 4.4 wt % and 37.5 wt % enrichment and include a variety of thermal, intermediate, and mixed neutron spectra. This work finds no strong evidence of an enrichment-dependent bias, but there are indications of a spectrum-dependent bias. There are, however, only a limited number of independent experiments included at this time to distinguish intermediate enrichment effects from intermediate spectrum effects. Further work in this area is ongoing for intermediate enrichment systems. In the long run, the addition of these models to the Verified, Archived Library and Inputs and Data is planned, as funding and staff availability allow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Productive Programming of Distributed Systems with the SHAD C++ Library

High-performance computing (HPC) is often perceived as a matter of making large-scale systems (e.g., clusters) run as fast as possible, regardless the required programming effort. However, the idea of "bringing HPC to the masses" has recently emerged. Inspired by this vision, we have designed SHAD, the Scalable High-performance Algorithms and Data-structures library. SHAD is open source software, written in C++, for C++ developers. Unlike other HPC libraries for distributed systems, which rely on SPMD models, SHAD adopts a shared-memory programming abstraction, to make C++ programmers feel at home. Underneath, SHAD manages tasking and data-movements, moving the computation where data resides and taking advantage of asynchrony to tolerate network latency. At the bottom of his stack, SHAD can interface with multiple runtime systems: this not only improves developer’s productivity, by hiding the complexity of such software and of the underlying hardware, but also greatly enhance code portability. Thanks to its abstraction layers, SHAD can indeed target different systems, ranging from laptops to HPC clusters, without any need for modifying the user-level code. We have prototyped and open-sourced the implementation of (a subset of) the C++ standard library (STL) targeting multi-node HPC clusters. Our work allows plain STL-based C++ code to scale on HPC systems, with no need for rewriting the code to exploit the complex hardware. SHAD is available under Apache v2 License at https://github.com/pnnl/SHAD. In this paper we overview the design of the SHAD library, depicting its main components: runtime systems abstractions for tasking; parallel and distributed data-structures; STL-compliant interfaces and algorithms.

Castellana, Vito G.↗

SCALE 6.2.4 Validation for Light Water Reactor Decay Heat Analysis

Energy release from the decay of radionuclides in nuclear fuel after its discharge from reactor is a critical parameter for design, safety, and licensing analyses of used nuclear fuel storage, transportation, and repository systems. Well-validated computational tools and nuclear data are essential for decay heat prediction. This paper summarizes the validation of the SCALE nuclear analysis code system version 6.2.4, used with ENDF/B-VII.1 libraries, for decay heat analysis of light water reactor used fuel. The experimental data used for validation include full-assembly decay heat measurements that cover assembly burnups of 5 to 51 GWd/tonne U, cooling times after discharge in the 2- to 27-year range, and initial fuel enrichments up to 4 wt% 235U. The comparison between calculated (C) and experimental (E) decay heat showed very good agreement, with an average C/E over all considered measurements of 1.006 (σ = 0.016) for pressurized water reactor and 0.984 (σ = 0.077) for boiling water reactor assembly measurements. The effect of using assembly-average versus axially varying modeling data on the calculated decay heat, important to thermal analyses for used fuel transportation and storage systems, is discussed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Asynchronous-many-task systems: Challenges and opportunities - Scaling an AMR astrophysics code on exascale machines using Kokkos and HPX

Dynamic and adaptive mesh refinement is pivotal in high-resolution, multi-physics, multi-model simulations, necessitating precise physics resolution in localized areas across expansive domains. Today’s supercomputers’ extreme heterogeneity presents a significant challenge for dynamically adaptive codes, highlighting the importance of achieving performance portability at scale. Our research focuses on astrophysical simulations, particularly stellar mergers, to elucidate early universe dynamics. Here, we present Octo-Tiger, leveraging Kokkos, HPX, and SIMD for portable performance at scale in complex, massively parallel adaptive multi-physics simulations. Octo-Tiger supports diverse processors, accelerators, and network backends. Experiments demonstrate exceptional scalability across several heterogeneous supercomputers including Perlmutter, Frontier, and Fugaku, encompassing major GPU architectures and x86, ARM, and RISC-V CPUs. Parallel efficiency of 47.59% (110,080 cores and 6880 hybrid A100 GPUs) on a full-system run on Perlmutter (26% HPCG peak performance) and 51.37% (using 32,768 cores and 2048 MI250X) on Frontier are achieved.

97 MATHEMATICS AND COMPUTING↗

Dose Rate Analysis of the WCS Consolidated Interim Storage Facility

This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.

61 RADIATION PROTECTION AND DOSIMETRY↗

Corrections to official ENDF/B Evaluation Releases for SCALE Nuclear Data

The official ENDF/B nuclear data evaluation releases in the past two decades have incrementally incorporated more detailed information, new nuclide evaluations, and very often have improved the accuracy of radiation transport codes when modeling shielding, fission reactors, criticality benchmarks, and fusion systems. However, like any large collaborative data compilation, these releases have all included small errors. This report documents the small corrections to ENDF/B releases ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1 that have been applied during nuclear data processing to produce data libraries for the SCALE code system.

Brown, Jesse M. [Oak Ridge National Laboratory (OR↗

Scaling NASA Applications to 1024 CPUs on Origin 3K

The long and highly successful joint SGI-NASA research effort in ever larger SSI systems was to a large degree the result of the successful development of the MLP scalable parallel programming paradigm developed at ARC: 1) MLP scaling in real production codes justified ever larger systems at NAS; 2) MLP scaling on 256p Origin 2000 gave SGl impetus to productize 256p; 3) MLP scaling on 512 gave SGI courage to build 1024p O3K; and 4) History of MLP success resulted in IBM Star Cluster based MLP effort.

Taft, Jim↗

MLEC-Sim: A Simulator for Evaluating Multi-Level Erasure Coding

We present MLEC-Sim, a sophisticated simulator for Multi-Level Erasure Coding (MLEC), developed in approximately 13 KLOC. The simulator is engineered to analyze the impact of various system configurations and erasure coding policies on system durability and network overhead. It supports a comprehensive range of parameters including disk capacity, disk I/O bandwidth, failure rates, network bandwidth, and system scale, accommodating various erasure coding approaches such as Single-Level Erasure Coding (SLEC), Multi-Level Erasure Coding (MLEC), and Local Reconstruction Codes (LRC). MLEC-Sim provides support for multiple chunk placement policies, including clustered parity and declustered parity, and encompasses a variety of repair methods like Repair-ALL, Repair-FCO, Repair-HYB, and Repair-MIN. It is capable of simulating disk failures through a variety of means, including distribution-based or trace-based mechanisms, and can handle complex multi-level (de)clustered placements and repair processes. A key feature of MLEC-Sim is its adoption of the splitting simulation method for evaluating system durabilities at extremely high levels, which are challenging to assess with traditional simulation approaches. This feature allows for a detailed evaluation of system resilience under a range of conditions, aiding in the selection of appropriate erasure coding solutions for enhancing system durability. MLEC-Sim contributes to the field of data storage and reliability by providing a tool for the detailed evaluation of the durability and efficiency of erasure coding configurations, intended for use by researchers and practitioners in the design and optimization of storage systems.

Wang, Meng↗

Application of the Denovo Discrete Ordinates Radiation Transport Code to Large-Scale Fusion Neutronics

Fusion energy systems pose unique challenges to the modeling and simulation community. These challenges must be met to ensure the success of the ITER experimental fusion reactor. ITER’s complex systems require detailed modeling that goes beyond the scale of comparable simulations to date. In this work, the Denovo radiation transport code was used to calculate neutron fluence and kerma for the JET streaming benchmark. This work was performed on the Titan supercomputer at the Oak Ridge Leadership Computing Facility. Denovo is a novel three-dimensional discrete ordinates transport code designed to be highly scalable. Sensitivity studies have been completed to examine the impact of several deterministic parameters. Furthermore, results were compared against experiment as well as the MCNP and Shift Monte Carlo codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SCALE 6.2.4: Validation Overview

SCALE is a multi-application code system with tools for reactor physics, criticality safety, radiation shielding, and spent fuel characterization for a range of nuclear systems. The last of the SCALE 6.2 series, version 6.2.4, was released in June 2020. Validation studies in nuclear criticality safety, reactor physics, and shielding applications were performed to confirm results, assess accuracy, and identify biases from SCALE simulations. Nuclear criticality safety applications yield small code biases for a wide range of systems, generally overestimating cross section uncertainty. Reactor physics applications improvements in nuclide inventory and full-core eigenvalue predictions, with similar performance in decay heat predictions. Radiation shielding applications yield sufficient agreement in comparisons to measurements and alternate simulations for eight shielding experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Comparative study with two reduced test facilities for air-cooled RCCS scaling law

Korea Atomic Energy Research Institute (KAERI) and Argonne National Laboratory (ANL) conducted steady-state tests at constant heat flux conditions to investigate the scale effect on heat removal behavior of the air-cooled Reactor Cavity Cooling System (RCCS). Two differently scaled-down test facilities were used in this study, KAERI’s ¼ scale with 4.0-m cavity height, and ANL’s ½ scale with 6.8-m cavity height. Two scaling laws were proposed to simulate the radiation across the cavity and the buoyancy-driven duct flow in the riser, respectively. The test matrix for KAERI and ANL test facilities focused on the scaling effect of the heated riser length with constant heat flux at the PMR200 RCCS design. The test results showed that mixed convection in the riser duct is an important factor for accurately extrapolating the thermo-fluid behavior in the prototype from the test results in the scale-down facilities. The system analysis code, GAMMA+, with improved heat transfer models predicted fairly well the air-cooled RCCS test data from two facilities. GAMMA+ analysis results showed that the predicted radiation fractions on the heated plate in the scale-down test conditions were larger than those on the reactor vessel in the prototype. The scaling law for air-cooled RCCS was improved by considering the mixed convection in the riser duct and the same radiation fraction on the heated plate. The mixed convection effect was calculated by the height ratio. The same radiation fraction provided a conservative extrapolation from the scale-down test results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗