CN113673006B - Physical thermal coupling visual complex geometric modeling method and system for whole reactor core - Google Patents

Abstract

The invention discloses a physical and thermal coupling visual complex geometric modeling method and system for a full reactor core, which relate to the technical field of nuclear reactor core design and have the technical scheme that: drawing to obtain cell geometry; automatically identifying the cell geometry to obtain a polygonal grid; performing physical and thermal attribute definition processing based on grids divided by a physical program to obtain corresponding grid cells; integrating a plurality of grid cells to complete the construction of the component and the reflecting layer structure, so as to obtain a component model; the sub-channel information of the thermal program is established based on the component model established by the physical program, and a component cell model is obtained; after integrating the plurality of component cell models, obtaining a full reactor core model; the mapping relation between grids of the physical program and the thermal program is output in a specific format. The invention can rapidly and simply carry out grid mapping and data transmission between the physical program and the thermal program, and reduces the coupling difficulty of the physical calculation program and the thermal calculation program.

Description

Physical thermal coupling visual complex geometric modeling method and system for whole reactor core

Technical Field

The invention relates to the technical field of nuclear reactor core design, in particular to a physical and thermal coupling visualization full-core complex geometric modeling method and system.

Background

The power calculated by the physical program can influence the fuel temperature, the coolant density and the coolant temperature calculated by the thermal program, and the thermal hydraulic parameters can influence the fission section so as to influence the power calculated by the physical program.

However, the physical program and the thermal calculation program have many differences in the size, shape and arrangement of the grids due to the characteristics of the respective physical and mathematical models, and if the grid shapes, the grid arrangement forms, the grid densities and the like established by the two types of programs are different, the grid mapping and the data transmission between the two types of programs are difficult in the coupling simulation of the two types of programs. Millions of physical program calculated full-core grids are used, and if manual modes are adopted for setting each geometric description, grid division and material property, the input modes are abnormally complex.

Therefore, how to research and design a rapid physical thermal coupling visual complex geometric modeling method and system for the whole reactor core is an urgent problem to be solved at present.

Disclosure of Invention

The invention aims to provide a physical thermodynamic coupling visual modeling method and system for complex geometry of a full reactor core, which are used for solving the problems of adaptability of complex geometry input of the transport calculation geometry of the full reactor core and convenience for users to use and reducing the coupling difficulty of a physical calculation program and a thermodynamic calculation program.

The technical aim of the invention is realized by the following technical scheme:

in a first aspect, a physical thermodynamic coupling visual complex geometry modeling method for a full core is provided, including the steps of:

adopting the basic function of AutoCAD to interactively draw the target shape to obtain the cell geometry;

automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid;

polygonal grids which are consistent with the physical program grid division and the thermal program grid division are subjected to physical and thermal attribute definition processing based on the grids of the physical program division, so that corresponding grid element grids are obtained; physical properties include material type, material temperature, material density; the thermal property comprises a channel type, and if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient property of the flux need to be defined;

adopting a modular modeling method to integrate a plurality of cell grids to complete the construction of the component and the reflecting layer structure, and obtaining a component model;

the sub-channel information of the thermal program is established based on the component model established by the physical program, and a component cell model is obtained; resetting the sub-channels by any combination of creating sub-channels, merging sub-channels, deleting sub-channels, adding grids, removing grids, generating according to grids and generating according to boundaries to finish modification;

adopting a modular modeling method to integrate and process a plurality of component cell models to obtain a full reactor core model;

the mapping relation between grids of the physical program and the thermal program is output in a specific format.

Furthermore, the modeling method is also used for preprocessing the intersecting line graph of the cell geometry based on the quadtree, and the specific process of preprocessing the intersecting line graph is as follows:

inserting all elements in the graph into corresponding quadtree nodes in a mode that the rectangular area of the maximum bounding box of the geometric graph of the grid element is continuously divided into four parts;

