Radiation transfer coefficient calculation method, system and device and medium

Abstract

The invention discloses a radiation transfer coefficient calculation method, system and device and a medium. According to the method, the radiation transfer coefficient is calculated through the Monte Carlo method, heat radiation energy of a certain surface is regarded as a large number of energy beams emitted by the surface, energy carried by the energy beams is equal, and the energy emitted by the radiation surface is directly determined by the number of the energy beams emitted by the surface. Secondly, the possible path of each energy beam is tracked until the energy beam is finally absorbed by a certain surface or enters the environment space; and tracking the strokes of all the energy beams, counting the number of the energy beams which are emitted from one surface, arrive at the other surface and are absorbed by the surface, and calculating the radiation transfer coefficient between the surfaces through a statistical value. According to the method, the complex geometric configuration can be processed when the radiation transfer coefficient is calculated, the radiation reflection effect between the surfaces in the geometric configuration is considered, and the heat transfer characteristic of the three-dimensional lattice sandwich structure can be truly described. The method is mainly used in the technical field of materials.

Description

Technical Field

[0001] The present invention relates to the field of material technology, and in particular to a method, system, device and medium for calculating a radiation transfer coefficient. Background Art

[0002] Ceramic-based lattice structure composite materials, as a functional structure that combines lightweight, high strength, and high temperature resistance, have been widely used in the thermal protection structure of hypersonic aircraft. Therefore, it is very necessary to deeply study the heat transfer mechanism of ceramic-based lattice structures. However, when traditional methods are used to analyze the heat transfer performance of ceramic-based lattice structures, due to their complex internal core configuration, the solution of the radiation heat transfer integral equation and the radiation transfer coefficient will become extremely difficult. Traditional analytical methods will no longer be applicable, and even common numerical calculation methods will be difficult to work. Therefore, the problems existing in related technologies need to be improved. Summary of the invention

[0003] The present invention provides a method, system, device and medium for calculating the radiation transfer coefficient to solve one or more technical problems existing in the prior art and at least provide a beneficial choice or create conditions.

[0004] A method for calculating the radiation transfer coefficient is provided, including: Step 1: Establish a unit cell model according to the microscopic configuration of the ceramic-based lattice material; Step 2: using the i-th imaginary surface of the unit cell model as the emission surface of the energy beam, obtaining the emission position and emission direction of the energy beam by Monte Carlo method, and counting the number of energy beams emitted from the i-th imaginary surface; Step 3, determining the absorption surface of the energy beam in the unit cell model according to the propagation path of the energy beam, and the absorption surface is recorded as the jth absorption surface; Step 4, counting the number of energy beams absorbed by the j-th absorption surface according to the shielding of the energy beams by the core; Step 5: Calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, where the mathematical model is: ; It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

[0005] Furthermore, the probability distribution model of the emission position of the energy beam in step 2 is: ; ; ; Wherein, x, y, z are the coordinates of the emission point of the energy beam, , , is the maximum coordinate point that can be taken in the area occupied by the unit cell model, , , is the minimum coordinate point that can be taken in the area occupied by the unit cell model, , , is a random number uniformly distributed in the interval [0, 1].

[0006] Furthermore, the probability distribution model of the emission direction of the energy beam in step 2 is: ; ; ; ; in, is the minimum zenith angle that each imaginary surface in the unit cell model can take, is the minimum azimuth angle that each imaginary surface in the unit cell model can take, =90°, , , , , =0°, =270°, =180°, =90°, and is a random number uniformly distributed in the interval [0, 1].

[0007] Furthermore, the method for calculating the propagation path of the energy beam in step 3 is: through the energy beam emission position and emission direction calculated in step 2, a positioning point is taken outside the unit cell model along its emission direction, and the line connecting the positioning point and the emission position of the energy beam is used as the motion trajectory of the energy beam, and the surface of the unit cell model intersecting with the line is the absorption surface.

[0008] Furthermore, the method for calculating the energy beam being blocked by the core in step 4 is: the motion trajectory of the energy beam is regarded as a line segment m, the core is regarded as a cylinder n, and the line segment m is represented by the centroid position vector and the unit vector along the axial direction Indicates that the cylinder n is represented by the center of mass position vector and the unit vector along the axial direction Indicates that the diameter of the cylinder n is d; the shortest distance between the line segment m and the cylinder n is The calculation formula is: , ; in, is perpendicular to and The unit vector of the line segment m; the distance between the point on the axis closest to the cylinder n , the distance between the center of mass of the cylinder n and the point on the axis closest to the line segment m The calculation formula is: ; ; When the following conditions are met, the line segment m intersects with the cylinder n, that is, the energy beam is blocked by the core: .

[0009] Furthermore, when the energy beam is not blocked by the core, it is considered to be absorbed by the j-th absorption surface. When the energy beam is blocked by the core, it is necessary to calculate whether the energy beam is absorbed by the core. The energy beam is calculated by the absorption ratio of the core surface. The calculation formula is: ; in, is a random number uniformly distributed in the interval [0, 1]. is the absorption ratio of the core surface.

[0010] Furthermore, when the energy beam is not absorbed by the core, it is regarded as reflected, and the calculation is returned to step 2.

[0011] In order to achieve the above-mentioned purpose, another aspect of the embodiment of the present application provides a radiation transfer coefficient calculation system, including: a modeling module, a transmission azimuth calculation module, a propagation path calculation module, a statistical module and a radiation transfer coefficient calculation module; Modeling module, used to establish unit cell model according to the microstructure of ceramic-based lattice materials; An emission direction calculation module, used for taking the i-th imaginary surface of the unit cell model as the emission surface of the energy beam, obtaining the emission position and emission direction of the energy beam by Monte Carlo method, and counting the number of energy beams emitted from the i-th imaginary surface; A propagation path calculation module, used to determine the absorption surface of the energy beam in the unit cell model according to the propagation path of the energy beam, wherein the absorption surface is recorded as the jth absorption surface; A statistical module, used for counting the number of energy beams absorbed by the j-th absorption surface according to the shielding of the energy beams by the core; The radiation transfer coefficient calculation module is used to calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, wherein the mathematical model is: ; It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

[0012] In order to achieve the above-mentioned purpose, another aspect of an embodiment of the present application proposes a radiation transfer coefficient calculation device, including a memory and a processor, the memory stores a computer program, and the processor implements any of the above-mentioned radiation transfer coefficient calculation methods when executing the computer program.

[0013] In order to achieve the above-mentioned purpose, another aspect of an embodiment of the present application proposes a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and is characterized in that when the computer program is executed by a processor, it implements any of the above-mentioned methods for calculating the radiation transfer coefficient.

[0014] The present invention at least includes the following beneficial effects: The present invention provides a method, system, device and medium for calculating the radiation transfer coefficient. In view of the fact that the existing traditional methods cannot calculate the radiation transfer coefficient of the complex three-dimensional porous structure, the present invention calculates its radiation transfer coefficient by the Monte Carlo method. That is, the thermal radiation energy of a certain surface is regarded as composed of a large number of energy beams emitted by the surface, and the energy carried by each energy beam is equal, then the energy emitted by the radiation surface will be directly determined by the number of energy beams emitted by this surface. Secondly, the possible path of each energy beam is tracked until the energy beam is finally absorbed by a certain surface or emitted into the ambient space. In this process, the emission position and emission direction of each energy beam are random. The travel of all energy beams is tracked, and the number of energy beams emitted from a certain surface and then reaching another surface and absorbed by the surface is counted. Finally, according to the above statistical values, the radiation transfer coefficient between each surface is calculated. The present application can handle complex geometric configurations and consider the radiation reflection effect between the surfaces in the geometric configuration when calculating the radiation transfer coefficient, and can truly describe the heat transfer characteristics of the three-dimensional lattice sandwich structure. The present invention is mainly used in the field of material technology. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The accompanying drawings are used to provide a further understanding of the technical solution of the present invention and constitute a part of the specification. Together with the embodiments of the present invention, they are used to explain the technical solution of the present invention and do not constitute a limitation on the technical solution of the present invention.

[0016] Figure 1 This is a schematic diagram of the structure where the energy beam is blocked by the core; Figure 2It is a schematic diagram of the intersection relationship between the energy beam and the core; Figure 3 It is a schematic diagram of the unit cell model structure of a tetrahedral pyramid lattice material; Figure 4 Schematic diagram of the unit cell model structure of the TPMS-Gyroid rod structure. DETAILED DESCRIPTION

[0017] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0018] With the rapid development of my country's aerospace industry, the research and development of hypersonic vehicles has received increasing attention in my country. Some complex geometric configurations, such as lattice sandwich structures, a new type of three-dimensional sandwich structure, can control radiation heat transfer by adjusting the density and shape of the lattice, thereby further improving the thermal insulation effect. This structure has high design flexibility and can be customized according to specific usage scenarios and needs to achieve optimal thermal performance. Its excellent thermal insulation and heat preservation effects can be achieved by optimizing the design and material selection. In addition, the lattice sandwich structure also has significant advantages in improving material performance, optimizing design flexibility and reducing weight, so it is widely used in various engineering fields, such as the field of composite materials for aerospace thermal protection systems.

[0019] The main technical purpose of the present invention is to calculate the radiation transfer coefficient in complex geometric configurations.

[0020] In order to achieve this technical purpose, the present invention provides a method, system, device and medium for calculating the radiation transfer coefficient.

[0021] Monte Carlo method is also called statistical simulation method and statistical test method. It is a numerical simulation method that takes probability phenomena as the research object. It is a calculation method that obtains statistical values ​​according to the sampling survey method to infer unknown characteristic quantities. Monte Carlo is a famous casino in Monaco. The method is named to indicate its random sampling nature. Therefore, it is suitable for computational simulation tests on discrete systems. In computational simulation, by constructing a probability model that is similar to the system performance and conducting random tests on a digital computer, the random characteristics of the system can be simulated.

[0022] refer to Figure 1 and Figure 2 , Figure 1 This is a schematic diagram of the structure where the energy beam H is blocked by the core P. Figure 2 It is a schematic diagram of the intersection relationship between the energy beam and the core.

[0023] Example 1: Four-sided pyramid lattice material refer to Figure 3 , Figure 3 It is a schematic diagram of the unit cell model structure of a tetrahedral pyramid lattice material, where A, B, C, D, E, and S represent six imaginary surfaces respectively.

[0024] Step 1: Establish a unit cell model according to the microscopic configuration of the tetrahedral pyramid lattice material. The unit cell structure of the tetrahedral pyramid lattice material has 6 imaginary faces and 4 cores.

[0025] Step 2: Use the i-th imaginary surface of the tetrahedral lattice material unit cell model as the emission surface of the energy beam, obtain the emission position and emission direction of the energy beam through the Monte Carlo method, and count the number of energy beams emitted from the i-th imaginary surface.

[0026] The probability distribution model of the launch position and launch direction is established according to the Monte Carlo method. The probability distribution model formula of the launch position coordinates is: (1); (2); (3); Among them, x, y, z are the coordinates of the emission point of the energy beam, , , is the maximum coordinate point that can be taken in the area occupied by the tetrahedral lattice material unit cell model, , , is the minimum coordinate point that can be obtained in the area occupied by the tetrahedral lattice material unit cell model, , , is a random number uniformly distributed in the interval [0, 1].

[0027] The probability distribution model of the emission direction is: (4); (5); (6); (7); in, It is the minimum zenith angle that each imaginary surface in the unit cell model of the tetrahedral lattice material can take. The zenith angle is the angle between the incident light, the incident light source and the ground normal. It is the minimum azimuth that each imaginary surface in the unit cell model of the tetrahedral lattice material can take. The azimuth is the horizontal angle from the north direction line of a certain point to the target direction line in a clockwise direction. =90°, , , , , =0°, =270°, =180°, =90°, and is a random number uniformly distributed in the interval [0, 1].

[0028] Step 3: Determine the absorption surface of the energy beam in the tetrahedral lattice material unit cell model according to the propagation path of the energy beam. According to the emission position and emission direction of the energy beam calculated in step 2, a positioning point is taken outside the tetrahedral lattice material unit cell model along its emission direction, and the line connecting the positioning point and the energy beam emission position is used as the motion trajectory of the energy beam. The imaginary surface intersecting with the line is the absorption surface, which is recorded as the jth absorption surface.

[0029] Step 4: Calculate the shielding of the energy beam according to the distribution of the cores. The energy beam that is not shielded by the core is absorbed by the jth absorption surface. The motion trajectory of the energy beam is regarded as a line segment m, and the core is regarded as a cylinder n. The line segment m is represented by the centroid position vector and the unit vector along the axial direction Indicates that the cylinder n is represented by the center of mass position vector and the unit vector along the axial direction Indicates that the diameter of cylinder n is d. The shortest distance between line segment m and cylinder n The calculation formula is: , (8); in, is perpendicular to and The unit vector of .

[0030] The distance between the line segment m and the point on the axis closest to the cylinder n , the distance between the center of mass of cylinder n and the point on the axis closest to line segment m The calculation formula is: (9); (10); When the condition of formula 11 is met, the line segment m intersects with the cylinder n, that is, the energy beam is blocked by the core: (11); When the energy beam is not blocked, it is considered to be absorbed by the jth absorption surface.

[0031] When the energy beam is blocked, the absorption ratio of the core surface is used to calculate whether the energy beam is absorbed by the core. The calculation formula is: ; (12); in, is a random number uniformly distributed in the interval [0, 1]. is the absorption ratio of the core surface. When formula 12 is satisfied, the energy beam is absorbed by the core. When formula 12 is not satisfied, the energy beam is reflected, and the process returns to step 2. A new emission position and emission direction are randomly selected from the core surface, and the calculation is performed again.

[0032] Step 5: Calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, where the mathematical model is: (13); in, It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

[0033] Example 2: TPMS-Gyroid rod-shaped structure lattice material refer to Figure 4 , Figure 4 It is a schematic diagram of the unit cell model structure of the TPMS-Gyroid rod-shaped structure lattice material, where A, B, C, D, E, and S represent six imaginary surfaces respectively.

[0034] Step 1: Establish a unit cell model based on the microscopic configuration of the TPMS-Gyroid rod-shaped structure lattice material, which can be expressed by the following formula: (14); in, is the volume fraction of the unit cell model that controls the TPMS-Gyroid rod-shaped structure lattice material, is the size of the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material. Due to the complex geometry of the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material, it needs to be simplified to a rod of equal diameter for calculation. The simplified formula is: (15); Wherein, D is the diameter of the simplified equal-diameter rod, and the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material is established based on this. The unit cell model of the TPMS-Gyroid rod-shaped structure lattice material has 6 imaginary surfaces and 15 cores.

[0035] Step 2: The i-th imaginary surface of the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material is used as the emission surface of the energy beam, the emission position and emission direction of the energy beam are obtained by the Monte Carlo method, and the number of energy beams emitted from the i-th imaginary surface is counted.

[0036] The probability distribution model of the launch position and launch direction is established according to the Monte Carlo method. The probability distribution model formula of the launch position coordinates is: (16); (17); (18); Among them, x, y, z are the coordinates of the emission point of the energy beam, , , is the maximum coordinate point that can be obtained in the area occupied by the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material. , , is the minimum coordinate point that can be obtained in the area occupied by the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material. , , is a random number uniformly distributed in the interval [0, 1].

[0037] The probability distribution model of the emission direction is: (19); (20); (twenty one); (twenty two); in, It is the minimum zenith angle that each imaginary surface can take in the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material. The zenith angle is the angle between the incident light, the incident light source and the ground normal. It is the minimum azimuth that each imaginary surface in the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material can take. The azimuth is the horizontal angle from the north direction line of a certain point to the target direction line in a clockwise direction. =90°, , , , , =0°, =270°, =180°, =90°, and is a random number uniformly distributed in the interval [0, 1].

[0038] Step 3: Determine the absorption surface of the energy beam in the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material according to the propagation path of the energy beam. According to the emission position and emission direction of the energy beam calculated in step 2, a positioning point is taken outside the unit cell model of the TPMS-Gyroid rod-shaped structure lattice material along its emission direction, and the line connecting the positioning point and the energy beam emission position is used as the motion trajectory of the energy beam. The imaginary surface intersecting with the line is the absorption surface, which is recorded as the jth absorption surface.

[0039] Step 4: Calculate the shielding of the energy beam according to the distribution of the cores. The energy beam that is not shielded by the core is absorbed by the jth absorption surface. The motion trajectory of the energy beam is regarded as a line segment m, and the core is regarded as a cylinder n. The line segment m is represented by the centroid position vector and the unit vector along the axial direction Indicates that the cylinder n is represented by the center of mass position vector and the unit vector along the axial direction Indicates that the diameter of cylinder n is d. The shortest distance between line segment m and cylinder n The calculation formula is: , (twenty three); in, is perpendicular to and The unit vector of the line segment m; the distance between the point on the axis closest to the cylinder n , the distance between the center of mass of cylinder n and the point on the axis closest to line segment m The calculation formula is: (twenty four); (25); When the condition of formula 26 is met, the line segment m intersects with the cylinder n, that is, the energy beam is blocked by the core: (26); When the energy beam is not blocked, it is considered to be absorbed by the jth absorption surface.

[0040] When the energy beam is blocked, the absorption ratio of the core surface is used to calculate whether the energy beam is absorbed by the core. The calculation formula is: (27); in, is a random number uniformly distributed in the interval [0, 1]. is the absorption ratio of the core surface. When formula 27 is satisfied, the energy beam is absorbed by the core. When formula 27 is not satisfied, the energy beam is reflected, and we return to step 2, randomly select a new emission position and emission direction from the core surface, and calculate again.

[0041] Step 5: Calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, where the mathematical model is: (28); in, It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

[0042] The embodiment of the present application also provides a radiation transfer coefficient calculation system, which can implement the above-mentioned radiation transfer coefficient calculation method, and the system includes: a modeling module, an emission azimuth calculation module, a propagation path calculation module, a statistical module and a radiation transfer coefficient calculation module; Modeling module, used to establish unit cell model according to the microstructure of ceramic-based lattice materials; An emission direction calculation module is used to use the i-th imaginary surface of the unit cell model as the emission surface of the energy beam, obtain the emission position and emission direction of the energy beam through the Monte Carlo method, and count the number of energy beams emitted from the i-th imaginary surface; A propagation path calculation module is used to determine the absorption surface of the energy beam in the unit cell model according to the propagation path of the energy beam, and the absorption surface is recorded as the jth absorption surface; A statistical module, used for counting the number of energy beams absorbed by the j-th absorption surface according to the shielding of the energy beams by the core; The radiation transfer coefficient calculation module is used to calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, wherein the mathematical model is: ; It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

[0043] It can be understood that the contents of the above method embodiments are all applicable to the present system embodiments, the functions specifically implemented by the present system embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0044] The embodiment of the present application also provides a device for calculating the radiation transfer coefficient, which includes a memory and a processor. The processor can be implemented in the form of a general-purpose CPU (Central Processing Unit), a microprocessor, an application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), or one or more integrated circuits, and the memory can be implemented in the form of a read-only memory (Read Only Memory, ROM), a static storage device, a dynamic storage device, or a random access memory (Random Access Memory, RAM). The memory stores a computer program, and when the processor executes the computer program, any of the above-mentioned radiation transfer coefficient calculation methods is implemented.

[0045] It can be understood that the contents of the above method embodiments are all applicable to the present device embodiments, the functions specifically implemented by the present device embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0046] An embodiment of the present application further provides a computer-readable storage medium, which stores a computer program. When the computer program is executed by a processor, the calculation method of the radiation transfer coefficient is implemented.

[0047] It can be understood that the contents of the above method embodiments are all applicable to the present storage medium embodiments, the functions specifically implemented by the present storage medium embodiments are the same as those of the above method embodiments, and the beneficial effects achieved are also the same as those achieved by the above method embodiments.

[0048] The memory, as a non-transient computer-readable storage medium, can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some embodiments, the memory may optionally include a memory remotely disposed relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0049] The embodiments described in the embodiments of the present application are intended to more clearly illustrate the technical solutions of the embodiments of the present application and do not constitute a limitation on the technical solutions provided in the embodiments of the present application. Those skilled in the art will appreciate that with the evolution of technology and the emergence of new application scenarios, the technical solutions provided in the embodiments of the present application are also applicable to similar technical problems.

[0050] Those skilled in the art will appreciate that the technical solutions shown in the figures do not constitute a limitation on the embodiments of the present application, and may include more or fewer steps than shown in the figures, or a combination of certain steps, or different steps.

[0051] The device embodiments described above are merely illustrative, and the units described as separate components may or may not be physically separated, that is, they may be located in one place or distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0052] Those skilled in the art will appreciate that all or some of the steps in the methods disclosed above, and the functional modules / units in the systems and devices may be implemented as software, firmware, hardware, or a suitable combination thereof.

[0053] The terms "first", "second", "third", "fourth", etc. (if any) in the specification of the present application and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable where appropriate, so that the embodiments of the present application described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device comprising a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0054] It should be understood that in the present application, "at least one (item)" means one or more, and "plurality" means two or more. "And / or" is used to describe the association relationship of associated objects, indicating that three relationships may exist. For example, "A and / or B" can mean: only A exists, only B exists, and A and B exist at the same time, where A and B can be singular or plural. The character " / " generally indicates that the objects associated before and after are in an "or" relationship. "At least one of the following" or similar expressions refers to any combination of these items, including any combination of single or plural items. For example, at least one of a, b or c can mean: a, b, c, "a and b", "a and c", "b and c", or "a and b and c", where a, b, c can be single or multiple.

[0055] In the several embodiments provided in the present application, it should be understood that the disclosed methods, systems and devices can be implemented in other ways. For example, the device embodiments described above are only schematic. For example, the division of the above units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0056] The units described above as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0057] In addition, each functional unit in each embodiment of the present application may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit. The above-mentioned integrated unit may be implemented in the form of hardware or in the form of software functional units.

[0058] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product, and the computer software product is stored in a storage medium, including multiple instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods of various embodiments of the present application. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (Read-Only Memory, referred to as ROM), random access memory (Random Access Memory, referred to as RAM), disk or optical disk and other media that can store programs.

[0059] Although the description of the present application has been quite detailed and specifically describes several described embodiments, it is not intended to be limited to any of these details or embodiments or any particular embodiment, but should be regarded as providing a broad possible interpretation of these claims by reference to the attached claims, taking into account the prior art, so as to effectively cover the intended scope of the present application. In addition, the above description of the present application is based on the embodiments foreseeable by the inventor, and its purpose is to provide a useful description, and those non-substantial changes to the present application that have not yet been foreseen may still represent equivalent changes to the present application.

Claims

1. A method for calculating a radiation transfer coefficient, characterized in that: include: Step 1: Establish a unit cell model according to the microscopic configuration of the ceramic-based lattice material; Step 2: using the i-th imaginary surface of the unit cell model as the emission surface of the energy beam, obtaining the emission position and emission direction of the energy beam by Monte Carlo method, and counting the number of energy beams emitted from the i-th imaginary surface; Step 3, determining the absorption surface of the energy beam in the unit cell model according to the propagation path of the energy beam, and the absorption surface is recorded as the jth absorption surface; Step 4, counting the number of energy beams absorbed by the j-th absorption surface according to the shielding of the energy beams by the core; Step 5: Calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, where the mathematical model is: ; It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

2. A method for calculating the radiation transfer coefficient according to claim 1, characterized in that: The probability distribution model of the emission position of the energy beam in step 2 is: ； ； ； Wherein, x, y, z are the coordinates of the emission point of the energy beam, , , is the maximum coordinate point that can be taken in the area occupied by the unit cell model, , , is the minimum coordinate point that can be taken in the area occupied by the unit cell model, , , is a random number uniformly distributed in the interval [0, 1].

