Study on the accelerated method for solving the CRAM burnup equations based on Gauss-Seidel method

To evaluate computational accuracy, the results of the SGE and GS methods were compared, thereby demonstrating the precision of the GS method. Furthermore, with the Transmutation Trajectory Analysis (TTA) as the reference solution, detailed analyses of the computational results of the Partial Fraction Decomposition (PFD) and Incomplete Partial Fractions (IPF) formulations based on the GS method were conducted. Finally, a comparative analysis of the computational efficiency of the SGE and GS methods was performed.

[Results]: The results show that the numerical precision of GS method is good for solving IPF and PFD burnup equations under different scale burnup databases, and can meet the accuracy requirements of actual burnup calculation. For the calculation of short-lived nuclides, the calculation accuracy of IPF is better than that of PFD. In terms of efficiency, the GS method significantly surpassed SGE, achieving up to 80.17% acceleration across the three databases.

[Conclusions]: This study recommends the adoption of the GS-accelerated IPF formalism for practical burnup calculations to effectively balance computational accuracy and efficiency.