CN115048810A - Three-dimensional neutron transport equation calculation method and system based on linear source acceleration - Google Patents
Three-dimensional neutron transport equation calculation method and system based on linear source acceleration Download PDFInfo
- Publication number
- CN115048810A CN115048810A CN202210815861.XA CN202210815861A CN115048810A CN 115048810 A CN115048810 A CN 115048810A CN 202210815861 A CN202210815861 A CN 202210815861A CN 115048810 A CN115048810 A CN 115048810A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- equation
- flux
- source
- neutron transport
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 71
- 230000001133 acceleration Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 45
- 230000004907 flux Effects 0.000 claims description 109
- 238000010168 coupling process Methods 0.000 claims description 11
- 230000004992 fission Effects 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 230000008878 coupling Effects 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 4
- 238000000265 homogenisation Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009257 reactivity Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
The invention discloses a three-dimensional neutron transport equation calculation method and a three-dimensional neutron transport equation calculation system based on linear source acceleration, which relate to the technical field of nuclear reactor core design and are improved on the basis of a traditional flat source approximation method; the plane source approximation is eliminated theoretically, the number of grids of the plane source area is reduced, and the calculation efficiency is improved.
Description
Technical Field
The invention relates to the technical field of nuclear reactor core design, in particular to a linear source acceleration-based three-dimensional neutron transport equation calculation method and system.
Background
The reactor physical analysis calculation is used as the basis of the nuclear reactor system analysis calculation, and the reactor core reactivity and the full-reactor fine power distribution are obtained by solving a neutron transport equation. In order to rapidly develop the research and development of the advanced nuclear power reactor core, advanced high-precision reactor physical design software needs to be researched and developed. In order to simulate the core with a complex structure, research on a 'one-step' reactor physical calculation method based on an accurate physical model and a fine geometric modeling is widely carried out at home and abroad. The neutron angular flux of the three-dimensional neutron transport equation contains 7 dependent variables (3 dimensions in space, 2 dimensions in angle, 1 dimension in energy and 1 dimension in time), and accurate numerical simulation is very difficult.
The method for directly solving the three-dimensional neutron transport equation by the one-step method has the advantages of large calculation amount, high memory consumption and difficulty in realization under the existing calculation condition, so that a two-dimensional one-dimensional method is provided, the direct three-dimensional solution is converted into axial one-dimensional solution and radial two-dimensional solution respectively, and the coupling is performed through a leakage item, so that the calculation requirement of the one-step method is reduced. In the two-dimensional one-dimensional coupling method, flat source approximation is adopted, in order to obtain an accurate calculation result, a fine flat source region grid needs to be divided, the calculation efficiency is reduced, and the calculation requirement for solving a three-dimensional neutron transport equation by a one-step method is reduced by coupling through a leakage item. However, in the conventional two-dimensional one-dimensional coupling method, in order to obtain an accurate calculation result, a fine flat source area grid needs to be divided, so that the calculation efficiency is linearly reduced.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: in the traditional two-dimensional one-dimensional coupling method, in order to obtain an accurate calculation result, a fine flat source area grid needs to be divided, so that the calculation efficiency is linearly reduced; the invention aims to provide a three-dimensional neutron transport equation calculation method and system based on linear source acceleration, which are improved on the traditional plain form, carry out flux expansion based on the linear source form to calculate three-dimensional neutron flux and high-order standard flux moment, and solve the problem of low calculation efficiency caused by the fact that a plain source region grid is divided densely in the traditional two-dimensional one-dimensional method.
The invention is realized by the following technical scheme:
the invention provides a linear source acceleration-based three-dimensional neutron transport equation calculation method, which comprises the following steps of:
establishing a three-dimensional neutron transport equation;
converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
respectively solving the one-dimensional equation and the two-dimensional equation to obtain a characteristic value and a three-dimensional neutron flux of the reactor core; and when the two-dimensional equation is solved, flux expansion is carried out on the basis of a linear source form to calculate three-dimensional neutron flux and high-order standard flux moment.
The working principle of the scheme is as follows: in the traditional two-dimensional one-dimensional coupling method, in order to obtain an accurate calculation result, a fine flat source region grid needs to be divided, so that the calculation efficiency is linearly reduced, the invention aims to provide a three-dimensional neutron transport equation calculation method and a system based on linear source acceleration, the method is improved in the traditional plain form, the three-dimensional neutron flux and the high-order standard flux moment are calculated by performing flux expansion based on the linear source form, the high-order expansion is performed on the space through the linear source form, the high-order standard flux moment is obtained at the same time, and the approximation introduced by the flat source method is eliminated, so that the problem of low calculation efficiency caused by the division of the denser flat source region grid in the traditional two-dimensional one-dimensional method is solved; the scheme theoretically eliminates the flat source approximation introduced in the two-dimensional one-dimensional method, reduces the number of grids in the flat source region and improves the calculation efficiency.
The further optimization scheme is that the three-dimensional neutron transport equation is as follows:
wherein m represents an angle, g represents an energy group,. phi g,m (x, y, z) represents angular flux, x, y, z represent x, y, z coordinates of the location in space, ξ, respectively m Representing the cosine of the angle of incidence of the azimuth and the x-axis, sigma t,g (r) represents the total cross-section, η represents the amplitude sine, μ represents the polar cosine, Q g (x, y, z) represents the total source term.
The further optimization scheme is that the method for acquiring the one-dimensional equation and the two-dimensional equation comprises the following steps:
and (3) integrating the radial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a one-dimensional equation:
and (3) integrating the axial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a two-dimensional equation:
wherein psi g,m,i,j (z) angular flux, Q, of the z-th layer representing the radial (i, j) position angle m energy group g g,i,j (z) a one-dimensional total source term, Q, representing the radial (i, j) position g (x, y) represents a two-dimensional total source term,
indicating a radial leakage term, # g,m (x, y) denotes radial angular flux, Σ t,g,i,j (z) represents a one-dimensional total cross section, Sigma t,g (x, y) represents a two-dimensional total cross section,an axial leakage term is indicated.
The further optimization scheme is that the method for acquiring the characteristic value and the three-dimensional neutron flux of the reactor core comprises the following steps:
calculating fission sources and scattering sources;
on the basis of a fission source and a scattering source, respectively solving a one-dimensional equation and a two-dimensional equation to obtain radial flow, axial flow and radial flux;
performing CMFD iterative updating based on the radial flow, the axial flow and the radial flux to obtain a reactor core characteristic value and a three-dimensional neutron flux;
the CMFD iterative update process comprises:
calculating three-dimensional neutron average flux, flow coupling factors and homogenization cross sections according to the radial flow, the axial flow and the radial flux;
and updating the three-dimensional neutron flux and the characteristic value based on the CMFD characteristic value iteration result, and calculating a leakage item.
The further optimization scheme is that when the two-dimensional equation is solved, linear source item calculation is firstly carried out, then linear source characteristic line scanning is carried out, and finally the fine network standard flux and the high-order standard flux distance are calculated.
The further optimization scheme is that the method for acquiring the characteristic value and the three-dimensional neutron flux of the reactor core further comprises the following steps:
and after obtaining the three-dimensional neutron flux and the characteristic value through CMFD iterative updating, judging whether the current three-dimensional neutron flux and the characteristic value are both converged, if so, outputting the current three-dimensional neutron flux and the characteristic value, and otherwise, recalculating the fission source and the scattering source.
The specific two-dimensional equation solving method comprises the following steps:
writing a two-dimensional equation into a characteristic linear equation;
deducing an analytical solution of the angular flux and the first-order integral angular flux according to a characteristic linear equation;
and extracting a high-order standard flux moment from the analytic solutions of the angular flux and the first-order integral angular flux.
The further optimization scheme is that a two-dimensional neutron transport equation with an introduced axial leakage term is solved by a linear source characteristic line method, and the equation can be written into a characteristic line formula equation, wherein the characteristic linear equation is expressed as:
wherein,is the angular flux in the i direction m of the flat source region, sigma t,i Is the total cross-section of the pipe,is the total source term of the order 0,is the 1 st order source term expansion moment, s i,m Is the characteristic line length, s i,m,k Is the length of the characteristic line segment k.
The further optimization scheme is that an analytic solution of the angular flux and the first-order integral angular flux is derived according to a characteristic linear equation, wherein the angular flux is expressed as:
the analytical solution for the first order integral angular flux is expressed as:
wherein,is the incident angle flux, F 1 (Σ t,i s)、F 2 (Σ t,i s)、G 1 (Σ t,i s k )、G 2 (Σ t,i s k ) Are all high order coefficients.
The further optimization scheme is that the higher-order scalar quantity moment is represented by a higher-order coefficient:
F 1 (Σ t,i s)=1-exp(-Σ t,i s)
F 2 (Σ t,i s)=2[Σ t,i s-F 1 (Σ t,i s)]-Σ t,i s i,m,k F 1 (Σ t,i s)
wherein s is k The length of the characteristic line segment k is indicated and s the characteristic line length.
The scheme provides a linear source acceleration-based three-dimensional neutron transport equation calculation system, which is applied to the linear source acceleration-based three-dimensional neutron transport equation calculation method and comprises the following steps: the device comprises a construction module, a conversion module and a calculation module;
the construction module is used for establishing a three-dimensional neutron transport equation;
the conversion module is used for converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
the calculation module is used for solving the one-dimensional equation and the two-dimensional equation respectively to obtain a reactor core characteristic value and three-dimensional neutron flux; the calculation module is also used for calculating the three-dimensional neutron flux and the high-order scalar flux moment based on the linear source form for flux expansion when solving the two-dimensional equation.
The present solution provides a non-transitory computer readable storage medium having stored thereon computer instructions for execution by a processor to perform the steps described for implementing a linear source acceleration based three-dimensional neutron transport equation calculation method.
Compared with the prior art, the invention has the following advantages and beneficial effects:
the three-dimensional neutron transport equation calculation method and system based on linear source acceleration provided by the invention are improved on the traditional plain form, flux expansion is carried out on the basis of the linear source form to calculate three-dimensional neutron flux and high-order standard flux moment, high-order expansion is carried out on the space through the linear source form, and the high-order standard flux moment is obtained at the same time, so that approximation introduced by the plain source method is eliminated, and the problem of low calculation efficiency caused by the fact that a dense plain source area grid is divided in the traditional two-dimensional one-dimensional method is solved; the method theoretically eliminates the flat source approximation, reduces the number of grids in the flat source region and improves the calculation efficiency.
Drawings
In order to more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the drawings that are required in the embodiments will be briefly described below, it should be understood that the following drawings only illustrate some embodiments of the present invention and therefore should not be considered as limiting the scope, and that those skilled in the art may also derive other related drawings based on these drawings without inventive effort. In the drawings:
FIG. 1 is a schematic flow chart of a method for calculating a three-dimensional neutron transport equation based on relaxation factors;
FIG. 2 is a schematic diagram of the two-dimensional/one-dimensional method basic principle;
FIG. 3 is a schematic diagram of a general flow of solving and calculating two-dimensional equations and one-dimensional equations;
fig. 4 is a schematic diagram of a two-dimensional equation scanning calculation solution flow.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
The three-dimensional neutron transport equation is directly solved by a one-step method, so that the calculation amount is large, the memory consumption is high, and the method is difficult to realize under the existing calculation condition. Therefore, a two-dimensional/one-dimensional method is provided, direct three-dimensional solution is converted into axial one-dimensional solution and radial two-dimensional solution respectively, and coupling is carried out through leakage items, so that the calculation requirement of directly solving a three-dimensional neutron transport equation by a one-step method is reduced; in the traditional two-dimensional one-dimensional coupling method, in order to obtain an accurate calculation result, a fine flat source area grid needs to be divided, so that the calculation efficiency is linearly reduced; the present invention provides the following embodiments to solve the above problems:
example 1
The embodiment provides a linear source acceleration-based three-dimensional neutron transport equation calculation method, as shown in fig. 1 and 2, including the steps of:
establishing a three-dimensional neutron transport equation;
converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
respectively solving the one-dimensional equation and the two-dimensional equation to obtain a reactor core characteristic value and three-dimensional neutron flux; and when the two-dimensional equation is solved, flux expansion is carried out on the basis of a linear source form to calculate three-dimensional neutron flux and high-order standard flux moment.
The three-dimensional neutron transport equation is:
wherein m represents an angle, g represents an energy group,. phi g,m (x, y, z) represents angular flux, x, y, z represent x, y, z coordinates of the location in space, ξ, respectively m Representing the cosine of the angle of incidence of the azimuth and the x-axis, sigma t,g (r) represents the total cross-section, η represents the amplitude sine, μ represents the polar cosine, Q g (x, y, z) represents the total source term.
The one-dimensional equation and the two-dimensional equation are obtained by the following method:
and (3) integrating the radial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a one-dimensional equation:
and (3) integrating the axial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a two-dimensional equation:
wherein psi g,m,i,j (z) angular flux, Q, of the z-th layer representing the radial (i, j) position angle m energy group g g,i,j (z) a one-dimensional total source term, Q, representing the radial (i, j) position g (x, y) represents a two-dimensional total source term,
indicating a radial leakage term, # g,m (x, y) denotes radial angular flux, Σ t,g,i,j (z) represents a one-dimensional total cross section, Sigma t,g (x, y) represents a two-dimensional total cross section,an axial leakage term is indicated.
The method for acquiring the characteristic value and the three-dimensional neutron flux of the reactor core comprises the following steps:
calculating fission sources and scattering sources;
on the basis of a fission source and a scattering source, respectively solving a one-dimensional equation and a two-dimensional equation to obtain radial flow, axial flow and radial flux;
performing CMFD iterative updating based on the radial flow, the axial flow and the radial flux to obtain a reactor core characteristic value and a three-dimensional neutron flux;
the CMFD iterative update process comprises:
calculating three-dimensional neutron average flux, flow coupling factors and homogenization cross sections according to the radial flow, the axial flow and the radial flux;
and updating the three-dimensional neutron flux and the characteristic value based on the CMFD characteristic value iteration result, and calculating a leakage item.
As shown in fig. 4, the two-dimensional equation solving method includes the steps of: and (3) performing linear source item calculation, performing linear source characteristic line scanning, and finally calculating the fine network standard flux and the high-order standard flux distance.
Specifically, a two-dimensional equation is written into a characteristic linear equation;
deducing an analytical solution of the angular flux and the first-order integral angular flux according to a characteristic linear equation;
and extracting a high-order standard flux moment from the analytic solutions of the angular flux and the first-order integral angular flux.
The characteristic linear equation is expressed as:
wherein,is the angular flux in the i direction m of the flat source region, sigma t,i Is the total cross-section of the pipe,is the total source term of the order 0,is the 1 st order source term expansion moment, s i,m Is the characteristic line length, s i,m,k Is the length of the characteristic line segment k.
The angular flux is expressed as:
the analytical solution for the first order integral angular flux is expressed as:
wherein,is the incident angle flux, F 1 (Σ t,i s)、F 2 (Σ t,i s)、G 1 (Σ t,i s k )、G 2 (Σ t,i s k ) Are all high order coefficients.
The higher order scalar flux moment is represented by a higher order coefficient:
F 1 (Σ t,i s)=1-exp(-Σ t,i s)
F 2 (Σ t,i s)=2[Σ t,i s-F 1 (Σ t,i s)]-Σ t,i s i,m,k F 1 (Σ t,i s)
wherein s is k The length of the characteristic line segment k is indicated and s the characteristic line length.
Example 2
The embodiment provides a linear source acceleration-based three-dimensional neutron transport equation calculation system, which is applied to the method in the previous embodiment, and is characterized by comprising the following steps: the device comprises a construction module, a conversion module and a calculation module;
the construction module is used for establishing a three-dimensional neutron transport equation;
the conversion module is used for converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
the calculation module is used for solving the one-dimensional equation and the two-dimensional equation respectively to obtain a reactor core characteristic value and three-dimensional neutron flux; the calculation module is also used for calculating the three-dimensional neutron flux and the high-order scalar flux moment based on the linear source form for flux expansion when solving the two-dimensional equation.
Example 3
The present embodiments provide a non-transitory computer readable storage medium having stored thereon computer instructions that, when executed by a processor, implement the steps of the method of embodiment 1.
The processor performs the steps further comprising:
(1) before iterative computation, geometric modeling is carried out, a characteristic line is generated, and linear source related parameters are computed;
(2) calculating fission source and scattering source items;
(3) and respectively carrying out one-dimensional equation solution and two-dimensional equation solution to obtain axial flow, radial flow and radial flux. The method comprises the following steps of solving a two-dimensional equation, wherein the two-dimensional equation is based on a linear source method, namely linear source item calculation, characteristic line scanning based on the linear source method, and standard flux and high-order standard flux moment calculation are carried out;
(4) carrying out coarse net effective differential Calculation (CMFD) according to two-dimensional and one-dimensional calculation results;
(5) and (5) judging whether the characteristic value and the crack rate are converged, finishing the calculation if the characteristic value and the crack rate are converged, and returning to the step (2) if the characteristic value and the crack rate are not converged.
In the step (3), the two-dimensional calculation by using the linear source method is the key of the embodiment. The method is improved on the basis of the traditional flat source approximation method, the high-order scalar flux moment is additionally calculated in the processes of source item calculation, characteristic line scanning, flux calculation and the like, high-order expansion is carried out in space, approximation introduced by the flat source method is eliminated theoretically, and therefore accurate calculation results can be obtained by adopting fewer flat source area grids, and calculation efficiency is improved.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are merely exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, 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. The linear source acceleration-based three-dimensional neutron transport equation calculation method is characterized by comprising the following steps of:
establishing a three-dimensional neutron transport equation;
converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
respectively solving the one-dimensional equation and the two-dimensional equation to obtain a reactor core characteristic value and three-dimensional neutron flux; and when the two-dimensional equation is solved, flux expansion is carried out on the basis of a linear source form to calculate three-dimensional neutron flux and high-order standard flux moment.
2. The linear source acceleration based three-dimensional neutron transport equation calculation method of claim 1, wherein the three-dimensional neutron transport equation is:
wherein m represents an angle, g represents an energy group,. phi g,m (x, y, z) represents angular flux, x, y, z represent x, y, z coordinates of the location in space, ξ, respectively m Representing the cosine of the angle of incidence of the azimuth and the x-axis, sigma t,g (r) represents the total cross-section, η represents the amplitude sine, μ represents the polar cosine, Q g (x, y, z) represents the total source term.
3. The linear source acceleration-based three-dimensional neutron transport equation calculation method according to claim 2, wherein the one-dimensional equation and the two-dimensional equation are obtained by:
and (3) integrating the radial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a one-dimensional equation:
and (3) integrating the axial direction in the area of each rod of each layer by using a three-dimensional neutron transport equation to obtain a two-dimensional equation:
wherein psi g,m,i,j (z) angular flux, Q, of the z-th layer representing the radial (i, j) position angle m energy group g g,i,j (z) a one-dimensional total source term, Q, representing the radial (i, j) position g (x, y) represents a two-dimensional total source term,
4. The linear source acceleration-based three-dimensional neutron transport equation calculation method according to claim 1, wherein the reactor core characteristic value and three-dimensional neutron flux acquisition method comprises the steps of:
calculating fission sources and scattering sources;
on the basis of a fission source and a scattering source, respectively solving a one-dimensional equation and a two-dimensional equation to obtain radial flow, axial flow and radial flux;
performing CMFD iterative updating based on the radial flow, the axial flow and the radial flux to obtain a reactor core characteristic value and a three-dimensional neutron flux;
the CMFD iterative update process comprises:
calculating three-dimensional neutron average flux, flow coupling factors and homogenization cross sections according to the radial flow, the axial flow and the radial flux;
and updating the three-dimensional neutron flux and the characteristic value based on the CMFD characteristic value iteration result, and calculating a leakage item.
5. The linear source acceleration-based three-dimensional neutron transport equation calculation method according to claim 3, wherein the two-dimensional equation solving method comprises the steps of:
writing a two-dimensional equation into a characteristic linear equation;
deducing an analytical solution of the angular flux and the first-order integral angular flux according to a characteristic linear equation;
and extracting a high-order standard flux moment from the analytic solutions of the angular flux and the first-order integral angular flux.
6. The linear source acceleration-based three-dimensional neutron transport equation calculation method of claim 5, wherein the characteristic linear equation is expressed as:
wherein,is the angular flux in the i direction m of the flat source region, sigma t,i Is the total cross-section of the pipe,is the total source term of the order 0,is the unfolding moment of the source term of order 1, s i,m Is the characteristic line length, s i,m,k Is the length of the characteristic line segment k.
7. The linear source acceleration-based three-dimensional neutron transport equation calculation method of claim 6, wherein the angular flux is expressed as:
the analytical solution for the first order integral angular flux is expressed as:
8. The linear-source-acceleration-based three-dimensional neutron transport equation calculation method of claim 7, wherein the higher-order scalar moments are represented by higher-order coefficients:
F 1 (Σ t,i s)=1-exp(-Σ t,i s)
F 2 (Σ t,i s)=2[Σ t,i s-F 1 (Σ t,i s)]-Σ t,i s i,m,k F 1 (Σ t,i s)
wherein s is k The length of the characteristic line segment k is indicated and s the characteristic line length.
9. The linear source acceleration-based three-dimensional neutron transport equation calculation system is applied to the method of any one of claims 1 to 8, and is characterized by comprising the following steps of: the device comprises a construction module, a conversion module and a calculation module;
the construction module is used for establishing a three-dimensional neutron transport equation;
the conversion module is used for converting a three-dimensional neutron transport equation into a one-dimensional equation and a two-dimensional equation;
the calculation module is used for solving the one-dimensional equation and the two-dimensional equation respectively to obtain a reactor core characteristic value and three-dimensional neutron flux; the calculation module is also used for calculating the three-dimensional neutron flux and the high-order scalar flux moment based on the linear source form for flux expansion when solving the two-dimensional equation.
10. A non-transitory computer readable storage medium having stored thereon computer instructions, characterized in that the instructions, when executed by a processor, implement the steps of the method of any one of claims 1-8.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210815861.XA CN115048810A (en) | 2022-07-12 | 2022-07-12 | Three-dimensional neutron transport equation calculation method and system based on linear source acceleration |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210815861.XA CN115048810A (en) | 2022-07-12 | 2022-07-12 | Three-dimensional neutron transport equation calculation method and system based on linear source acceleration |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115048810A true CN115048810A (en) | 2022-09-13 |
Family
ID=83165510
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210815861.XA Pending CN115048810A (en) | 2022-07-12 | 2022-07-12 | Three-dimensional neutron transport equation calculation method and system based on linear source acceleration |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115048810A (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112632822A (en) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | Nuclear reactor neutron flux obtaining method and device based on three-dimensional leakage item segmentation |
CN113672849A (en) * | 2021-08-26 | 2021-11-19 | 中国核动力研究设计院 | One-step transportation calculation method and system based on axial flux expansion |
-
2022
- 2022-07-12 CN CN202210815861.XA patent/CN115048810A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112632822A (en) * | 2020-12-21 | 2021-04-09 | 中国核动力研究设计院 | Nuclear reactor neutron flux obtaining method and device based on three-dimensional leakage item segmentation |
CN113672849A (en) * | 2021-08-26 | 2021-11-19 | 中国核动力研究设计院 | One-step transportation calculation method and system based on axial flux expansion |
Non-Patent Citations (2)
Title |
---|
CHEN ZHAO, XINGJIE PENG, HONGBO ZHANG, WENBO ZHAO, ZHANG CHEN, QING LI,: "A linear source scheme for the 2D/1D transport method in SHARK", ANNALS OF NUCLEAR ENERGY, vol. 161, pages 1 - 2 * |
张宏博,赵晨,彭星杰等: "数字化反应堆高保真中子学程序SHARK研发", 原子能科学技术, vol. 56, no. 02, pages 336 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103106301B (en) | A kind of method of the calculating radiation shield be coupled with characteristic line method based on Monte Carlo method | |
CN107273582B (en) | Calculation method for neutron transport burnup coupling analysis of fast neutron reactor | |
Smith et al. | Challenges in the development of high-fidelity LWR core neutronics tools | |
CN112906271B (en) | Reactor transient physical thermal full-coupling fine numerical simulation method and system | |
KR20060015760A (en) | Computer modelling method for the core of a nuclear reactor | |
Xu et al. | Development of a new parallel SN code for neutron-photon transport calculation in 3-D cylindrical geometry | |
Du et al. | Development of an assistant program for CAD-to-cosRMC modelling | |
CN112906272B (en) | Reactor steady-state physical thermal full-coupling fine numerical simulation method and system | |
CN103116667B (en) | Method of obtaining fusion reactor experimental covering module neutronics parameters | |
CN113672849A (en) | One-step transportation calculation method and system based on axial flux expansion | |
CN111814343B (en) | Reactor core power distribution online reconstruction method for comprehensive in-reactor and out-reactor detector measurement values | |
CN111651896A (en) | Wind power plant flow field calculation method based on actual wind speed and actual complex terrain | |
CN112632822B (en) | Nuclear reactor neutron flux obtaining method and device based on three-dimensional leakage item segmentation | |
CN106202865A (en) | A kind of calculate the method for arbitrary order coefficient in neutron transport discrete locking nub method | |
CN105512387A (en) | Fine power reconstruction method for fuel assembly | |
CN115048810A (en) | Three-dimensional neutron transport equation calculation method and system based on linear source acceleration | |
CN112507445A (en) | Geometrical structure reconstruction method of topological optimization result | |
CN113704996A (en) | Quasi-three-dimensional transport calculation method and system based on axial flux expansion | |
CN108694299A (en) | Two dimensional finite element neutronics stable state computational methods based on ICEM-CFD | |
Fiorina et al. | Creation of an OpenFOAM fuel performance class based on FRED and integration into the GeN-foam multi-physics code | |
CN114547952A (en) | Pressurized water reactor neutron transport method based on deep learning-interface flow coupling | |
CN114091135A (en) | Lisp language-based ship body curved surface outer plate weight gravity center calculation program | |
Kang et al. | The advanced multilevel predictor-corrector quasi-static method for pin-resolved neutron kinetics simulation | |
CN110020501B (en) | Equipment structure improvement method and equipment based on multiphase computational fluid mechanics interpolation | |
CN112380744A (en) | Method for accelerated solution of generalized conjugate neutron transport equation |
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 |