ELECTROMAGNETIC DEVICE DESIGN SYSTEM FOR FAST FREQUENCY SWEEP BASED ON FINITE ELEMENT METHOD

Abstract

The disclosure provides an electromagnetic device design system including: one or more processors of a machine; and computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device. The operations include: initiating physical parameters of the EM device, wherein the EM device comprises multiple ports; performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM); applying single-size matrix Padé via Lanczos (MPVL) method in fast frequency sweep and performing EM simulation for the EM device under excitation at each port to obtain field solution of the EM device at frequencies in a frequency range; calculating S-parameters for the multiple ports of the EM device; and calculating derivatives of the S-parameters with respect to one of the physical parameters in the frequency range.

Description

CROSS-REFERENCE TO RELAYED APPLICATIONS

Pursuant to 35 U.S.C.§ 119 and the Paris Convention Treaty, this application claims foreign priority to Chinese Patent Application No. 202110992300.2 filed Aug. 27, 2021, the contents of which, including any intervening amendments thereto, are incorporated herein by reference. Inquiries from the public to applicants or assignees concerning this document or the related applications should be directed to: Matthias Scholl P C., Attn.: Dr. Matthias Scholl Esq., 245 First Street, 18th Floor, Cambridge, Mass. 02142.

BACKGROUND

The disclosure relates to the technical field of electromagnetic device design. In particular, the disclosure addresses systems and methods for performing an electromagnetic sensitivity analysis by using fast frequency sweep based on the finite element method.

Electromagnetic sensitivity analysis is important for electromagnetic (EM) device design. An EM sensitivity is defined as a derivative of an EM response, e.g., a full scattering matrix for ports of an EM device, with respect to a physical parameter of the EM device including geometrical parameters and material parameters. EM device design, such as design optimization, what if analysis and yield-driven design, benefits greatly from the EM sensitivity analysis.

As one of the most efficient numerical techniques, the finite-element method (FEM) is used to perform the EM device design routinely. However, using FEM to solve a large linear system is very time-consuming. Especially when the equations are frequency dependent, the solution process for a single frequency must be repeated once at each different frequency.

Recent efforts have focused on Model order reduction (MORe) techniques, such as the asymptotic waveform evaluation (AWE), the Galerkin AWE (GAWE), the well-conditioned AWE (WCAWE), the projection via Arnoldi, the Padé via Lanczos (PVL), the matrix Padé via Lanczos (MPVL), and the adaptive Lanczos-Pade sweep (ALPS). Using MORe techniques, the EM equation is merely solved at a single value of a parameter and solutions at other values of the parameter are estimated approximately with a small computational overhead, thus speeding up the whole simulating process.

However, model order reduction techniques cannot be directly used in the EM sensitivity analysis.

SUMMARY

The disclosure provides a system for performing an electromagnetic sensitivity analysis by using fast frequency sweep based on finite element method (FEM). Adjoint/self-adjoint formulas based on fast frequency sweep are derived, predicting the EM sensitivity information for the entire frequency range, and thus the efficiency of the EM sensitivity analysis is improved.

Conventional commercial software typically computes electromagnetic response derivatives at multiple frequencies and are less efficient. The disclosure incorporates a model order reduction algorithm and aims to speed up the electromagnetic sensitivity analysis. Increasing the number of frequencies to be solved results in a reduction in the execution time of the disclosure (for example, the disclosure gets a 30 times speedup when calculating the electromagnetic response derivatives at 100 frequencies).

Specifically, in one aspect, the disclosure provides an electromagnetic device design system comprising:

one or more processors of a machine; and

computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising:

initiating physical parameters of the EM device, wherein the EM device comprises multiple ports;

performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);

applying single-size matrix Padé via Lanczos (MPVL) method in fast frequency sweep and performing EM simulation for the EM device under excitation at each port to obtain field solutions of the EM device in a frequency range;

calculating S-parameters for the multiple ports of the EM device;

obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;

selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with a selected solution;

performing EM simulation for the EM device at a frequency of the selected solution using FEM; and

determining the simulated result satisfies a physical specification of the EM device;

wherein:

applying single-size MPVL method in fast frequency sweep comprises:

transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents a number of elements in a field vector and the single-size system matrix is of a linear combination of global finite-element system matrices;

generating a first linear system using the double-size system matrix;

representing solving vectors of the first linear system with the global finite-element system matrices using the block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; and

solving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining field solution in a frequency range as the following:

x _k ≈x _k ^q =∥r ₀ ∥V _k ^q(I_q−(s−s₀)T_k^q)⁻¹e₁

where

s represents a frequency;

s₀ represents a pre-solution frequency;

q is a reduced order in MPVL method;

x_k is a vector of field solution under excitation at port k;

x_k^q, represents a vector of q_th order reduced field solution under excitation at port k;

r₀ represents a solution vector of the first linear system at 0^th MPVL iteration;

V_k^q is represented as V_k^q=[v_m]_m+1^q=[v₁ v₂ . . . v_q], where v_m is an orthonormal basis vector of the Krylov subspace for model order reduction;

I_q is an identity matrix with a dimension of q×q;

T_k^q is a reduced order matrix; and

e₁ is represented as e₁=[1 0 . . . 0]^T with a dimension of 1×q.

In a second aspect, an electromagnetic device design system comprises:

one or more processors of a machine; and

computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising:

initiating physical parameters of the EM device wherein the EM device comprises multiple ports;

performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);

solving a linear system under excitation at port j to obtain a field solution under excitation at port j by fast frequency sweep;

solving an adjoint representation of the linear system under excitation at port k to obtain an adjoint field solution between port k and port j by fast frequency sweep;

calculating derivatives of a S-parameter of port j and port k with respect to O_i based on an adjoint sensitivity formula;

obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;

selecting a solution and updating EM device design to replace the one of the physical parameters of the EM device with the selected solution;

performing EM simulation for the EM device at a frequency of the selected solution using FEM; and

determining the simulated result satisfies a physical specification of the EM device;

wherein the adjoint sensitivity formula is written as

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}=-{\hat{x}}_{k,j}^T\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j,$

where

x_j represents a vector of the field solution under excitation at port j;

{tilde over (G)}_i represents a derivative of G with respect to ϕ_i;

{tilde over (C)}_i represents a derivative of C with respect to ϕ_i;

s is a frequency;

{circumflex over (x)}_k,j^T represents a transpose vector of the adjoint field solution between port j and port k;

S_k,j represents the S-parameter of port j and port k; and

ϕ_i represents an i_th one of the physical parameters.

In a third aspect, an electromagnetic device design system comprises:

one or more processors of a machine; and

computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising:

initiating physical parameters of the EM device wherein the EM device comprises multiple ports;

performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);

solving a linear system under excitation at port j to obtain a field solution under excitation at port j by fast frequency sweep;

solving the linear system under excitation at port k to obtain a field solution under excitation at port k by fast frequency sweep;

calculating derivatives of a S-parameter of port j and port k with respect to ϕ_i based on a self-adjoint sensitivity formula;

obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;

selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;

performing EM simulation for the EM device at a frequency of the selected solution using FEM; and

determining the simulated result satisfies a physical specification of the EM device;

wherein

the self-adjoint sensitivity formula is written as

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}=-{\kappa }_{k,j}{x}_k^T\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j,$

where

x_j represents a vector of the field solution under excitation at port j;

{tilde over (G)}_i represents a derivative of G with respect to ϕ_i;

{tilde over (C)}_i represents a derivative of C with respect to ϕ_i;

s is a frequency;

x_k^T represents a transpose vector of the field solution under excitation at port k;

S_k,j represents the S-parameter of port j and port k;

ϕ_i represents an i_th one of the physical parameters; and

κ_k,j represents a correlated coefficient of port j and port k.

In a fourth aspect, an electromagnetic device design system comprises:

one or more processors of a machine; and

computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising:

initiating physical parameters of the EM device wherein the EM device comprises multiple ports;

performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);

solving a linear system under excitation at port j by performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port j in a frequency range;

solving an adjoint representation of the linear system under excitation at port k by performing fast frequency sweep with single-size MPVL method to obtain an order reduced adjoint field solution between port k and port j in the frequency range;

calculating derivatives of a S-parameter of port j and port k with respect to O_i based on an adjoint sensitivity formula;

obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;

selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;

performing EM simulation for the EM device at a frequency of the selected solution using FEM; and

determining the simulated result satisfies a physical specification of the EM device;

wherein:

performing fast frequency sweep with single-size MPVL method comprises:

transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents the number of elements in a field vector and the single-size system matrix comprises a linear combination of global finite-element system matrices;

generating a first linear system using the double-size system matrix;

representing solving vectors of the first linear system with the global finite-element system matrices using block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; and

solving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining a field solution in a frequency range by the following:

x _k ≈x _k ^q =∥r ₀ ∥V _k ^q(I_q−(s−s₀)T_k^q)⁻¹e₁;

and the adjoint sensitivity formula is written as:

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\hat{x}}_{k,j}^{\mathrm{Tq}}\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j^q,$

where

s represents a frequency;

s₀ represents a pre-solution frequency;

q is a reduced order in MPVL method;

x_k is a vector of field solution under excitation at port k;

x_k^q, represents a vector of q_th order reduced field solution under excitation at port k;

r₀ represents a solution vector of the first linear system at 0^th MPVL iteration;

V_k^q is represented as V_k^q=[v_m]_m=1^q=[v₁ v₂ . . . v_q], where v_m is an orthonormal basis vector of the Krylov subspace for model order reduction;

I_q is an identity matrix with a dimension of q×q;

T_k^q is a reduced order matrix;

e₁ is represented as e₁=[1 0 . . . 0]^T with a dimension of 1×q;

{tilde over (G)}_i represents a derivative of G with respect to ϕ_i;

{tilde over (C)}_i represents a derivative of C with respect to ϕ_i;

{circumflex over (x)}_k,j^qT represents a transpose vector of qth order reduced adjoint field solution between port j and port k;

S_k,j represents the S-parameter of port j and port k; and

ϕ_i represents an i_th one of the physical parameters.

In a fifth aspect, an electromagnetic device design system comprises:

one or more processors of a machine; and

computer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising:

initiating physical parameters of the EM device wherein the EM device comprises multiple ports;

performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);

performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port j of the EM device in a frequency range;

performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port k of the EM device in a frequency range;

calculating derivatives of a S-parameter of port j and port k with respect to ϕ_i based on a self-adjoint sensitivity formula;

obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;

selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;

performing EM simulation for the EM device at a frequency of the selected solution using FEM; and

determining the simulated result satisfies a physical specification of the EM device;

wherein:

performing fast frequency sweep with single-size MPVL method comprises:

transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents the number of elements in a field vector and the single-size system matrix comprises a linear combination of global finite-element system matrices;

generating a first linear system using the double-size system matrix;

representing solving vectors of the first linear system with the global finite-element system matrices using block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; and

solving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining field solution of the EM device in a frequency range as the following:

x _k ≈x _k ^q =∥r ₀ ∥V _k ^q(I_q−(s−s₀)T_k^q)⁻¹e₁,

and the self-adjoint sensitivity formula is written as

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\kappa }_k{x}_k^{\mathrm{qT}}\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j^q,$

where

s represents a frequency;

s₀ represents a pre-solution frequency;

q is a reduced order in MPVL method;

x_k is a vector of the field solution under excitation at port k;

x_k^q represents q_th order reduced solution vector under excitation at port k;

r₀ represents a solution vector of the first linear system at 0^th MPVL iteration;

V_k^q is represented as V_k^q=[v_m]_m=1^q=[v₁ v₂ . . . v_q], where v_m is an orthonormal basis vector of the Krylov subspace for model order reduction;

I_q is an identity matrix with a dimension of q×q;

T_k^q is a reduced order matrix;

e₁ is represented as e₁=[1 0 . . . 0]^T Or with a dimension of 1×q;

x_j^q represents q_th order reduced solution vector under excitation at port j;

{tilde over (G)}_i represents a derivative of G with respect to ϕ_i;

{tilde over (C)}_i represents a derivative of C with respect to ϕ_i;

x_k^qT represents a transpose vector of qth order reduced field solution under excitation at port k;

S_k,j represents the S-parameter of port j and port k;

ϕ_i represents an i_th one of the physical parameters; and

κ_k,j represents a correlated coefficient of port j and port k.

The following advantages are associated with the electromagnetic sensitivity analysis method of the disclosure:

The proposed single-size MPVL method transforms an EM linear system to be solved into a one of single-size matrices. And the adjoint/self-adjoint EM sensitivity analysis for fast frequency sweep are generated to calculate derivative of S-parameter of the multiple ports of the EM device, allowing for direct application of MPVL method for model order reduction calculation in EM sensitivity analysis. When fast frequency sweep using single-size MPVL method is performed for EM sensitivity analysis, repetitive solving process due to different frequencies is avoided and computational cost of the EM device design system is saved, leading to a sufficiently enhanced simulation efficiency.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flow chart for the operations performed by the invented device design system in the embodiment of this microwave filter.

