CN110704806B - Rapid online calculation method for one-dimensional cylindrical geometric collision probability - Google Patents
Rapid online calculation method for one-dimensional cylindrical geometric collision probability Download PDFInfo
- Publication number
- CN110704806B CN110704806B CN201910884705.7A CN201910884705A CN110704806B CN 110704806 B CN110704806 B CN 110704806B CN 201910884705 A CN201910884705 A CN 201910884705A CN 110704806 B CN110704806 B CN 110704806B
- Authority
- CN
- China
- Prior art keywords
- collision probability
- total cross
- macroscopic
- section
- macroscopic total
- 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.)
- Active
Links
Images
Classifications
-
- 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
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Abstract
A quick online calculation method for geometric collision probability of a one-dimensional cylinder is characterized in that a range from a minimum macroscopic total section to a maximum macroscopic total section value of a resonance ring is divided into n intervals according to an equal ratio sequence, the first macroscopic total section interval is dispersed, a parallel line cutting method is used for calculating collision probability, and a collision probability table is established; when the collision probability is calculated for each superfine group, if the macroscopic total cross sections of all the resonance rings are in the first macroscopic total cross section interval, linear interpolation is adopted, and if the parallel line cutting method is not used, the collision probability is calculated on line; compared with a method of directly adopting parallel line cutting, the method has the advantages that the precision is basically not lost, and the efficiency is obviously improved; compared with the mode that all superfine groups adopt collision probability interpolation, the collision probability table used by the invention can be calculated on line, is more convenient and has wider applicability, and saves the memory; and the first macroscopic total cross section interval is subjected to special dispersion, so that interpolation intervals of the macroscopic total cross sections are avoided being searched one by one.
Description
Technical Field
The invention relates to the field of nuclear reactor physical computation, in particular to a rapid online computing method for a one-dimensional cylindrical geometric collision probability.
Background
The resonance self-shielding calculation provides a plurality of group resonance group constants for nuclear reactor neutron calculation, the group constants correspond to average multi-group resonance sections in each fuel area, and the group constants are constant coefficient items in a neutron transport equation, so that accurate resonance calculation is the basis of accurate neutron calculation.
The superfine group resonance calculation method is to disperse the resonance energy interval into several superfine groups and solve the neutron moderation equation in each superfine group to obtain precise superfine group energy spectrum for merging sections of multiple groups. The superfine group resonance calculation method can simultaneously consider the spatial self-shielding effect and the resonance interference effect to obtain an accurate resonance self-shielding calculation section. Solving the neutron moderation equation of each ultrafine group requires first calculating the collision probability related to geometry and materials. The collision probability of tens of thousands of groups and even hundreds of thousands of groups of ultrafine groups is low in efficiency during on-line calculation, so that the use of the ultrafine group resonance calculation method is limited, and the actual resonance calculation requirements cannot be met. The collision probability interpolation table established in advance can improve the calculation efficiency of the collision probability through interpolation, but the manufacture of the collision probability table is often limited by geometry and materials, and the universality is not high.
Disclosure of Invention
In order to improve the efficiency of the superfine group resonance calculation method, the invention aims to provide a rapid online calculation method for the geometric collision probability of a one-dimensional cylinder, which couples the advantages of an online calculation method based on a parallel line cutting method and an interpolation calculation method based on a collision probability table and provides a high-precision and high-efficiency online calculation method for the collision probability.
In order to achieve the purpose, the invention adopts the following technical scheme to implement:
a fast online calculation method for geometric collision probability of a one-dimensional cylinder mainly comprises the steps of building and interpolating a collision probability table and performing online calculation based on a parallel line cutting method. The method comprises the following steps:
step 1: determining whether the nuclear species is a resonance ring or not according to the composition of each ring nuclear species of the one-dimensional cylindrical geometry; for each resonance ring, calculating the superfine group macroscopic total cross sections corresponding to all superfine groups, and determining the maximum value and the minimum value in all the superfine group macroscopic total cross sections;
step 2: dividing the range between the maximum value and the minimum value in all the superfine group macroscopic total cross sections of the 1 st resonance ring obtained in the step (1) into n macroscopic total cross section intervals according to an equal ratio sequence;
and step 3: dispersing the 1 st macroscopic total cross section interval of the first resonance ring obtained in the step 2 into a plurality of macroscopic total cross section points, and setting the interval between the dispersed macroscopic total cross section points and the points to be d cm-1(ii) a Repeating the operations from the step 2 to the step 3 for all the rest resonance rings; taking a typical macroscopic total cross section for the non-resonant circular rings, wherein the independent variable is the macroscopic total cross section of all the resonant circular rings, and the collision probability is taken as a dependent variable; calculating collision probability on line by using a parallel line cutting method, and establishing a multi-dimensional collision probability table, wherein the dimension is the total number of resonance rings;
and 4, step 4: based on the steps 1-3, the collision probability table is already manufactured on line; when solving the neutron moderation equation, the collision probability in each superfine group needs to be obtained; the collision probability is rapidly obtained as follows: for the first superfine group, judging whether the macroscopic total cross sections of all the resonance circular rings are in the range of the corresponding first macroscopic total cross section interval; if not, calculating the collision probability on line by using a parallel line cutting method, and if so, obtaining the collision probability by adopting linear interpolation of a collision probability table; and (4) repeatedly executing the operation of the step (4) on all the remaining superfine groups, thereby quickly obtaining the collision probability of all the superfine groups.
Compared with the online calculation of the collision probability based on the parallel line cutting method and the interpolation calculation based on the collision probability table, the method has the following innovation points:
1. compared with the on-line calculation of the collision probability based on the parallel line cutting method, the method has the characteristic of higher efficiency under the condition of keeping the precision.
2. Compared with the interpolation calculation method based on the collision probability table, the method has the advantages that the collision probability table is calculated on line aiming at specific problems, the convenience is realized, the applicability is higher, and the memory is saved.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
FIG. 2 is a schematic view of a macro total cross-section of an ultrafine cluster.
Detailed Description
The present invention will be described in further detail with reference to the accompanying drawings and specific embodiments.
The method disclosed by the invention is coupled with an online calculation collision probability based on a parallel line cutting method and an interpolation calculation method based on a collision probability table, wherein the collision probability in an ultra-fine group integral neutron transport equation is derived from the two methods. The invention relates to a rapid online calculation method for geometric collision probability of a one-dimensional cylinder, the flow of which is shown in figure 1, and the method comprises the following steps:
1. and determining whether the ring is a resonance ring according to the composition of each ring element of the one-dimensional cylindrical geometry. And for each resonance ring, calculating the superfine group macroscopic total cross sections corresponding to all superfine groups, and determining the maximum value and the minimum value in all the superfine group macroscopic total cross sections.
2. As shown in fig. 2, the range between the maximum value and the minimum value in the macroscopic total cross section of all the ultrafine groups of the 1 st resonant ring obtained in step 1Dividing the macro total section into n (n takes 2 or 3) sections according to the geometric series, wherein the geometric series ratio calculation formula of the 1 st resonance ring is as follows:
in the formula:
n-number of divisions of macroscopic Total Cross section
R1-geometric array ratio of the 1 st resonant ring;
3. The 1 st macroscopic total section interval obtained by the calculation in the step 2The dispersion is a plurality of macroscopic total cross section points, and the distance between the discrete macroscopic total cross section points and the points is set to be d cm-1(d is less than or equal to 0.1), and the first macroscopic total cross-sectional point value after dispersion is as follows:
in the formula:
floor-the operator of rounding down the value in parentheses;
the value of the j-th macroscopic total cross section point after the 1 st resonance ring is dispersed isRepeating the steps 2 to 3 for all the remaining resonant rings. And taking a typical macroscopic total cross section for the non-resonant circular rings, wherein the independent variable is the macroscopic total cross section of all the resonant circular rings, and the dependent variable is the collision probability. And calculating the collision probability on line by using a parallel line cutting method, and establishing a multi-dimensional collision probability table, wherein the dimension is the total number of the resonance rings.
4. Based on steps 1-3, the collision probability table is already manufactured on line. In solving the neutron moderation equation, the collision probability within each ultrafine cluster is needed. The collision probability is rapidly obtained as follows: for the first superfine group, judging whether the macroscopic total cross sections of all the resonance circular rings are in the first macroscopic total cross section intervalAnd (4) the following steps. If not, using a parallel line cutting method to calculate the collision probability on line, if so, adopting linear interpolation of a collision probability table, wherein the calculation formula of the interpolation coefficient corresponding to the ith resonance ring is as follows:
in the formula:
ai,bi-the ith resonance ring corresponds to a linear interpolation coefficient;
for a one-dimensional cylinder with two resonance rings, the linear interpolation calculation collision probability formula is as follows:
in the formula:
Pi→j-probability of collision from the ith to the jth ring
aj,bj-the jth resonance ring corresponds to a linear interpolation coefficient
ai,bi-the ith resonance ring corresponds to a linear interpolation coefficient
li-the ith resonance ring finds the intermediate index variable of the macroscopic total section interpolation interval
The i-th circular ring has a total macroscopic cross section ofThe total macroscopic cross section of the jth ring isProbability of collision from the ith to the jth ring
The i-th circular ring has a total macroscopic cross section ofThe total macroscopic cross section of the jth ring isProbability of collision from the ith to the jth ring
The i-th circular ring has a macroscopic total cross sectionThe total macroscopic cross section of the jth circular ring isProbability of collision from the ith to the jth ringThe i-th circular ring has a total macroscopic cross section ofThe total macroscopic cross section of the jth ring isThe probability of collision from the ith circle to the jth circle. And (4) repeatedly executing the operation of the step (4) on all the superfine groups, thereby quickly obtaining the collision probability of all the superfine groups.
Claims (3)
1. A fast online calculation method aiming at the geometric collision probability of a one-dimensional cylinder is characterized by comprising the following steps: the method mainly comprises the steps of establishing and interpolating a collision probability table and performing online calculation based on a parallel line cutting method; the method comprises the following steps:
step 1: determining whether the nuclear species is a resonance ring or not according to the composition of each ring nuclear species of the one-dimensional cylindrical geometry; for each resonance ring, calculating the superfine group macroscopic total cross sections corresponding to all superfine groups, and determining the maximum value and the minimum value in all the superfine group macroscopic total cross sections;
step 2: dividing the range between the maximum value and the minimum value in all the superfine group macroscopic total cross sections of the 1 st resonance ring obtained in the step (1) into n macroscopic total cross section intervals according to an equal ratio sequence;
and step 3: dispersing the 1 st macroscopic total cross section interval of the first resonance ring obtained in the step 2 into a plurality of macroscopic total cross section points, and setting the interval between the dispersed macroscopic total cross section points and the points to be d cm-1(ii) a Repeating the operations from the step 2 to the step 3 for all the rest resonance rings; taking a typical macroscopic total cross section for the non-resonant circular rings, wherein the independent variable is the macroscopic total cross section of all the resonant circular rings, and the collision probability is taken as a dependent variable; calculating collision probability on line by using a parallel line cutting method, and establishing a multi-dimensional collision probability table, wherein the dimension is the total number of resonance rings; and 4, step 4: based on the steps 1-3, the collision probability table is already manufactured on line; when solving the neutron moderation equation, the collision probability in each superfine group needs to be obtained; the collision probability is rapidly obtained as follows: for the first superfine group, judging whether the macroscopic total cross sections of all the resonance circular rings are in the range of the corresponding first macroscopic total cross section interval; if not, calculating the collision probability on line by using a parallel line cutting method, and if so, obtaining the collision probability by adopting linear interpolation of a collision probability table; and (4) repeatedly executing the operation of the step (4) on all the remaining superfine groups, thereby quickly obtaining the collision probability of all the superfine groups.
2. The fast online calculation method for the collision probability of the one-dimensional cylindrical geometry according to claim 1, characterized in that: the value of n in step 2 is 2 or 3.
3. The fast online calculation method for the collision probability of the one-dimensional cylindrical geometry according to claim 1, characterized in that: d in the step 3 is less than or equal to 0.1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910884705.7A CN110704806B (en) | 2019-09-19 | 2019-09-19 | Rapid online calculation method for one-dimensional cylindrical geometric collision probability |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910884705.7A CN110704806B (en) | 2019-09-19 | 2019-09-19 | Rapid online calculation method for one-dimensional cylindrical geometric collision probability |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110704806A CN110704806A (en) | 2020-01-17 |
CN110704806B true CN110704806B (en) | 2021-03-02 |
Family
ID=69194848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910884705.7A Active CN110704806B (en) | 2019-09-19 | 2019-09-19 | Rapid online calculation method for one-dimensional cylindrical geometric collision probability |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110704806B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114491906B (en) * | 2020-10-27 | 2023-10-31 | 中国核动力研究设计院 | Equivalent homogenization double non-uniformity calculation method based on cylinder random model |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177154A (en) * | 2013-02-05 | 2013-06-26 | 西安交通大学 | Method for acquiring nuclear fuel assembly resonance parameters |
CN105373667A (en) * | 2015-11-27 | 2016-03-02 | 西安交通大学 | Multi-group section perturbation method for uncertainty analysis of reactor physics calculation |
CN106096182A (en) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | A kind of reactor embedded Resonance self-shielding computational methods |
CN107145721A (en) * | 2017-04-20 | 2017-09-08 | 西安交通大学 | A kind of mixing computational methods for obtaining the few group cross-section parameter of fast neutron reactor |
CN109493924A (en) * | 2018-11-13 | 2019-03-19 | 西安交通大学 | A method of obtaining the effective multigroup cross section of FCM fuel |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107038293B (en) * | 2017-04-05 | 2020-05-15 | 西安交通大学 | Resonance self-screen calculation method for plate-shaped fuel based on equivalent one-dimensional plate model |
CN107092785B (en) * | 2017-04-05 | 2020-04-10 | 西安交通大学 | Method for obtaining resonance group constant for nuclear reactor dual heterogeneous fuel |
-
2019
- 2019-09-19 CN CN201910884705.7A patent/CN110704806B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103177154A (en) * | 2013-02-05 | 2013-06-26 | 西安交通大学 | Method for acquiring nuclear fuel assembly resonance parameters |
CN105373667A (en) * | 2015-11-27 | 2016-03-02 | 西安交通大学 | Multi-group section perturbation method for uncertainty analysis of reactor physics calculation |
CN106096182A (en) * | 2016-06-24 | 2016-11-09 | 西安交通大学 | A kind of reactor embedded Resonance self-shielding computational methods |
CN107145721A (en) * | 2017-04-20 | 2017-09-08 | 西安交通大学 | A kind of mixing computational methods for obtaining the few group cross-section parameter of fast neutron reactor |
CN109493924A (en) * | 2018-11-13 | 2019-03-19 | 西安交通大学 | A method of obtaining the effective multigroup cross section of FCM fuel |
Non-Patent Citations (4)
Title |
---|
Accurate resonance calculation method coupling simplified embedded self-shielding method and ultra-fine group method;Tiejun Zu et.al;《Annals of Nuclear Energy》;20180630;期刊第198-206页 * |
Multi-group effective cross section calculation method for Fully Ceramic Micro-encapsulated fuel;Wen Yin et.al;《Annals of Nuclear Energy》;20181028;期刊第26-37页 * |
RESEARCH ON ACCELERATION OF ULTRA-FINE-GROUP CALCULATION BASED ON METHOD OF CHARACTERISTIC;Shuai Qin et.al;《The Proceedings of the International Conference on Nuclear Engineering》;20190524;正文第1-9页 * |
二维任意几何耦合共振计算方法研究;贺雷等;《核动力工程》;20130731;第34卷;期刊第5-9页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110704806A (en) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111340359B (en) | Comprehensive evaluation method for multi-energy complementary distributed energy system | |
CN105046079B (en) | A kind of field mouthful test design method based on the optimal interior table designs of D- | |
CN111199016A (en) | DTW-based improved K-means daily load curve clustering method | |
CN110689195A (en) | Power daily load prediction method | |
CN110580935B (en) | Method for acquiring full-stack effective resonance self-shielding cross section | |
CN110704806B (en) | Rapid online calculation method for one-dimensional cylindrical geometric collision probability | |
CN103047959B (en) | A kind of flat form error detection method based on entropy theory towards Fine Boring | |
CN105846448B (en) | Method for determining reactive compensation capacity of power distribution network based on random matrix theory | |
CN107038293B (en) | Resonance self-screen calculation method for plate-shaped fuel based on equivalent one-dimensional plate model | |
CN110765703B (en) | Wind power plant aggregation characteristic modeling method | |
CN109544006A (en) | A method of evaluation electric automobile charging station economic benefit | |
CN108459993B (en) | Complex high-dimensional system optimization method based on rapid peak-tracking sampling | |
CN111553568A (en) | Line loss management method based on data mining technology | |
CN112308298A (en) | Multi-scenario performance index prediction method and system for semiconductor production line | |
CN108073074A (en) | A kind of assembling quality evaluating method based on Machine Tool Feeding System kinetic characteristic | |
CN113592538A (en) | Method and system for calculating approximate construction cost | |
CN110276478B (en) | Short-term wind power prediction method based on segmented ant colony algorithm optimization SVM | |
CN108446452A (en) | A kind of mixed-flow pump impeller Robust Optimal Design | |
CN109193791B (en) | Wind power convergence tendency state-based quantification method based on improved shape value | |
Wijuniamurti et al. | Agglomerative Nesting (AGNES) Method and Divisive Analysis (DIANA) Method For Hierarchical Clustering On Some Distance Measurement Concepts | |
CN107066741A (en) | Optimization Design for Aerodynamic method based on data mining | |
CN115982872A (en) | Design method for performance of rolling bearing | |
CN115173421A (en) | Probability optimal power flow calculation method based on progressive Latin hypercube sampling | |
CN110555470B (en) | Oil sprayer grouping method | |
CN114004525A (en) | Electric energy substitution benefit evaluation method facing user side based on fuzzy comprehensive evaluation method |
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 |