A method, device, equipment and storage medium for calculating diffuse angle coefficient in high temperature engineering field

Abstract

The present application discloses a method, device, equipment and storage medium for calculating a diffusion angle coefficient, which relates to the field of high temperature engineering. The method comprises: a processing device determines multiple radiation surfaces of a target space, selects an sth radiation surface and a kth radiation surface from the multiple radiation surfaces, determines an hth radiation surface from the multiple radiation surfaces, the hth radiation surface, the sth radiation surface and the kth radiation surface are different radiation surfaces, determines the occlusion between the sth radiation surface and the kth radiation surface by judging the positional relationship between the sth radiation surface, the kth radiation surface and the hth radiation surface, and calculates the diffusion angle coefficient according to the occlusion. The method can reduce the amount of calculation and improve the calculation efficiency.

Description

Technical Field

[0001] The present application relates to the field of high temperature engineering, and in particular to a method, device, equipment and storage medium for calculating a diffusion angle coefficient. Background Art

[0002] In the field of high-temperature engineering, especially in various industrial furnaces such as high-temperature combustion chambers and nuclear reactors, radiation heat transfer is one of the important heat transfer methods. When medium radiation is ignored, only radiation heat transfer between surfaces needs to be considered. In practical applications, accurate calculation of the diffuse angle coefficient between surfaces is crucial to solving surface thermal radiation transmission. The diffuse angle coefficient describes the degree of direct radiation heat exchange between two surfaces, which depends not only on the geometric shape and relative position of the surfaces, but also on the shielding effect.

[0003] With the increasing complexity of industrial design, traditional analytical methods or simplified assumptions can no longer meet actual engineering needs, especially when complex geometric shapes are involved. Numerical simulation methods have been widely used in surface radiation calculations due to their ability to handle arbitrarily complex shapes and shading effects.

[0004] However, existing methods for calculating the diffuse angle coefficient usually require dividing the surface polygonal mesh into triangular meshes when considering occlusion between surfaces. However, such methods have a large amount of processing and low computational efficiency when dealing with systems containing a large number of surfaces. Summary of the invention

[0005] The present application provides a method for calculating a diffusion angle coefficient, which can reduce the amount of calculation and improve the calculation efficiency when processing a system containing a large number of surfaces.

[0006] In order to achieve the above objectives, this application adopts the following technical solutions:

[0007] In a first aspect, the present application provides a method for calculating a diffusion angle coefficient applicable to a polygonal mesh, the method comprising:

[0008] Determine a plurality of radiation surfaces of the target space, wherein the radiation surfaces are in the shape of polygonal meshes;

[0009] Determine whether an s-th radiation surface and a k-th radiation surface among the plurality of radiation surfaces face each other, and obtain a first determination result;

[0010] If the first judgment result indicates that the s th radiation surface and the k th radiation surface are facing each other, determining the h th radiation surface from a plurality of radiation surfaces, and the h th radiation surface, the s th radiation surface and the k th radiation surface are different radiation surfaces;

[0011] Determine whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, and obtain a second determination result;

[0012] If the second judgment result indicates that the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, it is determined whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located, and a third judgment result is obtained, where the first center line is a line connecting the center point of the sth radiation surface and the center point of the kth radiation surface;

[0013] If the third judgment result indicates that the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, determine whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and obtain a fourth judgment result;

[0014] If the fourth judgment result indicates that the first intersection is located between the s th radiation surface and the k th radiation surface, determine whether the first intersection is located inside the h th radiation surface to obtain a fifth judgment result;

[0015] If the fifth judgment result indicates that the first intersection point is located inside the h th radiation surface, it is determined that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero.

[0016] In some possible implementations, the hth radiation surface is an m-polygonal grid, and the determining whether the first intersection point is located inside the hth radiation surface to obtain a fifth determination result includes:

[0017]

[0018] if to are both greater than or equal to 0, then a fifth judgment result is obtained that the first intersection is located inside the hth radiation surface; if to If any item of is less than 0, a fifth judgment result is obtained that the first intersection point is not located inside the hth radiation surface;

[0019] in, is the qth edge and the first intersection point The relationship indicator, Greater than or equal to 0, indicating the first intersection Located inside the qth edge, is the starting point of the qth edge, is the end point of the qth edge, is the starting point of the mth edge, is the end point of the mth edge, represents the normal vector of the h-th radiation surface, and q is an integer.

[0020] In some possible implementations, the determining whether a first intersection point of a plane where the hth radiation surface is located and a straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface to obtain a fourth determination result includes:

[0021]

[0022] If t is less than 0 or t is greater than 1, a fourth judgment result is obtained that the first intersection is not located between the s th radiation surface and the k th radiation surface; if t is greater than 0 and t is less than 1, a fourth judgment result is obtained that the first intersection is located between the s th radiation surface and the k th radiation surface;

[0023] Among them, t is the relationship index between the vector corresponding to the first center line and the first intersection point, represents the center point of the s-th radiating surface, represents any point within the hth radiation surface, represents the center point of the k-th radiating surface, Represents the normal vector of the h-th radiating surface.

[0024] In some possible implementations, the determining whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface to obtain a second determination result includes:

[0025] if or , then a second judgment result is obtained that the hth radiation surface is not within the radiation area of ​​the sth radiation surface or the kth radiation surface;

[0026] if and , then a second judgment result of the hth radiation surface being within the radiation area of ​​the sth radiation surface or the kth radiation surface is obtained;

[0027] in, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the h-th radiation surface, is the center point of the kth radiation surface.

[0028] In some possible implementations, the determining whether the s th radiation surface and the k th radiation surface among the plurality of radiation surfaces face each other to obtain a first determination result includes:

[0029] if or 0, then the first judgment result is obtained that the s-th radiation surface and the k-th radiation surface are not facing each other;

[0030] if and 0, then the first judgment result that the s-th radiation surface and the k-th radiation surface are facing each other is obtained;

[0031] in, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the kth radiation surface.

[0032] In some possible implementations, the method further includes:

[0033] If the first judgment result indicates that the s th radiation surface and the k th radiation surface are not facing each other, it is determined that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero.

[0034] In some possible implementations, the method further includes:

[0035] If the fifth judgment result indicates that the first intersection point is not located inside the hth radiation surface, the diffusion angle coefficient of the sth radiation surface to the kth radiation surface is calculated according to the following formula:

[0036]

[0037] in, represents the angle between the normal vector of the sth radiation surface and the line connecting the first center, represents the angle between the normal vector of the kth radiation surface and the line connecting the first center, is the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface, is the length of the line segment from the center point of the s th radiation surface to the center point of the k th radiation surface, represents the area of ​​the kth radiating surface.

[0038] In a second aspect, the present application provides a device for calculating a diffusion angle coefficient, the device comprising:

[0039] A determination module, used for determining a plurality of radiation surfaces of a target space, wherein the shape of the radiation surfaces is a polygonal mesh;

[0040] A judgment module is used to judge whether an s-th radiation surface and a k-th radiation surface among a plurality of radiation surfaces are facing each other, and obtain a first judgment result; if the first judgment result indicates that the s-th radiation surface and the k-th radiation surface are facing each other, determine an h-th radiation surface from a plurality of radiation surfaces, and the h-th radiation surface, the s-th radiation surface and the k-th radiation surface are different radiation surfaces; judge whether the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, and obtain a second judgment result; if the second judgment result indicates that the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, judge whether the plane where the h-th radiation surface is located is consistent with the plane where the h-th radiation surface is located. A third judgment result is obtained by determining whether the straight line where the first center line is located is parallel to the straight line where the first center line is located, and the first center line is the line connecting the center point of the s th radiation surface and the center point of the k th radiation surface; if the third judgment result indicates that the plane where the h th radiation surface is located is not parallel to the straight line where the first center line is located, it is determined whether the first intersection point of the plane where the h th radiation surface is located and the straight line where the first center line is located is located between the s th radiation surface and the k th radiation surface, and a fourth judgment result is obtained; if the fourth judgment result indicates that the first intersection point is located between the s th radiation surface and the k th radiation surface, it is determined whether the first intersection point is located inside the h th radiation surface, and a fifth judgment result is obtained;

[0041] A calculation module is configured to determine that a diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero if the fifth judgment result indicates that the first intersection point is located inside the h th radiation surface.

[0042] In a third aspect, the present application provides a computing device, including a memory and a processor;

[0043] One or more computer programs are stored in the memory, and the one or more computer programs include instructions; when the instructions are executed by the processor, the computing device executes the method as described in any one of the first aspects.

[0044] In a fourth aspect, the present application provides a computer-readable storage medium, wherein the computer-readable storage medium is used to store a computer program, and the computer program is used to execute the method as described in any one of the first aspects.

[0045] It can be seen from the above technical solution that the present application has at least the following beneficial effects:

[0046] In the present application, the processing device determines multiple radiation surfaces of the target space, wherein the shapes of the multiple radiation surfaces are polygonal meshes, selects the sth radiation surface and the kth radiation surface from the multiple radiation surfaces, and determines whether the sth radiation surface and the kth radiation surface are facing each other. If the sth radiation surface and the kth radiation surface are facing each other, it indicates that there may be occlusion between the sth radiation surface and the kth radiation surface, and it is necessary to further determine whether there is occlusion. Then the processing device determines the hth radiation surface from multiple radiation surfaces, where the hth radiation surface, the sth radiation surface and the kth radiation surface are different radiation surfaces, and determines whether the hth radiation surface is in the radiation area of ​​the sth radiation surface or the kth radiation surface. If so, it indicates that the hth radiation surface may be blocked between the sth radiation surface and the kth radiation surface, and then determines whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located. If the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, it indicates that the hth radiation surface intersects with the straight line where the first center line is located, and it is necessary to continue to determine whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and then determine whether the first intersection point is located inside the hth radiation surface to determine the shielding situation between the sth radiation surface and the kth radiation surface, and calculate the diffusion angle coefficient according to the shielding situation. If the hth radiation surface is blocked between the sth radiation surface and the kth radiation surface, the diffusion angle coefficient between the sth radiation surface and the kth radiation surface is determined to be zero. The traditional solution is to divide the polygonal mesh into multiple triangular meshes and calculate the diffuse angle coefficient on each triangular mesh. This method has a large amount of calculation in the calculation process and has low calculation efficiency. It can be seen that in this application, the occlusion situation can be directly determined by determining multiple radiation surfaces in the target space, and then the diffuse angle coefficient of the polygonal mesh can be calculated based on the judgment results. It can be seen that the solution of this application can reduce the amount of calculation and improve the calculation efficiency.

[0047] It should be understood that the description of technical features, technical solutions, beneficial effects or similar language in this application does not imply that all features and advantages can be realized in any single embodiment. On the contrary, it is understood that the description of features or beneficial effects means that specific technical features, technical solutions or beneficial effects are included in at least one embodiment. Therefore, the description of technical features, technical solutions or beneficial effects in this specification does not necessarily refer to the same embodiment. Furthermore, the technical features, technical solutions and beneficial effects described in the present embodiment can also be combined in any appropriate manner. Those skilled in the art will understand that the embodiment can be realized without one or more specific technical features, technical solutions or beneficial effects of a specific embodiment. In other embodiments, additional technical features and beneficial effects can also be identified in a specific embodiment that does not embody all embodiments. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] Figure 1 A flow chart of a method for calculating a diffusion angle coefficient provided in an embodiment of the present application;

[0049] Figure 2 A schematic diagram of a first judgment result provided in an embodiment of the present application;

[0050] Figure 3 A schematic diagram of a second first judgment result provided in an embodiment of the present application;

[0051] Figure 4 A schematic diagram of a first second judgment result provided in an embodiment of the present application;

[0052] Figure 5 A schematic diagram of a second second judgment result provided in an embodiment of the present application;

[0053] Figure 6 A schematic diagram of a third second judgment result provided in an embodiment of the present application;

[0054] Figure 7 A schematic diagram of the first third judgment result provided in an embodiment of the present application;

[0055] Figure 8 A schematic diagram of the second third judgment result provided in an embodiment of the present application;

[0056] Fig. 9 A schematic diagram of the first fourth judgment result provided in an embodiment of the present application;

[0057] Fig.10 A schematic diagram of the second fourth judgment result provided in an embodiment of the present application;

[0058] Fig.11 A schematic diagram of the third and fourth judgment results provided in the embodiment of the present application;

[0059] Fig.12 A schematic diagram of a fifth judgment result provided in an embodiment of the present application;

[0060] Fig.13 A schematic diagram of an h-th radiation surface provided in an embodiment of the present application;

[0061] Fig.14 An example diagram of a polygonal mesh provided in an embodiment of the present application;

[0062] Fig.15 A schematic diagram of calculating a diffusion angle coefficient provided in an embodiment of the present application;

[0063] Fig.16A schematic diagram of a device for calculating a diffusion angle coefficient provided in an embodiment of the present application;

[0064] Fig.17 A schematic diagram of a computing device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0065] The terms "first", "second", "third", etc. in the specification of this application and the accompanying drawings are used to distinguish different objects rather than to limit a specific order.

[0066] In the embodiments of the present application, words such as "exemplary" or "for example" are used to indicate examples, illustrations or descriptions. Any embodiment or design described as "exemplary" or "for example" in the embodiments of the present application should not be interpreted as being more preferred or more advantageous than other embodiments or designs. Specifically, the use of words such as "exemplary" or "for example" is intended to present related concepts in a specific way.

[0067] In order to make the description of the following embodiments clear and concise, a brief introduction to the related technology is first given:

[0068] Radiation is a phenomenon of energy transmission in the form of electromagnetic waves. Objects can emit radiation energy for various reasons. The process of emitting radiation energy due to heat is called thermal radiation. Thermal radiation is emitted in the form of electromagnetic waves and propagates in space. When it encounters another object, it is partially or completely absorbed and converted into thermal energy again. Radiation is not only the transfer of energy but also accompanied by the transformation of energy forms. The most important thermal radiation in industry is the mutual radiation between solids, and radiation can only become the main way of heat transfer at high temperatures.

[0069] The most important thing in surface thermal radiation calculation in the field of high temperature engineering is to accurately calculate the diffusion angle coefficient between surfaces. The diffusion angle coefficient refers to the share of the radiation energy emitted from one solid surface that reaches another solid surface. The diffusion angle coefficient is a coefficient that reflects the geometric shape and position relationship between different objects that radiate each other. The diffusion angle coefficient is a physical quantity that describes the relationship between radiation heat transfer between surfaces and factors such as the surface geometry and the relative position of each surface. The concept of the diffusion angle coefficient is of great significance in the field of thermal radiation heat transfer. The diffusion angle coefficient can help understand and calculate the radiation heat exchange process between objects.

[0070] At present, the traditional solution is to divide the polygonal mesh into multiple triangular meshes and calculate the diffuse angle coefficient on each triangular mesh. This method has a large amount of calculation in the calculation process and has low calculation efficiency.

[0071] In view of this, an embodiment of the present application provides a method for calculating a diffusion angle coefficient, in which a processing device determines multiple radiation surfaces of a target space, wherein the shapes of the multiple radiation surfaces are polygonal meshes, selects an s th radiation surface and a k th radiation surface from the multiple radiation surfaces, and determines whether the s th radiation surface and the k th radiation surface are facing each other. If the s th radiation surface and the k th radiation surface are facing each other, it indicates that there may be occlusion between the s th radiation surface and the k th radiation surface, and it is necessary to further determine whether there is occlusion. Then the processing device determines the hth radiation surface from multiple radiation surfaces, where the hth radiation surface, the sth radiation surface and the kth radiation surface are different radiation surfaces, and determines whether the hth radiation surface is in the radiation area of ​​the sth radiation surface or the kth radiation surface. If so, the hth radiation surface may be blocked between the sth radiation surface and the kth radiation surface, and then determines whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located. If the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, it indicates that the hth radiation surface intersects with the straight line where the first center line is located, and it is necessary to continue to determine whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and then determine whether the first intersection point is located inside the hth radiation surface to determine the shielding situation between the sth radiation surface and the kth radiation surface, and calculate the diffusion angle coefficient according to the shielding situation. If the hth radiation surface is blocked between the sth radiation surface and the kth radiation surface, the diffusion angle coefficient between the sth radiation surface and the kth radiation surface is determined to be zero. The traditional solution is to divide the polygonal mesh into multiple triangular meshes and calculate the diffuse angle coefficient on each triangular mesh. This method has a large amount of calculation in the calculation process and has low calculation efficiency. It can be seen that in this application, the occlusion situation can be directly determined by determining multiple radiation surfaces in the target space, and then the diffuse angle coefficient of the polygonal mesh can be calculated based on the judgment results. It can be seen that the solution of this application can reduce the amount of calculation and improve the calculation efficiency.

[0072] In order to make the technical solution of the present application clearer and easier to understand, a calculation method for the diffusion angle coefficient of a polygonal mesh provided in an embodiment of the present application is introduced. Figure 1 As shown in FIG. 1 , the figure is a flow chart of a method for calculating a diffusion angle coefficient provided in an embodiment of the present application. The method for calculating the diffusion angle coefficient of the polygonal mesh can be executed by a processing device, and the method includes:

[0073] S101. A processing device determines a plurality of radiation surfaces of a target space, where the radiation surfaces are in the shape of polygonal meshes.

[0074] Processing equipment refers to hardware or software systems used to perform specific tasks or process data. In scenarios involving radiation surface analysis, processing equipment may refer to computers, servers, or specific software platforms dedicated to calculating, simulating, or analyzing radiation characteristics.

[0075] The target space refers to a specific area or environment that needs to be analyzed, simulated or processed. For example, the target space can be a room or a building.

[0076] A radiating surface is a surface that can emit, absorb or reflect radiant energy. In physics and engineering, these surfaces can be walls and windows in a room or the exterior walls and roofs of buildings.

[0077] It should be noted that the present application does not limit the method for determining the multiple radiation surfaces of the target space, and personnel can choose the method for determining the multiple radiation surfaces based on actual needs. For example, when the processing device determines the multiple radiation surfaces, a grid calculated based on CFD (Computational Fluid Dynamics) can be determined as the radiation surface.

[0078] S102: The processing device determines whether the diffusion angle coefficient of the kth radiation surface to the sth radiation surface is known.

[0079] If the processing device determines that the diffusion angle coefficient of the kth radiation surface to the sth radiation surface is known, in this case S103 is executed.

[0080] If the processing device determines that the diffusion angle coefficient of the kth radiation surface with respect to the sth radiation surface is not known, in this case S104 is performed.

[0081] S103, calculating the diffusion angle coefficient of the s th radiation surface to the k th radiation surface.

[0082] Since the diffusion angle coefficient of the kth radiation surface to the sth radiation surface is known, the diffusion angle coefficient of the sth radiation surface to the kth radiation surface can be calculated by the following formula.

[0083]

[0084] in, is the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface, is the diffusion angle coefficient of the k-th radiation surface to the k-th radiation surface, represents the area of ​​the kth radiating surface, represents the area of ​​the sth radiating surface.

[0085] S104: The processing device determines whether the s th radiation surface and the k th radiation surface among the plurality of radiation surfaces face each other.

[0086] The processing device determines whether the s-th radiation surface and the k-th radiation surface among the multiple radiation surfaces are facing each other, and obtains a first judgment result. The first judgment result may be that the s-th radiation surface and the k-th radiation surface among the multiple radiation surfaces are facing each other, in which case S105 is executed, or the s-th radiation surface and the k-th radiation surface among the multiple radiation surfaces are not facing each other, in which case S111 is executed.

[0087] like Figure 2 As shown, this figure is a schematic diagram of the first judgment result provided in the embodiment of the present application, Figure 3 This is a schematic diagram of a second first judgment result provided in an embodiment of the present application, wherein the radiation direction of the s th radiation surface 201 is the direction indicated by the first arrow 202 , and the radiation direction of the k th radiation surface 203 is the direction indicated by the second arrow 204 .

[0088] Determine the normal vector of the s-th radiation surface 201 , center point S, normal vector of the kth radiation surface 203 , center point K, and determine whether the s-th radiation surface and the k-th radiation surface are facing each other based on the positional relationship between the normal vector and the line connecting the first center point.

[0089] Specifically, the following formula is used for judgment:

[0090] if or 0, then the first judgment result that the s-th radiation surface 201 and the k-th radiation surface 203 are not facing each other is obtained, such as Figure 2 As shown, it means that the s-th radiation surface 201 and the k-th radiation surface 203 are not in the radiation area of ​​each other, and the diffusion angle coefficient is zero.

[0091] if and 0, then the first judgment result that the s-th radiation surface 201 and the k-th radiation surface 203 face each other is obtained, such as Figure 3 As shown, it means that the s th radiation surface 201 and the k th radiation surface 203 are in the radiation area of ​​each other, and it is necessary to further determine whether there is any shielding between the s th radiation surface 201 and the k th radiation surface 203 .

[0092] in, is the center point of the s-th radiation surface 201, is the center point of the k-th radiation surface 203, is the vector from the center point S of the s-th radiation surface 201 to the center point K of the k-th radiation surface 203. The center point S of the s-th radiation surface 201 is The center point K of the kth radiation surface 203 is The end point, is a vector from the center point K of the k-th radiation surface 203 to the center point S of the s-th radiation surface 201. The center point K of the k-th radiation surface 203 is The starting point of the s radiation surface 201 is the center point S The end point.

[0093] S105. The processing device determines an hth radiation surface from a plurality of radiation surfaces.

[0094] The hth radiation surface 301, the sth radiation surface 201, and the kth radiation surface 203 are different radiation surfaces. In some examples, the plurality of radiation surfaces include the 1st radiation surface, the 2nd radiation surface ... the 10th radiation surface, wherein the sth radiation surface 201 may be the 1st radiation surface, the kth radiation surface 203 may be the 2nd radiation surface, and the hth radiation surface 301 may be one of the 3rd to 10th radiation surfaces.

[0095] S106. The processing device determines whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface.

[0096] The processing device determines whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, and obtains a second judgment result. The second judgment result may be that the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, in which case S107 is executed, or that the hth radiation surface is not within the radiation area of ​​the sth radiation surface or the kth radiation surface, in which case S110 is executed.

[0097] like Figure 4 As shown, this figure is a schematic diagram of the first second judgment result provided in an embodiment of the present application, Figure 5 A schematic diagram of a second second judgment result provided in an embodiment of the present application, Figure 6 A schematic diagram of a third second judgment result provided in an embodiment of the present application determines the hth radiation surface 301 from multiple radiation surfaces, wherein the hth radiation surface 301, the sth radiation surface 201 and the kth radiation surface 203 are different radiation surfaces.

[0098] Determine the center point H of the hth radiation surface 301, and according to the normal vectors of the hth radiation surface 301 and the sth radiation surface 201 The positional relationship between the h-th radiation surface 301 and the k-th radiation surface 203 The positional relationship is used to determine whether the h-th radiation surface 301 is within the radiation area of ​​the s-th radiation surface 201 or the k-th radiation surface 203.

[0099] Specifically, the following formula is used for judgment:

[0100] if or , then a second judgment result is obtained that the hth radiation surface 301 is not within the radiation area of ​​the sth radiation surface 201 or the kth radiation surface 203, such as Figure 4 and Figure 5 As shown, it means that there is no shielding of the hth radiation surface 301 between the sth radiation surface 201 and the kth radiation surface 203, and the next group of grids is used to judge the shielding between the sth radiation surface and the kth radiation surface.

[0101] if and , then a second judgment result is obtained that the hth radiation surface 301 is within the radiation area of ​​the sth radiation surface 201 or the kth radiation surface 203, such as Figure 6 As shown, the location of the hth radiation surface 301 needs to be further determined.

[0102] in, is the center point of the s-th radiation surface 201, is the center point of the k-th radiation surface 203, H is the center point of the h-th radiation surface 301, is the vector from the center point S of the s-th radiation surface 201 to the center point H of the h-th radiation surface 301. The center point S of the s-th radiation surface 201 is The starting point of the hth radiation surface 301 is The end point, is a vector from the center point K of the k-th radiation surface 203 to the center point H of the h-th radiation surface 301. The center point K of the k-th radiation surface 203 is The starting point of the hth radiation surface 301 is The end point.

[0103] S107. The processing device determines whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located.

[0104] The first center line is a line segment connecting the center point S of the s th radiation surface and the center point K of the k th radiation surface. The processing device determines whether the plane where the h th radiation surface is located is parallel to the straight line where the first center line is located, and obtains a third judgment result. The third judgment result may be that the plane where the h th radiation surface is located is not parallel to the straight line where the first center line is located, in which case S108 is executed, or the plane where the h th radiation surface is located is parallel to the straight line where the first center line is located, in which case S110 is executed.

[0105] like Figure 7 As shown, this figure is a schematic diagram of the first third judgment result provided in the embodiment of the present application, Figure 8A schematic diagram of the second third judgment result provided in an embodiment of the present application.

[0106] Determine the normal vector of the hth radiation surface 301 , according to the normal vector of the hth radiation surface 301 The positional relationship with the first center line is used to determine whether the plane where the hth radiation surface 301 is located is parallel to the straight line where the first center line is located.

[0107] Specifically, the following formula is used for judgment:

[0108] if , then a third judgment result is obtained that the plane where the hth radiation surface 301 is located is parallel to the straight line where the first center line is located, such as Figure 7 As shown, it means that there is no shielding of the hth radiation surface 301 between the sth radiation surface 201 and the kth radiation surface 203, and the next group of grids is used to judge the shielding between the sth radiation surface and the kth radiation surface.