performing two-by-two intersection calculation on all geometric elements in the cell geometric figure and performing figure breaking treatment to obtain grid line segments in a discrete mode;

and carrying out pairwise intersection calculation on the subtree where the quadtree nodes are located during the intersection calculation.

Further, the process of automatically identifying the cell geometry by the grid automatic identification algorithm specifically comprises the following steps:

grid recognition is carried out on grid line segments obtained after graph preprocessing is carried out on cell geometry, and a point-edge topological relation is established according to recognition results;

identifying a ring in the point-edge topological relation, and establishing a grid topological structure according to the ring;

the grid topological structure is expressed in a BREP form in a three-dimensional entity, grids are formed on the basis of rings, and a plurality of adjacent grids form a region or a coarse grid required by physical analysis to obtain a polygonal grid.

Further, the construction process of the point-edge topological relation specifically comprises the following steps:

constructing non-coincident vertexes according to the coincidence relation of the geometric element endpoints in the graph;

the points are indexed, and the edges form the index of the end points of the edges based on the vertex index numbers, and the number of the edges connected on the vertices is recorded as the degree of the vertices.

Further, the identifying process of the ring in the point-edge topological relation specifically comprises the following steps:

optionally selecting one edge as an initial edge, and searching the next edge with the smallest included angle with the initial edge in the clockwise direction by taking the starting point to the ending point of the initial edge as a searching direction;

and taking the next searched edge as a search edge, and continuing searching according to the direction of the searched edge until the initial edge is searched.

Further, the process of identifying the ring in the point-edge topological relation further comprises the following steps:

and searching the primary search mark +1 in the forward direction by the edge, and searching the primary search mark +2 in the reverse direction; when the search mark is 3, the corresponding side does not participate in the search any more;

when the search marks of all the edges are 3, stopping searching, and completing the automatic construction process of the ring.

Further, the loops identified in the point-edge topological relation are divided into an inner loop and an outer loop, and the specific process for judging the inner loop and the outer loop is as follows:

judging whether the two rings are intersected according to the bounding box;

if so, obtaining an intersection by Boolean operation;

comparing whether the inner ring area is equal to the intersection area; if the two parts are equal, the inner ring is arranged inside the outer ring;

traversing the outer rings of all the inner rings;

finding out the outer ring with the smallest area in the outer rings as the smallest father outer ring of the corresponding inner ring;

all inner rings form a grid with the corresponding smallest parent outer ring.

Further, the definition processing of the physical and thermal attributes adopts 2-layer nested attribute description, and the specific description process is as follows:

defining a grid layer of a region in the geometric construction process, wherein the region consists of a plurality of grids with the same attribute, and the region is used for the description of a graphical interface and the maintenance of attribute information;

describing attribute information by using grids as units in an input file for outputting and forming a neutron transport program;

for a physical program, grids of the same material and the same temperature are defined as the same area;

for the thermal process, the fuel rods, fluid channels, insulating material, and thermally conductive material are defined as distinct regions.

Further, the whole core information includes whole core overall information and detailed geometric information of the components;

the global core information includes distribution of components, global core size information, and component geometry reference categories;

the detailed geometric information of the component includes edges, grids, attributes, and coarse grids.

In a second aspect, a physical thermodynamic coupling visualization full core complex geometry modeling system is provided, comprising:

the geometric drawing module is used for interactively drawing the target shape by adopting the basic function of AutoCAD to obtain the cell geometry;

the grid identification module is used for automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid;

the attribute definition module is used for carrying out physical and thermal attribute definition processing on polygon grids which are consistent with the physical program grid division and the thermal program grid division on the basis of the grids of the physical program division to obtain corresponding grid element grids; physical properties include material type, material temperature, material density; the thermal property comprises a channel type, and if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient property of the flux need to be defined;

the first integration module is used for integrating a plurality of grid cells to complete the construction of the component and the reflecting layer structure by adopting a modular modeling method to obtain a component model;

