CN115295187A - Hexagonal assembly reactor core physical calculation coordinate conversion method - Google Patents

Hexagonal assembly reactor core physical calculation coordinate conversion method Download PDF

Info

Publication number
CN115295187A
CN115295187A CN202210800391.XA CN202210800391A CN115295187A CN 115295187 A CN115295187 A CN 115295187A CN 202210800391 A CN202210800391 A CN 202210800391A CN 115295187 A CN115295187 A CN 115295187A
Authority
CN
China
Prior art keywords
component
assembly
hexagonal
core
coordinate
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
CN202210800391.XA
Other languages
Chinese (zh)
Other versions
CN115295187B (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.)
Nuclear Power Institute of China
Original Assignee
Nuclear Power Institute of China
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 Nuclear Power Institute of China filed Critical Nuclear Power Institute of China
Priority to CN202210800391.XA priority Critical patent/CN115295187B/en
Publication of CN115295187A publication Critical patent/CN115295187A/en
Application granted granted Critical
Publication of CN115295187B publication Critical patent/CN115295187B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/06—Devices or arrangements for monitoring or testing fuel or fuel elements outside the reactor core, e.g. for burn-up, for contamination
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Plasma & Fusion (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

The invention belongs to the technical field of nuclear reactor core calculation, and particularly discloses a hexagonal assembly core physical calculation coordinate conversion method.

Description

Hexagonal assembly reactor core physical calculation coordinate conversion method
Technical Field
The invention belongs to the technical field of nuclear reactor core calculation, and particularly relates to a hexagonal assembly reactor core physical calculation coordinate conversion method.
Background
In traditional light water reactors such as pressurized water reactors and boiling water reactors, a reactor core is generally constructed by rectangular fuel assemblies, and when reactor core nuclear design is carried out by adopting reactor core physical computing software, the reactor core needs to be constructed by a coordinate system and divided into grids. Because the rectangular components are perfectly aligned with the rectangular coordinates, it is very easy to describe the position and order of each component or grid using a cartesian rectangular coordinate system.
And in the fourth generation reactors such as sodium-cooled fast reactors and lead-bismuth fast reactors, the reactor core is constructed by adopting hexagonal components, compared with rectangular components, due to the angle difference of the hexagonal components, the difficulty of constructing a coordinate system for the reactor core is increased when reactor core nuclear design is carried out by adopting reactor core physical calculation software, and the grid position and the sequence definition mode of the hexagonal components are also obviously different. At present, there are two typical methods for building a core coordinate system and defining a grid position of a hexagonal assembly: (1) A two-dimensional non-rectangular coordinate system is adopted, such as an oblique coordinate system with an included angle of 60 degrees, the upper left corner is a coordinate origin, a hexagonal component represented by a dotted line at the upper left corner does not exist really, the core coordinates are supplemented in a virtual component mode, the grid coordinates are described in (i, j), and a typical pressurized water reactor core and the grid coordinates thereof are shown in FIG. 1; (2) The grid coordinates are described by adopting one-dimensional numbering, and the numbering rule is as follows: the central component is 1, numbering is carried out anticlockwise from the second circle, and after one circle is numbered, the next circle is numbered continuously, so that a one-dimensional serial number coordinate system is formed, and grid coordinates of the one-dimensional serial number coordinate system are shown in figure 2.
Generally, only one coordinate system description mode exists in certain hexagonal assembly core physical calculation software, but in recent years, with the rapid development of hexagonal assembly core physical calculation software and the deepening of multi-professional coupling calculation, two sets of coordinate systems exist in the same software or software modules containing different calculation coordinate systems exist in one calculation system, and in order to better realize data transmission and parameter conversion between software or modules, it is necessary to establish a conversion method between different coordinate systems of a hexagonal assembly core.
Disclosure of Invention
In order to meet the parameter transmission between different coordinate systems in the hexagonal assembly core physical calculation software or between the software, the invention provides a hexagonal assembly core physical calculation coordinate conversion method.
In order to achieve the above purpose, the invention is realized by the following technical scheme:
a physical calculation coordinate conversion method for a reactor core of a hexagonal assembly converts a two-dimensional non-rectangular coordinate system into a one-dimensional serial number coordinate system, and implements the following steps:
s1, determining a coordinate A0 (i 0, j 0) of a hexagonal assembly reactor core central assembly in a two-dimensional non-rectangular coordinate system;
s2, calculating relative coordinates rA (ri, rj) of each component A (i, j) of the hexagonal component reactor core relative to a central component A0 (i 0, j 0), and calculating the distance L (i, j) from each component A (i, j) to the central component A0 (i 0, j 0) according to an oblique coordinate distance calculation formula;
s3, judging a circle layer range in which the distance L (i, j) from each assembly A (i, j) to the central assembly A0 (i 0, j 0) falls according to the distance range from each circle assembly of the hexagonal assembly reactor core to the central assembly A0 (i 0, j 0), and determining the number of circles C (i, j) where each assembly A (i, j) is located;
s4, determining that each component A (i, j) is located on the first edge of the circle layer, namely the edge number B (i, j), according to the relative coordinates rA (ri, rj) of each component A (i, j) obtained in the step S2 and the circle number C (i, j) of each component A (i, j) obtained in the step S3;
s5, determining the number G (i, j) of the assemblies A (i, j) as the number of the edges of the assemblies A (i, j) according to the relative coordinates rA (ri, rj) of the assemblies A (i, j) obtained in the step S2 and the number B (i, j) of the edges of the assemblies A (i, j) obtained in the step S4;
s6, calculating one-dimensional serial number coordinates S (i, j) of each component A (i, j) in the one-dimensional serial number coordinate system according to the number of turns C (i, j) determined in the step S3, the number of sides B (i, j) determined in the step S4 and the number G (i, j) determined in the step S5.
The scheme establishes a calculation algorithm for converting the two-dimensional non-rectangular coordinate system into the one-dimensional serial number coordinate system, provides a specific implementation process for data transmission between reactor core physical software modules or software of different coordinate systems, and can expand the application range and the calculation capacity of reactor core physical calculation software.
Further, in step S1: the coordinates A0 (i 0, j 0) of the hexagonal assembly core center assembly is determined to be A0 (R, R), and R is the total number of turns of the core.
Further, in step S2: relative coordinates rA (ri, rj) of each assembly A (i, j) of the hexagonal assembly core relative to the center assembly A0 (R, R) are calculated, where ri = i-R and rj = j-R.
Further, in step S2: the formula for calculating the distance L (i, j) from each assembly a (i, j) to the center assembly A0 (R, R) is: l (i, j) = (ri ^2+ rj ^2+ ri + rj ^ cos (pi/3)) ^0.5.
Further, in step S3: the distance range from the nth turn assembly to the central assembly A0 (R, R) is
Figure RE-GDA0003850186030000021
When the value of n is determined when the distance L (i, j) from each element a (i, j) to the center element A0 (R, R) falls within the range, a (i, j) is located in the nth turn, i.e., C (i, j) = n.
Further, in step S3: for the six corner point components of the hexagonal component core, when the number of turns is larger than or equal to 8, the number of turns C (i, j) = n +1 of the component A (i, j).
Further, in step S4:
when the relative coordinates rA (ri, rj) of component A (i, j) satisfy both (1): ri ≧ 0, (2) = C (i, j) -1, (3): ri ≠ -rj, component A (i, j) is located on the first side, i.e., B (i, j) =1;
when the relative coordinates rA (ri, rj) of component a (i, j) simultaneously satisfy (1) = -C (i, j) +1, (2) >0, component a (i, j) is located on the second side, i.e., B (i, j) =2;
when the relative coordinates rA (ri, rj) of the component a (i, j) satisfy both (1) = -C (i, j) +1, (2) = ri 0 or (1) </0 and (2) </0, the component a (i, j) is located on the third side, i.e., B (i, j) =3;
when the relative coordinates rA (ri, rj) of component a (i, j) simultaneously satisfy (1) <0, (2) = -C (i, j) +1, (3) ≠ -rj, then component a (i, j) is on the fourth side, i.e., B (i, j) =4;
when the relative coordinates rA (ri, rj) of the component a (i, j) satisfy both (1) ' ri < <0, (2) = -C (i, j) +1, (3) ' = -rj or (1) ' ri <0, (2) = -C (i, j) +1, (3);
when the relative coordinates rA (ri, rj) of component A (i, j) satisfy both (1) ri ≧ 0, (2) > <0, or both (1) ri <0, (2) = C (i, j) -1, (3) = ri =0, component A (i, j) is located on the sixth side, i.e., B (i, j) =6.
Further, in step S5:
when the number of sides B (i, j) =1 of the component a (i, j), the number G (i, j) of components a (i, j) located on the first side is equal to the absolute value of the difference between the relative coordinates rj and 1 in the y direction, i.e., G (i, j) = abs (rj-1);
when the number of sides B (i, j) =2 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the second side is equal to the absolute value of the difference between the number of turns C (i, j) of the component a (i, j) and the relative coordinate ri in the x direction, i.e., G (i, j) = abs (C (i, j) -ri);
when the number of sides B (i, j) =3 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the third side is equal to the absolute value of the difference between the relative coordinates ri in the x direction and 1, that is, G (i, j) = abs (ri-1);
when the number of sides B (i, j) =4 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the fourth side is equal to the absolute value of the sum of the relative coordinates rj and 1 in the y direction, i.e., G (i, j) = abs (rj + 1);
when the number of sides B (i, j) =5 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the fifth side is equal to the absolute value of the sum of the number of turns C (i, j) of the component a and the relative coordinate ri in the x direction, that is, G (i, j) = abs (C (i, j) + ri);
when the number of sides B (i, j) =6 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the sixth side is equal to the absolute value of the sum of the relative coordinates ri in the x direction and 1, that is, G (i, j) = abs (ri + 1).
Further, in step S6:
when the number of turns C (i, j) <3, the calculation formula for converting the two-dimensional non-rectangular coordinate into the one-dimensional serial number coordinate is as follows: s (i, j) =1+ (B (i, j) -1) × (C (i, j) -1) + G (i, j).
Further, in step S6:
when the number of turns C (i, j) is more than or equal to 3, the calculation formula for converting the two-dimensional non-rectangular coordinate into the one-dimensional serial number coordinate is as follows: s (i, j) =1+3 (C (i, j) -2) × (C (i, j) -1) + (B (i, j) -1) × (C (i, j) -1) + G (i, j).
In summary, compared with the prior art, the invention has the following advantages and beneficial effects:
1. aiming at the hexagonal component reactor core, a calculation algorithm for converting a two-dimensional non-rectangular coordinate system into a one-dimensional serial number coordinate system is established, a specific implementation process is provided for data transmission between software modules or software of different coordinate systems, and the method can expand the application range and the calculation capability of reactor core physical calculation software.
2. The invention can be used for the coordinate system conversion and parameter transmission of the reactor cores of lead-bismuth fast reactors and the like which typically adopt hexagonal assemblies during the physical calculation of the reactor cores.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings, a core composed of 3-turn hexagonal assemblies is illustrated as an example:
FIG. 1 shows a schematic diagram of the component coordinates of a hexagonal component core in a two-dimensional non-rectangular coordinate system.
FIG. 2 shows a schematic component coordinate diagram of a hexagonal component core in a one-dimensional serial number coordinate system.
FIG. 3 gives a numbered schematic of three key parameters of the hexagonal assembly in the core.
Figure 4 gives the coordinate distance versus the number of turns.
Detailed Description
To make the objects, technical solutions and advantages of the present invention more apparent, the following detailed description of the principles, features and the like of the present invention is made with reference to the following examples and accompanying drawings, and the exemplary embodiments and descriptions of the present invention are only used for explaining the present invention, and are not used as limiting the scope of the present invention.
Flowcharts or text are used in this specification to illustrate the operational steps performed in accordance with embodiments of the present application. It should be understood that the operational steps in the embodiments of the present application are not necessarily performed in the exact order recited. Rather, the various steps may be processed in reverse order or simultaneously, as desired. Meanwhile, other operations may be added to the processes, or a certain step or several steps of operations may be removed from the processes.
In order to meet the parameter transmission between different coordinate systems in the hexagonal assembly core physical calculation software or between different coordinate systems, the invention provides a hexagonal assembly core physical calculation coordinate conversion method. Specific embodiments of the present invention are given below:
examples
The embodiment provides a method for converting physical calculation coordinates of a core of a hexagonal assembly, which converts a two-dimensional non-rectangular coordinate system into a one-dimensional serial number coordinate system and implements the following steps:
s1, determining a coordinate A0 (i 0, j 0) of a hexagonal assembly reactor core central assembly in a two-dimensional non-rectangular coordinate system;
specifically, in step S1: the coordinates A0 (i 0, j 0) of the hexagonal assembly core center assembly is determined to be A0 (R, R), and R is the total number of turns of the core.
S2, calculating relative coordinates rA (ri, rj) of each component A (i, j) of the hexagonal component reactor core relative to a central component A0 (i 0, j 0), and calculating the distance L (i, j) from each component A (i, j) to the central component A0 (i 0, j 0) according to an oblique coordinate distance calculation formula;
specifically, in step S2: relative coordinates rA (ri, rj) of each assembly a (i, j) of the hexagonal assembly core with respect to the center assembly A0 (R, R) are calculated, where ri = i-R and rj = j-R.
Specifically, in step S2: and calculating the distance L (i, j) from each component A (i, j) to the central component A0 (R, R) according to an oblique coordinate distance calculation formula, wherein L (i, j) = (ri ^2+ rj ^2+ ri + rj cos (pi/3)) ^0.5.
In step S2, the distance from the center of each component to the center component is calculated by using an oblique coordinate formula, and the distance is calculated not by using the actual width of the component but by using the coordinate identifier.
S3, judging a circle layer range in which the distance L (i, j) from each assembly A (i, j) to the central assembly A0 (i 0, j 0) falls according to the distance range from each circle assembly of the hexagonal assembly reactor core to the central assembly A0 (i 0, j 0), and determining the number of circles C (i, j) where each assembly A (i, j) is located;
specifically, in step S3: the distance range from the nth turn assembly to the central assembly A0 (R, R) is
Figure RE-GDA0003850186030000051
When the value of n is determined when the distance L (i, j) from each component a (i, j) to the center component A0 (R, R) falls within this range, a (i, j) is located in the nth turn, i.e., C (i, j) = n. In fig. 4, the distance range between the nth circle component and the central component A0 (R, R) is the range between two circles, that is, the range between two circles
Figure RE-GDA0003850186030000052
Specifically, in step S3: for the six corner point assemblies of the hexagonal assembly core, when the number of turns is larger than or equal to 8, C (i, j) = n +1. In step S3, it is particularly clear whether the number of turns is greater than or equal to 8, and the formula at six angular points changes when the number of turns of the assembly is calculated after the number of turns exceeds 8.
S4, determining that each component A (i, j) is located on the first edge of the circle layer, namely the edge number B (i, j), according to the relative coordinates rA (ri, rj) of each component A (i, j) obtained in the step S2 and the circle number C (i, j) of each component A (i, j) obtained in the step S3;
specifically, in step S4:
when the relative coordinates rA (ri, rj) of component A (i, j) satisfy both (1): ri ≧ 0, (2) = C (i, j) -1, (3): ri ≠ -rj, component A (i, j) is located on the first side, i.e., B (i, j) =1;
when the relative coordinates rA (ri, rj) of component a (i, j) satisfy (1) = -C (i, j) +1, (2) >, at the same time, component a (i, j) is located on the second side, i.e., B (i, j) =2;
when the relative coordinates rA (ri, rj) of the component a (i, j) satisfy both (1) = -C (i, j) +1, (2) = ri 0 or (1) </0 and (2) </0, the component a (i, j) is located on the third side, i.e., B (i, j) =3;
when the relative coordinates rA (ri, rj) of component a (i, j) simultaneously satisfy (1) <0, (2) = -C (i, j) +1, (3) ≠ -rj, component a (i, j) is located on the fourth side, i.e., B (i, j) =4;
when the relative coordinates rA (ri, rj) of the component a (i, j) satisfy both (1) ' ri < <0, (2) = -C (i, j) +1, (3) ' = -rj or (1) ' ri <0, (2) = -C (i, j) +1, (3);
when the relative coordinates rA (ri, rj) of component A (i, j) satisfy both (1) ri ≧ 0, (2) rj >0, or both (1) ri <0, (2) rj = C (i, j) -1, (3) ri =0, component A (i, j) is located on the sixth side, i.e., B (i, j) =6.
S5, determining that each component A (i, j) is the number G (i, j) of the edge according to the relative coordinates rA (ri, rj) of each component A (i, j) acquired in the step S2 and the number B (i, j) of the edges of each component A (i, j) acquired in the step S4;
specifically, in step S5:
when the number of sides B (i, j) =1 of the component a (i, j), the number G (i, j) of components a (i, j) located on the first side is equal to the absolute value of the difference between the relative coordinates rj and 1 in the y direction, i.e., G (i, j) = abs (rj-1);
when the number of sides B (i, j) =2 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the second side is equal to the absolute value of the difference between the number of turns C (i, j) of the component a (i, j) and the relative coordinate ri in the x direction, i.e., G (i, j) = abs (C (i, j) -ri);
when the number of sides B (i, j) =3 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the third side is equal to the absolute value of the difference between the relative coordinates ri in the x direction and 1, that is, G (i, j) = abs (ri-1);
when the number of sides B (i, j) =4 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the fourth side is equal to the absolute value of the sum of the relative coordinates rj in the y direction and 1, that is, G (i, j) = abs (rj + 1);
when the number of sides B (i, j) =5 of the component a (i, j), the number G (i, j) of the component a (i, j) located on the fifth side is equal to the absolute value of the sum of the number of turns C (i, j) of the component a and the relative coordinate ri in the x direction, i.e., G (i, j) = abs (C (i, j) + ri);
when the number of sides B (i, j) =6 of the component a (i, j), the number G (i, j) of the components a (i, j) located on the sixth side is equal to the absolute value of the sum of the relative coordinates ri and 1 in the x direction, i.e., G (i, j) = abs (ri + 1).
In fig. 3, the three key parameters are the number of turns, the number of edges and the number of components in the reactor core respectively, the number of turns refers to the number of turns of the components in the reactor core, the number of edges refers to the number of edges of the components in the number of turns of the reactor core (the number of edges is 1 to 6, the upper right is 1, and the number increases counterclockwise), and the number refers to the number of components on a certain edge (the number increases counterclockwise, and the range is 1-turn-1).
S6, according to the number of turns C (i, j) determined in the step S3, the number of edges B (i, j) determined in the step S4 and the number G (i, j) determined in the step S5, calculating a one-dimensional serial number coordinate S (i, j) of each component A (i, j) in a one-dimensional serial number coordinate system;
specifically, in step S6:
when the number of turns C (i, j) <3, the calculation formula for converting the two-dimensional non-rectangular coordinate into the one-dimensional serial number coordinate is as follows: s (i, j) =1+ (B (i, j) -1) (C (i, j) -1) + G (i, j).
When the number of turns C (i, j) is more than or equal to 3, the calculation formula of converting the two-dimensional non-rectangular coordinate into the one-dimensional serial number coordinate is as follows: s (i, j) =1+3 (C (i, j) -2) × (C (i, j) -1) + (B (i, j) -1) × (C (i, j) -1) + G (i, j). And S6, when only 1 circle or 2 circles are calculated, the calculation formula is changed with the one-dimensional sequence number of more than 2 circles, and the one-dimensional sequence numbers are sorted according to a reverse clock.
The embodiment discloses a hexagonal component reactor core physical calculation coordinate conversion method which comprises the steps of firstly determining the center component coordinates of a hexagonal component reactor core and the two-dimensional coordinates of other components relative to the center component, then determining the number of turns, the number of edges and the number of the components, and finally calculating the one-dimensional serial number coordinates numbered in a counterclockwise sequence to realize the conversion between a two-dimensional coordinate system and a one-dimensional coordinate system, thereby providing a specific implementation process for data transmission between reactor core physical software modules or software of different coordinate systems and being capable of expanding the use range and the calculation capability of reactor core physical calculation software.
Taking the component a (5, 2) as an example, the two-dimensional non-rectangular coordinate a (5, 2) in fig. 1 is converted into the one-dimensional serial number coordinate S (5, 2) =9 in fig. 2, and the following steps are performed:
s1, determining the coordinates of a hexagonal assembly reactor core central assembly as A0 (3, 3);
s2, calculating the relative coordinates of the component A (5, 2) relative to the central component A0 (3, 3) as rA (2, -1); calculate the distance of component A (5, 2) to the center component A0 (3, 3)
Figure RE-GDA0003850186030000081
S3 distance of A (5, 2) to A0 (3, 3) central component
Figure RE-GDA0003850186030000082
The distance range from the 3 rd circle (n = 3) assembly to the center assembly A0 (3, 3)
Figure RE-GDA0003850186030000083
I.e. the number of turns C (5, 2) =3 of a (5, 2);
s4. The relative coordinates rA (2, -1) of the component A (5, 2) satisfy simultaneously (1): ri ≧ 0, (2) = C (5, 2) -1, (3): ri ≠ -rj, thus the number of edges B (5, 2) =1 of the component A (5, 2);
s5, the number B (5, 2) =1 of sides of the component A (5, 2), the number G (5, 2) of the components A (5, 2) on the first side is equal to the absolute value of the difference between the relative coordinates rj and 1 in the y direction, namely G (5, 2) =2;
s6, calculating according to a calculation formula of converting the two-dimensional non-rectangular coordinates into one-dimensional serial number coordinates when the number of turns C (5, 2) ≥ 3 to obtain S (5, 2) =9.
The above-mentioned embodiments, objects, technical solutions and advantages of the present invention are further described in detail, it should be understood that the above-mentioned embodiments are merely preferred embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (10)

1. A physical calculation coordinate conversion method for a reactor core of a hexagonal assembly is characterized in that a two-dimensional non-rectangular coordinate system is converted into a one-dimensional serial number coordinate system, and the following steps are implemented:
s1, determining a coordinate A0 (i 0, j 0) of a hexagonal assembly reactor core central assembly in a two-dimensional non-rectangular coordinate system;
s2, calculating relative coordinates rA (ri, rj) of each component A (i, j) of the hexagonal component reactor core relative to a central component A0 (i 0, j 0), and calculating the distance L (i, j) from each component A (i, j) to the central component A0 (i 0, j 0) according to an oblique coordinate distance calculation formula;
s3, judging a circle layer range in which the distance L (i, j) from each assembly A (i, j) to the central assembly A0 (i 0, j 0) falls according to the distance range from each circle assembly of the hexagonal assembly reactor core to the central assembly A0 (i 0, j 0), and determining the number of circles C (i, j) where each assembly A (i, j) is located;
s4, determining that each component A (i, j) is located on the first edge of the circle layer, namely the edge number B (i, j), according to the relative coordinates rA (ri, rj) of each component A (i, j) obtained in the step S2 and the circle number C (i, j) of each component A (i, j) obtained in the step S3;
s5, determining that each component A (i, j) is the number G (i, j) of the edge where each component A (i, j) is located according to the relative coordinates rA (ri, rj) of each component A (i, j) obtained in the step S2 and the number B (i, j) of the edges of each component A (i, j) obtained in the step S4;
s6, calculating one-dimensional serial number coordinates S (i, j) of each component A (i, j) in the one-dimensional serial number coordinate system according to the number of turns C (i, j) determined in the step S3, the number of sides B (i, j) determined in the step S4 and the number G (i, j) determined in the step S5.
2. The method for converting physical calculation coordinates of a hexagonal assembly core as claimed in claim 1, wherein in step S1: and determining the coordinates A0 (i 0, j 0) of the central assembly as A0 (R, R), wherein R is the total number of turns of the core.
3. The method for physical computation coordinate conversion of a hexagonal assembly core as claimed in claim 2, wherein in step S2: relative coordinates rA (ri, rj) of each component a (i, j) are calculated, where ri = i-R and rj = j-R.
4. The method for physical computation coordinate transformation of a hexagonal assembly core as claimed in claim 3, wherein in step S2: the formula for calculating the distance L (i, j) from each component a (i, j) to the center component A0 (R, R) is: l (i, j) = (ri ^2+ rj ^2+ ri + rj ^ cos (pi/3)) ^0.5.
5. The method for physical computation coordinate conversion of a hexagonal assembly core as claimed in claim 4, wherein in step S3: the distance range from the nth turn assembly to the central assembly A0 (R, R) is
Figure FDA0003737344220000011
The number of turns C (i, j) = n of the assemblies a (i, j) is judged as the value of n when the distance L (i, j) from each assembly a (i, j) to the center assembly A0 (R, R) falls within the range.
6. The method for physical computation coordinate conversion of a hexagonal assembly core as claimed in claim 5, wherein in step S3: for the six corner point components of the hexagonal component core, when the number of turns is larger than or equal to 8, the number of turns C (i, j) = n +1 of the component A (i, j).
7. The method for physical computation coordinate transformation of a hexagonal assembly core as claimed in claim 6, wherein in step S4:
when the relative coordinates rA (ri, rj) of the component A (i, j) simultaneously satisfy (1) ≧ 0, (2) = C (i, j) -1, (3) ≠ -rj, the number of edges B (i, j) =1 of the component A (i, j);
when the relative coordinates rA (ri, rj) of the component a (i, j) satisfy (1) = -C (i, j) +1, (2) >0 at the same time, the number of sides B (i, j) =2 of the component a (i, j);
when the relative coordinates rA (ri, rj) of the modules a (i, j) satisfy both (1) = -C (i, j) +1, (2) = ri 0 or (1) <0 and (2) <0, the number of sides B (i, j) =3 of the modules a (i, j);
when the relative coordinates rA (ri, rj) of component a (i, j) simultaneously satisfy (1) <0, (2) = -C (i, j) +1, (3) ≠ -rj, then the number of edges B (i, j) =4 for component a (i, j);
when the relative coordinates rA (ri, rj) of the component A (i, j) satisfy both (1) < (0), (2) = -C (i, j) +1, (3);
when the relative coordinates rA (ri, rj) of the component A (i, j) satisfy both (1) ≧ 0, (2) > <0, or both (1) < ≧ 0, (2) = C (i, j) -1, (3) = ri =0, the number of sides B (i, j) =6 of the component A (i, j).
8. The method for physical computation coordinate transformation of a hexagonal assembly core as claimed in claim 7, wherein in step S5:
when the number of sides B (i, j) =1 for the component a (i, j), the number G (i, j) = abs (rj-1) for the component a (i, j);
when the number of sides B (i, j) =2 for the component a (i, j), the number G (i, j) = abs (C (i, j) -ri for the component a (i, j);
when the number of sides B (i, j) =3 of the component a (i, j), the number G (i, j) = abs (ri-1) of the component a (i, j);
when the number of sides B (i, j) =4 for the component a (i, j), the number G (i, j) = abs (rj + 1) for the component a (i, j);
when the number of sides B (i, j) =5 for the component a (i, j), the number G (i, j) = abs (C (i, j) + ri) for the component a (i, j);
when the number of sides B (i, j) =6 for the component a (i, j), the number G (i, j) = abs (ri + 1) for the component a (i, j).
9. The method for physical computation coordinate conversion of a hexagonal assembly core of claim 8, wherein in step S6:
when the number of turns C (i, j) <3, the two-dimensional non-rectangular coordinate is converted into a one-dimensional serial number coordinate by the following formula: s (i, j) =1+ (B (i, j) -1) × (C (i, j) -1) + G (i, j).
10. The method for physical computation coordinate transformation of a hexagonal assembly core as claimed in claim 9, wherein in step S6:
when the number of turns C (i, j) is more than or equal to 3, the calculation formula for converting the two-dimensional non-rectangular coordinate into the one-dimensional serial number coordinate is as follows: s (i, j) =1+3 (C (i, j) -2) × (C (i, j) -1) + (B (i, j) -1) × (C (i, j) -1) + G (i, j).
CN202210800391.XA 2022-07-08 2022-07-08 Method for converting physical calculation coordinates of hexagonal assembly reactor core Active CN115295187B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210800391.XA CN115295187B (en) 2022-07-08 2022-07-08 Method for converting physical calculation coordinates of hexagonal assembly reactor core

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210800391.XA CN115295187B (en) 2022-07-08 2022-07-08 Method for converting physical calculation coordinates of hexagonal assembly reactor core

Publications (2)

Publication Number Publication Date
CN115295187A true CN115295187A (en) 2022-11-04
CN115295187B CN115295187B (en) 2024-03-12

Family

ID=83823208

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210800391.XA Active CN115295187B (en) 2022-07-08 2022-07-08 Method for converting physical calculation coordinates of hexagonal assembly reactor core

Country Status (1)

Country Link
CN (1) CN115295187B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1184060A (en) * 1997-09-01 1999-03-26 Toshiba Corp Method and device for calculating output distribution of fast reactor
JP2007278991A (en) * 2006-04-11 2007-10-25 Nuclear Fuel Ind Ltd Control method for three-dimensional mesh system for parallel calculation
US20110029294A1 (en) * 2008-02-18 2011-02-03 Areva Np Computer Implemented Method for Modelizing a Nuclear Reactor Core and a Corresponding Computer Program Product
US20110268239A1 (en) * 2010-04-30 2011-11-03 David Jerome Krieg Method of calibrating excore detectors in a nuclear reactor
US20150310946A1 (en) * 2012-11-30 2015-10-29 Electricite De France Improved determination of positions of fuel assembly elements
KR101653155B1 (en) * 2015-03-06 2016-09-02 한국원자력연구원 Test apparatus and method for fast reactor's core physics test
CN108647402A (en) * 2018-04-19 2018-10-12 西安交通大学 A method of it obtaining hexagonal geometry pressurized water reactor core and reflects layer parameter
CN110867262A (en) * 2019-11-21 2020-03-06 中国核动力研究设计院 Liquid metal cooling reactor based on improvement of fuel utilization rate and management method
CN111145924A (en) * 2019-12-31 2020-05-12 西安交通大学 Power density distribution statistical method for porous hexagonal prism type nuclear fuel element
CN111696197A (en) * 2020-05-27 2020-09-22 北京科技大学 Ray generation and tracking method for neutron transport calculation of hexagonal reactor core of fast reactor

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1184060A (en) * 1997-09-01 1999-03-26 Toshiba Corp Method and device for calculating output distribution of fast reactor
JP2007278991A (en) * 2006-04-11 2007-10-25 Nuclear Fuel Ind Ltd Control method for three-dimensional mesh system for parallel calculation
US20110029294A1 (en) * 2008-02-18 2011-02-03 Areva Np Computer Implemented Method for Modelizing a Nuclear Reactor Core and a Corresponding Computer Program Product
US20110268239A1 (en) * 2010-04-30 2011-11-03 David Jerome Krieg Method of calibrating excore detectors in a nuclear reactor
US20150310946A1 (en) * 2012-11-30 2015-10-29 Electricite De France Improved determination of positions of fuel assembly elements
KR101653155B1 (en) * 2015-03-06 2016-09-02 한국원자력연구원 Test apparatus and method for fast reactor's core physics test
CN108647402A (en) * 2018-04-19 2018-10-12 西安交通大学 A method of it obtaining hexagonal geometry pressurized water reactor core and reflects layer parameter
CN110867262A (en) * 2019-11-21 2020-03-06 中国核动力研究设计院 Liquid metal cooling reactor based on improvement of fuel utilization rate and management method
CN111145924A (en) * 2019-12-31 2020-05-12 西安交通大学 Power density distribution statistical method for porous hexagonal prism type nuclear fuel element
CN111696197A (en) * 2020-05-27 2020-09-22 北京科技大学 Ray generation and tracking method for neutron transport calculation of hexagonal reactor core of fast reactor

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
安萍;刘东;潘俊杰;赵文博;芦;: "稳态堆芯多物理耦合***CSSS V1.0的研发", 原子能科学技术, no. 05, pages 863 - 868 *
徐李等: "快堆三维六角形节块法输运计算研究", 原子能科学技术, 28 February 2013 (2013-02-28), pages 161 - 165 *

Also Published As

Publication number Publication date
CN115295187B (en) 2024-03-12

Similar Documents

Publication Publication Date Title
CN103150424B (en) Method for acquiring fine distribution of reactor core three dimensional neutron flux density of reactor
CN102208025B (en) Method for correcting geometric distortion of text image
Karian et al. Fitting the generalized lambda distribution to data: a method based on percentiles
CN111046505B (en) Axial flow pump spoke parameter optimization design method based on response surface model
CN110501733B (en) Self-adaptive grid VRS generation and service method
CN111209690B (en) Modeling method for random distribution of TRISO fuel particles in FCM fuel element
CN114387421A (en) Method for acquiring reactor core parameters of pebble bed reactor based on CMFD accelerated transport calculation
CN103218512A (en) Method for acquiring neutron angular flux density in nuclear fuel assembly
CN102385242A (en) Method and system for manufacturing mask
CN109858143B (en) Method for generating special-shaped door and window model on arc wall
CN105512387A (en) Fine power reconstruction method for fuel assembly
CN116011089A (en) Turbine guide vane laminate cooling structure design method and device, terminal and storage medium
CN115295187A (en) Hexagonal assembly reactor core physical calculation coordinate conversion method
CN113672849A (en) One-step transportation calculation method and system based on axial flux expansion
CN117610354A (en) Method and device for mapping parallelepiped grid
CN114399428A (en) Planar artificial element dimension reduction processing method
CN110542410B (en) Geodetic survey result conversion method based on grid model
CN114528719A (en) Online energy group compression method for pressurized water reactor based on two-dimensional reactor core
CN103544654A (en) Method for determining local minimal solution and searching global minimal solution of power grid economic dispatching
CN114970293B (en) Coarse mesh diffusion coefficient calculation method based on two-dimensional one-dimensional coupling Fourier analysis
Schiff et al. Many-fold moment method
CN111444655A (en) Static and dynamic reinforcement method for hydraulic building structure
CN115221813B (en) Nuclear-thermal-force cross-dimension coupling method
CN112632825B (en) Electrostatic field smooth finite element numerical algorithm based on finite element super-convergence
CN116205093A (en) Structural topology optimization design method considering load multi-peak uncertainty

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