Method and apparatus for simplifying large scale simulation

Abstract

The present disclosure relates to a method and device for simplifying large-scale simulation, the method comprising obtaining an analysis task for a target structure and performing structure partitioning on the target structure to obtain a plurality of sub-structures; determining a propagation space required for performing the analysis task according to the analysis task, and determining a plurality of first modes according to the size of the propagation space, the plurality of first modes comprising a plurality of propagation modes and a plurality of evanescent modes; performing electromagnetic simulation on each sub-structure using all the first modes of electromagnetic waves to obtain a sub-structure near-field scattering matrix of each sub-structure; and determining a target near-field scattering matrix of the target structure based on all the sub-structure near-field scattering matrices. The present disclosure accurately predicts the near-field scattering matrix of the complete target structure, and this way of decomposing a large-scale simulation problem into a plurality of small-scale simulation problems can simplify the electromagnetic simulation process of complex structures, greatly reduce the computing resource requirements of full-wave simulation, and thus cope with the near-field characteristic analysis task of complex structures.

Description

TECHNICAL FIELD

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

[0002] With the support of electromagnetic wave simulation technology, more and more novel electromagnetic properties can be achieved through complex structure design. However, as the structure becomes more and more complex, the scale of numerical simulation is significantly expanded accordingly, which leads to higher memory requirements and longer solution time, bringing great challenges to electromagnetic simulation. This trend may affect the efficiency of electromagnetic simulation and break the limit of hardware capacity, which may lead to certain simulations being unable to be solved. Although more powerful computer configurations and better numerical algorithms can alleviate these challenges to some extent, they cannot completely solve the fundamental problem of constantly upgrading computing requirements as the simulation scale continues to increase. SUMMARY

[0003] Therefore, the present disclosure proposes a method and device for simplifying large-scale simulation, which determines the near-field scattering matrix of the complete target structure by using the near-field scattering matrix of each substructure included in the target structure to be analyzed. This way of decomposing large-scale simulation problems 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 an aspect of the present disclosure, a method for simplifying large-scale simulation is provided, including: obtaining an analysis task for a target structure, and performing structure division on the target structure to obtain a plurality of substructures; determining a propagation space required for performing the analysis task according to the analysis task, and determining a plurality of first modes of electromagnetic waves in the propagation space according to the size of the propagation space, the plurality of first modes including a plurality of propagation modes and a plurality of evanescent modes; performing electromagnetic simulation of each of the substructures by using all the first modes of the electromagnetic waves, to obtain a substructure near-field scattering matrix of each of the substructures, the substructure near-field scattering matrix being used to indicate the scattering characteristics of the substructure in a near-field coupling process; determining a target near-field scattering matrix of the target structure based on all the substructure near-field scattering matrices, and taking the target near-field scattering matrix as an analysis result of the analysis task, the target near-field scattering matrix being used to indicate the scattering characteristics of the target structure in the near-field coupling process.

[0005] Thus, the near-field scattering matrix of the complete target structure can be accurately predicted by using the near-field scattering matrix of each substructure included in the complex target structure. This way of dividing a large-scale simulation problem into multiple small-scale simulation problems can simplify the electromagnetic simulation process of the complex target structure, greatly reduce the computational resource requirement of full-wave simulation, and thus cope with the near-field characteristic analysis task of the complex structure.

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

[0007] Thus, the more substructures the complete complex structure is divided into, the more computing time and memory can be saved. Whether parallel operation or serial operation is used, the computing resource can be greatly saved. Whether parallel operation or serial operation is used can be selected according to the computing time and memory occupation.

[0008] In a possible implementation, the electromagnetic simulation of each substructure is performed using all the first modes of the electromagnetic wave in parallel or sequentially to obtain the 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 the second mode of the electromagnetic wave formed by the reflection of the corresponding first mode on a reflection reference surface of the substructure and the third mode of the electromagnetic wave formed by the transmission of the corresponding first mode on a transmission reference surface of the substructure; determining the first reflection matrix corresponding to the substructure according to all the first modes and all the second modes, and determining the first transmission matrix corresponding to the substructure according to all the first modes and all the third modes; and inputting all the first modes of the electromagnetic wave into the substructure along a second direction opposite to the first direction respectively to obtain the fourth mode of the electromagnetic wave formed by the reflection of the corresponding first mode on the reflection reference surface and the fifth mode of the electromagnetic wave formed by the transmission of the corresponding first mode on the transmission reference surface; determining the second reflection matrix corresponding to the substructure according to all the first modes and all the fourth modes, and determining the second transmission matrix corresponding to the substructure according to all the first modes and all the fifth modes; and determining the substructure near-field scattering matrix according to the first reflection matrix, the first transmission matrix, the second reflection matrix, and the second transmission matrix.

[0009] Thus, the substructure near-field scattering matrix of each substructure is calculated, which provides a basis for accurately predicting the near-field scattering matrix of the complex target structure subsequently.

[0010] In a possible implementation, the target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the sub-structures, including: performing calculation of a transmission matrix of a corresponding sub-structure by using the first scattering matrix of each sub-structure to obtain a sub-structure transmission matrix of each sub-structure, each sub-structure transmission matrix being used to indicate a propagation characteristic of an electromagnetic wave in the sub-structure; calculating a target transmission matrix according to the sub-structure transmission matrices of all the sub-structures, the target transmission matrix being used to indicate a propagation characteristic of an electromagnetic wave in the entire target structure; and calculating the target near-field scattering matrix according to the target transmission matrix.

[0011] In a possible implementation, the sub-structure near-field scattering matrix of each sub-structure includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; wherein the calculation of the transmission matrix of the corresponding sub-structure by using the first scattering matrix of each sub-structure to obtain the sub-structure transmission matrix of each sub-structure includes: performing subtraction operation based on the first transmission matrix and a first result matrix to obtain a first sub-matrix, wherein the first result matrix is obtained by performing multiplication operation according to the second reflection matrix, an inverse matrix of the second transmission matrix, and the first reflection matrix; performing multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second sub-matrix; performing multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third sub-matrix; taking the inverse matrix of the second transmission matrix as a fourth sub-matrix; and determining the sub-structure transmission matrix of the sub-structure according to the first sub-matrix, the second sub-matrix, the third sub-matrix, and the fourth sub-matrix.

[0012] In a possible implementation, the target transmission matrix is calculated according to the sub-structure transmission matrices of all the sub-structures, including: determining a calculation sequence based on positions of each sub-structure in the entire target structure; and performing multiplication operation on all the sub-structure transmission matrices according to the calculation sequence to obtain the target transmission matrix.

[0013] In a possible implementation, the target transmission matrix comprises a first transmission sub-matrix, a second transmission sub-matrix, a third transmission sub-matrix, and a fourth transmission sub-matrix; and the target near-field scattering matrix is calculated according to the target transmission matrix, including: performing multiplication operation on the inverse matrix of the fourth transmission sub-matrix and the third transmission sub-matrix to obtain a first matrix; taking the inverse matrix of the fourth transmission sub-matrix as a second matrix; performing subtraction operation on the first transmission sub-matrix and a second resultant matrix to obtain a third matrix, where the second resultant matrix is obtained by performing multiplication operation on the second transmission sub-matrix, the inverse matrix of the fourth transmission sub-matrix, and the third transmission sub-matrix; performing multiplication operation 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, there is provided a device for simplifying large-scale simulation, comprising: an acquisition module configured to acquire an analysis task for a target structure, and perform structure partitioning on the target structure to obtain a plurality of sub-structures; a first determination module configured to determine a propagation space required for performing 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, the plurality of first modes comprising a plurality of propagation modes and a plurality of evanescent modes; a simulation module configured to perform electromagnetic simulation on each of the sub-structures by using all the first modes of the electromagnetic waves, to obtain a sub-structure near-field scattering matrix of each of the sub-structures, the sub-structure near-field scattering matrix being used to indicate the scattering characteristics of the sub-structure in a near-field coupling process; and a second determination module configured to determine a target near-field scattering matrix of the target structure based on all the sub-structure near-field scattering matrices, and take the target near-field scattering matrix as an analysis result of the analysis task, the target near-field scattering matrix being used to indicate the scattering characteristics of the target structure in the near-field coupling process.

[0015] In this way, the near-field scattering matrix of the complete target structure can be accurately predicted by using the near-field scattering matrix of each sub-structure included in the complex target structure, and the large-scale simulation problem can be divided into a plurality of small-scale simulation problems, so that the electromagnetic simulation process of the complex target structure is simplified, the calculation resource requirement of the full-wave simulation is greatly reduced, and the near-field characteristic analysis task of the complex structure is handled.

[0016] In one 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 of the substructures, including: performing electromagnetic simulation of each of the substructures 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.

[0017] In this way, the more substructures a complete complex structure is divided into, the greater the savings in computing time and memory. Both parallel and serial computing can significantly save computing resources. The choice between parallel and serial computing can be made based on the trade-off between computing time and memory usage.

