Method and device for simplifying large-scale simulation

Abstract

The invention relates to a method and device for simplifying large-scale simulation. The method comprises the steps that an analysis task for a target structure is obtained, and the target structure is subjected to structure division to obtain a plurality of substructures; determining a propagation space required for executing the analysis task according to the analysis task, and determining a plurality of first modes according to the size of the propagation space, wherein the plurality of first modes comprise a plurality of propagation modes and a plurality of evanescent modes; performing electromagnetic simulation of each substructure by using all first modes of the electromagnetic waves to obtain a substructure near-field scattering matrix of each substructure; and determining a target near-field scattering matrix of the target structure based on the near-field scattering matrixes of all the substructures. The near-field scattering matrix of the complete target structure is accurately predicted, the electromagnetic simulation process of the complex structure can be simplified in the mode that a large-scale simulation problem is divided into a plurality of small-scale simulation problems, the computing resource requirement of full-wave simulation is greatly reduced, and therefore the near-field characteristic analysis task of the complex structure is dealt with.

Description

Technical Field

[0001] The present disclosure relates to the technical field of electromagnetic simulation, and in particular to a method and device for simplifying large-scale simulation. Background Art

[0002] With the support of electromagnetic wave simulation technology, more and more novel electromagnetic properties can be realized through complex structural design. However, as the structure becomes more and more complex, the scale of numerical simulation has increased significantly, which leads to higher memory requirements and longer solution time, bringing huge challenges to electromagnetic simulation. This trend may affect the efficiency of electromagnetic simulation and break the limit of hardware capabilities, which may make some simulations unsolvable. Although more powerful computer configurations and better numerical algorithms can alleviate these challenges to a certain extent, they cannot completely solve the fundamental problem of escalating computing requirements as the scale of simulation continues to increase. Summary of the invention

[0003] In view of this, the present disclosure proposes a method and device for simplifying large-scale simulation, which uses the near-field scattering matrix of each substructure included in the target structure to be analyzed to determine the near-field scattering matrix of the complete target structure. This method of decomposing a large-scale simulation problem into multiple small-scale simulation problems can greatly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristic analysis task of complex structures.

[0004] According to one aspect of the present disclosure, a method for simplifying large-scale simulation is provided, comprising: obtaining an analysis task for a target structure, and performing structural division on the target structure to obtain a plurality of substructures; determining a propagation space required for executing the analysis task according to the analysis task, and determining a plurality of first modes of an electromagnetic wave in the propagation space according to a size of the propagation space, wherein the plurality of first modes include a plurality of propagation modes and a plurality of evanescent modes; performing electromagnetic simulation of each of the substructures using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each of the substructures, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure during near-field coupling; determining a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures, and using the target near-field scattering matrix as an analysis result of the analysis task, wherein the target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure during near-field coupling.

[0005] In this way, the near-field scattering matrix of each substructure included in the complex target structure can be used to accurately predict the near-field scattering matrix of the complete target structure. This method of decomposing large-scale simulation problems into multiple small-scale simulation problems can simplify the electromagnetic simulation process of complex target structures and significantly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristics analysis tasks of complex structures.

[0006] In a possible implementation, electromagnetic simulation of each of the substructures is performed using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure.

[0007] In this way, if the complete complex structure is divided into more substructures, the more computing time and memory can be saved. Whether parallel computing or serial computing, it can greatly save computing resources. The specific choice of parallel computing or serial computing can be made based on the trade-off between computing time and memory usage.

[0008] In a possible implementation, electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: inputting all the first modes of the electromagnetic wave into the substructure along a first direction respectively, to obtain a second mode of the electromagnetic wave formed on a reflection reference plane of the substructure corresponding to the first mode and a third mode of the electromagnetic wave formed on a transmission reference plane of the substructure corresponding to the first mode; determining a first reflection matrix corresponding to the substructure according to all the first modes and all the second modes, and determining a third mode of the electromagnetic wave corresponding to the substructure according to all the first modes and all the third modes. A first transmission matrix; and inputting all first modes of the electromagnetic wave into the substructure along a second direction opposite to the first direction, respectively, to obtain a fourth mode of the electromagnetic wave formed on the reflection reference plane corresponding to the first mode and a fifth mode of the electromagnetic wave formed on the transmission reference plane corresponding to the first mode; determining a second reflection matrix corresponding to the substructure according to all first modes and all fourth modes, and determining a second transmission matrix corresponding to the substructure according to all first modes and all fifth modes; determining a near-field scattering matrix of the substructure according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix.

[0009] In this way, calculating the substructure near-field scattering matrix of each substructure provides a basis for the subsequent accurate prediction of the near-field scattering matrix of the complex target structure.

[0010] In a possible implementation, a target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the substructures, including: using the first scattering matrix of each of the substructures to calculate the transmission matrix of the corresponding substructure to obtain a substructure transmission matrix of each substructure, each substructure transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the substructure; calculating a target transmission matrix based on the substructure transmission matrices of all substructures, the target transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the entire target structure; and calculating the target near-field scattering matrix based on the target transmission matrix.

[0011] In a possible implementation, the substructure near-field scattering matrix of each of the substructures includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; wherein, the first scattering matrix of each of the substructures is used to calculate the transmission matrix of the corresponding substructure to obtain the substructure transmission matrix of each of the substructures, including: performing a subtraction operation based on the first transmission matrix and the first result matrix to obtain a first submatrix, wherein the first result matrix is ​​obtained by multiplying the second reflection matrix, the inverse matrix of the second transmission matrix, and the first reflection matrix; performing a multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second submatrix; performing a multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third submatrix; using the inverse matrix of the second transmission matrix as the fourth submatrix; and determining the substructure transmission matrix of the substructure according to the first submatrix, the second submatrix, the third submatrix, and the fourth submatrix.

[0012] In a possible implementation, a target transmission matrix is ​​calculated based on substructure transmission matrices of all substructures, including: determining a calculation order based on the position of each substructure in the entire target structure; and multiplying all substructure transmission matrices according to the calculation order to obtain the target transmission matrix.

[0013] In a possible implementation, the target transmission matrix includes a first transmission sub-matrix, a second transmission sub-matrix, a third transmission sub-matrix, and a fourth transmission sub-matrix; wherein, the target near-field scattering matrix is ​​calculated according to the target transmission matrix, including: performing multiplication operations based on the inverse matrix of the fourth transmission sub-matrix and the third transmission sub-matrix to obtain a first matrix; using the inverse matrix of the fourth transmission sub-matrix as the second matrix; performing subtraction operations based on the first transmission sub-matrix and the second result matrix to obtain a third matrix, wherein the second result matrix is ​​obtained by multiplying the second transmission sub-matrix, the inverse matrix of the fourth transmission sub-matrix, and the third transmission sub-matrix; performing multiplication operations based on the second transmission sub-matrix and the inverse matrix of the fourth transmission sub-matrix to obtain a fourth matrix; and determining the target near-field scattering matrix according to the first matrix, the second matrix, the third matrix, and the fourth matrix.

[0014] According to another aspect of the present disclosure, a device for simplifying large-scale simulation is provided, including: an acquisition module, used to acquire an analysis task for a target structure, and structurally divide the target structure to obtain multiple substructures; a first determination module, used to determine the propagation space required for executing the analysis task according to the analysis task, and determine multiple first modes of the electromagnetic wave in the propagation space according to the size of the propagation space, wherein the multiple first modes include multiple propagation modes and multiple evanescent modes; a simulation module, used to perform electromagnetic simulation of each of the substructures using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each of the substructures, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure during the near-field coupling process; and a second determination module, used to determine a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures, and use the target near-field scattering matrix as the analysis result of the analysis task, wherein the target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure during the near-field coupling process.

[0015] In this way, the near-field scattering matrix of each substructure included in the complex target structure can be used to accurately predict the near-field scattering matrix of the complete target structure. This method of decomposing large-scale simulation problems into multiple small-scale simulation problems can simplify the electromagnetic simulation process of complex target structures and significantly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristics analysis tasks of complex structures.

[0016] In a possible implementation, electromagnetic simulation of each of the substructures is performed using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure.

[0017] In this way, if the complete complex structure is divided into more substructures, the more computing time and memory can be saved. Whether parallel computing or serial computing, it can greatly save computing resources. The specific choice of parallel computing or serial computing can be made based on the trade-off between computing time and memory usage.

[0018] In a possible implementation, electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: inputting all the first modes of the electromagnetic wave into the substructure along a first direction respectively, to obtain a second mode of the electromagnetic wave formed on a reflection reference plane of the substructure corresponding to the first mode and a third mode of the electromagnetic wave formed on a transmission reference plane of the substructure corresponding to the first mode; determining a first reflection matrix corresponding to the substructure according to all the first modes and all the second modes, and determining a third mode of the electromagnetic wave corresponding to the substructure according to all the first modes and all the third modes. A first transmission matrix; and inputting all first modes of the electromagnetic wave into the substructure along a second direction opposite to the first direction, respectively, to obtain a fourth mode of the electromagnetic wave formed on the reflection reference plane corresponding to the first mode and a fifth mode of the electromagnetic wave formed on the transmission reference plane corresponding to the first mode; determining a second reflection matrix corresponding to the substructure according to all first modes and all fourth modes, and determining a second transmission matrix corresponding to the substructure according to all first modes and all fifth modes; determining a near-field scattering matrix of the substructure according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix.

[0019] In this way, calculating the substructure near-field scattering matrix of each substructure provides a basis for the subsequent accurate prediction of the near-field scattering matrix of the complex target structure.

[0020] In a possible implementation, a target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the substructures, including: using the first scattering matrix of each of the substructures to calculate the transmission matrix of the corresponding substructure to obtain a substructure transmission matrix of each substructure, each substructure transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the substructure; calculating a target transmission matrix based on the substructure transmission matrices of all substructures, the target transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the entire target structure; and calculating the target near-field scattering matrix based on the target transmission matrix.

[0021] In a possible implementation, the substructure near-field scattering matrix of each of the substructures includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; wherein, the first scattering matrix of each of the substructures is used to calculate the transmission matrix of the corresponding substructure to obtain the substructure transmission matrix of each of the substructures, including: performing a subtraction operation based on the first transmission matrix and the first result matrix to obtain a first submatrix, wherein the first result matrix is ​​obtained by multiplying the second reflection matrix, the inverse matrix of the second transmission matrix, and the first reflection matrix; performing a multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second submatrix; performing a multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third submatrix; using the inverse matrix of the second transmission matrix as the fourth submatrix; and determining the substructure transmission matrix of the substructure according to the first submatrix, the second submatrix, the third submatrix, and the fourth submatrix.

[0022] In a possible implementation, a target transmission matrix is ​​calculated based on substructure transmission matrices of all substructures, including: determining a calculation order based on the position of each substructure in the entire target structure; and multiplying all substructure transmission matrices according to the calculation order to obtain the target transmission matrix.

[0023] In a possible implementation, the target transmission matrix includes a first transmission sub-matrix, a second transmission sub-matrix, a third transmission sub-matrix, and a fourth transmission sub-matrix; wherein, the target near-field scattering matrix is ​​calculated according to the target transmission matrix, including: performing multiplication operations based on the inverse matrix of the fourth transmission sub-matrix and the third transmission sub-matrix to obtain a first matrix; using the inverse matrix of the fourth transmission sub-matrix as the second matrix; performing subtraction operations based on the first transmission sub-matrix and the second result matrix to obtain a third matrix, wherein the second result matrix is ​​obtained by multiplying the second transmission sub-matrix, the inverse matrix of the fourth transmission sub-matrix, and the third transmission sub-matrix; performing multiplication operations based on the second transmission sub-matrix and the inverse matrix of the fourth transmission sub-matrix to obtain a fourth matrix; and determining the target near-field scattering matrix according to the first matrix, the second matrix, the third matrix, and the fourth matrix.

[0024] According to another aspect of the present disclosure, an electronic device is provided, comprising: a processor; and a memory for storing processor-executable instructions; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0025] According to another aspect of the present disclosure, a non-volatile computer-readable storage medium is provided, on which computer program instructions are stored, wherein the computer program instructions implement the above method when executed by a processor.

[0026] According to another aspect of the present disclosure, a computer program product is provided, including a computer-readable code, or a non-volatile computer-readable storage medium carrying the computer-readable code. When the computer-readable code runs in a processor of an electronic device, the processor in the electronic device executes the above method.

[0027] Further features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments with reference to the attached drawings.

[0028] The method and device for simplifying large-scale simulation provided by the embodiments of the present disclosure obtain an analysis task for a target structure, and perform structural division on the target structure to obtain a plurality of substructures; determine a propagation space required for executing the analysis task according to the analysis task, and determine a plurality of first modes of electromagnetic waves in the propagation space according to the size of the propagation space, wherein the plurality of first modes include a plurality of propagation modes and a plurality of evanescent modes; perform electromagnetic simulation of each substructure using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix for each substructure, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure in the near-field coupling process; based on the The target near-field scattering matrix of the target structure is determined based on the near-field scattering matrix of the substructure, and the target near-field scattering matrix is ​​used as the analysis result of the analysis task. The target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure in the near-field coupling process. In this way, the near-field scattering matrix of each substructure included in the target structure to be analyzed can accurately predict the near-field scattering matrix of the complete target structure. This method of decomposing large-scale simulation problems into multiple small-scale simulation problems can simplify the electromagnetic simulation process of complex structures and greatly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristic analysis tasks of complex structures. BRIEF DESCRIPTION OF THE DRAWINGS

[0029] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate exemplary embodiments, features, and aspects of the disclosure and, together with the description, serve to explain the principles of the disclosure.

[0030] Figures 1 to 4 A schematic diagram showing a method for simplifying large-scale simulation provided according to an embodiment of the present disclosure.

[0031] Figure 5 A schematic diagram of near-field coupling of a disordered metamaterial provided according to an embodiment of the present disclosure is shown.

[0032] Figure 6 A near-field scattering matrix and a near-field coupling schematic diagram of a chessboard-type metamaterial provided according to an embodiment of the present disclosure are shown.

[0033] Figure 7 A near-field scattering matrix and a near-field coupling schematic diagram of a metasurface provided according to an embodiment of the present disclosure are shown.

[0034] Figure 8 A schematic diagram showing near-field coupling of different types of metamaterials provided according to an embodiment of the present disclosure is shown.

[0035] Fig. 9 A schematic diagram of near-field coupling between a disordered metamaterial and air according to an embodiment of the present disclosure is shown.

[0036] Fig.10A block diagram of an apparatus for simplifying large-scale simulation according to an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0037] Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference numerals in the accompanying drawings represent elements with the same or similar functions. Although various aspects of the embodiments are shown in the accompanying drawings, the drawings are not necessarily drawn to scale unless otherwise specified.

[0038] The word “exemplary” is used exclusively herein to mean “serving as an example, example, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments.

[0039] In addition, in order to better illustrate the present disclosure, numerous specific details are given in the following specific embodiments. It should be understood by those skilled in the art that the present disclosure can also be implemented without certain specific details. In some examples, methods, means, components and circuits well known to those skilled in the art are not described in detail in order to highlight the subject matter of the present disclosure.

[0040] In order to facilitate those skilled in the art to understand the technical solution provided by the embodiments of the present disclosure, the technical environment in which the technical solution is implemented is described below.

[0041] The rapid development of metamaterials is closely related to the progress of full-wave simulation technology. Since most metamaterials lack analytical solutions, numerically solving Maxwell's equations through full-wave simulation is the most effective way to study metamaterials. The process involves discretizing the metamaterial structure in the solution space, analyzing the field distribution at each node in the grid, and post-processing the results to predict the scattering properties of the structure. With the support of electromagnetic wave simulation technology, more and more novel electromagnetic properties can be achieved through complex structural designs. However, with the development of metamaterials, their structural design has become more and more complex. Overall, the development trend is from simple periodic design to non-periodic design, from single layer to multilayer structure, and from ordered to disordered arrangement. Correspondingly, the scale of numerical simulation has been significantly expanded, which leads to higher memory requirements and longer solution time, bringing great challenges to electromagnetic simulation. As the complexity of metamaterial structures continues to increase, this trend may affect the efficiency of electromagnetic simulation and break the limit of hardware capabilities, which may make some simulations unsolvable. Although more powerful computer configurations and better numerical algorithms can alleviate these challenges to a certain extent, they cannot completely solve the fundamental problem of escalating computing requirements as the simulation scale continues to increase.

[0042] In order to solve the above technical problems, the embodiment of the present disclosure provides a method for simplifying large-scale simulation, obtaining an analysis task for a target structure, and structurally dividing the target structure to obtain multiple substructures; determining the propagation space required for executing the analysis task according to the analysis task, and determining multiple first modes of the electromagnetic wave in the propagation space according to the size of the propagation space, the multiple first modes including multiple propagation modes and multiple evanescent modes; using all the first modes of the electromagnetic wave to perform electromagnetic simulation of each substructure, to obtain a substructure near-field scattering matrix of each substructure, the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure during the near-field coupling process; determining a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all substructures, and using the target near-field scattering matrix as the analysis result of the analysis task, the target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure during the near-field coupling process. In this way, the near-field scattering matrix of each substructure included in the complex target structure can be used to accurately predict the near-field scattering matrix of the complete target structure. This method of decomposing large-scale simulation problems into multiple small-scale simulation problems can simplify the electromagnetic simulation process of complex target structures and significantly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristics analysis tasks of complex structures.

[0043] Now combined Figures 1 to 8 The method for simplifying large-scale simulation provided by the embodiment of the present disclosure is schematically illustrated by taking a metamaterial with a complex structure (referred to as target metamaterial) as the target structure. In addition to metamaterials, the method can also be applied to other types of complex structures, simplifying the large-scale simulation process of complex structures and saving computing resources.

[0044] like Figure 1 As shown, the method for simplifying large-scale simulation may include the following steps S101 to S104.

[0045] S101, obtaining an analysis task for a target structure, and dividing the target structure into multiple substructures.

[0046] From a geometric point of view, the complex large-scale structure in metamaterials can be regarded as the result of near-field interactions between simple structures. For example, the structure of the target metamaterial is Figure 2 The disordered complete metamaterial structure shown in (a) can be divided into four substructures, which are air 1, metamaterial A, metamaterial B, and air 2 on the right from left to right. The disclosed embodiment decomposes a complex system such as a metamaterial structure into multiple simple substructures for separate analysis. In this way, solving the analysis task of several simple structures is usually more efficient than solving a single complex structure. In addition, since the simple substructures are smaller and involve fewer grid points, the computational complexity is lower, which can significantly reduce the overall computational cost.

[0047] S102: Determine a propagation space required for performing the analysis task according to the analysis task, and determine a plurality of first modes of the electromagnetic wave in the propagation space according to the size of the propagation space.