FIG. 2 is a perspective view of a four-pole waveguide filter according to one embodiment of the disclosure;

FIG. 3 is another perspective view of a four-pole waveguide filter according to one embodiment of the disclosure;

FIG. 4 is graph depicting Derivatives of ReS₁₁ and ImS₁₁ with respect to h₁ at ϕ₁ for the four-pole waveguide filter with the frequency swept from 10.5 GHz-11.5 GHz;

FIG. 5 is graph depicting Derivatives of ReS₂₁ and ImS₂₁ with respect to h₁ at ϕ₁ for the four-pole waveguide filter with the frequency swept from 10.5 GHz-11.5 GHz;

FIG. 6 is graph depicting Derivatives of ReS₁₁ and ImS₁₁ with respect to h_c1 at ϕ₁ for the four-pole waveguide filter with the frequency swept from 10.5 GHz-11.5 GHz; and

FIG. 7 is graph depicting Derivatives of ReS₂₁ and ImS₂₁ with respect to h_c2 at ϕ₁ for the four-pole waveguide filter with the frequency swept from 10.5 GHz-11.5 GHz.

DETAILED DESCRIPTION

For an EM device design system, an adjoint formula and a self-adjoint formula that are respectively integrated with fast frequency sweeps are generated to improve EM sensitivity analysis for EM device designs. Further, a single-size MPVL method is proposed for model order reduction and is used in fast frequency sweep for the present EM sensitivity analysis. The present EM sensitivity analysis method obtains the same accuracy as the techniques of the prior art, while taking much less time by avoiding repetitively solving large systems of EM equations for different frequencies.

The single-size MPVL method incorporated in fast frequency sweep is derived as follows:

the Helmholtz equation for the full wave EM simulation is formulated as:

$\begin{array}{cc}\nabla \times \left(\frac{1}{\mu }\nabla \times E\right)-{\omega }^2ϵE=-j\omega J,& \left(1\right)\end{array}$

where ∈ and μ represent the permittivity and the permeability of the medium, respectively; ω represents the angular frequency; J represents the electric current source of the EM problem; E represents the electric field intensity to be solved; finite element method (FEM) is one of the commonly used methods to solve this Helmholtz equation; assuming that the EM simulation is performed on an EM device comprising multiple ports; the FEM equation in a general form is written as:

(K₀+sK₁+s²K₂)x=sb_j,  (2)

where s represents the complex frequency corresponding to ω; x represents the unknown field solution vector containing the unknown values of E for FEM; K₀, K₁, and K₂ represent global finite-element system matrices, which are dependent on the physical parameters of the EM device but independent on frequency; b_j represents a vector describing the excitation at the port j; to obtain the full scattering matrix of the multi-port structure, the S-parameters can be written as

$\begin{array}{cc}\begin{array}{c}{S}_{k,j}\left(s\right)={\kappa }_{k,j}{b}_k^{T}x-{\delta }_{k,j}\\ =s{\kappa }_{k,j}{b_k^T\left(K_0+{\mathrm{sK}}_1+s^2{K}_2\right)}^{-1}{b}_j-{\delta }_{k,j},\end{array}& \left(3\right)\end{array}$ $\mathrm{where}$ $\begin{array}{cc}{\delta }_{k,j}=\{\begin{array}{cc}1,& k=j\\ 0,& k\ne j\end{array},& \left(4\right)\end{array}$

κ_k,j is a correlated coefficient of port j and port k, which is constant for a constant EM device and is dependent on the powers incident upon the port k and port j of the EM device.

To use MPVL method on FEM system, the order of frequency terms contained in the EM equations would be no more than a first order. Therefore, Equation (3) is reformatted as:

$\begin{array}{cc}{S}_{k,j}\left(s\right)={{\hat{\beta }}_{k,j}^T\left(G+\mathrm{sC}\right)}^{-1}{\beta }_j,& \left(5\right)\end{array}$ $\mathrm{where}$ $\begin{array}{cc}{\hat{\beta }}_{k,j}=\left[\begin{array}{c}s{\kappa }_{k,j}{b}_k\\ 0_{N\times 1}\end{array}\right],G=\left[\begin{array}{cc}{0}_N& I_N\\ K_0& K_1\end{array}\right],C=\left[\begin{array}{cc}-I_N& 0_N\\ 0_N& K_2\end{array}\right],{\beta }_j=\left[\begin{array}{c}{0}_{N\times 1}\\ b_j\end{array}\right].& \left(6\right)\end{array}$

where I_N and 0_N represent the identity matrix and zero matrix in ^N×N, respectively. Here, N is the number of elements in x. Let ϕ represent a vector of the physical parameters of the EM device, where ϕ_i (i=1, . . . , p) represent the i_th element in ϕ, and p is the total number of physical parameters.

To use MPVL method to perform the fast frequency sweep, one specific frequency is selected and defined as the pre-solution frequency; with the suitable selection of the pre-solution frequency, the accurate EM solutions can be obtained over the entire frequency band by solving the large linear system once at the pre-solution frequency; let s₀ represent the pre-solution frequency for MPVL method; let K_s represent a single-size system matrix at the pre-solution frequency, which has a dimension of N×N and is written as

K _s =K ₀ +s ₀ K ₁ +s ₀ ² K ₂,  (7)

let A_s represent a double-size system matrix having a dimension of 2N×2N and is defined as

$\begin{array}{cc}{A}_s=G+s_{0}C=\left[\begin{array}{cc}-s_0{I}_N& I_N\\ K_0& K_1+s_0{K}_2\end{array}\right].& \left(8\right)\end{array}$

where A_s is calculated at the pre-solution frequency s₀; let q be the order of the reduced model using MPVL; let T_k^q be defined as the reduced order matrix with the elements t_l,m, i.e.,

$T_k^q={\left[t_{l,m}\right]}_{q\times q},$

where l=1, 2, . . . , q and m=1, 2, . . . , q; let {w_m}_m=1^q be defined as a set of unit vectors which are orthonormal to each other; let V: be defined as the matrix containing the orthonormal basis of the Krylov subspaces for the model order reduction, where V_k^q=[v_m]_m=1^q=[v₁ v₂ . . . v_q]; let r_m be defined as the vector for the calculation of V_k^q and T_k^q in the m_th MPVL iteration; r_m is calculated by solving a first linear system as:

$\begin{array}{cc}{A}_s{r}_m=\{\begin{array}{cc}{\beta }_k,& m=0,\\ -{\mathrm{Cv}}_m,& m=1,2,\dots ,q.\end{array}& \left(9\right)\end{array}$

using MPVL as the model order reduction technique to perform the fast frequency sweep; calculating the LU factors of the matrix A_s;

because the size of A_s is twice the size of K_s, the computational cost for the LU factorization of A_s takes much longer time than that of Ks;

therefore, a single-size MPVL method is derived to perform the LU factorization of the single-size system matrix K_s, instead of the double-size system matrix A_s; using the block matrix inversion method, Equation (9) is reformulated as:

$\begin{array}{cc}\begin{array}{c}{r}_0=\frac{1}{s_0}\left[\begin{array}{cc}{K}_s^{-1}{K}_0-I& s_0{K}_s^{-1}\\ s_0{K}_s^{-1}{K}_0& s_0^2{K}_s^{-1}\end{array}\right]\left[\begin{array}{c}{0}_{N\times 1}\\ b_k\end{array}\right]\\ =\left[\begin{array}{c}{I}_N\\ s_0{I}_N\end{array}\right]K_s^{-1}{b}_k,\end{array}& \left(10\right)\end{array}$ $\mathrm{and}$ $\begin{array}{cc}\begin{array}{c}{r}_m=\frac{1}{s_0}\left[\begin{array}{cc}{K}_s^{-1}{K}_0-I& s_0{K}_s^{-1}\\ s_0{K}_s^{-1}{K}_0& s_0^2{K}_s^{-1}\end{array}\right]\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& -K_2\end{array}\right]v_m\\ =\frac{1}{s_0}\left[\begin{array}{c}{I}_N\\ s_0{I}_N\end{array}\right]{K_s^{-1}\left[\begin{array}{c}{K}_0\\ -s_0{K}_2\end{array}\right]}^T{v}_m\\ -\frac{1}{s_0}\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& 0_N\end{array}\right]v_{m,}\end{array}& \left(11\right)\end{array}$

