Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multigroup”

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 163 records · Page 9

Parallel computation safety analysis irradiation targets fission product molybdenum in neutronic aspect using the successive over-relaxation algorithm

One of the research activities in support of commercial radioisotope production program is a safety research on target FPM (Fission Product Molybdenum) irradiation. FPM targets form a tube made of stainless steel which contains nuclear-grade high-enrichment uranium. The FPM irradiation tube is intended to obtain fission products. Fission materials such as Mo{sup 99} used widely the form of kits in the medical world. The neutronics problem is solved using first-order perturbation theory derived from the diffusion equation for four groups. In contrast, Mo isotopes have longer half-lives, about 3 days (66 hours), so the delivery of radioisotopes to consumer centers and storage is possible though still limited. The production of this isotope potentially gives significant economic value. The criticality and flux in multigroup diffusion model was calculated for various irradiation positions and uranium contents. This model involves complex computation, with large and sparse matrix system. Several parallel algorithms have been developed for the sparse and large matrix solution. In this paper, a successive over-relaxation (SOR) algorithm was implemented for the calculation of reactivity coefficients which can be done in parallel. Previous works performed reactivity calculations serially with Gauss-Seidel iteratives. The parallel method can be used to solve multigroup diffusion equation system and calculate the criticality and reactivity coefficients. In this research a computer code was developed to exploit parallel processing to perform reactivity calculations which were to be used in safety analysis. The parallel processing in the multicore computer system allows the calculation to be performed more quickly. This code was applied for the safety limits calculation of irradiated FPM targets containing highly enriched uranium. The results of calculations neutron show that for uranium contents of 1.7676 g and 6.1866 g (× 10{sup 6} cm{sup −1}) in a tube, their delta reactivities are the still within safety limits; however, for 7.9542 g and 8.838 g (× 10{sup 6} cm{sup −1}) the limits were exceeded.

07 ISOTOPE AND RADIATION SOURCES↗

SU-G-TeP1-15: Toward a Novel GPU Accelerated Deterministic Solution to the Linear Boltzmann Transport Equation

Purpose: To develop a Graphic Processor Unit (GPU) accelerated deterministic solution to the Linear Boltzmann Transport Equation (LBTE) for accurate dose calculations in radiotherapy (RT). A deterministic solution yields the potential for major speed improvements due to the sparse matrix-vector and vector-vector multiplications and would thus be of benefit to RT. Methods: In order to leverage the massively parallel architecture of GPUs, the first order LBTE was reformulated as a second order self-adjoint equation using the Least Squares Finite Element Method (LSFEM). This produces a symmetric positive-definite matrix which is efficiently solved using a parallelized conjugate gradient (CG) solver. The LSFEM formalism is applied in space, discrete ordinates is applied in angle, and the Multigroup method is applied in energy. The final linear system of equations produced is tightly coupled in space and angle. Our code written in CUDA-C was benchmarked on an Nvidia GeForce TITAN-X GPU against an Intel i7-6700K CPU. A spatial mesh of 30,950 tetrahedral elements was used with an S4 angular approximation. Results: To avoid repeating a full computationally intensive finite element matrix assembly at each Multigroup energy, a novel mapping algorithm was developed which minimized the operations required at each energy. Additionally, a parallelized memory mapping for the kronecker product between the sparse spatial and angular matrices, including Dirichlet boundary conditions, was created. Atomicity is preserved by graph-coloring overlapping nodes into separate kernel launches. The one-time mapping calculations for matrix assembly, kronecker product, and boundary condition application took 452±1ms on GPU. Matrix assembly for 16 energy groups took 556±3s on CPU, and 358±2ms on GPU using the mappings developed. The CG solver took 93±1s on CPU, and 468±2ms on GPU. Conclusion: Three computationally intensive subroutines in deterministically solving the LBTE have been formulated on GPU, resulting in two orders of magnitude speedup. Funding support from Natural Sciences and Engineering Research Council and Alberta Innovates Health Solutions. Dr. Fallone is a co-founder and CEO of MagnetTx Oncology Solutions (under discussions to license Alberta bi-planar linac MR for commercialization).

60 APPLIED LIFE SCIENCES↗

SCALE 6: Comprehensive Nuclear Safety Analysis Code System

Version 6 of the Standardized Computer Analyses for Licensing Evaluation (SCALE) computer software system developed at Oak Ridge National Laboratory, released in February 2009, contains significant new capabilities and data for nuclear safety analysis and marks an important update for this software package, which is used worldwide. This paper highlights the capabilities of the SCALE system, including continuous-energy flux calculations for processing multigroup problem-dependent cross sections, ENDF/B-VII continuous-energy and multigroup nuclear cross-section data, continuous-energy Monte Carlo criticality safety calculations, Monte Carlo radiation shielding analyses with automated three-dimensional variance reduction techniques, one- and three-dimensional sensitivity and uncertainty analyses for criticality safety evaluations, two- and three-dimensional lattice physics depletion analyses, fast and accurate source terms and decay heat calculations, automated burnup credit analyses with loading curve search, and integrated three-dimensional criticality accident alarm system analyses using coupled Monte Carlo criticality and shielding calculations.

99 GENERAL AND MISCELLANEOUS↗

Demonstration of MOOSE-based Griffin reactor physics, code for heterogeneous lead-cooled fast reactor analysis

The MOOSE-based reactor physics code Griffin was assessed on a heterogeneous pin-resolved model of a prototype lead-cooled fast reactor assembly. This model was developed in preparation for future use in MOOSE-based multiphysics calculations for computing hot channel factors. Heterogeneous multigroup cross sections were prepared using the fast reactor multi-group cross section processing code MC{sup 2}-3 using a two-step method. Griffin simulations were performed using the DFEM-SN solver on 576 cores on Argonne's LCRC cluster. Diffusion-based acceleration methods were applied (NDA and CMFD). Reference solutions were generated with continuous energy MCNP and the hybrid MOC/finite element solver PROTEUS-MOC for code-to-code comparison. Space-angle convergence studies were conducted to observe convergence in k-eigenvalue and axial pin power distributions. The fully resolved Griffin calculation was within 68 pcm of the MCNP eigenvalue and exhibited max 1.2% relative error in the axial pin power distribution. Griffin produced nearly identical results to PROTEUS-MOC when using the same 9-group multigroup cross-section set. Griffin demonstrated favorable scaling in wall-clock time and memory usage when using diffusion-based acceleration methods. Griffin is capable of simulating the pin-resolved heterogeneous LFR assembly with good accuracy and performance, and is suitable for future use in coupled high-fidelity hot channel factor simulations. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

COMBINE7.1 - A Portable ENDF/B-VII.0 Based Neutron Spectrum and Cross-Section Generation Program

COMBINE7.1 is a FORTRAN 90 computer code that generates multigroup neutron constants for use in the deterministic diffusion and transport theory neutronics analysis. The cross-section database used by COMBINE7.1 is derived from the Evaluated Nuclear Data Files (ENDF/B-VII.0). The neutron energy range covered is from 20 MeV to 1.0E-5 eV. The Los Alamos National Laboratory NJOY code is used as the processing code to generate a 167 fine-group cross-section library in MATXS format for Bondarenko self-shielding treatment. Resolved resonance parameters are extracted from ENDF/B-VII.0 File 2 for a separate library to be used in an alternate Nordheim self-shielding treatment in the resolved resonance energy range. The equations solved for energy dependent neutron spectrum in the 167 fine-group structure are the B3 or B1 zero-dimensional approximations to the transport equation. The fine group cross sections needed for the spectrum calculation are first prepared by Bondarenko self-shielding interpolation in terms of background cross section and temperature. The geometric lump effect, when present, is accounted for by augmenting the background cross section. Nordheim self-shielded fine group cross sections for a material having resolved resonance parameters overwrite correspondingly the existing self-shielded fine group cross sections when this option is used. COMBINE7.1 coalesces fine group cross sections into broad group macroscopic and microscopic constants. The coalescing is performed by utilizing fine-group fluxes and/or currents obtained by spectrum calculation as the weighting functions. The multigroup constants may be output in any of several standard formats including INL format, ANISN 14** free format, CCCC ISOTXS format, and AMPX working library format. ANISN-PC, a one-dimensional (1-D) discrete-ordinate transport code, is incorporated into COMBINE7.1. As an option, the 167 fine-group constants generated by zero-dimensional COMBINE portion in the program can be used to calculate regionwise spectra in the 1-D ANISN portion, all internally to reflect the 1-D transport correction. The regionwise spectra are then used to generate mutigroup regionwise neutron constants. The 1-D neutron transport can be performed up to three stages, e.g., from a TRISO fuel to PEBBLE to 1-D full core wedge. In addition, COMBINE7.1 has now the capability of adjoint flux calculation through the 1-D ANISN transport. Photon transport capability is also added. For this, a photon production and photo-atomic cross section library, MATNG.LIB, was generated in MATXS format through NJOY code. The photon production cross section matrix is of 167 neutron - 18 photon groups. Photo-atomic cross sections, including heating, are in 18 energy groups.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Use of Covariances in a Consistent Data Assimilation for Improvement of Basic Nuclear Parameters in Nuclear Reactor Applications: From Meters to Femtometers

Recognition of the key role physics plays in innovative reactor design was already made clear by E. Wigner in his warning: “It has been a bit forgotten that in all really creative thinking in reactor design, a working knowledge of nuclear reaction theory is required.” More recently, extensive sensitivity and uncertainty studies and the availability of new covariance data have allowed preliminary quantification of the impact of current nuclear data uncertainties on the design parameters of the Generation-IV systems (both fast neutron and thermal neutron systems). Similar quantifications can be done for the parameters that characterize the associated innovative fuel cycle to insure sustainability, waste minimization and drastic reduction of the proliferation risks. In parallel, there has been a growing and very significant trend to develop a new generation of reactor simulation tools that rely more and more on first principles. It is safe to say that the next generation of innovative nuclear systems will be assessed with a completely new set of tools, more science based, in sharp contrast with what has been done in the past. However, the studies mentioned above point out that the present uncertainties in nuclear data should be significantly reduced in order to fully benefit from advances in modeling and simulation. Only a parallel effort in advanced simulation and nuclear data improvement will be able to provide designers with more general and well validated calculational tools to meet new design target accuracies. One further consideration related to the development of new advanced simulation tools is that there should be a more explicit link with the more fundamental physics parameters underlying the models used to describe cross sections that would avoid the use of processed, application-oriented, data (like multigroup cross sections). Finally, current methodologies for reducing uncertainties coming from nuclear data, which rely on the use of integral experiment information to perform statistical multigroup (i.e., energy spectrum weighted) cross section data assimilation should evolve to embrace new frontiers in advanced simulation. One of the drawbacks of the classical methodologies is related to the energy group structure and the type of neutron energy spectrum that are adopted in the assimilation (adjustment). In fact, after such an adjustment is performed, neutronic designers are then tied to this energy group structure and neutron energy spectrum when carrying out further calculations. In reality, this can be quite a limiting factor in view of the complex spectral issues that are involved in most reactor physics. This work combines novel, but proven, methodologies for overcoming these limitations. In fact, this is the first attempt to build up a link between the wealth of precise integral experiments and a basic theory of nuclear reactions. Essential ingredients of such a procedure, denominated here as assimilation, are covariances for model parameters and sensitivity matrices. The latter provide direct link between reaction theory and integral experiments. The result is a consistent data assimilation performed directly on the basic nuclear physics parameters that are being used in a variety of nuclear reaction mechanisms. The resulting improvement in their performance will consequently reduce related uncertainties when employed in reactor calculations. By using integral reactor physics experiments (meter scale), information is propagated back to the nuclear physics level (femtometers) covering a range of more than 13 orders of magnitude. The assimilation procedure should result in more accurate and more reliable evaluated data files of universal validity rather than tailored to a particular application. In fact, after data assimilation is carried out, the basic nuclear data file can be processed by a dedicated code into any energy group structure that the reactor physicist deems to be useful. On the other hand, integral experiments used in the assimilation should provide additional, possibly quite strict, constraints on the parameters entering nuclear reaction modeling, as well as the reaction models themselves. This report describes three years of combined research by Brookhaben National Lab (BNL) and Idaho National Lab (INL) on establishing viable assimilation methodology. The emphasis of this paper is on the empire code calculations to prepare the priors, sensitivity matrices and covariances for the model parameters. For completeness, we also include a short summary of the assimilation results produced by INL and reported earlier in annual reports [6], [7], [8].

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Verification of the CENTRM Module for Adaptation of the SCALE Code to NGNP Prismatic and PBR Core Designs