the channel construction module is used for establishing sub-channel information of the thermal program based on the component model established by the physical program to obtain a component cell model; resetting the sub-channels by any combination of creating sub-channels, merging sub-channels, deleting sub-channels, adding grids, removing grids, generating according to grids and generating according to boundaries to finish modification;

the second integration module is used for integrating the plurality of component cell models by adopting a modular modeling method to obtain a full reactor core model;

and the information output module is used for outputting the whole core information and the grid mapping relation between the physical program and the thermal program in a specific format.

Compared with the prior art, the invention has the following beneficial effects:

1. the full-core grid is millions in scale, the modeling of the component and the modeling of the full-core are performed by adopting a modular modeling method, the modeling can be completed quickly, the modeling efficiency is high, and meanwhile, the visual modeling is performed, so that the user can check errors in time conveniently;

2. the grids of the full-reactor core model established by the physical program and the thermal program are completely consistent, so that grid mapping and data transmission between the physical program and the thermal program are easier, and the coupling difficulty of the physical calculation program and the thermal calculation program is reduced.

Drawings

The accompanying drawings, which are included to provide a further understanding of embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiments of the invention. In the drawings:

FIG. 1 is a flow chart in an embodiment of the invention;

FIG. 2 is a schematic diagram of preprocessing a quadtree intersection graph in an embodiment of the invention, a is a space structure diagram, and b is a node storage structure diagram;

FIG. 3 is a schematic diagram illustrating the judgment of the inner ring and the outer ring according to an embodiment of the present invention;

FIG. 4 is a schematic diagram showing classification of rings in an embodiment of the present invention, a is a nested ring, b is a parallel ring, and c is a combination ring;

FIG. 5 is a schematic view of region division in an embodiment of the present invention;

FIG. 6 is a diagram of a physical and thermal program grid mapping relationship in an embodiment of the invention;

fig. 7 is a system block diagram in an embodiment of the invention.

Detailed Description

For the purpose of making apparent the objects, technical solutions and advantages of the present invention, the present invention will be further described in detail with reference to the following examples and the accompanying drawings, wherein the exemplary embodiments of the present invention and the descriptions thereof are for illustrating the present invention only and are not to be construed as limiting the present invention.

Example 1: a physical thermal coupling visual complex geometric modeling method for a full reactor core is shown in fig. 1, and is specifically realized by the following steps.

And step one, carrying out interactive drawing on the target shape by adopting the basic function of AutoCAD to obtain the cell geometry.

Since AutoCAD has a strong geometric description capability, various complex geometric figures can be drawn. Therefore, the geometric shapes of any shapes can be interactively drawn through functions such as straight lines, circles, circular arcs and the like provided by AutoCAD.

And step two, automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid.

AutoCAD can draw arbitrarily complex geometric figures, but does not have the function of identifying grids. If the relation between the grid and the edge is identified manually, the process is very complicated, and the editing is easy to make mistakes when the grid scale is large. To improve the efficiency of grid preprocessing, quadtrees are employed herein to improve the efficiency of graphics preprocessing; and then forming a complete grid topological structure by the discrete grid edges according to the positions of the discrete grid edges at the intersecting points.

The method comprises the following specific processes of preprocessing the intersecting line graph on the basis of the quadtree to carry out the preprocessing of the intersecting line graph on the cell geometry: inserting all elements in the graph into corresponding quadtree nodes in a mode that the rectangular area of the maximum bounding box of the geometric graph of the grid element is continuously divided into four parts; performing two-by-two intersection calculation on all geometric elements in the cell geometric figure and performing figure breaking treatment to obtain grid line segments in a discrete mode; and carrying out pairwise intersection calculation on the subtree where the quadtree nodes are located during the intersection calculation.

As shown in FIG. 2, the four-way tree structure is established by recursively dividing the geometric figure in the plane into tree structures of different levels according to the positions of the spatial quadrants, as shown in FIG. 2 a, each quadrant position is divided into four equal subspaces, and the recursion is continued until the level of the tree meets the condition that the current node object is smaller than a set value, and the segmentation is stopped. The structure of the quadtree is simple, and has high intersection calculation efficiency when the spatial data objects are distributed uniformly. After the graph is represented by the quadtree structure, the computational complexity can be reduced from O (n×2) to O (n× lgn).

The automatic identification process of the grid geometry by the grid automatic identification algorithm specifically comprises the following steps: grid recognition is carried out on grid line segments obtained after graph preprocessing is carried out on cell geometry, and a point-edge topological relation is established according to recognition results; identifying a ring in the point-edge topological relation, and establishing a grid topological structure according to the ring; the grid topological structure is expressed in a BREP form in a three-dimensional entity, grids are formed on the basis of rings, and a plurality of adjacent grids form a region or a coarse grid required by physical analysis to obtain a polygonal grid.

The construction process of the point-edge topological relation specifically comprises the following steps: constructing non-coincident vertexes according to the coincidence relation of the geometric element endpoints in the graph; the points are indexed, and the edges form the index of the end points of the edges based on the vertex index numbers, and the number of the edges connected on the vertices is recorded as the degree of the vertices.

As shown in fig. 3, the process of identifying the ring in the point-edge topological relation specifically includes: optionally selecting one edge as an initial edge, and searching the next edge with the smallest included angle with the initial edge in the clockwise direction by taking the starting point to the ending point of the initial edge as a searching direction; and taking the next searched edge as a search edge, and continuing searching according to the direction of the searched edge until the initial edge is searched. And, searching for the first search mark +1 in the forward direction and searching for the first search mark +2 in the reverse direction; when the search mark is 3, the corresponding side does not participate in the search any more; when the search marks of all the edges are 3, stopping searching, and completing the automatic construction process of the ring.

The ring identified in the point-edge topological relation is divided into an inner ring and an outer ring, and after the ring identification is completed, the ring is automatically identified into a grid according to the relation between the outer ring and the inner ring and between the two rings. According to the positional relationship table between the rings, a nested ring, a parallel ring, and a combination ring can be formed.

Nested loops as shown in fig. 4 a, ring 3 comprises ring 2, and ring 2 comprises ring 1, in a nested relationship. Meanwhile, we call the parent ring of the ring 1 to have a ring 2 and a ring 3. The parent ring of the ring 2 has a ring 3, the child ring of the ring 2 has a ring 1, the ring 3 has no parent ring, but there are child rings 1, 2.

Parallel rings as shown in fig. 4 b, ring 4 comprises ring 1, ring 2, and ring 3, with ring 1, ring 2, and ring 3 being in a positional parallel relationship. The parent rings of ring 1, ring 2 and ring 3 are ring 4.

As shown in fig. 4 c, the combined ring is shown in a side-by-side relationship with ring 1 and ring 2, ring 3 is the common parent ring of ring 1 and ring 2, and ring 3 is shown in a side-by-side relationship with ring 4, and the common parent ring is ring 5.

The specific process for judging the inner ring and the outer ring is as follows: judging whether the two rings are intersected according to the bounding box; if so, obtaining an intersection by Boolean operation; comparing whether the inner ring area is equal to the intersection area; if the two parts are equal, the inner ring is arranged inside the outer ring; traversing the outer rings of all the inner rings; finding out the outer ring with the smallest area in the outer rings as the smallest father outer ring of the corresponding inner ring; all inner rings form a grid with the corresponding smallest parent outer ring.

And thirdly, performing physical and thermal attribute definition processing on polygon grids which are consistent with the physical program grid division and the thermal program grid division on the basis of the grids of the physical program division to obtain corresponding grid element grids.

Since the physical program and the thermal program cell structures are consistent, the grid divided by the physical program is finer than the thermal program, and thus the physical program and the thermal program can use uniform grid division and are based on the grid divided by the physical program. The physical program needs to define the material, temperature, etc. information of each grid, and the thermal program needs to define the properties of each grid as fuel rods, fluid channels, heat insulating materials, heat conducting materials. If the geometric figure and the grids are completely constructed, then each grid is filled with attributes, and when the number of the grids is large, the user is very difficult to use; when a user modifies a grid property, if there are a large number of grid properties that need to be modified, the user modification is very labor intensive. For this purpose, a two-layer nested attribute description mode is adopted, and a special grid layer of the region is defined in the geometric construction process. The area is composed of a plurality of grids with the same attribute, the grids with the same material and the same temperature are defined as the same area for a physical program, and the fuel rod, the fluid channel, the heat insulating material and the heat conducting material can be defined as different areas for a thermal program, so that the workload of a user for maintaining the attribute information of the grids is greatly reduced due to the use of the areas.

The thermal program needs to define sub-channel information, which generally includes a grid crossing cells, but if sub-channels are defined on the basis of a component grid model, since the component grid can reach thousands or tens of thousands of scales, the sub-channel information can be set when the cell model is built, the sub-channels required by the grid building program calculation belonging to the fluid channel attribute can be set, and the attribute of the sub-channels can be defined.

The physical attribute comprises information such as material type, material temperature, material density and the like, the thermal attribute comprises channel type, if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient attribute of flux are required to be defined, the physical and thermal are respectively described by adopting nested attributes, the efficiency of setting the grid attribute is improved, the sub-channel information is required to be set for the thermal, the sub-channel is set for the grid with the attribute of the fluid channel, and the attribute inherits the grid attribute.

As shown in fig. 5, the problem includes 80 grids, 3 kinds of attribute information, if the user grids are constructed, the grid attribute information is directly filled in, and the user grids are filled in 80 times, namely the number of grids which need to be filled in is the same; if the attribute setting function of the "region" provided herein is applied, grids with the same attribute are combined into the same region (as shown, the grids are combined into 3 regions), attribute filling is performed for the region, and the user only needs to fill 3 times. The nesting property description method can effectively improve the efficiency of filling in the property by the user.

And step four, integrating a plurality of grid cells by adopting a modular modeling method to complete the construction of the component and the reflecting layer structure, thereby obtaining a component model.

The grid scale of the component can reach thousands or tens of thousands, and the component consists of repeated cell structures, if the geometry, the grid and the attribute setting are constructed by adopting a mode of directly interactively drawing all edges in AutoCAD software, the workload on a user is very large, errors are easy to occur, and a modular modeling mode is adopted. Firstly, a single regular cell geometry is developed, a single cell model is built, only one time for repeated cell structures is needed to be built, and then a filling technology or a splicing technology is applied to be combined into a final component model. In the modular modeling process, the complete automatic inheritance function of grids, physical and thermal attributes and coarse grids is provided, and all the physical and thermal information after splicing is automatically inherited. Meanwhile, in order to facilitate the convenience of the user, after the structural modeling of the components, the reflecting layer and the like is completed, a modification editing function is provided, and after modification, a mesh regeneration function is applied to generate a new model.

And fifthly, establishing sub-channel information of the thermal program based on the component model established by the physical program to obtain a component cell model.

In the process of physical and thermal coupling design, the division of the grid by the physical program is finer than that of the thermal program, the thermal program needs to define sub-channel information at the component level, and the physical program is mainly defined at the cell level, so that the definition is performed on the basis of the component model established by the physical program during the thermal program design. The subchannel information required by a thermal program can be defined in two ways: (1) When no sub-channel definition is performed during cell modeling, a grid building program belonging to the attribute of the fluid channel can be used for calculating the required sub-channel at the step, and the attribute of the sub-channel is defined; (2) Because the component grid scale can reach thousands or tens of thousands, the subchannel definition is difficult for users, so the subchannel information can be defined in the cell modeling process, the operation only needs to combine the subchannels crossing cells, and the attribute of the subchannels can be directly inherited.

