41 results about "Finite-difference time-domain method" patented technology

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 application provides a transient ground potential rise induced coupling determination method through cable shielding layer, and belongs to the technical field of cable shielding layer.The application establishes a distributed parameter partial element equivalent circuit model of a grounding network, solves transient ground potential distribution by using a finite difference time domain method, obtains equivalent transfer impedance frequency characteristics of a braided shielding layer by using a four-port S parameter vector network analysis method and a variational asymptotic homogenization multi-scale algorithm, solves core line common mode voltage and differential mode voltage by establishing a multi-conductor transmission line equation group, constructs an artificial intelligence prediction model based on a topology perception graph neural network, solves a large-scale sparse admittance matrix by using a multilevel fast multipole method to obtain full network induced voltage distribution, and establishes a quantitative relationship model between ground potential rise and core line induced voltage, so that the technical problems that cable core line induced voltage caused by transient ground potential rise is difficult to accurately predict and effectively suppress are solved.

The invention provides a plasma flow field parameter inversion method and system based on a 1D-CNN, and relates to the technical field of electromagnetic inverse scattering, and the method comprises the steps: building a plasma flow field model of a hypersonic speed target, calculating the information of a scattered field under the action of electromagnetic waves through a current density convolution-time domain finite difference method, and calculating the information of the scattered field under the action of the electromagnetic waves; solving a far field electric field component according to a near-far field extrapolation formula; sample data of multiple groups of far-zone time domain electric field components are obtained by changing inversion parameters, and a data set is constructed after the sample data are preprocessed; training a pre-built 1D-CNN model according to the data set until a preset training condition is met, and obtaining a parameter inversion model; and performing inversion according to the parameter inversion model to obtain plasma flow field parameters of the hypersonic speed target.

This invention discloses a method for predicting electromagnetic radiation from cables. First, a one-dimensional conductor model of the target cable is established in a three-dimensional simulation space, and its axial direction is divided into multiple segments. Current data for each segment is solved using the finite-difference time-domain method. Then, the three-dimensional simulation space is discretized using a Yee grid, and virtual observation points corresponding to actual observation points are determined. Standard electric field and standard magnetic field update equations are established on the Yee grid based on Maxwell's equations. A field-circuit co-simulation process is established, and in each time step iteration, the current data is equivalently coupled to the standard electric field and standard magnetic field update equations to achieve synchronous updates of transmission line current and spatial electromagnetic field. Finally, this process is executed until a preset termination condition is met, and the time-domain electromagnetic field response at the virtual observation points is extracted as the time-domain electromagnetic radiation prediction result. The advantages are high prediction accuracy and high prediction efficiency.

This invention discloses a method for detecting aflatoxin B1 based on a BIC terahertz metamaterial sensor. Based on electromagnetic theory and the time-domain finite-difference method, a continuous-domain bound-state (BIC) terahertz metamaterial sensor was designed. Using this BIC terahertz metamaterial sensor, terahertz enhanced spectra of aflatoxin B1 solutions of different concentrations were acquired. Analysis revealed that the amplitude of the transmission peak of aflatoxin B1 around 1.5 THz gradually decreased with increasing aflatoxin B1 solution concentration, while the frequency of the transmission peak gradually shifted to lower frequencies with increasing aflatoxin B1 solution concentration. PLS and LS-SVM quantitative models for aflatoxin were established. Analysis showed that the LS-SVM model performed better, and the detection limit for aflatoxin B1 was found to be 2.04 × 10⁻⁶. ‑7 μg / ml. This invention enables rapid, non-destructive, and highly sensitive detection of aflatoxin B1 in grains.

The invention discloses a preparation method of a near-infrared achromatic super lens based on a nanoimprint technology, and relates to the technical field of micro-nano manufacturing. The preparation method comprises the following steps: 1, designing a target super lens by using finite element simulation software based on a time domain finite difference method, a particle swarm algorithm and a scalar diffraction algorithm; 2, preparing a silicon female die corresponding to the super lens according to the designed target super lens parameters by using electron beam lithography and deep reactive ion etching; 3, carrying out surface fluorination treatment on the silicon female die; and 4, preparing a corresponding sub-mold according to the silicon female mold, carrying out imprinting by using a nano imprinting technology and the sub-mold, and finally, etching the a-Si layer on the surface of the substrate through deep reactive ion etching by taking imprinting glue as an etching mask to obtain the super lens. According to the invention, the nanoimprint technology can be used for manufacturing the near-infrared achromatic superlens, so that the near-infrared achromatic superlens gets rid of dependence on a high-precision photoetching machine.

