CN115828643A - Method for simultaneously moving rod search thresholds of multiple control rod groups in consideration of overlapping steps - Google Patents

Method for simultaneously moving rod search thresholds of multiple control rod groups in consideration of overlapping steps Download PDF

Info

Publication number
CN115828643A
CN115828643A CN202310115637.4A CN202310115637A CN115828643A CN 115828643 A CN115828643 A CN 115828643A CN 202310115637 A CN202310115637 A CN 202310115637A CN 115828643 A CN115828643 A CN 115828643A
Authority
CN
China
Prior art keywords
control rod
height
calculation
search
rod group
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310115637.4A
Other languages
Chinese (zh)
Other versions
CN115828643B (en
Inventor
刘宙宇
苏鑫
吴宏春
曹良志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN202310115637.4A priority Critical patent/CN115828643B/en
Publication of CN115828643A publication Critical patent/CN115828643A/en
Application granted granted Critical
Publication of CN115828643B publication Critical patent/CN115828643B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The invention discloses a method for simultaneously operating search criticality by considering a plurality of control rod groups in a cascade manner, which comprises the following steps of firstly, respectively inserting or extracting all the control rod groups participating in the search criticality in a reactor core according to the initial effective value-added factors of the reactor core of the search criticality; then, calculating the whole reactor core to obtain a new effective value-added factor of the reactor core; according to the total moving length of the control rod groups and the change of the effective value-added factors of the reactor core, total height interpolation and rod position height calculation are carried out to obtain the new rod position of each control rod group when the reactor core is critical; adjusting the corresponding control rod group to a new rod position and carrying out the transport calculation of the whole reactor core again; the above process is repeated until the core reaches the critical state. The method can not only carry out the search-critical calculation on a plurality of control rod groups together, reduce the total reactor core transport calculation times in the traditional search-critical calculation process, but also consider the condition that the control rod groups are stacked to carry out the search-critical calculation. The method can be used for the search-critical calculation of the numerical reactor, and improves the calculation speed and efficiency.

Description

Method for simultaneously moving rod search thresholds of multiple control rod groups in consideration of overlapping steps
Technical Field
The invention relates to the technical field of nuclear reactor core design and safety, in particular to a method for simultaneously performing rod search on a plurality of control rod groups in consideration of overlapping steps.
Background
The critical process of reactor core search is an important working condition in the reactor operation process, and the reactor core is required to be in a critical state before the calculation of the reactor core burnup point and the calculation of the reactor transient working condition. Controlling core reactivity by varying the position of the control rod sets is an important means of core criticality.
When a large pressurized water reactor is operated, the requirements on the reactor core are in a critical state. However, as the reactor operating time increases, the burn-up depth of the burnable material increases, so that the core cannot be maintained critical, and therefore, the core criticality is achieved by moving the control rod groups inside the core at different burn-up points. In the transient condition of the nuclear reactor, the reactor core must be in a critical state before the transient condition occurs. For the reactor physical procedure, if the transient operating condition needs to be calculated, the position of the control rod needs to be continuously changed, the calculation is repeatedly carried out, the transient calculation can not be carried out until the reactor core reaches the critical state, and the process is extremely complicated and complex.
The current critical search method mainly used in the numerical reactor process is the conventional linear interpolation method. And critical search is carried out on the control rod groups one by one through the specified control rod group moving rod sequence until the reactor core reaches the critical state when one control rod group reaches a certain position. However, when the number of rod groups participating in the search is too large, the method is inefficient because the reactor is moved at least once for all the cores before reaching the critical value, which consumes a lot of time, and the search is difficult to consider the cascade movement of the control rod groups.
Disclosure of Invention
In order to overcome the problems in the prior art, the present invention provides a method for simultaneously performing a search for a plurality of control rod groups in consideration of overlapping steps, which does not use the conventional method for searching for the control rod groups one by one, and solves the total height of the control rod groups by interpolation using all the control rod groups participating in the search as a whole, thereby obtaining the height of each rod group. Compared with the traditional method, the method greatly reduces the number of times of calculating the whole reactor core of the search criticality under a plurality of control rod groups, reduces the time required by the search criticality calculation, and can realize the process of moving the rod search criticality in a cascading manner.
In order to achieve the purpose, the invention adopts the following technical scheme to implement:
a method for simultaneously operating a plurality of control rod groups in consideration of overlapping steps, comprising the steps of:
step 1: reading geometric information, material information, boundary conditions and control rod group information of a critical core to be simulated;
step 2: obtaining information of each control stick group participating in the search threshold according to the control stick group information obtained in the step 1, wherein the information comprises the sequence of each control stick group participating in the search threshold, the minimum and maximum movable heights and the initial position of each control stick group; then, performing neutron transport calculation on the whole reactor core to obtain the effective reactor core multiplication coefficient in an initial state;
and step 3: respectively processing according to the difference between the effective reactor core multiplication coefficient calculated in the step 2 and the absolute difference value of 1; if the absolute difference value of the effective reactor core multiplication coefficient and 1 is greater than the convergence criterion and the effective reactor core multiplication coefficient is greater than 1, inserting all control rod groups of the reactor core participating in the critical search into the reactor core by moving the control rods; if the absolute difference value of the effective reactor core multiplication coefficient and 1 is greater than the convergence criterion and the effective reactor core multiplication coefficient is less than 1, all control rod groups of the reactor core participating in the critical search are lifted out of the reactor core by moving the control rods; if the absolute difference value between the effective multiplication coefficient of the reactor core and 1 is smaller than the convergence criterion, the reactor core reaches a critical state at the moment, and the calculation is finished;
and 4, step 4: carrying out the transport calculation of the whole reactor core again according to the new rod position of the control rod group to obtain the effective value-added factors of the reactor core under the new rod position of the control rod group; calculating to obtain the difference of the total heights of the control rod groups in two adjacent whole core calculations through a formula (4.1), then calculating to obtain the average differential value of the control rod groups in the two whole core calculations through a formula (4.2), and finally carrying out total height interpolation calculation through a formula (4.3) and a formula (4.4) to obtain the total height of a new control rod group;
Figure SMS_1
(4.1)
Figure SMS_2
(4.2)
Figure SMS_3
(4.3)
Figure SMS_4
(4.4)
wherein, the first and the second end of the pipe are connected with each other,
Figure SMS_5
-the difference between the total height of the control rod groups participating in the search criticality in the jth full core calculation and the total height of the control rod groups participating in the search criticality in the last full core calculation
Figure SMS_6
Figure SMS_7
-the number of control stick groups participating in the search;
Figure SMS_8
-the height of the ith control rod group at the jth full core calculation;
Figure SMS_9
-the height of the ith control rod group at the time of j-1 th full core calculation;
Figure SMS_10
after j th full core calculation, comparing with the last full core calculationThe average differential value of the control rod groups participating in the search;
Figure SMS_11
obtaining the effective increment coefficient of the reactor core after the jth total reactor core calculation;
Figure SMS_12
the effective incremental coefficient of the reactor core obtained after j-1 th calculation of the whole reactor core;
Figure SMS_13
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure SMS_14
-total height of control rod groups participating in the search for criticality at jth full core calculation;
Figure SMS_15
the total height of the control rod groups participating in the criticality search during the 1 st full core calculation;
Figure SMS_16
-height of ith control rod group at 1 st full core calculation;
and 5: if the overlapping movement of the control rod groups is not considered, numbering the N control rod groups according to the searching critical sequence of the control rod groups and the maximum height capable of moving of the control rod groups obtained in the step 2 into 1-N, and supposing that the heights of m control rod groups are the maximum at the moment; calculating according to a formula (4.5) and a formula (4.6) to obtain the height of each control rod group in the next total reactor core calculation;
Figure SMS_17
(4.5)
Figure SMS_18
(4.6)
wherein, the first and the second end of the pipe are connected with each other,
h i+1,j+1 -the height of the (i + 1) th control rod group at the time of the (j + 1) th total core calculation;
Figure SMS_19
-the maximum height that the ith control stick set can move;
if the stacking movement of the control rod groups is considered, assuming that the stacking height of each control rod group during the stacking movement is d, the maximum height of the control rod groups capable of moving is equal;
when the control bar group numbered 1 is moved in steps to the maximum height, it is assumed that it is common at this time
Figure SMS_20
The height of each control rod group is more than 0; the height sum of each control rod group is obtained through the formula (4.7)
Figure SMS_21
Then the total height of the control rod group at the moment is calculated by the formula (4.8);
Figure SMS_22
(4.7)
Figure SMS_23
(4.8)
wherein the content of the first and second substances,
Figure SMS_24
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
Figure SMS_25
1 st controlThe maximum height that the rod making group can move;
Figure SMS_26
-stack height between control rod groups;
Figure SMS_27
-the number of control stick groups participating in the search;
Figure SMS_28
when the control rod group with the number of 1 moves to the maximum height in a stacking mode, the total height of the control rod group participating in searching the critical area is increased;
Figure SMS_29
when the control rod group with the number of 1 moves to the maximum height in the stacking step, the height of the control rod group is more than 0;
when in use
Figure SMS_30
If so, assuming that the height of the n control rod groups is more than 0; the height of the control rod group with the number from n to 1 is shown in a formula (4.9), and the size of n is calculated according to the formula (4.10);
Figure SMS_31
(4.9)
when the temperature is higher than the set temperature
Figure SMS_32
When N = N;
otherwise, the formula (4.10) is satisfied
Figure SMS_33
(4.10)
Calculated according to equation (4.10):
Figure SMS_34
since n is a positive integer and the value range of n in the calculation result of the formula (4.10) is less than 1, only one integer n meets the condition; after n is obtained through calculation, the height of each control rod group is obtained through calculation according to a formula (4.11) and a formula (4.9);
Figure SMS_35
(4.11)
wherein the content of the first and second substances,
Figure SMS_36
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure SMS_37
the jth total core calculation, the th
Figure SMS_38
The height of each control rod set;
Figure SMS_39
-stack height between control rod groups;
Figure SMS_40
-the number of control rod groups having a height greater than 0;
Figure SMS_41
-the number of control stick groups participating in the search;
Figure SMS_42
-the number of control rod groups with number 1, the control rod group height being greater than 0, when the control rod groups move in steps to the maximum height;
Figure SMS_43
when the control rod group with the number of 1 moves to the maximum height in a stacking mode, the total height of the control rod group participating in searching the critical area is increased;
Figure SMS_44
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
when H is present j+1 >H critical And is
Figure SMS_45
Figure SMS_46
When the height of the control rod group reaches the maximum height, the size of m is calculated according to a formula (4.12);
Figure SMS_47
(4.12)
calculated according to equation (4.12):
Figure SMS_48
since m is a positive integer and the value range of m in the calculation result of the formula (4.12) is less than 1, only one integer m meets the condition, and after m is obtained through calculation, the heights of the rest bar groups are obtained through calculation of the formula (4.10), the formula (4.11) and the formula (4.9);
when H is present j+1 >H critical
Figure SMS_49
If the number of the control rod groups with the rod position height smaller than the maximum value is k, calculating to obtain the new rod position of each control rod group according to the following formula (4.13);
Figure SMS_50
(4.13)
calculated according to equation (4.13):
Figure SMS_51
k is a positive integer, and the value range of k in the calculation result of the formula (4.13) is less than 1, so that only one integer k meets the condition; then calculating the size of m by m = N-k, and finally calculating the height of the rest control rod groups by (4.11) and a formula (4.9);
wherein the content of the first and second substances,
Figure SMS_52
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure SMS_53
-the height of the ith control rod group at the jth full core calculation;
Figure SMS_54
-the maximum height that the ith control bar set can move;
Figure SMS_55
-stack height between control rod groups;
Figure SMS_56
-the number of control rod groups having a height greater than 0;
Figure SMS_57
-the number of control stick groups participating in the search;
Figure SMS_58
-the number of control rod groups with number 1, the control rod group height being greater than 0, when the control rod groups move in steps to the maximum height;
Figure SMS_59
-the total height of the control rod group participating in the search for criticality when the control rod group numbered 1 moves in steps to the maximum height;
Figure SMS_60
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
Figure SMS_61
-the number of control rod groups having a rod position height smaller than the maximum value;
step 6: repeating the processes in the step 4 and the step 5 until the absolute difference value between the effective reactor core multiplication coefficient obtained by the calculation of the whole reactor core and 1 is less than the convergence criterion under a certain control rod group height; at the moment, the critical calculation of the reactor core is finished, and the finally obtained height of each control rod group is the critical rod position of the reactor core control rod group.
The convergence criterion of the step 3 and the step 6 considers the tine effect and the related test basis in the control rod moving process, and the convergence criterion is set not to exceed
Figure SMS_62
Compared with the prior art, the invention has the following outstanding advantages:
compared with the traditional method, the method disclosed by the invention has the advantages that the problem of the core search criticality of a plurality of control rod groups is avoided being calculated by adopting a traditional single rod group linear interpolation mode, the total heights of all the control rod groups participating in the search criticality are interpolated, the heights of all the control rod groups are respectively solved, the calculation times of the whole core are reduced, the calculation time of the core search criticality is greatly reduced, and the function of moving the core search criticality in a stacked step is realized.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a schematic view of a 3 × 3 single component.
Detailed Description
The invention is described in further detail below with reference to the drawings and the detailed description.
The present invention has been tested in the high fidelity physical computing program NECP-X. The specific implementation steps are shown in fig. 1. The simulated reactor core is first modeled using the NECP-X program based on the geometric information, material information, boundary conditions and control rod set information for the core in the input card. By processing the information in the input card in the NECP-X program, information of each control stick group participating in the search threshold is obtained, including the sequence of participation of each control stick group in the search threshold, the minimum and maximum heights of the movement possible, and the initial position of each control stick group. A first full core transport calculation is then made based on the initial position of the control rods. And carrying out full interpolation or full extraction treatment on the control rod groups participating in the search criticality according to the calculated result, then carrying out full stack calculation again, obtaining the rod position of a new control rod group through total height interpolation and rod position calculation, and carrying out full core transport calculation again according to the rod position of the new control rod group. And repeating the process until the absolute difference value between the effective value-added coefficient of the reactor core obtained by the whole-reactor calculation under a certain rod position and 1 is smaller than the convergence criterion, and ending the critical calculation, wherein the rod position of the control rod group at the moment is the critical rod position. The following describes the specific computation steps using a 3 × 3 single component problem:
step 1: a 3 x 3 single component as shown in fig. 2, the component consists of 9 cells, with 9 rods in each cell. B, C, D and E are four different control rod groups, A is a guide tube, MOX fuel is adopted as the material of the fuel rods in the fourth grid cell, and UO is adopted as the material of the fuel rods in the other grid cells 2 The fuel and the rest materials except the fuel rods in the grid cells are all water.
Step 2: there are a total of four control rod groups in the core: b, C, D, E; wherein the control rod groups participating in the search are sequentially B, C and D. Setting the minimum height moved by each bar group to be 0 (full interpolation), the maximum height to be 60 steps (1 step =0.714 cm), the initial positions of the four bar groups to be full interpolation, and the convergence criterion to be 30 × 10 -5 . Then, the neutron transport calculation of the whole reactor core is carried out to obtainUntil the initial effective multiplication coefficient of the reactor core is 0.896115;
and step 3: since the effective core multiplication coefficient calculated in step 2 is less than 1 and the absolute difference between the effective core multiplication coefficient and 1 is greater than the convergence criterion, the NECP-X program moves all the control rod groups bc D of which the core participates in the search criticality to the maximum height by the control rod movement by 60 steps;
and 4, step 4: and performing the whole core calculation again according to the new rod position of the B C D control rod group. The new effective value-added factor of the reactor core is 1.026925. Calculating by the formula (4.1) to obtain that the total length of the control rod group in the two adjacent whole core calculations is 180 steps, and then calculating by the formula (4.2) to obtain that the average differential value of the control rod group in the two whole core calculations is 72.67 × 10 -5 Step one, finally, carrying out total height interpolation calculation through a formula (4.3) to obtain the total height of the new control rod group in step 143;
and 5: because there is no overlapping motion between control rod groups. Calculating according to the formulas (4.5) and (4.6) to obtain a new control rod group height, B: step 60, C: step 60, D:23, step (2);
step 6: repeating the processes in the step 4 and the step 5 until the rod position of the control rod group is B: step 60, C: step 60, D: in 26 steps, the effective reactor core multiplication coefficient obtained by the calculation of the whole reactor core is 1.000212, and the absolute difference value of the result and 1 is 21 multiplied by 10 -5 And less than 30 x 10 of the convergence criterion -5 . At this time, the critical calculation of the reactor core is finished, and the rod position of the control rod group obtained finally is the critical rod position of the control rod group obtained by the calculation of the reactor core.
Since the tine effect exists during the control rod movement, the calculation results were verified by explicitly modeling the control rod set rod positions calculated above using the NECP-X program, by comparing B: step 60, C: step 60, D: the height of the rod position in the 26 steps is explicitly modeled and calculated, and the obtained final effective multiplication coefficient is 1.000430, and the difference from the critical value is 43 multiplied by 10 -5 . To ensure that the rod position is indeed the optimal critical rod position, the rod positions for the control rod group are B: step 60, C: step 60, D: the 25-step rod position is also subjected to explicit modeling calculation to obtain the effective multiplication coefficient of 0.999297 and the critical difference of 70.3 multiplied by 10 -5 . Thus, the deviceThe calculated rod position is determined as the optimal critical rod position.
Table 1 is a comparison of the total core count times for the conventional search and criticality approach when the number of control rod sets participating in the search is different. Compared with the traditional method, the method has obvious advantages in processing the search-critical calculation with more control rod groups, greatly reduces the total reactor core calculation times in the calculation process and improves the calculation efficiency.
TABLE 1
Figure SMS_63
In addition, the method can also consider the condition that the overlapped walking rods exist among the control rod groups. Taking the core as an example, the stacking height between B, C, D control rod groups is set to 15 steps.
Wherein step 1, step 2, step 3 and step 4 are the same as above.
And 5: at this time H critical In the case of the number of steps of = 135,
Figure SMS_64
= 3。H critical less than the total height of the control rod set calculated in step four
Figure SMS_65
Equal to the total number of control rod groups. K =2,m =1 is calculated according to the formula (4.13), and then the height size of each bar group is calculated by the formulas (4.11) and (4.9) as B: step 34, C: step 49, D:60 steps;
wherein the content of the first and second substances,
Figure SMS_66
when the control rod group with the number of 1 moves to the maximum height in the stacking step, the height of the control rod group is more than 0;
Figure SMS_67
control rod set numbered 1 moves in cascade to maximum heightIn the time of searching, the total height of the control rod group participating in the critical searching;
m-when obtaining the new control rod group rod position, the rod position height of the control rod group with the number of 1-m is the maximum;
Figure SMS_68
-the number of control rod groups having a rod position height smaller than the maximum value;
step 6: repeating the processes in the step 4 and the step 5 until the rod position of the control rod group is B: step 29, C: step 44, D: in step 59, the effective multiplication coefficient of the reactor core obtained by the calculation of the whole reactor core is 1.000172, and the absolute difference value of the result and 1 is 17 multiplied by 10 -5 And less than 30 x 10 of the convergence criterion -5 . At this time, the critical calculation of the reactor core is finished, and the rod position of the control rod group obtained finally is the critical rod position of the control rod group obtained by the calculation of the reactor core.
The control rod group position obtained by the calculation is subjected to explicit modeling by using a NECP-X program to verify the calculation result, and the final effective multiplication coefficient is 0.999222, and the critical difference is 77.8 multiplied by 10 -5 . To ensure that the rod position is indeed the optimal rod position, the rod positions for the control rod group are B: step 30, C: step 45, D: the bar position of 60 steps is also subjected to explicit modeling calculation to obtain the effective multiplication coefficient of 1.001162 and the critical difference of 116.2 multiplied by 10 -5 . Therefore, the calculated rod position is determined as the optimal critical rod position. The method can be well applied to the calculation of the search critical of the control rod stacking step motion rod, and has higher precision.

Claims (2)

1. A method for simultaneously performing a search by multiple control rod groups in consideration of overlapping steps, comprising: the method comprises the following steps:
step 1: reading geometric information, material information, boundary conditions and control rod group information of a critical core to be simulated;
and 2, step: obtaining information of each control stick group participating in the search threshold according to the control stick group information obtained in the step 1, wherein the information comprises the sequence of each control stick group participating in the search threshold, the minimum and maximum movable heights and the initial position of each control stick group; then, performing neutron transport calculation on the whole reactor core to obtain the effective reactor core multiplication coefficient in an initial state;
and 3, step 3: respectively processing according to the difference between the effective reactor core multiplication coefficient calculated in the step 2 and the absolute difference value of 1; if the absolute difference value of the effective reactor core multiplication coefficient and 1 is greater than the convergence criterion and the effective reactor core multiplication coefficient is greater than 1, inserting all control rod groups of the reactor core participating in the critical search into the reactor core by moving the control rods; if the absolute difference value of the effective reactor core multiplication coefficient and 1 is greater than the convergence criterion and the effective reactor core multiplication coefficient is less than 1, all control rod groups of the reactor core participating in the critical search are lifted out of the reactor core by moving the control rods; if the absolute difference value between the effective multiplication coefficient of the reactor core and 1 is smaller than the convergence criterion, the reactor core reaches a critical state at the moment, and the calculation is finished;
and 4, step 4: carrying out the transport calculation of the whole reactor core again according to the new rod position of the control rod group to obtain the effective value-added factors of the reactor core under the new rod position of the control rod group; calculating to obtain the difference of the total heights of the control rod groups in two adjacent whole core calculations through a formula (4.1), then calculating to obtain the average differential value of the control rod groups in the two whole core calculations through a formula (4.2), and finally carrying out total height interpolation calculation through a formula (4.3) and a formula (4.4) to obtain the total height of a new control rod group;
Figure QLYQS_1
(4.1)
Figure QLYQS_2
(4.2)
Figure QLYQS_3
(4.3)
Figure QLYQS_4
(4.4)
wherein the content of the first and second substances,
Figure QLYQS_5
-the difference between the total height of the control rod groups participating in the search criticality in the jth full core calculation and the total height of the control rod groups participating in the search criticality in the last full core calculation
Figure QLYQS_6
Figure QLYQS_7
-the number of control stick groups participating in the search;
Figure QLYQS_8
-the height of the ith control rod group at the jth full core calculation;
Figure QLYQS_9
-the height of the ith control rod group at j-1 time of the full core calculation;
Figure QLYQS_10
after j times of total core calculation, comparing the average differential value of the control rod groups participating in critical search with the average differential value of the control rod groups obtained in the previous total core calculation;
Figure QLYQS_11
the effective incremental coefficient of the reactor core obtained after the jth total reactor core calculation;
Figure QLYQS_12
the effective incremental coefficient of the reactor core obtained after j-1 th calculation of the whole reactor core;
Figure QLYQS_13
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure QLYQS_14
-total height of control rod groups participating in critical search in jth full core calculation;
Figure QLYQS_15
-total height of control rod groups participating in the search for criticality at the 1 st full core calculation;
Figure QLYQS_16
-height of the ith control rod group at 1 st full core calculation;
and 5: if the overlapping movement of the control rod groups is not considered, numbering N control rod groups according to the sequence of the search critical participation of each control rod group obtained in the step 2 and the maximum height capable of moving of each control rod group as 1-N, and assuming that the heights of m control rod groups reach the maximum at the moment; calculating according to a formula (4.5) and a formula (4.6) to obtain the height of each control rod group in the next total reactor core calculation;
Figure QLYQS_17
(4.5)
Figure QLYQS_18
(4.6)
wherein the content of the first and second substances,
h i+1,j+1 -the height of the (i + 1) th control rod group at the time of the (j + 1) th total core calculation;
Figure QLYQS_19
-the maximum height that the ith control stick set can move;
if the step-stacking movement of the control rod groups is considered, assuming that the step-stacking height of each control rod group during the step-stacking movement is d, the maximum height of the control rod groups capable of moving is equal;
when the control bar group numbered 1 is moved in steps to the maximum height, it is assumed that it is common at this time
Figure QLYQS_20
The height of each control rod group is more than 0; the height sum of each control rod group is obtained through the formula (4.7)
Figure QLYQS_21
Then the total height of the control rod group at the moment is calculated by the formula (4.8);
Figure QLYQS_22
(4.7)
Figure QLYQS_23
(4.8)
wherein the content of the first and second substances,
Figure QLYQS_24
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
Figure QLYQS_25
-the maximum height that the 1 st control stick set can move;
Figure QLYQS_26
-stack height between control rod groups;
Figure QLYQS_27
-the number of control stick groups participating in the search;
Figure QLYQS_28
-the total height of the control rod group participating in the search for criticality when the control rod group numbered 1 moves in steps to the maximum height;
Figure QLYQS_29
-the number of control rod groups with number 1, the control rod group height being greater than 0, when the control rod groups move in steps to the maximum height;
when in use
Figure QLYQS_30
If so, assuming that the height of the n control rod groups is more than 0; the height of the control rod group with the number from n to 1 is shown in a formula (4.9), and the size of n is calculated according to the formula (4.10);
Figure QLYQS_31
(4.9)
when in use
Figure QLYQS_32
When N = N;
otherwise, the formula (4.10) is satisfied
Figure QLYQS_33
(4.10)
Calculated according to equation (4.10):
Figure QLYQS_34
since n is a positive integer and the value range of n in the calculation result of the formula (4.10) is less than 1, only one integer n meets the condition; computingAfter n is obtained, calculating the height of each control rod group through a formula (4.11) and a formula (4.9);
Figure QLYQS_35
(4.11)
wherein the content of the first and second substances,
Figure QLYQS_36
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure QLYQS_37
the jth total core calculation, the th
Figure QLYQS_38
The height of each control rod group;
Figure QLYQS_39
-stack height between control rod groups;
Figure QLYQS_40
-the number of control rod groups having a height greater than 0;
Figure QLYQS_41
-the number of control stick groups participating in the search;
Figure QLYQS_42
when the control rod group with the number of 1 moves to the maximum height in the stacking step, the height of the control rod group is more than 0;
Figure QLYQS_43
when the control rod group with the number of 1 moves to the maximum height in a stacking mode, the total height of the control rod group participating in searching the critical area is increased;
Figure QLYQS_44
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
when H is present j+1 >H critical And is
Figure QLYQS_45
Figure QLYQS_46
When the height of the control rod group reaches the maximum height, the size of m is calculated according to a formula (4.12);
Figure QLYQS_47
(4.12)
calculated according to equation (4.12):
Figure QLYQS_48
since m is a positive integer and the value range of m in the calculation result of the formula (4.12) is less than 1, only one integer m meets the condition, and after m is obtained through calculation, the heights of the rest bar groups are obtained through calculation of the formula (4.10), the formula (4.11) and the formula (4.9);
when H is present j+1 >H critical
Figure QLYQS_49
If the number of the control rod groups with the rod position height smaller than the maximum value is k, calculating to obtain the new rod position of each control rod group according to the following formula (4.13);
Figure QLYQS_50
(4.13)
calculated according to equation (4.13):
Figure QLYQS_51
k is a positive integer, and the value range of k in the calculation result of the formula (4.13) is less than 1, so that only one integer k meets the condition; then calculating the size of m by m = N-k, and finally calculating the heights of the rest control rod groups by a formula (4.11) and a formula (4.9);
wherein the content of the first and second substances,
Figure QLYQS_52
-total height of control rod group participating in critical search in j +1 th total core calculation;
Figure QLYQS_53
-the height of the ith control rod group at the jth full core calculation;
Figure QLYQS_54
-the maximum height that the ith control stick set can move;
Figure QLYQS_55
-stack height between control rod groups;
Figure QLYQS_56
-the number of control bar sets having a height greater than 0;
Figure QLYQS_57
-the number of control stick groups participating in the search;
Figure QLYQS_58
-the number of control rod groups with number 1, the control rod group height being greater than 0, when the control rod groups move in steps to the maximum height;
Figure QLYQS_59
-the total height of the control rod group participating in the search for criticality when the control rod group numbered 1 moves in steps to the maximum height;
Figure QLYQS_60
-the height of the ith control rod group when the control rod group numbered 1 moves in steps to the maximum height;
Figure QLYQS_61
-the number of control rod groups having a rod position height smaller than the maximum value;
step 6: repeating the processes in the step 4 and the step 5 until the absolute difference value between the effective reactor core multiplication coefficient obtained by the calculation of the whole reactor core and 1 is less than the convergence criterion under a certain control rod group height; at the moment, the critical calculation of the reactor core is finished, and the finally obtained height of each control rod group is the critical rod position of the reactor core control rod group.
2. The method of claim 1, wherein the step-by-step simultaneous motion of the plurality of control rod groups is considered as follows: the convergence criterion of the steps 3 and 6 considers the tine effect and the relevant test basis in the control rod moving process, and the convergence criterion is set not to exceed
Figure QLYQS_62
CN202310115637.4A 2023-02-15 2023-02-15 Method for simultaneously moving multiple control rod groups to search critical by considering overlapping steps Active CN115828643B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310115637.4A CN115828643B (en) 2023-02-15 2023-02-15 Method for simultaneously moving multiple control rod groups to search critical by considering overlapping steps

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310115637.4A CN115828643B (en) 2023-02-15 2023-02-15 Method for simultaneously moving multiple control rod groups to search critical by considering overlapping steps

Publications (2)

Publication Number Publication Date
CN115828643A true CN115828643A (en) 2023-03-21
CN115828643B CN115828643B (en) 2023-05-02

Family

ID=85521470

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310115637.4A Active CN115828643B (en) 2023-02-15 2023-02-15 Method for simultaneously moving multiple control rod groups to search critical by considering overlapping steps

Country Status (1)

Country Link
CN (1) CN115828643B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000162365A (en) * 1998-11-26 2000-06-16 Hitachi Ltd Device for calculating reactor core performance of reactor
JP2008026151A (en) * 2006-07-21 2008-02-07 Toshiba Corp Core monitoring device and method
CN107122546A (en) * 2017-04-27 2017-09-01 西安交通大学 A kind of coupling of multiple physics method that presurized water reactor stable state is calculated
CN110427706A (en) * 2019-08-05 2019-11-08 中国核动力研究设计院 A kind of dynamic promotes the search reactor core critical rod position calculation method of fitting order
CN111523233A (en) * 2020-04-23 2020-08-11 西安交通大学 Neutron transport calculation method for three-dimensional pressurized water reactor core
CN112233824A (en) * 2020-09-08 2021-01-15 中广核研究院有限公司 Nuclear power reactor core control rod cascade movement method and nuclear power reactor core
CN113345605A (en) * 2021-04-29 2021-09-03 广西防城港核电有限公司 Fast critical control method for refueling and starting of nuclear reactor
CN114974631A (en) * 2022-04-15 2022-08-30 华能山东石岛湾核电有限公司 First charging and critical method for high-temperature gas cooled reactor
CN115331844A (en) * 2022-09-08 2022-11-11 中国核动力研究设计院 Method and system for measuring value of control rod under subcritical condition of nuclear reactor

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000162365A (en) * 1998-11-26 2000-06-16 Hitachi Ltd Device for calculating reactor core performance of reactor
JP2008026151A (en) * 2006-07-21 2008-02-07 Toshiba Corp Core monitoring device and method
CN107122546A (en) * 2017-04-27 2017-09-01 西安交通大学 A kind of coupling of multiple physics method that presurized water reactor stable state is calculated
CN110427706A (en) * 2019-08-05 2019-11-08 中国核动力研究设计院 A kind of dynamic promotes the search reactor core critical rod position calculation method of fitting order
CN111523233A (en) * 2020-04-23 2020-08-11 西安交通大学 Neutron transport calculation method for three-dimensional pressurized water reactor core
CN112233824A (en) * 2020-09-08 2021-01-15 中广核研究院有限公司 Nuclear power reactor core control rod cascade movement method and nuclear power reactor core
CN113345605A (en) * 2021-04-29 2021-09-03 广西防城港核电有限公司 Fast critical control method for refueling and starting of nuclear reactor
CN114974631A (en) * 2022-04-15 2022-08-30 华能山东石岛湾核电有限公司 First charging and critical method for high-temperature gas cooled reactor
CN115331844A (en) * 2022-09-08 2022-11-11 中国核动力研究设计院 Method and system for measuring value of control rod under subcritical condition of nuclear reactor

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
DAO YU HU; JIE SHENG; ZHU YONG LI; ZHI YONG HONG; ZHI JIAN JIN: "Experimental and numerical study on an HTS air core transformer with pancake structure" *
卢宗健;刘同先;王金雨;吴磊;于颖锐;: "VVER-1000型反应堆停堆及临界模式优化研究" *
郑继业;鲍杰;肖红;付浩;王一鸣;付陟玮;: "AP1000机组首次临界试验验证计算" *
马党伟;刘宙宇;赵晨;贺涛;张文鑫;曹良志;吴宏春;柴晓明;: "NECP-X的多重并行策略及效率优化" *

Also Published As

Publication number Publication date
CN115828643B (en) 2023-05-02

Similar Documents

Publication Publication Date Title
CN104102782B (en) Method for controlling temperatures of reactors of pressurized water reactor nuclear power stations by aid of RBF (radial basis function) neural networks
CN113591954B (en) Filling method of missing time sequence data in industrial system
CN113076684B (en) Intelligent calculation method for transient parameters in rod adjusting process of nuclear reactor core
CN110427706B (en) Method for calculating critical rod position of search reactor core by dynamically increasing fitting order
Stout et al. Determination of optimum fuel loadings in pressurized water reactors using dynamic programming
Tayefi et al. Using Hopfield neural network to optimize fuel rod loading patterns in VVER/1000 reactor by applying axial variation of enrichment distribution
CN115828643A (en) Method for simultaneously moving rod search thresholds of multiple control rod groups in consideration of overlapping steps
CN111274624B (en) Multi-working-condition special-shaped node topology optimization design method based on RBF proxy model
Motoda Optimization of Control Rod Programming and Loading Pattern in Multiregion Nuclear Reactor by the Method of Approximation Programming
Kubiński et al. Optimization of the loading pattern of the PWR core using genetic algorithms and multi-purpose fitness function
CN112987665A (en) Flow shop scheduling method based on reinforcement learning
Dzianisau et al. Macroscopic cross-section generation for nodal code RAST-K using artificial neural network
CN113095596B (en) Photovoltaic power prediction method based on multi-stage Gate-SA-TCN
US6161056A (en) Placement method and apparatus
CN113656852B (en) Method for rapidly generating fine river terrain
CN106611082A (en) Efficient radial basis function supporting point simplification method
Charles et al. Mixed oxide LWR assembly design optimization using differential evolution algorithms
Nam et al. Development of deep convolutional neural network for prediction of cycle maximum pin power peaking factor in pressurized water reactor
Jessee et al. Coupled bundle-core design using fuel rod optimization for boiling water reactors
US20230197302A1 (en) Core performance calculation apparatus
CN110427681B (en) Method for parameterizing shape factor of pressurized water reactor assembly
CN110705054B (en) Method for obtaining resonance group constant for neutron strong absorber
An et al. Nuclear Fuel Batch Optimization for Extra Long Equilibrium Core Design of APR-1400
Ougouag Why and how to Adapt the Two-step Method of LWRS for the Analysis of PBRS
Zaki et al. Modified CANDLE burnup scheme and its application for long life Pb-Bi cooled fast reactor with natural uranium as fuel cycle input

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant