Multi-body Modeling and Simulation Mesh Stitching Method and Device Based on Physical Conservation

Abstract

The present application relates to a multi-body modeling and simulation grid splicing method and device based on physical conservation. The method includes: performing simulation modeling on the physical space region of the reactor thermal-hydraulic system to obtain two different sub-domains; the two different sub-domains form two interfaces; calculating the fluxes on the two interfaces; deriving based on the non-matching interface flux conservation law, performing coefficient allocation on the fluxes on the two interfaces, and calculating the interface field quantity values in the overlapping region of the two interfaces; performing weighted combination on the interface field quantity values in the overlapping region of the two interfaces, updating the weighted interface field quantity values, and reorganizing the updated interface field quantity values based on the flux conservation law to obtain an interpolation result; realizing simulation grid splicing according to the interpolation result. Using this method can improve the efficiency of reactor thermal-hydraulic analysis.

Description

Technical Field

[0001] This application relates to the technical field of reactor thermal-hydraulic system simulation, and particularly to a multi-body modeling and simulation grid splicing method and device based on physical conservation. Background Art

[0002] In the field of reactor thermal-hydraulic analysis, the application of CFD (Computational Fluid Dynamics) technology is extremely crucial. The internal structure of a reactor is intricate, especially the multi-body geometric structure part, which includes numerous independent regions such as fuel assemblies, control rod assemblies, coolant channels, etc. There are complex contacts and interfaces between these regions, such as between fuel rods and the wall of the coolant channel, between control rods and surrounding components, etc. When conducting CFD simulations, geometric modeling and mesh generation are fundamental steps. However, due to the complexity of the multi-body geometric structure inside the reactor, it is difficult to fully match the meshes of different component regions. For example, in the fuel assembly region, in order to accurately capture the flow and heat transfer characteristics of the coolant on the surface of the fuel rods, fine meshes need to be generated; while in some relatively simple spatial regions, the meshes can be appropriately sparse. In this case, the problem of mesh mismatch will occur.

[0003] In traditional simulation methods, non-conservative interpolation algorithms are often used to handle these mismatched mesh interfaces. However, in high-demand application scenarios such as reactor thermal-hydraulic analysis, the defects of non-conservative interpolation algorithms are fully exposed. On the one hand, the flow conditions of the coolant inside the reactor are complex. In some high-flow-rate regions (similar to high Reynolds number flows), non-conservative interpolation is difficult to ensure the strict conservation of physical quantities such as mass, momentum, and energy, resulting in large errors in key parts such as near the fuel rods and in key flow channels in the simulation results, seriously affecting the accurate evaluation of the reactor's thermal-hydraulic performance. On the other hand, during the operation of the reactor, the flow and heat transfer of the coolant change rapidly, with a small time scale and a large spatial gradient. The numerical stability of the non-conservative interpolation algorithm is poor, and the calculation results are prone to oscillation or divergence phenomena, making the simulation results unreliable. Moreover, non-conservative interpolation also requires additional post-processing steps to correct errors, which greatly reduces the efficiency for the work that needs to quickly obtain simulation results to guide the safe operation and optimized design of the reactor.

[0004] In recent years, although various grid splicing methods based on conservation principles have been proposed, they still face challenges when applied to reactor thermal-hydraulic analysis. For example, the method of modifying the grid topology structure will greatly increase the workload for such a complex multi-body geometric structure of the reactor, and may reduce the local grid quality, affecting the simulation accuracy; while the method of constructing a super grid has excessive computational and storage overheads due to the huge scale of reactor simulation, and it is difficult to be effectively implemented in actual reactor engineering analysis.

[0005] Therefore, in the field of reactor thermal-hydraulic analysis, how to simplify the algorithm complexity while maintaining flux conservation to achieve more efficient and accurate CFD simulation has become a key issue that needs to be urgently addressed in the field of reactor safe operation and optimized design. Summary of the Invention

[0006] Based on this, it is necessary to provide a multi-body modeling and simulation grid splicing method and device based on physical conservation that can improve the efficiency of reactor thermal-hydraulic analysis to address the above technical problems.

[0007] A multi-body modeling and simulation grid splicing method based on physical conservation, the method comprising:

[0008] The physical space of the reactor thermal hydraulic system is simulated and modeled to obtain two different subdomains; the two different subdomains form two interfaces; the fluxes on the two interfaces are calculated;

[0009] Based on the law of conservation of flux of non-matching interfaces, the fluxes on the two interfaces are coefficient-distributed, and the interface field magnitude in the overlapping region of the two interfaces is calculated; the interface field magnitudes in the overlapping region of the two interfaces are weightedly combined to obtain the weighted interface field magnitude;

[0010] The weighted interface field values are updated to obtain updated interface field values; the updated interface field values are reorganized based on the flux conservation law to obtain interpolation results; and simulation grid splicing is achieved based on the interpolation results.

[0011] A multi-body modeling and simulation grid splicing device based on physical conservation, the device comprising:

[0012] The flux calculation module is used to simulate and model the physical space area of the reactor thermal hydraulic system to obtain two different subdomains; the two different subdomains form two interfaces; and the flux on the two interfaces is calculated;

[0013] An interface field magnitude calculation module is used to derive the flux conservation law of the non-matching interface, distribute the coefficients of the flux on the two interfaces, and calculate the interface field magnitude in the overlapping area of the two interfaces; perform a weighted combination of the interface field magnitudes in the overlapping area of the two interfaces to obtain the weighted interface field magnitude;

[0014] The field value reorganization module is used to update the weighted interface field value to obtain the updated interface field value; reorganize the updated interface field value based on the flux conservation law to obtain the interpolation result; and realize the simulation grid splicing according to the interpolation result.

[0015] The above multi-body modeling and simulation grid splicing method based on physical conservation can, first of all, reasonably define and abstract complex physical phenomena and regions in the reactor thermal-hydraulic analysis by establishing sub-domains, providing a clear spatial range for subsequent numerical calculations, enabling the analysis to focus on specific regions, avoiding unnecessary calculations and interferences, and helping to improve the calculation efficiency. Among them, the sub-domains can form interfaces, which is conducive to the regional treatment of complex reactor systems. Different sub-domains may have different physical characteristics or flow states. Such a division can adopt more appropriate calculation methods or parameter settings according to the characteristics of different sub-domains. At the same time, the existence of the interfaces also facilitates the handling of the interactions between sub-domains, making the calculation more organized and improving the accuracy and efficiency of the analysis. By accurately calculating the flux, the transmission laws of mass, momentum, energy, etc. in the reactor thermal-hydraulic process can be better understood, providing basic data for subsequent analysis based on the conservation laws, helping to improve the accuracy and reliability of the analysis, and thus improving the overall analysis efficiency. Then, based on the flux conservation law, coefficients are derived and distributed to calculate the interface field quantity values. On the one hand, it ensures the conservation of physical quantities, which is crucial for reactor thermal-hydraulic analysis because conservation is the basic requirement for accurately describing physical processes and can avoid incorrect results caused by non-conservation of physical quantities, improving the accuracy of the analysis; on the other hand, reasonable coefficient distribution can more accurately calculate the interface flux according to the characteristics of different interfaces and the requirements of physical processes, providing accurate data for subsequent calculations and helping to improve the efficiency of the entire analysis process. By weighted combination of the interface field quantity values, the information of two interfaces can be comprehensively considered to generate comprehensive and consistent results. The weighting method can be reasonably set according to factors such as the importance of the interface and the distribution of physical quantities, so as to more accurately reflect the actual physical situation in the reactor, avoid the one-sidedness of single-interface information, improve the reliability and accuracy of the analysis results, and thus improve the analysis efficiency. Updating the weighted interface field quantity values can continuously adjust and optimize the field quantity values according to the real-time information and physical laws in the calculation process to make them closer to the real situation. Mapping based on the flux conservation law to obtain the interpolation result once again ensures the conservation of physical quantities throughout the calculation process, ensuring that the analysis results conform to physical reality. At the same time, the interpolation result can provide more accurate boundary conditions or initial conditions for subsequent calculations, enabling the calculation to converge more stably and quickly, improving the calculation efficiency, and providing reliable numerical solutions for reactor thermal-hydraulic analysis, helping to more deeply understand and analyze the thermal-hydraulic phenomena in the reactor, and thus improving the overall analysis efficiency. Description of the Drawings