let u_m be defined as the solution vector for a second linear system as

$\begin{array}{cc}{K}_s{u}_m=\{\begin{array}{cc}{b}_k,& m=0\\ {\frac{1}{s_0}\left[\begin{array}{c}{K}_0\\ -s_0{K}_2\end{array}\right]}^T{v}_m,& m=1,2,\dots ,q.\end{array}& \left(12\right)\end{array}$

substituting (12) into (10) and (11),

$\begin{array}{cc}{r}_0=\left[\begin{array}{c}{u}_0\\ s_0{u}_0\end{array}\right],& \left(13\right)\end{array}$ $\mathrm{and}$ $\begin{array}{cc}{r}_m=\left[\begin{array}{c}{u}_m\\ s_0{u}_m\end{array}\right]-\frac{1}{s_0}\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& 0_N\end{array}\right]v_m,& \left(14\right)\end{array}$

where m represents the MPVL iteration index, i.e., m=1, 2, . . . , q; let x_k^q represent the q_th order reduced vector (q×1 vector), the solution vector x_k under excitation at port k of the device can be formulated as,

x _k ≈x _k ^q =∥r ₀ ∥V _k ^q(I_q−(s−s₀)T_k^q)⁻¹e₁,  (15)

where e₁[1 0 . . . 0]^T∈^q, is the unit vector in ^q; I_q is an identity matrix in ^q×q.

according to Equation (3), the S-parameters can be obtained by substituting the solution vector x_k with the Equation (15) and be formulated as

S _k,j(S)=∥r₀∥{circumflex over (β)}_k,j^TV_k^q(I_q−(s−s₀)T_k^q)⁻¹e₁−δ_k,j,  (16)

where, S-parameters corresponding to a frequency range can be evaluated by performing fast frequency sweep integrated with the single-size MPVL method.

The following is derivation of the adjoint EM sensitivity analysis using fast frequency sweep model order reduction technique.

Letting {tilde over (G)}_i and {tilde over (C)}_i represent the derivatives of G and C with respect to a physical parameter of the device ϕ_i, respectively. {tilde over (G)}_i and {tilde over (C)}_i can then be derived as,

$\begin{array}{cc}{\tilde{G}}_i=\left[\begin{array}{cc}{0}_N& 0_N\\ {\tilde{K}}_0^i& {\tilde{K}}_1^i\end{array}\right],{\tilde{C}}_i=\left[\begin{array}{cc}{0}_N& 0_N\\ 0_N& {\tilde{K}}_2^i\end{array}\right].& \left(17\right)\end{array}$

where {tilde over (K)}₀ⁱ, {tilde over (K)}₁ⁱ, and {tilde over (K)}₂ⁱ represent the sensitivity of K₀, K₁, and K₂ with respect to the physical parameter ϕ_i, respectively; to calculate the sensitivity of S_{i, j}, the sensitivities of the global finite-element system matrices w.r.t. physical parameters {tilde over (K)}₀ⁱ, {tilde over (K)}₁ⁱ, and {tilde over (K)}₂ⁱ need to be firstly obtained;

obtaining the vector b_j from the EM excitation and/or the boundary conditions on the j_th port; the structures of the ports normally do not change during the EM design; thus b_j is a constant vector w.r.t. to a physical parameter ϕ_i in the EM problem, i.e.,

$\begin{array}{cc}\frac{\partial b_j}{\partial {\varphi }_i}=0_{N\times 1}.& \left(18\right)\end{array}$

based on (5), formulating the derivatives of S_k,j as,

$\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}=-{{\hat{\beta }}_{k,j}^T\left(G+\mathrm{sC}\right)}^{-1}\left({\tilde{G}}_i+s{\tilde{C}}_i\right){\left(G+\mathrm{sC}\right)}^{-1}{\beta }_j.& \left(19\right)\end{array}$

letting x_j represent the solution vector for a linear system that is written as

(G+sC)x_j=β_j.  (20)

calculating the solution vector x_j in (20) by MORe techniques such as MPVL method;

defining {circumflex over (x)}_k,j as a vector of adjoint field solution of the adjoint representation of the linear system (20), the adjoint representation is written as follows,

(G+sC)^T{circumflex over (x)}_k,j={circumflex over (β)}_k,j.  (21)

the adjoint solution vector {circumflex over (x)}_k,j in (21) is calculated similarly to x_j by MORe techniques such as MPVL method;

substituting (20) and (21) into (19) to obtain the adjoint sensitivity formula of S_k,j w.r.t. ϕ_i as,

$\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}={\hat{x}}_{k,j}^T\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j,& \left(22\right)\end{array}$

