• Nuclear Engineering and Technology
• /
• 제46권3호
• /
• pp.313-342
• /
• 2014

A Nuclear Energy Agency (NEA), Organization for Economic Co-operation and Development (OECD) benchmark for Uncertainty Analysis in Modeling (UAM) is defined in order to facilitate the development and validation of available uncertainty analysis and sensitivity analysis methods for best-estimate Light water Reactor (LWR) design and safety calculations. The benchmark has been named the OECD/NEA UAM-LWR benchmark, and has been divided into three phases each of which focuses on a different portion of the uncertainty propagation in LWR multi-physics and multi-scale analysis. Several different reactor cases are modeled at various phases of a reactor calculation. This paper discusses Phase I, known as the "Neutronics Phase", which is devoted mostly to the propagation of nuclear data (cross-section) uncertainty throughout steady-state stand-alone neutronics core calculations. Three reactor systems (for which design, operation and measured data are available) are rigorously studied in this benchmark: Peach Bottom Unit 2 BWR, Three Mile Island Unit 1 PWR, and VVER-1000 Kozloduy-6/Kalinin-3. Additional measured data is analyzed such as the KRITZ LEU criticality experiments and the SNEAK-7A and 7B experiments of the Karlsruhe Fast Critical Facility. Analyzed results include the top five neutron-nuclide reactions, which contribute the most to the prediction uncertainty in keff, as well as the uncertainty in key parameters of neutronics analysis such as microscopic and macroscopic cross-sections, six-group decay constants, assembly discontinuity factors, and axial and radial core power distributions. Conclusions are drawn regarding where further studies should be done to reduce uncertainties in key nuclide reaction uncertainties (i.e.: $^{238}U$ radiative capture and inelastic scattering (n, n') as well as the average number of neutrons released per fission event of $^{239}Pu$).

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2004년도 추계학술발표회 발표논문집
• /
• pp.75-76
• /
• 2004

• Nuclear Engineering and Technology
• /
• 제52권11호
• /
• pp.2422-2430
• /
• 2020

In order to improve calculational efficiency of the CAPP code in the analysis of the hexagonal reactor core, we have tried to implement a refined AFEN method with transverse gradient basis functions and interface flux moments in the hexagonal geometry. The numerical scheme for the refined AFEN method adopted here is the response matrix method that uses the interface partial currents as nodal unknowns instead of the interface fluxes used in the original AFEN method. Since the response matrix method is single-node based, it has good properties such as good calculational efficiency and parallel computing affinity. Because a refined AFEN method equivalent nonlinear FDM response matrix method tried first could not provide a numerically stable solution, a direct formulation of the refined AFEN response matrix were developed. To show the numerical performance of this response matrix method against the original AFEN method, the numerical error analyses were performed for several benchmark problems including the VVER-440 LWR benchmark problem and the MHTGR-350 HTGR benchmark problem. The results showed a more than three times speedup in computing time for the LWR and HTGR benchmark problems due to good convergence and excellent calculational efficiency of the refined AFEN response matrix method.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1999년도 춘계학술발표회요약집
• /
• pp.24-24
• /
• 1999

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 2005년도 춘계학술발표회
• /
• pp.147-148
• /
• 2005

• Nuclear Engineering and Technology
• /
• 제50권7호
• /
• pp.1024-1036
• /
• 2018

In this paper, McCARD code was verified using various models listed in the NUREG/CR-6361 benchmark guide, which provides specifications for single pin-cells, single assemblies, and the whole core classified depending on the nuclear properties and structural characteristics. McCARD code was verified by comparing its results with those of SCALE code for single pin-cell and single assembly benchmark problems. The difference in the multiplication factor obtained through the two codes did not exceed 90 pcm. The benchmark guide treats a total of 173 whole core experiments. The experiments are categorized as simple lattices, separator plates, reflecting walls, reflecting walls and separator plates, burnable absorber fuel rods, water holes, poison rods, and borated moderator. As a result of numerical simulation using McCARD, the mean value of the multiplication factors is 1.00223 and the standard deviation of the multiplication factors is 285 pcm. The difference between the multiplication factors and the experimental value is in the range of -665 pcm to + 1609 pcm. In addition, statistics of results for experiments categorized by reactor shape, additional structure, burnable poison, etc., are detailed in the main text.

• Nuclear Engineering and Technology
• /
• 제46권3호
• /
• pp.343-352
• /
• 2014

Uncertainty and sensitivity analyses with respect to nuclear data are performed with depletion calculations for BWR and PWR fuel assemblies specified in the framework of the UAM-LWR Benchmark Phase II. For this, the GRS sampling based tool XSUSA is employed together with the TRITON depletion sequences from the SCALE 6.1 code system. Uncertainties for multiplication factors and nuclide inventories are determined, as well as the main contributors to these result uncertainties by calculating importance indicators. The corresponding neutron transport calculations are performed with the deterministic discrete-ordinates code NEWT. In addition, the Monte Carlo code KENO in multi-group mode is used to demonstrate a method with which the number of neutron histories per calculation run can be substantially reduced as compared to that in a calculation for the nominal case without uncertainties, while uncertainties and sensitivities are obtained with almost the same accuracy.

• Nuclear Engineering and Technology
• /
• 제52권8호
• /
• pp.1626-1637
• /
• 2020

This work covers an important point of the benchmark released by the expert group on Uncertainty Analysis in Modeling of Light Water Reactors. This ambitious benchmark aims to determine the uncertainty in light water reactors systems and processes in all stages of calculation, with emphasis on multi-physics (coupled) and multi-scale simulations. The Gesellschaft für Anlagen und Reaktorsicherheit methodology is used to propagate the thermal-hydraulic uncertainty of macroscopic parameters through TRACE5.0p3/PARCSv3.0 coupled code. The main innovative points achieved in this work are i) a new thermal-hydraulic model is developed with a highly-accurate 3D core discretization plus an iterative process is presented to adjust the 3D bypass flow, ii) a control rod insertion occurrence -which data is obtained from a real PWR test- is used as a transient simulation, iii) two approaches are used for the propagation process: maximum response where the uncertainty and sensitivity analysis is performed for the maximum absolute response and index dependent where the uncertainty and sensitivity analysis is performed at each time step, and iv) RESTING MATLAB code is developed to automate the model generation process and, then, propagate the thermal-hydraulic uncertainty. The input uncertainty information is found in related literature or, if not found, defined based on expert judgment. This paper, first, presents the Gesellschaft für Anlagen und Reaktorsicherheit methodology to propagate the uncertainty in thermal-hydraulic macroscopic parameters and, then, shows the results when the methodology is applied to a PWR reactor.