Because the setting of the sub-channel information of the cell is already carried out during the modeling of the cell, the sub-channel information is inherited by adopting a modular modeling method, and only the sub-channel crossing the cell is needed to be reset in this step, and thousands of grids of the component are not needed to be operated.

And step six, integrating the plurality of component cell models by adopting a modular modeling method to obtain a full reactor core model.

The method is characterized in that the grid scale of the whole reactor core is millions, if geometry, grids and attribute setting are constructed in AutoCAD software in a 'direct interactive drawing mode', the whole reactor core is unrealizable for users, as the whole reactor core is composed of a plurality of components with repeated structures, a modular modeling mode is adopted, namely, different components are modeled through the previous four steps, the components with repeated structures only need to be modeled once, then a filling technology or a splicing technology is applied to be combined into a final whole reactor core model, and as the information of complete physics and thermal engineering, namely, the information of grid attributes, coarse grids, sub-channels and the like is defined in the modeling process of the components, all the information after splicing can be automatically inherited, so that users only need to set physical and thermal engineering boundary information of the whole reactor core, and users can check the correctness of information such as functional checking materials, temperature and the like through the attributes after splicing.

And step seven, outputting the whole reactor core information and the mapping relation between grids of the physical program and the thermal program in a specific format.

If the grid scale of the whole reactor core is millions, if the grid scale of the whole reactor core directly outputs information such as edges, grids, attributes and the like of the whole reactor core, the output quantity is larger, the efficiency is lower, the grid scale is suitable for a high-precision transportation calculation program of the whole reactor core, the output information comprises two parts of components, whole reactor core overall information and detailed geometric information of the components, the whole reactor core overall information comprises distribution conditions of the components, whole reactor core size information and component geometric reference types, the detailed geometric information of the components comprises information such as edges, grids, attributes, coarse grids and the like, the information quantity is larger, but only one time of geometric information is required to be output for the components with repeated structures, and the output information quantity is greatly reduced.

Because the full core model calculated by the physical program and the full core model calculated by the thermal program are established on the basis of a unified grid model, the grid mapping relation between structures such as fuel rods, fluid channels and sub-channels in the thermal program and the physical program can be output, and the grid mapping and data transmission in the coupling calculation of the physical program and the thermal program are facilitated.

As shown in fig. 6, the physical procedure is generally centered on the fuel rod core, a square space with a width of the rod core pitch is a computational cell, and the grid division is finer by dividing concentric circles of different radii from the midpoint outwards with the center line. The thermodynamic subchannel program uses a coolant flow channel surrounded by four fuel rods as a calculation subchannel, and the grid division of the calculation subchannel is thicker. Because the thermal program sub-channel information is established on the basis of the physical program grid division, the program can automatically output the grid mapping relation between the physical program grid and the thermal program, and meanwhile, because the thermal program sub-channel is used for corresponding to the grids of a plurality of physical programs, the grid crossing the sub-channel does not exist, the temperature information calculated by the thermal program can be directly applied to the corresponding physical program grid, and the data transmission between the physical program and the thermal program is greatly simplified.

Example 2: a physical thermal coupling visual complex geometric modeling system for a full reactor core is shown in fig. 7, and comprises a geometric drawing module, a grid identification module, an attribute definition module, a first integration module, a channel construction module and a second integration module and good information output module.