3. The method for calculating the radiation transfer coefficient according to claim 1, characterized in that: The probability distribution model of the emission direction of the energy beam in step 2 is: ； ； ； ； in, is the minimum zenith angle that each imaginary surface in the unit cell model can take, is the minimum azimuth angle that each imaginary surface in the unit cell model can take, =90°, , , , , =0°, =270°, =180°, =90°, and is a random number uniformly distributed in the interval [0, 1].

4. The method for calculating the radiation transfer coefficient according to claim 1, characterized in that: The method for calculating the propagation path of the energy beam in step 3 is: through the energy beam emission position and emission direction calculated in step 2, a positioning point is taken outside the unit cell model along its emission direction, and the line connecting the positioning point and the emission position of the energy beam is used as the movement trajectory of the energy beam, and the surface of the unit cell model intersecting with the line is the absorption surface.

5. The method for calculating the radiation transfer coefficient according to claim 1, characterized in that: The method for calculating the energy beam being blocked by the core in step 4 is: the motion trajectory of the energy beam is regarded as a line segment m, the core is regarded as a cylinder n, and the line segment m is represented by the centroid position vector and the unit vector along the axial direction Indicates that the cylinder n is represented by the center of mass position vector and the unit vector along the axial direction Indicates that the diameter of the cylinder n is d; the shortest distance between the line segment m and the cylinder n is The calculation formula is: ， ； in, is perpendicular to and The unit vector of the line segment m; the distance between the point on the axis closest to the cylinder n , the distance between the center of mass of the cylinder n and the point on the axis closest to the line segment m The calculation formula is: ； ； When the following conditions are met, the line segment m intersects with the cylinder n, that is, the energy beam is blocked by the core: 。 6. A method for calculating the radiation transfer coefficient according to claim 5, characterized in that: When the energy beam is not blocked by the core, it is considered to be absorbed by the jth absorption surface. When the energy beam is blocked by the core, it is necessary to calculate whether the energy beam is absorbed by the core. The energy beam is calculated by the absorption ratio of the core surface. The calculation formula is: ； in, is a random number uniformly distributed in the interval [0, 1]. is the absorption ratio of the core surface.

7. A method for calculating the radiation transfer coefficient according to claim 6, characterized in that: When the energy beam is not absorbed by the core, it is regarded as reflected, and the calculation is returned to step 2.

8. A radiation transfer coefficient calculation system, characterized in that: include: Modeling module, emission azimuth calculation module, propagation path calculation module, statistics module and radiation transfer coefficient calculation module; Modeling module, used to establish unit cell model according to the microstructure of ceramic-based lattice materials; An emission direction calculation module, used for taking the i-th imaginary surface of the unit cell model as the emission surface of the energy beam, obtaining the emission position and emission direction of the energy beam by Monte Carlo method, and counting the number of energy beams emitted from the i-th imaginary surface; A propagation path calculation module, used to determine the absorption surface of the energy beam in the unit cell model according to the propagation path of the energy beam, wherein the absorption surface is recorded as the jth absorption surface; A statistical module, used for counting the number of energy beams absorbed by the j-th absorption surface according to the shielding of the energy beams by the core; The radiation transfer coefficient calculation module is used to calculate the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface through a mathematical model, wherein the mathematical model is: ; It is expressed as the radiation transfer coefficient between the i-th imaginary surface and the j-th absorption surface, N is expressed as the number of energy beams emitted from the i-th imaginary surface, and M is expressed as the number of energy beams absorbed by the j-th absorption surface.

9. A radiation transfer coefficient calculation device, characterized in that: include: A memory and a processor, wherein the memory stores a computer program, and the processor implements the method for calculating the radiation transfer coefficient according to any one of claims 1 to 7 when executing the computer program.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the method for calculating the radiation transfer coefficient according to any one of claims 1 to 7 is implemented.

Images