Method and device for generating two-dimensional hybrid grid for electromagnetic simulation based on region division

Abstract

The present invention discloses a method and device for generating two-dimensional hybrid grids for electromagnetic simulation based on region division, which includes step one, preprocessing of electromagnetic simulation region division; step two, generating different electromagnetic simulation sub-regions; step three, post-processing each of the divided electromagnetic simulation sub-regions; step four, generating a hybrid grid for electromagnetic simulation and outputting grid information. For the grid division method of narrow and long path-shaped metal conductors, the method of the present invention not only improves the grid generation speed and electromagnetic simulation calculation efficiency, but also maintains the calculation accuracy of electromagnetic simulation, which has obvious advantages over the existing methods.

Description

Technical Field

[0001] The present invention belongs to the technical field of electromagnetic simulation, and in particular relates to a method and a device for generating a two-dimensional hybrid grid for electromagnetic simulation based on region division. Background Art

[0002] Computational electromagnetics (CEM) is an important field in modern electromagnetic research, and electronic design automation (EDA) software is an effective tool for CEM, among which mesh generation technology is a key element.

[0003] According to the grid type classification, the computational grid generally used in modern CEM research is an unstructured grid. The nodes and units of the unstructured grid are formed randomly, which is more flexible when dealing with shapes and more suitable for complex shapes. According to the different unit geometric structures, the two-dimensional unstructured grid can be divided into triangular grids, quadrilateral grids and hybrid grids.

[0004] In computational electromagnetics simulation, triangular meshes can adapt to various geometric shapes and have strong flexibility in simulating complex electromagnetic problems in the real world; at the same time, triangular meshes can also provide relatively high numerical solution accuracy; but compared with rectangles or quadrilaterals, triangular meshes have lower electromagnetic computing efficiency. Quadrilateral mesh generation is relatively simple, with a regular structure, and the connection relationship between mesh cells is relatively simple, which makes matrix operations and computational operations more efficient in the process of solving electromagnetic simulation numerical solutions, but quadrilateral meshes have limited modeling capabilities for electromagnetic simulation areas containing complex geometric shapes or curved boundaries. Rectangular meshes are a type of quadrilateral mesh, but rectangular meshes have the advantage of controllable accuracy over quadrilateral meshes, that is, the mesh cells in the rectangular mesh can be refined or coarsened according to the requirements of the problem, but rectangular meshes, like quadrilateral meshes, cannot adapt to geometric simulation areas with complex shapes.

[0005] Hybrid grids are widely used in solving large-scale problems in computational electromagnetics. The purpose of the present invention is to improve the computational efficiency of electromagnetic simulation, and the improvement of the computational efficiency of electromagnetic simulation will be reflected in the significant reduction of operation time when solving large-scale problems. Rectangular and triangular hybrid grid generation is one of the relatively common methods in computational electromagnetic simulation; this hybrid grid generation technology combines the advantages of rectangular grids and triangular grids, which can better adapt to complex geometric shapes and improve the efficiency of numerical calculations. By selecting appropriate grid types in different areas, triangular grids can be used in areas that require higher precision, and rectangular grids can be used in larger and regular areas, so as to obtain better results in electromagnetic simulation. Hybrid grid generation methods are widely used in electromagnetic finite element method, finite difference method and other methods for numerically solving electromagnetic problems.

[0006] The present invention proposes a new hybrid grid generation method for path-shaped metals (such as wires, planar spiral inductors, etc.) in the field of electromagnetic simulation. The grid generation method can generate grids of path-shaped simulation areas more quickly, and the generated grids can further improve the efficiency of path-shaped simulation areas during electromagnetic calculations. Summary of the invention

[0007] The purpose of the present invention is to solve the problem that the mesh generated by the existing mesh generation method applied to the path-shaped electromagnetic simulation area cannot achieve a high electromagnetic calculation efficiency, and proposes a two-dimensional hybrid mesh generation method and device for electromagnetic simulation based on regional division. The core idea of ​​the present invention is to reasonably divide the target simulation area and adopt appropriate mesh generation strategies for different simulation sub-areas. For the path-shaped simulation area, this method can quickly generate a mesh suitable for more efficient electromagnetic calculations, overcoming the bottleneck of traditional methods in this case. The present invention divides the two-dimensional electromagnetic simulation area into a number of triangular sub-areas and rectangular sub-areas, and uses the Delaunay algorithm to generate triangular meshes for the triangular sub-areas, while the rectangular sub-areas directly draw rectangular meshes according to the boundary points.

[0008] In a first aspect, the present invention provides a method for generating a two-dimensional hybrid grid for electromagnetic simulation based on region division, comprising the following steps:

[0009] Step 1: Preprocessing of electromagnetic simulation area division

[0010] S1.1 obtain electromagnetic simulation data, and use the electromagnetic simulation data to establish an electromagnetic simulation plane region; then remove useless vertices on the boundary of the electromagnetic simulation plane region;

[0011] S1.2 calculates the normal vector of the boundary of the electromagnetic simulation plane area after processing in S1.1;

[0012] S1.3 Insert area division auxiliary points into all opposite sides on the boundary of electromagnetic simulation area;

[0013] Step 2: Generate different electromagnetic simulation sub-areas;

[0014] Step 3: Post-process each divided electromagnetic simulation sub-area;

[0015] Step 4: Generate hybrid mesh for electromagnetic simulation

[0016] S4.1 Use the set P and set Q constructed in step 2.1; insert mesh generation auxiliary points on the boundaries of the post-processed electromagnetic simulation plane sub-regions in the set P and set Q according to the mesh size required by the electromagnetic simulation design;

[0017] S4.2 For the electromagnetic simulation plane sub-region processed in step S4.1, if it comes from the set P, the electromagnetic simulation plane sub-region is triangularly meshed, all boundary vertices after the mesh generation auxiliary points are inserted in step 4.1 are used as boundary conditions, and the post-processed electromagnetic simulation plane sub-region is triangulated through Constrained Delaunay; if it comes from the set Q, the electromagnetic simulation plane sub-region is rectangularly meshed, and a rectangular mesh is generated for the post-processed electromagnetic simulation plane sub-region;

[0018] S4.3 merges all electromagnetic simulation plane sub-regions processed in step 4.2 to obtain an overall hybrid mesh, and then removes duplicate edges and points.

[0019] In a second aspect, the present invention provides a two-dimensional hybrid grid generation system for electromagnetic simulation, characterized in that it includes:

[0020] A preprocessing module, which obtains electromagnetic simulation data, establishes an electromagnetic simulation plane area using the electromagnetic simulation data, and preprocesses the electromagnetic simulation plane area;

[0021] A partitioning module generates different electromagnetic simulation sub-areas for the preprocessed electromagnetic simulation plane area;

[0022] A post-processing module is used to post-process each divided electromagnetic simulation sub-area;

[0023] The hybrid grid generation module generates a hybrid grid for electromagnetic simulation for each electromagnetic simulation sub-area after post-processing.

[0024] In a third aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method described.

[0025] In a fourth aspect, the present invention provides a computing device, comprising a memory and a processor, wherein the memory stores executable code, and when the processor executes the executable code, the described method is implemented.

[0026] The beneficial effects of the present invention are:

[0027] For the meshing method of narrow and long path-shaped metal conductors, the present invention is based on the idea of ​​divide and conquer, and with the help of the opposite edge relationship of each edge on the boundary of the electromagnetic simulation plane area, proposes the electromagnetic simulation plane area division method described in steps S1-S3; the electromagnetic simulation plane area division method described in steps S1-S3 not only improves the electromagnetic simulation calculation efficiency, but also maintains the calculation accuracy of the electromagnetic simulation, which has obvious advantages over the existing methods; applying the method described in steps S1-S3 can improve the quality of the meshing unit in the mesh generation process and speed up the entire electromagnetic simulation calculation process.

[0028] At the same time, the region division preprocessing method proposed in step S1 makes the subsequent segmentation process simpler and faster, making the electromagnetic simulation plane region division method described in S1-S3 more efficient than the existing method. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] Figure 1 This is a schematic diagram of inserting auxiliary points into opposite edges during the preprocessing process.

[0030] Figure 2 The initial electromagnetic simulation domain is input, a planar spiral inductor.

[0031] Figure 3 Schematic diagram of the process of inserting auxiliary points on the boundary of the simulation area during preprocessing.

[0032] Figure 4 It is the region boundary of the preprocessed planar spiral inductor.

[0033] Figure 5 This is a schematic diagram of the electromagnetic simulation sub-areas after the planar spiral inductor simulation area is divided.

[0034] Figure 6 The hybrid grid for planar spiral inductors is ultimately used in electromagnetic simulation calculations.

[0035] Figure 7 The initial electromagnetic simulation domain is input, and a virtual inductor is implemented using other circuit elements.

[0036] Figure 8 is the sub-region boundary of the preprocessed virtual inductor.

[0037] Fig. 9 The virtual inductor is divided into simulation sub-areas after the simulation area is divided.

[0038] Fig.10 Each simulation sub-area after post-processing of the virtual inductor.

[0039] Fig.11 It is a flow chart of the present invention. DETAILED DESCRIPTION

[0040] In order to more clearly illustrate the problem solved by the present invention and the technical solution adopted, the specific implementation mode of the present invention is described below with reference to the drawings.

[0041] A two-dimensional hybrid grid generation method for electromagnetic simulation based on region division is implemented in the following scenarios:

[0042] For metal layer simulation only, there are paths in the simulation area (such as Figure 2 The planar spiral inductor shown and Figure 7 In the case of the virtual inductor shown in the figure, there may be both regular simulation areas and irregular simulation areas in the local area; the use of pure triangular partitioning can adapt well to the irregular simulation area, but it also sacrifices the calculation efficiency of electromagnetic simulation; the use of rectangular partitioning has higher electromagnetic simulation calculation efficiency, but cannot adapt well to irregular simulation areas.

[0043] For the above implementation scenario, among the currently disclosed fast hybrid grid generation methods, no method can provide more efficient electromagnetic simulation calculations while maintaining high calculation accuracy. However, the method of the present invention can effectively solve this problem.

[0044] It should be noted that although the present invention is a method for generating a mixed grid of rectangles and triangles, when the method of the present invention is used for segmentation Figure 7 When the overall regular virtual inductor is shown in the geometry, the resulting mesh is a pure rectangular mesh.

[0045] When applying the two-dimensional hybrid grid generation method for electromagnetic simulation based on region partitioning, Fig.11 It mainly includes the following steps:

[0046] Step 1: Preprocessing of electromagnetic simulation area division

[0047] S1.1 obtain electromagnetic simulation data, and use the electromagnetic simulation data to establish an electromagnetic simulation plane region; then remove useless vertices on the boundary of the electromagnetic simulation plane region;

[0048] The removal of useless vertices specifically includes: determining and removing duplicate points and multiple collinear points;

[0049] S1.2 calculates the normal vector of the boundary of the electromagnetic simulation plane area after processing in S1.1; specifically, calculates the normal vector of each edge pointing to the inside of the electromagnetic simulation plane area according to the order of arrangement of vertices (clockwise or counterclockwise);

[0050] S1.3 Insert area division auxiliary points on all opposite sides of the electromagnetic simulation area boundary; specifically:

[0051] S1.3.1 Let C be the set of all edges on the boundary of the electromagnetic simulation plane region after S1.1. The two edges E satisfying formulas (1)-(3) i 、E j is the required opposite side;

[0052]

[0053]

[0054]

[0055] Where E i represents the i-th edge on the boundary of the electromagnetic simulation plane region after step 1.1, E j represents the jth edge on the boundary of the electromagnetic simulation plane region after step 1.1, E k represents the kth edge on the boundary of the electromagnetic simulation plane region after processing in step 1.1; It is side E i The normal vector of It is side E j The normal vector of is the zero vector; D(E i ,E j ) is a parallel side E i and E j The distance, D(E i ,E k ) is a parallel side E i and E k The distance, D(E k ,E j ) is a parallel side E k and E j distance;

[0056] S1.3.2 Insert area division auxiliary points on all opposite sides:

[0057] Use the opposite side E obtained in step 1.3.1 i and E j , for E i The two endpoints are related to E j Take two perpendicular lines to E jThe intersection of E j The two endpoints are related to E i Take two perpendicular lines to E i The intersection of is used as the auxiliary point insertion position; the region division auxiliary point is directly inserted at the auxiliary point insertion position, and finally the electromagnetic simulation plane region after the region division auxiliary point is inserted is obtained;

[0058] like Figure 1 As shown, find a pair of parallel sides E1 and E2, draw a perpendicular line through the two endpoints of E1, and the intersection of the perpendicular line and E2 is the auxiliary point to be inserted; E2 should also do the same operation on E1; the process of inserting auxiliary points for planar spiral inductors is shown in Figure 3 Finally, the electromagnetic simulation plane area is obtained after inserting the area division auxiliary points; the plane spiral inductor simulation area boundary is as follows Figure 4 shown.

[0059] Step 2: Generate different EM simulation sub-areas

[0060] S2.1 constructs a stack T, and initializes the top simulation area of ​​T to the electromagnetic simulation plane area preprocessed in step 1; constructs a set P to store the electromagnetic simulation sub-areas finally generated for triangular mesh partitioning; constructs a set Q to store the electromagnetic simulation sub-areas finally generated for rectangular mesh partitioning.

[0061] S2.2 Take out the top simulation area from T, record it as R, and find opposite sides of equal length on the boundary of R. A pair of opposite sides constitute a rectangular simulation area;

[0062] S2.3 determines whether each rectangular simulation area is legal, specifically:

[0063] Determine whether the rectangular simulation area exceeds the boundary of the electromagnetic simulation plane area preprocessed in step 1; if so, the current rectangular simulation area is considered illegal, and then the simulation area R is placed in the set P, and then step 2.6 is executed; otherwise, the current rectangular simulation area is considered legal, and then step 2.4 is executed;

[0064] like Figure 8 The virtual inductor simulation sub-area after preprocessing shown in the figure has two marked edges that just form a rectangular simulation area, but the rectangular simulation area formed by it has exceeded the boundary of the initial simulation area, that is, it is an illegal rectangular simulation area;

[0065] S2.4 determines whether the size of the current rectangular simulation area is large enough relative to the electromagnetic simulation plane area preprocessed in step 1; the criteria for determining whether the size of the rectangular simulation area is large enough are as follows:

[0066] (1) The ratio of the area of ​​the current rectangular simulation region to the area of ​​the electromagnetic simulation plane region after preprocessing in step 1;

[0067] (2) the ratio of the radius of the circumscribed circle of the current rectangular simulation area to the radius of the circumscribed circle of the electromagnetic simulation plane area preprocessed in step 1;

[0068] If any one of the two ratios is greater than the threshold value A (0.13 is used in this embodiment), it is considered that the size of the current rectangular simulation area is large enough, and then step S2.5 is executed; otherwise, it is considered that the size of the current rectangular simulation area is too small, and then the simulation area R is put into the set P, and then step S2.6 is executed;

[0069] S2.5 splits the current rectangular simulation region from the simulation region R, and then stores the rectangular simulation region as an electromagnetic simulation sub-region into the cache queue; if other regions remain after the simulation region R is split, the other remaining regions after the split are placed on the top of the stack of T, and then returns to step 2.2; otherwise, proceed to step 2.6;

[0070] S2.6 determines whether there are other simulation areas in stack T. If yes, it returns to step 2.2. Otherwise, it goes to step 3. The simulation sub-areas after the planar spiral inductor is divided are as follows: Figure 5 shown.

[0071] Step 3: Post-process each electromagnetic simulation sub-area

[0072] S3.1 remove multiple collinear vertices on the boundaries of all electromagnetic simulation sub-regions;

[0073] S3.2 builds a stack and initializes the stack to be empty; pushes all electromagnetic simulation sub-regions in the cache queue into the stack;

[0074] S3.3 Take out the top simulation area R from the stack, set the edge E to be any edge on the boundary of R; set the point N to be any point on the boundary of the electromagnetic simulation plane area after preprocessing in step 1; search according to the post-processing strategy;

[0075] The post-processing strategy is specifically:

[0076] (1) If point N is on edge E but not on the two endpoints of E, draw a perpendicular line to point N about edge E to divide the current rectangular simulation area into two, obtain two rectangular simulation areas, push the two rectangular simulation areas into the stack, and then execute step S3.5;

[0077] (2) If point N is not on edge E or at the two endpoints of E, put the simulation area R into the set Q constructed in step 2.1, and then go to step 3.5;

[0078] S3.5 determines whether the stack is empty; if it is empty, the post-processed electromagnetic simulation plane area is obtained and the process goes to step 4; if it is not empty, the process returns to S3.3;

[0079] At once Figure 7 For the virtual inductor shown, the simulation sub-area without post-processing is as follows Fig. 9 As shown in the figure, the simulation sub-area after post-processing is shown in Fig.10 shown.

[0080] Step 4: Generate hybrid mesh for electromagnetic simulation

[0081] S4.1 Use the set P and set Q constructed in step 2.1; insert mesh generation auxiliary points on the boundaries of the post-processed electromagnetic simulation plane sub-regions in the set P and set Q according to the mesh size required by the electromagnetic simulation design;

[0082] The insertion of auxiliary points is specifically: let any edge on the boundary of the post-processed electromagnetic simulation plane sub-area in the set P and the set Q be E, the side length of E be L(E), let the grid size designed for electromagnetic simulation requirements be Z, then the number X of grid generation auxiliary points that need to be inserted on the edge E is obtained by formula (4):

[0083]

[0084] Auxiliary points are generated by uniformly inserting the grid into the edge E according to the insertion number X;

[0085] S4.2 Take out the electromagnetic simulation plane sub-region from the set P and the set Q used in step 4.1; if the electromagnetic simulation plane sub-region is taken out from the set P, the electromagnetic simulation plane sub-region is triangularly meshed, and all boundary vertices after the mesh generation auxiliary points are inserted in step 4.1 are used as boundary conditions, and the post-processed electromagnetic simulation plane sub-region is triangulated through Constrained Delaunay; if the electromagnetic simulation plane sub-region is taken out from the set Q, the electromagnetic simulation plane sub-region is rectangularly meshed, and a rectangular mesh is generated for the post-processed electromagnetic simulation plane sub-region;

[0086] The specific process of generating a rectangular grid is:

[0087] Draw perpendicular lines about the boundary for auxiliary points of mesh generation on the boundary of the post-processed electromagnetic simulation plane sub-region to obtain a rectangular grid;

[0088] S4.3 merges all the electromagnetic simulation plane sub-regions processed in step 4.2 to obtain an overall hybrid mesh, and then removes duplicate edges and points. The final hybrid mesh of the planar spiral inductor is as follows: Figure 6 As shown;

[0089] S4.4 generates grid information for EDA calculation;

[0090] The grid information includes the positions of grid points and the corresponding relationships among unit points, lines, and surfaces.

Claims

1. A two-dimensional hybrid grid generation method for electromagnetic simulation based on region division, characterized in that: The following steps are involved: Step 1: Preprocessing of electromagnetic simulation area division S1.1 obtain electromagnetic simulation data, and use the electromagnetic simulation data to establish an electromagnetic simulation plane region; then remove useless vertices on the boundary of the electromagnetic simulation plane region; S1.2 calculates the normal vector of the boundary of the electromagnetic simulation plane area after processing in S1.1; S1.3 Insert area division auxiliary points on all opposite sides of the electromagnetic simulation area boundary; specifically: S1.3.1 Let C be the set of all edges on the boundary of the electromagnetic simulation plane region after S1.

1. The two edges E satisfying formulas (1)-(3) i 、E j is the required opposite side; Where E i represents the i-th edge on the boundary of the electromagnetic simulation plane region after step 1.1, E j represents the jth edge on the boundary of the electromagnetic simulation plane region after step 1.1, E k represents the kth edge on the boundary of the electromagnetic simulation plane region after processing in step 1.1; It is side E i The normal vector of It is side E j The normal vector of is the zero vector; D(E i ,E j ) is a parallel side E i and E j The distance, D(E i ,E k ) is a parallel side E i and E k The distance, D(E k ,E j ) is a parallel side E k and E j distance; S1.3.2 Insert area division auxiliary points on all opposite sides: Use the opposite side E obtained in step 1.3.1 i and E j , for E i The two endpoints are related to E j Take two perpendicular lines to E j The intersection of E j The two endpoints are related to E i Take two perpendicular lines to E i The intersection of is inserted as an auxiliary point; Directly insert the region division auxiliary point at the auxiliary point insertion position, and finally obtain the electromagnetic simulation plane region after the region division auxiliary point is inserted; Step 2: Generate different EM simulation sub-areas; specifically: S2.1 constructs a stack T, initializes the top simulation region of T to be the electromagnetic simulation plane region preprocessed in step 1; constructs a set P to store the electromagnetic simulation sub-regions finally generated for triangular mesh partitioning; constructs a set Q to store the electromagnetic simulation sub-regions finally generated for rectangular mesh partitioning; S2.2 Take out the top simulation area from T, record it as R, and find opposite sides of equal length on the boundary of R. A pair of opposite sides constitute a rectangular simulation area; S2.3 judge whether each rectangular simulation area is legal one by one; S2.4 determines whether the size of the current rectangular simulation area is large enough relative to the electromagnetic simulation plane area preprocessed in step 1; S2.5 splits the current rectangular simulation region from the simulation region R, and then stores the rectangular simulation region as an electromagnetic simulation sub-region into the cache queue; if other regions remain after the simulation region R is split, the other remaining regions after the split are placed on the top of the stack of T, and then returns to step 2.2; otherwise, proceed to step 2.6; S2.6 determines whether there are other simulation areas in stack T, if yes, returns to step 2.2, otherwise goes to step 3; Step 3: Post-process each divided electromagnetic simulation sub-area; Step 4: Generate hybrid mesh for electromagnetic simulation S4.1 inserting mesh generation auxiliary points on the boundary of the post-processed electromagnetic simulation plane sub-region according to the mesh size required by the electromagnetic simulation design; Specifically: S4.1 Use the set P and set Q constructed in step 2.1; insert mesh generation auxiliary points on the boundaries of the post-processed electromagnetic simulation plane sub-regions in the set P and set Q according to the mesh size required by the electromagnetic simulation design; S4.2 For the electromagnetic simulation plane sub-region processed in step S4.1, if it comes from the set P, the electromagnetic simulation plane sub-region is triangularly meshed, all boundary vertices after the mesh generation auxiliary points are inserted in step 4.1 are used as boundary conditions, and the post-processed electromagnetic simulation plane sub-region is triangulated through Constrained Delaunay; if it comes from the set Q, the electromagnetic simulation plane sub-region is rectangularly meshed, and a rectangular mesh is generated for the post-processed electromagnetic simulation plane sub-region; The insertion grid generation auxiliary points are specifically: Assume that any edge on the boundary of the post-processed electromagnetic simulation plane sub-region is E, the side length of E is L(E), and the grid size designed for electromagnetic simulation requirements is Z. Then, the number X of grid generation auxiliary points that need to be inserted on the edge E is obtained by formula (4): Auxiliary points are generated by uniformly inserting the grid into the edge E according to the insertion number X; S4.2, according to the meshing type marked on each electromagnetic simulation plane sub-region obtained in step 4.1, determine whether to perform triangular meshing or rectangular meshing; if triangular meshing is performed, all boundary vertices after the meshing auxiliary points are inserted in step 4.1 are used as boundary conditions, and the post-processed electromagnetic simulation plane sub-region is triangulated through Constrained Delaunay; otherwise, rectangular meshing is used to generate a rectangular mesh for the post-processed electromagnetic simulation plane sub-region; S4.3 merges all electromagnetic simulation plane sub-regions processed in step 4.2 to obtain an overall hybrid mesh, and then removes duplicate edges and points.

2. The method according to claim 1, characterized in that Step S1.2 specifically includes: calculating the normal vector of each edge pointing to the inside of the electromagnetic simulation plane area according to the order in which the vertices are arranged.

3. The method according to claim 1, characterized in that Step 3 specifically is: S3.1 remove multiple collinear vertices on the boundaries of all electromagnetic simulation sub-regions; S3.2 builds a stack and initializes the stack to be empty; pushes all electromagnetic simulation sub-regions in the cache queue into the stack; S3.3 Take out the top simulation area R from the stack, set the edge E to be any edge on the boundary of R; set the point N to be any point on the boundary of the electromagnetic simulation plane area after preprocessing in step 1; search according to the post-processing strategy; S3.5 determines whether the stack is empty; if it is empty, the post-processed electromagnetic simulation plane area is obtained and the process goes to step 4; If it is not empty, return to S3.

3.

4. A two-dimensional hybrid grid generation system for electromagnetic simulation implementing the method described in any one of claims 1 to 3, characterized in that: include: A preprocessing module, which acquires electromagnetic simulation data and uses the electromagnetic simulation data to establish an electromagnetic simulation plane area; Preprocessing of electromagnetic simulation plane area; A partitioning module generates different electromagnetic simulation sub-areas for the preprocessed electromagnetic simulation plane area; A post-processing module is used to post-process each divided electromagnetic simulation sub-area; The hybrid grid generation module generates a hybrid grid for electromagnetic simulation for each electromagnetic simulation sub-area after post-processing.

5. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method according to any one of claims 1 to 3.

6. A computing device, comprising a memory and a processor, wherein the memory stores executable code, and when the processor executes the executable code, the method according to any one of claims 1 to 3 is implemented.

Images