Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order

The goal of high-fidelity burnup calculations is to accurately describe the nuclide composition of the reactor over its entire lifetime. As the core physics calculation method evolves to high-fidelity, the number of burnup regions will increase thousands of times from the traditional decimeter level...

Full description

Bibliographic Details
Main Author: ZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao2
Format: Article
Language:English
Published: Editorial Board of Atomic Energy Science and Technology 2023-10-01
Series:Yuanzineng kexue jishu
Subjects:
Online Access:https://www.aest.org.cn/CN/10.7538/yzk.2022.youxian.0904
_version_ 1797635149232340992
author ZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao2
author_facet ZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao2
author_sort ZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao2
collection DOAJ
description The goal of high-fidelity burnup calculations is to accurately describe the nuclide composition of the reactor over its entire lifetime. As the core physics calculation method evolves to high-fidelity, the number of burnup regions will increase thousands of times from the traditional decimeter level to the centimeter level, and the cost and complexity of the burnup calculation will increase steeply while the accuracy of the burnup calculation increases. In this paper, the fuel consumption calculation method based on the proper orthogonal decomposition (POD) reduced-order method was studied. By establishing a parametric reduced-order model (ROM), a set of complete orthogonal basis functions was used to construct the burnup matrix in each burnup region, and further based on the offline-online calculation strategy to achieve efficient numerical solution. In the offline phase, a high-fidelity burnup calculation was performed to obtain the snapshot matrices of reaction rate and neutron flux. Based on singular value decomposition (SVD), a set of basis functions of the burnup matrix and neutron flux were constructed respectively. In the online phase, the neutron flux was obtained from the real-time transport calculation, the parametric ROM and basis functions in the offline stage were used to quickly construct the burnup matrix and then completed the full-life burnup calculation. Based on the benchmark problem of MOX fuel cell released by JAEA, the numerical results verify the correctness and feasibility of this method. For the calculation of kinf, the relative deviation from the full-order model (FOM) is less than 0.5% with only the first-order ROM. For the calculation of nuclear density of heavy nuclides 235U and 239Pu, the first-, fifth- and tenth-order ROMs were used for calculation. Comparing with the FOM, the relative deviations are within 5%, 0.1% and 0.005%, respectively. However, for the calculation of 155Gd and 149Sm, which have larger neutron absorption cross sections, the calculation accuracy of the first-order ROM is poor, the relative deviation of the fifth-order ROM is within 5%, and the relative deviation of the tenth-order ROM reaches within 0.1%. From the calculation results, it can be seen that a ROM of about fifth-order can be used to construct the burnup matrix for 3 fuel regions at every burnup step, which can obtain good calculation accuracy without repeatedly calculating and storing the burnup matrix. Comparing with the traditional burnup calculation method, the acceleration ratio obtained by different reduced-order models is up to 18-23. Compared with the traditional burnup calculation, the acceleration ratio of the 5-order reduced-order model can reach 17-26, which improves the calculation efficiency. This method can be used as a reference to effectively reduce the cost and improve the efficiency of high-fidelity burnup calculations.
first_indexed 2024-03-11T12:18:13Z
format Article
id doaj.art-fc67e28312644c158ff09837df9d14a5
institution Directory Open Access Journal
issn 1000-6931
language English
last_indexed 2024-03-11T12:18:13Z
publishDate 2023-10-01
publisher Editorial Board of Atomic Energy Science and Technology
record_format Article
series Yuanzineng kexue jishu
spelling doaj.art-fc67e28312644c158ff09837df9d14a52023-11-07T07:16:32ZengEditorial Board of Atomic Energy Science and TechnologyYuanzineng kexue jishu1000-69312023-10-0157101938194810.7538/yzk.2022.youxian.0904Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-orderZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao201.College of Mechanical & Power Engineering, China Three Gorges University, Yichang 443002, China;2.College of Science, China Three Gorges University, Yichang 443002, ChinaThe goal of high-fidelity burnup calculations is to accurately describe the nuclide composition of the reactor over its entire lifetime. As the core physics calculation method evolves to high-fidelity, the number of burnup regions will increase thousands of times from the traditional decimeter level to the centimeter level, and the cost and complexity of the burnup calculation will increase steeply while the accuracy of the burnup calculation increases. In this paper, the fuel consumption calculation method based on the proper orthogonal decomposition (POD) reduced-order method was studied. By establishing a parametric reduced-order model (ROM), a set of complete orthogonal basis functions was used to construct the burnup matrix in each burnup region, and further based on the offline-online calculation strategy to achieve efficient numerical solution. In the offline phase, a high-fidelity burnup calculation was performed to obtain the snapshot matrices of reaction rate and neutron flux. Based on singular value decomposition (SVD), a set of basis functions of the burnup matrix and neutron flux were constructed respectively. In the online phase, the neutron flux was obtained from the real-time transport calculation, the parametric ROM and basis functions in the offline stage were used to quickly construct the burnup matrix and then completed the full-life burnup calculation. Based on the benchmark problem of MOX fuel cell released by JAEA, the numerical results verify the correctness and feasibility of this method. For the calculation of kinf, the relative deviation from the full-order model (FOM) is less than 0.5% with only the first-order ROM. For the calculation of nuclear density of heavy nuclides 235U and 239Pu, the first-, fifth- and tenth-order ROMs were used for calculation. Comparing with the FOM, the relative deviations are within 5%, 0.1% and 0.005%, respectively. However, for the calculation of 155Gd and 149Sm, which have larger neutron absorption cross sections, the calculation accuracy of the first-order ROM is poor, the relative deviation of the fifth-order ROM is within 5%, and the relative deviation of the tenth-order ROM reaches within 0.1%. From the calculation results, it can be seen that a ROM of about fifth-order can be used to construct the burnup matrix for 3 fuel regions at every burnup step, which can obtain good calculation accuracy without repeatedly calculating and storing the burnup matrix. Comparing with the traditional burnup calculation method, the acceleration ratio obtained by different reduced-order models is up to 18-23. Compared with the traditional burnup calculation, the acceleration ratio of the 5-order reduced-order model can reach 17-26, which improves the calculation efficiency. This method can be used as a reference to effectively reduce the cost and improve the efficiency of high-fidelity burnup calculations.https://www.aest.org.cn/CN/10.7538/yzk.2022.youxian.0904burnup calculationreduced-order modelproper orthogonal decompositionnumerical validation
spellingShingle ZHANG Binhang1;BI Yanzhao1;GONG Hanyuan1;ZHANG Yonghong1,*;YUAN Xianbao2
Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
Yuanzineng kexue jishu
burnup calculation
reduced-order model
proper orthogonal decomposition
numerical validation
title Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
title_full Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
title_fullStr Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
title_full_unstemmed Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
title_short Study on Burnup Calculation Method Based on Proper Orthogonal Decomposition Reduced-order
title_sort study on burnup calculation method based on proper orthogonal decomposition reduced order
topic burnup calculation
reduced-order model
proper orthogonal decomposition
numerical validation
url https://www.aest.org.cn/CN/10.7538/yzk.2022.youxian.0904
work_keys_str_mv AT zhangbinhang1biyanzhao1gonghanyuan1zhangyonghong1yuanxianbao2 studyonburnupcalculationmethodbasedonproperorthogonaldecompositionreducedorder