• Nuclear Engineering and Technology
• /
• 제50권3호
• /
• pp.327-339
• /
• 2018

This article presents a new global-local hybrid coarse-mesh finite difference (HCMFD) method for efficient parallel calculation of pin-by-pin heterogeneous core analysis. In the HCMFD method, the one-node coarse-mesh finite difference (CMFD) scheme is combined with a nodal expansion method (NEM)-based two-node CMFD method in a nonlinear way. In the global-local HCMFD algorithm, the global problem is a coarse-mesh eigenvalue problem, whereas the local problems are fixed source problems with boundary conditions of incoming partial current, and they can be solved in parallel. The global problem is formulated by one-node CMFD, in which two correction factors on an interface are introduced to preserve both the surface-average flux and the net current. Meanwhile, for accurate and efficient pin-wise core analysis, the local problem is solved by the conventional NEM-based two-node CMFD method. We investigated the numerical characteristics of the HCMFD method for a few benchmark problems and compared them with the conventional two-node NEM-based CMFD algorithm. In this study, the HCMFD algorithm was also parallelized with the OpenMP parallel interface, and its numerical performances were evaluated for several benchmarks.

• Nuclear Engineering and Technology
• /
• 제46권6호
• /
• pp.783-796
• /
• 2014

The coarse mesh finite difference (CMFD) method is applied to the discontinuous finite element method based discrete ordinate calculation for source convergence acceleration. The three-dimensional (3-D) DFEM-Sn code FEDONA is developed for general geometry applications as a framework for the CMFD implementation. Detailed methods for applying the CMFD acceleration are established, such as the method to acquire the coarse mesh flux and current by combining unstructured tetrahedron elements to rectangular coarse mesh geometry, and the alternating calculation method to exchange the updated flux information between the CMFD and DFEM-Sn. The partial current based CMFD (p-CMFD) is also implemented for comparison of the acceleration performance. The modified p-CMFD method is proposed to correct the weakness of the original p-CMFD formulation. The performance of CMFD acceleration is examined first for simple two-dimensional multigroup problems to investigate the effect of the problem and coarse mesh sizes. It is shown that smaller coarse meshes are more effective in the CMFD acceleration and the modified p-CMFD has similar effectiveness as the standard CMFD. The effectiveness of CMFD acceleration is then assessed for three-dimensional benchmark problems such as the IAEA (International Atomic Energy Agency) and C5G7MOX problems. It is demonstrated that a sufficiently converged solution is obtained within 7 outer iterations which would require 175 iterations with the normal DFEM-Sn calculations for the IAEA problem. It is claimed that the CMFD accelerated DFEM-Sn method can be effectively used in the practical eigenvalue calculations involving general geometries.

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

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

• Journal of Electrical Engineering and Technology
• /
• 제6권1호
• /
• pp.128-132
• /
• 2011

In this paper, two automatic mesh generation algorithms are presented. The methods seek to optimize mesh density with regard to geometries exhibiting both fine and coarse physical structures. When generating meshes, the algorithms attempt to satisfy the conditions on the maximum mesh spacing and the maximum grading ratio simultaneously. Both algorithms successfully produce non-uniform meshes that satisfy the requirements for finite-difference time-domain simulations of microwave components. Additionally, an algorithm successfully generates a minimum number of grid points while maintaining the simulation accuracy.

• Nuclear Engineering and Technology
• /
• 제52권10호
• /
• pp.2162-2172
• /
• 2020

Due to its computational efficiency and geometrical flexibility, the Method of Characteristics (MOC) has been widely used for light water reactor lattice physics analysis. Usually acceleration methods are necessary for MOC to achieve acceptable convergence on practical reactor physics problems. Among them, Coarse Mesh Finite Difference (CMFD) is very popular and can drastically reduce the number of transport iterations. In OpenMOC, CMFD acceleration was implemented but had the limitation of supporting only a uniform CMFD mesh, which would often lead to splitting MOC source regions, thus creating an unnecessary increase in computation and memory use. In this study, CMFD acceleration with a non-uniform Cartesian mesh is implemented into OpenMOC. We also propose a quadratic fit based CMFD prolongation method in the axial direction to further improve the acceleration when multiple MOC source regions are contained in one CMFD coarse mesh. Numerical results are presented to demonstrate the improvement of the CMFD acceleration capability in OpenMOC in terms of both efficiency and stability.