[0109] if , then a third judgment result is obtained that the plane where the hth radiation surface 301 is located is not parallel to the straight line where the first center line is located, such as Figure 8 As shown, the location of the hth radiation surface 301 needs to be further determined.

[0110] S108. The processing device determines whether a first intersection point of a plane where the hth radiation surface is located and a straight line where the first center connecting line is located is located between the sth radiation surface and the kth radiation surface.

[0111] The processing device determines whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and obtains a fourth judgment result. The fourth judgment result can be that the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, in which case S109 is executed, or the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is not located between the sth radiation surface and the kth radiation surface, in which case S110 is executed.

[0112] like Fig. 9 As shown, this figure is a schematic diagram of the first fourth judgment result provided in the embodiment of the present application, Fig.10 A schematic diagram of the second fourth judgment result provided in an embodiment of the present application, Fig.11 A schematic diagram of the third and fourth judgment results provided in an embodiment of the present application, wherein the plane where the h-th radiation surface 301 is located intersects with the straight line where the first center line is located to obtain a first intersection point P; the position of the h-th radiation surface 301 is judged based on the positional relationship of the first intersection point P with respect to the straight line where the first center line is located.

[0113] Specifically, the following formula is used for judgment:

[0114]

[0115]

[0116] If t is less than 0 or t is greater than 1, a fourth judgment result is obtained that the first intersection is not located between the s th radiation surface 201 and the k th radiation surface 203, such as Fig. 9 and Fig.10 As shown, it means that there is no shielding of the hth radiation surface 301 between the sth radiation surface 201 and the kth radiation surface 203, and the next group of grids is used to judge the shielding between the sth radiation surface and the kth radiation surface.

[0117] If t is greater than 0 and t is less than 1, a fourth judgment result is obtained in which the first intersection is located between the s th radiation surface 201 and the k th radiation surface 203, such as Fig.11 As shown, it is necessary to further determine the location of the hth radiation surface 301;

[0118] Among them, O is the coordinate origin, is the vector from the coordinate origin O to the first intersection point P, is the vector from the coordinate origin O to the center point S of the s-th radiation surface 201, is the coordinate origin O to point h on the hth radiation surface 301 a t is the relationship index between the vector corresponding to the first center line and the first intersection point. represents the center point of the s-th radiation surface 201, represents any point within the hth radiation surface 301, represents the center point of the k-th radiation surface 203, The distance from the center point S of the s-th radiation surface 201 to the point h of the h-th radiation surface 301 a The center point S of the s radiation surface 201 is The starting point of the hth radiation surface 301 is point h a yes The end point, represents the normal vector of the h-th radiation surface 301 .

[0119] S109: The processing device determines whether the first intersection point is located inside the h th radiation surface.

[0120] The processing device determines whether the first intersection is located inside the hth radiation surface to obtain a fifth determination result. The fifth determination result may be that the first intersection is located inside the hth radiation surface, in which case S111 is executed, or that the first intersection is not located inside the hth radiation surface, in which case S110 is executed.

[0121] like Fig.12 As shown in FIG. 1 , this figure is a schematic diagram of a fifth judgment result provided by an embodiment of the present application. In order to more clearly indicate whether the h-th radiation surface 301 is blocked between the s-th radiation surface 201 and the k-th radiation surface 203, the present application flattens the h-th radiation surface 301, as shown in FIG. Fig.13 As shown, this figure is a schematic diagram of an h-th radiation surface provided in an embodiment of the present application.

[0122] It can be seen in the figure that m points are randomly selected on the hth radiation surface 301 and connected end to end. , according to the first intersection point P and , it is determined whether the first intersection is located inside the h th radiation surface 301 , and further whether there is any shielding between the s th radiation surface 201 and the k th radiation surface 203 .

[0123] Specifically, the following formula is used for judgment:

[0124]

[0125] if to are greater than or equal to 0, then the fifth judgment result is obtained that the first intersection is located inside the hth radiation surface 301, indicating that there is an obstruction between the sth radiation surface 201 and the kth radiation surface 203, and the diffusion angle coefficient is zero. If to If any item in is less than 0, the fifth judgment result is obtained that the first intersection is not located inside the hth radiation surface 301, indicating that there is no occlusion between the sth radiation surface 201 and the kth radiation surface 203 by the hth radiation surface, and the next group of grids is used to judge the occlusion between the sth radiation surface and the kth radiation surface.

[0126] in, is the qth edge and the first intersection point The relationship indicator, Greater than or equal to 0, indicating the first intersection Located on the inside of the qth edge, for arrive The vector of for To the first intersection The vector of for arrive The vector of for To the first intersection The vector of is the starting point of the qth edge, is the end point of the qth edge, is the starting point of the mth edge, is the end point of the mth edge, and q is an integer.

[0127] S110, traversal completed.

[0128] If the traversal of all radiation surfaces is completed, S112 is executed in this case; if the traversal of all radiation surfaces is not completed, S105 is executed in this case.

[0129] S111. The diffusion angle coefficient of the s th radiation surface to the k th radiation surface is zero.

[0130] The processing device determines the shielding situation between the s-th radiation surface and the k-th radiation surface based on S104~S109, and concludes that there is shielding of the h-th radiation surface between the s-th radiation surface and the k-th radiation surface, so the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface is zero.

[0131] S112. Calculate the diffusion angle coefficient of the s th radiation surface to the k th radiation surface.

[0132] The processing device determines the shielding situation between the s-th radiation surface and the k-th radiation surface based on S104~S109, and concludes that there is no shielding of the h-th radiation surface between the s-th radiation surface and the k-th radiation surface, and therefore calculates the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface.

[0133] The diffusion angle coefficient of the sth radiation surface to the kth radiation surface is calculated according to the following formula:

[0134]

[0135] in, represents the angle between the normal vector of the sth radiation surface and the line connecting the first center, represents the angle between the normal vector of the kth radiation surface and the line connecting the first center, is the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface, is the length of the line segment from the center point of the s th radiation surface to the center point of the k th radiation surface, represents the area of ​​the kth radiating surface.

[0136] Based on the above description, an embodiment of the present application provides a method for calculating a diffuse angle coefficient suitable for a polygonal mesh, in which a processing device determines multiple radiation surfaces of a target space, wherein the shapes of the multiple radiation surfaces are polygonal meshes, selects an s th radiation surface and a k th radiation surface from the multiple radiation surfaces, and determines whether the s th radiation surface and the k th radiation surface are facing each other. If the s th radiation surface and the k th radiation surface are facing each other, it indicates that there may be occlusion between the s th radiation surface and the k th radiation surface, and it is necessary to further determine whether there is occlusion. Then the processing device determines the hth radiation surface from multiple radiation surfaces, where the hth radiation surface, the sth radiation surface and the kth radiation surface are different radiation surfaces, and determines whether the hth radiation surface is in the radiation area of ​​the sth radiation surface or the kth radiation surface. If so, it indicates that the hth radiation surface may be blocked between the sth radiation surface and the kth radiation surface, and then determines whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located. If the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, it indicates that the hth radiation surface intersects with the straight line where the first center line is located, and it is necessary to continue to determine whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and then determine whether the first intersection point is located inside the hth radiation surface to determine the shielding situation between the sth radiation surface and the kth radiation surface, and calculate the diffusion angle coefficient according to the shielding situation. If the hth radiation surface is blocked between the sth radiation surface and the kth radiation surface, the diffusion angle coefficient between the sth radiation surface and the kth radiation surface is determined to be zero. The traditional solution is to divide the polygonal mesh into multiple triangular meshes and calculate the diffuse angle coefficient on each triangular mesh. This method has a large amount of calculation in the calculation process and has low calculation efficiency. It can be seen that in this application, the occlusion situation can be directly determined by determining multiple radiation surfaces in the target space, and then the diffuse angle coefficient of the polygonal mesh can be calculated based on the judgment results. It can be seen that the solution of this application can reduce the amount of calculation and improve the calculation efficiency.

[0137] According to the above-mentioned calculation method of the diffusion angle coefficient applicable to polygonal meshes, a test example is provided to verify the effectiveness of the above-mentioned method. Fig.14 As shown, this figure is an example diagram of a polygonal mesh provided by an embodiment of the present application. Fig.14 In the calculation, the diffusion angle coefficient of the i-th grid of the bottom plane of the cavity to other planes of the cavity and a spherical surface in the cavity is calculated, wherein the cavity surface is divided into quadrilateral grids, and the spherical surface is divided into triangular grids and quadrilateral grids. The diffusion angle coefficient of the i-th grid of the bottom plane of the cavity to other planes of the cavity and a spherical surface in the cavity in this test example is calculated according to the above-mentioned calculation method of the diffusion angle coefficient applicable to polygonal grids. The calculation process is as follows:

[0138] The processing device determines multiple polygonal meshes to obtain data information of polygonal meshes, selects the jth mesh, wherein the ith mesh and the jth mesh are different meshes, and determines whether there is occlusion between the ith mesh and the jth mesh according to steps S101 to S111, and calculates the diffusion angle coefficient according to the occlusion situation. It should be noted that the jth mesh and the ith mesh are not facing each other, and the non-occlusion situation can be quickly determined, and the diffusion angle coefficient is zero, and there is no need to traverse the third mesh, which improves the judgment rate. If the jth mesh and the ith mesh are facing each other, it can also be determined whether there is a third mesh between the radiation area of ​​the jth mesh and the ith mesh. If it is not between the radiation area of ​​the jth mesh and the ith mesh, the diffusion angle coefficient is zero, and there is no need to calculate the intersection position of the third mesh and the first center line, which improves the calculation rate, and there is no need to divide the original mesh into triangular meshes to determine the occlusion situation, which reduces the calculation time.

[0139] According to the above test example, we can further obtain the calculation results, such as Fig.15 As shown in FIG. 1 , this figure is a schematic diagram of calculating a diffusion angle coefficient provided in an embodiment of the present application. Fig.15 It can be seen that the area on the spherical surface that is invisible to the i-th grid, the areas on both sides of the cavity blocked by the spherical surface, and the area on the bottom plane of the cavity are all effectively judged to be blocked with the i-th grid, and the diffuse angle coefficient of the area on the spherical surface that is invisible to the i-th grid, the areas on both sides of the cavity blocked by the spherical surface, and the area on the bottom plane of the cavity with the i-th grid is zero. Except for the area on the spherical surface that is invisible to the i-th grid, the areas on both sides of the cavity blocked by the spherical surface, and the area on the bottom plane of the cavity, all other areas are judged to have no occlusion with the i-th grid, and the diffuse angle coefficient can be calculated by the diffuse angle coefficient calculation formula.

[0140] Exemplarily, the single-core calculation time of this test example is about 0.08s; if steps S104 to S106 are not executed, and S107 and subsequent judgments are directly executed, the calculation time is about 0.23s; if steps S104 to S106 are executed, the mesh is triangulated, and then subsequent steps are performed, the calculation time is about 28s; if steps S104 to S106 are not executed, and the mesh is triangulated, and then subsequent steps are performed, the calculation time is about 79s.

[0141] The above calculation results are only an example of calculation speed. In actual situations, the calculation time is affected by the equipment and various factors, and the calculation results will be different. The calculation efficiency is further explained below.

[0142] In this test case, the calculation efficiency of the diffuse angle coefficient calculated by a single core is improved compared to the case where the diffuse angle coefficient is calculated without judging whether the i-th grid and the j-th grid are facing each other and whether the third grid is between the radiation areas of the i-th grid and the j-th grid. times; the single-core calculation of the diffuse angle coefficient is compared with judging whether the i-th grid and the j-th grid are facing each other, whether the third grid is between the radiation area of ​​the i-th grid and the j-th grid, and triangulating the grid to calculate the diffuse angle coefficient. The calculation efficiency is improved by times; the single-core calculation of the diffuse angle coefficient is compared with not judging whether the i-th grid and the j-th grid are facing each other, whether the third grid is between the radiation area of ​​the i-th grid and the j-th grid, and triangulating the grid to calculate the diffuse angle coefficient, the calculation efficiency is improved by Times. The decimal part of the calculation efficiency is rounded to three significant figures.

[0143] The present application also provides a schematic diagram of a device for calculating a diffusion angle coefficient, such as Fig.16 As shown, the figure is a schematic diagram of a calculation device for a diffusion angle coefficient provided in an embodiment of the present application, and the device includes a determination module 901, a judgment module 902 and a calculation module 903.

[0144] A determination module 901 is used to determine a plurality of radiation surfaces of a target space, wherein the shape of the radiation surface is a polygonal mesh;