[0048] The multiple first modes include multiple propagation modes and multiple evanescent modes. After confirming the specific size of the propagation space, such as the size of the cross section, the wavelength or frequency of the electromagnetic wave is determined, so that the number of propagation modes and evanescent modes of the electromagnetic wave at the frequency or wavelength in the propagation space at the specific size can be determined. Taking the metamaterial in the waveguide as an example, for a waveguide of a certain size, the number of propagation modes of the electromagnetic wave therein is limited, and the number of propagation modes is marked as N; although the number of evanescent modes is infinite in theory, once the near-field reference surface to be calculated is determined, the evanescent waves exceeding the near-field reference surface will be completely attenuated and will not affect the near-field interaction. Therefore, the number of evanescent wave modes (abbreviated as evanescent modes) is also limited, and the maximum number of evanescent modes is marked as N'. In addition to waveguides, the propagation space can also be free space. The main difference between free space and waveguides is that the physical meaning of different modes in free space scattering is different scattering directions, and the calculation of free space is similar to that of waveguides.

[0049] The effectiveness of this method of decomposing a large-scale simulation problem of a complex structure into multiple small-scale simulation problems provided by the disclosed embodiment depends on whether the scattering characteristics of a simple substructure can accurately predict the scattering characteristics of a complex system. For metamaterials, this depends on achieving effective near-field coupling and decoupling of metamaterials. However, current theoretical models, such as effective medium theory, transmission line theory, equivalent circuit method, multiple reflection model and coupled mode theory, cannot support near-field decoupling and recoupling. The unpredictability of near-field interactions usually stems from the complex and unknown resonances that occur between different metamaterial components during the near-field coupling process. The disclosed embodiment quantifies this resonance and studies the interaction of the complex target structure with the propagation mode electromagnetic wave and the evanescent electromagnetic wave for complete description. This scattering characteristic is called the near-field scattering characteristic of the target structure, and then the near-field scattering matrix of the substructure is used to predict the near-field scattering matrix of the complex target structure. For better explanation, the near-field scattering matrix of the substructure is referred to as the near-field scattering matrix of the substructure, and the near-field scattering matrix of the target structure is referred to as the target near-field scattering matrix.

[0050] S103, using all first modes of electromagnetic waves to perform electromagnetic simulation on each substructure, to obtain a substructure near-field scattering matrix of each substructure.

[0051] The substructure near-field scattering matrix can be used to indicate the scattering characteristics of the substructure during the near-field coupling process. The embodiment of the present disclosure obtains the near-field scattering matrix of the substructure through electromagnetic simulation to ensure that effective near-field coupling of the target structure can be achieved.

[0052] Step S103 may include: using all first modes of electromagnetic waves to perform electromagnetic simulations of each substructure in parallel or in sequence to obtain a substructure near-field scattering matrix of each substructure. The disclosed embodiment can decouple a complex metamaterial structure into a simple metamaterial structure. It has been verified that the computational resource consumption of computing all simple structures is much greater than the computational resource consumption of the complete structure. Figure 3 As shown, (a) is a schematic diagram of the decoupling of a large-scale complex structure into a small-scale simple structure. In this example, the complex structure is decoupled into two equal-sized structures and four equal-sized structures, and the time consumption and memory usage of calculating all structures are studied. In order to demonstrate the relationship between the scale of different grids and the consumption of computing resources during the full-wave simulation analysis process, the grid scale is adjusted by modifying the minimum grid size in the full-wave simulation software to demonstrate the calculation of large-scale structures. The calculation time and memory usage of the complete structure, the complete structure divided into two parts, and the complete structure divided into four parts are studied when the grid scale is different multiples of the number of basic grids. The specific comparison results of the calculation time can be found in Figure 3 For the specific comparison results of memory usage in (b), please refer to Figure 3 (c). The results show that, compared with directly simulating the complete structure, regardless of the scale of the meshing, dividing the complex structure into two equal-sized substructures can reduce the simulation time by 16% and memory usage by 56% on average when the substructures are run serially. When all substructures are studied in parallel, memory usage is reduced by about 7% and run time is reduced by about 58%. If the structure is divided into four equal-sized substructures, about 27% of time and 79% of memory can be saved when running serially, and 17% of memory usage and 82% of average run time can be reduced when operating in parallel. In this way, the more substructures a complete complex structure is divided into, the more computing time and memory can be saved. Whether it is parallel or serial operation, it can greatly save computing resources. The specific choice of parallel operation or serial operation can be weighed based on the calculation time and memory usage.

[0053] Using all the first modes of the electromagnetic wave to perform electromagnetic simulation of each substructure in parallel or in sequence to obtain the substructure near-field scattering matrix of each substructure may include: inputting all the first modes of the electromagnetic wave to the substructure along the first direction respectively, to obtain the second mode of the electromagnetic wave formed by the reflection reference surface of the substructure corresponding to the first mode and the third mode of the electromagnetic wave formed by the transmission reference surface of the substructure corresponding to the first mode, wherein the second mode and the third mode can be obtained simultaneously by inputting all the first modes once, or the second mode and the third mode can be obtained respectively by inputting all the first modes twice, and the specific method depends on the situation. The fourth mode and the fifth mode described below are similar to this, and the input of each first mode can be performed simultaneously or sequentially, and the specific method also depends on the situation. The embodiments of the present disclosure are for This is not limited; the first reflection matrix of the corresponding substructure is determined according to all the first modes and all the second modes, and the first transmission matrix of the corresponding substructure is determined according to all the first modes and all the third modes; and all the first modes of the electromagnetic wave are input to the substructure along the second direction opposite to the first direction, respectively, to obtain the fourth mode of the electromagnetic wave formed by the corresponding first mode on the reflection reference surface and the fifth mode of the electromagnetic wave formed by the corresponding first mode on the transmission reference surface; the second reflection matrix of the corresponding substructure is determined according to all the first modes and all the fourth modes, and the second transmission matrix of the corresponding substructure is determined according to all the first modes and all the fifth modes; the near-field scattering matrix of the substructure is determined according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix. The embodiment of the present disclosure performs electromagnetic simulation on the interaction of the propagation mode electromagnetic wave and the evanescent electromagnetic wave based on each substructure of the target structure to obtain the substructure near-field scattering matrix of each substructure, which provides a basis for the subsequent accurate prediction of the near-field scattering matrix of the complex target structure.

