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 397 records · Page 22

Performance of the block-Krylov energy group solvers in Jaguar

A new method of coupling the inner and outer iterations for deterministic transport problems is proposed. This method is termed the Multigroup Energy Blocking Method (MEBM) and has been implemented in the deterministic transport solver Jaguar, which is currently under development at KAPL. The method is derived for both fixed-source and eigenvalue problems. The method is then applied to a PWR pin cell model, both in fixed-source mode and eigenvalue mode. The results show that the MEBM improves the convergence of both types of problems when applied to the thermal (up-scattering) groups. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

HTTR criticality calculations with SCALE6: Studies of various geometric and unit-cell options in modeling

The fuel element of the High Temperature Engineering Test Reactor (HTTR) presents a doubly heterogeneous geometry, where tiny TRISO fuel particles dispersed in a graphite matrix form the fuel region of a cylindrical fuel rod, and a number of fuel rods together with moderator or reflector then constitute the lattice design of the core. In this study, a series of full-core HTTR criticality calculations were performed with the SCALE6 code system using various geometric and unit-cell options in order to systematically investigate their effects on neutronic analysis. Two geometric descriptions (ARRAY or HOLE) in SCALE6 can be used to construct a complicated and repeated model. The result shows that eliminating the use of HOLE in the HTTR geometric model can save the computation time by a factor of 4. Four unit-cell treatments for resonance self-shielding corrections in SCALE6 were tested to create problem-specific multigroup cross sections for the HTTR core model. Based on the same ENDF/B-VII cross-section library, their results were evaluated by comparing with continuous-energy calculations. The comparison indicates that the INFHOMMEDIUM result overestimates the system multiplication factor (k{sub eff}) by 55 mk, whereas the LATTICECELL and MULTIREGION treatments predict the k{sub eff} values with similar biases of approximately 10 mk overestimation. The DOUBLEHET result shows a more satisfactory agreement, about 4.2 mk underestimation in the k{sub eff} value. In addition, using cell-weighted cross sections instead of an explicit modeling of TRISO particles in fuel region can further reduce the computation time by a factor of 5 without sacrificing accuracy. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development of a statistical sampling method for uncertainty analysis with SCALE

A new statistical sampling sequence called Sampler has been developed for the SCALE code system. Random values for the input multigroup cross sections are determined by using the XSUSA program to sample uncertainty data provided in the SCALE covariance library. Using these samples, Sampler computes perturbed self-shielded cross sections and propagates the perturbed nuclear data through any specified SCALE analysis sequence, including those for criticality safety, lattice physics with depletion, and shielding calculations. Statistical analysis of the output distributions provides uncertainties and correlations in the desired responses. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hybrid parallel code acceleration methods in full-core reactor physics calculations

When dealing with nuclear reactor calculation schemes, the need for three dimensional (3D) transport-based reference solutions is essential for both validation and optimization purposes. Considering a benchmark problem, this work investigates the potential of discrete ordinates (Sn) transport methods applied to 3D pressurized water reactor (PWR) full-core calculations. First, the benchmark problem is described. It involves a pin-by-pin description of a 3D PWR first core, and uses a 8-group cross-section library prepared with the DRAGON cell code. Then, a convergence analysis is performed using the PENTRAN parallel Sn Cartesian code. It discusses the spatial refinement and the associated angular quadrature required to properly describe the problem physics. It also shows that initializing the Sn solution with the EDF SPN solver COCAGNE reduces the number of iterations required to converge by nearly a factor of 6. Using a best estimate model, PENTRAN results are then compared to multigroup Monte Carlo results obtained with the MCNP5 code. Good consistency is observed between the two methods (Sn and Monte Carlo), with discrepancies that are less than 25 pcm for the k{sub eff}, and less than 2.1% and 1.6% for the flux at the pin-cell level and for the pin-power distribution, respectively. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Verification of the shift Monte Carlo code with the C5G7 reactor benchmark

Shift is a new hybrid Monte Carlo/deterministic radiation transport code being developed at Oak Ridge National Laboratory. At its current stage of development, Shift includes a parallel Monte Carlo capability for simulating eigenvalue and fixed-source multigroup transport problems. This paper focuses on recent efforts to verify Shift's Monte Carlo component using the two-dimensional and three-dimensional C5G7 NEA benchmark problems. Comparisons were made between the benchmark eigenvalues and those output by the Shift code. In addition, mesh-based scalar flux tally results generated by Shift were compared to those obtained using MCNP5 on an identical model and tally grid. The Shift-generated eigenvalues were within three standard deviations of the benchmark and MCNP5-1.60 values in all cases. The flux tallies generated by Shift were found to be in very good agreement with those from MCNP. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Potential of pin-by-pin SPN calculations as an industrial reference

This paper aims at analysing the potential of pin-by-pin SP{sub n} calculations to compute the neutronic flux in PWR cores as an alternative to the diffusion approximation. As far as pin-by-pin calculations are concerned, a SPH equivalence is used to preserve the reactions rates. The use of SPH equivalence is a common practice in core diffusion calculations. In this paper, a methodology to generalize the equivalence procedure in the SP{sub n} equations context is presented. In order to verify and validate the equivalence procedure, SP{sub n} calculations are compared to 2D transport reference results obtained with the APOLL02 code. The validation cases consist in 3x3 analytical assembly color sets involving burn-up heterogeneities, UOX/MOX interfaces, and control rods. Considering various energy discretizations (up to 26 groups) and flux development orders (up to 7) for the SP{sub n} equations, results show that 26-group SP{sub 3} calculations are very close to the transport reference (with pin production rates discrepancies < 1%). This proves the high interest of pin-by-pin SP{sub n} calculations as an industrial reference when relying on 26 energy groups combined with SP{sub 3} flux development order. Additionally, the SP{sub n} results are compared to diffusion pin-by-pin calculations, in order to evaluate the potential benefit of using a SP{sub n} solver as an alternative to diffusion. Discrepancies on pin-production rates are less than 1.6% for 6-group SP{sub 3} calculations against 3.2% for 2-group diffusion calculations. This shows that SP{sub n} solvers may be considered as an alternative to multigroup diffusion. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An improved resonance self-shielding method for heterogeneous fast reactor assembly and core calculations

An improved resonance self-shielding method has been developed to accurately estimate the effective multigroup cross sections for heterogeneous fast reactor assembly and core calculations. In the method, the heterogeneity effect is considered by the use of isotopic escape cross sections while the resonance interference effect is accounted for through the narrow resonance approximation or slowing-down calculations for specific compositions. The isotopic escape cross sections are calculated by solving fixed-source transport equations with the method of characteristics for the whole problem domain. This method requires no pre-calculated resonance integral tables or parameters that are typically necessary in the subgroup method. Preliminary results for multi pin-cell fast reactor problems show that the escape cross sections estimated from the explicit-geometry fixed source calculations produce more accurate eigenvalue and self-shielded effective cross sections than those from conventional one-dimensional geometry models. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupling the core analysis program DeCART to the fuel performance application BISON

The 3D neutron transport and core analysis program DeCART was coupled to the fuels performance application BISON to provide a higher fidelity tool for fuel performance simulation. This project is motivated by the desire to couple a high fidelity core analysis program (based on the method of characteristics) to a high fidelity fuel performance program, both of which can simulate 3D problems. DeCART provides sub-pin level resolution of the multigroup neutron flux, with resonance treatment, during burnup or a fast transient. BISON implicitly solves coupled thermomechanical equations for the fuel on a sub-millimeter level finite element mesh. A method was developed for mapping the fission rate density and fast neutron flux from DeCART to BISON. Multiple depletion cases were run with one-way data transfer from DeCART to BISON. The one-way data transfer of fission rate density is shown to agree with the fission rate density obtained from an internal Lassman-style model in BISON. One-way data transfer was also demonstrated in a 3D case in which azimuthal asymmetry was induced in the fission rate density profile of a fuel rod modeled in DeCART. Two-way data transfer was established by mapping the temperature distribution from BISON to DeCART. A Picard iterative algorithm was developed for the loose coupling with two-way data transfer. (authors)

97 MATHEMATICS AND COMPUTING↗

Spatial and spectral effects in subcritical system pulsed experiments

Accurate neutronic models are needed for the interpretation of pulsed experiments in subcritical systems. In this work, the extent of spatial and spectral effects in the pulse propagation phenomena is investigated and the analysis is applied to the GUINEVERE experiment. The multigroup cross section data is generated by the Monte Carlo SERPENT code and the neutronic evolution following the source pulse is simulated by a kinetic diffusion code. The results presented show that important spatial and spectral aspects need to be properly accounted for and that a detailed energy approach may be needed to adequately capture the physical features of the system to the pulse injection. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Improvements of the boundary projection acceleration technique applied to the discrete-ordinates transport solver in XYZ geometries

In this paper we discuss the latest upgrades of the Boundary Projection Acceleration (BPA) applied to the XYZ transport solver of APOLLO3, namely IDT. The acceleration method is a well-known effective technique for the speed-up of the source iterations of the discrete-ordinates method. The BPA in IDT has been improved in three aspects: the taking into account of the residue on boundary conditions as a boundary source for the acceleration problem, the extension of the method to higher order angular moments in the case of anisotropic scattering and, finally, the application of the method to the multigroup iterations for the acceleration of the fission source and k-effective. The spectrum of the method has been Fourier-analyzed to explore the effectiveness. The 3D mock-up geometry of the ZPPR is presented as final study to test the performances of the acceleration on a realistic whole-core 3D calculation. (authors)

97 MATHEMATICS AND COMPUTING↗

Calculating infinite-medium {alpha}-eigenvalue spectra with a transition rate matrix method

The time-dependent behavior of the energy spectrum in neutron transport was investigated with a formulation, based on continuous-time Markov processes, for computing {alpha}-eigenvalues and eigenvectors in an infinite medium. For this, a research Monte Carlo code called TORTE was created and used to estimate elements of a transition rate matrix. TORTE is capable of using both multigroup and continuous-energy nuclear data, and verification was performed. Eigenvalue spectra for infinite homogeneous mixtures were obtained and an eigenfunction expansion was used to investigate transient behavior of the neutron energy spectrum. (authors)

97 MATHEMATICS AND COMPUTING↗

ENDF/B-VIII.0 testing with AMPX and SCALE

As the new ENDF/B-VIII.0 data has continued development from its initial beta release into the final release version, it has subsequently undergone processing and testing at Oak Ridge National Laboratory with AMPX and SCALE. Each of the beta releases has been processed with AMPX to create both multigroup (MG) and continuous-energy (CE) libraries. The MG libraries were all processed using the SCALE 252-group structure. These libraries were then used with the VALID suite for benchmarking the new data library's performance against the currently distributed ENDF/B-VII.1 (VII.1)-based SCALE libraries. Here, we compare the MG and CE VALID results for the two current ENDF releases, versions VII.1 and the most recent beta release of ENDF/B-VIII.0 as of this writing (VIII.0.β4). Generally, both the MG and CE VALID aggregate results for VIII.0.β4 compare reasonably well with their VII.1-derived counterparts. As can be expected of any new ENDF release, some of the individual cases have substantially different computed neutron multiplicity (k{sub eff}) values between the two releases, although typically the differences are less than 100 percent-milli (pcm). However, changes to the plutonium cross sections have induced much larger changes in the ratio of the computed k{sub eff} to the experimental value (C/E) for thermal spectrum plutonium solution critical systems. This paper details the comparative performance of the VII.1 and VIII.0.β4 data for the broad range of test cases covered by the VALID suite. A companion paper at this conference provides details of the baseline validation of SCALE 6.2, including the first release of the VII.1 data. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

An Atomic Mix Closure for Stochastic Media Transport Problems

The accurate modeling of stochastic media in radiation transport calculations remains on ongoing research challenge. The Levermore-Pomraning (LP) closure for the statistical transport equation is simple and inexpensive, but it is often inaccurate. More accurate closures have been proposed, but they can be prohibitively expensive. In this work we present a relatively inexpensive closure that we have observed to be more accurate than LP in almost all cases studied. Obtaining solutions to stochastic transport problems can be quite difficult. Transport calculations on an ensemble of explicit realizations generated by Monte Carlo sampling can be prohibitively expensive. Atomic mix or LP calculations are relatively inexpensive, but they can be inaccurate. The present work illustrates an approach that is less expensive than Monte Carlo sampling yet more accurate than atomic mix or LP approximations. We still need to analyze the reasons why the atomic mix closures derived for thicker problems lead to source iteration instabilities. It is not clear if the closure itself induces an ill-posed problem or if it is merely the iterative process that is problematic. We hope to examine alternative source terms to drive the subsidiary calculations. We also want to extend the work to multigroup and multidimensional problems. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nuclear Data Uncertainty Quantification for Adjoint-Weighted Quantities via XGPT in Monte Carlo Codes