Equation (22) is the adjoint sensitivity formula which is well fit for fast frequency sweep. MORe technique can be used to calculate x_j and {circumflex over (x)}_k,j; since the format of (22) is derived using MPVL, x_j and {circumflex over (x)}_k,j can be formulated using MPVL method,

$\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\hat{x}}_{k,j}^{qT}\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)x_j^q.& \left(23\right)\end{array}$

The following is derivation of the self-adjoint EM sensitivity analysis using fast frequency sweep model order reduction technique.

The disclosure further provides a self-adjoint EM sensitivity analysis. Because G and C are formulated as 2×2 block matrices, by substituting (6) into (20) and using the 2×2 block matrix inversion method, we can obtain,

$\begin{array}{cc}\begin{array}{c}{x}_k=\frac{1}{s}\left[\begin{array}{cc}{K}^{-1}{K}_0-I_N& s{K}^{-1}\\ s{K}^{-1}{K}_0& s^2{K}^{-1}\end{array}\right]\left[\begin{array}{c}{0}_{N\times 1}\\ b_k\end{array}\right]\\ =\left[\begin{array}{c}{K}^{-1}{b}_k\\ s{K}^{-1}{b}_k\end{array}\right],\end{array}\mathrm{where}& \left(24\right)\end{array}$ $\begin{array}{cc}K=K_0+s{K}_1+s^2{K}_2.& \left(25\right)\end{array}$

because K₀, K₁, and K₂ are symmetrical matrices, substituting (6) into (21) and using the 2×2 block matrix inversion method to obtain,

$\begin{array}{cc}\begin{array}{c}{\hat{x}}_{k,j}=\frac{1}{s}\left[\begin{array}{cc}{K}^{-1}{K}_0-I_N& s{K}^{-1}{K}_0\\ s{K}^{-1}& s^2{K}^{-1}\end{array}\right]\left[\begin{array}{c}s{\kappa }_{k,j}{b}_k\\ 0_{N\times 1}\end{array}\right]\\ ={\kappa }_{k,j}\left[\begin{array}{c}{K}^{-1}{K}_0{b}_k-b_k\\ s{K}^{-1}{b}_k\end{array}\right].\end{array}& \left(26\right)\end{array}$

defining ξ_k as the vector of weighted differences between {circumflex over (x)}_k,j and x_k, formulated as,

$\begin{array}{cc}{\xi }_k=\frac{{\hat{x}}_{k,j}}{{\kappa }_{k,j}}-x_k=\left[\begin{array}{c}{K}^{-1}\left(K_0{b}_k-b_k\right)-b_k\\ 0_{N\times 1}\end{array}\right].& \left(27\right)\end{array}$

formulating {circumflex over (x)}_k,j by the weighted summation of x_k and ξ_k as,

{circumflex over (x)} _k,j=κ_k,jx_k+κ_k,jξ_k,  (28)

substituting (17) into (22) to obtain,

$\begin{array}{cc}\begin{array}{c}\frac{\partial S_{k,j}}{\partial {\varphi }_i}=-{\kappa }_{k,j}\left(x_k^T+{\xi }_k^T\right)\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)x_j\\ =-{\kappa }_{k,j}\left[x_k^T\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)+{\xi }_k^T\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)\right]x_j.\end{array}& \left(29\right)\end{array}$

because

$\begin{array}{cc}\begin{array}{c}{\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)}^T{\xi }_k=\left[\begin{array}{cc}{0}_N& {\stackrel{\sim }{K}}_0^i\\ 0_N& {\stackrel{\sim }{K}}_1^i+s{\stackrel{\sim }{K}}_2^i\end{array}\right]\left[\begin{array}{c}{K}^{-1}\left(K_0{b}_k-b_k\right)-b_k\\ 0_{N\times 1}\end{array}\right]\\ =0_{2N\times 1},\end{array}& \left(30\right)\end{array}$

simplifying the formulation for calculating the derivative (29) as,

$\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}=-{\kappa }_{k,j}{x}_k^T\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)x_j.& \left(31\right)\end{array}$ $\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\kappa }_{k,j}{x}_k^{qT}\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)x_j^q.& \left(32\right)\end{array}$

Equation (20) is the self-adjoint sensitivity formula which is well fit for fast frequency sweep; MORe techniques can be used to calculate x_j; since the format of (31) is derived using MPVL, x_j can be formulated using MPVL method.

By using Equations (12)-(14), the disclosure performs the LU factorization and forward/backward substitutions of the original system K_s (single-size system matrix) instead of A_s (double-size system matrix), to reduce calculation time r_m, m=0, 1, . . . , q. Because the MPVL method takes a long time to calculate r_m, it is simplified for efficient x_k calculation. Note that, according to Equations (12)-(14), the disclosure uses only one LU factorization of the original system K_s, whereas the MPVL method uses q+1 forward/backward substitutions.

Because Equation (31) is derived using MPVL, x_j and {circumflex over (x)}_k,j can also be formulated using the MPVL method,

$\begin{array}{cc}\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\kappa }_{k,j}{x}_k^{qT}\left({\stackrel{\sim }{G}}_i+s{\tilde{C}}_i\right)x_j^q,& \left(32\right)\end{array}$

where x_j^q is calculated using fast frequency sweep with MORe technique as in equation (15). Note that, in the single-size MPVL method, {w_m}_m=1^q makes each λ_m equal to non-zero. When performing MPVL method, Equation (32) calculates the derivatives of S-parameters for the self-adjoint EM sensitivity analysis. Fast frequency sweep uses a reduced-order model of order q to estimate the S-parameters of a p port structure and the derivatives thereof, and involves only one LU factorization and p (q+1) forward/backward substitutions. The computational time of the self-adjoint EM sensitivity analysis is not related to the number of physical parameters and increases slightly with the number of frequencies to be solved.

The comparison of the different methods of EM sensitivity analysis of a four-pole waveguide filter is shown in Table 1. A flow chart for the operations performed by the invented device design system in the embodiment of this microwave filter is shown in FIG. 1. The detailed number of LU factorizations and FB substitutions of the original FEM matrix for the proposed adjoint EM sensitivity analysis using fast frequency sweep is illustrated in Table 1. The disclosure further provides details of the self-adjoint EM sensitivity analysis using discrete frequency sweep, the sensitivity analysis using the finite difference method (i.e., perturbing variables, re-simulating EM and computing the solution difference) for both fast frequency sweep and discrete frequency sweep. For illustration purpose, the disclosure further provides the details of EM simulation using both fast frequency sweep and discrete frequency sweep. In Table 1, n represents the number of physical parameters; n_f represents the number of frequencies; q represents the order for MPVL method; p represents the number of ports of the microwave filter. From Table 1, if the EM application only simulates at 1 or 2 frequencies, the sensitivity analysis using self-adjoint EM sensitivity analysis with discrete frequency sweep is performed. If the EM application only has 1 or 2 physical parameters, the sensitivity analysis using finite difference method with fast frequency sweep is performed. If the EM application has more physical parameters and needs to simulate at more frequencies, the adjoint EM sensitivity analysis using fast frequency sweep obtains the results faster than the existing sensitivity analysis methods. The larger number of physical parameters and frequencies the EM application has, the more significant advantage the proposed method will have.