[0054] Still Figure 2 In (a), the four substructures of air 1, metamaterial A, metamaterial B, and air 2 on the left are shown as examples. The near-field scattering matrix of the substructure of air 1 can be recorded as S air1 , the substructure near-field scattering matrix of metamaterial A can be recorded as S A , the substructure near-field scattering matrix of metamaterial B can be recorded as S B , the near-field scattering matrix of the substructure of air 2 can be recorded as S ari2 , the substructure near-field scattering matrix of each substructure can be expressed by the following equation 1:

[0055]

[0056] In formula 1, S represents the substructure near-field scattering matrix, r represents the first reflection matrix, r′ represents the second reflection matrix, t represents the first transmission matrix, and t′ represents the second transmission matrix. The substructure near-field scattering matrix of each substructure includes the corresponding first reflection matrix, first transmission matrix, second reflection matrix, and second transmission matrix. For details, see Figure 2 It should be noted that all parameters shown in Formula 1 and the subsequent formulas in the embodiment of the present disclosure are in matrix form rather than numerical form used in related technologies, which helps to obtain more accurate near-field characteristic analysis results.

[0057] For a waveguide of fixed size, some of its eigenmodes are transmission modes that can be transmitted, and some are evanescent modes that cannot be transmitted. All elements in the near-field scattering matrix are calculated by electromagnetic full-wave simulation, that is, all first modes (such as Figure 4 The input shown is the scattering situation under the excitation source. Figure 4 For the values ​​of the i-th column in the first transmission matrix t shown, first set the i-th first mode as the excitation source and solve it through full-wave simulation. After the solution, the field distribution of the transmission reference surface will be obtained. The proportion of each second mode in the transmission field is obtained by modal decomposition of the transmission field distribution, and the proportions are filled into the i-th column of the first transmission matrix in order from low to high according to the modal order, where i is determined according to N and N′. That is to say, each column of the first transmission matrix represents the components (or coefficients) of each second mode contained in the transmission field after an electromagnetic wave of a first mode interacts with the substructure in the corresponding metamaterial. When all possible first modes are used as incident waves for excitation, the responses of all first mode electromagnetic waves and substructures are obtained (for example Figure 4 After the output shown in FIG. 1 is obtained, the complete first transmission matrix can be obtained. Similar operations are performed on the first reflection matrix, the second reflection matrix, and the second transmission matrix to further obtain the target near-field scattering matrix of the target structure.

[0058] S104. Determine a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all substructures, and use the target near-field scattering matrix as an analysis result of the analysis task.

[0059] Determining the target near-field scattering matrix of the target structure based on the near-field scattering matrices of all substructures in S104 may include: using the first scattering matrix of each substructure to calculate the transmission matrix of the corresponding substructure to obtain the substructure transmission matrix of each substructure, and each substructure transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the substructure; calculating the target transmission matrix based on the substructure transmission matrices of all substructures, and the target transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the entire target structure; and calculating the target near-field scattering matrix based on the target transmission matrix.

[0060] Among them, using the first scattering matrix of each substructure to calculate the transmission matrix of the corresponding substructure to obtain the substructure transmission matrix of each substructure may include: performing a subtraction operation based on the first transmission matrix and the first result matrix to obtain the first submatrix, wherein the first result matrix is ​​obtained by multiplying the second reflection matrix, the inverse matrix of the second transmission matrix, and the first reflection matrix; performing a multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain the second submatrix; performing a multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain the third submatrix; using the inverse matrix of the second transmission matrix as the fourth submatrix; and determining the substructure transmission matrix of the substructure based on the first submatrix, the second submatrix, the third submatrix, and the fourth submatrix. Now still with Figure 2 In (a), the four substructures of air 1, metamaterial A, metamaterial B, and air 2 on the left are used as examples to illustrate the calculation process of the substructure transmission matrix. The substructure transmission matrix of air 1 can be denoted as M air1 , the substructure transmission matrix of metamaterial A can be denoted as M A , the substructure transmission matrix of metamaterial B can be denoted as M B , the substructure transmission matrix of air 2 can be recorded as M air2 Each substructure transmission matrix of each substructure may include four submatrices, and the four submatrices are the first submatrix M 11 , the second sub-matrix M 12 , the third sub-matrix M 21 , the fourth sub-matrix M 22 For the substructure transmission matrix M of metamaterial A A , M A The first sub-matrix M included 11 The formula M 11 =tr′t′ -1 r determines (where r′t′ -1 r represents the first result matrix), M A The second sub-matrix M included 12 The formula M 12 = r′t′ -1 OK, M A The third sub-matrix M included 21 The formula M 21 = -t′ -1 r is confirmed, M A The fourth submatrix M included 22 The formula M 12 =t′ -1 Determine, where r is the substructure near-field scattering matrix S of metamaterial A A The first reflection matrix in , r′ is the substructure near-field scattering matrix S of metamaterial A AThe second reflection matrix in , t is the substructure near-field scattering matrix S of metamaterial A A The first transmission matrix in , t′ is the substructure near-field scattering matrix S of metamaterial A A The calculation process of the substructure transmission matrices of the remaining substructures is similar to that of metamaterial A, and will not be repeated in this paper for the sake of brevity.

[0061] Calculating the target transmission matrix based on the substructure transmission matrices of all substructures may include: determining the calculation order based on the position of each substructure in the entire target structure; performing multiplication operations on all substructure transmission matrices according to the calculation order to obtain the target transmission matrix. The target transmission matrix includes a first transmission submatrix, a second transmission submatrix, a third transmission submatrix, and a fourth transmission submatrix. Figure 2 In (a), the four substructures of air 1, metamaterial A, metamaterial B, and air 2 on the left are used as examples to illustrate the calculation process of the target transmission matrix. After the substructure transmission matrices of all substructures are obtained, the target transmission matrix of the complete target metamaterial can be calculated by the calculation sequence shown in the following formula 2:

[0062] M comp =M air2 *M B *M A *M air1 Formula 2

[0063] In formula 2, M comp represents the target transmission matrix, M air2 represents the substructure transmission matrix of air 2, M B represents the substructure transmission matrix of metamaterial B, M A represents the substructure transmission matrix of metamaterial A, M air1 represents the substructure transmission matrix of air 1. Similarly to the substructure transmission matrix, the target transmission matrix also includes four submatrices, which are the first transmission submatrix Mc 11 , the second transmission sub-matrix Mc 12 , the third transmission sub-matrix Mc 21 , the fourth transmission sub-matrix Mc 22 .

[0064] Among them, calculating the target near-field scattering matrix according to the target transmission matrix may include: performing multiplication operations based on the inverse matrix of the fourth transmission submatrix and the third transmission submatrix to obtain a first matrix; using the inverse matrix of the fourth transmission submatrix as the second matrix; performing subtraction operations based on the first transmission submatrix and the second result matrix to obtain a third matrix, wherein the second result matrix is ​​obtained by multiplying the second transmission submatrix, the inverse matrix of the fourth transmission submatrix, and the third transmission submatrix; performing multiplication operations based on the second transmission submatrix and the inverse matrix of the fourth transmission submatrix to obtain a fourth matrix; determining the target near-field scattering matrix according to the first matrix, the second matrix, the third matrix, and the fourth matrix. The target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure during the near-field coupling process. It is still based on Figure 2 (a) shows the four substructures of air 1 on the left, metamaterial A, metamaterial B, and air 2 on the right as examples to illustrate the calculation process of the target near-field scattering matrix. The first matrix r c The formula can be used to calculate r c =-Mc 22 -1 Mc 21 Determine, the second matrix t c ′ can be obtained by formula t c ′=Mc 22 -1 Determine, the third matrix t c The formula t c =Mc 11 -Mc 12 Mc 22 -1 Mc 21 Determine (where Mc 12 Mc 22 -1 Mc 21 represents the second result matrix), the fourth matrix r c ′ can be obtained by formula r c ′=Mc 12 Mc 22 -1 OK, among them, Mc 22 , Mc 21 , Mc 12 , Mc 11 See the above explanation. After the first matrix, the second matrix, the third matrix, and the fourth matrix are determined, the target near-field scattering matrix of the complete target metamaterial can be obtained by the following formula 3:

[0065]

[0066] In Equation 3, Sc represents the target near-field scattering matrix, r c represents the first matrix, t c′ represents the second matrix, t c represents the third matrix, r c ′ represents the fourth matrix. In fact, in the target near-field scattering matrix Sc, the first matrix r c It can be understood as a reflection matrix. The fourth matrix r c ′ can be understood as the reflection matrix corresponding to the electromagnetic wave transmitted in the reverse direction. The third matrix t c It can be understood as the transmission matrix, the second matrix t c ' can be understood as the transmission matrix corresponding to the electromagnetic wave transmitted in the reverse direction.

[0067] In this way, the target near-field scattering matrix of the entire target structure is calculated by integrating the substructure near-field scattering matrices of all substructures, thereby achieving a comprehensive prediction of the electromagnetic characteristics of the complex structure. Compared with the existing scheme, the simplified large-scale simulation method (or the characteristic analysis method for complex structures) provided in the embodiment of the present disclosure helps to accurately describe the propagation characteristics of electromagnetic waves inside the complex target structure. This method avoids directly simulating the complex target structure as a whole, but instead splits the complex target structure into multiple substructures, calculates their scattering matrices and transmission matrices respectively, and then integrates them to obtain the target near-field scattering matrix. This splitting-integration method significantly improves the efficiency and accuracy of the characteristic analysis of the complex target structure, and is particularly suitable for processing large-scale and complex structures to achieve near-field characteristic analysis tasks for complex structures.

[0068] Currently, the commonly used full-wave simulation methods include finite element method, time-domain finite difference method, frequency-domain finite difference method, etc. Here, the frequency-domain finite difference method is selected as the full-wave simulation to demonstrate its technical effect. Figure 5 Near-field coupling of the disordered metamaterial shown, Figure 5 (a) shows the structure of the disordered metamaterial in the waveguide. Figure 5 (b) shows the solution space and near-field reference cross section of disordered metamaterials. Figure 5 (c) and Figure 5 (d) shows two different disordered metamaterials and the near-field scattering matrix obtained based on the finite difference frequency domain (FDFD) method. Figure 5 (e) shows a schematic diagram of the structure after near-field coupling of two substructures (metamaterial A and metamaterial B), a scattering matrix calculated based on the frequency domain finite difference method, a scattering matrix after near-field coupling based on the scattering parameters of the substructures, and the difference between the two calculation methods. It can be determined that the calculation results obtained by the near-field decoupling scheme in the method provided in the embodiment of the present disclosure match the results obtained by direct simulation of the complete structure very well.

[0069] The method for simplifying large-scale simulation provided by the embodiments of the present disclosure is not only applicable to disordered metamaterials, but also to other types of metamaterials. Figure 6 The near-field scattering matrix and near-field coupling of the chessboard metamaterial are shown. Figure 7 The near-field scattering matrix and near-field coupling of the metasurface shown in FIG. The simplified large-scale simulation method provided by the embodiment of the present disclosure can also be applied to the near-field coupling of different types of metamaterials, such as Figure 8 The near-field coupling of three types of metamaterials, namely disordered metamaterials, chessboard metamaterials, and metasurfaces, is shown in FIG. Fig. 9 The process of converting near-field scattering into far-field scattering is essentially to couple the metamaterial and air in the near field.

[0070] The embodiment of the present disclosure also provides a device for simplifying large-scale simulation, including: an acquisition module, used to acquire an analysis task for a target structure, and structurally divide the target structure to obtain multiple substructures; a first determination module, used to determine the propagation space required for executing the analysis task according to the analysis task, and determine multiple first modes of the electromagnetic wave in the propagation space according to the size of the propagation space, wherein the multiple first modes include multiple propagation modes and multiple evanescent modes; a simulation module, used to perform electromagnetic simulation of each of the substructures using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each of the substructures, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure during the near-field coupling process; a second determination module, used to determine a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures, and use the target near-field scattering matrix as the analysis result of the analysis task, wherein the target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure during the near-field coupling process.

[0071] In a possible implementation, electromagnetic simulation of each of the substructures is performed using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure.

[0072] In a possible implementation, electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, including: inputting all the first modes of the electromagnetic wave into the substructure along a first direction respectively, to obtain a second mode of the electromagnetic wave formed on a reflection reference plane of the substructure corresponding to the first mode and a third mode of the electromagnetic wave formed on a transmission reference plane of the substructure corresponding to the first mode; determining a first reflection matrix corresponding to the substructure according to all the first modes and all the second modes, and determining a third mode of the electromagnetic wave corresponding to the substructure according to all the first modes and all the third modes. A first transmission matrix; and inputting all first modes of the electromagnetic wave into the substructure along a second direction opposite to the first direction, respectively, to obtain a fourth mode of the electromagnetic wave formed on the reflection reference plane corresponding to the first mode and a fifth mode of the electromagnetic wave formed on the transmission reference plane corresponding to the first mode; determining a second reflection matrix corresponding to the substructure according to all first modes and all fourth modes, and determining a second transmission matrix corresponding to the substructure according to all first modes and all fifth modes; determining a near-field scattering matrix of the substructure according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix.

[0073] In a possible implementation, a target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the substructures, including: using the first scattering matrix of each of the substructures to calculate the transmission matrix of the corresponding substructure to obtain a substructure transmission matrix of each substructure, each substructure transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the substructure; calculating a target transmission matrix based on the substructure transmission matrices of all substructures, the target transmission matrix is ​​used to indicate the propagation characteristics of electromagnetic waves in the entire target structure; and calculating the target near-field scattering matrix based on the target transmission matrix.

[0074] In a possible implementation, the substructure near-field scattering matrix of each of the substructures includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; wherein, the first scattering matrix of each of the substructures is used to calculate the transmission matrix of the corresponding substructure to obtain the substructure transmission matrix of each of the substructures, including: performing a subtraction operation based on the first transmission matrix and the first result matrix to obtain a first submatrix, wherein the first result matrix is ​​obtained by multiplying the second reflection matrix, the inverse matrix of the second transmission matrix, and the first reflection matrix; performing a multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second submatrix; performing a multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third submatrix; using the inverse matrix of the second transmission matrix as the fourth submatrix; and determining the substructure transmission matrix of the substructure according to the first submatrix, the second submatrix, the third submatrix, and the fourth submatrix.

[0075] In a possible implementation, a target transmission matrix is ​​calculated based on substructure transmission matrices of all substructures, including: determining a calculation order based on the position of each substructure in the entire target structure; and multiplying all substructure transmission matrices according to the calculation order to obtain the target transmission matrix.

[0076] In a possible implementation, the target transmission matrix includes a first transmission sub-matrix, a second transmission sub-matrix, a third transmission sub-matrix, and a fourth transmission sub-matrix; wherein, the target near-field scattering matrix is ​​calculated according to the target transmission matrix, including: performing multiplication operations based on the inverse matrix of the fourth transmission sub-matrix and the third transmission sub-matrix to obtain a first matrix; using the inverse matrix of the fourth transmission sub-matrix as the second matrix; performing subtraction operations based on the first transmission sub-matrix and the second result matrix to obtain a third matrix, wherein the second result matrix is ​​obtained by multiplying the second transmission sub-matrix, the inverse matrix of the fourth transmission sub-matrix, and the third transmission sub-matrix; performing multiplication operations based on the second transmission sub-matrix and the inverse matrix of the fourth transmission sub-matrix to obtain a fourth matrix; and determining the target near-field scattering matrix according to the first matrix, the second matrix, the third matrix, and the fourth matrix.

[0077] In some embodiments, the functions or modules included in the device provided by the embodiments of the present disclosure can be used to execute the method described in the above method embodiments. The specific implementation can refer to the description of the above method embodiments, and for the sake of brevity, it will not be repeated here.

[0078] The embodiment of the present disclosure also provides a computer-readable storage medium on which computer program instructions are stored, and the computer program instructions implement the above method when executed by a processor. The computer-readable storage medium can be a volatile or non-volatile computer-readable storage medium.

[0079] An embodiment of the present disclosure further proposes an electronic device, comprising: a processor; and a memory for storing instructions executable by the processor; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0080] The embodiments of the present disclosure also provide a computer program product, including a computer-readable code, or a non-volatile computer-readable storage medium carrying the computer-readable code. When the computer-readable code runs in a processor of an electronic device, the processor in the electronic device executes the above method.

[0081] Fig.10 1900 is a block diagram of a device for simplifying large-scale simulation according to an embodiment of the present disclosure. For example, the device 1900 may be provided as a server or a terminal device. Fig.10 , the apparatus 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932 for storing instructions, such as an application, that can be executed by the processing component 1922. The application stored in the memory 1932 may include one or more modules, each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute instructions to perform the above method.

[0082] The device 1900 may also include a power supply component 1926 configured to perform power management of the device 1900, a wired or wireless network interface 1950 configured to connect the device 1900 to a network, and an input / output interface 1958 (I / O interface). The device 1900 may operate based on an operating system stored in the memory 1932, such as Windows Server 2000. TM , MacOS X TM , Unix TM ,Linux TM , FreeBSD TM or similar.

[0083] In an exemplary embodiment, a non-volatile computer-readable storage medium is also provided, such as a memory 1932 including computer program instructions, which can be executed by the processing component 1922 of the device 1900 to perform the above method.

[0084] The present disclosure may be a system, a method and / or a computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

[0085] A computer-readable storage medium may be a tangible device that can hold and store instructions used by an instruction execution device. A computer-readable storage medium may be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing. More specific examples of computer-readable storage media (a non-exhaustive list) include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a static random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanical encoding device, such as a punch card or a raised structure in a groove on which instructions are stored, and any suitable combination of the foregoing. As used herein, a computer-readable storage medium is not to be interpreted as a transient signal per se, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagating through a waveguide or other transmission medium (e.g., a light pulse through a fiber optic cable), or an electrical signal transmitted through a wire.

[0086] The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, optical fiber transmissions, wireless transmissions, routers, firewalls, switches, gateway computers, and / or edge servers. The network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing / processing device.

[0087] The computer program instructions for performing the operation of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-related instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages, such as Smalltalk, C++, etc., and conventional procedural programming languages, such as "C" language or similar programming languages. Computer-readable program instructions may be executed completely on a user's computer, partially on a user's computer, as an independent software package, partially on a user's computer, partially on a remote computer, or completely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer via any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., using an Internet service provider to connect via the Internet). In some embodiments, an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), may be customized by utilizing the state information of the computer-readable program instructions, and the electronic circuit may execute the computer-readable program instructions, thereby realizing various aspects of the present disclosure.

[0088] Various aspects of the present disclosure are described herein with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer-readable program instructions.

[0089] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device that implements the functions / actions specified in one or more boxes in the flowchart and / or block diagram is generated. These computer-readable program instructions can also be stored in a computer-readable storage medium, and these instructions cause the computer, programmable data processing device, and / or other equipment to work in a specific manner, so that the computer-readable medium storing the instructions includes a manufactured product, which includes instructions for implementing various aspects of the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0090] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operating steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0091] The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system, method and computer program product according to multiple embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a part of a module, program segment or instruction, and a part of the module, program segment or instruction includes one or more executable instructions for realizing the specified logical function. In some alternative implementations, the function marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square boxes can actually be executed substantially in parallel, and they can sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of special hardware and computer instructions.

[0092] The embodiments of the present disclosure have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The selection of terms used herein is intended to best explain the principles of the embodiments, practical applications, or technical improvements in the market, or to enable other persons of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for simplifying large-scale simulation, characterized in that: include: Acquire an analysis task for a target structure, and perform structural division on the target structure to obtain a plurality of substructures; Determine a propagation space required for performing the analysis task according to the analysis task, and determine a plurality of first modes of the electromagnetic wave in the propagation space according to the size of the propagation space, wherein the plurality of first modes include a plurality of propagation modes and a plurality of evanescent modes; Performing electromagnetic simulation on each of the substructures using all first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix for each of the substructures, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure during near-field coupling; The target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the substructures, and the target near-field scattering matrix is ​​used as the analysis result of the analysis task. The target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure in the near-field coupling process.

2. The method according to claim 1, characterized in that Using all the first modes of the electromagnetic wave to perform electromagnetic simulation on each of the substructures, a substructure near-field scattering matrix of each of the substructures is obtained, including: Electromagnetic simulation of each of the substructures is performed in parallel or sequentially using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each of the substructures.

3. The method according to claim 2, characterized in that Using all the first modes of the electromagnetic wave to perform electromagnetic simulation of each of the substructures in parallel or in sequence to obtain a substructure near-field scattering matrix of each of the substructures, including: Inputting all first modes of the electromagnetic wave to the substructure along the first direction respectively, obtaining a second mode of the electromagnetic wave formed on a reflection reference plane of the substructure corresponding to the first mode and a third mode of the electromagnetic wave formed on a transmission reference plane of the substructure corresponding to the first mode; determining a first reflection matrix corresponding to the substructure according to all first modes and all second modes, and determining a first transmission matrix corresponding to the substructure according to all first modes and all third modes; and Inputting all the first modes of the electromagnetic wave to the substructure along a second direction opposite to the first direction, respectively, to obtain a fourth mode of the electromagnetic wave formed on the reflection reference plane corresponding to the first mode and a fifth mode of the electromagnetic wave formed on the transmission reference plane corresponding to the first mode; determining a second reflection matrix corresponding to the substructure according to all the first modes and all the fourth modes, and determining a second transmission matrix corresponding to the substructure according to all the first modes and all the fifth modes; The substructure near-field scattering matrix is ​​determined according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix.

4. The method according to claim 1, characterized in that: Determining a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures includes: Calculating the transmission matrix of the corresponding substructure using the first scattering matrix of each substructure to obtain the substructure transmission matrix of each substructure, wherein each substructure transmission matrix is ​​used to indicate the propagation characteristics of the electromagnetic wave in the substructure; Calculating a target transmission matrix based on substructure transmission matrices of all substructures, wherein the target transmission matrix is ​​used to indicate propagation characteristics of electromagnetic waves in the entire target structure; The target near-field scattering matrix is ​​calculated according to the target transmission matrix.

5. The method according to claim 4, characterized in that The substructure near-field scattering matrix of each of the substructures includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; The calculation of the transmission matrix of the corresponding substructure using the first scattering matrix of each substructure to obtain the substructure transmission matrix of each substructure includes: Performing a subtraction operation based on the first transmission matrix and the first result matrix to obtain a first submatrix, wherein the first result matrix is ​​obtained by performing a multiplication operation based on the second reflection matrix, the inverse matrix of the second transmission matrix, and the first reflection matrix; Perform multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second sub-matrix; Perform multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third sub-matrix; Using the inverse matrix of the second transmission matrix as the fourth sub-matrix; A substructure transmission matrix of the substructure is determined according to the first submatrix, the second submatrix, the third submatrix, and the fourth submatrix.