• Nuclear Engineering and Technology
• /
• 제55권6호
• /
• pp.2230-2245
• /
• 2023

KANT (KAIST Advanced Nuclear Tachygraphy) is a PWR core simulator recently developed at Korea Advance Institute of Science and Technology, which solves three-dimensional steady-state and transient multigroup neutron diffusion equations under Cartesian geometries alongside the incorporation of thermal-hydraulics feedback effect for multi-physics calculation. It utilizes the standard Nodal Expansion Method (NEM) accelerated with various Coarse Mesh Finite Difference (CMFD) methods for neutronics calculation. For thermal-hydraulics (TH) calculation, a single-phase flow model and a one-dimensional cylindrical fuel rod heat conduction model are employed. The time-dependent neutronics and TH calculations are numerically solved through an implicit Euler scheme, where a detailed coupling strategy is presented in this paper alongside a description of nodal equivalence, macroscopic depletion, and pin power reconstruction. For validation of the steady, transient, and depletion calculation with pin power reconstruction capacity of KANT, solutions for various benchmark problems are presented. The IAEA 3-D PWR and 4-group KOEBERG problems were considered for the steady-state reactor benchmark problem. For transient calculations, LMW (Lagenbuch, Maurer and Werner) LWR and NEACRP 3-D PWR benchmarks were solved, where the latter problem includes thermal-hydraulics feedback. For macroscopic depletion with pin power reconstruction, a small PWR problem modified with KAIST benchmark model was solved. For validation of the multi-physics analysis capability of KANT concerning large-sized PWRs, the BEAVRS Cycle1 benchmark has been considered. It was found that KANT solutions are accurate and consistent compared to other published works.

• 한국원자력학회:학술대회논문집
• /
• 한국원자력학회 1996년도 춘계학술발표회논문집(1)
• /
• pp.131-136
• /
• 1996

An analytic nodal expansion method has been derived for the multigroup neutron diffusion equation in 2-D cylindrical(R-Z) coordinate. In this method we used the second order Legendre polynomials for source, and transverse leakage, and then the diffusion eqaution was solved analytically. This formalism has been applied to 2-D LWR model. $textsc{k}$$_{eff}$, power distribution, and computing time have been compared with those of ADEP code(finite difference method). The benchmark showed that the analytic nodal expansion method in R-Z coordinate has good accuracy and quite faster than the finite difference method. This is another merit of using R-Z coordinate in that the transverse integration over surfaces is better than the linear integration over length. This makes the discontinuity factor useless.s.