The geometric drawing module is used for interactively drawing the target shape by adopting the basic function of AutoCAD to obtain the cell geometry; the grid identification module is used for automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid; the attribute definition module is used for carrying out physical and thermal attribute definition processing on polygon grids which are consistent with the physical program grid division and the thermal program grid division on the basis of the grids of the physical program division to obtain corresponding grid element grids; physical properties include material type, material temperature, material density; the thermal property comprises a channel type, and if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient property of the flux need to be defined; the first integration module is used for integrating a plurality of grid cells to complete the construction of the component and the reflecting layer structure by adopting a modular modeling method to obtain a component model; the channel construction module is used for establishing sub-channel information of the thermal program based on the component model established by the physical program to obtain a component cell model; resetting the sub-channels by any combination of creating sub-channels, merging sub-channels, deleting sub-channels, adding grids, removing grids, generating according to grids and generating according to boundaries to finish modification; the second integration module is used for integrating the plurality of component cell models by adopting a modular modeling method to obtain a full reactor core model; and the information output module is used for outputting the whole core information and the grid mapping relation between the physical program and the thermal program in a specific format.

It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The foregoing detailed description of the invention has been presented for purposes of illustration and description, and it should be understood that the invention is not limited to the particular embodiments disclosed, but is intended to cover all modifications, equivalents, alternatives, and improvements within the spirit and principles of the invention.

Claims (10)

1. A physical thermal coupling visual complex geometric modeling method for a whole reactor core is characterized by comprising the following steps:

adopting the basic function of AutoCAD to interactively draw the target shape to obtain the cell geometry;

automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid;

the grids of the full reactor core model established by the physical program and the thermal program are completely consistent, polygonal grids which are consistent in the physical program grid division and the thermal program grid division are subjected to physical and thermal attribute definition processing based on the grids of the physical program division, and corresponding grid element grids are obtained; physical properties include material type, material temperature, material density; the thermal property comprises a channel type, and if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient property of the flux need to be defined;

adopting a modular modeling method to integrate a plurality of cell grids to complete the construction of the component and the reflecting layer structure, and obtaining a component model;

the sub-channel information of the thermal program is established based on the component model established by the physical program, and a component cell model is obtained; resetting the sub-channels by any combination of creating sub-channels, merging sub-channels, deleting sub-channels, adding grids, removing grids, generating according to grids and generating according to boundaries to finish modification;

adopting a modular modeling method to integrate and process a plurality of component cell models to obtain a full reactor core model;

the mapping relation between grids of the physical program and the thermal program is output in a specific format.

2. The physical thermal coupling visualization full-core complex geometric modeling method according to claim 1, wherein the modeling method is further based on a quadtree to perform intersection line graph preprocessing on cell geometry, and the specific process of the intersection line graph preprocessing is as follows:

inserting all elements in the graph into corresponding quadtree nodes in a mode that the rectangular area of the maximum bounding box of the geometric graph of the grid element is continuously divided into four parts;

performing two-by-two intersection calculation on all geometric elements in the cell geometric figure and performing figure breaking treatment to obtain grid line segments in a discrete mode;

and carrying out pairwise intersection calculation on the subtree where the quadtree nodes are located during the intersection calculation.

3. The physical thermal coupling visualization full core complex geometry modeling method according to claim 1 or 2, wherein the process of automatically identifying the cell geometry by the grid automatic identification algorithm is specifically as follows:

grid recognition is carried out on grid line segments obtained after graph preprocessing is carried out on cell geometry, and a point-edge topological relation is established according to recognition results;

identifying a ring in the point-edge topological relation, and establishing a grid topological structure according to the ring;

the grid topological structure is expressed in a BREP form in a three-dimensional entity, grids are formed on the basis of rings, and a plurality of adjacent grids form a region or a coarse grid required by physical analysis to obtain a polygonal grid.

4. The physical thermodynamic coupling visual modeling method for the complex geometry of the whole reactor core according to claim 3, wherein the construction process of the point-side topological relation is specifically as follows:

constructing non-coincident vertexes according to the coincidence relation of the geometric element endpoints in the graph;

the points are indexed, and the edges form the index of the end points of the edges based on the vertex index numbers, and the number of the edges connected on the vertices is recorded as the degree of the vertices.

5. The physical thermodynamic coupling visualization full core complex geometry modeling method of claim 3, wherein the process of identifying the loops in the point-to-edge topology is specifically as follows:

optionally selecting one edge as an initial edge, and searching the next edge with the smallest included angle with the initial edge in the clockwise direction by taking the starting point to the ending point of the initial edge as a searching direction;

and taking the next searched edge as a search edge, and continuing searching according to the direction of the searched edge until the initial edge is searched.

6. The method for modeling a complex geometry of a physical thermodynamic coupling visualization total core of claim 5, wherein the process of identifying the loops in the point-to-edge topology further comprises the steps of:

and searching the primary search mark +1 in the forward direction by the edge, and searching the primary search mark +2 in the reverse direction; when the search mark is 3, the corresponding side does not participate in the search any more;

when the search marks of all the edges are 3, stopping searching, and completing the automatic construction process of the ring.

7. The physical thermodynamic coupling visual complex geometry modeling method of the full reactor core according to claim 3, wherein the identified rings in the point-side topological relation are divided into an inner ring and an outer ring, and the specific process for judging the inner ring and the outer ring is as follows:

judging whether the two rings are intersected according to the bounding box;

if so, obtaining an intersection by Boolean operation;

comparing whether the inner ring area is equal to the intersection area; if the two parts are equal, the inner ring is arranged inside the outer ring;

traversing the outer rings of all the inner rings;

finding out the outer ring with the smallest area in the outer rings as the smallest father outer ring of the corresponding inner ring;

all inner rings form a grid with the corresponding smallest parent outer ring.

8. The physical and thermal coupling visualization full core complex geometry modeling method of claim 1, wherein the defining process of the physical and thermal attributes adopts 2-layer nested attribute description, and the specific description process is as follows:

defining a grid layer of a region in the geometric construction process, wherein the region consists of a plurality of grids with the same attribute, and the region is used for the description of a graphical interface and the maintenance of attribute information;

describing attribute information by using grids as units in an input file for outputting and forming a neutron transport program;

for a physical program, grids of the same material and the same temperature are defined as the same area;

for the thermal process, the fuel rods, fluid channels, insulating material, and thermally conductive material are defined as distinct regions.

9. The method of claim 1, wherein the total core information comprises total core population information and detailed geometric information of the components;

the global core information includes distribution of components, global core size information, and component geometry reference categories;

the detailed geometric information of the component includes edges, grids, attributes, and coarse grids.

10. A physical thermal coupling visual complex geometry modeling system for a full reactor core is characterized by comprising:

the geometric drawing module is used for interactively drawing the target shape by adopting the basic function of AutoCAD to obtain the cell geometry;

the grid identification module is used for automatically identifying the cell geometry through a grid automatic identification algorithm to obtain a polygonal grid;

the attribute definition module is used for completely conforming grids of the full-reactor core model established by the physical program and the thermal program, polygonal grids which are consistent in the physical program grid division and the thermal program grid division are used for carrying out physical and thermal attribute definition processing on the basis of the grids divided by the physical program to obtain corresponding grid element grids; physical properties include material type, material temperature, material density; the thermal property comprises a channel type, and if the channel is a fluid channel, the weight factor, the inlet flow ratio, the inlet temperature ratio and the resistance coefficient property of the flux need to be defined;

the first integration module is used for integrating a plurality of grid cells to complete the construction of the component and the reflecting layer structure by adopting a modular modeling method to obtain a component model;

the channel construction module is used for establishing sub-channel information of the thermal program based on the component model established by the physical program to obtain a component cell model; resetting the sub-channels by any combination of creating sub-channels, merging sub-channels, deleting sub-channels, adding grids, removing grids, generating according to grids and generating according to boundaries to finish modification;

the second integration module is used for integrating the plurality of component cell models by adopting a modular modeling method to obtain a full reactor core model;

and the information output module is used for outputting the whole core information and the grid mapping relation between the physical program and the thermal program in a specific format.