The Total Monte Carlo method (TMC) is an innovative approach to nuclear data uncertainty propagation. In TMC, a large number of independent ENDF files are randomly generated starting from nuclear models parameters and their uncertainties. These files are then processed with NJOY or similar codes to produce a set of formatted continuous-energy cross sections (i.e., ACE files), and independent Monte Carlo neutron transport simulations are run with these different ACE files as input. Finally the distributions of the considered response functions are obtained directly from the results of the Monte Carlo runs. Despite several improvements of the TMC method, this approach is often considered an inefficient way of propagating nuclear data uncertainties, as it requires several computationally expensive Monte Carlo runs. Covariance-based methods present important advantages for several applications, especially when different sources of uncertainties are to be investigated. For this reason, sensitivity analysis and uncertainty quantification capabilities are being implemented in several Monte Carlo neutron transport codes. On the other hand, the fact that no artificial multi-group covariance matrices are required and that higher moments of the responses distributions can be easily obtained, are definite advantages of the TMC method. In an effort to overcome the limitations of both approaches, a new method, called eXtended Generalized Perturbation Theory (XGPT) is under development. When performing nuclear data uncertainty propagation, such XGPT approach does not require generating discretized (i.e., multigroup) covariance matrices, and allows to estimate higher moments of the response probability distributions, via linear perturbation theory. One of the promising applications of the XGPT method involves its adoption in conjunction with the Proper Orthogonal Decomposition (POD) of nuclear data, to build a Reduced Order Model (ROM) for the approximation of thousands of independent TMC results in a single criticality source simulation. In the present work, this approach is employed for the first time to evaluate the impact of uncertainty in nuclear data on adjoint-weighted quantities. The method has been implemented in an extended version of the Serpent-2.1.24 Monte Carlo code and is demonstrated on two fast spectrum benchmarks, focusing on {sup 208}Pb nuclear data. A new method for propagating nuclear data uncertainties in criticality source Monte Carlo simulations has been developed. This eXtended Generalized Perturbation Theory (XGPT) approach allows to accurately reproduce reference Total Monte Carlo (TMC) uncertainty estimates in a single criticality source simulation. The present work aims to demonstrate that the proposed method can quantify nuclear data uncertainty based on adjoint-weighted quantities. The XGPT approach was applied to evaluate the impact of {sup 208}Pb nuclear data uncertainties on the effective prompt neutron lifetime (l{sub eff}) in two fast spectrum systems. The obtained XGPT results show a good agreement with the reference TMC results, obtained from independent Monte Carlo runs with {sup 208}Pb nuclear data from random TENDL-2013 evaluations. The XGPT method, implemented in a extended version of the Serpent Monte Carlo code and adopted along with the Proper Orthogonal Decomposition (POD) of the nuclear data, is able to provide an accurate reduced order model to approximate thousands TMC results within a single run, ensuring a large saving in computational requirements. The effect of {sup 208}Pb nuclear data uncertainties appears to be rather large. This is consistent with previous studies that adopted the same input random evaluations data from TENDL. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Finite Volume Method Based Neutronic Solvers for Steady and Transient Analysis of Molten Salt Reactors

In recent years, the coupled neutronics and thermal hydraulics analysis becomes more and more important, because it helps to identify the most relevant safety issues without conservative assumptions. In addition, such coupled simulations are essential for the generation IV advanced nuclear power systems such as molten salt reactor (MSR), because the presence of strong coupling between different phenomena, like the transport of delayed neutron precursors by the fluid flow?. Due to the different discretization methods used in reactor physics analysis and thermal-hydraulics analysis, it is hard to couple different kinds of codes. This work presents an attempt to solve the multigroup neutron diffusion equation with the finite volume method which is quite widely used for the numerical simulation of a variety of applications involving fluid flow and heat and mass transfer. The unstructured finite volume method (FVM) based on open source C++ library OpenFOAM was chosen as the development platform?. The objective of this paper is to concisely present the main features of the developed codes in terms of the theoretical background and numerical approach and as can be seen in the following sections, the calculation results obtained by the NDSFoam which solves the steady-state neutron diffusion equation and NDTFoam which solves the transient-state neutron diffusion equation are discussed and compared to the benchmark problems. In this paper, two new neutronic solvers using finite volume method have been developed to perform the neutron physics analysis. The solvers are based on OpenFOAM, and have been tested with several steady state and transient-state benchmark problems. By verifications of the codes, it is proved that the development of NDSFoam and NDTFoam are correct. In our future work, both of these two solvers will be coupled with the precompiled CFD solvers in OpenFOAM to analyze the molten salt fast reactors. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

