CN-115879276-B - Electromagnetic characteristic analysis method, device and equipment of target object and medium

CN115879276BCN 115879276 BCN115879276 BCN 115879276BCN-115879276-B

Abstract

The application discloses a method and a device for analyzing electromagnetic characteristics of a target object, electronic equipment and a computer readable storage medium, wherein the method comprises the steps of performing simulation modeling according to boundary condition information of the target object to obtain a simulation model of the target object; the method comprises the steps of carrying out grid subdivision on a simulation model by adopting a preset structure, establishing a Maxwell equation set time-varying electromagnetic field based on the simulation model after grid subdivision, carrying out implicit solving on the Maxwell equation set time-varying electromagnetic field based on a physical time and virtual time double iteration mode, and outputting an electromagnetic characteristic analysis result of the target object when the Maxwell equation set time-varying electromagnetic field converges. The electromagnetic characteristic analysis method of the target object improves the electromagnetic characteristic analysis efficiency of the target object.

Inventors

XU YONG
Gao Tiesuo
CHEN JIANQIANG
SUN JUNFENG
JIANG TAO
DONG WEIZHONG
Ding Mingsong

Assignees

中国空气动力研究与发展中心计算空气动力研究所

Dates

Publication Date: 20260512
Application Date: 20221103

Claims (7)

1. A method of analyzing electromagnetic properties of a target object, comprising: Performing simulation modeling according to boundary condition information of a target object to obtain a simulation model of the target object; performing grid subdivision on the simulation model by adopting a preset structure, and establishing a Maxwell equation set time-varying electromagnetic field based on the simulation model after grid subdivision; Implicitly solving the time-varying electromagnetic field of the Maxwell equation set based on a physical time and virtual time double iteration mode; Outputting an electromagnetic characteristic analysis result of the target object when the time-varying electromagnetic field of the Maxwell equation set converges; The implicit solving of the time-varying electromagnetic field of the Maxwell equation set based on the physical time and the virtual time in a double iterative mode comprises the following steps: acquiring electromagnetic parameters and control parameters, and implicitly solving the time-varying electromagnetic field of the Maxwell equation set by utilizing the electromagnetic parameters and the control parameters based on a physical time and virtual time double iteration mode; the implicit solving of the maxwell equation set time-varying electromagnetic field based on a physical time and virtual time double iteration mode by utilizing the electromagnetic parameters and the control parameters comprises the following steps: The physical time step iteration loops until the physical time step iteration converges; In each physical time step iteration process, the virtual time step iteration loops until the virtual time step iteration converges; in each virtual time sub-iteration process, implicitly solving the time-varying electromagnetic field of the Maxwell equation set based on a physical time step and a virtual time step by utilizing the electromagnetic parameters and the control parameters so as to update the value of the conservation electromagnetic field of the next virtual time sub-iteration step number; wherein the maxwell equations change into electromagnetic fields: ; Wherein, the For the physical time step size of the device, For a virtual time step-size, Is the j, k, l grid cell The electromagnetic conservation variable of the virtual time step, Is the j, k, l grid cell The electromagnetic conservation variable of the virtual time step, Is the first The electromagnetic conservation variable of the physical time step, Is the first The electromagnetic conservation variable of the physical time step, Is the first The electromagnetic conservation variable of the physical time step, Is an implicit control parameter that is used to control the device, 、、 Respectively corresponding to a curve coordinate system 、、 A differential operator of the direction of the light, 、、 Respectively corresponding to a curve coordinate system 、、 A matrix of jacobian coefficients of the directional electromagnetic flux, Is the difference between adjacent physical time fringe field conservation variables, Is the first Physical time step spatial flux residual, RHS is the total spatial flux residual of the last physical time step plus physical time step and virtual time step correction.
2. The electromagnetic property analysis method according to claim 1, wherein the meshing of the simulation model with a predetermined structure includes: and meshing the simulation model by adopting a preset structure quadrilateral structure or hexahedral structure.
3. The electromagnetic property analysis method according to claim 1, wherein a grid density at a target position in the simulation model is inversely related to a distance between the target position and a wall surface, and inversely related to a distance between the target position and a geometric singular point.
4. The electromagnetic property analysis method according to claim 1, wherein outputting the electromagnetic property analysis result of the target object when the maxwell's equations time-varying electromagnetic field converges, comprises: And outputting any one or a combination of any two of time distribution, spatial distribution, surface induced current and radar scattering cross section spatial distribution data of the target object when the time-varying electromagnetic field of the Maxwell equation set is converged.
5. An electromagnetic property analysis apparatus for a target object, comprising: the modeling module is used for carrying out simulation modeling according to boundary condition information of the target object to obtain a simulation model of the target object; The establishing module is used for carrying out mesh division on the simulation model by adopting a preset structure, and establishing a Maxwell equation set time-varying electromagnetic field based on the simulation model after mesh division; the solving module is used for implicitly solving the time-varying electromagnetic field of the Maxwell equation set based on a physical time and virtual time double iteration mode; the analysis module is used for outputting an electromagnetic characteristic analysis result of the target object when the time-varying electromagnetic field of the Maxwell equation set is converged; the solving module is specifically used for acquiring electromagnetic parameters and control parameters, and implicitly solving the time-varying electromagnetic field of the Maxwell equation set based on a physical time and virtual time double iteration mode by utilizing the electromagnetic parameters and the control parameters; The solving module is specifically used for carrying out implicit solving on the time-varying electromagnetic field of the Maxwell equation set based on the physical time steps and the virtual time steps by utilizing the electromagnetic parameters and the control parameters in each virtual time sub-iteration process so as to update the value of the conservation electromagnetic field of the next virtual time sub-iteration step number; wherein the maxwell equations change into electromagnetic fields: ; Wherein, the For the physical time step size of the device, For a virtual time step-size, Is the j, k, l grid cell The electromagnetic conservation variable of the virtual time step, Is the j, k, l grid cell The electromagnetic conservation variable of the virtual time step, Is the first The electromagnetic conservation variable of the physical time step, Is the first The electromagnetic conservation variable of the physical time step, Is the first The electromagnetic conservation variable of the physical time step, Is an implicit control parameter that is used to control the device, 、、 Respectively corresponding to a curve coordinate system 、、 A differential operator of the direction of the light, 、、 Respectively corresponding to a curve coordinate system 、、 A matrix of jacobian coefficients of the directional electromagnetic flux, Is the difference between adjacent physical time fringe field conservation variables, Is the first Physical time step spatial flux residual, RHS is the total spatial flux residual of the last physical time step plus physical time step and virtual time step correction.
6. An electronic device, comprising: A memory for storing a computer program; a processor for implementing the steps of the electromagnetic property analysis method of the target object according to any one of claims 1 to 4 when executing the computer program.
7. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the steps of the electromagnetic property analysis method of a target object according to any one of claims 1 to 4.

Description

Electromagnetic characteristic analysis method, device and equipment of target object and medium Technical Field The present application relates to the field of time domain electromagnetics, and more particularly, to a method and apparatus for analyzing electromagnetic characteristics of a target object, and an electronic device and a computer readable storage medium. Background The time domain method can compatibly simulate complex phenomena such as scattering, multiple scattering, hole penetration, cavity excitation and the like, and can obtain system broadband information only through one calculation, so that the time domain method is widely applied to electromagnetic field problems, can accurately simulate time histories, and does not need to provide special processing methods for special components and special electromagnetic phenomena such as geometric edge diffraction unlike the traditional high-frequency progressive method. The time-varying electromagnetic field in the time domain satisfies the time domain maxwell's equations, which are directly solved as computer technology develops. The hyperbolic mathematical features, which are the same as the euler equation, facilitate the application of computational fluid dynamics (Computational Fluid Dynamics, CFD) techniques in electromagnetic field computation, with time-Domain finite difference methods (FINITE DIFFERENCE TIME Domain, FDTD) and time-Domain finite volume methods (Finite Volume Time Domain, FVTD) being the most notable. K.S.Yee published a precursor time domain finite difference algorithm in the 60 th century of 20 th, directly and differentially calculated a time-varying Max Wei Weifen equation set, successfully simulated the time domain response of the electromagnetic pulse and the action of an ideal conductor, and opened a new electromagnetic field time domain calculation method. In the Yee algorithm, first, a cartesian orthogonal grid is generated in a region of interest (a target and a certain space around the target), and electric field and magnetic field components are placed in a crossing manner at a value point of a grid space, so that the periphery of each electric field component is surrounded by the magnetic field component on each coordinate plane, and the periphery of each magnetic field component is surrounded by the electric field component, so that the electromagnetic field configuration meets the requirements of faraday induction law and ampere loop law, and the grid is commonly called a Yee grid. The traditional time domain finite difference method and the time domain finite volume method adopt a 2-order central difference or space-time coupling Lax-Wendroff format and a multi-step Runge-Kutta method, and the common point is an explicit format of time calculation. The explicit method represented by the Runge-Kutta method has the advantages of simple programming and easy realization of time high precision, and is a reliable time discrete method for calculating the time electromagnetic field. However, the time display method has a biggest defect that the time step is limited by stability, the whole calculation space must adopt a uniform minimum global calculation step, a patch encryption grid generated for simulating the severe change of geometric shape (for example, the severe change of electromagnetic field gradient caused by the geometric singularities of the edges of the front and rear edges of wings requires the encryption grid to be carefully simulated and the electromagnetic multiscale problem) can bring a small global time step, and a large grid unit needs more time steps to transmit information in the unit, so that a stable time-varying electromagnetic field needs longer calculation time, particularly when solving the electromagnetic scattering time domain problem of a high-frequency and electric large-size target, the small time step limited by stability brings remarkable increase of the calculation amount of the time-domain electromagnetic field, and a large amount of calculation resources are consumed. On the other hand, the implicit calculation method can relax the stability limit of the calculation step length, but the time precision is reduced and the encryption causes the increase of the dimension of the coefficient matrix to improve the matrix inversion operation difficulty. ADI (ALTERNATING DIRETION IMPLICIT, ADI) overcomes the limitation of the number of CFLs (Courant FRIEDRICHS LEWY, CFL) so that the selection of time steps is determined only by the calculation accuracy independent of Courant stability conditions, and thus the selection of time steps can be multiplied and the calculation time is also reduced. Namiki et al 1999 proposed an alternating direction implicit time domain difference (ADI-FDTD) method, which was subsequently demonstrated to be unconditionally stable. In 2009, cooke et al obtained a simplified single-step alternating direction implicit time domain finite difference (Leapfrog ADI-FD