57 results about "Computational electromagnetics" patented technology

The invention discloses an ocean-ionized layer ultralow-frequency electromagnetic target detection simulation method based on ADI-FDTD and SIBC, and belongs to the field of computational electromagnetics. The ocean-ionized layer ultralow-frequency electromagnetic target detection simulation method comprises the following steps: introducing a frequency domain surface impedance boundary condition SIBC into an ADI-FDTD model; the ocean-ionosphere simulation model is used for carrying out ocean-ionosphere ultra-low frequency electromagnetic wave attenuation analysis on the ocean-ionosphere ultra-low frequency electromagnetic wave attenuation. According to the method, the frequency domain surface impedance boundary condition (SIBC) is combined with the ADI-FDTD method for the first time, so that the calculation space for analyzing the influence of the seabed on electromagnetic wave propagation is reduced. The calculation requirement is obviously reduced, and the efficiency is improved. In addition, the SIBC model based on the ADI-FDTD method can be used for simulating the scattering effect of electromagnetic waves on a high-conductivity and low-conductivity area interface.

The invention belongs to the technical field of computational electromagnetics, and discloses a numerical calculation method for simulating electromagnetic wave propagation characteristics in time-varying plasma, which is a three-dimensional Newton-HIE-FDTD method for efficiently simulating propagation of electromagnetic waves in time-varying plasma based on a hybrid explicit-implicit finite difference time domain method and Newton's second law of motion. The method inherits the advantages of the Newton-ADE-FDTD method, the time-varying electron density is accurately represented through the time derivative, and the abundant dynamic characteristics of the time-varying plasma can be accurately captured. On the basis, the method also fuses the advantages of an HIE-FDTD method, so that the time step length of the algorithm is only restricted by the sizes of grids in two directions, and compared with other simulation methods based on time-varying plasma, the method provided by the invention is more excellent in calculation efficiency and is more suitable for large-scale popularization and application. The method is especially suitable for analysis and simulation of electromagnetic wave propagation characteristics in time-varying plasma.

The invention discloses a target echo classification method and system based on multi-modal data fusion, relates to the technical field of computational electromagnetics, and solves the technical problems that the classification accuracy of a model is reduced in a noise environment, the generalization ability is limited, and the performance fluctuation is large when modal data distribution changes. The time sequence dependency is modeled through the Bi LSTM, the weight of each feature is analyzed, after the features are removed, global information is reserved through cross-modal fusion without removing the features, key discrimination features are focused through cross fusion after the features are removed, and the information integrity and the calculation efficiency can be both considered through combination of the two. And dynamically adjusting the contribution of different modal features to a classification result through the hidden state of the GRU. If certain modal data is missing, the original fusion path can still utilize other modal information, and the fusion path screening can depend on the key features of the remaining modals. And the sensitivity of the model to a specific noise mode can be reduced through the balance of redundant feature reservation and key feature screening.

The invention aims to provide an electromagnetic parameter extraction method based on a non-conformal partial element equivalent circuit method, and belongs to the technical field of computational electromagnetism. According to the method, the multi-scale electromagnetic structure is divided into different areas according to the size difference, each area is subdivided by adopting a triangular prism grid with the corresponding size, and meanwhile, an h-SPV primary function is defined at the non-conformal junction, so that the continuity of current at the non-conformal junction is ensured. According to the method, in the calculation process of the multi-scale structure problem, a large number of transition grid units at the interfaces of different areas are avoided through non-conformal grid subdivision, and the calculation complexity is reduced while the calculation precision is guaranteed by establishing the current continuity boundary conditions at the interface ports.

The present invention belongs to the field of computational electromagnetics and provides a multi-state simulation method for a reconfigurable intelligent surface to achieve the purpose of shortening the calculation time and reducing the memory consumption. First, obtain the unit model of the reconfigurable intelligent surface, expand the unit model to form a 3×3 characteristic sub-unit basic array, and regard it as an invariant geometric structure item, and regard the two-dimensional equivalent circuit model as a variable state description item. Then, extract the characteristics of the invariant geometric structure item and the variable state description item to obtain the system matrix of the reconfigurable intelligent surface and construct a system equation. Finally, solve the system equation to obtain the unknown electromagnetic field quantities to be solved, so as to calculate the electromagnetic characteristics of the reconfigurable intelligent surface and complete the multi-state simulation of the reconfigurable intelligent surface. The present invention effectively optimizes the repeated calculation part in the multi-state simulation process of the reconfigurable intelligent surface, significantly reduces the memory consumption, and greatly shortens the simulation time under the condition of ensuring the solution accuracy.

The application provides a method for calculating interface polarization charge and current distribution based on reflection coefficient, and relates to the technical field of computational electromagnetics. The method comprises the following steps: establishing a reflection coefficient-dielectric parameter mapping, directly correlating the reflection characteristics of electromagnetic waves with the intrinsic parameters of the medium through a reflection coefficient formula, and laying a physical foundation for polarization analysis; interface electromagnetic field reconstruction, using boundary conditions and constitutive relations to determine the electric field distribution on both sides of the interface and the corresponding polarization intensity. The distribution rule of the electromagnetic field in the medium is determined, the macroscopic electric field is converted into the spatial description of the polarization intensity, and an intermediate variable is provided for charge and current calculation; polarization dynamic response quantification, through the spatial and temporal derivatives of the polarization intensity, the interface polarization charge density and the polarization current density are derived. The theoretical model is converted into actual physical quantities, and the dynamic response mechanism of charges and currents when electromagnetic waves interact with the medium is revealed.

The application discloses a kind of electromagnetic response prediction methods of nuclear explosion pulse excitation under power plant rotating equipment, it is related to the technical field of computational electromagnetics.Nuclear explosion pulse signal is decomposed, and spectrum signal under multiple frequencies is obtained;According to the spectrum signal, generate the electric field excitation source of power plant rotating equipment at each frequency;Establish the multi-conductor transmission line model of power plant rotating equipment;At each frequency, the electric field excitation source of frequency is regarded as input, and the multi-conductor transmission line model is solved, and the electromagnetic response of power plant rotating equipment at frequency is obtained;According to the electromagnetic response at each frequency, determine the electromagnetic response of power plant rotating equipment, and the electromagnetic response of power plant rotating equipment under nuclear explosion pulse excitation can be accurately predicted by the method.

The invention discloses a triangular mesh-based high-order electromagnetic field simulation method, and belongs to the field of computational electromagnetism. In order to solve the problems of significant numerical dispersion and phase error accumulation when an existing triangular mesh time domain finite difference method is used for processing a complex geometric structure, the method comprises the following steps: firstly, introducing a first additional loop on a basic triangular mesh, and combining with a first optimization coefficient to construct a high-order model; and a second additional loop is further introduced, and an optimization model is constructed in combination with a second optimization coefficient, so that the phase error is effectively inhibited under the coarse mesh, and high-precision electromagnetic simulation of a large-scale polygon structure is realized.

The invention relates to the technical field of computational electromagnetism, in particular to a metasurface scattering behavior analysis method and device, equipment and a storage medium. The method comprises the following steps: tracking an incident wave emitted into the metasurface based on geometrical optics; for the tracked incident wave, determining a first relationship between the equivalent surface current and the scattered field based on physical optics; determining an equivalent surface current based on the generalized boundary condition and a surface polarizability model; and determining a scattering field by combining the equivalent surface current and the first relationship by adopting Gordon integral. According to the metasurface scattering behavior analysis method provided by the embodiment of the invention, the surface polarizability model is combined with physical optics and geometrical optics, so that the scattering behavior related to the metasurface target can be effectively analyzed, and meanwhile, the calculation amount is reduced.