[0145] The judgment module 902 is used to judge whether the s-th radiation surface and the k-th radiation surface among the multiple radiation surfaces are facing each other, and obtain a first judgment result; if the first judgment result indicates that the s-th radiation surface and the k-th radiation surface are facing each other, determine the h-th radiation surface from the multiple radiation surfaces, and the h-th radiation surface, the s-th radiation surface and the k-th radiation surface are different radiation surfaces; judge whether the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, and obtain a second judgment result; if the second judgment result indicates that the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, judge the plane where the h-th radiation surface is located A third judgment result is obtained by determining whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located, and the first center line is the line connecting the center point of the sth radiation surface and the center point of the kth radiation surface; if the third judgment result indicates that the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, determine whether the first intersection point between the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and obtain a fourth judgment result; if the fourth judgment result indicates that the first intersection point is located between the sth radiation surface and the kth radiation surface, determine whether the first intersection point is located inside the hth radiation surface, and obtain a fifth judgment result;

[0146] The calculation module 903 is configured to determine that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero if the fifth judgment result indicates that the first intersection point is located inside the h th radiation surface.

[0147] Optionally, the h-th radiation surface is an m-polygonal grid, and the judgment module 902 is specifically configured to determine the fifth judgment result in the following manner:

[0148]

[0149] if to are both greater than or equal to 0, then a fifth judgment result is obtained that the first intersection is located inside the hth radiation surface; if to If any item of is less than 0, a fifth judgment result is obtained that the first intersection point is not located inside the hth radiation surface;

[0150] in, is the qth edge and the first intersection point The relationship indicator, Greater than or equal to 0, indicating the first intersection Located inside the qth edge, is the starting point of the qth edge, is the end point of the qth edge, is the starting point of the mth edge, is the end point of the mth edge, represents the normal vector of the h-th radiation surface, and q is an integer.

[0151] Optionally, the judgment module 902 is specifically configured to obtain a fourth judgment result in the following manner:

[0152]

[0153] If t is less than 0 or t is greater than 1, a fourth judgment result is obtained that the first intersection is not located between the s th radiation surface and the k th radiation surface; if t is greater than 0 and t is less than 1, a fourth judgment result is obtained that the first intersection is located between the s th radiation surface and the k th radiation surface;

[0154] Among them, t is the relationship index between the vector corresponding to the first center line and the first intersection point, represents the center point of the s-th radiating surface, represents any point within the hth radiation surface, represents the center point of the k-th radiating surface, Represents the normal vector of the h-th radiating surface.

[0155] Optionally, the judgment module 902 is specifically configured to obtain a second judgment result in the following manner:

[0156] if or , then a second judgment result is obtained that the hth radiation surface is not within the radiation area of ​​the sth radiation surface or the kth radiation surface; if and , then a second judgment result is obtained that the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface; wherein, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the h-th radiation surface, is the center point of the kth radiation surface.

[0157] Optionally, the judgment module 902 is specifically configured to obtain a first judgment result in the following manner:

[0158] if or 0, then the first judgment result is obtained that the s-th radiation surface and the k-th radiation surface are not facing each other;

[0159] if and 0, then the first judgment result that the s-th radiation surface and the k-th radiation surface are facing each other is obtained;

[0160] in, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the kth radiation surface.

[0161] Optionally, the calculation module 903 is further configured to determine that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero if the first judgment result indicates that the s th radiation surface and the k th radiation surface are not facing each other.

[0162] Optionally, the calculation module 903 is further configured to calculate the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface using the following formula:

[0163]

[0164] in, represents the angle between the normal vector of the sth radiation surface and the line connecting the first center, represents the angle between the normal vector of the kth radiation surface and the line connecting the first center, is the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface, is the length of the line segment from the center point of the s th radiation surface to the center point of the k th radiation surface, represents the area of ​​the kth radiating surface.

[0165] The present application also provides a computing device. Fig.17 As shown, this figure is a schematic diagram of a computing device provided in an embodiment of the present application, and the computing device 1000 includes a bus 1001, a processor 1002, a communication interface 1003 and a memory 1004. The processor 1002, the memory 1004 and the communication interface 1003 communicate through the bus 1001.

[0166] The bus 1001 may be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of representation, Fig.17 Only one thick line is used in the diagram, but this does not mean that there is only one bus or only one type of bus.

[0167] The processor 1002 may be any one or more of a central processing unit (CPU), a graphics processing unit (GPU), a microprocessor (MP), or a digital signal processor (DSP).

[0168] The communication interface 1003 is used for communicating with the outside.

[0169] The memory 1004 may include a volatile memory, such as a random access memory (RAM). The memory 1004 may also include a non-volatile memory, such as a read-only memory (ROM), a flash memory, a hard disk drive (HDD), or a solid state drive (SSD).

[0170] The memory 1004 stores executable codes, and the processor 1002 executes the executable codes to perform the above method.

[0171] The embodiment of the present application also provides a computer-readable storage medium. The computer-readable storage medium can be any available medium that can be stored by a computing device or a data storage device such as a data center containing one or more available media. The available medium can be a magnetic medium (such as a floppy disk, a hard disk, a tape), an optical medium (such as a DVD), or a semiconductor medium (such as a solid-state hard disk). The computer-readable storage medium includes instructions that instruct the computing device to execute the above method.

[0172] The embodiment of the present application further provides a computer program product, which includes one or more computer instructions. When the computer instructions are loaded and executed on a computing device, the process or function described in the embodiment of the present application is generated in whole or in part.

