220 GHz共焦波導(dǎo)回旋振蕩管的理論研究與模擬仿真

李雷雷+寧智勇+王平+劉大剛+鄢陽

摘 要： 隨著毫米波波段回旋管的研究深入，傳統(tǒng)的封閉式圓波導(dǎo)諧振腔的弊端越發(fā)明顯。根據(jù)回旋管非線性理論，設(shè)計了一只采用開放式準(zhǔn)光諧振腔作為高頻結(jié)構(gòu)的回旋振蕩管，其工作頻率為220 GHz、工作電壓為60 kV、電流為3 A、橫縱速度比α為1.5、工作模式為HE06模。通過采用自主研發(fā)的三維粒子模擬軟件CHIPIC對其進行數(shù)值模擬研究，分析其工作特性，并進一步優(yōu)化參數(shù)。仿真結(jié)果表明：所設(shè)計的回旋管在磁場為8.57 T的條件下工作，可獲得36 kW的峰值功率輸出，輸出功率效率可達20%。

關(guān)鍵詞： 回旋振蕩管； 共焦波導(dǎo)； 三維粒子模擬； CHIPIC； 衍射損耗

中圖分類號： TN12?34； O441.4 文獻標(biāo)識碼： A 文章編號： 1004?373X（2014）22?0117?04

Theoretical research and simulation of 0.22 THz confocal gyrotron oscillator

LI Lei?lei， NING Zhi?yong， WANG Ping， LIU Da?gang， YAN Yang

（School of Physical Electronics， University of electronic Science and Technology of China， Chengdu 610054， China）

Abstract： With the research of gyrotron working in millimeter wave band， the disadvantages of traditional closed circular waveguide resonator reveals more clearly. According to the gyrotron nonlinear theory， a gyrotron oscillator was designed， which adopts an open type quasi?optical resonant cavity as its HF structure. Its work frequency is 220 GHz， voltage 60 kV， current 3 A， velocity ratio 1.5 and working mode HE06. The independently developed three?dimensional particle simulation software CHIPIC is used to carry out numerical simulation， analyze its operating characteristics and optimize its parameters. The simulation result shows that the gyrotron oscillator can obtain 36 kW peak power output under the condition of the magnetic field of 8.57 T. The efficiency of its output power can reach 20%.

Keywords： gyrotron oscillator； confocal waveguide； three?dimensional particle simulation； CHIPIC； diffraction loss

0 引 言

電子回旋脈塞作為一種重要的毫米波輻射源器件，在等離子體加熱、受控核聚變、等離子體診斷、毫米波雷達、通信及電子對抗、工業(yè)加工等方面有著十分重要的應(yīng)用前景，電子回旋脈塞器件主要有回旋振蕩管、回旋速調(diào)管、回旋行波管以及回旋返波管等類型，其中回旋振蕩管作為一種具有大功率輻射源器件，應(yīng)用潛力十分遠大。但當(dāng)器件的工作頻率逐漸升高時，為了解決尺寸共渡效應(yīng)以及功率容量等問題，往往需要使用高次模作為工作模式，這就引入了模式競爭問題。在毫米波頻段，由于器件的尺寸很小，模式競爭也就更加難以抑制，因此調(diào)諧性能嚴(yán)重下降，所以需研究新型結(jié)構(gòu)的適用于毫米波頻段的回管。近年來， 很多學(xué)者都在努力尋求解決的辦法。一個目前認(rèn)為有前途的研究方向是采用準(zhǔn)光諧腔。準(zhǔn)光腔與電子回旋脈塞不穩(wěn)定性相結(jié)合， 是解決這個難點的一個有用途徑。

普通圓波導(dǎo)諧振腔工作于高頻率時，模式競爭問題嚴(yán)重，模式密度很大，模式分割度較小，并且非常難采取有效的模式抑制措施。而準(zhǔn)光波導(dǎo)不存在上述問題。在準(zhǔn)光腔內(nèi)部，同樣存在無限個分立的振蕩模式。模式密度會隨頻率的增大而增大。由于這種準(zhǔn)光腔的側(cè)壁不存在，導(dǎo)致一部分模式的場從側(cè)壁輻射出去。如果把全部的側(cè)壁去掉，那么會有很多模式的場從側(cè)壁輻射出去，因而相當(dāng)一部分模式的場不存在，模式密度減小。留下的是能夠穩(wěn)定存在的場，這些場的模式不會因為側(cè)壁不存在而消失。模式密度大體上保持不變。從而達到有效的模式抑制效果。

美國MIT提出采用一種新型準(zhǔn)光結(jié)構(gòu)作為回旋管的高頻結(jié)構(gòu)，并采用該結(jié)構(gòu)成功研制出0.14 THz回旋單腔管與回旋行波管[1?2]，國內(nèi)對此結(jié)構(gòu)的研究仍處于起步階段，本文采用這種準(zhǔn)光結(jié)構(gòu)設(shè)計了一只工作頻率0.22 THz，工作模式為TE06模的準(zhǔn)光波導(dǎo)回旋振蕩管，并利用自主研發(fā)的三維粒子模擬軟件CHIPIC對設(shè)計進行了模擬分析。

1 回旋管原理

回旋振蕩管（Gyro?Backward Wave Oscillator，Gyro?BWO）是在外加磁場引導(dǎo)的回旋電子注與高頻結(jié)構(gòu)中電磁波的返波相互作用的一類電子回旋脈塞器件?；匦袷幑艿碾娮踊匦C振條件如下：

[ω-kzνz-sΩc?0] （1）

[Ωc=em01γ0B] （2）

式中：s是諧波數(shù)；kz是縱向波數(shù)；[ω]是波的角頻率；vz是電子的軸向速度；[γ0]=（1-β2）-1/2是相對論因子。

當(dāng)回旋電子注通過高頻結(jié)構(gòu)時，由于受到TE模Eφ場的作用，電子速度會發(fā)生變化，從而導(dǎo)致一部分電子處于加速狀態(tài)，一部分電子處于減速狀態(tài)，電子回旋頻率也會隨之變化。這種回旋頻率的變化進一步導(dǎo)致回旋軌道中相位的群聚。所以當(dāng)波的多普勒漂移頻率略大于電子回旋頻率Ω或其諧波頻率sΩ時，高頻場使得大部分電子處于減速區(qū)。經(jīng)過一段時間，回旋電子注在角向形成群聚，電子注交出能量，高頻場獲得能量。就形成了電子注與場的能量交換。

2 準(zhǔn)光腔結(jié)構(gòu)設(shè)計與分析

選取準(zhǔn)光波導(dǎo)結(jié)構(gòu)的本征模式作為工作模式的情況下，可以將其視為一個腔體結(jié)構(gòu)。在此條件下利用準(zhǔn)光學(xué)方法和衍射損耗的數(shù)值計算法，對本征模式的場分布進行求解，得出共焦波導(dǎo)結(jié)構(gòu)在本征頻率下類高斯模式的電磁場分布方程如下：[ExHz∝exp[-iky-12x2ω2（y）-i12kx2R（y）+i12arctanykω20]] （3）

式中：k為電磁波的波數(shù)；高斯電子注的束半徑公式：

[ω（y）=ω01+（ykω20）2] （4）

[ω0]為半徑y(tǒng)=0時的束半徑值；波前曲率半徑公式：

[R（y）=y1+kω20y2] （5）

采用CST對準(zhǔn)光腔結(jié)構(gòu)進行了3維模型的建立和冷腔仿真，仿真的準(zhǔn)光波導(dǎo)結(jié)構(gòu)的本征模式TE05，TE06模式的場分布圖如圖1所示，與理論公式（3）求解出的場分布完全一致。

圖1 CST中準(zhǔn)光腔模型電場分布圖

當(dāng)[y=L2]，準(zhǔn)光腔的曲率半徑與高斯束的相前曲率半徑相匹配：

[Rc=RL2=L21+2kw20L2] （6）

因此，可以得到：

[w20=L2k2Rc-LL] （7）

在鏡面的邊界條件為Ex=0，根據(jù)式（1）～式（5）可以推導(dǎo)出準(zhǔn)光腔橫向的諧振條件如下：

[ktL-arctanL2RcL-L2=nπ] （8）

由于設(shè)計的準(zhǔn)光腔結(jié)構(gòu)為更為特殊的共焦準(zhǔn)光結(jié)構(gòu)（[Rc=L]），所以式（8）可以簡化為：

[kt=（n+14）πL] （9）

因此，準(zhǔn)光腔HEmn的本征頻率可簡化為：

[f=c2L（n+14）] （10）

[K2=K2t+K2z] （11）

準(zhǔn)光波導(dǎo)結(jié)構(gòu)如圖2所示，這是一種開放式結(jié)構(gòu)。上下兩鏡面半徑等于鏡面間距（Rc=L）。整個腔體結(jié)構(gòu)選取較為簡單的三段式單腔結(jié)構(gòu)，其中L2為主要注波互作用區(qū)域，L1段半徑逐漸縮小，與電子槍區(qū)域相接， 波被截止，L3段半徑逐漸擴大，與輸出段相接，波由此輸出?；匦艿母鱾€尺寸根據(jù)反復(fù)優(yōu)化得到，其中Rc=4.25 mm，L1=9 mm，L2=L3=19 mm，傾角[θ1]=3°，[θ2]=1°。

圖2 準(zhǔn)光腔回旋管結(jié)構(gòu)

根據(jù)式（10）可以得出準(zhǔn)光腔中TE模式的色散曲線如圖3所示，可看出準(zhǔn)光波導(dǎo)里面的模式密度較稀疏。

圖3 準(zhǔn)光腔色散曲線圖

3 準(zhǔn)光腔結(jié)構(gòu)模式選擇性分析

準(zhǔn)光波導(dǎo)之所以具有模式選擇性，是由于它的橫向開敞結(jié)構(gòu)能使他產(chǎn)生衍射損耗，不同模式的衍射損耗不同。為了能夠?qū)崿F(xiàn)損耗的定量加載，從而更好地解決回旋管中模式競爭的問題，這里對準(zhǔn)光波導(dǎo)中衍射損耗進行定量的分析。

首先引入一個菲涅爾參數(shù)CF，用它來描述準(zhǔn)光波導(dǎo)結(jié)構(gòu)損耗特性的一個參數(shù)，

[CF=Kta2Rc] （12）

式中[a]為準(zhǔn)光波導(dǎo)鏡面的半寬度。由于衍射損耗的存在，準(zhǔn)光波導(dǎo)的橫向波速實際上是一個復(fù)數(shù)，即（[kt=ktr2+kti2] ） ，前面得到的[kt] 即為這里的[ktr]，使用[kti] 來對衍射損耗進行描述。Boyd和Gordon兩人對準(zhǔn)光波導(dǎo)衍射損耗進行了研究得到

[Λ=2ktiRckti=1Rclnπ2CF1R（1）0，m（CF，1）] （13）

式中[R（1）0，m（CF，1）]為橢球坐標(biāo)系中的徑向波函數(shù)。考慮到用徑向波函數(shù)計算量較大，Baker等人對數(shù)據(jù)進行了擬合，對HEmn模將準(zhǔn)光波導(dǎo)用來作為模式選擇結(jié)構(gòu)設(shè)計的關(guān)鍵點在于利用衍射損耗對可能的寄生振蕩進行抑制，從圖可以看出，對于HEmn?？芍?，在鏡面寬度X一定的情況下，對于m>0的模式衍射損耗比m=0的模式要大得多，所以應(yīng)選取m=0的模式為工作模式。對于m=0時，n越大，衍射損耗越小，但是，n越大也會造成模式競爭的情況更加復(fù)雜（競爭模式為比工作模式低的模式），模式與電子注的耦合也更弱，所以選取工作模式應(yīng)該更加綜合的考慮。由于衍射損耗的存在，這使準(zhǔn)光波導(dǎo)即使不采用吸收體也具有很好的分布衰減特性，并且這種衰減對不同的模式是不用的，這就使得準(zhǔn)光波導(dǎo)具有良好的模式選擇性。另一方面，回旋脈塞器件的實用化需要輸出模式為高斯束，對于準(zhǔn)光波導(dǎo)來說，其中的模式場分布在橫向上已經(jīng)是高斯分布，這時準(zhǔn)光波導(dǎo)轉(zhuǎn)換器的設(shè)計被大大簡化。

圖4 模式與衍射損耗的關(guān)系

4 CHIPIC模擬結(jié)果分析

CHIPIC是由電子科技大學(xué)自主研發(fā)是一款三維帶電粒子與電磁場互作用的粒子模擬軟件，可對微波電真空器件進行粒子模擬仿真。為更準(zhǔn)確的研究所設(shè)計的回旋管的模式競爭問題，采用CHIPIC三維模擬軟件進行仿真研究，采用其中的回旋發(fā)射模型等效代替磁注入式電子槍。由于準(zhǔn)光腔器件是開放式結(jié)構(gòu)，本文在CHIPIC模擬時處理為在橫切面方向設(shè)置4個吸收邊界（即開放邊界）。仿真電參數(shù)根據(jù)回旋管非線性理論計算取電壓V=60 kV，電流=3 A，引導(dǎo)半徑1.1 mm，速度橫縱比a=1.5進行模擬。根據(jù)式（1）～式（3）可以算出。當(dāng)引導(dǎo)磁場B=7.3與B=8.6左右時，由準(zhǔn)光腔色散曲線可以看出，HE05與HE06模式能夠起振。圖5、圖6為CHIPIC仿真得到的HE05和模式HE06模式等位圖。

圖5 HE05模式的等位圖

從圖5中可以看出HE05模式的波瓣寬度要HE06模式的波瓣寬度大許多，這也能說明，在相同的準(zhǔn)光波導(dǎo)橫向?qū)挾葪l件下，HE05模式的功率更容易從兩邊策側(cè)壁衍射出去，從而得到較小的輸出功率。通過優(yōu)化引導(dǎo)磁場大小B，當(dāng)選取B=7.28與B=8.57時，回旋管工作在HE05、HE06模式下，它的輸出功率分別達到最大。如圖7、8所示，HE05和HE06模式的輸出功率分別為18 kW和35.6 kW，輸出頻率分別為185 GHz和220.1 GHz。圖中可以看出HE06模式的輸出功率要比HE05模式的輸出功率大上2倍。這也證明了上述所說的正確性，此時HE06工作模式的效率約為20%。

圖6 HE06模式的等位圖

圖7 HE05模式的頻率與輸出功率圖

圖8 HE06模式的頻率與輸出功率圖

本文針對大氣窗口頻率220 GHz，以及模擬得到的參數(shù)，選擇了HE06的工作模式進行以后的實驗。

5 結(jié) 論

本文設(shè)計了一只中心工作頻率在0.22 THz的準(zhǔn)光腔振蕩管，根據(jù)回旋管非線性理論選取了工作參數(shù)，并采用自主研發(fā)的三維粒子模擬軟件CHIPIC對其進行了模擬分析。模擬結(jié)果表明：在B=8.57，I=3 A，a=1.5，Lz=10 mm，V=60 kV時，可獲得36 kW的峰值功率輸出和20%的效率。本文為準(zhǔn)光腔回旋器件的實驗與研發(fā)奠定了一定的基礎(chǔ)。

參考文獻

[1] JOYE C D， SHAPIRO M A. Demonstration of a 140 GHz 1 kW confocal gyro?traveling?wave amplifier [J]. IEEE Transactions on Electron Devices， 2009， 56（5）： 818?827.

[2] HU Wen， SHAPIRO M A， KREISCHER K E. 140?GHz gyrotron experiments based on a confocal cavity [J]. IEEE Transactions on Plasma Sci，1998， 26（3）： 366?374.

[3] HONG K D， BRAND G F， IDEHARA T. A 150?600 GHz step?tunable gyrotron [J]. Journal of Applied Phys， 1993， 74（8）： 5250?5258.

[4] CHU K R. The electron cycltron maser [J]. Rev Modern Phy， 2004， 76（2）： 489?540.

[5] 劉盛綱.準(zhǔn)光諧振腔電子回旋脈塞的動力學(xué)理論[J].電子學(xué)報，1984，12（1）：13?27.

[6] BASTEN M A， GUSS W C， KREISCHER K E， et al. Experimental investigation of a 140 GHz gyrotron?backward wave oscillator [J]. International Journal of Infrared Millimeter Waves， 1995，16（5）：880?905.

[7] 鄢揚，劉盛綱，黎曉云，等.0.22 THz回旋管振蕩器的研制與實驗[J].科學(xué)通報，2009，54（4）：522?526.

[8] NUSINOVICH G S， DUMBRAJS O. Theory of gyro?backward wave oscillators with tapered magnetic field and waveguide cross section [J]. IEEE Transactions on Plasma Sci， 1996， 24（3） ： 620?629.

[9] 王文祥.微波工程技術(shù)[M].北京：國防工業(yè)出版社，2009.

[10] 王麗.回旋管電子光學(xué)系統(tǒng)的研究[D].成都：電子科技大學(xué)，2006.

圖5 HE05模式的等位圖

從圖5中可以看出HE05模式的波瓣寬度要HE06模式的波瓣寬度大許多，這也能說明，在相同的準(zhǔn)光波導(dǎo)橫向?qū)挾葪l件下，HE05模式的功率更容易從兩邊策側(cè)壁衍射出去，從而得到較小的輸出功率。通過優(yōu)化引導(dǎo)磁場大小B，當(dāng)選取B=7.28與B=8.57時，回旋管工作在HE05、HE06模式下，它的輸出功率分別達到最大。如圖7、8所示，HE05和HE06模式的輸出功率分別為18 kW和35.6 kW，輸出頻率分別為185 GHz和220.1 GHz。圖中可以看出HE06模式的輸出功率要比HE05模式的輸出功率大上2倍。這也證明了上述所說的正確性，此時HE06工作模式的效率約為20%。

圖6 HE06模式的等位圖

圖7 HE05模式的頻率與輸出功率圖

圖8 HE06模式的頻率與輸出功率圖

本文針對大氣窗口頻率220 GHz，以及模擬得到的參數(shù)，選擇了HE06的工作模式進行以后的實驗。

5 結(jié) 論

本文設(shè)計了一只中心工作頻率在0.22 THz的準(zhǔn)光腔振蕩管，根據(jù)回旋管非線性理論選取了工作參數(shù)，并采用自主研發(fā)的三維粒子模擬軟件CHIPIC對其進行了模擬分析。模擬結(jié)果表明：在B=8.57，I=3 A，a=1.5，Lz=10 mm，V=60 kV時，可獲得36 kW的峰值功率輸出和20%的效率。本文為準(zhǔn)光腔回旋器件的實驗與研發(fā)奠定了一定的基礎(chǔ)。

參考文獻

[1] JOYE C D， SHAPIRO M A. Demonstration of a 140 GHz 1 kW confocal gyro?traveling?wave amplifier [J]. IEEE Transactions on Electron Devices， 2009， 56（5）： 818?827.

[2] HU Wen， SHAPIRO M A， KREISCHER K E. 140?GHz gyrotron experiments based on a confocal cavity [J]. IEEE Transactions on Plasma Sci，1998， 26（3）： 366?374.

[3] HONG K D， BRAND G F， IDEHARA T. A 150?600 GHz step?tunable gyrotron [J]. Journal of Applied Phys， 1993， 74（8）： 5250?5258.

[4] CHU K R. The electron cycltron maser [J]. Rev Modern Phy， 2004， 76（2）： 489?540.

[5] 劉盛綱.準(zhǔn)光諧振腔電子回旋脈塞的動力學(xué)理論[J].電子學(xué)報，1984，12（1）：13?27.

[6] BASTEN M A， GUSS W C， KREISCHER K E， et al. Experimental investigation of a 140 GHz gyrotron?backward wave oscillator [J]. International Journal of Infrared Millimeter Waves， 1995，16（5）：880?905.

[7] 鄢揚，劉盛綱，黎曉云，等.0.22 THz回旋管振蕩器的研制與實驗[J].科學(xué)通報，2009，54（4）：522?526.

[8] NUSINOVICH G S， DUMBRAJS O. Theory of gyro?backward wave oscillators with tapered magnetic field and waveguide cross section [J]. IEEE Transactions on Plasma Sci， 1996， 24（3） ： 620?629.

[9] 王文祥.微波工程技術(shù)[M].北京：國防工業(yè)出版社，2009.

[10] 王麗.回旋管電子光學(xué)系統(tǒng)的研究[D].成都：電子科技大學(xué)，2006.

圖5 HE05模式的等位圖

從圖5中可以看出HE05模式的波瓣寬度要HE06模式的波瓣寬度大許多，這也能說明，在相同的準(zhǔn)光波導(dǎo)橫向?qū)挾葪l件下，HE05模式的功率更容易從兩邊策側(cè)壁衍射出去，從而得到較小的輸出功率。通過優(yōu)化引導(dǎo)磁場大小B，當(dāng)選取B=7.28與B=8.57時，回旋管工作在HE05、HE06模式下，它的輸出功率分別達到最大。如圖7、8所示，HE05和HE06模式的輸出功率分別為18 kW和35.6 kW，輸出頻率分別為185 GHz和220.1 GHz。圖中可以看出HE06模式的輸出功率要比HE05模式的輸出功率大上2倍。這也證明了上述所說的正確性，此時HE06工作模式的效率約為20%。

圖6 HE06模式的等位圖

圖7 HE05模式的頻率與輸出功率圖

圖8 HE06模式的頻率與輸出功率圖

本文針對大氣窗口頻率220 GHz，以及模擬得到的參數(shù)，選擇了HE06的工作模式進行以后的實驗。

5 結(jié) 論

本文設(shè)計了一只中心工作頻率在0.22 THz的準(zhǔn)光腔振蕩管，根據(jù)回旋管非線性理論選取了工作參數(shù)，并采用自主研發(fā)的三維粒子模擬軟件CHIPIC對其進行了模擬分析。模擬結(jié)果表明：在B=8.57，I=3 A，a=1.5，Lz=10 mm，V=60 kV時，可獲得36 kW的峰值功率輸出和20%的效率。本文為準(zhǔn)光腔回旋器件的實驗與研發(fā)奠定了一定的基礎(chǔ)。

參考文獻

[1] JOYE C D， SHAPIRO M A. Demonstration of a 140 GHz 1 kW confocal gyro?traveling?wave amplifier [J]. IEEE Transactions on Electron Devices， 2009， 56（5）： 818?827.

[2] HU Wen， SHAPIRO M A， KREISCHER K E. 140?GHz gyrotron experiments based on a confocal cavity [J]. IEEE Transactions on Plasma Sci，1998， 26（3）： 366?374.

[3] HONG K D， BRAND G F， IDEHARA T. A 150?600 GHz step?tunable gyrotron [J]. Journal of Applied Phys， 1993， 74（8）： 5250?5258.

[4] CHU K R. The electron cycltron maser [J]. Rev Modern Phy， 2004， 76（2）： 489?540.

[5] 劉盛綱.準(zhǔn)光諧振腔電子回旋脈塞的動力學(xué)理論[J].電子學(xué)報，1984，12（1）：13?27.

[6] BASTEN M A， GUSS W C， KREISCHER K E， et al. Experimental investigation of a 140 GHz gyrotron?backward wave oscillator [J]. International Journal of Infrared Millimeter Waves， 1995，16（5）：880?905.

[7] 鄢揚，劉盛綱，黎曉云，等.0.22 THz回旋管振蕩器的研制與實驗[J].科學(xué)通報，2009，54（4）：522?526.

[8] NUSINOVICH G S， DUMBRAJS O. Theory of gyro?backward wave oscillators with tapered magnetic field and waveguide cross section [J]. IEEE Transactions on Plasma Sci， 1996， 24（3） ： 620?629.

[9] 王文祥.微波工程技術(shù)[M].北京：國防工業(yè)出版社，2009.

[10] 王麗.回旋管電子光學(xué)系統(tǒng)的研究[D].成都：電子科技大學(xué)，2006.