The embodiment of the invention provides an electromagnetic scattering calculation method based on a heterogeneous GPU cluster. The method is applied to the field of computational electromagnetism, and comprises the following steps: constructing a hierarchical bounding volume geometric acceleration structure of a target object and adding an edge index, and meanwhile, completely copying and distributing data of the acceleration structure to a video memory of each computational node in a heterogeneous GPU cluster; generating a random ray path starting from an emission source, modeling an electromagnetic scattering process as a ray path set in a probability space, and generating a ray direction through an importance sampling strategy; determining a static scheduling strategy for the ray batch based on the calculation cost of pre-sampling and geometric enhancement, and carrying out non-uniform distribution on task tiles according to the heterogeneous GPU calculation power; and monitoring the execution state of the heterogeneous GPU cluster in real time, and dynamically adjusting the distribution strategy of the residual ray tasks based on execution feedback. According to the method, the calculation overhead can be remarkably reduced while the calculation precision is ensured, and the efficiency and expandability of complex target electromagnetic scattering calculation are improved.

The present invention discloses a multi-physics numerical calculation method for an on-chip device loaded with a ferrite thin film, belonging to the fields of integrated circuit technology and computational electromagnetics. First, a hybrid grid meshing is performed on the on-chip device loaded with the ferrite thin film; then, each field component of the three physical fields is arranged in the discrete space; finally, according to the coupling relationship of the three physical fields, the physical field equations are solved and the parameter transfer between the equations is carried out. According to the structural characteristics of the on-chip device, the present invention combines two grid encryption technologies in the conventional FDTD space, significantly improving the modeling accuracy; according to the multi-physics phenomena faced by the on-chip device, on the basis of solving the Maxwell equation by FDTD, the LLG equation and the heat conduction equation are coupled and calculated, and reasonable parameter transfer is carried out, completely restoring the electromagnetic field phenomena, the magnetization and demagnetization processes of the ferrite thin film and their thermal effects when the on-chip device is working, and effectively improving the simulation accuracy of the calculation method for the target.

This invention belongs to the field of computational electromagnetics numerical simulation and provides a simulation method for electromagnetic scattering of thin-coated targets based on Schur complement preconditions. First, for the IBC scattering target to be solved, equivalent surface current and equivalent surface magnetic current are introduced as independent unknowns. The electromagnetic scattering solution equation is established using SDIE and then discretized. Next, the impedance matrix is ​​decomposed into near-field and far-field matrices using an octree structure. For the near-field matrix, key operator sub-blocks are explicitly assembled and extracted. Then, approximate operators in block diagonal form and Schur complement operators are constructed and approximately inverted. Simultaneously, the diagonal product of intermediate blocks is pre-calculated and stored based on the approximate operators. Finally, the constructed preconditions are applied to the GMRES iterative solver. When the iteration error reaches the preset convergence criterion, the current and magnetic current distributions on the target surface are output. Further calculations of radiation characteristic parameters such as the target's radar cross section are then performed, achieving low-cost, high-efficiency, and easily parallelizable electromagnetic simulation solutions.

The invention belongs to the field of computational electromagnetism numerical simulation, and particularly relates to a non-uniform conformal FDTD method based on a sub-grid technology. According to the method, non-uniform conformal geometric correction is introduced into the sub-grids, so that high-precision fitting of a curved surface boundary is realized, step approximation errors are reduced, and the calculation precision of field quantities near the boundary is improved; according to the method, bidirectional coupling transfer which simultaneously meets line integral conservation and flux conservation is established on a main-sub grid interface, and a conservation relation required by a non-uniform conformal FDTD method is prevented from being damaged, so that interface numerical value reflection is effectively inhibited, energy non-physical errors are reduced, and long-term stability is improved; and fine grids are only adopted in the local fine structure neighborhood and are matched with time sub-cycle propulsion, so that global grid encryption and limitation of global time step length to a minimum unit are avoided, and the storage and calculation time overhead is remarkably reduced while the precision is ensured. According to the method, geometric modeling precision, interface stability and calculation efficiency are considered, and the method has high engineering application value.

The invention discloses a chaff cloud cluster electromagnetic scattering calculation method based on a binary-voxel method, and belongs to the technical field of computational electromagnetism and electronic countermeasures. According to the method, an incidence voxel coordinate system and a scattering voxel coordinate system are respectively established through two times of dynamic coordinate rotation transformation, voxel division, chaff concentration statistics and attenuation factor calculation are independently executed on a bidirectional path, and accurate modeling of incidence and scattering path attenuation is realized in combination with a linear moment method basic scattering field solution. According to the method, the scattered field calculation precision of the non-uniform chaff cloud cluster in the complex electromagnetic environment is remarkably improved, meanwhile, the calculation complexity is reduced through voxel-level discretization, and the method can be applied to electronic countermeasure efficiency evaluation.

The invention discloses a design method and system of a high-temperature-resistant bearing wave-absorbing / wave-transmitting structure, and belongs to the technical field of calculation of electromagnetic structural strength crossing. Comprising the following steps: carrying out mechanical, thermal and electromagnetic performance characterization on a bearing core layer, and constructing an equivalent transmission line model of the bearing core layer; constructing an equivalent circuit model based on the equivalent transmission line model; importing the structure parameters into the equivalent circuit model, calculating the electrical performance, and researching the influence of the structure parameters on the electromagnetic performance to obtain key design parameters; establishing an optimization model by taking the key design parameters as optimization objects, the electromagnetic performance indexes as optimization targets and the force and thermal performance indexes of the structure as constraint conditions; key design parameters are initialized and imported into the optimization model, deformation and temperature performance indexes of the structure are calculated, and electromagnetic performance indexes under all temperature detection points are obtained; judging whether the thermal performance meets a constraint condition or not and whether the electromagnetic performance index meets an optimization target or not; and if not, modifying the key design parameters until the conditions are met.

The invention discloses a coated target electromagnetic scattering modeling method based on multi-core heterogeneous parallel computing, and belongs to the field of computational electromagnetism and high-performance computing. The method comprises the steps that an MPI + OpenMP + DSP three-level heterogeneous parallel framework is constructed, and optimal configuration of computing resources is achieved through one-to-one binding of an MPI process and a DSP acceleration card and one-to-one mapping of an OpenMP thread and a CPU core; dispersing a to-be-solved coating target by adopting a triangular grid, and defining an RWG primary function to expand surface current; substituting the surface current approximate expansion into an electric field integral equation based on impedance boundary conditions; carrying out impedance matrix calculation by adopting a process-based parallel filling strategy, and carrying out parallel solution on a matrix equation based on a ScaLAPACK library; and finally, calculating a far-region scattering field and a radar scattering sectional area according to the solved current coefficient. According to the method, through an innovative heterogeneous parallel architecture and a calculation task scheduling strategy, the problems of insufficient memory and low calculation efficiency during processing of an electrically large coated target in a traditional method are effectively solved, and an efficient and reliable technical means is provided for electromagnetic stealth design of a complex target.

The present application discloses a flux linkage identification method and device for a six-phase permanent magnet synchronous motor, relating to the technical field of computational electromagnetics. The method includes: constructing a motor model of a target motor in a synchronous speed coordinate system; constructing a model reference adaptive system of the target motor based on the motor model; performing inductance identification on the parameter to be identified based on the model reference adaptive system to obtain an accurate inductance value; determining an inductance error factor based on the accurate inductance value, and further obtaining a flux linkage error factor, and substituting the flux linkage error factor into a flux linkage identification model to obtain an accurate flux linkage value, realizing the flux linkage identification of the motor based on inductance identification transfer and the model reference adaptive system, effectively improving the parameter identification efficiency and accuracy, enhancing the robustness of the motor model to parameter changes, and improving the control performance of the permanent magnet synchronous motor.

