A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method
Coarse-mesh finite difference (CMFD) method is a widely used numerical acceleration method. However, the stability of CMFD method is not good for the problems with optically thick regions. In this paper, a stability rule named the “sign preservation rule” in the field of numerical heat transfer is e...
Main Authors: | , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Frontiers Media S.A.
2022-05-01
|
Series: | Frontiers in Energy Research |
Subjects: | |
Online Access: | https://www.frontiersin.org/articles/10.3389/fenrg.2022.836363/full |
_version_ | 1818213626187612160 |
---|---|
author | Zhitao Xu Hongchun Wu Youqi Zheng Qian Zhang |
author_facet | Zhitao Xu Hongchun Wu Youqi Zheng Qian Zhang |
author_sort | Zhitao Xu |
collection | DOAJ |
description | Coarse-mesh finite difference (CMFD) method is a widely used numerical acceleration method. However, the stability of CMFD method is not good for the problems with optically thick regions. In this paper, a stability rule named the “sign preservation rule” in the field of numerical heat transfer is extended to the scheme of CMFD. It is required that the disturbance of neutron current is positively correlated with that of the negative value of flux gradient. A necessary condition for stability of the CMFD method is derived, an adaptive diffusion coefficient equation is proposed to improve the stability of CMFD method, and the corresponding revised CMFD method is called the rCMFD method. With a few modifications of the code, the rCMFD method was implemented in the hexagonal-Z nodal SN (discrete-ordinates) solver in the NECP-SARAX code system. The rCMFD method and other similar acceleration methods were tested by three fast reactor problems which were obtained by modifying the hexagonal pitches of a benchmark problem. The numerical results indicated that the rCMFD method showed better stability than the traditional CMFD method and the artificially diffusive CMFD (adCMFD) method and a better convergence rate than the adCMFD method and the optimally diffusive CMFD (odCMFD) method for these fast reactor problems. |
first_indexed | 2024-12-12T06:07:17Z |
format | Article |
id | doaj.art-af66c44dd62f45bca70c2ba86ae1946d |
institution | Directory Open Access Journal |
issn | 2296-598X |
language | English |
last_indexed | 2024-12-12T06:07:17Z |
publishDate | 2022-05-01 |
publisher | Frontiers Media S.A. |
record_format | Article |
series | Frontiers in Energy Research |
spelling | doaj.art-af66c44dd62f45bca70c2ba86ae1946d2022-12-22T00:35:16ZengFrontiers Media S.A.Frontiers in Energy Research2296-598X2022-05-011010.3389/fenrg.2022.836363836363A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference MethodZhitao Xu0Hongchun Wu1Youqi Zheng2Qian Zhang3School of Electric Power, North China University of Water Resources and Electric Power, Zhengzhou, ChinaSchool of Nuclear Science and Technology, Xi’an Jiaotong University, Xi’an, ChinaSchool of Nuclear Science and Technology, Xi’an Jiaotong University, Xi’an, ChinaSchool of Physics, Zhejiang University, Hangzhou, ChinaCoarse-mesh finite difference (CMFD) method is a widely used numerical acceleration method. However, the stability of CMFD method is not good for the problems with optically thick regions. In this paper, a stability rule named the “sign preservation rule” in the field of numerical heat transfer is extended to the scheme of CMFD. It is required that the disturbance of neutron current is positively correlated with that of the negative value of flux gradient. A necessary condition for stability of the CMFD method is derived, an adaptive diffusion coefficient equation is proposed to improve the stability of CMFD method, and the corresponding revised CMFD method is called the rCMFD method. With a few modifications of the code, the rCMFD method was implemented in the hexagonal-Z nodal SN (discrete-ordinates) solver in the NECP-SARAX code system. The rCMFD method and other similar acceleration methods were tested by three fast reactor problems which were obtained by modifying the hexagonal pitches of a benchmark problem. The numerical results indicated that the rCMFD method showed better stability than the traditional CMFD method and the artificially diffusive CMFD (adCMFD) method and a better convergence rate than the adCMFD method and the optimally diffusive CMFD (odCMFD) method for these fast reactor problems.https://www.frontiersin.org/articles/10.3389/fenrg.2022.836363/fullstable conditionadaptive diffusion coefficientsrCMFDacceleration methodstabilityIFDF |
spellingShingle | Zhitao Xu Hongchun Wu Youqi Zheng Qian Zhang A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method Frontiers in Energy Research stable condition adaptive diffusion coefficients rCMFD acceleration method stability IFDF |
title | A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method |
title_full | A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method |
title_fullStr | A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method |
title_full_unstemmed | A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method |
title_short | A Stable Condition and Adaptive Diffusion Coefficients for the Coarse-Mesh Finite Difference Method |
title_sort | stable condition and adaptive diffusion coefficients for the coarse mesh finite difference method |
topic | stable condition adaptive diffusion coefficients rCMFD acceleration method stability IFDF |
url | https://www.frontiersin.org/articles/10.3389/fenrg.2022.836363/full |
work_keys_str_mv | AT zhitaoxu astableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT hongchunwu astableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT youqizheng astableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT qianzhang astableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT zhitaoxu stableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT hongchunwu stableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT youqizheng stableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod AT qianzhang stableconditionandadaptivediffusioncoefficientsforthecoarsemeshfinitedifferencemethod |