CN105717371B - 一种慢波结构慢波色散特性的获取方法 - Google Patents

一种慢波结构慢波色散特性的获取方法 Download PDF

Info

Publication number
CN105717371B
CN105717371B CN201610061581.9A CN201610061581A CN105717371B CN 105717371 B CN105717371 B CN 105717371B CN 201610061581 A CN201610061581 A CN 201610061581A CN 105717371 B CN105717371 B CN 105717371B
Authority
CN
China
Prior art keywords
slow
wave
electron beam
wave structure
beam group
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.)
Expired - Fee Related
Application number
CN201610061581.9A
Other languages
English (en)
Other versions
CN105717371A (zh
Inventor
王媛媛
魏彦玉
李大治
张鲁齐
丁冲
郭彍
宫玉彬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610061581.9A priority Critical patent/CN105717371B/zh
Publication of CN105717371A publication Critical patent/CN105717371A/zh
Application granted granted Critical
Publication of CN105717371B publication Critical patent/CN105717371B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16—Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Particle Accelerators (AREA)

Abstract

本发明公开了一种慢波结构慢波色散特性的获取方法，在慢波结构纵向的一端设置一个发射极，使之可以沿纵向发射一个速度固定的短电子束团，同时，在慢波结构的纵向中间位置设置一个观测点，用于记录短电子束团所激励起的电磁信号随时间的变化；在此基础上，改变所发射的短电子束团的速度，对不同电子束团速度下记录的时间变化的电磁信号进行频谱分析，得到不同短电子束团速度下所激励起的电磁波的频率，从而得到慢波色散特性。实验表明，本发明对于封闭的周期型慢波结构、开放的周期性慢波结构以及非周期性的慢波结构都具有良好的效果，获取的慢波色散特性与理论(实际)的基本一致，实现对开放的或非周期型的慢波结构色散特性的获取。

Description

一种慢波结构慢波色散特性的获取方法
技术领域
本发明属于电子器件技术领域，更为具体地讲，涉及一种慢波结构慢波色散特性的获取方法。
背景技术
O型真空电子器件如行波管和返波震荡管在通信、雷达以及空间应用等领域发挥着重大作用。这类器件主要利用慢波结构中的电磁慢波与直线型电子注在同步条件下得到互作用进而放大电磁场，作为注波互作用的场所，慢波结构是其的核心部件之一。所谓慢波是指相速小于光速的电磁波。在设计慢波结构的过程中，首要任务是要研究其电磁波的相速与频率的关系，即色散特性。
早期，主要采用实验法来研究慢波结构的慢波色散特性：实际制造一个慢波结构，将其两端短路，然后通过实验，测其谐振频率，进而得到慢波结构的慢波色散关系。然而，此方法费时费力费钱，尤其是对于目前作为研究热点的尺寸极小、极难加工的太赫兹波段的慢波结构。
随着计算机技术的发展，人们开始用数值计算和模拟的方法来研究慢波结构的慢波色散特性。常用的模拟方法有两种：谐振频率法和准周期法。其中，谐振频率法与实验法原理类似，由实验制造慢波结构后测试得到谐振频率，变为在模拟中建立虚拟模型，用模拟的方式获得谐振频率，这种方法只能用来研究封闭的慢波结构。准周期法是在模拟软件中如CST MWS(CST MICROWAVE STUDIO，是德国CST(Computer Simulation Technology)公司推出的高频三维电磁场仿真软件)、HFSS(High Frequency Structure Simulator)建立模型，利用准周期边界，在一个周期长度或半个周期长度上设置一个相移，相应的频率可以通过本征模式计算得到，从而得到相移和频率的对应关系。但是这种方法只能用于封闭的周期型慢波结构。
发明内容
本发明的目的在于克服现有技术的不足，提供一种慢波结构慢波色散特性的获取方法，以实现对开放的或非周期型的慢波结构色散特性的获取。
为实现上述发明目的，本发明慢波结构慢波色散特性的获取方法，其特征在于，包括以下步骤：
(1)、在慢波结构纵向的一端设置一个发射极，使之可以沿纵向发射一个速度固定的短电子束团；
(2)、在慢波结构的纵向中间位置设置一个观测点，用于记录短电子束团所激励起的电磁信号随时间的变化；
(3)、改变所发射的短电子束团的速度，对不同短电子束团速度下对记录的时间变化的电磁信号进行频谱分析，得到不同短电子束团速度下所激励起的电磁波的频率；
这样，得到的短电子束团速度和所激励起的电磁波的频率的对应关系为电磁波相速vp和频率f的关系(慢波色散特性)。
本发明的目的是这样实现的。
本发明慢波结构慢波色散特性的获取方法，在慢波结构纵向的一端设置一个发射极，使之可以沿纵向发射一个速度固定的短电子束团，同时，在慢波结构的纵向中间位置设置一个观测点，用于记录短电子束团所激励起的电磁信号随时间的变化；在此基础上，改变所发射的短电子束团的速度，对不同电子束团速度下记录的时间变化的电磁信号进行频谱分析，得到不同短电子束团速度下所激励起的电磁波的频率；得到的电子束团速度和所激励起的电磁波的频率的对应关系为电磁波相速vp和频率f的关系即慢波色散特性。
实验表明，本发明慢波结构慢波色散特性的获取方法对于封闭的周期型慢波结构、开放的周期性慢波结构以及非周期性的慢波结构都具有良好的效果，获取的慢波色散特性与理论(实际)的基本一致，实现对开放的或非周期型的慢波结构色散特性的获取。
附图说明
图1是本发明慢波结构慢波色散特性的获取方法中慢波结构为正弦波纹波导的结构示意图；
图2是图1所示正弦波纹波导的慢波色散特性曲线准周期法与本发明的比较图；
图3是本发明慢波结构慢波色散特性的获取方法中慢波结构为开放光栅的结构示意图；
图4是图3所示开放光栅的慢波色散特性曲线理论与本发明的比较图；
图5本发明慢波结构慢波色散特性的获取方法中慢波结构为介质波导的结构示意图；
图6是图5所示介质波导的慢波色散特性曲线理论与本发明的比较图。
具体实施方式
下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。
一、基本原理
当一个短电子束团(理想状态下近似为点电荷)沿着某个慢波结构的纵向匀速运动时，会在慢波结构中激励起电磁波，所激励的电磁波的相速度与短电子束团(点电荷)的速度相同。通过计算或模拟的方法，可以得到电子束团(点电荷)在慢波结构中激励的电磁波的频率，如此便得到频率和相速度的对应关系，这便是本发明电子束团激励法得到慢波结构色散的基本原理。
二、理论依据
首先研究点电荷在无界空间中匀速运动时所产生的电磁场。假设一个电量为q的点电荷以速度(为x方向单位矢量)在无界空间中匀速运动，相应的波动方程如下两式所示。
上式中E为电场强度，H为磁场强度，J为电流密度，ρ为电荷密度，ε和μ为介质的介电常数和磁导率。x，y和z是位置坐标，t为时间，
为了简化方程，对上述各个变量做傅里叶变换，如下式所示。
傅里叶逆变换为
上两式中表示E，H，J和ρ；表示和为各变量的傅里叶变换。β和γ分别是y方向和z方向的波数。i是虚数单位，ω是角频率。
对于一个匀速运动的点电荷，其电荷密度和电流密度矢量分别可以表示为
将上面两式带入到式(3)中，可以得到相应变量的傅里叶变换如下两式所示。
上式中α＝ω/ve。
对式(1)和(2)同样做傅里叶变换，并把式(5)和(6)带入其中。此处给出纵向x方向的场分量的波动方程：
上两式的解可以分别写为：
式中，是z方向的波数，是介质中的光速。Ai(i＝1，2，3)是待确定的系数。其他场分量可以由上两式根据纵向场法得到。将上两式带入到(4)中即可得到时域场的表达式。
由上两式可知，点电荷所激励起的电磁波在电荷运动方向上的相速为：
vp＝ω/α＝ve (13)
即，点电荷所激励起的电磁波在电荷运动方向上的相速等于点电荷的速度。当点电荷以某一速度匀速运动时，所激励起的场为幅值不同频率不同的电磁波的叠加，但是这些电磁波在纵向的相速与电子的速度相同。点电荷沿着慢波结构的纵向运动时，只有满足由慢波结构所构成的边界条件的某一频率的电磁波才能被激励起来，其纵向相速仍等于电子的速度。我们已知点电荷的速度，即电磁波在纵向的相速，通过模拟的方式则可以得到在慢波结构中所被激励起的电磁波的频率。由此可以得到电磁波频率和相速的关系，即慢波色散关系。
三、具体实施步骤(模型仿真)
(1)、建立一个纵向足够长的慢波结构模型，慢波结构的边界条件与实际结构相同；
(2)、在慢波结构纵向的一端设置一个发射极，使之可以沿纵向发射一个速度固定的短电子束团；
(3)、在慢波结构纵向中间的位置设置一个观测点；
(4)、求解模型，记录短电子束团所激励起的电磁信号随时间的变化；
(5)、对步骤(3)中所记录的电磁信号进行频谱分析，记录所激励起的频率；
(6)、改变所发射的短电子束团的速度，重复步骤4)和5)；
(7)、得到短电子束团速度和频率的对应关系，即为电磁波相速vp和频率f的关系(色散关系)。
对于电磁波的相速，我们有如下关系：
vp＝ω/β＝2πf/β (14)
式中β为波数。所以电磁波相速vp和频率f的关系也可以通过上式变为f和β的对应关系或者β和ω的对应关系，这些对应关系是色散关系的不同表现形式。
四、物理条件
1、点电荷是理想状态。在本发明中，用横向和纵向尺寸均远远小于慢波结构中电磁波波长的短电子束团来代替；
2、短电子束团所带的电荷量需尽可能少。当电荷量增大时，会增强与电磁场的互作用，短电子束团的速度变化会增大，会影响本发明的准确性。一般情况下，其电荷量量级小于等于10-15C；
3、要求所建立的慢波结构模型在纵向相对较长，以便减弱边缘效应的影响。在分析周期型慢波结构的慢波色散特性时，可用50个周期长度；
4、本发明不适用于计算快波色散特性。快波是指相速大于光速的波。由于所发射的短电子束团速度不能超过光速，所以本发明只能用于获取慢波结构的慢波色散特性。
五、实例
在实例中，所有实例均采用CST Particle Studio进行模拟。
1、封闭的周期型慢波结构：正弦波纹波导
在本实施例中，如图1所示，正弦波纹波导的半径可以由下式表述
R＝R0[1+h cos(2πx/L)] (15)
在公式(15)中R0为波纹波导的平均半径，h为归一化的波纹幅值，L为波导周期长度。此处所用的尺寸参数如下R0＝1.5cm，L＝1.67cm，h＝0.273。由准周期和电子束团激励法得到的色散关系如图2所示，此处将相速和频率的对应关系根据式(13)处理为f和β的对应关系，由于一个周期上的相移等于βL，所以进而处理成为相移βL和频率的对应关系。两种结果吻合的非常好。但是由于本发明即短电子束团激励法不能获得快波色散关系，所以快波区域没有电子束团激励法的结果。
2、开放的周期性慢波结构：开放光栅
如图3所示，在本实施例中，开放光栅的尺寸如下：周期L＝0.729mm，光栅槽的深度h＝0.447mm，光栅槽的宽度d＝0.533mm。
简化的开放光栅的色散方程如下
上式中βn＝β+2nπ/L是x方向的波数，是y方向的衰减常数，k＝ω/c是真空中波的传播常数。用电子束团激励法和理论所得到的色散曲线如下图所示，所画色散曲线仍是相移和频率的对应关系。如图4所示，由于理论计算与本发明模拟结果趋势相同，但有一定差距。这是由于理论是二维的简化理论，默认z方向无限大，而模拟时z方向不能取无限大。此外，在理论推导的过程中，对光栅槽里的场表达式做了近似处理，这也是区别产生的原因。
3、非周期性的慢波结构：介质波导
如图5所示，在本实施例中，平板介质波导尺寸：h0＝2mm，L＝100mm，w＝50mm。所用介质的参数为相对介电常数εr＝10，相对磁导率μr＝1。
介质波导的简化二维理论色散方程为：
εrα＝τ tan τh0 (17)
上式中，和分别为介质波导外和介质波导内的电磁场在x方向的传播常数。此处β为纵向波数，ω＝2πf为角频率。由理论得到和电子束团激励法得到的色散曲线如下图所示。由于结构的非周期性，不能计算其周期长度上的相移，所以此处我们选择了β和f的关系来表征色散关系。图6表明，理论和电子束团激励法的结果比较吻合，但是由于理论是二维理论，而模拟是三维的，所以两条色散曲线有一定偏离。
尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

Claims (2)

1.一种慢波结构慢波色散特性的获取方法，其特征在于，包括以下步骤：
(1)、在慢波结构纵向的一端设置一个发射极，使之可以沿纵向发射一个速度固定的短电子束团；
(2)、在慢波结构的纵向中间位置设置一个观测点，用于记录短电子束团所激励起的电磁信号随时间的变化；
(3)、改变所发射的短电子束团的速度，对不同短电子束团速度下对记录的时间变化的电磁信号进行频谱分析，得到不同短电子束团速度下所激励起的电磁波的频率；
这样，得到的短电子束团速度和所激励起的电磁波的频率的对应关系为电磁波相速vp和频率f的关系即慢波色散特性。
2.根据权利要求1所述的获取方法，其特征在于，所述的短电子束团电荷量量级小于等于10-15C。
CN201610061581.9A 2016-01-29 2016-01-29 一种慢波结构慢波色散特性的获取方法 Expired - Fee Related CN105717371B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610061581.9A CN105717371B (zh) 2016-01-29 2016-01-29 一种慢波结构慢波色散特性的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610061581.9A CN105717371B (zh) 2016-01-29 2016-01-29 一种慢波结构慢波色散特性的获取方法

Publications (2)

Publication Number Publication Date
CN105717371A CN105717371A (zh) 2016-06-29
CN105717371B true CN105717371B (zh) 2018-05-29

Family

ID=56155316

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610061581.9A Expired - Fee Related CN105717371B (zh) 2016-01-29 2016-01-29 一种慢波结构慢波色散特性的获取方法

Country Status (1)

Country Link
CN (1) CN105717371B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106680684B (zh) * 2017-03-24 2023-06-02 中国工程物理研究院应用电子学研究所 一种折叠波导行波管高频系统用测量装置及方法
CN107144739B (zh) * 2017-05-26 2019-10-18 电子科技大学 螺旋线慢波结构色散特性自动测试系统及方法
CN110909515B (zh) * 2019-10-31 2022-03-25 电子科技大学 一种获取慢波结构色散特性和耦合阻抗的方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011099600A1 (ja) * 2010-02-15 2011-08-18 国立大学法人九州大学 被験体状態解析用信号のピーク周波数測定システム
CN102255121A (zh) * 2011-05-11 2011-11-23 东南大学 基于圆柱线波导激励的宽带慢波系统
CN103414000A (zh) * 2013-08-26 2013-11-27 天津大学 局部分级深度材料光栅超慢太赫兹调制方法及波导
CN103412989A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6196793B2 (ja) * 2013-03-28 2017-09-13 アンリツ株式会社 ミリ波帯スペクトラム解析装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011099600A1 (ja) * 2010-02-15 2011-08-18 国立大学法人九州大学 被験体状態解析用信号のピーク周波数測定システム
CN102255121A (zh) * 2011-05-11 2011-11-23 东南大学 基于圆柱线波导激励的宽带慢波系统
CN103412989A (zh) * 2013-08-01 2013-11-27 电子科技大学 基于参数化降阶模型周期结构的三维电磁场仿真模拟方法
CN103414000A (zh) * 2013-08-26 2013-11-27 天津大学 局部分级深度材料光栅超慢太赫兹调制方法及波导

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
正弦波导及其应用的研究;许雄;《中国博士学位论文全文数据库信息科技辑》;20121215(第12期);第13页第2段、第17页第2段以及图2-7 *