TABLE 1
Comparison of the Number of LU Factorizations and F/B
Substitutions of FEM Matrix for Different Methods
		No. of LU	No. of F/B
	Method	Factorizations	Substitutions
	Proposed Self-Adjoint	1	p(q + 1)
	EM Sensitivity Analysis
	Using Fast Frequency
	Sweep
	Proposed Adjoint EM	1	p(2q + 2)
	Sensitivity Analysis
	Using Fast Frequency
	Sweep
	Self-Adjoint EM	nf	p · nf
	Sensitivity Analysis
	Using Discrete Frequency
	Sweep
	Finite-Difference	n + 1	p(q + l)(n + 1)
	Sensitivity Using Fast
	Frequency Sweep
	Finite-Difference	nf(n + 1)	p · nf(n + 1)
	Sensitivity Analysis
	Using Discrete Frequency
	Sweep
	Fast Frequency Sweep	1	p(q + 1)
	Simulation
	Discrete Frequency	nf	p • nf
	Sweep Simulation

n—Number of Physical Parameters;

n_f- Number of Frequencies at a Frequency range;

q—Order for MPVL method;

p—Number of ports of the microwave filter.

Remarks: Note that, the self-adjoint EM sensitivity method takes the same number of LU decomposition and F/B substitutions as EM simulation using fast frequency sweep using finite element method. In another word, after performing EM simulation using fast frequency sweep, no extra LU or F/B time is needed to calculate the EM sensitivities wr.t. all variables. Only a small amount of extra time of matrix multiplications is needed for EM sensitivity analysis after EM simulation.

4. Sensitivity Analysis for Four-Pole Waveguide Filter

The example under consideration is the EM sensitivity analysis of a four-pole waveguide filter. The tuning elements are penetrating posts of square cross section placed at the center of each cavity and each coupling window, shown in FIG. 2. Let h₁, h₂, h₃ represent the heights of the tuning posts in the coupling windows. Let h_c1 and h_c2 represent the heights of the square cross section placed in the center of the resonator cavities. Let w_c and w_p represent the widths of the square cross section placed at the center of the resonant cavity and the coupling window respectively. Let d₁₂ and d₂₃ represent the diameters of resonant cavities. The iris widths w₁, w₂, w₃ and w₄ are pre-tuned to achieve the desired coupling bandwidth for each section of the filter. The input and output waveguides, as well as the resonant cavities, are standard WR-75 waveguides (a=19.05 mm and b=9.525 mm). For this example, the design space vector is ϕ=[d₁₂ d₂₃ h₁ h₂ h₃ h_c1 h_c2 w₁ w₂ w₃ W_p W_c]^T The sensitivity analysis is performed at the nominal design point, i.e., (P_i=[10.5 13.3 3.39 4.15 3.62 3.30 2.98 8.7 5.1 5.1 2.0 4.0]^T mm. The frequency range for this filter is selected as 10.5 GHz-11.5 GHz with the step of 20 MHz.

The derivatives of S-parameters for the microwave filter device is calculated by the proposed self-adjoint EM sensitivity analysis using fast frequency sweep for the EM simulation of this filter example. The order for the reduced order model is 16, i.e., q=16. Through the instant EM device design system disclosed in this embodiment for the microwave filter, physical specifications, such as passband bandwidth, return loss, and passband ripple can be satisfied.

For comparison purposes, the disclosure further provides the self-adjoint EM sensitivity analysis using discrete frequency analysis, the sensitivity analysis using the finite difference method for fast frequency sweep and discrete frequency analysis. Three cases for the sensitivity analysis with different number of frequency points are performed for the comparison: Case 1 with 11 frequency points and Case 2 with 51 frequency points. The comparison of the different methods of sensitivity analysis for the four-pole waveguide filter example is shown in Table 2.

TABLE 2
Comparison of CPU of different sensitivity analysis methods for the four-
pole waveguide filter
		No. of LU	No. of F/B	Time of LU	Total
	Sensitivity method	factorizations	substitutions	Factorizations	Time
11 Frequencies	Proposed Self-	1	34	2.1 min	4.1
	adjoint EM				min
	sensitivity analysis
	Using Fast
	frequency sweep
	Proposed Adjoint	1	68	2.1 min	6.2
	EM Sensitivity				min
	Analysis Using Fast
	frequency sweep
	Self-adjoint EM	11	22	22.3 min	23.9
	sensitivity analysis				min
	Using Discrete
	Frequency Analysis
	Finite-Difference	13	442	26.4 min	41.2
	Sensitivity Analysis				min
	Using Fast
	frequency sweep
	Finite-Difference	143	286	4.7 h	4.9 h
	Sensitivity Analysis
	Using Discrete
	Frequency Sweep
101 Frequencies	Proposed Self-	1	34	2.1 min	4.1
	adjoint EM				min
	sensitivity analysis
	Using Fast
	Frequency Sweep
	Proposed Adjoint	1	68	2.1 min	6.2
	Sensitivity Analysis				min
	Using Fast
	frequency sweep
	Self-adjoint EM	102	204	3.3 h	3.4 h
	sensitivity analysis
	Using Discrete
	Frequency Sweep
	Finite-Difference	13	442	26.4 min	41.2
	Sensitivity Analysis				min
	Using Fast
	Frequency Sweep
	Finite-Difference	1326	2652	44.3 h	45.1
	Sensitivity Analysis				h
	Using Discrete
	Frequency Sweep

Because 12 physical parameters are used in this example, 13 complete evaluations of the S-parameters (once at nominal values and 12 perturbations for 12 different physical parameters) are needed for the sensitivity analysis using the finite different method. For the self-adjoint EM sensitivity analysis using discrete frequency analysis, the cost increases as the number of frequency points for the EM design increases. From Table 2, the self-adjoint EM sensitivity analysis uses only one LU factorization and 34 forward/backward substitutions to obtain the derivatives of the S-parameters w.r.t 12 different physical parameters. The LU factorizations and forward/backward substitutions are the most time-consuming part during the overall sensitivity analysis process. The number of LU factorizations and forward/backward substitutions used by the proposed method is much fewer than that used by the existing methods listed in Table 2, therefore, the disclosed method takes much less time than the existing methods.

FIG. 4 and FIG. 5 show detailed comparisons of the derivatives of S₁₁ and S₂₁ with respect to h₁ for the four-pole waveguide filter with the frequency swept from 10.5 GHz-11.5 GHz, respectively. As shown in FIG. 3 and FIG. 4, the derivatives calculated using the self-adjoint EM sensitivity analysis match very well with the derivatives calculated using the existing self-adjoint method and finite difference methods. To further verify the accuracy of the proposed self-adjoint EM sensitivity analysis method, the disclosure further provides the calculation of the derivatives of S-parameters w.r.t. other physical parameters. Detailed comparisons of the derivatives of S₁₁ and S₂₁ with respect to h_c2 for this four-pole waveguide filter are shown in FIG. 6 and FIG. 7, respectively. From the comparisons, the self-adjoint EM sensitivity analysis using fast frequency sweep achieves the correct solution with a very low computational cost.

Appendix 1: Single-Size MPVL Method for EM sensitivity Analysis

1) Define δ_k,j, κ_k,j, e₁, I_q, K₀, K₁, K₂, b_k, b_j. Define frequency s. Define pre-solution frequency s₀.