6. The method according to claim 5, characterized in that The target transfer matrix is ​​calculated based on the substructure transfer matrices of all substructures, including: Determining a calculation order based on the position of each of the substructures in the entire target structure; All substructure transmission matrices are multiplied according to the calculation order to obtain the target transmission matrix.

7. The method according to any one of claims 4 to 6, characterized in that The target transmission matrix includes a first transmission sub-matrix, a second transmission sub-matrix, a third transmission sub-matrix, and a fourth transmission sub-matrix; The step of calculating the target near-field scattering matrix according to the target transmission matrix includes: Performing a multiplication operation based on the inverse matrix of the fourth transmission submatrix and the third transmission submatrix to obtain a first matrix; Using the inverse matrix of the fourth transmission submatrix as the second matrix; Performing a subtraction operation based on the first transmission submatrix and the second result matrix to obtain a third matrix, wherein the second result matrix is ​​obtained by performing a multiplication operation based on the second transmission submatrix, the inverse matrix of the fourth transmission submatrix, and the third transmission submatrix; Performing a multiplication operation based on the second transmission submatrix and the inverse matrix of the fourth transmission submatrix to obtain a fourth matrix; The target near-field scattering matrix is ​​determined according to the first matrix, the second matrix, the third matrix, and the fourth matrix.

8. A device for simplifying large-scale simulation, characterized in that: include: An acquisition module is used to acquire an analysis task for a target structure and divide the target structure into multiple substructures; A first determination module is used to determine a propagation space required for performing the analysis task according to the analysis task, and determine a plurality of first modes of the electromagnetic wave in the propagation space according to a size of the propagation space, wherein the plurality of first modes include a plurality of propagation modes and a plurality of evanescent modes; A simulation module, used for performing electromagnetic simulation of each of the substructures using all the first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each of the substructures, wherein the substructure near-field scattering matrix is ​​used to indicate the scattering characteristics of the substructure in the near-field coupling process; The second determination module is used to determine the target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures, and use the target near-field scattering matrix as the analysis result of the analysis task. The target near-field scattering matrix is ​​used to indicate the scattering characteristics of the target structure in the near-field coupling process.

9. An electronic device, characterized in that: include: processor; a memory for storing processor-executable instructions; The processor is configured to implement the method described in any one of claims 1 to 7 when executing the instructions stored in the memory.

10. A non-volatile computer-readable storage medium having computer program instructions stored thereon, characterized in that: When the computer program instructions are executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images