use of cross-section covariance information in uncertainty analysis at modeling HTGR neutron physics - 18551

This work is aimed at improvement of high temperature gas cooled reactors (HTGR) neutron physics design calculations by application of uncertainty analysis with the use of cross-section covariance information. Methodology and codes for preparation of multigroup libraries of covariance information for individual isotopes from the basic 44-group binary library of SCALE-6 code system were developed. A 69-group library of covariance information in a special format for main isotopes and elements typical for HTGR was generated. This library can be used for estimation of uncertainties associated with nuclear data in analysis of HTGR neutron physics with design codes. As an example, calculations of one-group cross-section uncertainties for fission and capture reactions for main isotopes of the MHTGR-350 benchmark, as well as uncertainties of the multiplication factor (k{sub ∞}) for the MHTGR-350 fuel compact cell model and fuel block model were performed. These uncertainties were estimated by the developed technology with the use of WIMS-D code (technology SUW) and modules of SCALE-6 code system, namely, by TSUNAMI, KENO-VI and SAMS. The eight most important reactions on isotopes for MHTGR-350 benchmark were identified, namely: {sup 10}B(capt), {sup 238}U(n,γ), ν{sub 5}, {sup 235}U(n,γ), {sup 238}U(el), {sup nat}C(el), {sup 235}U(fiss)-{sup 235}U(n,γ) and {sup 235}U(fiss). It is shown that the main contribution in the uncertainty of multiplication factor of the fuel block model gives the reaction {sup 10}B(capt) ((n,α)+(n,γ)). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Localized Radial Basis Functions for Radiation Transport in Slab Geometry

Meshing the spatial domain for the neutron transport equation can be computationally expensive and labor-intensive. Mesh-free methods, such as radial basis function (RBF) approximations, may alleviate these concerns by using a pointwise solution of the transport equation instead of the standard volume-integrated solution using cells or elements. This allows for the use of any geometric description that can provide the material at solution points and the boundary normal directions for boundary points. RBFs give spectral convergence when used to solve differential equations, but only if they are infinitely smooth. This occurs when the RBFs are effectively flat, creating an ill-conditioned, dense system of equations. One simple way to alleviate this issue is to use localized RBFs. Global RBFs have been applied previously to the monoenergetic radiative transport equation in X-Y geometry using source iteration. In this summary, localized RBFs are applied to the multigroup neutron transport equation in slab geometry using Krylov iteration. The results show agreement with benchmark solutions, provided the coefficients for the RBF basis are properly chosen. (authors)

97 MATHEMATICS AND COMPUTING↗

Performance Improvements to the Cross Section Calculation in MPACT

The MPACT code is the primary deterministic neutron transport solver available in the Virtual Environment for Reactor Applications (VERA) as part of the Consortium for Advanced Simulation of Light Water Reactors (CASL). In contrast to the conventional two-step method using homogenization techniques, MPACT employs a direct transport methodology, where the transport calculation is performed with the explicit reactor core geometry, material compositions, and temperature distribution. As the direct whole-core transport calculation requires significant computational resources, performance improvements have been recently conducted in various components of MPACT. The initial MPACT profiling was performed on the EOS cluster at the Oak Ridge Leadership Computing Facility. Table I shows that the cross section processing accounted for 74% and 33% of the total runtime for VERA Progression Problems 5 and 7. Starting with a problem-independent multigroup (MG) library, two primary functionalities are performed in MPACT to obtain the problem-dependent MG cross sections: (1) the resonance self-shielding calculation, which currently MPACT is approaching with the subgroup method; (2) use of equivalence cross section from resonance calculation to compute the microscopic cross sections, and thus, the macroscopic cross sections of material regions. As shown in Table I, the time distribution in cross section processing varies from case to case. For Problem 5, a 2-D quarter core depletion case, the macroscopic XS dominates due to the increased number of isotopes in the depleted fuel. For Problem 7, a 3-D quarter core PWR at HFP and BOL, more subgroup calculations are performed to obtain updated self-shielded cross sections during the iteration. In this summary, the approaches to reduce the computing time in cross section processing are presented. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