The invention belongs to the technical field of lunar and deep space exploration, and relates to a simulation method, device, equipment, medium and program for a lunar seismic exploration task, which are used for realizing seismic source-propagation-observation full-link closed-loop simulation, and the method comprises the following steps: establishing a two-dimensional non-uniform lunar soil wave field propagation numerical model based on a finite difference time domain method; the whole process that lunar seismic waves are generated from a seismic source and spread in different circle layer media in the moon is simulated, and a theoretical acceleration signal at a preset station position is output; establishing an MEMS lunar seismograph system model for simulating the working characteristics of the MEMS lunar seismograph in the extreme environment of the lunar surface; the theoretical acceleration signal is input into an MEMS lunar seismometer system model for simulation calculation, and a simulation output voltage signal containing noise and instrument response characteristics is obtained; and calculating the cross correlation coefficient of the signals under different lunar earthquake station layout schemes based on the simulation output voltage signal, and further quantitatively evaluating the detection efficiency of the different station layout schemes.

The application discloses a complex ground surface ground penetrating radar amplitude-preserving imaging method based on an orthogonal skin mesh, and comprises the following steps: 1) constructing a ground surface model according to ground fluctuation data; constructing an electromagnetic wave underground propagation speed model according to geological exploration data; 2) generating an orthogonal skin mesh by using a Ryskin-Leal method; 3) discretizing the orthogonal skin mesh by using a finite difference time domain method; 4) performing reverse extrapolation on the received ground penetrating radar signal, and reversely propagating the electromagnetic wave field to an initial time; in the reverse extrapolation process, the spatial position of a reflection interface is calculated through a cross-correlation imaging condition, and an underground imaging profile is obtained; 5) calculating a target function gradient, and solving an inversion problem by using a conjugate gradient method to iteratively update the electromagnetic wave underground propagation speed model; and 6) generating an imaging profile of an underground structure according to the iteratively converged electromagnetic wave underground propagation speed model. The application significantly improves the precision, stability and anti-noise performance of underground imaging.

The application discloses a mine electromagnetic wave analysis method, which solves the problem that the ray tracing method cannot calculate the mine electromagnetic wave characteristics when the transmitting and receiving antennas are close to the roadway sides at the same time, and solves the problem that the finite difference time domain method cannot calculate the characteristics of the mine electromagnetic field above GHz in a large-section long roadway. Different calculation regions are divided in the simulation roadway, a data transmission interface between the ray tracing method and the finite difference time domain method is designed according to the law of conservation of energy, and the accuracy and efficiency of the mine electromagnetic wave propagation characteristics when the transmitting and receiving antennas are close to the roadway sides at the same time are improved. Guidance is provided for mine wireless communication base station planning and layout and installation position of the transmitting and receiving antennas, and the wireless transmission distance of the wireless communication system is improved under the premise that the transmitting power and the receiving sensitivity are unchanged.

The application discloses a near-bit formation detection method based on while-drilling acoustic forward-looking, and belongs to the technical field of formation detection, and comprises the following steps: S1, acquiring a while-drilling acoustic signal generated by rock breaking of a drill bit in a drilling process, performing time-frequency analysis on the while-drilling acoustic signal, and determining a while-drilling acoustic source feature in combination with theoretical analysis; S2, performing independent component analysis on the while-drilling acoustic signal, and combining with ensemble empirical mode decomposition (CEEMD) and wavelet transform to identify effective while-drilling acoustic signals from complex background noise; and S3, establishing a geophysical model of a bad geological body containing a lithology boundary surface, a fault fracture zone, a karst cave and a water-rich formation, and obtaining a propagation rule and characteristic parameters of the while-drilling acoustic signal in the model by using a time-domain finite difference numerical simulation method. Through multi-algorithm collaborative processing, the while-drilling acoustic signal with high signal-to-noise ratio and obvious characteristics is extracted from the complex background noise, and the reliability and accuracy of signal processing are improved.

The invention discloses a three-phase circuit lightning stroke identification method and device, relates to the technical field of power system positioning, and aims to accurately judge the lightning stroke type and the lightning stroke position, conveniently identify the line accident reason and quickly position the line fault caused by the lightning stroke. The method comprises the following steps: fitting a relationship between a lightning peak value on a line and a transmission distance according to a finite difference time domain method; a transient magnetic field generated by lightning stroke is induced through a PCB differential ring, and induced electromotive forces generated at the two ends of the transient magnetic field are sampled; when the induced electromotive force is greater than a set threshold value, recording current waveform data of a three-phase line; and judging an output lightning stroke type according to the current waveform data and whether the polarities of the three-phase current head waves are the same or not, and judging a lightning stroke position in combination with a relationship between a lightning peak value and a transmission distance. The lightning stroke type and the lightning stroke position can be accurately judged, and line accident reasons can be conveniently identified and line faults caused by lightning stroke can be rapidly positioned.

The invention belongs to the technical field of electronic materials, and particularly relates to an ultrahigh-temperature-resistant high-entropy ceramic carbide composite wave-absorbing material and a preparation method thereof. The composite wave-absorbing material has a three-layer structure, namely a boron nitride high-entropy ceramic layer, a low-carbon-content high-entropy ceramic carbide layer and a high-carbon-content high-entropy ceramic carbide layer. And screening the thickness and the wave-absorbing efficiency of the composite wave-absorbing material by combining an impedance matching principle and a time domain finite difference method. And putting the screened composite wave-absorbing material with the three-layer structure into a high-temperature furnace, heating to be more than or equal to 1600 DEG C in an air atmosphere, and cooling to obtain the ultrahigh-temperature-resistant high-entropy ceramic carbide composite wave-absorbing material. The average wave-absorbing efficiency of the material in the frequency band of 2-18 GHz is smaller than or equal to-8 dB, and material support is provided for electromagnetic protection requirements in the aerospace extreme high-temperature aerobic environment.

The invention discloses a field line coupling time domain modeling method, device and equipment for a bent shielded cable and a medium, and relates to the field of electromagnetic compatibility and electromagnetic field numerical calculation, and the method comprises the steps: decomposing the bent shielded cable into an inner transmission line system and an outer transmission line system according to a cable shielding layer; the shielding cable structure is removed, and FDTD is used for simulating electromagnetic field distribution generated in the area where the shielding cable is located; dividing the shielding cable according to FDTD grids, and establishing an outer transmission line equation; performing differential discretization on the outer transmission line equation by adopting an FDTD method to obtain transient voltage and current response of each point along an outer layer; extracting current along a shielding outer layer, and establishing an inner transmission line equation; and performing differential discretization on the internal transmission line equation by adopting an FDTD method to obtain transient voltage and current responses of each point along the core wire and the terminating load. According to the method, a time domain finite difference method, a transmission line equation theory and an interpolation technology are combined, and synchronous calculation of space electromagnetic field radiation and transient response of a shielding cable bending shielding outer layer and a core wire is achieved.

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 discloses a hybrid implicit-explicit finite-difference time domain electromagnetic simulation method based on a body-centered square grid, belongs to the technical field of electromagnetic simulation, and aims at solving the problem that an existing method cannot give consideration to calculation efficiency and numerical isotropy at the same time. A mixed implicit-explicit time marching strategy is introduced, implicit updating is carried out on one space direction, explicit updating is carried out on other directions, and field quantity iteration is completed by solving a column-by-column three-diagonal linear system. According to the method, high precision and isotropy are kept, meanwhile, time step length limitation is remarkably relaxed, and the calculation efficiency of large-scale electromagnetic simulation is improved.

The invention discloses an electromagnetic pulse simulator field time domain waveform rapid estimation method based on a neural network. The technical problems that an existing calculation method is long in consumed time and large in calculation resource, and calculation results cannot be effectively utilized are solved. According to the method, a full-connection neural network model is constructed, and the mapping relation of different input voltage waveforms in the aspects of time and voltage amplitude is learned, so that an agent model is established; an electromagnetic field time-domain waveform of the electromagnetic pulse simulator at a diagnosis position under voltage excitation is calculated by adopting a time-domain finite difference method, and a calculation result is used as input of an agent model, so that rapid estimation of the field time-domain waveform of the electromagnetic pulse simulator under different voltage waveforms is realized. And the calculation speed of the field time domain waveform of the electromagnetic pulse simulator is greatly improved.

The invention discloses an electromagnetic simulation method of a microstrip line and electronic equipment. According to the microstrip line electromagnetic simulation method provided by the invention, the frequency domain polarization model is constructed by decomposing a complex dispersive medium model into superposition of a plurality of Debye relaxation terms, and the frequency domain response of each unit is converted into an independent time domain update relationship by using an auxiliary differential equation method; therefore, the high calculation overhead caused by the traditional time domain convolution operation is thoroughly avoided; furthermore, the updating relations and Maxwell equations are coupled to construct a joint solving model, discretization is carried out through a time domain finite difference method, an explicit field quantity updating operator is generated, and the electric field, the magnetic field and the polarization current component can be efficiently and synchronously updated in time domain iterative calculation. According to the method, high-precision modeling of the broadband dispersion characteristic of the PCB microstrip line is ensured, and meanwhile, the simulation calculation efficiency is remarkably improved.

The application discloses a kind of ground electrode in time domain finite difference simulation soil ionization effect characterization method, including obtaining the geometry of the ground electrode to be characterized, soil characteristics, environment and excitation source parameters, construct electrode in FTDT simulation model;The model is meshed, the absorption boundary condition and simulation total length are determined and initialized;Current density and electric field intensity around electrode are calculated, and soil resistivity is updated;The total resistance of electrode under the distribution of current soil resistivity is calculated, and the equivalent radius of electrode is calculated correspondingly;The correction factor of time domain finite difference is calculated using the equivalent radius of electrode, and then the vertical electric field component parameter of electrode is corrected;When the current simulation length reaches the simulation total length, the ground electrode top voltage, resistance, equivalent radius characterizing soil ionization effect are output;Method can accurately characterize soil ionization effect under coarse grid, greatly improve the calculation efficiency.

This invention discloses a method for fabricating a large-aperture convex grating with a high-depth-width blazed groove shape. The method includes: confirming the grating line density using dispersion theory; determining the grating's operating temperature under corresponding wavelength radiation backgrounds; determining the grating substrate material; performing a rough calculation of the blaze angle using scalar diffraction theory; performing a precise calculation of the blaze angle using the finite-difference time-domain method; determining the film material and thickness; designing the consistency of TE and TM diffraction efficiencies; analyzing the grating groove shape difference; confirming the grating groove shape parameters; electron beam exposure and development of the blazed groove shape; designing, exposing, and developing the exposure dose compensation; stitching together the large-aperture grating exposure areas; ion beam etching of the grating; and coating the grating with a metal reflective film. This fabrication method ensures the stability of the diffraction efficiency of the grating substrate and film materials under varying temperatures from room temperature to deep cryogenics, as well as the consistency of TE and TM diffraction efficiencies, thus guaranteeing the imaging quality and response stability of spectroscopic instruments.

Methods, systems, and non-transitory computer readable media for bypassing the bottleneck of graphics processing unit (GPU) parallelization during three-dimensional (3D) finite difference time domain (FDTD) simulation. Prior to running the 3D FDTD simulation, a central processing unit (CPU) creates a 3D array and a mirror image set of the 3D array and stores the mirror image set on the GPU. During 3D FDTD simulation run, the CPU sends an instruction to the GPU to update an array of interest of a mirror set on the GPU, while updating array elements of a 3D electromagnetic field component of each of the arrays of interest. And after 3D FDTD simulation is completed, the CPU copies the updated array of interest from the GPU to a corresponding mirror array pointed on the RAM by indicating the GPU, copies the updated array of interest from the GPU to the CPU, and writes the updated array of interest into an output file.

This invention discloses a non-uniform hexahedral mesh generation method based on the regional density distribution law, belonging to the field of software development technology, including the following steps: S1: For the solution target, triangular facets are used for mesh generation, and regions are divided according to different regional densities to obtain the maximum and minimum mesh size corresponding to each region; S2: After executing step S1, n regions are obtained. In the Cartesian coordinate system, mesh lines are generated sequentially along the three dimensions of x, y, and z. Optimization starts from the region with the smallest maximum mesh size. The advantages of this invention are: through modules such as the identification and calculation of the density distribution of triangular facet regions, the generation of three-dimensional mesh lines in the Cartesian coordinate system, the reconstruction of ray-intersecting hexahedral mesh, and mesh output, the conversion from triangular facet mesh to non-uniform hexahedral mesh is realized, solving the mesh generation and geometric target modeling problems of the finite-difference time-domain (FDTD) method.

The invention discloses an embedded framework-oriented lightning fire partial differential equation solving method. The method comprises the following steps of: inputting a form, an initial condition, a boundary condition, differential grid division and a neural network structure of a lightning fire equation; using a finite difference method and a time domain finite difference method to carry out discrete solution on the lightning equation based on a preset grid to obtain a discrete numerical solution on a definition domain; solving the lightning equation by adopting a physical information neural network, and performing error correction by embedding the discrete numerical solution to obtain a complete numerical solution in a time-space domain; and outputting a visual result of the complete numerical solution. According to the invention, high-precision approximation and solution of the whole partial differential equation system are realized.

The invention discloses a method for determining a coating area of a conical structure of a wave-absorbing material based on multi-objective optimization. The method comprises the following steps: acquiring geometric parameters and material attribute parameters of the conical structure through a data acquisition terminal; the modeling terminal establishes an axial layering three-dimensional model of a conical structure; the grid division terminal performs curvature self-adaptive unstructured tetrahedral grid division on each layer; dynamically optimizing the grid units by calculating the curvature of the vertex of the grid; the data analysis terminal adopts a time domain finite difference method to calculate electromagnetic response data of each grid unit in the target frequency band; the optimization decision terminal uses a non-dominated sorting genetic algorithm with an elitist strategy to iteratively search a solution space, optimizes electromagnetic performance and economical efficiency, and outputs a Pareto leading solution; the comprehensive evaluation value of each Pareto solution is calculated according to a preset weight coefficient, and a maximum value scheme is selected to generate a coating area coordinate mapping table, so that the triangular contradiction of'precision-efficiency-cost 'in coating of the cone structure wave-absorbing material is solved.

This application discloses an electromagnetic simulation method, apparatus, device, and readable storage medium for diodes. The method includes: constructing an electromagnetic simulation calculation region, which contains a simulated diode and one or more simulated objects; the electromagnetic simulation calculation region is the area of ​​interest when performing electromagnetic simulation on the simulated diode and each simulated object; based on the simulated diode as a lumped parameter element, performing finite-difference time-domain (FDTD) meshing on the electromagnetic simulation calculation region to obtain a set of orthogonal meshes; using the FDTD algorithm, iterating continuously on the electric field vector and magnetic field vector of each mesh and the simulated diode until a preset condition is reached, wherein the electric field of the electromagnetic simulation calculation region is composed of the electric field vectors of each mesh and the simulated diode, and the magnetic field of the electromagnetic simulation calculation region is composed of the magnetic field vectors of each mesh and the simulated diode. Therefore, this application can systematically provide a method for electromagnetic simulation of diodes.

This invention relates to a method for calculating the electromagnetic transient processes of high-voltage transmission lines based on the FD-WENO scheme. This method employs a high-order, inherently oscillatory-free differential FD-WENO scheme for numerical flux reconstruction in the spatial domain, and a total variation decay (TVD-RK) scheme for time-domain numerical integration, progressively solving for the time-varying curves of electrical quantities at various points in the transmission line. Compared with existing technologies, this invention effectively solves the numerical oscillation problem encountered when using traditional finite-difference time-domain (FDTD) methods to solve high-voltage transmission lines, and the spatial discretization accuracy of high-voltage transmission lines can reach arbitrary precision of order 2 or higher.

The invention discloses a surface wave absorber based on a silicon-based thermo-optical metasurface and a temperature sensor. The wave absorber comprises a plurality of resonator structures on the top layer, a middle dielectric layer and a bottom reflecting layer which are sequentially stacked from top to bottom. The top resonator structure is a cross-shaped notch circular ring coupling resonator based on amorphous silicon, the middle dielectric layer is silicon nitride, and the bottom reflecting layer is gold. Through simulation and optimization of a finite difference time domain method, the structure realizes perfect absorption of an ultra-narrow band with the absorption rate exceeding 99% and the full width at half maximum of only 7nm in a near-infrared band, and shows excellent light field regulation and control capability. In the aspect of temperature sensing, when the temperature is increased to 120 DEG C from 20 DEG C, the peak value of an absorption peak is shifted to 1539.27 nm from 1536.95 nm, the thermo-optical sensitivity is 80 pm / DEG C, the temperature measurement precision can reach + / -0.15 DEG C, and the huge potential of the super-structure surface based on the thermo-optical effect in the field of temperature sensing is shown.

The application discloses a block storage method suitable for a time domain finite difference method electric large-size grid, and relates to the technical field of software development.The method comprises the following steps: classifying a grid file, and acquiring grid data by using a hexahedral grid subdivision method; dividing a three-dimensional array of cell entities along three dimensions of x, y and z, and storing block files in an HDF5 format on a hard disk; mapping the three-dimensional array in a main grid file, mapping cell entity data dispersed in each block file into a continuous global array; and storing conformal data according to a conformal scheme.The block storage method of the electric large-size FDTD grid enables the FDTD solver to randomly read and write all the cells like accessing a single data set, while retaining the I / O advantage of parallel loading of sub-files, and meanwhile, the grid file is standardized and has expandability.

The application discloses a method for determining a coating area of a wave-absorbing material cone structure based on multi-target optimization, and comprises the following steps: collecting geometric parameters and material attribute parameters of the cone structure through a data acquisition terminal; establishing an axial layered three-dimensional model of the cone structure through a modeling terminal; performing curvature self-adaptive unstructured tetrahedral mesh division on each layer through a mesh division terminal; dynamically optimizing a grid unit by calculating the curvature of a grid vertex; calculating electromagnetic response data of each grid unit in a target frequency band through a data analysis terminal by using a finite difference time domain method; iteratively searching a solution space while optimizing electromagnetic performance and economy by using a non-dominated sorting genetic algorithm with an elitist strategy through an optimization decision terminal, and outputting a Pareto front solution; and calculating a comprehensive evaluation value of each Pareto solution according to a preset weight coefficient, selecting a highest value scheme to generate a coating area coordinate mapping table, and solving the triangular contradiction of "precision-efficiency-cost" in coating of the cone structure wave-absorbing material.

The invention discloses a method for designing an X-shaped crossed waveguide based on reverse design. The method comprises the following steps: an initial structure of the X-shaped crossed waveguide sequentially comprises two input waveguides (IN1 and IN2), a square region to be optimized and two output waveguides (O1 and O2) from left to right; the optimized area is divided into 16 * 16 pixel points; setting particle swarm algorithm parameters and initializing a particle swarm; processing information carried by the particles; simulating the particle swarm through a three-dimensional time domain finite difference method, and solving an FOM value; judging whether the FOM requirement is met or not, and outputting an optimal structure if the requirement is met; if the total number of iterations is reached, stopping optimization; if yes, outputting the current optimal structure; performing particle swarm optimization to obtain new particles; repeating simulation; and when the optimal structure is output, the simulation is finished, namely, the required X-shaped crossed waveguide is designed, and the method has the advantages of small size, high normalization power, large bandwidth and good performance of tolerance to manufacturing process deviation.

The application discloses a kind of based on conformal technology's core particle through silicon via interconnection structure signal integrity simulation excitation method, using time domain finite difference method analysis core particle interconnection structure application scene, to through silicon via structure model using conformal technology is analyzed;Overall calculation region is divided into excitation region and structure calculation region two sub-regions;Using total field / scattering field boundary technology, the result of excitation region is introduced as incident wave into structure calculation region, as excitation source to through silicon via structure is excited;The surface power integration at port in total field and scattering field is calculated, and by Fourier transform is converted into frequency domain calculation each port S parameter.The present application reduces the step approximation error of time domain finite difference method using conformal technology, by dividing the calculation region into two sub-regions solve the high-order mode crosstalk problem existing in typical coaxial structure such as through silicon via, and greatly improve the calculation precision of through silicon via structure signal integrity.