• Nuclear Engineering and Technology
• /
• v.53 no.7
• /
• pp.2084-2094
• /
• 2021

Delta-form-based methods for solving high order spatial discretization schemes are introduced into the reactor S_N transport equation. Due to the nature of the delta-form, the final numerical accuracy only depends on the residuals on the right side of the discrete equations and have nothing to do with the parts on the left side. Therefore, various high order spatial discretization methods can be easily adopted for only the transport term on the right side of the discrete equations. Then the simplest step or other robust schemes can be adopted to discretize the increment on the left hand side to ensure the good iterative convergence. The delta-form framework makes the sweeping and iterative strategies of various high order spatial discretization methods be completely the same with those of the traditional S_N codes, only by adding the residuals into the source terms. In this paper, the flux limiter method and weighted essentially non-oscillatory scheme are used for the verification purpose to only show the advantages of the introduction of delta-form-based solving methods and other high order spatial discretization methods can be also easily extended to solve the S_N transport equations. Numerical solutions indicate the correctness and effectiveness of delta-form-based solving method.

• Nuclear Engineering and Technology
• /
• v.55 no.4
• /
• pp.1350-1364
• /
• 2023

2D/1D coupling method is an important neutron transport calculation method due to its high accuracy and relatively low computation cost. However, 2D/1D coupling method may diverge especially in small axial mesh size. To analyze the convergence behavior of 2D/1D coupling method, a Fourier analysis for k-eigenvalue neutron transport problems is implemented. The analysis results present the divergence problem of 2D/1D coupling method in small axial mesh size. Several common attempts are made to solve the divergence problem, which are to increase the number of inner iterations of the 2D or 1D calculation, and two times 1D calculations per outer iteration. However, these attempts only could improve the convergence rate but cannot deal with the divergence problem of 2D/1D coupling method thoroughly. Moreover, the choice of axial solvers, such as DGFEM SN and traditional SN, and its effect on the convergence behavior are also discussed. The results show that the choice of axial solver is a key point for the convergence of 2D/1D method. The DGFEM SN based 2D/1D method could converge within a wide range of optical thickness region, which is superior to that of traditional SN method.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.485-498
• /
• 2020

A diffusion synthetic acceleration (DSA) technique for the S_N transport equation discretized with the linear discontinuous expansion method with subcell balance (LDEM-SCB) on unstructured tetrahedral meshes is presented. The LDEM-SCB scheme solves the transport equation with the discrete ordinates method by using the subcell balances and linear discontinuous expansion of the flux. Discretized DSA equations are derived by consistently discretizing the continuous diffusion equation with the LDEM-SCB method, however, the discretized diffusion equations are not fully consistent with the discretized transport equations. In addition, a fine mesh rebalance (FMR) method is devised to accelerate the discretized diffusion equation coupled with the preconditioned conjugate gradient (CG) method. The DSA method is applied to various test problems to show its effectiveness in speeding up the iterative convergence of the transport equation. The results show that the DSA method gives small spectral radii for the tetrahedral meshes having various minimum aspect ratios even in highly scattering dominant mediums for the homogeneous test problems. The numerical tests for the homogeneous and heterogeneous problems show that DSA with FMR (with preconditioned CG) gives significantly higher speedups and robustness than the one with the Gauss-Seidel-like iteration.

• Nuclear Engineering and Technology
• /
• v.22 no.2
• /
• pp.151-158
• /
• 1990

The new discrete elements method (DEM) is applied to the one-group neutron transport equation in one-dimensional slab geometry. The fixed source and the criticality problems are treated and three spatial differencing schemes (the DD, the SC, -and the LC schemes) are tested to determine the most computationally efficient in the DEM. In all cases, the accuracy of the results obtained from the DEM shows an improvement over that obtained from the standard discrete ordinates calculations. And the LC scheme gives the most accurate results in the DEM.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1165-1171
• /
• 2017

The traditional form of parallelism in Monte Carlo particle transport simulations, wherein each individual particle history is considered a unit of work, does not lend itself well to data-level parallelism. Event-based algorithms, which were originally used for simulations on vector processors, may offer a path toward better utilizing data-level parallelism in modern computer architectures. In this study, a simple model is developed for estimating the efficiency of the event-based particle transport algorithm under two sets of assumptions. Data collected from simulations of four reactor problems using OpenMC was then used in conjunction with the models to calculate the speedup due to vectorization as a function of the size of the particle bank and the vector width. When each event type is assumed to have constant execution time, the achievable speedup is directly related to the particle bank size. We observed that the bank size generally needs to be at least 20 times greater than vector size to achieve vector efficiency greater than 90%. When the execution times for events are allowed to vary, the vector speedup is also limited by differences in the execution time for events being carried out in a single event-iteration.

• Nuclear Engineering and Technology
• /
• v.54 no.9
• /
• pp.3478-3487
• /
• 2022

In this study, we incorporate an anisotropic scattering scheme involving spherical harmonics into the method of characteristics (MOC). The neutron transport solution in a light water reactor can be significantly improved because of the impact of an anisotropic scattering source with the MOC flat source approximation. Several problems are selected to verify the proposed scheme and investigate its effects and accuracy. The MOC anisotropic scattering source is based on the expansion of spherical harmonics with Legendre polynomial functions. The angular flux, scattering source, and cross section are expanded in terms of the surface spherical harmonics. Later, the polynomial is expanded to achieve the odd and even parity of the source components. Ultimately, the MOC angular and scalar fluxes are calculated from a combination of two sources. This paper presents various numerical examples that represent the hot and cold conditions of a reactor core with boron concentration, burnable absorbers, and control rod materials, with and without a reflector or baffle. Moreover, a small critical core problem is considered which involves significant neutron leakage at room temperature. We demonstrate that an anisotropic scattering source significantly improves solution accuracy for the small core high-leakage problem, as well as for practical large core analyses.

• Nuclear Engineering and Technology
• /
• v.55 no.3
• /
• pp.1071-1083
• /
• 2023

The on-the-fly energy release per fission (OTFK) model is implemented in STREAM to continuously update the Kappa values during the depletion calculation. The explicit neutron and photon energy distribution, which has not been considered in previous STREAM versions, is incorporated into the existing on-the-fly model. The impacts of the modified OTFK model with explicit neutron and photon heating in STREAM on the power distribution, fuel temperature, and other core parameters during depletion with feedback calculations are studied using several problems from the VERA benchmark suit. Overall, the explicit heating calculation provides a better power map for the feedback calculations particularly when strong gamma emitters are present. Generally, the fuel temperature decreases when neutron and photon heating is employed because fission neutrons and gamma rays are transported away from their points of generation. This energy release model in STREAM indicates that gamma energy accounts for approximately 9.5%-10% of the total energy released, and approximately 2.4%-2.6% of the total energy released will be deposited in the coolant for the VERA 5, NuScale, and Yonggwang Unit 3 2D cores.

• Nuclear Engineering and Technology
• /
• v.24 no.3
• /
• pp.319-328
• /
• 1992

The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation.

• Nuclear Engineering and Technology
• /
• v.49 no.6
• /
• pp.1240-1249
• /
• 2017

An essential component of the neutron transport solver is the resonance self-shielding calculation used to determine equivalence cross sections. The neutron transport code, MPACT, is currently using the subgroup self-shielding method, in which the method of characteristics (MOC) is used to solve purely absorbing fixed-source problems. Recent efforts incorporating multigroup kernels to the MOC solvers in MPACT have reduced runtime by roughly $2{\times}$. Applying the same concepts for self-shielding and developing a novel lumped parameter approach to MOC, substantial improvements have also been made to the self-shielding computational efficiency without sacrificing any accuracy. These new multigroup and lumped parameter capabilities have been demonstrated on two test cases: (1) a single lattice with quarter symmetry known as VERA (Virtual Environment for Reactor Applications) Progression Problem 2a and (2) a two-dimensional quarter-core slice known as Problem 5a-2D. From these cases, self-shielding computational time was reduced by roughly $3-4{\times}$, with a corresponding 15-20% increase in overall memory burden. An azimuthal angle sensitivity study also shows that only half as many angles are needed, yielding an additional speedup of $2{\times}$. In total, the improvements yield roughly a $7-8{\times}$ speedup. Given these performance benefits, these approaches have been adopted as the default in MPACT.

• Nuclear Engineering and Technology
• /
• v.56 no.2
• /
• pp.673-684
• /
• 2024

STREAM, developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST), is a deterministic neutron- and photon-transport code primarily designed for light water reactor (LWR) analysis. Initially, the photon module in STREAM did not account for fluorescence and bremsstrahlung photons. This article presents recent developments regarding the integration of atomic relaxation and bremsstrahlung models into the existing photon module, thus allowing for the transport of secondary photons. The photon flux and photon heating computed with the newly incorporated models is compared to results obtained with the Monte Carlo code MCS. The incorporation of secondary photons has substantially improved the accuracy of photon flux calculations, particularly in scenarios involving strong gamma emitters. However, it is essential to note that despite the consideration of secondary photon sources, there is no noticeable improvement in the photon heating for LWR problems when compared to the photon heating obtained with the previous version of STREAM.