CN108536929B - 一种应用arpack求解波导结构色散特性的方法 - Google Patents
一种应用arpack求解波导结构色散特性的方法 Download PDFInfo
- Publication number
- CN108536929B CN108536929B CN201810243782.XA CN201810243782A CN108536929B CN 108536929 B CN108536929 B CN 108536929B CN 201810243782 A CN201810243782 A CN 201810243782A CN 108536929 B CN108536929 B CN 108536929B
- Authority
- CN
- China
- Prior art keywords
- waveguide
- characteristic value
- solving
- modes
- arpack
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Optical Integrated Circuits (AREA)
Abstract
Description
技术领域
本发明属于计算电磁学技术领域。在微波传输系统仿真设计过程中,涉及一种波导结构色散特性数值求解的方法,具体为一种应用ARPACK求解波导结构色散特性的方法。
背景技术
波导是一种微波定向传输结构,主要应用于通信和微波能量传输等领域。在设计波导传输系统的过程中,需要求解波导的色散特性。实际工程应用中,一般以波导传播的主模和前面几个高次特征模式作为工作模式,即截止频率最小的几个模式。波导色散特性的求解方法有有限差分法、有限元方法、矩量法等。在有限差分法中,有一种称为二维四分量频域有限差分法,其具体实施步骤如下。
步骤1、设定求解波导色散特性的频率范围,在求解的频率范围内等间距选取N个频率,并求出N个频率所对应的自由空间波数k0,n,n=1,2,…,N。设定所要研究的模式为波导前M个模式。
步骤2、对波导横截面进行二维网格划分,索引编号(u,v)网格的电磁场分量位置及索引编号如图1所示。
步骤3、将求解频率范围内某一频率的归一化频域麦克斯韦方程组离散化,消去纵向电磁场,就得到四个横向场分量方程,组合所有网格的四个横向场分量方程,得到特征值问题
[A]·[X]=λ[X] (1)
其中,[A]表示由所有网格的归一化频域离散麦克斯韦方程电磁场分量系数组成的矩阵,[X]表示所有网格四个横向电磁场分量组成的列向量,λ为特征值,λ=β/k0,n,k0,n为给定频率的自由空间波数,β为传播常数。
步骤4、求解特征值问题[A]·[X]=λ[X],得到所有特征值,在所有特征值中选取截止频率最小的M个模式的特征值,求得对应的传播常数。
步骤5、对所研究频率范围内的其余频率进行步骤3、4,得到所有频率的传播常数,即得到了该波导前M个模式的色散特性曲线。
在上述方法的步骤4中,需要求解特征值问题[A]·[X]=λ[X]的全部特征值,这种处理方法求解效率很低,特别是当网格数目较多时,其求解所花费时间是实际工程应用所不能接受的。ARPACK是一款求解特征值问题的开源软件,可以求解部分特征值。ARPACK有六种策略去求解前几个特征值:求解实部最大的前几个特征值,求解实部最小的前几个特征值,求解虚部最大的前几个特征值,求解虚部最小的前几个特征值,求解幅值最大的前几个特征值,求解幅值最小的前几个特征值。波导前M个模式所对应的截止波数大小关系为
kc,1≤kc,2≤…≤kc,M (2)
当所设定的频率大于波导模式截止频率,波导模式所对应的特征值为实数;当所设定的频率小于波导模式截止频率,波导模式所对应的特征值为纯虚数。如果给定频率所对应的自由空间波数k0满足关系
kc,m≤k0≤kc,m+1 (3)
那么,波导前M个模式所对应的特征值满足
λ1≥λ2≥…≥λm≥0且0≤|λm+1|≤|λm+2|≤…≤|λM| (4)
λ1~λm为实数且大小递减,λm+1~λM为纯虚数且幅值递增,可以看出波导前M个模式所对应的特征值不满足ARPACK六种求解模式,所以ARPACK不适用于上述方法特征值问题的求解。
发明内容
针对上述存在问题或不足,为了能使用ARPACK高效求解波导结构的相移常数且适用于衰减常数求解,本发明提供了一种应用ARPACK求解波导结构色散特性的方法。
具体技术方案如下:
为了描述简单,设波导横截面位于x-y平面,波导沿z方向传输。波导内传输的电磁场表示为
步骤1、设定求解波导色散特性的频率范围,最小频率fmin,最大频率fmax。在求解的频率范围内等间距选取N个频率,即Δf=(fmax-fmin)/(N-1),求出N个频率所对应的自由空间波数k0,n=2π[fmin+Δf(n-1)]/c0,n=1,2,…,N,c0为光速。设定所要研究的模式为波导前M个模式。
步骤2、对波导横截面进行二维网格划分,并对网格及电磁场分量进行编号(如图1所示);
其中,k0,n表示给定频率所对应的自由空间波数,er表示电介质相对介电常数,dx、dy分别表示x、y方向网格长度,u、v分别为x、y方向索引编号。
组合所有网格的(6)至(9)四个横向场分量方程,得到特征值问题[A]·[X]=λ[X],其中,[A]表示由所有网格(6)至(9)式中等号右边电磁场分量系数所组成的矩阵,[X]表示所有网格四个横向电磁场分量组成的列向量,即[X]=[Ex|Ey|Hx|Hy]T,λ=-jγ,j为虚数单位。在等式两边同乘以λ,进而得到处理后的特征值问题
[A]·[A]·[X]=λ2[X] (10)
步骤4、应用ARPACK求解特征值问题前2*M个实部最大的特征值,然后得到前M个模式的传播常数。
步骤5、对所研究频率范围内的其余频率进行步骤3、4,得到所有频率的传播常数,即得到了该波导前M个模式的色散特性曲线。
本发明通过在特征值问题[A]·[X]=λ[X]等号两边同时乘以λ的方法,对特征值方程进行了等价改造,原特征值问题的特征值λ改为λ2,波导特征值问题[A]·[A]·[X]=λ2[X]所对应的前M个模式的特征值满足关系
综上所述,本发明通过特征值方程的改造,使得波导特征值问题可以采用ARPACK求解最大实部特征值的策略去求解,不仅适用于波导相移常数求解,而且适用于波导衰减常数的求解,并极大地提高了计算效率。
附图说明
图1为索引编号为(u,v)的二维网格电磁场分量分布及编号示意图;
图2为实施实例标准矩形波导横截面示意图;
图3为矩形波导相移常数曲线与理论值对比图;
图4为矩形波导衰减常数曲线与理论值对比图。
具体实施方式
下面以一个标准矩形波导为例,其宽边长度为a=20.86mm,窄边长度为b=10.16mm,如图2所示,通过具体实施过程对本发明作进一步详细说明。
步骤1、设定求解频率范围为fmin=1GHz,fmax=30GHz。在求解的频率范围内等间距选取N=30个频率,即Δf=1GHz,求出N个频率所对应的自由空间波数k0,n=2π[fmin+Δf(n-1)]/c0,n=1,2,…,30,c0为光速。
步骤2、对波导横截面进行二维网格划分。波导横截面x、y方向网格数目分别为nx=40、ny=20,网格大小为dx=0.5715mm、dy=0.508mm。
步骤3、给定求解频率范围内某一频率,组合所有网格的(6)至(9)频域离散麦克斯韦方程,得到特征值问题[A]·[X]=λ[X],在等式两边同乘以λ,进而得到处理后的特征值问题[A]·[A]·[X]=λ2[X]
步骤4、应用ARPACK求解特征值问题[A]·[A]·[X]=λ2[X]的前10个实部最大的特征值,然后得到前5个模式的传播常数。
步骤5、对所研究频率范围内的其余频率进行步骤3、4,得到所有频率的传播常数,即得到了该波导前5个模式的色散特性曲线。
表1
通过计算我们得到了图2所示波导的色散特性,即相移常数曲线和衰减常数曲线,如图3、4所示。从图3、4可以看出,求解得到的前5个模式(TE10、TE20、TE01、TE11、TM11)的色散曲线与理论值完全吻合。
表1展示的是求解特征值问题[A]·[X]=λ[X]全部特征值和采用ARPACK求解[A]·[A]·[X]=λ2[X]前10个实部最大的特征值花费时间的对比。从表中可以看出,采用ARPACK求解波导特征值问题是非常高效的,且求解全部特征值在网格数较多的情况下所花费的时间对于实际工程是无法接受的。
综上可见,本发明通过特征值方程的改造,使得波导特征值问题可以采用ARPACK求解最大实部特征值的策略去求解,不仅适用于波导相移常数求解,而且适用于波导衰减常数的求解,并极大地提高了计算效率。
Claims (1)
1.一种应用ARPACK求解波导结构色散特性的方法,具体步骤如下:
设波导横截面位于x-y平面,波导沿z方向传输,波导内传输的电磁场表示为
步骤1、设定求解波导色散特性的频率范围,最小频率fmin,最大频率fmax;在求解的频率范围内等间距选取N个频率,即Δf=(fmax-fmin)/(N-1),求出N个频率所对应的自由空间波数k0,n=2π[fmin+Δf(n-1)]/c0,n=1,2,…,N,c0为光速,设定所要研究的模式为波导前M个模式;
步骤2、对波导横截面进行二维网格划分,并对网格及电磁场分量进行编号;
其中,k0,n表示给定频率所对应的自由空间波数,er表示电介质相对介电常数,dx、dy分别表示x、y方向网格长度,u、v分别为x、y方向索引编号;
组合所有网格的(2)至(5)四个横向场分量方程,得到特征值问题[A]·[X]=λ[X],其中,[A]表示由所有网格(2)至(5)式中等号右边电磁场分量系数所组成的矩阵,[X]表示所有网格四个横向电磁场分量组成的列向量,即[X]=[Ex Ey Hx Hy]T,λ=-jγ,j为虚数单位;然后在等式两边同乘以λ,进而得到处理后的特征值问题
[A]·[A]·[X]=λ2[X] (6)
步骤4、应用ARPACK求解特征值问题前2*M个实部最大的特征值,然后得到前M个模式的传播常数;
步骤5、对所研究频率范围内的其余频率进行步骤3、4,得到所有频率的传播常数,即得到了该波导前M个模式的色散特性曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810243782.XA CN108536929B (zh) | 2018-03-23 | 2018-03-23 | 一种应用arpack求解波导结构色散特性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810243782.XA CN108536929B (zh) | 2018-03-23 | 2018-03-23 | 一种应用arpack求解波导结构色散特性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108536929A CN108536929A (zh) | 2018-09-14 |
CN108536929B true CN108536929B (zh) | 2021-09-24 |
Family
ID=63483618
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810243782.XA Active CN108536929B (zh) | 2018-03-23 | 2018-03-23 | 一种应用arpack求解波导结构色散特性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108536929B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117494533B (zh) * | 2024-01-02 | 2024-03-12 | 电子科技大学 | 一种涉及磁边界的四分量波导端口特性求解方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556623A (zh) * | 2008-04-09 | 2009-10-14 | 中国科学院微电子研究所 | 基于方波导的2×2单面双鳍线阵的制作方法 |
CN101728054A (zh) * | 2009-12-25 | 2010-06-09 | 中国科学院光电技术研究所 | 一种应用于定向辐射调制的表面人工电磁材料 |
CN102411676A (zh) * | 2011-10-31 | 2012-04-11 | 中国舰船研究设计中心 | 一种适用于不同形状波导方向图计算的面激励方法 |
CN102592057A (zh) * | 2012-01-17 | 2012-07-18 | 电子科技大学 | 周期结构指定频率的本征分析方法 |
CN103412989A (zh) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法 |
CN104335679A (zh) * | 2011-12-28 | 2015-02-04 | 王子控股株式会社 | 有机发光二极管、有机发光二极管的制造方法、图像显示装置及照明装置 |
CN104584326A (zh) * | 2012-05-09 | 2015-04-29 | 杜克大学 | 超材料设备及使用该超材料设备的方法 |
US10310186B2 (en) * | 2015-03-13 | 2019-06-04 | California Institute Of Technology | Method of engineering the dispersion of whispering gallery mode resonators and the resonators with dispersion engineered by the method |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1885757A (zh) * | 2006-07-05 | 2006-12-27 | 东南大学 | 多向传输接口高性能平面光波光路谐振腔滤波器 |
US9791618B2 (en) * | 2013-08-22 | 2017-10-17 | The University Of North Carolina At Charlotte | Waveguides incorporating novel metamaterials and associated methods of manufacture |
CN107255755B (zh) * | 2017-06-27 | 2019-09-20 | 华北电力大学 | 一种计算开孔矩形屏蔽体内部场分布的分区解析模型的方法 |
-
2018
- 2018-03-23 CN CN201810243782.XA patent/CN108536929B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101556623A (zh) * | 2008-04-09 | 2009-10-14 | 中国科学院微电子研究所 | 基于方波导的2×2单面双鳍线阵的制作方法 |
CN101728054A (zh) * | 2009-12-25 | 2010-06-09 | 中国科学院光电技术研究所 | 一种应用于定向辐射调制的表面人工电磁材料 |
CN102411676A (zh) * | 2011-10-31 | 2012-04-11 | 中国舰船研究设计中心 | 一种适用于不同形状波导方向图计算的面激励方法 |
CN104335679A (zh) * | 2011-12-28 | 2015-02-04 | 王子控股株式会社 | 有机发光二极管、有机发光二极管的制造方法、图像显示装置及照明装置 |
CN102592057A (zh) * | 2012-01-17 | 2012-07-18 | 电子科技大学 | 周期结构指定频率的本征分析方法 |
CN104584326A (zh) * | 2012-05-09 | 2015-04-29 | 杜克大学 | 超材料设备及使用该超材料设备的方法 |
CN103412989A (zh) * | 2013-08-01 | 2013-11-27 | 电子科技大学 | 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法 |
US10310186B2 (en) * | 2015-03-13 | 2019-06-04 | California Institute Of Technology | Method of engineering the dispersion of whispering gallery mode resonators and the resonators with dispersion engineered by the method |
Non-Patent Citations (1)
Title |
---|
电磁问题中有限元快速求解方法的研究;丁卫营;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20140715(第07期);C042-7 * |
Also Published As
Publication number | Publication date |
---|---|
CN108536929A (zh) | 2018-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Polstyanko et al. | Fast frequency sweep technique for the efficient analysis of dielectric waveguides | |
Ludick et al. | A tracking algorithm for the eigenvectors calculated with characteristic mode analysis | |
Gao et al. | A research on the electromagnetic properties of plasma photonic crystal based on the symplectic finite-difference time-domain method | |
CN108536929B (zh) | 一种应用arpack求解波导结构色散特性的方法 | |
Sadrpour et al. | A 3-D hybrid implicit explicit single-field finite-difference time-domain method | |
CN108090296B (zh) | 基于高阶辛紧致格式的波导全波分析方法 | |
Mishakin et al. | Analysis of dispersion and losses in helically corrugated metallic waveguides by 2-D vector finite-element method | |
Aghabarati et al. | Multilevel methods for $ p $-adaptive finite element analysis of electromagnetic scattering | |
Jiang et al. | Mixed finite element method for 2D vector Maxwell's eigenvalue problem in anisotropic media | |
Ahmed et al. | A time-domain adjoint variable method for materials with dispersive constitutive parameters | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
Zemlyakov | The Band-pass SIW-filter Based on L-ridged Rectangular Waveguide. | |
Brazhnyi et al. | Management of matter waves in optical lattices by means of the Feshbach resonance | |
Kaufmann | The meshless radial point interpolation method for electromagnetics | |
Tavallaee et al. | Finite-element modeling of evanescent modes in the stopband of periodic structures | |
Czarniewska et al. | Local mesh morphing technique for parametrized macromodels in the finite element method | |
Erden et al. | Excitation by a transient signal of the real-valued electromagnetic fields in a cavity | |
Hrizi et al. | Improving the convergence of the wave iterative method by filtering techniques | |
Shiraiwa et al. | Plasma wave simulation based on versatile FEM solver on Alcator C‐mod | |
Khalef et al. | An unconditionally stable radial point interpolation meshless method based on the Crank–Nicolson scheme solution of wave equation | |
Alexandru et al. | Efficient large electromagnetic problem solving by hybrid tlm and modal approach on grid computing | |
Hebermehl et al. | Improved numerical methods for the simulation of microwave circuits | |
Kulygin | Calculation of dispersion characteristics of a circular waveguide with deep helical corrugation by the FDTD method | |
Wright et al. | Coupling an ICRF core spectral solver to an edge FEM code | |
Gusev et al. | Symbolic-Numerical calculations of high-| m| rydberg states and decay rates in strong magnetic fields |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |