A method and system for calculating fuel consumption of a fuel assembly that moves with a control rod
By using an adaptive meshing method to dynamically divide and merge axial meshes, the matching problem in the calculation of fuel consumption of continuously moving control rod fuel assemblies is solved, improving the calculation accuracy and efficiency, and making it suitable for fuel consumption calculation in nuclear reactors.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- SHANGHAI NUCLEAR ENGINEERING RESEARCH & DESIGN INSTITUTE CO LTD
- Filing Date
- 2023-10-24
- Publication Date
- 2026-06-30
AI Technical Summary
In nuclear reactors, the continuous movement of control rods makes it difficult to match the fuel assembly burnup calculation grid with the neutronics calculation grid, resulting in low computational efficiency and the easy generation of control rod tooth effect. Traditional equidistant axial meshing methods cannot effectively handle the cumulative distribution effect of fuel assemblies.
An adaptive meshing method is adopted to dynamically divide and merge axial meshes based on the positional changes of the control rods and their connected fuel assemblies, calculate the fuel consumption weight, ensure the matching of the fuel consumption calculation mesh with the neutronics calculation mesh, and consider the cumulative distribution effect of the fuel assemblies.
It improves the accuracy and efficiency of fuel consumption calculation, and can more accurately reflect the fuel consumption changes of fuel components without reducing the computational grid, thus solving the computational difficulties caused by control rod movement.
Smart Images

Figure CN117436150B_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the field of nuclear reactor physics numerical calculation technology, and in particular relates to a method and system for calculating the burnup of a fuel assembly that moves with the control rods. Background Technology
[0002] In most nuclear reactors, the fuel assemblies are fixed in position within the core, and control rods move up and down within the fuel assemblies via drive mechanisms to achieve safe control of the reactor. However, in some research reactors, to achieve reactivity control in a compact core, a design that follows the fuel assembly is used. In this design, the fuel assembly and control rod are axially connected, and the up-and-down movement of the control rod also moves the fuel assembly. Unlike step-by-step control rods, these control rods use a continuously moving mechanism, and the positions of the control rod and the following fuel assembly are not fixed. During reactor operation, fuel is continuously consumed, a phenomenon known as burnup. Unlike parameters such as power and neutron flux, which only reflect the current state of the reactor, burnup is a cumulative effect. Fuel assembly burnup is crucial for the safe operation of the reactor and is also an important target parameter in numerical simulation calculations.
[0003] The inventors discovered that, due to computational efficiency requirements, the number of axial layers in reactor neutronics calculations is limited. Traditional methods for calculating the burnup of movable fuel assemblies use equally spaced axial grids. For step-type control rods, the grid spacing equals the step size of the control rod, and the burnup calculation grid and the neutronics calculation grid can be easily matched as the fuel assembly moves with the control rod. However, for continuously moving control rods, a very small calculation grid is required; otherwise, the control rod tooth effect will occur. The smaller the calculation grid, the lower the computational efficiency. Adaptive grid methods are an effective means of handling the control rod tooth effect. During reactor operation, the control rod positions constantly change to compensate for reactivity. When the fuel assembly moves with the control rod, the axial position of the interface between different materials also changes, leading to inconsistent axial layering. Burnup has a cumulative distribution effect, and the continuous change in axial layering makes it difficult to establish an effective mapping between the fuel assembly burnup calculation grid and the neutronics calculation grid. Summary of the Invention
[0004] To address the aforementioned problems, this invention proposes a method and system for calculating the fuel consumption of a fuel assembly that moves with the control rod. This invention ensures the matching of the fuel consumption calculation grid and the neutronics calculation grid, and can consider the cumulative distribution effect of fuel consumption of the fuel assembly moving with the control rod, thereby improving the accuracy of fuel consumption calculation.
[0005] To achieve the above objectives, the present invention is implemented through the following technical solution:
[0006] In a first aspect, the present invention provides a method for calculating the fuel consumption of a fuel assembly that moves following a control rod, comprising:
[0007] A calculation model is established based on the relevant parameters of the reactor core;
[0008] The computational model is meshed to obtain an initial mesh;
[0009] Based on the initial position of the control rod and its connected fuel assembly, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh before the control rod position is changed. Based on the mesh before the control rod position is changed, the fuel consumption value of each mesh in the mesh before the control rod position is changed is calculated according to the given fuel consumption point.
[0010] Based on the position of the control rod and its connected fuel assembly when the control rod position changes, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh after the control rod position changes. The mesh after the control rod position changes is offset according to the amount of change in the axial position of the control rod to obtain the offset mesh.
[0011] Calculate the fuel consumption weight of each grid in the grid before the control rod position change to each grid in the grid after the offset; based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position change, calculate the fuel consumption value of each grid in the grid after the offset.
[0012] Furthermore, the relevant parameters of the reactor core include geometric dimensions and material properties.
[0013] Furthermore, the initial mesh remains unchanged throughout the computation.
[0014] Furthermore, the formulas for calculating the power and fuel consumption values for each grid are as follows:
[0015]
[0016] Where i is the grid number; The fuel consumption value of the previous fuel consumption step is given by k, which is a constant representing the fuel consumption step. This represents the fuel consumption value; Δt k ρ is the time step; i This refers to the initial uranium loading. Power.
[0017] Furthermore, based on the difference between the position of the control rod after the change and the position before the change, the change in the axial position of the control rod is obtained; the change in the axial position of the control rod is subtracted from the upper and lower boundaries of each axial segment of the grid after the change in the control rod position to obtain the offset grid after the change in the control rod position and axial offset.
[0018] Furthermore, based on the percentage of volume occupied by each grid in the grid before the control rod position change to each grid in the grid after the offset, the fuel consumption weight of each grid in the grid before the control rod position change to each grid in the grid after the offset is calculated.
[0019] Furthermore, the mesh of the offset mesh corresponds one-to-one with the mesh of the mesh after the control rod position is changed.
[0020] Secondly, the present invention also provides a fuel consumption calculation system for a fuel assembly that moves following a control rod, comprising:
[0021] The model building module is configured to: build a computational model based on the relevant parameters of the reactor core;
[0022] The initial mesh generation module is configured to: generate an initial mesh from the computational model;
[0023] The first calculation module is configured to: determine the interface of the axial material segment based on the initial position of the control rod and its connected fuel assembly; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh before the control rod position changes; and calculate the fuel consumption value of each mesh in the mesh before the control rod position changes according to the given fuel consumption point based on the mesh before the control rod position changes.
[0024] The offset module is configured to: determine the interface of the axial material segment based on the position of the control rod and its connected fuel assembly when the control rod position changes; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh after the control rod position change; and offset the mesh after the control rod position change according to the amount of change in the axial position of the control rod to obtain the offset mesh.
[0025] The second calculation module is configured to: calculate the fuel consumption weight of each grid in the grid before the control rod position is changed, and the fuel consumption value of each grid in the grid after the offset, based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position is changed.
[0026] Thirdly, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the fuel consumption calculation method for moving the fuel assembly following the control rod as described in the first aspect.
[0027] Fourthly, the present invention also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the program to implement the steps of the fuel consumption calculation method for following the control rod to move the fuel assembly as described in the first aspect.
[0028] Compared with the prior art, the beneficial effects of the present invention are as follows:
[0029] In this invention, firstly, based on the initial position of the control rod and its connected fuel assembly, the interface of the axial material segment is determined. An axial mesh is then divided at the interface, and the divided mesh is merged with the initial mesh to obtain the mesh before the control rod position change. Based on the mesh before the control rod position change, a fuel consumption value is calculated according to a given fuel consumption point to obtain the fuel consumption value of each mesh in the mesh before the control rod position change. Then, based on the position of the control rod and its connected fuel assembly when the control rod position changes, the interface of the axial material segment is determined. An axial mesh is then divided at the interface, and the divided mesh is merged with the initial mesh to obtain the mesh after the control rod position change. The control rod is then adjusted according to the amount of axial position change. After the control rod position is changed, the grid is offset to obtain the offset grid. Finally, the fuel consumption weight of each grid in the grid before the control rod position change is calculated, and the fuel consumption value of each grid in the offset grid is calculated based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position change. The grid generation process fully considers the change of the axial position of the control rod and the relationship between the grids before and after the axial position change of the control rod. It can ensure the matching of the fuel consumption calculation grid and the neutronics calculation grid, and can take into account the cumulative distribution effect of fuel assembly fuel consumption as the control rod moves. It improves the calculation accuracy of fuel consumption without reducing the calculation grid and ensures the calculation efficiency. Attached Figure Description
[0030] The accompanying drawings, which form part of this embodiment, are used to provide a further understanding of this embodiment. The illustrative embodiments and their descriptions are used to explain this embodiment and do not constitute an improper limitation of this embodiment.
[0031] Figure 1 This is a schematic diagram of the axial grid relationship in Embodiment 1 of the present invention;
[0032] Figure 2 This is a schematic diagram of the M0 grid in Embodiment 2 of the present invention;
[0033] Figure 3 This is a schematic diagram of the M1 mesh in Embodiment 2 of the present invention;
[0034] Figure 4 This is a schematic diagram of the M2 mesh in Embodiment 2 of the present invention;
[0035] Figure 5 This is a schematic diagram of the mesh relationship between M1, M2 and M3 in Embodiment 2 of the present invention. Detailed Implementation
[0036] The present invention will be further described below with reference to the accompanying drawings and embodiments.
[0037] It should be noted that the following detailed descriptions are exemplary and intended to provide further explanation of this application. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application pertains.
[0038] The serrated effect refers to the phenomenon where, when the end of the control rod is inserted into a homogeneous layer along the axial direction, the strong absorption effect of the control rod on neutrons increases the axial equivalent homogenization error, resulting in a non-physical "serrated" appearance on the differential value curve of the control rod.
[0039] Example 1:
[0040] For continuously moving control rods, a very small computational grid is required; otherwise, the control rod tooth effect will occur. However, the smaller the computational grid, the lower the computational efficiency. Adaptive grid methods are an effective means of handling the control rod tooth effect. During reactor operation, the control rod positions continuously change to compensate for reactivity. As the fuel assembly moves with the control rods, the axial position of the interfaces between different materials also changes, leading to inconsistent axial stratification. Burnup exhibits a cumulative distribution effect, and the continuous change in axial stratification makes it difficult to establish an effective mapping between the fuel assembly burnup computational grid and the neutronics computational grid. Therefore, designing a burnup calculation method suitable for fuel assemblies that move with the control rods is essential.
[0041] To address the above problems, this embodiment provides a method for calculating the fuel consumption of a fuel assembly that moves with a control rod, including:
[0042] A calculation model is established based on the relevant parameters of the reactor core;
[0043] The computational model is meshed to obtain an initial mesh;
[0044] Based on the initial position of the control rod and its connected fuel assembly, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh before the control rod position is changed. Based on the mesh before the control rod position is changed, the fuel consumption value of each mesh in the mesh before the control rod position is changed is calculated according to the given fuel consumption point.
[0045] Based on the position of the control rod and its connected fuel assembly when the control rod position changes, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh after the control rod position changes. The mesh after the control rod position changes is offset according to the amount of change in the axial position of the control rod to obtain the offset mesh.
[0046] Calculate the fuel consumption weight of each grid cell in the grid before the control rod position is changed, and the fuel consumption value of each grid cell in the grid after the offset. Based on the fuel consumption weight and the fuel consumption value of each grid cell in the grid before the control rod position is changed, calculate the fuel consumption value of each grid cell in the grid after the offset.
[0047] In this embodiment, the mesh generation process fully considers the changes in the axial position of the control rod and the relationship between the meshes before and after the axial position of the control rod. This ensures the matching of the burnup calculation mesh and the neutronics calculation mesh, and takes into account the cumulative distribution effect of fuel assembly burnup as the control rod moves. This improves the accuracy of burnup calculation without reducing the computational mesh size, thus ensuring computational efficiency. In this embodiment, physical parameters with cumulative distribution effects, including but not limited to burnup, cumulative neutron flux, and nucleon density, can all be calculated using the method described in this embodiment. The specific steps of the method in this embodiment are as follows:
[0048] S1. Based on the obtained geometry, material, and other parameters of the target core, complete the core modeling to obtain the computational model. Perform axial meshing on the model according to the specified grid, and use the resulting mesh as the initial mesh M0. The initial mesh M0 remains unchanged throughout the entire calculation process.
[0049] Understandably, the computational model can be understood as a finite element model or other simulation model used for computation; axial meshing according to a specified mesh can be understood as meshing according to a preset mesh type, preset size and / or preset mesh ratio, etc.
[0050] S2. Based on the initial position of the control rod and its connected fuel assembly, determine the interface of the axial material segment, divide the axial grid at the interface position, and merge the divided grid with the initial grid M0 in step S1 as the grid M1 before the control rod position change, and record the control rod position Z1 at this time.
[0051] S3. Using grid M1 before the control rod position change as the neutronics and fuel consumption calculation grid, perform neutronics and fuel consumption calculations according to the given fuel consumption point to obtain the power and fuel consumption values for each grid:
[0052]
[0053] Where i is the grid number; The fuel consumption value (MWd / tU) of the previous fuel consumption step is given, where k is a constant representing the fuel consumption step. This represents the fuel consumption value; Δt k ρ is the time step (s); i Initial uranium loading (kg); Power (MW).
[0054] S4. Repeat step S3 until the position of the control rod changes.
[0055] S5. Based on the position of the control rod and its connected fuel assembly at this time, determine the interface of the axial material segment, divide the axial grid at the interface position, and merge the grid with the initial grid M0 in step S1 as the grid M2 after the control rod position is changed, and record the control rod position Z2 at this time.
[0056] S6. Calculate the axial position change of the control rod ΔZ = Z2 - Z1. Subtract ΔZ from the upper boundary B1 and lower boundary B2 of each axial segment of the grid M2 after the control rod position change in step S5. This gives the position of each axial segment of M2 before the control rod movement. Use this grid as the offset grid M3 after the control rod position change and axial offset by -ΔZ. The offset grid M3 and the grid M1 before the control rod position change may not correspond one-to-one; there may be overlaps or inclusions. Figure 1 As shown, after the offset, grid i in grid M3 intersects with grid j+1 and grid j in grid M1 before the control rod position is changed, and grid i+2 in grid M3 after the offset is completely contained by grid j+2 in grid M1 before the control rod position is changed.
[0057] S7. Based on the crossover and misalignment between grid M1 before the control rod position change and grid M3 after the offset, and according to the volume percentage occupied by each grid in grid M1 before the control rod position change in grid M3 after the offset, calculate the fuel consumption weight P of each grid j in grid M1 before the control rod position change to each grid i in grid M3 after the offset. ij .
[0058] S8. Calculate the fuel consumption of each grid in the offset grid M3. Since the grids of grid M3 after the offset correspond one-to-one with the grids of grid M2 after the control rod position change, the fuel consumption value of each grid in grid M2 after the control rod position change can be obtained. The same approach can be used to calculate physical parameters with cumulative distribution effects, such as cumulative neutron flux and nucleon density.
[0059] S9. Take the current axial grid as the new M1 and the control bar position as the new Z1.
[0060] S10. Repeat steps S4 to S9 until all fuel consumption points have been calculated.
[0061] Example 2:
[0062] This embodiment provides a method for calculating the fuel consumption of a fuel assembly that moves following a control rod. It is a further explanation and supplement to the method in Embodiment 1, specifically as follows:
[0063] S1. Based on the obtained geometry, material, and other parameters of the target reactor core, complete the core modeling. The problem model is then meshed axially according to a specified grid, and this mesh is used as the initial grid M0. The initial grid M0 remains unchanged throughout the calculation process. Figure 2 As shown, the initial mesh M0 is generated without considering the influence of the control rod positions, and the initial mesh M0 is indicated by dashed lines. The illustration only shows the mesh for the active region; the mesh for the inactive region is generated in the same way.
[0064] S2. Based on the initial positions of the control rod and the fuel assembly connected to the control rod, determine the interface of the axial material segment. Divide the axial mesh at the interface position and merge this mesh with the initial mesh M0 from step S1, using it as the mesh M1 before the control rod position change. Record the control rod position Z1 at this time. Figure 3 As shown, the grid before the control rod position is changed is given by dashed lines, and the Z1 position is also marked in the figure.
[0065] S3. Using the grid M1 before the control rod position change as the grid for neutronics and fuel consumption calculations, perform neutronics and fuel consumption calculations according to the given fuel consumption point. Figure 3 The fuel assembly, which moves along with the control rod, is divided into 6 segments before its axial position is changed by the control rod. The fuel consumption value of each segment can be calculated by the following formula:
[0066]
[0067] Where i is the grid number, The fuel consumption value (MWd / tU) of the previous fuel consumption step, Δt k Let ρ be the time step (s). i Initial uranium loading (kg) Power (MW).
[0068] S4. Repeat step S3 until the position of the control rod changes.
[0069] S5. Based on the current positions of the control rod and the fuel assembly connected to the control rod, determine the interface of the axial material segment. Divide the axial mesh at the interface position and merge this mesh with the initial mesh M0 in step S1 to form the mesh M2 after the control rod position change. Record the current control rod position Z2. Figure 4 As shown, after the position of the control rod is changed, the grid M2 is given by the dashed line, and the position of Z2 is also marked in the figure.
[0070] S6. Calculate the axial position change of the control rod ΔZ = Z2 - Z1. Subtract ΔZ from the upper and lower boundaries of each axial segment of grid M2 after the control rod position change in step S5 to obtain the position of each axial segment of M2 before the control rod moves. Use this grid as the offset grid M3. Figure 5 As shown, the fuel assembly connected to the control rod in offset grid M3 is divided into 7 segments axially. M3 and M1 cannot be matched one-to-one, and there are overlapping grids. Therefore, the grids of M1 and M3 are merged, dividing the fuel assembly into 12 segments axially. The length of each segment is s. i exist Figure 4 The bid was successful.
[0071] S7. Based on the crossover and misalignment of M1 and M3, use the volume occupied by each grid in M1 within each grid in M3 as a weight to calculate the fuel consumption of each grid in M3.
[0072]
[0073]
[0074]
[0075]
[0076]
[0077]
[0078]
[0079] S8. Since the grids of M3 correspond one-to-one with those of M2, the fuel consumption value of each grid in M2 can be obtained.
[0080] S9. Take the current axial grid as the new M1 and the control bar position as the new Z1.
[0081] Repeat steps S4 to S9 until all flammability points have been calculated.
[0082] Example 3:
[0083] This embodiment provides a fuel consumption calculation system for a fuel assembly that moves following a control rod, including:
[0084] The model building module is configured to: build a computational model based on the relevant parameters of the reactor core;
[0085] The initial mesh generation module is configured to: generate an initial mesh from the computational model;
[0086] The first calculation module is configured to: determine the interface of the axial material segment based on the initial position of the control rod and its connected fuel assembly; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh before the control rod position changes; and calculate the fuel consumption value of each mesh in the mesh before the control rod position changes according to the given fuel consumption point based on the mesh before the control rod position changes.
[0087] The offset module is configured to: determine the interface of the axial material segment based on the position of the control rod and its connected fuel assembly when the control rod position changes; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh after the control rod position change; and offset the mesh after the control rod position change according to the amount of change in the axial position of the control rod to obtain the offset mesh.
[0088] The second calculation module is configured to: calculate the fuel consumption weight of each grid in the grid before the control rod position is changed, and the fuel consumption value of each grid in the grid after the offset, based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position is changed.
[0089] The operating method of the system is the same as the fuel consumption calculation method of the fuel assembly moving with the control rod in Embodiment 1, and will not be repeated here.
[0090] Example 4:
[0091] This embodiment provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the fuel consumption calculation method for moving a fuel assembly following a control rod as described in Embodiment 1.
[0092] Example 5:
[0093] This embodiment provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, it implements the steps of the fuel consumption calculation method for moving a fuel assembly following a control rod as described in Embodiment 1.
[0094] The above description is merely a preferred embodiment of this practice and is not intended to limit the scope of this practice. Various modifications and variations can be made to this practice by those skilled in the art. Any modifications, equivalent substitutions, or improvements made within the spirit and principles of this practice should be included within the protection scope of this practice.
Claims
1. A method for calculating the fuel consumption of a fuel assembly that moves with a control rod, characterized in that, include: A calculation model is established based on the relevant parameters of the reactor core; The computational model is meshed to obtain an initial mesh; Based on the initial position of the control rod and its connected fuel assembly, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh before the control rod position is changed. Based on the mesh before the control rod position is changed, the fuel consumption value of each mesh in the mesh before the control rod position is changed is calculated according to the given fuel consumption point. Based on the position of the control rod and its connected fuel assembly when the control rod position changes, the interface of the axial material segment is determined. An axial mesh is divided at the interface position, and the divided mesh is merged with the initial mesh to obtain the mesh after the control rod position changes. The mesh after the control rod position changes is offset according to the amount of change in the axial position of the control rod to obtain the offset mesh. Calculate the fuel consumption weight for each grid cell in the grid before the control rod position is changed, and the weight for each grid cell in the grid after the offset. Based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position is changed, the fuel consumption value of each grid in the grid after the offset is calculated.
2. The method for calculating fuel consumption of a fuel assembly following a control rod as described in claim 1, characterized in that, The relevant parameters of the reactor core include geometric dimensions and material properties.
3. The method for calculating fuel consumption of a fuel assembly following a control rod as described in claim 1, characterized in that, The initial mesh remains unchanged throughout the calculation process.
4. The method for calculating fuel consumption of a fuel assembly following a control rod as described in claim 1, characterized in that, The formulas for calculating the power and fuel consumption values for each grid are as follows: Where i is the grid number; The fuel consumption value of the previous fuel consumption step is given by k, which is a constant representing the fuel consumption step. This represents the fuel consumption value; Δt k ρ is the time step; i This represents the initial uranium loading. Power.
5. The method for calculating fuel consumption of a fuel assembly moving with a control rod as described in claim 1, characterized in that, The change in the axial position of the control rod is obtained by comparing its position after the change with its original position. The change in the axial position of the control rod is then subtracted from the upper and lower boundaries of each axial segment of the grid after the change in the control rod position to obtain the offset grid after the change in the control rod position and axial displacement.
6. The method for calculating fuel consumption of a fuel assembly moving with a control rod as described in claim 1, characterized in that, Calculate the fuel consumption weight of each grid in the grid before the control rod position change to each grid in the grid after the offset, based on the percentage of volume occupied by each grid in the grid before the control rod position change.
7. The method for calculating fuel consumption of a fuel assembly following a control rod as described in claim 1, characterized in that, The mesh of the offset grid corresponds one-to-one with the mesh of the grid after the control rod position is changed.
8. A fuel consumption calculation system for a fuel assembly that moves following a control rod, characterized in that, include: The model building module is configured to: build a computational model based on the relevant parameters of the reactor core; The initial mesh generation module is configured to: generate an initial mesh from the computational model; The first calculation module is configured to: determine the interface of the axial material segment based on the initial position of the control rod and its connected fuel assembly; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh before the control rod position changes; and calculate the fuel consumption value of each mesh in the mesh before the control rod position changes according to the given fuel consumption point based on the mesh before the control rod position changes. The offset module is configured to: determine the interface of the axial material segment based on the position of the control rod and its connected fuel assembly when the control rod position changes; divide the axial mesh at the interface position; merge the divided mesh with the initial mesh to obtain the mesh after the control rod position change; and offset the mesh after the control rod position change according to the amount of change in the axial position of the control rod to obtain the offset mesh. The second calculation module is configured to: calculate the fuel consumption weight of each grid in the grid before the control rod position is changed, and the weight of each grid in the grid after the offset. Based on the fuel consumption weight and the fuel consumption value of each grid in the grid before the control rod position is changed, the fuel consumption value of each grid in the grid after the offset is calculated.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by the processor, the program implements the steps of the fuel consumption calculation method for moving the fuel assembly following the control rod as described in any one of claims 1-7.
10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the program, it implements the steps of the fuel consumption calculation method for moving the fuel assembly following the control rod as described in any one of claims 1-7.