[0173] The computer instructions may be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions may be transmitted from one website, computer or data center to another website, computer or data center via wired (e.g., coaxial cable, optical fiber, digital subscriber line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means.

[0174] When the computer program product is executed by a computer, the computer executes any of the above methods. The computer program product may be a software installation package, and when any of the above methods is needed, the computer program product may be downloaded and executed on the computer.

[0175] The descriptions of the processes or structures corresponding to the above-mentioned figures have different emphases. For parts that are not described in detail in a certain process or structure, please refer to the relevant descriptions of other processes or structures.

[0176] The above description is only a specific implementation of the present application, but the protection scope of the present application is not limited thereto. Any changes or substitutions within the technical scope disclosed in the present application should be included in the protection scope of the present application.

Claims

1. A method for calculating the diffusion angle coefficient in the field of high temperature engineering, used to calculate the radiation heat exchange between objects, characterized in that: The method comprises: Determine a plurality of radiation surfaces of the target space, wherein the radiation surfaces are in the shape of polygonal meshes; Determine whether an s-th radiation surface and a k-th radiation surface among the plurality of radiation surfaces face each other, and obtain a first determination result; If the first judgment result indicates that the s th radiation surface and the k th radiation surface are facing each other, determining the h th radiation surface from a plurality of radiation surfaces, and the h th radiation surface, the s th radiation surface and the k th radiation surface are different radiation surfaces; Determine whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, and obtain a second determination result; If the second judgment result indicates that the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface, it is determined whether the plane where the hth radiation surface is located is parallel to the straight line where the first center line is located, and a third judgment result is obtained, where the first center line is a line connecting the center point of the sth radiation surface and the center point of the kth radiation surface; If the third judgment result indicates that the plane where the hth radiation surface is located is not parallel to the straight line where the first center line is located, determine whether the first intersection point of the plane where the hth radiation surface is located and the straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface, and obtain a fourth judgment result; If the fourth judgment result indicates that the first intersection is located between the s th radiation surface and the k th radiation surface, determine whether the first intersection is located inside the h th radiation surface to obtain a fifth judgment result; If the fifth judgment result indicates that the first intersection point is located inside the h th radiation surface, it is determined that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero.

2. The method according to claim 1, characterized in that The h-th radiation surface is an m-polygonal grid, and the determining whether the first intersection point is located inside the h-th radiation surface to obtain a fifth determination result includes: if to are both greater than or equal to 0, then a fifth judgment result is obtained that the first intersection is located inside the hth radiation surface; if to If any item of is less than 0, a fifth judgment result is obtained that the first intersection point is not located inside the hth radiation surface; in, is the qth edge and the first intersection point The relationship indicator, Greater than or equal to 0, indicating the first intersection Located inside the qth edge, is the starting point of the qth edge, is the end point of the qth edge, is the starting point of the mth edge, is the end point of the mth edge, represents the normal vector of the h-th radiation surface, and q is an integer.

3. The method according to claim 1, characterized in that: The determining whether a first intersection point of a plane where the hth radiation surface is located and a straight line where the first center line is located is located between the sth radiation surface and the kth radiation surface to obtain a fourth determination result includes: If t is less than 0 or t is greater than 1, a fourth judgment result is obtained that the first intersection is not located between the s th radiation surface and the k th radiation surface; if t is greater than 0 and t is less than 1, a fourth judgment result is obtained that the first intersection is located between the s th radiation surface and the k th radiation surface; Among them, t is the relationship index between the vector corresponding to the first center line and the first intersection point, represents the center point of the s-th radiation surface, represents any point within the hth radiation surface, represents the center point of the k-th radiating surface, Represents the normal vector of the h-th radiating surface.

4. The method according to claim 1, characterized in that The determining whether the hth radiation surface is within the radiation area of ​​the sth radiation surface or the kth radiation surface to obtain a second determination result includes: if or , then a second judgment result is obtained that the hth radiation surface is not within the radiation area of ​​the sth radiation surface or the kth radiation surface; if and , then a second judgment result of the hth radiation surface being within the radiation area of ​​the sth radiation surface or the kth radiation surface is obtained; in, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the h-th radiation surface, is the center point of the kth radiation surface.

5. The method according to claim 1, characterized in that The step of determining whether the s th radiation surface and the k th radiation surface among the plurality of radiation surfaces face each other to obtain a first determination result includes: if or 0, then the first judgment result is obtained that the s-th radiation surface and the k-th radiation surface are not facing each other; if and 0, then the first judgment result that the s-th radiation surface and the k-th radiation surface are facing each other is obtained; in, represents the normal vector of the s-th radiating surface, represents the normal vector of the k-th radiating surface, is the center point of the s-th radiation surface, is the center point of the kth radiation surface.

6. The method according to claim 1, characterized in that The method further comprises: If the first judgment result indicates that the s th radiation surface and the k th radiation surface are not facing each other, it is determined that the diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero.

7. The method according to any one of claims 1 to 6, characterized in that: The method further comprises: If the fifth judgment result indicates that the first intersection point is not located inside the hth radiation surface, the diffusion angle coefficient of the sth radiation surface to the kth radiation surface is calculated according to the following formula: in, represents the angle between the normal vector of the sth radiation surface and the line connecting the first center, represents the angle between the normal vector of the kth radiation surface and the line connecting the first center, is the diffusion angle coefficient of the s-th radiation surface to the k-th radiation surface, is the length of the line segment from the center point of the s th radiation surface to the center point of the k th radiation surface, represents the area of ​​the kth radiating surface.

8. A device for calculating the diffusion angle coefficient in the field of high temperature engineering, used for calculating the radiation heat exchange between objects, characterized in that: The device comprises: A determination module, used for determining a plurality of radiation surfaces of a target space, wherein the shape of the radiation surfaces is a polygonal mesh; A judgment module is used to judge whether an s-th radiation surface and a k-th radiation surface among a plurality of radiation surfaces are facing each other, and obtain a first judgment result; if the first judgment result indicates that the s-th radiation surface and the k-th radiation surface are facing each other, determine an h-th radiation surface from a plurality of radiation surfaces, and the h-th radiation surface, the s-th radiation surface and the k-th radiation surface are different radiation surfaces; judge whether the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, and obtain a second judgment result; if the second judgment result indicates that the h-th radiation surface is within the radiation area of ​​the s-th radiation surface or the k-th radiation surface, judge whether the plane where the h-th radiation surface is located is consistent with the plane where the h-th radiation surface is located. A third judgment result is obtained by determining whether the straight line where the first center line is located is parallel to the straight line where the first center line is located, and the first center line is the line connecting the center point of the s th radiation surface and the center point of the k th radiation surface; if the third judgment result indicates that the plane where the h th radiation surface is located is not parallel to the straight line where the first center line is located, it is determined whether the first intersection point of the plane where the h th radiation surface is located and the straight line where the first center line is located is located between the s th radiation surface and the k th radiation surface, and a fourth judgment result is obtained; if the fourth judgment result indicates that the first intersection point is located between the s th radiation surface and the k th radiation surface, it is determined whether the first intersection point is located inside the h th radiation surface, and a fifth judgment result is obtained; A calculation module is configured to determine that a diffusion angle coefficient between the s th radiation surface and the k th radiation surface is zero if the fifth judgment result indicates that the first intersection point is located inside the h th radiation surface.

9. A computing device, characterized in that including memory and processor; One or more computer programs are stored in the memory, and the one or more computer programs include instructions; when the instructions are executed by the processor, the computing device executes the method as claimed in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium is used to store a computer program, and the computer program is used to execute the method according to any one of claims 1 to 7.

Images