Simulation method for equivalent thermal conductivity of fuel pellets

Abstract

This application proposes a simulation calculation method for the equivalent thermal conductivity of a fuel core, which includes creating a representative volume element geometric model of the fuel core. The representative volume element geometric model includes a matrix and multiple randomly distributed fuel particles. The volume content of the multiple fuel particles in the matrix is ​​a preset volume content. The outer surface of each fuel particle is covered by a reaction layer of preset thickness. The representative volume element geometric model is meshed. Based on the meshed representative volume element geometric model, according to the material parameters of the representative volume element geometric model, the first temperature and the second temperature applied to the relative first and second sides of the representative volume element geometric model, the average heat flux density of the representative volume element geometric model is simulated and calculated. Based on the average heat flux density, the first temperature, the second temperature, and the distance between the first and second sides, the equivalent thermal conductivity of the fuel core is calculated.

Description

Technical Field

[0001] This application mainly relates to the field of nuclear power, and in particular to a simulation calculation method for the equivalent thermal conductivity of fuel cores, electronic equipment, and computer storage media. Background Technology

[0002] U3Si2-Al plate-type fuel elements use high-density U3Si2 as the fuel phase, dispersed within a high thermal conductivity, easily machinable aluminum alloy matrix. They possess advantages such as high uranium density, good irradiation stability, and high thermal conductivity, and are widely used in high-flux experimental and research reactors. During reactor operation, under strong neutron irradiation, U3Si2 particles chemically react with the aluminum matrix, forming a reaction layer dominated by the Al-Si-U ternary phase. The thermal conductivity of this reaction layer is less than 10% of that of pure aluminum, and it has weak fission gas retention capacity, easily leading to bubble aggregation. As burnup progresses, the reaction layer thickens, and the highly thermally conductive aluminum matrix is ​​continuously consumed, causing a gradual decrease in the overall equivalent thermal conductivity of the fuel core. This, in turn, leads to increased fuel temperature, uneven thermal stress distribution, and exacerbates risks such as fuel swelling and cladding creep, affecting element safety and structural integrity.

[0003] In addition, due to factors such as strong radioactivity and sample preparation and testing environment, it is extremely difficult to directly measure the thermal conductivity of fuel cores under different burn-out conditions. Existing methods are unable to systematically measure the variation of thermal conductivity with key parameters such as burn-out and reaction layer thickness. Summary of the Invention

[0004] To address the aforementioned issues, this application proposes a simulation calculation method, electronic device, and computer storage medium for the equivalent thermal conductivity of a fuel core, in order to more accurately obtain the thermal conductivity of a fuel core with a reaction layer.

[0005] In a first aspect, this application proposes a simulation calculation method for the equivalent thermal conductivity of a fuel core, comprising: creating a representative volume element geometric model of the fuel core, the representative volume element geometric model including a matrix and a plurality of randomly distributed fuel particles, the volume content of the plurality of fuel particles in the matrix being a preset volume content, and the outer surface of each fuel particle being covered by a reaction layer of preset thickness; meshing the representative volume element geometric model; based on the meshed representative volume element geometric model, and according to the material parameters of the representative volume element geometric model, a first temperature and a second temperature respectively applied to a first side and a second side of the representative volume element geometric model, simulating and calculating the average heat flux density of the representative volume element geometric model; and calculating the equivalent thermal conductivity of the fuel core according to the average heat flux density, the first temperature, the second temperature, and the distance between the first side and the second side.

[0006] In some embodiments, the step of creating a representative volumetric geometry model of the fuel core includes: S1: Create the base; S2: Based on the preset average radius and preset standard deviation of the radius of the fuel particles, a number of fuel particles and their corresponding center coordinates are randomly generated in the matrix, and the reaction layer is generated. S3: Determine whether the number of times S2 has been executed is greater than 1. If the result is yes, then execute S4; if the result is no, then execute S2. S4: Determine whether the fuel particles generated in the current execution interfere with the fuel particles generated in the previous execution. If the determination result is yes, return to S2. If the determination result is no, execute S5. S5: Determine whether the reaction layer generated in the current execution interferes with the reaction layer generated in the previous execution. If the determination result is yes, return to S2. If the determination result is no, execute S6. S6: Determine whether the volume content of the plurality of fuel particles in the matrix reaches the preset volume content. If the determination result is no, then execute S2. If the determination result is yes for all of them, then end the process.

[0007] In some embodiments, the substrate is a cube, wherein the side length of the cube is at least 10 times the preset average radius.

[0008] In some embodiments, during the step of meshing the representative volumetric geometry model, the mesh size of the reaction layer is no greater than one-third of the preset thickness, and the mesh size of the matrix and the mesh size of the fuel particles are no greater than the preset thickness.

[0009] In some embodiments, the material parameters of the representative volumetric geometry model include the thermal conductivity, specific heat and density of the matrix, the thermal conductivity, specific heat and density of the fuel particles, and the thermal conductivity, specific heat and density of the reaction layer.

[0010] In some embodiments, the equivalent thermal conductivity is the average of the equivalent thermal conductivity obtained from multiple simulation calculations.

[0011] In some embodiments, the preset thickness is calculated according to the following formula:

[0012] Where Y is the preset thickness, f is the fission rate, T is the reaction temperature, R is the gas constant, A is the reaction layer thickness coefficient, Q is the interdiffusion activation energy, p is the fission rate correlation index, and t is the reaction time.

[0013] In some embodiments, the fuel core is plate-shaped, the fuel particles are made of U3Si2 material, and the matrix is ​​made of aluminum.

[0014] Secondly, this application also proposes an electronic device, comprising: a memory for storing instructions executable by a processor; and a processor for executing the instructions to implement the simulation calculation method for the equivalent thermal conductivity of a fuel core as described above.

[0015] Thirdly, this application also proposes a computer storage medium storing computer program code, which, when executed by a processor, implements the simulation calculation method for the equivalent thermal conductivity of a fuel core as described above.

[0016] Compared with the prior art, the beneficial effects of this application are as follows: This application uses simulation software to create geometric models of fuel particles, reaction layers and matrix, and sets their physical parameters to simulate and calculate the equivalent thermal conductivity of the fuel core. It can efficiently and accurately calculate the equivalent thermal conductivity of the fuel core, and has the effects of being easy to operate and repeatable. Attached Figure Description

[0017] The accompanying drawings are included to provide a further understanding of this application; they are incorporated into and constitute a part of this application. The drawings illustrate embodiments of this application and, together with this specification, serve to explain the principles of this application. In the drawings: Figure 1 This is a flowchart illustrating the simulation calculation method for the equivalent thermal conductivity of the fuel core provided in this application embodiment; Figure 2 This is a schematic diagram of a representative volume element geometric model of the fuel core provided in the embodiments of this application; Figure 3 This is a cross-sectional schematic diagram of the composite fuel particles provided in the embodiments of this application; Figure 4 This is a flowchart illustrating the method for creating a representative volumetric geometric model of a fuel core according to an embodiment of this application; Figure 5 This is a schematic diagram of the structure of the electronic device provided in the embodiments of this application.

[0018] Reference numerals: Representative volumetric geometric model 200, composite fuel particle 210, fuel particle 211, reaction layer 212, substrate 220, system 400, communication bus 401, processor 402, read-only memory 403, random access memory 404, communication port 405, hard disk 406. Detailed Implementation

[0019] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are merely some examples or embodiments of this application. For those skilled in the art, these drawings can be applied to other similar scenarios without creative effort. Unless obvious from the context or otherwise specified, the same reference numerals in the drawings represent the same structures or operations.

[0020] As indicated in this application, unless the context clearly indicates otherwise, the words "a," "an," "an," and / or "the" do not specifically refer to the singular and may also include the plural. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of explicitly identified steps and elements, which do not constitute an exclusive list, and the method or apparatus may also include other steps or elements.

[0021] Unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps described in these embodiments do not limit the scope of this application. It should also be understood that, for ease of description, the dimensions of the various parts shown in the drawings are not drawn to actual scale. Techniques, methods, and devices known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and devices should be considered part of the specification. In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limitations. Therefore, other examples of exemplary embodiments may have different values. It should be noted that similar reference numerals and letters in the following drawings denote similar items; therefore, once an item is defined in one drawing, it need not be further discussed in subsequent drawings.

[0022] In the description of this application, it should be understood that the orientation or positional relationship indicated by directional terms such as "front, back, up, down, left, right", "horizontal, vertical, horizontal" and "top, bottom" is usually based on the orientation or positional relationship shown in the accompanying drawings, and is only for the convenience of describing this application and simplifying the description. Unless otherwise stated, these directional terms do not indicate or imply that the device or element referred to must have a specific orientation or be constructed and operated in a specific orientation, and therefore should not be construed as a limitation on the scope of protection of this application; the directional terms "inner" and "outer" refer to the inner and outer contours relative to the outline of each component itself.

[0023] For ease of description, spatial relative terms such as "above," "on top of," "on the upper surface of," "above," etc., are used herein to describe the spatial positional relationship of a device or feature as shown in the figures to other devices or features. It should be understood that spatial relative terms are intended to encompass different orientations in use or operation beyond the orientation of the device as described in the figures. For example, if the device in the figures were inverted, a device described as "above" or "on top of" other devices or structures would subsequently be positioned as "below" or "under" other devices or structures. Thus, the exemplary term "above" can include both "above" and "below." The device may also be positioned in other different ways (rotated 90 degrees or in other orientations), and the spatial relative descriptions used herein will be interpreted accordingly.

[0024] Furthermore, it should be noted that the use of terms such as "first" and "second" to define components is merely for the purpose of distinguishing the corresponding components. Unless otherwise stated, these terms have no special meaning and therefore should not be construed as limiting the scope of protection of this application. In addition, although the terminology used in this application is selected from commonly known and used terms, some terms mentioned in this application's specification may have been chosen by the applicant according to his or her judgment, and their detailed meanings are explained in the relevant sections of this description. Moreover, this application should be understood not only through the actual terms used, but also through the meaning implied by each term.

[0025] Flowcharts are used in this application to illustrate the operations performed by the system according to embodiments of this application. It should be understood that the preceding or following operations are not necessarily performed in exact order. Instead, various steps can be processed in reverse order or simultaneously. Furthermore, other operations may be added to these processes, or one or more steps may be removed from these processes.

[0026] The simulation calculation method for the equivalent thermal conductivity of the fuel core proposed in this application will be explained next through specific embodiments.

[0027] The simulation calculation method for the equivalent thermal conductivity of the fuel core in one embodiment of this application includes: Figure 1 Steps S101 to S104 are explained below.

[0028] refer to Figure 2 and Figure 3In step S101, a representative volumetric geometric model 200 of the fuel core is created. The representative volumetric geometric model 200 includes a substrate 220 and multiple composite fuel particles 210 randomly distributed within the substrate 220. Each composite fuel particle 210 includes a fuel particle 211 and a reaction layer 212 encapsulating the fuel particle 211. The volume fraction of the multiple fuel particles 211 in the substrate 220 is a preset volume fraction, and the thickness of the reaction layer 212 is a preset thickness. The reaction layer is generated by the reaction between the substrate 220 and the fuel particles 211.

[0029] The representative volumetric geometric model 200 is a physical model created in simulation software based on the physical properties of the fuel core. The representative volumetric geometric model 200 is established in simulation software, such as ABAQUS.

[0030] In some embodiments, a preset volume is set according to the actual fuel core being simulated, and the preset volume may be related to the type of uranium-containing material of the fuel particles 211 in the actual fuel core.

[0031] In some embodiments, the uranium-containing material in the fuel particles 211 is U3Si2, and the preset volume is calculated using the following formula: (1) in, This refers to the volume content of U3Si2 fuel particles 211 in the core. The density of metallic uranium (U) per unit volume of the core, expressed in g / cm³. 3 , This represents the mass percentage of metallic uranium in U3Si2. The density constant of U3Si2 per unit volume of the core is expressed in g / cm³. 3 .

[0032] In some embodiments, the preset thickness of the reaction layer 212 is calculated based on the actual working environment of the simulated actual fuel core, using the following formula: (2) Where Y is the preset thickness, f is the fission rate, T is the reaction temperature, R is the gas constant, A is the thickness coefficient of reaction layer 212, Q is the interdiffusion activation energy, p is the fission rate correlation index, and t is the reaction time. From formula (2), it can be concluded that the simulation algorithm of this application can truly reflect the thickness of reaction layer 212 under different reaction times and different reaction temperatures, that is, it can calculate the thickness of reaction layer 212 under multiple operating conditions, thus improving the accuracy of calculating the equivalent thermal conductivity of fuel core.

[0033] In some embodiments, the step of creating a representative volumetric geometry model 200 of the fuel core includes Figure 4 Steps S301 to S306 are detailed below.

[0034] refer to Figure 2 and Figure 3 In step S301, base 220 is created.

[0035] In this step, the substrate 220 can be created using simulation software. Figure 2 The base 220 is a cube. It should be noted that this application does not limit the shape of the base 220; it can also be a cuboid, etc.

[0036] In step S302, based on the preset average radius and the preset standard deviation of the radius of the fuel particles 211, a number of fuel particles 211 and their corresponding center coordinates are randomly generated in the substrate 220, and a reaction layer 212 is generated.

[0037] In this step, a preset average radius and a preset standard deviation of the radius can be set in the simulation software according to actual needs for simulating fuel particles 211. The simulation software randomly generates the radius corresponding to each fuel particle 211 based on the preset average radius and the preset standard deviation of the radius. Understandably, fuel particles 211 made of different materials have different preset average radii and preset standard deviations of the radius. The coordinates of the center of each fuel particle 211 determine its spatial position in the matrix 220, and the radius of each fuel particle 211 determines its volume in the matrix 220.

[0038] In some embodiments, when the base 220 is a cube, the side length of the cube is at least 10 times the preset average radius.

[0039] In some embodiments, a fuel particle 211 and its corresponding center coordinates are randomly generated at a time in the substrate 220. In other embodiments, multiple fuel particles 211 and their corresponding center coordinates can be randomly generated at a time.

[0040] In step S303, it is determined whether the current cumulative number of executions in step S302 is greater than 1. If the determination result is yes, then step S304 is executed; if the determination result is no, then step S302 is executed.

[0041] In step S304, it is determined whether the fuel particle 211 generated in the current step S302 interferes with the fuel particle 211 generated in the previous step S302. If interference occurs, the fuel particle 211 generated this time is deleted, and the process returns to step S302. If no interference occurs, step S305 is executed.

[0042] In step S305, it is determined whether the reaction layer 212 generated in the current step S302 interferes with the reaction layer 212 generated in the previous step S302. If interference occurs, the fuel particles 211 and the reaction layer 212 generated this time are deleted, and the process returns to step S302. If no interference occurs, step S306 is executed.

[0043] In step S306, it is determined whether the volume content of the multiple generated fuel particles 211 in the matrix 220 reaches the preset volume content. If the determination result is no, the process returns to step S302. If the determination result is yes, the process ends.

[0044] In some embodiments, reference Figure 2 and Figure 3 In the simulation software, Boolean operations are used to merge all fuel particles 211, reaction layer 212 and matrix 220. For fuel particles 211 and their surface reaction layer 212, only the surface of reaction layer 212 is retained, thereby obtaining composite fuel particles 210 containing reaction layer 212. Finally, a representative volume element geometric model 200 in which composite fuel particles 210 are randomly distributed in matrix 220 is obtained.

[0045] In step S102, the representative volumetric geometric model 200 is meshed.

[0046] Specifically, firstly, the representative volumetric geometry model 200 is divided into independent components according to the boundaries of fuel particles 211 and reaction layer 212. Then, a relatively large base mesh size is set for the entire representative volumetric geometry model 200, while a smaller local mesh size is set for the thinner reaction layer 212 to ensure that its internal temperature gradient can be accurately resolved.