[0018] In a possible implementation, all first modes of the electromagnetic wave are used to perform electromagnetic simulations of each substructure in parallel or sequentially to obtain a substructure near-field scattering matrix of each substructure, including: inputting all first modes of the electromagnetic wave into the substructure along a first direction, respectively, to obtain a second mode of the electromagnetic wave formed by the first mode on the reflection reference surface of the substructure and a third mode of the electromagnetic wave formed by the first mode on the transmission reference surface of the substructure; determining a first reflection matrix corresponding to the substructure according to all first modes and all second modes, and determining a third mode corresponding to the substructure according to all first modes and all 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 by the corresponding first mode on the reflection reference surface and a fifth mode of the electromagnetic wave formed by the corresponding first mode on the transmission reference surface; 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, the target near-field scattering matrix of the target structure is determined based on the near-field scattering matrices of all the sub-structures, including: performing calculation of a transmission matrix of a corresponding sub-structure by using the first scattering matrix of each sub-structure to obtain a sub-structure transmission matrix of each sub-structure, each sub-structure transmission matrix being used to indicate a propagation characteristic of an electromagnetic wave in the sub-structure; calculating a target transmission matrix according to the sub-structure transmission matrices of all the sub-structures, the target transmission matrix being used to indicate a propagation characteristic of an electromagnetic wave in the entire target structure; and calculating the target near-field scattering matrix according to the target transmission matrix.

[0021] In a possible implementation, the sub-structure near-field scattering matrix of each sub-structure includes a corresponding first reflection matrix, a first transmission matrix, a second reflection matrix, and a second transmission matrix; wherein the calculation of the transmission matrix of the corresponding sub-structure by using the first scattering matrix of each sub-structure to obtain the sub-structure transmission matrix of each sub-structure includes: performing subtraction operation based on the first transmission matrix and a first result matrix to obtain a first sub-matrix, wherein the first result matrix is obtained by performing multiplication operation according to the second reflection matrix, an inverse matrix of the second transmission matrix, and the first reflection matrix; performing multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second sub-matrix; performing multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third sub-matrix; taking the inverse matrix of the second transmission matrix as a fourth sub-matrix; and determining the sub-structure transmission matrix of the sub-structure according to the first sub-matrix, the second sub-matrix, the third sub-matrix, and the fourth sub-matrix.

[0022] In a possible implementation, the target transmission matrix is calculated according to the sub-structure transmission matrices of all the sub-structures, including: determining a calculation sequence based on positions of each sub-structure in the entire target structure; and performing multiplication operation on all the sub-structure transmission matrices according to the calculation sequence 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; and the target near-field scattering matrix is calculated according to the target transmission matrix, including: performing multiplication operation based on an inverse matrix of the fourth transmission sub-matrix and the third transmission sub-matrix to obtain a first matrix; taking the inverse matrix of the fourth transmission sub-matrix as a second matrix; performing subtraction operation based on the first transmission sub-matrix and a second result matrix to obtain a third matrix, where the second result matrix is obtained by performing multiplication operation based on the second transmission sub-matrix, the inverse matrix of the fourth transmission sub-matrix, and the third transmission sub-matrix; performing multiplication operation 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, including: a processor; 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 having stored thereon computer program instructions is provided, wherein the computer program instructions, when executed by a processor, implement the above method.

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

[0027] Other features and aspects of the present disclosure will become apparent from the following detailed description of example embodiments, taken in conjunction with the accompanying 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 structure division on the target structure to obtain a plurality of sub-structures; determine a propagation space required for performing 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 sub-structure by using all the first modes of electromagnetic waves, to obtain a sub-structure near-field scattering matrix of each sub-structure, wherein the sub-structure near-field scattering matrix is used to indicate the scattering characteristics of the sub-structure in a near-field coupling process; determine a target near-field scattering matrix of the target structure based on all the sub-structure near-field scattering matrices, and take 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 in the near-field coupling process, so that the near-field scattering matrix of the complete target structure can be accurately predicted by using the near-field scattering matrix of each sub-structure included in the target structure to be analyzed, and the electromagnetic simulation process of the complex structure can be simplified by decomposing the large-scale simulation problem into a plurality of small-scale simulation problems, and the computing resource requirement of the full-wave simulation is greatly reduced, so as to cope with the near-field characteristic analysis task of the complex structure. BRIEF DESCRIPTION OF DRAWINGS

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