Also Published As

Publication number Publication date
CN105717371A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
Jingzhu et al. Electromagnetic vibration noise analysis of transformer windings and core
Jiao et al. A fast higher-order time-domain finite element-boundary integral method for 3-D electromagnetic scattering analysis
CN105717371B (zh) 一种慢波结构慢波色散特性的获取方法
Staker et al. Alternating-direction implicit (ADI) formulation of the finite-difference time-domain (FDTD) method: Algorithm and material dispersion implementation
CN104123408B (zh) 基于有限元仿真分析的智能电表防外界电磁场干扰方法
Zhan et al. Stabilized DG-PSTD method with nonconformal meshes for electromagnetic waves
CN102385655B (zh) 一种电子设备辐射耦合电磁敏感性仿真方法
CN103218487B (zh) 旋转对称天线罩和抛物面天线一体化电磁散射仿真方法
Zheng et al. An efficient method to reduce the numerical dispersion in the ADI-FDTD
Gupta et al. Development of the virtual lab module for understanding the concepts of electric and magnetic field patterns in rectangular waveguides and cavities.
CN105054932A (zh) 一种基于磁纳米粒子一次谐波幅值的成像方法
CN103605633B (zh) 一种粗网格大时间步时域有限差分方法
CN104679920B (zh) 波导器件微波气体放电数值模拟方法
Kang et al. An efficient 2-D compact precise-integration time-domain method for longitudinally invariant waveguiding structures
CN105277927A (zh) 飞行器编队瞬态电磁特性时域阶数步进分析方法
Wu et al. Modeling of bandpass GPR problem by HIE procedure with enhanced absorption
CN103514143B (zh) 一种粗网格快速时域有限差分方法
CN103294871A (zh) 一种扬声器音圈电感非线性特性的数值仿真分析方法
Khamaru et al. Dynamics of a toroidal pure electron plasma using 3D PIC simulations
Munir Computational approach for resonant frequency calculation of coaxial cavity resonator using cylindrical coordinate system-based FDTD method
Wang et al. Stochastic FDTD for stochastic problems with axial symmetry
Daiffallah Simulation of f-mode propagation through a cluster of small identical magnetic flux tubes
Wu et al. Bandpass HIE-PML algorithm with improved absorption and efficiency for Lorentz medium
Mao et al. The CPML for hybrid implicit-explicit FDTD method based on auxiliary differential equation
He et al. Novel Scheme for GPU-Accelerated Finite-Difference Time-Domain Simulation of Electromagnetic Wave Interaction with Magnetic Plasma

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180529

Termination date: 20210129

CF01 Termination of patent right due to non-payment of annual fee