• Nuclear Engineering and Technology
• /
• 제49권6호
• /
• pp.1326-1338
• /
• 2017

The continuous energy Monte Carlo neutron transport code, MC21, was coupled to the CTF subchannel thermal-hydraulics code using a combination of Consortium for Advanced Simulation of Light Water Reactors (CASL) tools and in-house Python scripts. An MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 6 demonstrated good agreement with MC21/COBRA-IE and VERA solutions. The MC21/CTF solution for VERA Core Physics Benchmark Progression Problem 7, Watts Bar Unit 1 at beginning of cycle hot full power equilibrium xenon conditions, is the first published coupled Monte Carlo neutronics/subchannel T-H solution for this problem. MC21/CTF predicted a critical boron concentration of 854.5 ppm, yielding a critical eigenvalue of $0.99994{\pm}6.8E-6$ (95% confidence interval). Excellent agreement with a VERA solution of Problem 7 was also demonstrated for integral and local power and temperature parameters.

• Nuclear Engineering and Technology
• /
• 제55권9호
• /
• pp.3388-3400
• /
• 2023

Detailed analysis of the neutron pathway through matter inside the nuclear reactor core is exceedingly needed for safety and economic considerations. Due to the constant development of high-performance computing technologies, neutronics analysis using computer codes became more effective and efficient to perform sophisticated neutronics calculations. In this work, a commercial pressurized water reactor (PWR) presented by Virtual Environment for Reactor Applications (VERA) Core Physics Benchmark are modeled and simulated using a high-fidelity simulation of OpenMC code in terms of criticality and fuel pin power distribution. Various problems have been selected from VERA benchmark ranging from a simple two-dimension (2D) pin cell problem to a complex three dimension (3D) full core problem. The development of the code capabilities for reactor physics methods has been implemented to investigate the accuracy and performance of the OpenMC code against VERA SCALE codes. The results of OpenMC code exhibit excellent agreement with VERA results with maximum Root Mean Square Error (RMSE) values of less than 0.04% and 1.3% for the criticality eigenvalues and pin power distributions, respectively. This demonstrates the successful utilization of the OpenMC code as a simulation tool for a whole core analysis. Further works are undergoing on the accuracy of OpenMC simulations for the impact of different fuel types and burnup levels and the analysis of the transient behavior and coupled thermal hydraulic feedback.

• Nuclear Engineering and Technology
• /
• 제54권7호
• /
• pp.2670-2689
• /
• 2022

STREAM - a lattice transport calculation code with method of characteristics for the purpose of light water reactor analysis - has been developed by the Computational Reactor Physics and Experiment laboratory (CORE) of the Ulsan National Institute of Science and Technology (UNIST). Recently, efforts have been taken to develop a photon module in STREAM to assess photon heating and the influence of gamma photon transport on power distributions, as only neutron transport was considered in previous STREAM versions. A multi-group photon library is produced for STREAM based on the ENDF/B-VII.1 library with the use of the library-processing code NJOY. The developed photon solver for the computation of 2D and 3D distributions of photon flux and energy deposition is based on the method of characteristics like the neutron solver. The photon library and photon module produced and implemented for STREAM are verified on VERA pin and assembly problems by comparison with the Monte Carlo code MCS - also developed at UNIST. A short analysis of the impact of photon transport during depletion and thermal hydraulics feedback is presented for a 2D core also from the VERA benchmark.

• Nuclear Engineering and Technology
• /
• 제53권11호
• /
• pp.3543-3562
• /
• 2021

The methods and performance of a 3D pin-by-pin neutronics code based on the 2D/1D decoupling method are presented. The code was newly developed as an effort to achieve enhanced accuracy and high calculation performance that are sufficient for the use in practical nuclear design analyses. From the 3D diffusion-based finite difference method (FDM) formulation, decoupled planar formulations are established by treating pre-determined axial leakage as a source term. The decoupled axial problems are formulated with the radial leakage source term. To accelerate the pin-by-pin calculation, the two-level coarse mesh finite difference (CMFD) formulation, which consists of the multigroup node-wise CMFD and the two-group assembly-wise CMFD is implemented. To enhance the accuracy, both the discontinuity factor method and the super-homogenization (SPH) factor method are examined for pin-wise cross-section homogenization. The parallelization is achieved with the OpenMP package. The accuracy and performance of the pin-by-pin calculations are assessed with the VERA and APR1400 benchmark problems. It is demonstrated that pin-by-pin 2D/1D alternating calculations within the two-level 3D CMFD framework yield accurate solutions in about 30 s for the typical commercial core problems, on a parallel platform employing 32 threads.

• Nuclear Engineering and Technology
• /
• 제50권7호
• /
• pp.1043-1050
• /
• 2018

For an efficient Monte Carlo (MC) burnup analysis, an accurate high-order depletion scheme to consider the nonlinear flux variation in a coarse burnup-step interval is crucial accompanied with an accurate depletion equation solver. In a Seoul National University MC code, McCARD, the high-order depletion schemes of the quadratic depletion method (QDM) and the linear extrapolation/quadratic interpolation (LEQI) method and a depletion equation solver by the Chebyshev rational approximation method (CRAM) have been newly implemented in addition to the existing constant extrapolation/backward extrapolation (CEBE) method using the matrix exponential method (MEM) solver with substeps. In this paper, the quadratic extrapolation/quadratic interpolation (QEQI) method is proposed as a new high-order depletion scheme. In order to examine the effectiveness of the newly-implemented depletion modules in McCARD, four problems in the VERA depletion benchmarks are solved by CEBE/MEM, CEBE/CRAM, LEQI/MEM, QEQI/MEM, and QDM for gadolinium isotopes. From the comparisons, it is shown that the QEQI/MEM predicts ${k_{inf}}^{\prime}s$ most accurately among the test cases. In addition, statistical uncertainty propagation analyses for a VERA pin cell problem are conducted by the sensitivity and uncertainty and the stochastic sampling methods.

• Nuclear Engineering and Technology
• /
• 제52권5호
• /
• pp.878-888
• /
• 2020

This paper presents the validation of UNIST in-house Monte Carlo code MCS used for the high-fidelity simulation of commercial pressurized water reactors (PWRs). Its focus is on the accurate, spatially detailed neutronic analyses of startup physics tests for the initial core of the Watts Bar Nuclear 1 reactor, which is a vital step in evaluating core phenomena in an operating nuclear power reactor. The MCS solutions for the Consortium for Advanced Simulation of Light Water Reactors (CASL) Virtual Environment for Reactor Applications (VERA) core physics benchmark progression problems 1 to 5 were verified with KENO-VI and Serpent 2 solutions for geometries ranging from a single-pin cell to a full core. MCS was also validated by comparing with results of reactor zero-power physics tests in a full-core simulation. MCS exhibits an excellent consistency against the measured data with a bias of ±3 pcm at the initial criticality whole-core problem. Furthermore, MCS solutions for rod worth are consistent with measured data, and reasonable agreement is obtained for the isothermal temperature coefficient and soluble boron worth. This favorable comparison with measured parameters exhibited by MCS continues to broaden its validation basis. These results provide confidence in MCS's capability in high-fidelity calculations for practical PWR cores.

• Nuclear Engineering and Technology
• /
• 제55권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
• /
• 제53권1호
• /
• pp.30-43
• /
• 2021

In the 2D/1D method, a global adjoint CMFD based on the generalized equivalence theory is built to synthesize the 2D radial MOC adjoint and 1D axial NEM adjoint calculation and also to accelerate the iteration convergence of 3D whole-core adjoint transport calculation. Even more important, an advanced yet accurate two-level (TL) CMFD acceleration technique is proposed, in which an equivalent one-group adjoint CMFD is established to accelerate the multi-group adjoint CMFD and then to accelerate the 3D whole-core adjoint transport calculation efficiently. Based on these method, a new code is developed to perform 3D adjoint neutron flux calculation. Then a set of VERA and C5G7 benchmark problems are chosen to verify the capability of the 3D adjoint calculations and the effectiveness of TL CMFD acceleration. The numerical results demonstrate that acceptable accuracy of 2D/1D adjoint calculations and superior acceleration of TL CMFD are achievable.