Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method
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 metho...
| Main Author: | |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Elsevier
2022-08-01
|
| Series: | Nuclear Engineering and Technology |
| Subjects: | |
| Online Access: | http://www.sciencedirect.com/science/article/pii/S1738573322000766 |
| _version_ | 1798038146987851776 |
|---|---|
| author | Xiafeng Zhou |
| author_facet | Xiafeng Zhou |
| author_sort | Xiafeng Zhou |
| collection | DOAJ |
| description | 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/UO2 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. |
| first_indexed | 2024-04-11T21:36:15Z |
| format | Article |
| id | doaj.art-ef4e133bddca431d872fffaf40a5bb05 |
| institution | Directory Open Access Journal |
| issn | 1738-5733 |
| language | English |
| last_indexed | 2024-04-11T21:36:15Z |
| publishDate | 2022-08-01 |
| publisher | Elsevier |
| record_format | Article |
| series | Nuclear Engineering and Technology |
| spelling | doaj.art-ef4e133bddca431d872fffaf40a5bb052022-12-22T04:01:44ZengElsevierNuclear Engineering and Technology1738-57332022-08-0154830593072Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion methodXiafeng Zhou0Department of Nuclear Engineering and Technology, School of Energy and Power Engineering, Huazhong University of Science and Technology, Wuhan, 430074, ChinaA 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/UO2 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.http://www.sciencedirect.com/science/article/pii/S1738573322000766Jacobian-free Newton-KrylovTwo-node coarse mesh finite differenceNodal expansion methodModified incomplete LU preconditioner |
| spellingShingle | Xiafeng Zhou Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method Nuclear Engineering and Technology Jacobian-free Newton-Krylov Two-node coarse mesh finite difference Nodal expansion method Modified incomplete LU preconditioner |
| title | Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method |
| title_full | Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method |
| title_fullStr | Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method |
| title_full_unstemmed | Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method |
| title_short | Jacobian-free Newton Krylov two-node coarse mesh finite difference based on nodal expansion method |
| title_sort | jacobian free newton krylov two node coarse mesh finite difference based on nodal expansion method |
| topic | Jacobian-free Newton-Krylov Two-node coarse mesh finite difference Nodal expansion method Modified incomplete LU preconditioner |
| url | http://www.sciencedirect.com/science/article/pii/S1738573322000766 |
| work_keys_str_mv | AT xiafengzhou jacobianfreenewtonkrylovtwonodecoarsemeshfinitedifferencebasedonnodalexpansionmethod |