• Nuclear Engineering and Technology
• /
• 제56권3호
• /
• pp.772-784
• /
• 2024

The hexagonal nodal code RENUS has been enhanced to handle irregularly deformed hexagonal assemblies. The underlying RENUS methods involving triangle-based polynomial expansion nodal (T-PEN) and corner point balance (CPB) were extended in a way to use line and surface integrals of polynomials in a deformed hexagonal geometry. The nodal calculation is accelerated by the coarse mesh finite difference (CMFD) formulation extended to unstructured geometry. The accuracy of the unstructured nodal solution was evaluated for a group of 2D SFR core problems in which the assembly corner points are arbitrarily displaced. The RENUS results for the change in nuclear characteristics resulting from fuel deformation were compared with those of the reference McCARD Monte Carlo code. It turned out that the two solutions agree within 18 pcm in reactivity change and 0.46% in assembly power distribution change. These results demonstrate that the proposed unstructured nodal method can accurately model heterogeneous thermal expansion in hexagonal fueled cores.

• Nuclear Engineering and Technology
• /
• 제54권8호
• /
• pp.3059-3072
• /
• 2022

A Jacobian-Free Newton Krylov Two-Nodal Coarse Mesh Finite Difference algorithm based on Nodal Expansion Method (NEM_TNCMFD_JFNK) is successfully developed and proposed to solve the three-dimensional (3D) and multi-group reactor physics models. In the NEM_TNCMFD_JFNK method, the efficient JFNK method with the Modified Incomplete LU (MILU) preconditioner is integrated and applied into the discrete systems of the NEM-based two-node CMFD method by constructing the residual functions of only the nodal average fluxes and the eigenvalue. All the nonlinear corrective nodal coupling coefficients are updated on the basis of two-nodal NEM formulation including the discontinuity factor in every few newton steps. All the expansion coefficients and interface currents of the two-node NEM need not be chosen as the solution variables to evaluate the residual functions of the NEM_TNCMFD_JFNK method, therefore, the NEM_TNCMFD_JFNK method can greatly reduce the number of solution variables and the computational cost compared with the JFNK based on the conventional NEM. Finally the NEM_TNCMFD_JFNK code is developed and then analyzed by simulating the representative PWR MOX/UO₂ core benchmark, the popular NEACRP 3D core benchmark and the complicated full-core pin-by-pin homogenous core model. Numerical solutions show that the proposed NEM_TNCMFD_JFNK method with the MILU preconditioner has the good numerical accuracy and can obtain higher computational efficiency than the NEM-based two-node CMFD algorithm with the power method in the outer iteration and the Krylov method using the MILU preconditioner in the inner iteration, which indicates the NEM_TNCMFD_JFNK method can serve as a potential and efficient numerical tool for reactor neutron diffusion analysis module in the JFNK-based multiphysics coupling application.

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

• Nuclear Engineering and Technology
• /
• 제34권3호
• /
• pp.250-258
• /
• 2002

This Paper is to apply the well-established coarse mesh finite difference(CMFD) method to the method of characteristics(MOC) transport calculation as an acceleration scheme. The CMFD problem is first formulated at the pin-cell level with the multi-group structure To solve the cell- based multi-group CMFD problem efficiently, a two-group CMFD formulation is also derived from the multi-group CMFD formulation. The performance of the CMFD acceleration is examined for three test problems with different sizes including a realistic quarter core PWR problem. The CMFD formulation provides a significant reduction in the number of ray tracings and thus only about 9 ray tracing iterations are enough for the realistic problem. In computing time, the CMFD accelerated case is about two or three times faster than the coarse-mesh rebalancing(CMR) accelerated case.