The generation of multigroup cross sections lies at the heart of the very high temperature reactor (VHTR) core design, whether the prismatic (block) or pebble-bed type. The design process, generally performed in three steps, is quite involved and its execution is crucial to proper reactor physics analyses. The primary purpose of this project is to develop the CENTRM cross-section processing module of the SCALE code package for application to prismatic or pebble-bed core designs. The team will include a detailed outline of the entire processing procedure for application of CENTRM in a final report complete with demonstration. In addition, they will conduct a thorough verification of the CENTRM code, which has yet to be performed. The tasks for this project are to: Thoroughly test the panel algorithm for neutron slowing down; Develop the panel algorithm for multi-materials; Establish a multigroup convergence 1D transport acceleration algorithm in the panel formalism; Verify CENTRM in 1D plane geometry; Create and test the corresponding transport/panel algorithm in spherical and cylindrical geometries; and, Apply the verified CENTRM code to current VHTR core design configurations for an infinite lattice, including assessing effectiveness of Dancoff corrections to simulate TRISO particle heterogeneity.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

(U) Feynman-Y calculations using PARTISN

A prescription for computing the Feynman Y as a function of coincidence gate width using a deterministic multigroup neutron transport code has been published and the results compared favorably with measurements of the BeRP ball. In this paper, we report on our project to implement the method and reproduce the results. There are several clarifications and corrections of the published prescription. We show results using two multigroup cross section libraries compared with measurements and with Monte Carlo results. Deterministic simulations of the mean count rates compare very favorably with previously published Monte Carlo results, and deterministic simulations of the Feynman Y asymptote compare somewhat favorably. In Feynman beta plots, the deterministic simulations reached the asymptotic value much sooner than did a fit to the measured data.

42 ENGINEERING↗

Generation of the Cross Section Library for PROTEUS

The cross sections for a high-fidelity neutron transport code PROTEUS can be generated in multiple ways. Multigroup cross sections in the ISOTXS format can be generated using the external codes such as MC 2 -3 and Monte Carlo codes such as Serpent and OpenMC. For the Monte Carlo codes, their multigroup cross section outputs should be converted to the ISOTXS format using a post processor named GenISOTXS developed by ANL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) LNK3DNT Capability in SENSMG

The SENSMG multigroup neutron sensitivity code has been extended to read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Not all of the capability of SENSMG is available for redoin/lnk3dnt inputs. The only responses allowed are k eff or α eigenvalues. Fixed-source problems are not allowed, and there is no mechanism to input reaction rates or reaction-rate ratios. The next section of this report describes the user interface. Section III describes the verification of the capability. Section IV is a summary and conclusions. The input files used for verification are listed in the appendix.

97 MATHEMATICS AND COMPUTING↗

(U) Adjoint-Weighted Leakage Time in SENSMG

It is desired to compute the adjoint weighted leakage time $τ_1$ for an α-eigenvalue problem. The reciprocal of the leakage time is $\frac{1}{τ_1}$ = $\frac{\langle ψ^*, Ω∙∇ψ \rangle}{\langle ψ^*, \frac{1}{v}ψ \rangle}$ where $ψ^*$ and $ψ$ are the adjoint and forward α-eigenfunctions, respectively; $v$ is neutron speed; and the angle brackets indicate an inner product over all independent variables, in this case space, energy, and angle. This memo discusses the implementation of Eq. (1) in the multigroup neutron sensitivity code SENSMG (Refs. 2 and 3), which relies for the neutron transport on the multigroup discrete ordinates code PARTISN (Ref. 4).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) LNK3DNT Capability in SENSMG Extended

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. This report describes the user interface for this capability and the code verification that was done. Reference 5 described an initial version of this capability. Two important extensions have been made: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) Transport Corrections in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections . The PARTISN multigroup discrete ordinates neutron transport code has several transport correction options. These options can be used in the SENSMG multigroup neutron sensitivity code, which relies on PARTISN for the neutron transport. However, with one exception, the sensitivities computed by SENSMG do not include the effect of transport corrections. Thus, when certain transport corrections are used, SENSMG results are inconsistent with the transport solutions. This report discusses the issues surrounding computing sensitivities when transport corrections are used. Problems and possible solutions are identified.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) Transport Corrections Implemented in SENSMG

Truncating the spherical harmonics expansion of the neutron scattering source in the Boltzmann transport equation leads to truncation error. Methods for approximately correcting for this truncation error are called transport corrections. The PARTISN multigroup discrete ordinates neutron transport code has three transport correction options. These options can be specified in the SENSMG multigroup neutron sensitivity code and passed to PARTISN for the neutron transport.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

(U) LNK3DNT Capability in SENSMG Tested with Transport Correction

The SENSMG multigroup neutron sensitivity code can read materials and geometry from redoin and lnk3dnt files written by the PARTISN multigroup discrete ordinates code or the MCNP6 Monte Carlo code. Reference 5 described an initial version of this capability. Reference 6 introduced two important extensions: 1) Reaction-rate ratios are now allowed as responses with redoin/lnk3dnt inputs; and 2) Sensitivities of reaction-rate ratios to mass densities and interface locations are now calculated (for all problems, not just redoin/lnk3dnt problems). In this report, we document the testing of the redoin/lnk3dnt capability when used with transport correction. PARTISN’s transport corrections are described in Ref. 3, and Refs. 7 and 8 discuss the calculation of cross-section sensitivities when transport corrections are applied. References 7 and 8 verify SENSMG for problems with transport correction using central differences. This report verifies the capability with redoin/lnk3dnt inputs by comparison with regular SENSMG inputs. Unlike Refs. 5 and 6, in this report we use a spherical test problem because the sensitivity of the response with respect to the highest-order scattering is zero in cylinders due to the use of flux moments in the inner products.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Recent Updates in ETOE-2 and MC 2 -3

The MC 2 -3 code is a multigroup cross section generation code for fast reactor applications, developed by Argonne National Laboratory (ANL) under the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In this work, first, the cross section libraries for ENDF/B-VII.1 and ENDF/B-VIII.0 data were generated using the updated ETOE-2 code. Preliminary verification and validation tests of the ENDF/B-VII.1 MC 2 -3 library were performed with the selected fast reactor benchmark problems. Processing and verifying the libraries suggested that significant efforts would be required for thoroughly verifying the ENDF/B-VII.1 MC 2 -3 library and successfully processing the ENDF/B-VIII.0 MC 2 -3 library. Secondly, the cross section generation capability of MC 2 -3 was updated with the intermediate group lattice calculation and the equivalent Dancoff-factor cell (EDC) method to significantly improve the performance of a twodimensional assembly calculation using the method of characteristics (MOC). This effort was made useful in implementing and verifying the EDC method in the Griffin cross section API. Finally, we analyzed a SFR problem, for which unusually large deviations in core eigenvalues from Monte Carlo solutions were reported, by generating multigroup cross sections with MC 2 -3, performing core calculations with DIF3D-VARIANT, and analyzing cross section and eigenvalue results against Monte Carlo solutions. The analysis confirmed that the MC 2 -3/DIF3D solutions were in good agreement with Monte Carlo solutions, providing the correct process of accurately generating broad-group cross sections with MC 2 -3

97 MATHEMATICS AND COMPUTING↗

Transport @ LANL: An Overview [Slides]

Deterministic neutron transport requires discretization of the Boltzmann transport equation in space, angle, energy, and time. Discretization in energy is typically performed with the multigroup approximation, which tries to preserve continuous energy cross section information by using a weighted average over discrete energy groups. There are two important requirements for accurate multigroup cross sections: appropriate energy group boundaries and an accurate weight function. The former can be selected with appropriate optimizers; the latter either requires flux solutions from comparable simulations or expert knowledge.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling, Performance Assessment, and Nodal Data Analysis of TRISO-Fueled Systems with Shift

This technical report documents several enhancements to the Shift Monte Carlo (MC) code under the US Department of Energy (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in fiscal year (FY) 2022. Performance enhancements were added to Shift specifically for tristructural isotropic (TRISO)–fueled reactor systems and guided based on performance analysis in FY 2021. For the pebble performance model developed in previous studies, the runtime improved by ~ 91× compared to the original model and ~ 2× compared to the user-optimized model. Compared to Serpent, Shift is ~ 3× slower if Serpent delta-tracking is enabled but ~ 2× faster when delta-tracking is disabled. The multigroup cross section generation was improved through simplifying tally input definitions, porting several post-processing tally operations from Python scripts into the Shift code base, and accounting for production reactions in the scattering multiplicity. Progress was also made on two emerging capabilities: (1) the development of Titan (a Shift reactor physics user interface) and (2) initial investigation into path-length tallies for computing multigroup scattering matrices.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