Define {w_m}_m=1^q. Set K_s=K₀+s₀K₁+s₀²K₂

2) Obtain u₀ by solving the linear system K_su₀=b₁. Set

$r_0=\left[\begin{array}{c}{u}_0\\ s_0{u}_0\end{array}\right]\mathrm{and}r=r_0.$

For m=1, 2, q do:

3) Set

$v_m=\frac{r}{r}\mathrm{and}{\lambda }_m=w_m^T{v}_m.$

4) If m>1, set

t _m,m−1 =∥r∥.

5) Solving the linear system to obtain um,

$K_s{u}_m={\frac{1}{S_0}\left[\begin{array}{c}{K}_0\\ -s_0{K}_2\end{array}\right]}^T{v}_m.$

Calculate

$r=r_m=\left[\begin{array}{c}{u}_m\\ s_0{u}_m\end{array}\right]-\frac{1}{S_0}\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& 0_N\end{array}\right]v_m.$

6) For l=1, 2, . . . , m do: Set

$t_{l,m}=\frac{w_l^{T}r}{{\lambda }_l}\mathrm{and}r=r-t_{l,m}{v}_l.$

End For

End For

7) Set

T _j ^q=[t_l,m]_q×q,V_j^q=[v₁v₂ . . . v_q].

x _j ≈∥r ₀ ∥V _j ^q(I_q−(s−s₀)T_j^q)⁻¹e₁.

8) Calculate

$S_{k,j}\left(s\right)\approx sr_0{{\kappa }_{k,j}\left[\begin{array}{c}{b}_k\\ 0_{N\times 1}\end{array}\right]}^T{V_j^q\left(I_q-\left(s-s_0\right)T_j^q\right)}^{-1}{e}_1-{\delta }_{k,j}.$

Appendix 2: Adjoint Electromagnetic Sensitivity Analysis using Single-Size MPVL Method

1) Define κ_k,j, e₁, I_q, {tilde over (G)}_i, {tilde over (C)}_i, K₀, K₁, K₂, {circumflex over (β)}_k,j, and b_j.

Define frequency s. Define pre-solution frequency s₀.

Define {w_m}_m=1^q. Set K_s=K₀+s₀K₁+s₀²K₂.

2) Obtain u₀ by solving the linear system K_su₀=b₁.

Set

$r_0=\left[\begin{array}{c}{u}_0\\ s_0{u}_0\end{array}\right]\mathrm{and}r=r_0.$

For m=1, 2, . . . , q do:

3) Set

$v_m=\frac{r}{r}\mathrm{and}{\lambda }_m=w_m^T{v}_m.$

4) If m>1, set

t _m,m-1 =∥r∥.

5) Solving the linear system to obtain um,

$K_s{u}_m={\frac{1}{s_0}\left[\begin{array}{c}{K}_0\\ -s_0{K}_2\end{array}\right]}^T{v}_m.$

Calculate

$r=r_m=\left[\begin{array}{c}{u}_m\\ s_0{u}_m\end{array}\right]-\frac{1}{s_0}\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& 0_N\end{array}\right]v_{m.}$

6) For l=1, 2, . . . , m do: Set

$t_{l,m}=\frac{w_l^{T}r}{{\lambda }_l}\mathrm{and}r=r-t_{l,m}{v}_l.$

End For

End For

7) Set

T _k ^q=[t_l,m]_q−q,V_k^q=[v₁v₂ . . . v_q],

x _j ^q =∥r ₀ ∥V _j ^q(I_q−(s−s₀)T_j^q)⁻¹e₁.

8) Repeat 2)-7) by changing b_j to {circumflex over (β)}_k,j to calculate {circumflex over (x)}_k,j^qT

9) Calculate

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\hat{x}}_{k,j}^{q^T}\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j^q.$

Appendix 3: Self-Adjoint Electromagnetic Sensitivity Analysis using Single-Size MPVL Method

1) Define κ_k,j, e₁, I_q, {tilde over (G)}_i, {tilde over (C)}_i, K₀, K₁, K₂, b_k, and b_j.

Define frequency s. Define pre-solution frequency s₀.

Define {w_m}_m=1^q. Set K_s=K₀+s₀K₁+s₀²K₂.

2) Obtain u₀ by solving the linear system K_su₀=b_j.

Set

$r_0=\left[\begin{array}{c}{u}_0\\ s_0{u}_0\end{array}\right]\mathrm{and}r=r_0.$

For m=1, 2, . . . , q do:

3) Set

$v_m=\frac{r}{r}\mathrm{and}{\lambda }_m=w_m^T{v}_m.$

4) If m>1, set

t _m,m-1 =∥r∥.

5) Solving the linear system to obtain um,

$K_s{u}_m={\frac{1}{s_0}\left[\begin{array}{c}{K}_0\\ -s_0{K}_2\end{array}\right]}^T{v}_m.$

Calculate

$r=r_m=\left[\begin{array}{c}{u}_m\\ s_0{u}_m\end{array}\right]-\frac{1}{s_0}\left[\begin{array}{cc}{I}_N& 0_N\\ 0_N& 0_N\end{array}\right]v_m.$

6) For l=1, 2, . . . , m do: Set

$t_{l,m}=\frac{w_l^{T}r}{{\lambda }_l}\mathrm{and}r=r-t_{l,m}{v}_l.$

End For

End For