[0030] Figures 1 to 4 A schematic diagram of a method for simplifying large-scale simulation is shown.

[0031] Figure 5 A schematic diagram of near-field coupling of disordered metamaterials is shown.

[0032] Figure 6 A schematic diagram of near-field scattering matrix and near-field coupling of chessboard-type metamaterials is shown.

[0033] Figure 7 A schematic diagram of near-field scattering matrix and near-field coupling of metasurfaces is shown.

[0034] Figure 8 A schematic diagram of near-field coupling of different kinds of metamaterials is shown.

[0035] Figure 9 A schematic diagram of near-field coupling of disordered metamaterials and air is shown.

[0036] Figure 10A block diagram illustrating an apparatus for simplifying large-scale simulation is shown. DETAILED DESCRIPTION

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

[0038] The word "exemplary" is used herein to mean "serving as an example, instance, or illustration." Any implementation described herein as "exemplary" is not necessarily to be construed as preferred or advantageous over other implementations.

[0039] In addition, for the purpose of convenience and brevity, detailed descriptions of well-known functions and structures incorporated in the present disclosure can be omitted. It will be appreciated that the present disclosure can be practiced with the exact details as other implementations can omit some of the described elements, not to complicate the present disclosure with the while.

[0040] In order to facilitate the understanding of the technical solutions provided by the embodiments of the present disclosure, the technical environment for implementing the technical solutions will be described first.

[0041] The rapid development of metamaterials is closely related to the progress of full-wave simulation techniques. As most metamaterials lack analytical solutions, numerically solving Maxwell's equations through full-wave simulation is the most effective method for studying metamaterials. This process involves discretizing the metamaterial structure within 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 techniques, more and more novel electromagnetic properties can be achieved through complex structural design. However, as metamaterials develop, their structural design becomes increasingly complex. Overall, the trend is from simple periodic design to non-periodic design, from single-layer to multi-layer structure, and from ordered to disordered arrangement. Correspondingly, the scale of numerical simulation has significantly expanded, which leads to higher memory requirements and longer solution times, posing great challenges to electromagnetic simulation. As the complexity of metamaterial structures continues to increase, this trend may affect the efficiency of electromagnetic simulation and exceed the limits of hardware capabilities, potentially leading to certain simulations that cannot be solved. Although more powerful computer configurations and better numerical algorithms can alleviate these challenges to some extent, they cannot completely solve the fundamental problem of escalating computing demands as the scale of simulation continues to increase.

[0042] In order to solve the above technical problems, an 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 to perform 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; based on the near-field scattering matrices of all substructures, a target near-field scattering matrix of the target structure is determined, 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 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 greatly reduce the computing resource requirements of full-wave simulation, thereby coping with the near-field characteristic analysis task of complex structures.

[0043] Now combined Figures 1 to 8 The method for simplifying large-scale simulation provided by the embodiments of the present disclosure is schematically illustrated using a metamaterial with a complex structure (referred to as the target metamaterial) as the target structure. In addition to metamaterials, this 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. Obtain an analysis task for a target structure, and divide the target structure into multiple substructures.

[0046] From a geometric perspective, the complex large-scale structures in metamaterials can be viewed 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: Air 1, Metamaterial A, Metamaterial B, and Air 2 on the right, from left to right. The disclosed embodiments decompose complex systems, such as metamaterial structures, into multiple simple substructures for separate analysis. Solving the analysis of several simple structures is often more efficient than solving a single complex structure. Furthermore, since simple substructures are smaller and involve fewer grid points, the computational complexity is lower, significantly reducing overall computational costs.

[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 dimensions of the propagation space, such as the dimensions 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 that frequency or wavelength in the propagation space of that specific dimension can be determined. Taking the metamaterial in a waveguide as an example, for a waveguide of a certain dimension, the number of propagation modes of the electromagnetic wave therein is finite, and the number of propagation modes is marked as N. Although the number of evanescent modes is theoretically infinite, once the near-field reference surface to be calculated is determined, the evanescent waves exceeding the near-field reference surface will completely attenuate and will not affect the near-field interaction. Therefore, the number of evanescent wave modes (abbreviated as evanescent modes) is also finite, 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. The calculation of free space is similar to that of waveguides.

[0049] The effectiveness of the method provided by the embodiments of the present disclosure to decompose the large-scale simulation problem of a complex structure into multiple small-scale simulation problems 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 embodiments of the present disclosure quantify this resonance and study the interaction of a complex target structure with propagating mode electromagnetic waves and evanescent electromagnetic waves to fully describe it. This scattering characteristic is called the near-field scattering characteristic of the target structure, and then the near-field scattering matrix of the complex target structure is predicted by the near-field scattering matrix of the substructure. For better explanation, the near-field scattering matrix of the substructure will be referred to as the substructure near-field scattering matrix, and the near-field scattering matrix of the target structure will be referred to as the target near-field scattering matrix.

[0050] S103 , performing electromagnetic simulation on each substructure using all first modes of electromagnetic waves to obtain a substructure near-field scattering matrix for 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. In the embodiment of the present disclosure, the near-field scattering matrix of the substructure is obtained 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 the 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 that 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, 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 respectively 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 In (b), the specific comparison results of memory usage can be found in Figure 3 (c) The results show that, compared to directly simulating the complete structure, regardless of the scale of the meshing, dividing the complex structure into two equal-sized substructures can reduce simulation time by an average of 16% and memory usage by 56% when the substructures are run serially. When all substructures are studied in parallel, memory usage is reduced by approximately 7% and runtime by approximately 58%. If the structure is divided into four equal-sized substructures, approximately 27% of time and 79% of memory can be saved when running serially, and 17% of memory usage and 82% of average runtime can be reduced when running in parallel. Thus, the more substructures a complete complex structure is divided into, the greater the computational time and memory savings can be. Both parallel and serial operations can significantly save computing resources. The choice of parallel or serial operation can be made based on the trade-off between computational time and memory usage.

[0053] Performing electromagnetic simulations of each substructure in parallel or sequentially using all first modes of electromagnetic waves to obtain a substructure near-field scattering matrix of each substructure may include: inputting all first modes of electromagnetic waves into the substructure along a first direction, respectively, to obtain a second mode of the electromagnetic wave formed by the first mode on the reflection reference surface of the substructure and a third mode of the electromagnetic wave formed by the first mode on the transmission reference surface of the substructure, wherein the second mode and the third mode may be obtained simultaneously by inputting all first modes once, or the second mode and the third mode may be obtained respectively by inputting all first modes twice, and the specific method depends on the situation. The fourth mode and the fifth mode described below are similarly applied hereto, and the input of each first mode may 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; a first reflection matrix of the corresponding substructure is determined based on all first modes and all second modes, and a first transmission matrix of the corresponding substructure is determined based on all first modes and all third modes; and all first modes of the electromagnetic wave are input to the substructure along a second direction opposite to the first direction, respectively, to obtain a fourth mode of the electromagnetic wave formed by the corresponding first mode on the reflection reference surface and a fifth mode of the electromagnetic wave formed by the corresponding first mode on the transmission reference surface; a second reflection matrix of the corresponding substructure is determined based on all first modes and all fourth modes, and a second transmission matrix of the corresponding substructure is determined based on all first modes and all fifth modes; a near-field scattering matrix of the substructure is determined based on 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 propagation mode electromagnetic waves and evanescent electromagnetic waves based on each substructure of the target structure to obtain the substructure near-field scattering matrix of each substructure, which provides a basis for subsequent accurate prediction of the near-field scattering matrix of complex target structures.

[0054] Still Figure 2 In (a), the four substructures of air 1 on the left, metamaterial A, metamaterial B, and air 2 on the right are shown as examples. The near-field scattering matrix of the substructure of air 1 can be recorded as S air1 , the near-field scattering matrix of the substructure of metamaterial A can be recorded as S A , the near-field scattering matrix of the substructure 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 formula 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 matrices rather than numerical values ​​used in related arts, 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 of the near-field scattering matrix are calculated by electromagnetic full-wave simulation, that is, the calculation of 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 by full-wave simulation. After the solution is completed, 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 of each second mode in the transmission field are filled in 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 the 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 FIG1 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 , 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 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 the electromagnetic wave 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 the electromagnetic wave in the entire target structure; and calculating the target near-field scattering matrix based on the target transmission matrix.

[0060] The calculation of the transmission matrix of the corresponding substructure by using the first scattering matrix of each substructure to obtain the substructure transmission matrix of each substructure can include: performing 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 multiplication operation according to the second reflection matrix, the inverse matrix of the second transmission matrix and the first reflection matrix; performing multiplication operation based on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second submatrix; performing multiplication operation based on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third submatrix; taking the inverse matrix of the second transmission matrix as a 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. Still taking the left air 1, the metamaterial A, the metamaterial B and the right air 2 as an example, the calculation process of the substructure transmission matrix is shown in (a). Figure 2 The substructure transmission matrix of the air 1 can be denoted as M air1 , the substructure transmission matrix of the metamaterial A can be denoted as M A , the substructure transmission matrix of the metamaterial B can be denoted as M B , and the substructure transmission matrix of the air 2 can be denoted as M air2 . The substructure transmission matrix of each substructure can include four submatrices, which are a first submatrix M 11 , a second submatrix M 12 , a third submatrix M 21 and a fourth submatrix M 22 . For the substructure transmission matrix M A of the metamaterial A, the first submatrix M 11 included in M A can be determined by the formula M 11 =t-r′t′ -1 r (wherein r′t′ -1 r represents the first result matrix), the second submatrix M 12 included in M A can be determined by the formula M 12 =r′t′ -1 , the third submatrix M 21 included in M A can be determined by the formula M 21 =-t′ -1 r, and the fourth submatrix M 22 included in M A can be determined by the formula M 12 =t′ -1 , wherein r is the first reflection matrix in the substructure near-field scattering matrix S A of the metamaterial A, and r′ is the first reflection matrix in the substructure near-field scattering matrix S 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. For the sake of brevity, this article will not go into details.

[0061] Calculating the target transmission matrix based on the substructure transmission matrices of all substructures may include: determining a calculation order based on the position of each substructure in the entire target structure; and 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 Figure (a) shows the four substructures (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 transmission matrix. Once the substructure transmission matrices of all substructures are obtained, the target transmission matrix of the complete target metamaterial can be calculated using the calculation sequence shown in Equation 2 below:

[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. Similar 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 in 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 an example 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 previous explanation. After determining the first matrix, the second matrix, the third matrix, and the fourth matrix, 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 cdenotes a second matrix, t c denotes a third matrix, r c denotes a fourth matrix. In fact, in the target near-field scattering matrix Sc, the first matrix r c can be understood as a reflection matrix, the fourth matrix r c ' can be understood as a reflection matrix corresponding to the electromagnetic wave in the opposite direction, the third matrix t c can be understood as a transmission matrix, the second matrix t c ' can be understood as a transmission matrix corresponding to the electromagnetic wave in the opposite 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, realizing 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 by the embodiments of the present disclosure helps to accurately describe the propagation characteristics of electromagnetic waves inside the complex target structure. The present method avoids directly simulating the entire complex target structure, but instead splits the complex target structure into multiple substructures, calculates the scattering matrix and transmission matrix of each substructure, and then integrates them to obtain the target near-field scattering matrix. This splitting-integrating method significantly improves the efficiency and accuracy of the characteristic analysis of complex target structures, and is particularly suitable for processing large-scale and complex structures to achieve near-field characteristic analysis tasks for complex structures.

[0068] Currently 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 effects. As shown in Figure 5 the near-field coupling of disordered metamaterials, Figure 5 Fig. (a) shows the structure of disordered metamaterials in a waveguide, Figure 5 Fig. (b) shows the solution space and near-field reference section of disordered metamaterials, Figure 5 Fig. (c) and Figure 5 Fig. (d) respectively show two different disordered metamaterials and the near-field scattering matrix obtained based on the frequency domain finite difference (FDFD) method, Figure 5 Fig. (e) shows the structure of two substructures (metamaterial A and metamaterial B) after near-field coupling and the scattering matrix calculated based on the frequency domain finite difference method, the scattering matrix after near-field coupling based on the substructure scattering parameters, 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 by the embodiments of the present disclosure match the results obtained by directly simulating 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 applicable to other types of metamaterials. For details, refer to Figure 6 The near-field scattering matrix and near-field coupling of the chessboard-type metamaterial are shown, Figure 7 The near-field scattering matrix and near-field coupling of the super surface are shown. The method for simplifying large-scale simulation provided by the embodiments of the present disclosure is also applicable to the near-field coupling of different types of metamaterials, such as Figure 8 The near-field coupling of the three types of metamaterials, i.e., disordered metamaterials, chessboard-type metamaterials, and super surfaces, is shown. For example, Figure 9 The process of converting the near-field scattering into far-field scattering is shown. The essence is to perform near-field coupling between the metamaterial and air.

[0070] The embodiments of the present disclosure also provide a device for simplifying large-scale simulation, including: an acquisition module configured to acquire an analysis task for a target structure, and perform structure division on the target structure to obtain a plurality of sub-structures; a first determination module configured to determine a propagation space required for performing 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, the plurality of first modes including a plurality of propagation modes and a plurality of evanescent modes; a simulation module configured to perform electromagnetic simulation on each of the sub-structures by using all the first modes of the electromagnetic waves, to obtain a sub-structure near-field scattering matrix of each of the sub-structures, the sub-structure near-field scattering matrix being used to indicate scattering characteristics of the sub-structure in a near-field coupling process; and a second determination module configured to determine a target near-field scattering matrix of the target structure based on all the sub-structure near-field scattering matrices, and take the target near-field scattering matrix as an analysis result of the analysis task, the target near-field scattering matrix being used to indicate scattering characteristics of the target structure in the near-field coupling process.

[0071] In a possible implementation, the electromagnetic simulation on each of the sub-structures by using all the first modes of the electromagnetic waves to obtain the sub-structure near-field scattering matrix of each of the sub-structures includes: performing the electromagnetic simulation on each of the sub-structures by using all the first modes of the electromagnetic waves in parallel or sequentially to obtain the sub-structure near-field scattering matrix of each of the sub-structures.

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

[0073] In a possible implementation, the target near-field scattering matrix of the target structure is determined based on all the sub-structure near-field scattering matrices, including: performing calculation of a transmission matrix of each of the sub-structures using the first scattering matrix of the sub-structure, to obtain a sub-structure transmission matrix of each of the sub-structures, each of the sub-structure transmission matrices being used to indicate a propagation characteristic of the electromagnetic wave in the sub-structure; calculating a target transmission matrix according to the sub-structure transmission matrices of all the sub-structures, the target transmission matrix being used to indicate a propagation characteristic of the electromagnetic wave in the entire target structure; and calculating the target near-field scattering matrix according to 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; and the calculation of the transmission matrix of the corresponding substructure by using the first scattering matrix of each of the substructures to obtain the substructure transmission matrix of each of the substructures includes: performing subtraction operation on the first transmission matrix and a first result matrix to obtain a first submatrix, wherein the first result matrix is obtained by performing multiplication operation on the second reflection matrix, an inverse matrix of the second transmission matrix, and the first reflection matrix; performing multiplication operation on the second reflection matrix and the inverse matrix of the second transmission matrix to obtain a second submatrix; performing multiplication operation on the inverse matrix of the second transmission matrix and the first reflection matrix to obtain a third submatrix; taking the inverse matrix of the second transmission matrix as a 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, the target transmission matrix is calculated according to the substructure transmission matrices of all the substructures, including: determining a calculation sequence based on the positions of the substructures in the entire target structure; and performing multiplication operation on all the substructure transmission matrices according to the calculation sequence to obtain the target transmission matrix.

[0076] In a possible implementation, the target transmission matrix includes a first transmission submatrix, a second transmission submatrix, a third transmission submatrix, and a fourth transmission submatrix; and the target near-field scattering matrix is calculated according to the target transmission matrix, including: performing multiplication operation on the inverse matrix of the fourth transmission submatrix and the third transmission submatrix to obtain a first matrix; taking the inverse matrix of the fourth transmission submatrix as a second matrix; performing subtraction operation on the first transmission submatrix and a second result matrix to obtain a third matrix, wherein the second result matrix is obtained by performing multiplication operation on the second transmission submatrix, the inverse matrix of the fourth transmission submatrix, and the third transmission submatrix; performing multiplication operation on the second transmission submatrix and the inverse matrix of the fourth transmission submatrix 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 apparatus provided by the embodiments of the present disclosure has functions or includes modules that can be used to perform the methods described in the above method embodiments, and the specific implementation can be referred to the description of the above method embodiments. For briefness, details are not described herein.

[0078] The embodiment of the present disclosure further provides a computer readable storage medium, which has computer program instructions stored thereon. The computer program instructions are executed by a processor to implement the method described above. The computer readable storage medium can be a volatile or non-volatile computer readable storage medium.

[0079] The embodiment of the present disclosure further provides an electronic device, which comprises a processor, and a memory for storing processor-executable instructions. The processor is configured to implement the method described above when executing the instructions stored in the memory.

[0080] The embodiment of the present disclosure further provides a computer program product, which comprises computer readable codes or a non-volatile computer readable storage medium carrying the computer readable codes. When the computer readable codes are run in a processor of an electronic device, the processor in the electronic device implements the method described above.

[0081] Figure 10 A block diagram of an apparatus for simplifying large-scale simulation is shown according to an embodiment of the present disclosure. For example, the apparatus 1900 can be provided as a server or a terminal device. Referring to Figure 10 , the apparatus 1900 comprises a processing component 1922, which further comprises one or more processors, and a memory resource represented by a memory 1932, for storing instructions executable by the processing component 1922, such as an application program. The application program stored in the memory 1932 can comprise one or more than one module each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute the instructions to perform the method described above.

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

[0083] In an exemplary embodiment, a non-volatile computer readable storage medium is also provided, such as the memory 1932 comprising computer program instructions executable by the processing component 1922 of the apparatus 1900 to complete the method described above.

[0084] The present disclosure can be a system, a method, and / or a computer program product. The computer program product can include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present disclosure.

[0085] The computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device. The computer readable storage medium can be, for example, but is not limited to, an electronic 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 (a non-exhaustive list) of the computer readable storage medium include the following: a portable computer diskette, 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 disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device such as punch-cards or punched tape, a

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

[0087] Computer readable program instructions for carrying out operations of the present disclosure can be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Smalltalk, C++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The computer readable program instructions can execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the latter scenario, the remote computer can be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection can be made to an external computer (for example, through the Internet using an Internet Service Provider). In some embodiments, electronic circuitry including, for example, programmable logic circuitry, field-programmable gate array (FPGA), or programmable logic array (PLA) can execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present disclosure.

[0088] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0089] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide processes for implementing the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0090] The computer readable program instructions can also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus or other device to produce a computer implemented process such that the instructions which execute on the computer, other programmable data processing apparatus, or other device implement the functions / acts specified in the flowchart and / or block diagram block or blocks.

[0091] The flow diagrams and the block diagrams in the drawings are presented to illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flow diagrams and the block diagrams can represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logic functions. In some alternative implementations, the functions noted in the blocks can occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks can sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and / or flow diagrams, and combinations thereof, can be implemented by special purpose hardware-based systems that perform the specified functions or acts, or combinations of special purpose hardware and

[0092] Embodiments of the present disclosure have been described above, and the description is intended to be illustrative of the embodiments and not restrictive. Many modifications and variations of the described embodiments are possible and are within the scope of the disclosure. The selection of terms is intended to best describe the principles of the embodiments, practical application, or technical improvements in the art, or to enable others skilled in the art to understand the embodiments disclosed herein.

Claims

1. A method for simplifying large-scale simulation, characterized in that: include: Acquiring 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 performing the analysis task according to the analysis task, and determining a plurality of first modes of the electromagnetic wave in the propagation space according to a size of the propagation space, the plurality of first modes including 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 substructure, wherein the substructure near-field scattering matrix is ​​used to indicate a scattering characteristic of the substructure during near-field coupling; A 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 during the near-field coupling process.

2. The method according to claim 1, characterized in that 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 includes: Electromagnetic simulations of the substructures are performed in parallel or sequentially using all first modes of the electromagnetic waves to obtain a substructure near-field scattering matrix of each substructure.

3. The method according to claim 2, characterized in that Performing electromagnetic simulations of the substructures in parallel or sequentially using all first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure includes: Inputting all 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 based on all first modes and all second modes, and determining a first transmission matrix corresponding to the substructure based on all first modes and all third modes; 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 based on all first modes and all fourth modes, and determining a second transmission matrix corresponding to the substructure based on all first modes and all 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, wherein Determining a target near-field scattering matrix of the target structure based on the near-field scattering matrices of all the substructures includes: Calculating a transmission matrix of the corresponding substructure using the first scattering matrix of each substructure to obtain a substructure transmission matrix of each substructure, wherein each substructure transmission matrix is ​​used to indicate propagation characteristics of electromagnetic waves 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 substructure 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 a 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 sub-matrix; Performing a 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 a 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 substructure in the entire target structure; Perform multiplication operations on all substructure transfer matrices according to the calculation order to obtain the target transfer 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; 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 multiplying 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 sub-matrix and the inverse matrix of the fourth transmission sub-matrix 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 determining module, configured to determine, according to the analysis task, a propagation space required for performing the analysis task, and determine, according to a size of the propagation space, a plurality of first modes of the electromagnetic wave in the propagation space, the plurality of first modes including a plurality of propagation modes and a plurality of evanescent modes; a simulation module, configured to perform electromagnetic simulation of each of the substructures using all first modes of the electromagnetic wave to obtain a substructure near-field scattering matrix of each substructure, wherein the substructure near-field scattering matrix is ​​used to indicate a scattering characteristic of the substructure during a 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 during 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 according to 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