The invention belongs to the crossing field of computational electromagnetism and artificial intelligence, and particularly provides a three-dimensional metal target electromagnetic simulation method based on a neural network and an equivalence principle. The method is used for solving the technical problems that in three-dimensional metal target electromagnetic scattering characteristic analysis, an equivalent principle algorithm is high in calculation complexity and large in memory consumption in operator matrix construction and scattering characteristic solving stages, and a traditional numerical integration method is low in efficiency. According to the method, the deep neural network is combined with multi-frequency triangular features and a space coupling perception activation mechanism, a prediction model of an equivalent principle operator matrix is constructed, high prediction precision is maintained, computing resources and time required for operator matrix construction and scattering characteristic solving are remarkably reduced, and the method is suitable for large-scale popularization and application. And an efficient, stable and extensible rapid solution scheme is provided for the electromagnetic scattering problem of the complex three-dimensional metal target.

The application discloses a three-dimensional space separation heterogeneous region scattering field simulation method combined with a neural network, and belongs to the technical field of computational electromagnetics and artificial intelligence. The method divides a solution domain into a background domain and an observation domain, and the background domain and the observation domain are separated from each other in space. Cross-domain numerical Green functions and in-domain numerical Green functions are extracted through virtual probe scanning. Physical characteristic vectors corresponding to field source relations are extracted, and training sample sets are constructed in combination with the numerical Green functions. Neural network models are trained by using the training sample sets to respectively fit the cross-domain and in-domain numerical Green functions. The trained neural network models are used to reconstruct numerical Green functions corresponding to to-be-solved foreground scatterers in real time, integral equation solutions of equivalent currents of the foreground scatterers are established, and three-dimensional heterogeneous region scattering fields are calculated in the observation domain. The method can realize fast and accurate calculation of numerical Green functions in a three-dimensional heterogeneous region scattering scene.

The invention relates to the cross technical field of computational electromagnetism and precision permanent magnet device design, in particular to a method for improving the nonuniformity of a magnetic disk space magnetic field, and aims to solve the problems that a method for improving the uniformity of a magnetic field in the prior art is low in optimization efficiency and difficult to break through the uniformity bottleneck caused by angular harmonic waves and a structural coupling effect. According to the scheme, the method comprises the steps that an accurate finite element model containing seven circles of fan-shaped magnets, a parallel magnetic disk, an eddy current disk and an iron support is established in Opera3D, and meshing and solving are carried out; in the first stage, the structural height of the iron stand is optimized through parameterized scanning, and the optimal height parameter is determined; and in the second stage, on the basis of the optimal iron stand height, the outermost two circles of the magnetic disk are set as NdFeB. The non-uniformity deviation (delta B) of the magnetic field of the target spherical area can be reduced by 13%, and the method is suitable for high-precision scenes such as nuclear magnetic resonance and synchrotron radiation light sources and has the advantages of being efficient in design and high in repeatability.

The application belongs to the technical field of computational electromagnetics, and discloses a numerical calculation method for simulating electromagnetic wave propagation characteristics in time-varying plasma, which is a three-dimensional Newton-HIE-FDTD method for efficiently simulating electromagnetic wave propagation in time-varying plasma based on a mixed explicit-implicit finite-difference time-domain method and Newton's second law of motion. The method inherits the advantages of the Newton-ADE-FDTD method, accurately represents the time-varying electron density through the time derivative, and can accurately capture the rich dynamic characteristics of the time-varying plasma. On this basis, the method also integrates the advantages of the HIE-FDTD method, so that the time step of the algorithm is only restricted by the grid size in two directions. Compared with other simulation methods based on time-varying plasma, the method is more efficient in calculation, and is particularly suitable for the analysis and simulation of electromagnetic wave propagation characteristics in time-varying plasma.

The invention relates to the field of computational electromagnetics, and provides a surface wave impedance rapid calculation method, device and equipment, a medium and a program product, the method comprises the following steps: step 1, inputting a dielectric constant of a medium material uniformly coated above an infinite metal surface; 2, calculating the surface wave impedance of the input dielectric constant by adopting an explicit solution expression of the surface wave impedance; and step 3, outputting the calculated wave impedance of the surface wave. The calculation precision of the surface wave impedance is higher than that of an existing approximation method.

The present invention belongs to the field of computational electromagnetics and relates to a method for interpolating and sweeping frequencies of transmission scattering parameters with adaptive sampling. Based on rational function fitting and according to the results of the previous step, the present invention uses the difference between the integrals of two fitting functions within a sub-interval, which better reflects the degree of change of the transmission scattering parameters within the interval and reduces the number of fittings; effectively improves the selection efficiency of sampling points, reduces the insertion of redundant sampling points, and improves the fitting accuracy for severely changing sections. In severely changing sections where the number of sampling points is insufficient, the difference between the two fittings is large, and the adaptive method of the present invention will naturally tend to insert new sampling points into this section; conversely, in gently changing sections, the difference between the two fittings before and after is small, and there is no need to insert new sampling points; thus improving the efficiency and accuracy of interpolating and sweeping frequencies. The present invention effectively solves the problems of poor adaptive sampling effect and complexity in existing methods.

The invention discloses an electromagnetic field prediction method based on a deep learning network model, belongs to the technical field of calculation of absorption boundaries in electromagnetics, constructs and trains an Autoformer model, applies a sequence decomposition mechanism and a sparse autocorrelation mechanism to a WCS-FDTD framework based on an electromagnetic model so as to update electromagnetic field data more quickly and accurately, and improves the accuracy of the electromagnetic field prediction. And inputting the obtained historical internal field domain and the electromagnetic field data at the corresponding absorption boundary into a trained neural network model for prediction, and calculating the electromagnetic field data of the internal field domain at the next moment based on the predicted electromagnetic field data at the absorption boundary at the next moment and the electromagnetic field data of the internal field domain at the current moment. And combining the obtained predicted value of the electromagnetic field data at the absorption boundary corresponding to the next moment with the electromagnetic field data of the internal field domain at the next moment, generating input data at the next moment, carrying out iterative prediction, and finally obtaining the electromagnetic field data of the computational domain at each moment.

The invention belongs to the technical field of block iteration matrix solving in computational electromagnetics, and particularly relates to a block iteration matrix solving method based on matrix dimensionality reduction and a preprocessing technology. According to the method, two technologies of matrix dimension reduction and matrix preprocessing are applied, a right-end item matrix is processed by adopting a matrix dimension reduction mode, and through low-rank approximation, the solution precision can be effectively ensured, the calculation overhead is remarkably reduced, and a multi-right-end linear system can be conveniently processed; moreover, the BCOCG algorithm is applied to a preprocessing technology, so that the performance of the matrix is further improved, and the convergence and stability of the algorithm are improved. According to a system equation generated when multiple physical excitation sources or multiple scanning angles exist in the electromagnetic problem, memory expenditure is shortened from the perspective of a short recursion Krylov subspace iteration method, multiple columns of right-end term vectors can be calculated at the same time through a block iteration BCOCG algorithm, and the problem of matrix morbidity caused by correlation between columns is solved.

The invention discloses a broadband uniaxial anisotropic complete matching layer absorption boundary method, and belongs to the field of computational electromagnetism. The method comprises the following steps: initializing electromagnetic field component coefficients (Ex, Ey, Hz), an excitation source and W-UPML parameters; constructing a matrix equation through an AH domain differential operator to solve a magnetic field component coefficient, and converting an AH domain field value into a time domain field value in combination with a domain inverse transformation formula; and high-efficiency calculation and boundary absorption of a broadband electromagnetic field are realized by utilizing a Hermite orthogonal basis function and an attenuation factor b. Compared with a traditional method, the time step size stability limitation is avoided through AH domain transformation, and the calculation efficiency is remarkably improved; the low-frequency wave absorption performance is enhanced by the introduced attenuation factor, and the reflection error is reduced by 20dB compared with that of the traditional UPML. The method is suitable for the fields of electromagnetic simulation, seismic wave simulation and the like, and effective technical support is provided for high-precision and broadband electromagnetic field analysis.

The invention discloses an unconditionally stable UPML electromagnetic simulation method based on Dosin FDTD, and belongs to the field of computational electromagnetics and numerical analysis, and the method comprises the steps: introducing an anisotropic damping term under a Dosin finite difference time domain (MS-FDTD) framework, embedding the UPML absorption characteristic into a Dosin Maxwell system, and carrying out the UPML electromagnetic simulation. And the consistency of the main domain and the absorption layer in geometric structure and energy transfer is realized. By adopting a Preissmann implicit discrete format and a local one-dimensional operator splitting strategy, a numerical value updating form irrelevant to a time step and a space step is obtained, so that unconditional stable solution is realized. According to the method, the electromagnetic field quantity is decoupled through the intermediate variable of the electric displacement vector and the magnetic induction intensity, and the solving efficiency and the numerical precision are improved. The provided Dosin UPML method can keep energy conservation and structural stability in complex media and long-time simulation, and the electromagnetic wave absorption precision and simulation reliability are remarkably improved.

This invention belongs to the field of numerical simulation in computational electromagnetics, specifically a non-uniform conformal FDTD method based on submesh technology. This invention introduces non-uniform conformal geometric corrections within the submesh to achieve high-precision fitting of curved surface boundaries, reducing step approximation errors and improving the accuracy of field calculations near the boundaries. At the master-submesh interface, a bidirectional coupling transfer is established that simultaneously satisfies line integral conservation and flux conservation, avoiding disruption of the conservation relationships required by the non-uniform conformal FDTD method. This effectively suppresses numerical reflections at the interface, reduces energy non-physical errors, and improves long-term stability. Furthermore, by using fine meshes only in the neighborhood of local microstructures and coordinating with time sub-cycle progression, global mesh refinement and the global time step being limited by the smallest unit are avoided, significantly reducing storage and computation time overhead while maintaining accuracy. This invention balances geometric modeling accuracy, interface stability, and computational efficiency, and has high engineering application value.

The invention relates to a mixed absorption boundary electromagnetic wave numerical simulation method, belongs to the technical field of computational electromagnetics, and solves the technical problem that a Mur absorption boundary condition cannot be directly used for a coordinate telescopic medium, and the simulation method comprises the following steps: establishing a mixed absorption boundary electromagnetic wave numerical calculation model; setting CPML (Convolution Complete Matching Layer) parameters; calculating electromagnetic field distribution in the model; calculating a parameter I A, a parameter II B and a parameter III C required by the mixed boundary; calculating the electric field of the boundary of the outermost layer of the CPML; calculating a magnetic field at the next moment; and a step of iteration. The Mur absorption boundary is seamlessly combined with the CPML, and an electromagnetic field in a coordinate telescoping medium is solved by using time domain finite difference FDTD. The method is used for mixed absorption boundary electromagnetic wave numerical simulation in the coordinate telescopic medium.

The present application belongs to the technical field of block iterative matrix solution in computational electromagnetics, in particular to a block iterative matrix solution method based on matrix dimension reduction and preprocessing technology. The present application uses two technologies of matrix dimension reduction and matrix preprocessing: the right end item matrix is processed by using the matrix dimension reduction, through low rank approximation, which can effectively ensure the accuracy of the solution, significantly reduce the calculation cost, and facilitate the processing of multiple right end linear systems; and the BCOCG algorithm is applied to the preprocessing technology to further improve the matrix behavior and improve the convergence and stability of the algorithm. For the system equation generated by multiple physical excitation sources or multiple scanning angles in electromagnetic problems, the memory cost is shortened from the perspective of short recursive Krylov subspace iteration method, and the block iteration BCOCG algorithm can calculate multiple column right end item vectors at the same time, avoiding the matrix ill-conditioning problem caused by the correlation between columns and columns.