7) Set

T _j ^q=[t_l,m]_q×q,V_j^q=[v₁v₂ . . . v_q],

x _j ^q =∥r ₀ ∥V _j ^q(I_q−(s−s₀)T_j^q)⁻¹e₁.

8) Repeat 2)-7) by changing b_j to b_k to calculate x_k^q.

9) Calculate

$\frac{\partial S_{k,j}}{\partial {\varphi }_i}\approx -{\kappa }_{k,j}{x}_k^{q^T}\left({\tilde{G}}_i+s{\tilde{C}}_i\right)x_j^q.$

Claims

1. An electromagnetic device design system comprising: one or more processors of a machine; andcomputer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising: initiating physical parameters of the EM device, wherein the EM device comprises multiple ports;performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);applying single-size matrix Padé via Lanczos (MPVL) method in fast frequency sweep and performing EM simulation for the EM device under excitation at each port to obtain field solutions of the EM device in a frequency range;calculating S-parameters for the multiple ports of the EM device;obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with a selected solution;performing EM simulation for the EM device at a frequency of the selected solution using FEM; anddetermining the simulated result satisfies a physical specification of the EM device;wherein: applying single-size MPVL method in fast frequency sweep comprises:transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents a number of elements in a field vector and the single-size system matrix is of a linear combination of global finite-element system matrices;generating a first linear system using the double-size system matrix;representing solving vectors of the first linear system with the global finite-element system matrices using the block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; andsolving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining field solution in a frequency range as the following: xk≈xkq=∥r0∥vkq(Iq−(s−s0)Tkq)−1e1 wheres represents a frequency;s0 represents a pre-solution frequency;q is a reduced order in MPVL method;xk is a vector of field solution under excitation at port k;xkq represents a vector of qth order reduced field solution under excitation at port k;r0 represents a solution vector of the first linear system at 0th MPVL iteration;Vkq is represented as Vkq=[vm]m=1q=[v1 v2 . . . vq], where vm is an orthonormal basis vector of the Krylov subspace for model order reduction;Iq is an identity matrix with a dimension of q×q;Tkq, is a reduced order matrix; ande1 is represented as e1=[1 0 . . . 0]T with a dimension of 1×q.

2. An electromagnetic device design system comprising: one or more processors of a machine; andcomputer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising: initiating physical parameters of the EM device wherein the EM device comprises multiple ports;performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);solving a linear system under excitation at port j to obtain a field solution under excitation at port j by fast frequency sweep;solving an adjoint representation of the linear system under excitation at port k to obtain an adjoint field solution between port k and port j by fast frequency sweep;calculating derivatives of a S-parameter of port j and port k with respect to ϕi based on an adjoint sensitivity formula;obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;selecting a solution and updating EM device design to replace the one of the physical parameters of the EM device with the selected solution;performing EM simulation for the EM device at a frequency of the selected solution using FEM; anddetermining the simulated result satisfies a physical specification of the EM device;wherein the adjoint sensitivity formula is written as

3. An electromagnetic device design system comprising: one or more processors of a machine; andcomputer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising: initiating physical parameters of the EM device wherein the EM device comprises multiple ports;performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);solving a linear system under excitation at port j to obtain a field solution under excitation at port j by fast frequency sweep;solving the linear system under excitation at port k to obtain a field solution under excitation at port k by fast frequency sweep;calculating derivatives of a S-parameter of port j and port k with respect to ϕi based on a self-adjoint sensitivity formula;obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;performing EM simulation for the EM device at a frequency of the selected solution using FEM; anddetermining the simulated result satisfies a physical specification of the EM device;whereinthe self-adjoint sensitivity formula is written as

4. An electromagnetic device design system comprising: one or more processors of a machine; andcomputer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising: initiating physical parameters of the EM device wherein the EM device comprises multiple ports;performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);solving a linear system under excitation at port j by performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port j in a frequency range;solving an adjoint representation of the linear system under excitation at port k by performing fast frequency sweep with single-size MPVL method to obtain an order reduced adjoint field solution between port k and port j in the frequency range;calculating derivatives of a S-parameter of port j and port k with respect to ϕi based on an adjoint sensitivity formula;obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;performing EM simulation for the EM device at a frequency of the selected solution using FEM; anddetermining the simulated result satisfies a physical specification of the EM device;wherein: performing fast frequency sweep with single-size MPVL method comprises:transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents the number of elements in a field vector and the single-size system matrix comprises a linear combination of global finite-element system matrices;generating a first linear system using the double-size system matrix;representing solving vectors of the first linear system with the global finite-element system matrices using block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; andsolving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining a field solution in a frequency range by the following: xk≈xkq=∥r0∥vkq(Iq−(s−s0)Tkq)−1e1;and the adjoint sensitivity formula is written as:

5. An electromagnetic device design system comprising: one or more processors of a machine; andcomputer-storage medium storing instructions, which when executed by the machine, cause the machine to perform operations for EM sensitivity analysis for an electromagnetic (EM) device, the operations comprising: initiating physical parameters of the EM device wherein the EM device comprises multiple ports;performing EM simulation for the EM device at a pre-solution frequency using the finite-element method (FEM);performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port j of the EM device in a frequency range;performing fast frequency sweep with single-size MPVL method to obtain an order reduced field solution under excitation at port k of the EM device in a frequency range;calculating derivatives of a S-parameter of port j and port k with respect to ϕi based on a self-adjoint sensitivity formula;obtaining derivatives of a full scattering matrix for the multiple ports with respect to the physical parameters of the EM device in the frequency range;selecting a solution and updating EM device design to replace initial values of the physical parameters of the EM device with the selected solution;performing EM simulation for the EM device at a frequency of the selected solution using FEM; anddetermining the simulated result satisfies a physical specification of the EM device;wherein: performing fast frequency sweep with single-size MPVL method comprises:transforming a single-size system matrix with a dimension of N×N into a double-size system matrix with a dimension of 2N×2N for omitting second order terms of frequency, wherein N represents the number of elements in a field vector and the single-size system matrix comprises a linear combination of global finite-element system matrices;generating a first linear system using the double-size system matrix;representing solving vectors of the first linear system with the global finite-element system matrices using block matrix inversion method and transforming the first linear system into a second linear system, wherein the second linear system is of the single-size system matrix; andsolving the second linear system by performing fast frequency sweep incorporated with MPVL method and obtaining field solution of the EM device in a frequency range as the following: xk≈xkq=∥r0∥Vkq(Iq−(s−s0)Tkq)−1e1,and the self-adjoint sensitivity formula is written as