[0047] In some embodiments, the mesh size of the reaction layer 212 in the representative volumetric geometry model 200 is no greater than one-third of the preset thickness, and the mesh sizes of the substrate 220 and the fuel particles 211 are no greater than the preset thickness. This allows the simulation results to more closely approximate reality and enables more accurate capture and calculation of the characteristics of the reaction layer 212.

[0048] In step S103, based on the representative volume element geometric model 200 after mesh generation, and according to the material parameters of the representative volume element geometric model 200, the first temperature T1 and the second temperature T1 applied to the relative first side 221 and second side 222 of the representative volume element geometric model 200 respectively, the average heat flux density of the representative volume element geometric model 200 along the first direction D1 is calculated by simulation. .

[0049] In this step, a first temperature T1 and a second temperature T2 are applied to opposite sides of the representative volumetric geometric model 200. For example, if the representative volumetric geometric model 200 is a cube, the first side 221 and the second side 222 are opposite faces of the cube, where the side length of the cube is L. The first side 221 can be the face with Z=0 in the first direction D1, and the second side 222 is the face with Z=L. The first temperature T1 and the second temperature T2 are applied to the opposite faces, i.e., the faces with Z=0 and Z=L. In some embodiments, the material parameters of the representative volumetric geometric model 200 include the thermal conductivity, specific heat, and density of the matrix 220, the thermal conductivity, specific heat, and density of the fuel particles 211, and the thermal conductivity, specific heat, and density of the reaction layer 212. The above parameters are set in the simulation software, and then heat transfer calculations are performed according to the algorithm in the simulation software to obtain the temperature distribution of the representative volumetric geometric model 200, including the average heat flux density along the first direction D1 on the Z=L face. .

[0050] In step S104, based on the average heat flux density Calculate the equivalent thermal conductivity of the fuel core based on the first temperature T1, the second temperature T2, and the distance between the first side 221 and the second side 222.

[0051] Specifically, based on the average heat flux density obtained in step S103 And according to Fourier's law, the method for calculating the equivalent thermal conductivity of the fuel core is as follows: (3) in, The equivalent thermal conductivity of the fuel core, Let T1 be the average heat flux density along the first direction D1 on the Z=L surface, T2 be the first temperature, T2 be the second temperature, and L be the side length of the cube.

[0052] In some embodiments, the equivalent thermal conductivity is the average of the equivalent thermal conductivity obtained from multiple simulation calculations. Specifically, the fuel particles 211 in the representative volumetric geometry model 200 are randomly generated each time, and their distribution within the matrix 220 is different each time. By repeating steps S101 to S104 multiple times to calculate the equivalent thermal conductivity of the fuel core, and then averaging the results of multiple simulation calculations, a more accurate equivalent thermal conductivity of the fuel core can be obtained. In some embodiments, the equivalent thermal conductivity is the average of the equivalent thermal conductivity obtained from more than 10 simulation calculations.

[0053] Thus, this application uses simulation software to create geometric models of fuel particles 211, reaction layer 212 and substrate 220, and sets their physical parameters to simulate and calculate the equivalent thermal conductivity of the fuel core. This method can efficiently and accurately calculate the equivalent thermal conductivity of the fuel core, and it is easy to operate and repeatable. Furthermore, this simulation calculation method can also be used to predict the equivalent thermal conductivity of the fuel core under various operating conditions.

[0054] In some embodiments, the fuel core is a plate-shaped fuel core, wherein the material of the fuel particles 211 in the plate-shaped fuel core includes U3Si2, and the material of the matrix 220 includes aluminum (Al). When the fuel particles 211 include U3Si2 and the matrix 220 includes aluminum, the corresponding reaction layer 212 is composed of an amorphous or microcrystalline Al-Si-U ternary phase.

[0055] Figure 5 This is a schematic diagram of a simulation calculation method for the equivalent thermal conductivity of a fuel core and an electronic device according to an embodiment of this application. (Reference) Figure 5 As shown, system 400 is used to implement Figure 1 The method shown includes an internal communication bus 401, a processor 402, a read-only memory (ROM) 403, a random access memory (RAM) 404, a communication port 405, and a hard disk 406. The internal communication bus 401 enables data communication between components of the system 400. The processor 402 can make judgments and issue prompts. In some embodiments, the processor 402 may consist of one or more processors. The communication port 405 enables data communication between the system 400 and external systems. In some embodiments, the system 400 can send and receive information and data from a network through the communication port 405. The system 400 may also include different forms of program storage units and data storage units, such as the hard disk 406, the read-only memory (ROM) 403, and the random access memory (RAM) 404, capable of storing various data files used for computer processing and / or communication, as well as possible program instructions executed by the processor 402. The processor executes these instructions to implement the main part of the method. The above-mentioned simulation calculation method for the equivalent thermal conductivity of the fuel core can be implemented as a computer program, stored in the hard disk 406, and loaded into the processor 402 for execution.

[0056] This application also includes a computer-readable medium storing computer program code that, when executed by a processor, implements the aforementioned simulation calculation method for the equivalent thermal conductivity of a fuel core.

[0057] The basic concepts have been described above. Obviously, for those skilled in the art, the above disclosure is merely illustrative and does not constitute a limitation of this application. Although not explicitly stated herein, those skilled in the art may make various modifications, improvements, and corrections to this application. Such modifications, improvements, and corrections are suggested in this application, and therefore remain within the spirit and scope of the exemplary embodiments of this application.

[0058] Furthermore, this application uses specific terms to describe embodiments of the application. For example, "an embodiment," "one embodiment," and / or "some embodiments" refer to a particular feature, structure, or characteristic related to at least one embodiment of the application. Therefore, it should be emphasized and noted that "an embodiment," "one embodiment," or "an alternative embodiment" mentioned twice or more in different locations in this specification do not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of the application can be appropriately combined.

[0059] Some aspects of this application can be executed entirely by hardware, entirely by software (including firmware, resident software, microcode, etc.), or by a combination of hardware and software. The aforementioned hardware or software may be referred to as a "data block," "module," "engine," "unit," "component," or "system." The processor may be one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DAPDs), programmable logic devices (PLDs), field-programmable gate arrays (FPGAs), processors, controllers, microcontrollers, microprocessors, or combinations thereof. Furthermore, aspects of this application may manifest as computer products residing in one or more computer-readable media, including computer-readable program code. For example, computer-readable media may include, but are not limited to, magnetic storage devices (e.g., hard disks, floppy disks, magnetic tapes, etc.), optical discs (e.g., compressed CDs, digital multifunction DVDs, etc.), smart cards, and flash memory devices (e.g., cards, sticks, key drives, etc.).

[0060] A computer-readable medium may contain a propagated data signal containing computer program code, for example, on baseband or as part of a carrier wave. This propagated signal may take various forms, including electromagnetic, optical, and so on, or suitable combinations thereof. A computer-readable medium can be any computer-readable medium other than a computer-readable storage medium, which can be connected to an instruction execution system, apparatus, or device to enable communication, propagation, or transmission of a program for use. The program code located on the computer-readable medium can be propagated through any suitable medium, including radio, cable, fiber optic cable, radio frequency signals, or similar media, or any combination of the above media.

[0061] Similarly, it should be noted that, in order to simplify the description of the present application and thus aid in the understanding of one or more embodiments of the invention, the foregoing description of the embodiments of the present application sometimes combines multiple features into a single embodiment, drawing, or description thereof. However, this disclosure method does not imply that the subject matter of the application requires more features than those mentioned in the claims. In fact, the embodiments contain fewer features than all the features of the single embodiments disclosed above.

[0062] In some embodiments, numbers describing the quantity of components and attributes are used. It should be understood that such numbers used in the description of embodiments are modified in some examples with the terms "approximately," "approximately," or "generally." Unless otherwise stated, "approximately," "approximately," or "generally" indicates that the numbers are allowed to vary by ±20%. Accordingly, in some embodiments, the numerical parameters used in the specification and claims are approximate values, which may be changed depending on the characteristics required by individual embodiments. In some embodiments, numerical parameters should take into account specified significant digits and employ a general method of digit reservation. Although the numerical ranges and parameters used to confirm their breadth of scope in some embodiments of this application are approximate values, in specific embodiments, such values ​​are set as precisely as feasible.

[0063] Although this application has been described with reference to specific embodiments, those skilled in the art should recognize that the above embodiments are only used to illustrate this application, and various equivalent changes or substitutions can be made without departing from the spirit of this application. Therefore, any changes or modifications to the above embodiments within the essential spirit of this application will fall within the scope of the claims of this application.

Claims

1. A simulation calculation method for the equivalent thermal conductivity of a fuel core, characterized in that, include: Create a representative volume element geometric model of the fuel core. The representative volume element geometric model includes a matrix and multiple randomly distributed fuel particles. The volume content of the multiple fuel particles in the matrix is ​​a preset volume content. The outer surface of each fuel particle is covered by a reaction layer of preset thickness. The representative volumetric geometric model is meshed; Based on the representative volume element geometric model after mesh generation, the average heat flux density of the representative volume element geometric model is calculated by simulation according to the material parameters of the representative volume element geometric model and the first and second temperatures applied to the relative first and second sides of the representative volume element geometric model. The equivalent thermal conductivity of the fuel core is calculated based on the average heat flux density, the first temperature, the second temperature, and the distance between the first side and the second side.

2. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, The steps for creating a representative volumetric geometry model of the fuel core include: S1: Create the base; S2: Based on the preset average radius and preset standard deviation of the radius of the fuel particles, a number of fuel particles and their corresponding center coordinates are randomly generated in the matrix, and the reaction layer is generated. S3: Determine whether the number of times S2 has been executed is greater than 1. If the result is yes, then execute S4; if the result is no, then execute S2. S4: Determine whether the fuel particles generated in the current execution interfere with the fuel particles generated in the previous execution. If the determination result is yes, return to S2. If the determination result is no, execute S5. S5: Determine whether the reaction layer generated in the current execution interferes with the reaction layer generated in the previous execution. If the determination result is yes, return to S2. If the determination result is no, execute S6. S6: Determine whether the volume content of the plurality of fuel particles in the matrix reaches the preset volume content. If the determination result is no, then execute S2. If the determination result is yes for all of them, then end the process.

3. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 2, characterized in that, The substrate is a cube, wherein the side length of the cube is at least 10 times the preset average radius.

4. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, In the step of meshing the representative volumetric geometric model, the mesh size of the reaction layer is no greater than one-third of the preset thickness, and the mesh size of the matrix and the mesh size of the fuel particles are no greater than the preset thickness.

5. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, The material parameters of the representative volumetric geometry model include the thermal conductivity, specific heat and density of the matrix, the thermal conductivity, specific heat and density of the fuel particles, and the thermal conductivity, specific heat and density of the reaction layer.

6. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, The equivalent thermal conductivity is the average value of the equivalent thermal conductivity obtained from multiple simulation calculations.

7. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, The preset thickness is calculated using the following formula: Where Y is the preset thickness, f is the fission rate, T is the reaction temperature, R is the gas constant, A is the reaction layer thickness coefficient, Q is the interdiffusion activation energy, p is the fission rate correlation index, and t is the reaction time.

8. The simulation calculation method for the equivalent thermal conductivity of a fuel core as described in claim 1, characterized in that, The fuel core is plate-shaped, the fuel particles are made of U3Si2, and the matrix is ​​made of aluminum.

9. An electronic device, comprising: Memory is used to store instructions that can be executed by the processor; as well as A processor for executing the instructions to implement the simulation calculation method for the equivalent thermal conductivity of a fuel core as described in any one of claims 1-8.

10. A computer storage medium storing computer program code, which, when executed by a processor, implements the simulation calculation method for the equivalent thermal conductivity of a fuel core as described in any one of claims 1-8.

Images