[0016] Figure 1 It is a schematic flow chart of a multi-body modeling and simulation grid splicing method based on physical conservation in an embodiment;

[0017] Figure 2Schematic diagram of multi - sub - domain discretization for non - matching interfaces in an embodiment;

[0018] Figure 3 Schematic diagram of flux coefficient distribution and recombination in an embodiment;

[0019] Figure 4 Block diagram of a multi - body modeling and simulation grid splicing device based on physical conservation in an embodiment. Detailed implementation manners

[0020] In order to make the objectives, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application, and are not used to limit the present application.

[0021] In one embodiment, as Figure 1 shown, a multi - body modeling and simulation grid splicing method based on physical conservation is provided, including the following steps:

[0022] Step 102, perform simulation modeling on the physical space region of the reactor thermal - hydraulic system to obtain two different sub - domains; the two different sub - domains form two interfaces; calculate the fluxes on the two interfaces.

[0023] During the process of reactor thermal - hydraulic analysis, a sub - domain refers to the physical space region selected for numerical calculation and analysis in the reactor thermal - hydraulic system under study. It generally includes: the core region and the coolant system. Among them, the core region contains fuel assemblies, control rods, coolant channels, etc. Nuclear fission reactions occur in the fuel in the core, generating a large amount of heat. The coolant flows in the channels, taking the heat out of the core. The computational domain needs to accurately describe the fuel arrangement in the core, the coolant flow path, and related physical parameters to accurately analyze processes such as heat transfer, fluid flow, and neutron physics in the core. The coolant system includes reactor coolant pumps, steam generators, pressurizers, and connecting pipes, etc. The coolant circulates in these components, transferring the heat generated in the core to the steam generator to generate steam for power generation or other purposes. The computational domain should cover the entire process of the coolant system to study changes in parameters such as coolant flow rate, pressure, and temperature, as well as the stability and safety of the system.

[0024] By establishing subdomains, complex physical phenomena and regions in reactor thermal-hydraulic analysis can be reasonably defined and abstracted, providing a clear spatial range for subsequent numerical calculations. This enables the analysis to focus on specific regions, avoiding unnecessary calculations and interference, and helping to improve calculation efficiency. The formation of subdomain interfaces is conducive to the regional treatment of complex reactor systems. Different subdomains may have different physical characteristics or flow states. Such a division allows for the adoption of more appropriate calculation methods or parameter settings according to the characteristics of different subdomains. At the same time, the existence of interfaces also facilitates the handling of interactions between subdomains, making the calculation more organized and improving the accuracy and efficiency of the analysis. By accurately calculating the flux, the transport laws of mass, momentum, energy, etc. in the thermal-hydraulic process within the reactor can be better understood, providing basic data for subsequent analysis based on conservation laws, helping to improve the accuracy and reliability of the analysis, and thus enhancing the overall analysis efficiency.

[0025] Step 104: Based on the non-matching interface flux conservation law, conduct derivations, allocate coefficients to the fluxes on the two interfaces, and calculate the interface field quantity values in the overlapping region between the two interfaces; perform a weighted combination of the interface field quantity values in the overlapping region between the two interfaces to obtain the weighted interface field quantity values.

[0026] Derive and allocate coefficients according to the flux conservation law, calculate the interface field quantity values on the two interfaces, and combine them based on weighted averaging to generate a comprehensive and consistent result. Subsequently, the weighted result will be redistributed to the corresponding regions to ensure the precise conservation and balance of the flux or field quantity between the interfaces.

[0027] Step 106: Update the weighted interface field quantity values to obtain the updated interface field quantity values; map the updated interface field quantity values based on the flux conservation law to obtain the interpolation result; achieve simulation grid stitching according to the interpolation result.

[0028] Updating the weighted interface field quantity values can continuously adjust and optimize the field quantity values according to real-time information and physical laws during the calculation process, making them closer to the actual situation. Mapping based on the flux conservation law to obtain the interpolation result once again ensures the conservation of physical quantities throughout the calculation process, ensuring that the analysis results conform to physical reality. At the same time, the interpolation result can provide more accurate boundary conditions or initial conditions for subsequent calculations, enabling the calculation to converge more stably and quickly, improving the calculation efficiency, and providing a reliable numerical solution for reactor thermal-hydraulic analysis, helping to more deeply understand and analyze the thermal-hydraulic phenomena within the reactor, thereby enhancing the overall analysis efficiency. Through verification in actual engineering cases (such as reactor pressure vessel analysis), this application highlights its great value in improving simulation accuracy and solving non-matching interface problems in complex multi-body geometry situations, greatly improving the analysis efficiency.

[0029] The above multi-body modeling and simulation grid splicing method based on physical conservation can, first of all, reasonably define and abstract complex physical phenomena and regions in the reactor thermal-hydraulic analysis by establishing sub-domains, providing a clear spatial range for subsequent numerical calculations, enabling the analysis to focus on specific regions, avoiding unnecessary calculations and interferences, and helping to improve the calculation efficiency. Among them, the sub-domains can form interfaces, which is conducive to the regional treatment of complex reactor systems. Different sub-domains may have different physical characteristics or flow states. Such a division can adopt more appropriate calculation methods or parameter settings according to the characteristics of different sub-domains. At the same time, the existence of the interfaces also facilitates the handling of the interactions between sub-domains, making the calculation more orderly and improving the accuracy and efficiency of the analysis. By accurately calculating the flux, the transport laws of mass, momentum, energy, etc. in the reactor thermal-hydraulic process can be better understood, providing basic data for subsequent analysis based on the conservation laws, helping to improve the accuracy and reliability of the analysis, and further improving the overall analysis efficiency. Then, based on the flux conservation law, coefficients are derived and distributed to calculate the interface field quantity values. On the one hand, it ensures the conservation of physical quantities, which is crucial for reactor thermal-hydraulic analysis because conservation is the basic requirement for accurately describing physical processes and can avoid incorrect results caused by non-conservation of physical quantities, improving the accuracy of the analysis; on the other hand, reasonable coefficient distribution can more accurately calculate the interface flux according to the characteristics of different interfaces and the requirements of physical processes, providing accurate data for subsequent calculations and helping to improve the efficiency of the entire analysis process. By weighted combination of the interface field quantity values, the information of the two interfaces can be comprehensively considered to generate comprehensive and consistent results. The weighting method can be reasonably set according to factors such as the importance of the interface and the distribution of physical quantities, so as to more accurately reflect the actual physical situation in the reactor, avoid the one-sidedness of single-interface information, improve the reliability and accuracy of the analysis results, and further improve the analysis efficiency. Updating the weighted interface field quantity values can continuously adjust and optimize the field quantity values according to the real-time information and physical laws during the calculation process to make them closer to the real situation. Mapping based on the flux conservation law to obtain the interpolation result once again ensures the conservation of physical quantities throughout the calculation process, ensuring that the analysis results conform to physical reality. At the same time, the interpolation result can provide more accurate boundary conditions or initial conditions for subsequent calculations, enabling the calculation to converge more stably and quickly, improving the calculation efficiency, and providing reliable numerical solutions for reactor thermal-hydraulic analysis, helping to more deeply understand and analyze the thermal-hydraulic phenomena in the reactor, thus improving the overall analysis efficiency.

[0030] In one of the embodiments, calculating the fluxes on two interfaces includes:

[0031] The fluxes calculated on two interfaces are:

[0032] ;

[0033] Among them, represents the flux of the surface . represents the value of the interface field quantity on the interface . , represents the value of the interface field quantity of the interface . , represents the value of the interface field quantity on the interface . , represents the value of the interface field quantity of the interface . , represents the weight of each in . represents the weight of each in . represents the flux of the surface . respectively represent the values of the relevant variables, such as velocity and temperature, of the adjacent cells and of the interfaces and . represents the local number of each cell face on the interface . represents the local number of each cell face on the interface .

[0034] In a specific embodiment, as shown in Figure 2 , when the computational domain is divided into two subdomains and , these two subdomains are independently discretized, and a pair of non-matching interfaces will be formed between them, denoted as and , which are respectively composed of the surfaces of the subdomains and . These interfaces ensure the correct communication between the subdomains and are crucial for the accuracy of numerical simulation.

[0035] Assume that the interfaces , are respectively composed of the surfaces , . Based on the assumption or method of local connectivity: Each surface on an interface is coupled only with some or all of the surfaces on the opposite interface that overlap with it. In other words, only when there is a geometric relationship (such as overlap) between the surfaces on two interfaces, is numerical coupling allowed between them. The connectivity data between the interfaces is defined using two lists:

[0036] (1)

[0037] Among them, the operator 、 returns the area of the intersection surface.

[0038] In finite volume discretization, the flux is usually defined as the area expansion quantity related to the face F and can be calculated according to the values defined by the adjacent cells of the face. Specifically, the adjacent cells can be expressed as and , corresponding to the subdomains and respectively. The flux can be expressed as:

[0039] (2)

[0040] Among them, respectively represent the values of the relevant variables on the cells and , such as velocity, temperature, etc., represents the coefficient of finite volume discretization.

[0041] In equation (2), the terms on the right side can be respectively expressed as:

[0042] (3)

[0043] Therefore, equation (2) can be rewritten as:

[0044] (4)

[0045] Based on equation (4), for the faces and on the interfaces 、 , there is:

[0046] (5)

[0047] Combined with equation (3), it can be obtained that:

[0048] (6)

[0049] Due to the reason of grid splicing, based on the local connectivity assumption, there may be a situation where one grid cell corresponds to multiple grid cells on the opposite side at the splicing interface, then there will be:

[0050] (7)

[0051] Since the grid element mapping is not one-to-one, each in has a certain weight , and vice versa, each in has a certain weight , that is to say, equation (7) can be refined to:

[0052] (8)

[0053] The core idea of this application is to eliminate the geometric discontinuity between non-matching interfaces through interpolation technology while maintaining computational efficiency and accuracy. It does not require modifying the grid, but uses weighted interpolation to transfer numerical values from the relative interface to the current interface. Specifically, the calculation of the flux on the interface is similar to that of the internal face flux. By introducing an interpolation formula, the matrix of the finite volume problem can be directly assembled on the original interface. This method transforms the connectivity problem of non-matching interfaces into a unified numerical framework through interpolation technology, greatly reducing the computational complexity. This not only improves the computational efficiency but also simplifies the implementation process while maintaining consistency with traditional internal flux calculations. In summary, this application solves the geometric discontinuity problem of non-matching interfaces through interpolation technology, avoids complex grid adjustments, and achieves efficient and accurate numerical coupling.

[0054] In one of the embodiments, the derivation is based on the flux conservation law of non-matching interfaces, including:

[0055] Based on the flux conservation law of non-matching interfaces, the total fluxes on the two interfaces are equal, that is:

[0056] ;

[0057] where represents the flux of face , represents the flux of face ;

[0058] Further, by equating the coefficients, we get:

[0059] ;

[0060] where represents the value of the interface field quantity on interface , represents the value of the interface field quantity of interface , represents the value of the interface field quantity on interface , Display interface The interface field quantity value , Indicates each In The weight in Indicates each In The weight in

[0061] In a specific embodiment, the non-matching interface flux conservation law is: The total fluxes on the interfaces And Must be equal, that is:

[0062] (9)

[0063] Since the flux has a direction, Equation (9) can be discretized as:

[0064] (10)

[0065] Combined with Equation (8), there is:

[0066] (11)

[0067] To make Equation (11) hold, the corresponding terms on both sides of the equation And Can be made equal respectively, that is to say,

[0068] (12)

[0069] Furthermore, by making the coefficients equal, we can obtain:

[0070] (13)

[0071] In one of the embodiments, as Figure 3 Shown, calculate the interface field quantity values , On the overlapping regions of the two interfaces respectively And , and refine and allocate the two according to the area relationship of the overlapping region, we can obtain:

[0072] ;

[0073] Among them, Indicates the interface field quantity value of the surface , Indicates the interface field quantity value of the surface , Indicates each In the weights in represent each in the weights, the weights are expressed as:

[0074] ;

[0075] In one embodiment, the interfacial field quantity values on two interfaces are weighted and combined to obtain the weighted interfacial field quantity value, including:

[0076] For the refined part shared by the interfaces , assuming that , satisfy the local connectivity assumption, the field quantities of these two interfaces are weighted and averaged by a weighting coefficient which reflects the contribution of each interface to the shared part. After determining their respective contributions, the weighted combination is calculated as follows:

[0077] ;

[0078] where represents the weighting coefficient, represents the interfacial field quantity value of surface , represents the interfacial field quantity value of surface .

[0079] In one embodiment, the finally obtained weighted value represents the balanced contributions of the two interfaces in the shared overlapping region. Subsequently, this value is redistributed back to their respective corresponding regions, namely and , so as to ensure that the numerical solution remains consistent and the required physical quantities are conserved between the interfaces. At this time, the interfacial field quantity values and are updated as follows: [[ID=�9]]

[0080] ;

[0081] where represents the weighted interfacial field quantity value, represents the set that satisfies the local connectivity assumption with interface , represents the set that satisfies the local connectivity assumption with interface , represents the interfacial field quantity value of surface , represents the interfacial field quantity value of surface .

[0082] In one embodiment, mapping the updated interfacial field quantity values based on the flux conservation law includes:

[0083] In the finite volume discretization, There is usually a specific mapping relationship. Similarly, There is also the same mapping relationship. For example, for the diffusion term, the mapping relationship can be expressed as:

[0084] ;

[0085] where, represents the interfacial field quantity value of the interface , represents the interfacial field quantity value on the interface , represents the interfacial field quantity value of the interface , represents the interfacial field quantity value on the interface , represents the weight of each in , represents the weight of each in .

[0086] In one embodiment, mapping the updated interfacial field quantity values based on the flux conservation law to obtain an interpolation result includes:

[0087] Mapping the updated interfacial field quantity values based on the flux conservation law, the obtained interpolation result is:

[0088] .

[0089] In a specific embodiment, the interpolation result shows that the recombining of the weighting values ensures strict conservation of the flux, which is crucial for maintaining numerical stability and can accurately reflect the interaction between the two regions in the final solution.

[0090] In one embodiment, as Figure 4 shown, a multi-body modeling and simulation grid stitching device based on physical conservation, the device includes:

[0091] A flux calculation module 402, configured to perform simulation modeling on the physical space region of the reactor thermal-hydraulic system to obtain two different sub-domains; the two different sub-domains form two interfaces; calculate the fluxes on the two interfaces;

[0092] The interface field quantity value calculation module 404 is configured to perform derivation based on the non-matching interface flux conservation law, allocate coefficients to the fluxes on the two interfaces, and calculate the interface field quantity values in the overlapping region of the two interfaces; perform weighted combination on the interface field quantity values in the overlapping region of the two interfaces to obtain the weighted interface field quantity values;

[0093] The field quantity value recombination module 406 is configured to update the weighted interface field quantity values to obtain the updated interface field quantity values; recombine the updated interface field quantity values based on the flux conservation law to obtain the interpolation result; and implement simulation grid stitching according to the interpolation result.

[0094] It should be understood that although Figure 1 each step in the flowchart of Figure 1 is shown sequentially according to the indication of the arrow, these steps do not necessarily have to be executed sequentially in the order indicated by the arrow. Unless there is a clear indication in this article, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover,

[0095] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope recorded in this specification.

[0096] The above-described embodiments merely represent several implementation manners of the present application. The description is relatively specific and detailed, but it should not be construed as a limitation to the scope of the invention. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several deformations and improvements can still be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the appended claims.

Claims

1. A multi-body modeling and simulation grid splicing method based on physical conservation, characterized in that, The method includes: Performing simulation modeling on the physical space region of the reactor thermal-hydraulic system to obtain two different sub-domains; the two different sub-domains form two interfaces; calculating the fluxes on the two interfaces; Deriving based on the non-matching interface flux conservation law, performing coefficient assignment on the fluxes on the two interfaces, calculating the interface field quantity values in the overlapping region of the two interfaces; performing weighted combination on the interface field quantity values in the overlapping region of the two interfaces to obtain the weighted interface field quantity values; Updating the weighted interface field quantity values to obtain the updated interface field quantity values; mapping the updated interface field quantity values based on the flux conservation law to obtain the interpolation result; realizing simulation grid splicing according to the interpolation result.

2. The method according to claim 1, characterized in that, Calculating the fluxes on the two interfaces includes: Calculating the fluxes on the two interfaces as: Among them, represents the flux of surface . represents the value of the interface field quantity on interface , represents the value of the interface field quantity of interface , represents the value of the interface field quantity on interface , represents the value of the interface field quantity of interface , represents the weight of each in , represents the weight of each in , represents the flux of surface , respectively represent the values of the relevant variables on the adjacent cells and of interfaces , and , and the values of the relevant variables include velocity, temperature, i represents the local numbers of the cell faces on interface , j represents the local numbers of the cell faces on interface .​​​​ 3. The method according to claim 2, wherein Deriving based on the non-matching interface flux conservation law includes: Deriving based on the non-matching interface flux conservation law, the total fluxes on the two interfaces are equal, that is: Among them, represents the flux of surface , represents the flux of surface ; Further, by the equality of coefficients, we get: Among them, represents the interface field quantity value on the interface , represents the interface field quantity value of the interface , represents the interface field quantity value on the interface , represents the interface field quantity value of the interface , represents the weight of each in , represents the weight of each in .

4. The method according to claim 2, wherein Calculating the interface field quantity values in the overlapping region of the two interfaces includes: Calculating the interface field quantity values in the overlapping region of the two interfaces respectively as: Among them, represents the interfacial field quantity value of surface . represents the interfacial field quantity value of surface . represents the weight of each in . represents the weight of each in . The weight is expressed as: 。 5. The method according to claim 2, wherein Performing weighted combination on the interface field quantity values in the overlapping region of the two interfaces to obtain the weighted interface field quantity values includes: Performing weighted combination on the interface field quantity values in the overlapping region of the two interfaces to obtain the weighted interface field quantity values as: Among them, represents a weighting coefficient.

6. The method according to claim 2, characterized in that The way to update the weighted interface field quantity values is: Among them, represents the weighted interface field quantity value, represents the set that satisfies the local connectivity assumption with the interface ; represents the set that satisfies the local connectivity assumption with the interface ; represents the interface field quantity value of the surface ; represents the interface field quantity value of the surface ; It should be noted that there seems to be some inaccuracies or incompleteness in the original Chinese text, especially in the description of the relationship between different symbols and concepts. The above translation is based on the literal meaning as much as possible.

7. The method according to claim 2, wherein Mapping the updated interface field quantity values based on the flux conservation law includes: The relationship of the mapping is: Among them, represents the interface of the interface field quantity value , represents the interface of the interface field quantity value on it, represents the interface of the interface field quantity value , represents the interface of the interface field quantity value on it, represents each in the weight, represents each in the weight.

8. The method according to claim 7, wherein Mapping the updated interface field quantity values based on the flux conservation law to obtain the interpolation result includes: Mapping the updated interface field quantity values based on the flux conservation law to obtain the interpolation result as: 。 9. A multi-body modeling and simulation grid splicing device based on physical conservation, characterized in that, The device includes: A flux calculation module for performing simulation modeling on the physical space region of the reactor thermal-hydraulic system to obtain two different sub-domains; the two different sub-domains form two interfaces; calculating the fluxes on the two interfaces; An interface field quantity value calculation module for deriving based on the non-matching interface flux conservation law, performing coefficient assignment on the fluxes on the two interfaces, calculating the interface field quantity values in the overlapping region of the two interfaces; performing weighted combination on the interface field quantity values in the overlapping region of the two interfaces to obtain the weighted interface field quantity values; A field quantity value recombination module for updating the weighted interface field quantity values to obtain the updated interface field quantity values; mapping the updated interface field quantity values based on the flux conservation law to obtain the interpolation result; realizing simulation grid splicing according to the interpolation result.

Images