CN113887027A - Multi-scale coupling method and device for layered burnup of dispersed fuel and toxic particles - Google Patents
Multi-scale coupling method and device for layered burnup of dispersed fuel and toxic particles Download PDFInfo
- Publication number
- CN113887027A CN113887027A CN202111092461.2A CN202111092461A CN113887027A CN 113887027 A CN113887027 A CN 113887027A CN 202111092461 A CN202111092461 A CN 202111092461A CN 113887027 A CN113887027 A CN 113887027A
- Authority
- CN
- China
- Prior art keywords
- model
- fuel
- microscopic
- effective
- poison
- 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
Links
- 239000002245 particle Substances 0.000 title claims abstract description 108
- 239000000446 fuel Substances 0.000 title claims abstract description 104
- 238000010168 coupling process Methods 0.000 title claims abstract description 23
- 231100000331 toxic Toxicity 0.000 title abstract description 11
- 230000002588 toxic effect Effects 0.000 title abstract description 11
- 239000002574 poison Substances 0.000 claims abstract description 79
- 231100000614 poison Toxicity 0.000 claims abstract description 78
- 238000000034 method Methods 0.000 claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 44
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 230000004907 flux Effects 0.000 claims abstract description 20
- 238000002156 mixing Methods 0.000 claims abstract description 14
- 230000008878 coupling Effects 0.000 claims abstract description 10
- 238000005859 coupling reaction Methods 0.000 claims abstract description 10
- 230000009257 reactivity Effects 0.000 claims description 13
- 238000005253 cladding Methods 0.000 claims description 11
- 239000000463 material Substances 0.000 claims description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000002485 combustion reaction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 5
- 239000010410 layer Substances 0.000 description 33
- 230000008569 process Effects 0.000 description 9
- 239000000919 ceramic Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- CMIHHWBVHJVIGI-UHFFFAOYSA-N gadolinium(III) oxide Inorganic materials [O-2].[O-2].[O-2].[Gd+3].[Gd+3] CMIHHWBVHJVIGI-UHFFFAOYSA-N 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 230000035755 proliferation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 239000011248 coating agent Substances 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 239000006185 dispersion Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000265 homogenisation Methods 0.000 description 2
- 238000009462 micro packaging Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000011247 coating layer Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005538 encapsulation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000008187 granular material Substances 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000012088 reference solution Substances 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
The invention discloses a multi-scale coupling method and a multi-scale coupling device for layered burnup of dispersed fuel and toxic particles, wherein the method comprises the following steps: constructing a microscopic fine model according to the geometric parameters of the rod beam grid cell model; performing neutron transport calculation on the microscopic fine model to obtain an effective value-added factor of the microscopic fine model; mixing the poison particle balls with the SiC matrix according to the effective value-added factors to obtain a microscopic equivalent uniform model and the effective share of the poison particle balls; taking the rod bundle grid cell model as a macroscopic model, and mixing the poison particle balls in the fuel rods and the SiC matrix according to effective shares to obtain effective multiplication factors and average flux of the fuel rods; and performing fuel consumption calculation according to the effective multiplication factors and the average flux, and counting N, n, phi and other important fuel consumption data to realize multi-scale coupled fuel consumption calculation. According to the invention, the macroscopic rod bundle grid cell model is subjected to mixing processing according to the microscopic equivalent homogeneous model, so that the calculation efficiency is improved, and the method can be widely applied to the field of nuclear engineering.
Description
Technical Field
The invention relates to the field of nuclear engineering, in particular to a multi-scale coupling method and a multi-scale coupling device for layered burnup of dispersed fuel and toxic particles.
Background
In the safe operation process of the reactor, the control of the reactivity is particularly important; the setting of the burnable poison plays an important role in controlling the reactivity, is beneficial to reducing the initial reactivity of the reactor at the beginning of the service life, and plays an important role in unattended control of the reactor.
The full ceramic micro-packaging dispersion fuel is a novel accident-resistant fuel, and the coating fuel particles are dispersed in a matrix, so that the fuel has the advantages of higher burning depth, better containment and heat transfer performance and the like. Wherein, the dispersed burnable poison particles are added into the FCM fuel, and the reactivity process can be flexibly controlled for a long time under the condition of not causing power distortion. Furthermore, from a production process point of view, FCM is well compatible with particulate burnable poison.
Different from fuel particles, the burnable poison particles have extremely large neutron absorption cross sections, and the spatial self-shielding effect caused by the special spatial structure of the burnable poison particles causes a phenomenon of layered combustion in the burnup process, namely: the outer layer of the granule is preferentially consumed while the inner layer maintains higher integrity; the fuel particles do not experience such significant stratification. In the existing method, to accurately reduce the layering phenomenon of the toxic particles, the structure of the toxic particles needs to be further subdivided, so that huge grid calculation amount is caused, and the calculation load is large, which is not realizable in large-scale full-stack calculation; the fuel particles dispersed in the matrix are different from the poison particles, so that the layering phenomenon does not exist, the fuel particles do not need to be divided too finely, but have double heterogeneity and need to be treated separately.
The matrix is mixed with various complex structures such as cladding fuel particles, burnable poison particles and the like and corresponding characteristics thereof, so that an efficient and high-fidelity solving method is absent at present in the process of processing the layered burnup effect of the burnable poison particles in the all-ceramic micro-encapsulation dispersed Fuel (FCM).
Disclosure of Invention
In order to solve at least one of the technical problems in the prior art to a certain extent, the invention aims to provide a multi-scale coupling method and a multi-scale coupling device for layered burnup of dispersed fuel and toxic particles, which are used for solving the layered burnup phenomenon of the dispersed combustible toxic particles in the fully ceramic micro-encapsulated dispersed fuel (FCM fuel).
The technical scheme adopted by the invention is as follows:
a multi-scale coupling method of stratified burnup of dispersed fuel and poison particles, comprising the steps of:
s1, constructing a microscopic fine model according to the geometric parameters of the rod beam grid cell model; wherein the microscopic fine model comprises poison particle balls, fuel particles and a SiC matrix;
s2, performing neutron transport calculation on the microscopic fine model to obtain an effective value-added factor of the microscopic fine model;
s3, mixing the poison particle balls with the SiC matrix according to the effective value-added factors to obtain a microscopic equivalent uniform model with the same reactivity as the microscopic fine model and the effective share of the poison particle balls;
s4, taking the rod bundle grid cell model as a macro model, and mixing the poison particle balls and the SiC matrix in the fuel rod according to the effective share to obtain the effective multiplication factor and the average flux of the fuel rod;
s5, performing fuel consumption calculation according to the effective multiplication factors and the average flux to obtain the new nuclear density N of the fuel particles;
s6, in a microscopic fine model, performing burnup calculation according to the average flux to obtain the new nucleus density n of the poison particle ball;
s7, returning the new nucleus density N and the new nucleus density N to the step S1, and executing the steps S1-S6 to realize the multi-scale coupled burnup calculation.
Further, the microscopic fine model also comprises a coating layer and a moderator layer;
the sizes of the poison particle balls, the fuel particles, the SiC matrix, the cladding layer and the moderator layer are all set according to the volume proportion of each part of the fuel rod bundle grid element in the macro, and the material properties of each ball layer are consistent with those of the macro grid element.
Further, in the microscopic fine model, the poison particle spheres appear as a single geometric structure without double heterogeneity treatment; the fuel particles are dispersed in the SiC matrix and need to be subjected to double heterogeneous treatment;
the neutron transport calculation is carried out on the microscopic fine model to obtain an effective value-added factor, and the method comprises the following steps:
and performing neutron transport calculation on the microscopic fine model with double heterogeneity by adopting a collision probability method to obtain effective value-added factors of the microscopic fine model.
Further, the reactivity is the same, which means that the effective proliferation factors corresponding to the microscopic fine model and the microscopic equivalent homogeneous model are equal.
Further, the effective proliferation factors corresponding to the microscopic fine model and the microscopic equivalent homogeneous model are equalized by the following steps:
uniformly mixing the poison particle balls into the SiC matrix, and adjusting the nuclear density of the SiC matrix to keep the total amount of the SiC matrix unchanged;
adjusting the effective average nuclear density of the poison particle balls to enable the effective multiplication factor of the microscopic equivalent homogeneous model to be equal to that of the microscopic fine model;
wherein the effective average nucleus density obtained by the adjustment is used as the effective share of the poison particle balls.
Further, the obtaining of the effective multiplication factor and the average flux of the fuel rod comprises:
and calling a grid cell program for neutron physical analysis, and performing neutron transport calculation on the whole rod bundle grid cell by adopting a collision probability method to obtain the effective multiplication factor and the average flux of the fuel rod.
Further, the poison particle balls are divided into 10 layers.
Further, in the microscopic equivalent homogeneous model, the clad layer and the moderator layer remain unchanged.
Furthermore, light water is adopted in the pressurized water reactor of the moderator layer.
The other technical scheme adopted by the invention is as follows:
a multi-scale coupling device for stratified burnup of dispersed fuel and poison particles, comprising:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method described above.
The invention has the beneficial effects that: according to the invention, the macroscopic rod bundle grid cell model is subjected to mixing processing according to the microscopic equivalent homogeneous model, so that the technical problems of over-dense grids, over-large calculated amount and the like are avoided in the fuel consumption solving process, and the calculation efficiency is greatly improved while the fuel consumption performance of fuel particles and poison particles can be truly and comprehensively reflected.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following description is made on the drawings of the embodiments of the present invention or the related technical solutions in the prior art, and it should be understood that the drawings in the following description are only for convenience and clarity of describing some embodiments in the technical solutions of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of a multi-scale coupling method for stratified burnup of dispersed fuel and poison particles in an embodiment of the present disclosure;
FIG. 2 is a schematic representation of an FCM fuel and dispersed burnable poison particles in an embodiment of the present invention;
FIG. 3 is a schematic illustration of a microscopic fine model in an embodiment of the present invention;
FIG. 4 is a schematic illustration of a microscopic equivalent uniformity model in an embodiment of the present invention;
FIG. 5 is an error map of a multi-scale coupling algorithm in an embodiment of the invention.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like or similar reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are illustrative only for the purpose of explaining the present invention, and are not to be construed as limiting the present invention. The step numbers in the following embodiments are provided only for convenience of illustration, the order between the steps is not limited at all, and the execution order of each step in the embodiments can be adapted according to the understanding of those skilled in the art.
In the description of the present invention, it should be understood that the orientation or positional relationship referred to in the description of the orientation, such as the upper, lower, front, rear, left, right, etc., is based on the orientation or positional relationship shown in the drawings, and is only for convenience of description and simplification of description, and does not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention.
In the description of the present invention, the meaning of a plurality of means is one or more, the meaning of a plurality of means is two or more, and larger, smaller, larger, etc. are understood as excluding the number, and larger, smaller, inner, etc. are understood as including the number. If the first and second are described for the purpose of distinguishing technical features, they are not to be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
In the description of the present invention, unless otherwise explicitly limited, terms such as arrangement, installation, connection and the like should be understood in a broad sense, and those skilled in the art can reasonably determine the specific meanings of the above terms in the present invention in combination with the specific contents of the technical solutions.
As shown in fig. 1, the present embodiment provides a multi-scale coupling method for stratified burnup of dispersed fuel and poison particles, comprising the steps of:
s101, constructing a microscopic fine model according to the geometric parameters of the rod beam grid cell model; wherein the microscopic fine model comprises poison particle balls, fuel particles and a SiC matrix.
As shown in fig. 2, a microscopic fine sphere layer model was created based on the corresponding geometric parameters of the full ceramic microencapsulated dispersed fuel and burnable poison particle rod bundle cell model, as shown in fig. 3. The innermost part is a burnable poison particle ball which is finely divided into 10 layers; the space of the outer layer is filled with SiC matrix dispersed with coating fuel particles; the secondary outer layer is a cladding layer, and the outermost layer is a moderator layer (light water is adopted in a common pressurized water reactor). The size of each ball layer is set according to the volume proportion of each part of the macroscopic fuel bundle grid element, and the material property of each ball layer is consistent with that of the macroscopic grid element.
S102, performing neutron transport calculation on the microscopic fine model to obtain an effective value-added factor of the microscopic fine model.
In the microscopic model, a poison particle sphere appears as a single geometric structure, is not included in the dispersion-matrix, and does not need double heterogeneous treatment; and the fuel particles are dispersed in the SiC matrix, and then double heterogeneous treatment is needed.
The existing collision probability method can efficiently solve the problem of double heterogeneity of fuel particles and uniform SiC matrixes, a neutron physical analysis grid element program DRAGON is called, neutron transport calculation is carried out on the microscopic fine model with the double heterogeneity by adopting the collision probability method, and an effective value-added factor K-ref of the microscopic fine model is obtained.
S103, mixing the poison particle balls with the SiC matrix according to the effective value-added factors to obtain a microscopic equivalent uniform model with the same reactivity as the microscopic fine model and the effective share of the poison particle balls.
And equivalently homogenizing the microscopic fine model by adopting a reactive equivalence strategy according to the microscopic fine model to obtain a microscopic equivalent uniform model, as shown in figure 4.
In an equivalent model, an outer moderator layer and a cladding layer are kept unchanged, the spherical structure of the poison particles is scattered and uniformly mixed into the SiC matrix, and the structure and the components of the dispersed fuel particles in the matrix are kept unchanged. Among them, it is necessary to adjust a series of nuclear density changes due to volume changes of the structure according to the volume ratio.
The strategy of equivalence is to adjust the effective average nuclear density (effective share is fv) of the poison particle ball, so that the new uniform ball model and the original fine ball model have the same reactivity; namely, the effective multiplication factor K of the microscopic equivalent homogeneous model is ensured to be equal to the effective multiplication factor K-ref of the microscopic fine model. When neutron transport calculation is carried out, double heterogeneity processing is still carried out on dispersed fuel particles by using a collision probability method.
S104, taking the rod bundle grid cell model as a macroscopic model, and mixing the poison particle balls in the fuel rod and the SiC matrix according to the effective share to obtain the effective multiplication factor and the average flux of the fuel rod.
And S105, performing burnup calculation according to the effective multiplication factor and the average flux to obtain the new nuclear density N of the fuel particles.
In this example, the FCM fuel and the rod bundle cell model embedded with dispersed burnable poison particles are used as a macroscopic model. The effective share fv of the burnable poison obtained in the micro homogeneous model is used as the average effective share of the burnable poison mixed in the macro model.
In a macroscopic model, a SiC matrix in the fuel rod and dispersed poison particle balls are uniformly mixed, the nuclear density of the matrix is changed according to the change of the volume proportion, and the new nuclear density of the poison is determined according to the effective nuclear fraction fv of the combustible poison.
Neutron transport calculation is carried out on the whole rod bundle grid element by using a neutron physical analysis grid element program DRAGON and adopting a collision probability method (wherein double heterogeneity treatment is carried out on dispersed fuel particles), and an integral effective multiplication factor K and an average flux phi are obtained; and performing burnup calculation on the whole rod bundle grid cell to obtain the new nuclear density N of the fuel part.
And S106, in the microscopic fine model, performing burnup calculation according to the average flux to obtain the new nucleus density n of the poison particle ball.
And finally, in a microscopic fine model, performing one-time layered burnup calculation under the level of the average flux phi to obtain the multilayer new nuclear density n of the toxic particles.
And S107, returning the new kernel density N and the new kernel density N to the step S101, and executing the steps S101-S106 to realize the multi-scale coupled burn-up calculation.
And (4) updating the poison layer nucleus density N and the fuel layer nucleus density N of the micro model again, and circularly repeating the steps S101-S107 to finish the multi-scale coupled fuel consumption calculation in the whole service life.
The method of the present embodiment is described in further detail below with reference to the following embodiments and the accompanying drawings, but the embodiment of the present embodiment is not limited thereto.
This example will calculate a fuel bundle cell model of a fully ceramic microencapsulated dispersed Fuel (FCM) loaded with dispersed burnable poison particles; wherein, the moderator grid element is a square with the side length of 1.6cm, the outer radius of the fuel cladding is 0.55cm, and the radius of the whole fuel rod bundle is 0.5 cm. In the fuel rod region, fuel particles and poison particles are dispersed in the SiC matrix; wherein the fuel particles account for 30% of the total fuel bundle volume fraction and the particulate poisons account for 1% of the total fuel bundle volume fraction. The fuel adopts UO2 with the enrichment degree of 20%, the radius of fuel particles is 800um, and the outer diameter of the fuel particles is 980 um; the burnable poison is Gd2O3 with the particle radius of 500 um. The technical scheme is as follows:
s201, establishing a microscopic fine model.
And establishing a microscopic fine spherical layer model according to the corresponding geometric parameters of the macroscopic rod bundle grid cell model of the all-ceramic micro-packaging dispersion fuel and the burnable poison particles. The size of each ball layer is set according to the volume proportion of each part of the macroscopic fuel bundle grid element, and the material property of each ball layer is consistent with that of the macroscopic grid element.
In this case, the burnable poison particle sphere with a diameter of 500um was divided into 10 parts by equal volume, the outer layer was a mixed layer of fuel particles and a SiC matrix with a radius of 0.2321cm, the cladding was made of Zr-4 alloy, the outer diameter of the cladding was 0.2473cm, and the outer diameter of the moderator layer (water) was 0.3435 cm.
And S202, calculating neutron transport of the microscopic fine model.
Calling a neutron physical analysis cell program DRAGON, and performing neutron transport calculation on the microscopic fine model with double heterogeneity by adopting a collision probability method to obtain an effective incremental factor K-ref of 1.086904.
And S203, microscopic equivalent homogenization.
And (5) according to the microscopic fine model in the S201, equivalently homogenizing the microscopic fine model by adopting a reactive equivalence strategy to obtain a microscopic equivalent uniform model. In the equivalent model, the external moderator layer and the cladding layer are kept unchanged, the spherical structure of the poison particles is scattered and uniformly mixed into the SiC matrix, and the structure and the components of the dispersed fuel particles in the matrix are kept unchanged. Among them, it is necessary to adjust a series of nuclear density changes due to volume changes of the structure according to the volume ratio.
Adjusting the nuclear density of the SiC matrix according to the mixing proportion; the volume fraction of burnable poison is 1% here, so that the space occupied by the SiC molecules is expanded from 69% to 70%, and therefore the nucleus density needs to be multiplied by 69/70 to keep the total amount constant. While the volume fraction of 30% occupied by the fuel particles remains unchanged.
Combining a collision probability method, performing neutron transport calculation on the microstructure containing the dispersed fuel particles and the new SiC-poison matrix, and obtaining the reactivity of the microstructure; changing the reactivity of the model by continuously adjusting the effective average nuclear density (effective share is fv) of the poison material, so that the new uniform sphere model and the original fine sphere model have the same reactivity; that is, the effective proliferation factor (K) is ensured to be equal. The effective fraction fv is recorded.
In the initial state, the effective fraction fv is 0.0106 and the corresponding K is 1.086698.
And S204, calculating micro-macro coupling.
In the macroscopic model, the SiC matrix and the dispersed poison particle balls in the fuel rod are uniformly mixed, and the nuclear density of the matrix is changed according to the change of the volume proportion:
ND(SiC)＝ND(SiC)*69/70；
ND(Gd2O3)＝ND(Gd2O3)*1/70；
where ND denotes the nuclear density (number-density) of the corresponding nuclide.
Introducing the fv factor in S203 as the average effective share of the combustible poison in the macroscopic model; and multiplying the nucleus density of the diluted poison by the effective share to obtain the effective nucleus density of the poison under uniform mixing.
ND-eff(Gd2O3)＝ND(Gd2O3)*fv；
Calling a neutron physical analysis grid cell program DRAGON again, and performing neutron transport calculation on the whole rod bundle grid cell (wherein double heterogeneity processing is performed on dispersed fuel particles) by adopting a collision probability method to obtain an integral effective multiplication factor K and an average flux phi; and performing burnup calculation on the whole rod bundle grid cell to obtain the new nuclear density N of the fuel part.
In the initial state, K — 1.095868; the average flux level after normalization, phi, is 3.4065721E + 00.
And S205, macro-micro coupling calculation.
In the microscopic fine model, a stratified burnup calculation was made at a constant flux φ (obtained in S4) to obtain a multilayered new nucleus density n of poison particles. And (4) updating the poison layer nucleus density N and the fuel layer nucleus density N of the micro model again, and repeating the steps S201-S205 in a circulating manner, so that the multi-scale coupled fuel consumption calculation in the whole service life can be completed. And f, K, phi, N, n and other important fuel consumption data in the whole calculation process are counted.
S206, verifying and concluding.
The embodiment is based on a multi-scale coupling method, and has the inherent advantages of high efficiency, small calculation amount and the like in more complicated and large-scale core calculation. And solving the whole fuel consumption process by completely adopting a collision probability method to serve as a reference solution of the calculation precision. By comparison, in the whole layered burnup process of calculating and solving dispersed fuel and toxic particles, the overall error level is within 250pcm, the initial error level is about 100pcm, and the final error is higher due to a larger error caused by too low nuclide density at the end of the service life and an accumulated error along with burnup to a certain extent, but is kept within 300 pcm. Therefore, the method is proved to have higher precision, and the error graph is shown in figure 5.
In summary, compared with the prior art, the method of the embodiment has the following beneficial effects: the method combines the equivalent homogenization and the collision probability method, avoids the technical problems of over-dense grids, over-large calculated amount and the like in the fuel consumption solving process, and greatly improves the calculation efficiency while truly and comprehensively reflecting the fuel consumption performance of the fuel particles and the toxic particles. Wherein, under larger scale computing conditions (such as component level and full stack level computing), more excellent computing performance can be obtained.
The present embodiments also provide a multi-scale coupling device for stratified burnup of dispersed fuel and poison particles, comprising:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method described above.
The multi-scale coupling device for the layered burnup of dispersed fuel and poison particles can execute the multi-scale coupling method for the layered burnup of dispersed fuel and poison particles, which is provided by the method embodiment of the invention, can execute any combination implementation steps of the method embodiment, and has corresponding functions and beneficial effects of the method.
In alternative embodiments, the functions/acts noted in the block diagrams may occur out of the order noted in the operational illustrations. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality/acts involved. Furthermore, the embodiments presented and described in the flow charts of the present invention are provided by way of example in order to provide a more thorough understanding of the technology. The disclosed methods are not limited to the operations and logic flows presented herein. Alternative embodiments are contemplated in which the order of various operations is changed and in which sub-operations described as part of larger operations are performed independently.
Furthermore, although the present invention is described in the context of functional modules, it should be understood that, unless otherwise stated to the contrary, one or more of the described functions and/or features may be integrated in a single physical device and/or software module, or one or more functions and/or features may be implemented in a separate physical device or software module. It will also be appreciated that a detailed discussion of the actual implementation of each module is not necessary for an understanding of the present invention. Rather, the actual implementation of the various functional modules in the apparatus disclosed herein will be understood within the ordinary skill of an engineer, given the nature, function, and internal relationship of the modules. Accordingly, those skilled in the art can, using ordinary skill, practice the invention as set forth in the claims without undue experimentation. It is also to be understood that the specific concepts disclosed are merely illustrative of and not intended to limit the scope of the invention, which is defined by the appended claims and their full scope of equivalents.
In the foregoing description of the specification, reference to the description of "one embodiment/example," "another embodiment/example," or "certain embodiments/examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present invention have been shown and described, it will be understood by those of ordinary skill in the art that: various changes, modifications, substitutions and alterations can be made to the embodiments without departing from the principles and spirit of the invention, the scope of which is defined by the claims and their equivalents.
While the preferred embodiments of the present invention have been illustrated and described, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims.
Claims (10)
1. A multi-scale coupling method for stratified burnup of dispersed fuel and poison particles, comprising the steps of:
s1, constructing a microscopic fine model according to the geometric parameters of the rod beam grid cell model; wherein the microscopic fine model comprises poison particle balls, fuel particles and a SiC matrix;
s2, performing neutron transport calculation on the microscopic fine model to obtain an effective value-added factor of the microscopic fine model;
s3, mixing the poison particle balls with the SiC matrix according to the effective value-added factors to obtain a microscopic equivalent uniform model with the same reactivity as the microscopic fine model and the effective share of the poison particle balls;
s4, taking the rod bundle grid cell model as a macro model, and mixing the poison particle balls and the SiC matrix in the fuel rod according to the effective share to obtain the effective multiplication factor and the average flux of the fuel rod;
s5, performing fuel consumption calculation according to the effective multiplication factors and the average flux to obtain the new nuclear density N of the fuel particles;
s6, in a microscopic fine model, performing burnup calculation according to the average flux to obtain the new nucleus density n of the poison particle ball;
s7, returning the new nucleus density N and the new nucleus density N to the step S1, and executing the steps S1-S6 to realize the multi-scale coupled burnup calculation.
2. The method of claim 1, wherein the microscopic fine model further comprises a cladding layer and a moderator layer;
the sizes of the poison particle balls, the fuel particles, the SiC matrix, the cladding layer and the moderator layer are all set according to the volume proportion of each part of the fuel rod bundle grid element in the macro, and the material properties of each ball layer are consistent with those of the macro grid element.
3. The multi-scale coupling method for the stratified combustion of dispersed fuel and poison particles as claimed in claim 1, wherein in the microscopic fine model, the poison particle spheres appear as a single geometric structure without double heterogeneity; the fuel particles are dispersed in the SiC matrix and need to be subjected to double heterogeneous treatment;
the neutron transport calculation is carried out on the microscopic fine model to obtain an effective value-added factor, and the method comprises the following steps:
and performing neutron transport calculation on the microscopic fine model with double heterogeneity by adopting a collision probability method to obtain effective value-added factors of the microscopic fine model.
4. The method of claim 1, wherein said reactivity being the same means that the microscopic fine model and the microscopic equivalent homogeneous model correspond to an effective multiplication factor that is equal.
5. The method of claim 4, wherein the effective multiplication factors corresponding to the microscopic fine model and the microscopic equivalent homogeneous model are equalized by:
uniformly mixing the poison particle balls into the SiC matrix, and adjusting the nuclear density of the SiC matrix to keep the total amount of the SiC matrix unchanged;
adjusting the effective average nuclear density of the poison particle balls to enable the effective multiplication factor of the microscopic equivalent homogeneous model to be equal to that of the microscopic fine model;
wherein the effective average nucleus density obtained by the adjustment is used as the effective share of the poison particle balls.
6. The method of claim 1, wherein said obtaining an effective multiplication factor and an average flux of fuel rods comprises:
and calling a grid cell program for neutron physical analysis, and performing neutron transport calculation on the whole rod bundle grid cell by adopting a collision probability method to obtain the effective multiplication factor and the average flux of the fuel rod.
7. The multi-scale coupling method for the stratified combustion of dispersed fuel and poison particles as claimed in claim 1, wherein said poison particle spheres are divided into 10 layers.
8. The method of claim 1, wherein the cladding layer and the moderator layer remain unchanged in the microscopic equivalent homogenous model.
9. The method of claim 1, wherein light water is used in a pressurized water reactor of the moderator layer.
10. A multi-scale coupling device for stratified burnup of dispersed fuel and poison particles, comprising:
at least one processor;
at least one memory for storing at least one program;
when executed by the at least one processor, cause the at least one processor to implement the method of any one of claims 1-9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111092461.2A CN113887027B (en) | 2021-09-17 | 2021-09-17 | Multi-scale coupling method and device for layered burnup of dispersed fuel and poison particles |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111092461.2A CN113887027B (en) | 2021-09-17 | 2021-09-17 | Multi-scale coupling method and device for layered burnup of dispersed fuel and poison particles |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113887027A true CN113887027A (en) | 2022-01-04 |
| CN113887027B CN113887027B (en) | 2024-03-26 |
Family
ID=79009489
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111092461.2A Active CN113887027B (en) | 2021-09-17 | 2021-09-17 | Multi-scale coupling method and device for layered burnup of dispersed fuel and poison particles |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113887027B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110120270A (en) * | 2019-04-17 | 2019-08-13 | 中山大学 | A method of for the important operating parameter of presurized water reactor of collecting evidence of tracing to the source |
| CN115455796A (en) * | 2022-09-22 | 2022-12-09 | 中国核动力研究设计院 | Transport equivalence method for dispersed boron particles, electronic device and storage medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150170767A1 (en) * | 2013-11-07 | 2015-06-18 | Francesco Venneri | FULLY CERAMIC MICRO-ENCAPSULATED (FCM) FUEL FOR CANDUs AND OTHER REACTORS |
| CN106782684A (en) * | 2017-02-09 | 2017-05-31 | 中科瑞华原子能源技术有限公司 | For the fixed fuel assembly of wrapping wire and lead base reactor of high power density reactor core |
| CN109493924A (en) * | 2018-11-13 | 2019-03-19 | 西安交通大学 | A method of obtaining the effective multigroup cross section of FCM fuel |
| US20200089826A1 (en) * | 2018-09-14 | 2020-03-19 | Northwestern University | Integrated process-structure-property modeling frameworks and methods for design optimization and/or performance prediction of material systems and applications of same |
| CN113270206A (en) * | 2021-03-29 | 2021-08-17 | 中国核电工程有限公司 | Small-sized prismatic annular gas-cooled micro-reactor core system with densely arranged coolant channels |
-
2021
- 2021-09-17 CN CN202111092461.2A patent/CN113887027B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150170767A1 (en) * | 2013-11-07 | 2015-06-18 | Francesco Venneri | FULLY CERAMIC MICRO-ENCAPSULATED (FCM) FUEL FOR CANDUs AND OTHER REACTORS |
| CN106782684A (en) * | 2017-02-09 | 2017-05-31 | 中科瑞华原子能源技术有限公司 | For the fixed fuel assembly of wrapping wire and lead base reactor of high power density reactor core |
| US20200089826A1 (en) * | 2018-09-14 | 2020-03-19 | Northwestern University | Integrated process-structure-property modeling frameworks and methods for design optimization and/or performance prediction of material systems and applications of same |
| CN109493924A (en) * | 2018-11-13 | 2019-03-19 | 西安交通大学 | A method of obtaining the effective multigroup cross section of FCM fuel |
| CN113270206A (en) * | 2021-03-29 | 2021-08-17 | 中国核电工程有限公司 | Small-sized prismatic annular gas-cooled micro-reactor core system with densely arranged coolant channels |
Non-Patent Citations (1)
| Title |
|---|
| 娄磊;姚栋;柴晓明;于颖锐;彭星杰;刘琨;田超;王晨琳;刘勇;谢运利;: "弥散燃料与可燃毒物双重非均匀性分析", 中国核电, no. 01, 15 February 2020 (2020-02-15) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110120270A (en) * | 2019-04-17 | 2019-08-13 | 中山大学 | A method of for the important operating parameter of presurized water reactor of collecting evidence of tracing to the source |
| CN115455796A (en) * | 2022-09-22 | 2022-12-09 | 中国核动力研究设计院 | Transport equivalence method for dispersed boron particles, electronic device and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113887027B (en) | 2024-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113887027B (en) | Multi-scale coupling method and device for layered burnup of dispersed fuel and poison particles | |
| CN103366836B (en) | Fuel ball, method for making and nuclear reactor | |
| CN103189925B (en) | The nuclear fuel cluster of thoriated and comprise the nuclear reactor of this nuclear fuel cluster | |
| CN114242283A (en) | Fuel loading method for 24-month refueling period of pressurized water reactor core | |
| Tran et al. | Feasibility of using Gd2O3 particles in VVER-1000 fuel assembly for controlling excess reactivity | |
| Dobisesky | Reactor physics considerations for implementing silicon carbide cladding into a PWR environment | |
| JP4813876B2 (en) | Core design support device and program | |
| CN113963759B (en) | Multi-scale coupling calculation method and device for layered burnup of dispersed burnup poison particles | |
| JP2008170454A (en) | Mox fuel assembly for pressurized water reactor | |
| CN114048664A (en) | Double-heterogeneous-system reactivity equivalent physical transformation calculation method and module | |
| Chukbar | Verification of statistical method CORN for modeling of microfuel in the case of high grain concentration | |
| Girieud et al. | Science Version 2: the most recent capabilities of the Framatome 3-D nuclear code package | |
| Reiss et al. | Full-core SSCWR calculations applying a fast computational method | |
| CN117133490B (en) | Method and system for shortening establishment process of pebble-bed high-temperature air-cooled primary-loading core | |
| Yang et al. | Effects of the central graphite column dimension and pebble size on power density distribution in annular core pebble‐bed HTR | |
| Nagai et al. | Small PWRs using coated particle fuel for district heating, PFPWR50 | |
| Yuan et al. | Double heterogeneity effect analysis for plate fuel with dispersed particles | |
| Shaukat et al. | Research Article Multiobjective Core Reloading Pattern Optimization of PARR-1 Using Modified Genetic Algorithm Coupled with Monte Carlo Methods | |
| CN115455796A (en) | Transport equivalence method for dispersed boron particles, electronic device and storage medium | |
| Thanh et al. | Neutronics feasibility of using Gd₂O₃ particles in VVER-1000 fuel assembly | |
| Zhong et al. | Neutronic analysis of fuel assembly loading patterns in small integrated pressurized water reactor core and fuel shifting strategy concept design | |
| US20110158375A1 (en) | Fuel assembly | |
| Greenquist | UO2-liquid metal suspension fuel concept for enhanced passive safety of LWRs: A heat pipe case study | |
| Kozlowski et al. | Pressurised Water Reactor MOX/UO2 Core Transient Benchmark Final Report | |
| Fetuata | Optimization of Burnable Poison Pattern Design |
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 |