CN107992696B - Improved exponential time integral construction method in complex dispersion medium - Google Patents
Improved exponential time integral construction method in complex dispersion medium Download PDFInfo
- Publication number
- CN107992696B CN107992696B CN201711324894.XA CN201711324894A CN107992696B CN 107992696 B CN107992696 B CN 107992696B CN 201711324894 A CN201711324894 A CN 201711324894A CN 107992696 B CN107992696 B CN 107992696B
- Authority
- CN
- China
- Prior art keywords
- matrix
- vector
- grid cell
- unknown
- magnetic field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an improved exponential time integral construction method in a complex dispersion medium, which is applied to the field of three-dimensional electromagnetic analysis numerical values; the method comprises the steps of carrying out three-dimensional modeling on a multi-scale target to be analyzed and containing a complex dispersion medium, and establishing a corresponding geometric structure model; dividing the established geometric structure model by using a tetrahedral mesh; firstly, respectively integrating all local field components to respective corresponding global unknown quantities to construct a global semi-discrete format; then, the components of the electric field and the magnetic field are integrated into an unknown quantity to obtain an ordinary differential equation, namely an update equation of the electric field and the magnetic field; the updating of the electric field and the magnetic field only needs vector operation, and does not need to introduce a matrix index to realize updating like the traditional index time integration method, so that the dimension of the matrix index to be solved is obviously reduced, the calculation efficiency of the matrix index is effectively reduced, and the simulation efficiency of the time domain intermittent Galerkin method based on the index time integration method is improved.
Description
Technical Field
The invention belongs to the field of three-dimensional electromagnetic analysis, and particularly relates to a three-dimensional electromagnetic analysis numerical solving technology.
Background
With the development of stealth technology, internet of things electromagnetic compatibility, electromagnetic shielding and other technologies, various complex dispersion materials are more and more emphasized. Various stealth coating materials aiming at different application environments such as radar, infrared and even visible light in military stealth technology are widely adopted in the fields of anti-electromagnetic interference and electromagnetic shielding design of civil internet of things equipment, solar thin film batteries and the like. With the increasing complexity of electromagnetic environments, the demand for multi-spectral properties of these materials is increasing. Especially for a multi-scale target in a complex electromagnetic environment, the target itself has typical multi-scale characteristics in a geometric structure, and usually contains materials with dispersion or anisotropy properties, such as precious metals, transition metals or composite material coatings, and the coating materials not only inherit the multi-scale characteristics of the target itself in the geometric structure, but also have to face the multi-scale characteristics between different materials and between materials and metals (the precise simulation of the coating materials needs to perform grid encryption on material interfaces, thereby generating the so-called multi-scale characteristics on materials). These characteristics make a multi-scale target containing a complex dispersion medium to generate complex electromagnetic effects in a complex electromagnetic environment, and pose a serious challenge to the electromagnetic characteristics, stability and reliability of the target itself, and even to battlefield viability. Therefore, it is of great importance to accurately obtain the electromagnetic response characteristics of such a complex electromagnetic environment.
To truly and accurately simulate the electromagnetic response characteristics of a complex electromagnetic environment of such a multi-scale target, it is necessary to consider the dispersion and anisotropy characteristics of materials in simulation, and in addition, the materials have multi-spectral characteristics, and the accurate time-domain broadband electromagnetic simulation analysis of the materials is very challenging for the traditional computational electromagnetic method. The time domain discontinuous Galerkin method which is started in recent years is very flexible in terms of space dispersion and time dispersion, has high parallelism, and is particularly suitable for three-dimensional simulation analysis of such complex problems. However, the theory and the technology of the current time domain discontinuous Galerkin method have a great difference from the true realization of high-precision and high-performance electromagnetic simulation analysis of a multi-scale target containing complex dispersion media. For such large problems with typical multi-scale features, adaptive mesh encryption is one of the most important ways to guarantee computational accuracy and reduce computational overhead. However, the local mesh encryption makes the stability problem of the explicit time format more severe. If a global explicit time format is adopted, a normalized global time step will bring huge computational resource consumption, because the maximum time step must satisfy the stability condition of the minimum grid, resulting in a very small global time step and a large increase in the number of time iteration steps. This time step, determined by the minimum grid, is not necessary for large size grids. Therefore, in the local mesh encryption problem, the global time step will bring a lot of unnecessary computation overhead in the large-size mesh area. The fully implicit time format seems to solve this problem, and its unconditionally stable nature enables the time-domain discontinuous galois method to achieve the same accuracy of results with larger or even much larger time steps than the explicit time format. However, the cost is that a global linear equation set needs to be solved in each time iteration or a coefficient matrix of the global linear equation set needs to be inverted, and the huge computing resource consumption makes the full-implicit time-domain discontinuous Galerkin method difficult to be used for three-dimensional electromagnetic simulation analysis of actual large problems such as multi-scale complex electromagnetic environment problems. Therefore, there is an urgent need for a more efficient time format for time domain intermittent galileo analysis research of electromagnetic response characteristics of a complex electromagnetic environment for multi-scale targets.
In simulation analysis of a multi-scale target containing a complex dispersion medium, a traditional exponential time integration method starts from a local semi-discrete format derived by a time-domain discontinuous Galerkin method, and constructs a global semi-discrete format in the following form by respectively integrating all local field components into respective corresponding global unknown quantities
Wherein, the global diagonal block matrix is formed by the quality matrixes of grid cells; diagonal blocks of the global block matrix K correspond to the sum of the rigidity matrix and the self-acting flux matrix of each grid, and non-diagonal blocks correspond to the flux matrix on the interface;andrespectively, electric field, magnetic field, and auxiliary polarization current (introduced by the dispersion model) unknowns;∞,ωdand gammadIs a dispersion model parameter, is a constant. Then the field is dividedAndintegrated into an unknown vector u, the global semi-discrete format can be converted into the following ordinary differential equation form
Wherein
In the matrix index time integration method, it is necessary to solve a matrix index with respect to the coefficient matrix C. The calculation cost of the matrix index is almost exponentially increased along with the matrix dimension, and particularly for a multi-scale problem containing complex dispersion media, the matrix index dimension is huge, the solution cost is very large, and the performance improvement obtained by removing the stability limit of a fine grid in an index time integration format is obviously weakened.
Disclosure of Invention
In order to solve the problem that a traditional exponential time integration method is difficult to perform efficient and high-precision three-dimensional electromagnetic simulation analysis on a multi-scale target containing a complex dispersion medium, the invention provides an improved exponential time integration construction method in the complex dispersion medium.
The technical scheme adopted by the invention is as follows: a method of improved exponential time integral construction in a complex dispersive medium, comprising:
s1, carrying out three-dimensional modeling on a multi-scale target to be analyzed and containing a complex dispersion medium, and establishing a corresponding geometric structure model;
s2, subdividing the established geometric structure model by using tetrahedral meshes to obtain a plurality of tetrahedral meshes;
s3, constructing edge value problems, and deducing discontinuous Galerkin weak forms of local edge value problems in each tetrahedral mesh according to the Galerkin process;
s4, discretizing the discontinuous Galerkin weak form by adopting a basis function to obtain a discontinuous Galerkin semi-discrete format;
s5, integrating all field components into an unknown vector, and constructing a global semi-discrete format;
s6, integrating the electric field component and the magnetic field component into an unknown vector, and constructing an ordinary differential equation according to the global semi-discrete format obtained in the step S5;
s7, dividing the tetrahedral meshes obtained in the step S2 into fine meshes and coarse meshes according to the sizes of the meshes, and separating the unknown quantity of the ordinary differential equation according to the fine meshes and the coarse meshes;
s8, introducing a new exponential form unknown vector related by a fine grid to replace the unknown vector in the step S7, and obtaining a local unconditional stable exponential time integral format.
Further, the step S3 is to construct an edge value problem, specifically: according to a time domain Maxwell equation, by introducing a polarization current vectorObtaining a mixed Maxwell-Drude equation; describing metal boundary conditions by adopting an ideal electric wall; intercepting the calculation area by using a Silver-Muller absorption boundary condition;
wherein the content of the first and second substances,is an auxiliary potential shift vector, and t is a time variable.
Furthermore, the discontinuous galaojinweak form of the local edge value problem in each tetrahedral mesh specifically is:
wherein t is a time variable,∞is the relative dielectric constant at infinite frequency, omegadIs the plasma frequency, gammadAs the frequency of the impact is the frequency of the impact,an outer normal unit vector on the given plane is given for the ith grid cell,is the magnetic field vector inside the ith grid cell,is the electric field vector inside the ith grid cell,is the polarization current vector inside the ith grid cell,an outer normal unit vector on the given plane for the jth grid cell,is the electric field inside the jth grid cell, FhFor the set of all the faces, the number of faces,for the basis function vector corresponding to the ith grid cell,0for the vacuum dielectric constant, μ is the permeability of the medium in the calculated area.
Further, in step S4, the discontinuous galaojin semi-discrete format specifically includes:
wherein the content of the first and second substances,the unknown electric field vector inside the ith grid cell,the unknown magnetic field vector inside the ith grid cell,the unknown polarization current vector inside the ith grid cell,∞is the relative dielectric constant at infinite frequency, omegadIs the plasma frequency, gammadAs collision frequency, μ is the permeability of the medium in the calculation region,0is the dielectric constant in vacuum, t is the time variable, ViRepresenting a tetrahedron KiThe set of numbers for the mid-interface surface,a matrix of electric field-dependent qualities is represented,representing the magnetic field-dependent mass matrix, SiiPresentation interfaceijSelf-acting flux matrix of (S)ijRepresenting an interaction flux matrix of an ith grid cell and a jth grid cell, whereinijRepresenting the interface of the ith and jth grid cells.
Further, the ordinary differential equation expression in step S5 is:wherein C is a coefficient matrix and t is a time variable.
The invention has the beneficial effects that: the improved index time integral construction method in the complex dispersion medium utilizes the time domain intermittent Galerkin method based on the index time integral method to improve the index time integral construction method for the sparse characteristic of the matrix index generated during the fine simulation analysis of the multi-scale problem containing the complex dispersion medium, obviously reduces the dimension of the matrix index to be solved, effectively reduces the calculation efficiency of the matrix index and further improves the simulation efficiency of the time domain intermittent Galerkin method based on the index time integral method.
Drawings
FIG. 1 is a flow chart of the scheme of the invention.
Detailed Description
In order to facilitate the understanding of the technical contents of the present invention by those skilled in the art, the present invention will be further explained with reference to the accompanying drawings.
As shown in fig. 1, a scheme flow chart of the present invention is provided, and the technical scheme of the present invention is as follows: a method of improved exponential time integral construction in a complex dispersive medium, comprising:
and S1, carrying out three-dimensional modeling on the multi-scale target to be analyzed, which contains the complex dispersion medium, and establishing a corresponding geometric structure model. Multi-scale targets containing complex dispersive media to be analyzed, such as military vehicles, fighters and drones, are selected. And establishing a three-dimensional model according to the geometric characteristics and the material attributes of the target, and simultaneously adding a truncation absorption boundary condition to form a geometric structure model of the calculation region.
S2, subdividing the established geometric structure model by using tetrahedral meshes to obtain a plurality of tetrahedral meshes; specifically, the method comprises the following steps: a known procedure in the computational electromagnetics method is used for the geometric structure model created in the tetrahedral mesh partitioning step S1, and therefore this step will not be described in detail. The subdivided computing region is divided into a number of tetrahedral meshes, so that the continuous geometry space is converted into a discrete mesh space.
S3, constructing edge value problems, and deducing discontinuous Galerkin weak forms of local edge value problems in each tetrahedral mesh according to the Galerkin process;
the construction boundary value problem specifically includes: first, the following time domain Maxwell equation is used:
for complex dispersive media, embodiments of the present invention are described using the Drude dispersion model. In the dispersion model, the complex electric displacement can be expressed as
Wherein
Here, the∞Is the relative dielectric constant at infinite frequency, omegadIs the plasma frequency, gammadFor collision frequency, all three are Drude model parameters, which are directly related to material properties. Fourier transform is carried out on the formula (6) and the Fourier transform is brought into a Maxwell equation of the time domain to obtain the equation
The Fourier transform is obtained from the formula (7)
wherein the content of the first and second substances,is a vector of the electric field strength,is the vector of the magnetic field strength, mu is the magnetic permeability of the medium in the calculation area,an auxiliary potential shift vector is introduced for the Drude dispersion model, operator ^ represents rotation operation, and t is a time variable;
the embodiment of the invention adopts an ideal electric wall to describe the metal boundary condition:
truncating the calculation region by using a Silver-Muller absorption boundary condition:
wherein the content of the first and second substances,andthe incident excitation electric and magnetic fields, respectively, and the impedance Z and admittance Y, respectively, can be expressed as:
thus, equations (10), (11) and (12) together constitute the marginal problem of the present method.
Is provided withIs a unit of tetrahedronKiAdjacent grid cells, i represents the ith grid cell, then K can be deduced according to the Galerkin processiInterrupted galaojin weak form of medium local boundary problem:
wherein the content of the first and second substances,∞is the relative dielectric constant at infinite frequency, omegadIs the plasma frequency, gammadAs the frequency of the impact is the frequency of the impact,an outer normal unit vector on the given plane is given for the ith grid cell,is the magnetic field vector inside the ith grid cell,is the electric field vector inside the ith grid cell,is the polarization current vector inside the ith grid cell,an outer normal unit vector on the given plane for the jth grid cell,is the electric field inside the jth grid cell, FhFor the set of all the faces, the number of faces,for the basis function vector corresponding to the ith grid cell,0to calculate the relative vacuum permittivity of the medium in the region, μ is the permeability of the medium in the calculated region. The Galerkin process is a well-known process, is widely applied to a finite element method and a time-domain intermittent Galerkin method, and is not used any moreAnd (6) describing in detail.
S4, discretizing the discontinuous Galerkin weak form by adopting a basis function to obtain a discontinuous Galerkin semi-discrete format; the method specifically comprises the following steps: the discontinuous Galerkin method supports multiple types of basis functions, here high-order stacked vector basis functions are taken as an example. Local electromagnetic field inside tetrahedral mesh cellsCan be made of phiilThe linear combination of the basis functions is expressed as
Wherein d isiFor the number of local unknowns, determined by the order and type of the basis function, eilAnd hilAre the coefficients of the basis functions. Then, the discontinuous galaojin weak form obtained in step S4 is discretized according to the above formula to obtain tetrahedral mesh KiSemi-discrete format of local discontinuity Galerkin
Wherein, ViRepresenting a tetrahedron KiThe set of numbers for the mid-interface surface,a matrix of electric field-dependent qualities is represented,representing the magnetic field-dependent mass matrix, SiiPresentation interfaceijSelf-acting flux matrix of (S)ijRepresenting an interaction flux matrix of the ith grid cell and the jth grid cell; in addition, the method can be used for producing a composite material
ijRepresenting the interface of the ith and jth grid cells.
S5, integrating all field components into an unknown vector, and constructing a global semi-discrete format; the method specifically comprises the following steps: local unknownsRespectively integrated into global unknownsAndin the above-mentioned local semi-discrete format, the above-mentioned local semi-discrete format can be converted into the following global semi-discrete format
Wherein each diagonal block corresponds to a local quality matrix(electric field dependent mass matrix) or(magnetic field dependent mass matrix); the matrix K is a block matrix with diagonal blocks ofThe non-diagonal block is
S6, dividing the field Integration into an unknown vectorObtaining the following ordinary differential equation form, namely the update equation of the electric field and the magnetic field
Wherein C is a coefficient matrix;
here, the polarization current-dependent auxiliary equation introduced by the Drude dispersion model in equation (16)Only vector operation is needed, and updating is realized without introducing matrix indexes like a traditional exponential time integration method, so that the dimension of the matrix indexes to be required can be remarkably reduced.
S7, dividing the tetrahedral meshes obtained in the step S2 into fine meshes and coarse meshes according to the sizes of the meshes, and separating the unknown quantity of the ordinary differential equation according to the fine meshes and the coarse meshes; the method specifically comprises the following steps:
the meshes obtained in step S2 are classified into two types, fine meshes and coarse meshes, according to the mesh size. Then, based on this classification, the unknowns of the ordinary differential equations in step S6 can be separated into
Where the matrix P is a diagonal matrix with diagonal elements of 0 or 1, where 1 is used to mark the fine grid partial correlation unknowns and I is the identity matrix. The above ordinary differential equation (17) can be separated into
Wherein C isfCP and CcC (I-P) consists of fine and coarse grid related unknowns, respectively.
S8, introducing a new exponential form unknown vector related to the fine grid to replace the unknown vector in the step S7, removing the explicit dependence of the fine grid part on the ordinary differential equation time format, and obtaining a local unconditional stable exponential time integral format;
introducing new exponential form unknown vectors correlated by a fine grid
And the right end item can also be converted into
After simplification, the exponential time integral format with stable local unconditional condition is obtained
To this end, the time format-to-fine grid partial correlation matrix C is removed by introducing a new variable (22)fExplicit dependence of, CfNo separate time iteration is required. And the number of the first and second electrodes,andthe product of (a) is an identity matrix, so the matrix in equation (25)Is a coarse grid partial correlation matrix CcHave the same characteristic spectrum. In other words, the fine mesh portion no longer affects the stability of equation (25). Then after further time dispersion of equation (25), the resulting time format is unconditionally stable for the fine grid part.
By comparing the formula (3) and the formula (18), it can be known that the dimension of the coefficient matrix C in the improved exponential time format is reduced by 33% compared with the dimension of the coefficient matrix C in the conventional exponential time format derived matrix, and the corresponding generation matrix index is obtainedAlso reduced by 33%. The computational overhead of the matrix index increases almost exponentially with the matrix dimension, which is large especially for multi-scale problems involving complex dispersive media, becauseThe improved exponential time integration method can obviously improve the simulation analysis efficiency of the problems.
It will be appreciated by those of ordinary skill in the art that the embodiments described herein are intended to assist the reader in understanding the principles of the invention and are to be construed as being without limitation to such specifically recited embodiments and examples. Various modifications and alterations to this invention will become apparent to those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the claims of the present invention.
Claims (2)
1. A method of improved exponential time integral construction in a complex dispersive medium, comprising:
s1, carrying out three-dimensional modeling on a multi-scale target to be analyzed and containing a complex dispersion medium, and establishing a corresponding geometric structure model;
s2, subdividing the established geometric structure model by using tetrahedral meshes to obtain a plurality of tetrahedral meshes;
s3, constructing edge value problems, and deducing discontinuous Galerkin weak forms of local edge value problems in each tetrahedral mesh according to the Galerkin process; the discontinuous Galerkin weak form of the local edge value problem in each tetrahedral mesh specifically comprises:
wherein, KiRepresenting tetrahedral units, t is a time variable,∞is the relative dielectric constant at infinite frequency, omegadIs the plasma frequency, gammadAs the frequency of the impact is the frequency of the impact,an outer normal unit vector on the given plane is given for the ith grid cell,is the magnetic field vector inside the ith grid cell,the magnetic field vector inside the jth grid cell,is the electric field vector inside the ith grid cell,is the polarization current vector inside the ith grid cell,an outer normal unit vector on the given plane for the jth grid cell,is the electric field inside the jth grid cell, FhFor the set of all the faces, the number of faces,for the basis function vector corresponding to the ith grid cell,0is the vacuum dielectric constant, μ is the permeability of the medium in the calculation region;
s4, discretizing the discontinuous Galerkin weak form by adopting a basis function to obtain a discontinuous Galerkin semi-discrete format; step S4, where the discontinuous galaojin semi-discrete format specifically includes:
wherein the content of the first and second substances,the unknown electric field vector inside the ith grid cell,the unknown magnetic field vector inside the ith grid cell,the unknown magnetic field vector inside the jth grid cell,unknown polarization current vector, ω, inside the ith grid celldFor the frequency of the plasma to be,0is the dielectric constant in vacuum, t is the time variable, ViRepresenting a tetrahedron KiThe set of numbers for the mid-interface surface,a matrix of electric field-dependent qualities is represented,representing the magnetic field-dependent mass matrix, SiiPresentation interfaceijThe self-action matrix of (a) is,ijrepresenting the interface of the ith and jth grid cells, SijRepresenting an interaction matrix of the ith grid unit and the jth grid unit, and mu is the magnetic permeability of the medium in the calculation area;
s5, integrating all field components into an unknown vector, and constructing a global semi-discrete format; global semi-discrete format
Wherein each diagonal block corresponds to an electric field related quality matrixOr magnetic field dependent mass matrixThe matrix K is a block matrix with diagonal blocks ofThe non-diagonal block is
S6, integrating the electric field component and the magnetic field component into an unknown vectorConstructing an ordinary differential equation according to the global semi-discrete format obtained in the step S5; the ordinary differential equation expression described in step S6 is:wherein C is the corresponding coefficient matrix, t is the time variable,
s7, dividing the tetrahedral meshes obtained in the step S2 into fine meshes and coarse meshes according to the sizes of the meshes, and separating the unknown quantity of the ordinary differential equation according to the fine meshes and the coarse meshes; the ordinary differential equation is separated into
Wherein, CfCP and CcC (I-P) consists of fine and coarse grid related unknowns, respectively, P is a diagonal matrix with diagonal elements of 0 or 1;
s8, introducing a new exponential form unknown vector related to the fine grid to replace the unknown quantity in the step S7, and obtaining a local unconditional stable exponential time integral format; locally unconditionally stable exponential time integral format
2. The method according to claim 1, wherein the step S3 of constructing the boundary value problem is: according to a time domain Maxwell equation, by introducing a polarization current vectorObtaining a mixed Maxwell-Drude equation; describing metal boundary conditions by adopting an ideal electric wall; intercepting the calculation area by using a Silver-Muller absorption boundary condition;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711324894.XA CN107992696B (en) | 2017-12-13 | 2017-12-13 | Improved exponential time integral construction method in complex dispersion medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711324894.XA CN107992696B (en) | 2017-12-13 | 2017-12-13 | Improved exponential time integral construction method in complex dispersion medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107992696A CN107992696A (en) | 2018-05-04 |
CN107992696B true CN107992696B (en) | 2020-12-29 |
Family
ID=62037632
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711324894.XA Active CN107992696B (en) | 2017-12-13 | 2017-12-13 | Improved exponential time integral construction method in complex dispersion medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107992696B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109657288B (en) * | 2018-11-28 | 2022-07-26 | 电子科技大学 | Three-dimensional display and hiding time domain electromagnetic numerical method |
CN109684740B (en) * | 2018-12-27 | 2022-10-14 | 电子科技大学 | Electromagnetic multi-scale calculation method based on mixed grid and time step |
CN113536626B (en) * | 2021-06-25 | 2023-05-26 | 北京航空航天大学 | DGTD electromagnetic transient simulation method based on Cole-Cole model |
CN114076738B (en) * | 2021-11-22 | 2023-04-18 | 浙江大学 | Method for estimating soil salinity of residual straw farmland by using remote sensing construction index |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130041633A1 (en) * | 2011-08-11 | 2013-02-14 | Conocophillips Company | Heat front capture in thermal recovery simulations of hydrocarbon reservoirs |
CN103412989A (en) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | Parameterized reduced model based three-dimensional electromagnetic field simulation method of periodic structure |
CN103412988A (en) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | Three-dimensional electromagnetic field simulation method of periodic structure based on phase-shift reduced model |
CN104699909A (en) * | 2015-03-25 | 2015-06-10 | 厦门大学 | Variable step size multistep method time discrete algorithm keeping strong stability |
-
2017
- 2017-12-13 CN CN201711324894.XA patent/CN107992696B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130041633A1 (en) * | 2011-08-11 | 2013-02-14 | Conocophillips Company | Heat front capture in thermal recovery simulations of hydrocarbon reservoirs |
CN103412989A (en) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | Parameterized reduced model based three-dimensional electromagnetic field simulation method of periodic structure |
CN103412988A (en) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | Three-dimensional electromagnetic field simulation method of periodic structure based on phase-shift reduced model |
CN104699909A (en) * | 2015-03-25 | 2015-06-10 | 厦门大学 | Variable step size multistep method time discrete algorithm keeping strong stability |
Non-Patent Citations (2)
Title |
---|
A New Family of Exponential-Based High-Order DGTD Methods for Modeling 3-D Transient Multiscale Electromagnetic Problems;Hao Wang 等;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20171130;第65卷(第11期);第5960-5974页 * |
基于Rosenbrock型指数积分的一维间断Galerkin有限元方法;陈业飞 等;《应用数学和力学》;20130715;第34卷(第7期);第697-703页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107992696A (en) | 2018-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108052738B (en) | High-order local unconditionally stable time domain discontinuous Galerkin analysis method for dispersion medium | |
CN107992696B (en) | Improved exponential time integral construction method in complex dispersion medium | |
Deng et al. | High fidelity discontinuity-resolving reconstruction for compressible multiphase flows with moving interfaces | |
Musco et al. | Primordial black hole formation in the early universe: critical behaviour and self-similarity | |
Hildebrandt et al. | Novel formulation of nonlocal electrostatics | |
Pétri | The pulsar force-free magnetosphere linked to its striped wind: time-dependent pseudo-spectral simulations | |
Liu et al. | Mixed spectral-element method for 3-D Maxwell's eigenvalue problem | |
CN109684740B (en) | Electromagnetic multi-scale calculation method based on mixed grid and time step | |
CN111079278B (en) | Processing method for three-dimensional time domain hybridization discontinuous Galerkin method with additional electromagnetic source item | |
Andrade et al. | On the stability of gravity with Dirichlet walls | |
CN113158527B (en) | Method for calculating frequency domain electromagnetic field based on implicit FVFD | |
CN106294894B (en) | Finite element boundary integration method for rapidly analyzing electromagnetic scattering characteristics of non-uniform target | |
Vogman et al. | Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r, vr, vθ) phase space coordinates | |
CN107908903B (en) | Matrix exponential time decoupling method of exponential time integration method in complex dispersion medium | |
W. Hess et al. | A reduced basis method for microwave semiconductor devices with geometric variations | |
CN109657288B (en) | Three-dimensional display and hiding time domain electromagnetic numerical method | |
CN116776695A (en) | One-dimensional electromagnetic calculation method, system and equipment based on ultra-high-order finite element technology | |
Hendijani et al. | Constrained locally corrected Nyström method | |
CN107943755B (en) | Matrix index dimension reduction method of index time integration method in complex dispersion medium | |
Dolean et al. | Domain decomposition methods for electromagnetic wave propagation problems in heterogeneous media and complex domains | |
Beauregard et al. | Casimir effect in the presence of external fields | |
Fatalov | Brownian motion on with linear drift, reflected at zero: exact asymptotics for ergodic means | |
Lu et al. | Hybrid boundary integral-generalized (partition of unity) finite-element solvers for the scalar Helmholtz equation | |
Cheamsawat et al. | Free energy dependence on spatial geometry for (2+ 1)-dimensional QFTs | |
CN111931353A (en) | Scattered field solving method applied to simulated FSS structure |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |