新聞中心

EEPW首頁 > 模擬技術 > 設計應用 > 低能電子對單面介質加載微波部件微放電的影響*

低能電子對單面介質加載微波部件微放電的影響*

作者:倪鑫榮(湘潭大學 自動化與電子信息學院,湖南 湘潭 411105)時間:2021-11-22來源:電子產品世界收藏
編者按:在微波設備微放電效應的研究中,低能電子的作用通常被忽略。本文重點研究了低能電子對單面介質加載微波部件微放電的影響,基于蒙特卡洛算法,建立了單面介質加載微波部件的微放電的模型,模擬了單面介質加載微波部件微放電的規律,利用數值的方法對波導中靜電場進行求解并對電子的軌跡進行跟蹤。通過對比經典的Vaughan二次電子發射模型(考慮低能電子)和Rice模型(沒有考慮低能電子)的模擬結果,獲得了低能電子對單面介質加載微波部件微放電的影響。

*基金項目:湖南省自然科學基金(2017JJ3314)

本文引用地址:http://www.recipesarchive.com/article/202111/429820.htm

0   引言

效應(也被稱為二次電子倍增效應)是在真空條件下運行在微波和毫米波子系統中,在微波設備和微波元件中經常能觀察到的一種高功率共振電子放電,存在于各種不同的高功率微波場景中,如衛星通信有效載荷的無源部件、行波管、高功率射頻窗口以及粒子加速器等[1]

在許多現代射頻系統中,二次電子倍增效應是在真空或接近真空條件下工作的一個重要失效機制[2]。這種現象表現為自由電子在真空器件內的雪崩式增長,這是在射頻電場的作用下,被加速的高能電子撞擊器件表面時,從器件表面激發出新的二次電子發射引起的。在不同的射頻系統中,如星載通信、射頻加速器和高功率微波發生器,二次電子倍增擊穿問題變得越來越嚴重,這是為了滿足高性能的需求,微波功率不斷增加,并且盡可能減小制造微波裝置的尺寸所導致的。由于二次電子倍增放電可以顯著限制甚至損壞射頻系統,所以研究射頻系統中的二次電子倍增放電,對減少射頻系統的損壞有至關重要的作用[3-6]。這項工作將對快速發展的二次電子倍增效應的研究做出重大貢獻,將使研究人員能夠更好地設計射頻系統,同時,也使得二次電子倍增效應成為了空間研究領域重點關注的熱點課題之一[7]。在過去幾十年中,存在于金屬表面和介質表面之間的二次電子倍增效應已經引起了極大的關注。介質表面的二次電子倍增與金屬表面不同,在射頻電場的作用下,二次電子不會在金屬表面積累電荷,但會在介質表面積累相應的正電荷或負電荷[8],這些積累在介質表面的電荷會形成一個靜電場,影響電子的運動[9]

本文將通過對比經典的Vaughan 模型和Rice 模型在單面微波部件中二次電子倍增效應的差異,以獲得對單面微波部件中的影響。

1   模型介紹

圖1 展示了在下表面加載了介質的金屬平行板波導的示意圖,其中金屬平行板的長度和間距分別為a 和d。加載的介質層厚度為h,介質表面和金屬表面之間的間隙為H = d-h。在模型中設置了余弦變化的射頻電場(其中,E0 的取值為2×106 V/m),工作頻率為f(對應不同的諧振階數時,f 的取值有所不同),其對應的仿真時間步長設置為RF 周期的1/1 000,在每個射頻周期內電子運動1 000 個時間步,以保證仿真的精確度。設置初始電子的數目為N0 = 1 000,在前500 個時間步內,每個時間步從介質表面發射出2 個初始電子。每個初始電子的權重都設置為1,為了更加高效地進行計算,將電子數目的閾值上限設置為Nth = 1×106。當電子的數目超過Nth 時,會在下一步的模擬中將電子的數目隨機減少一半,同時使剩余電子的電荷和質量增加1 倍。在相反的情況下,當電子數目少于Nth 的10% 時,會在下一步的模擬中將每個電子分成2 個電子,并將電子的電荷和質量減半[10]。以這樣的方法進行計算,既保證了計算可行性,同時又保證了參與計算的電子數目。當激發出新的二次電子時,等效的正電荷會在介質表面積聚,形成非均勻分布的靜電場。隨著正電荷的不斷積聚,靜電場的強度亦隨之增強,雙面倍增可能演變為單面倍增并自我維持或熄滅,這容易受到許多其他因素的影響。當單面倍增發生時,電子云將在每個整數射頻周期內沖擊介質表面。

設置材料的最大二次電子產額系數為δmax0 = 2.36,與之對應的入射能量為Wmax0 = 300 eV[11],兩種模型的SEY 隨入射能量的變化如圖2 所示。從圖中可知,經典的Vaughan 模型將能量低于12.5 eV 的電子視為無效電子,沒有考慮其與金屬或介質表面的碰撞,而Rice模型則將其考慮在內并視為彈性散射電子[12]。經典的Vaughan 二次電子發射模型如下所示:

1637565788111589.png

Rice 模型對經典的Vaughan 模型進行了修改,其模型如下所示:

1637565845162364.png

其中,v = (Wim-Wth) / (Wmax-Wth),參數δmax、Wmax 和Wth 由材料和表面處理所決定。

初始電子從介質表面被隨機發射出去,初始電子的發射能量W0 滿足Maxwell-Boltzmann 分布,如式(3)所示:

1637565931273551.png

初始電子的發射角度服從正弦分布,如式(4)所示:

1637565990679154.png

1637566010796938.png

圖1 單面平行板波導模型

1637566036216504.png

圖2 SEY曲線

2   靜電場的計算

當電子與介質表面發生碰撞時,會在介質表面積聚電荷,然而電子與金屬表面碰撞則不會在金屬表面積聚電荷,兩者之間存在著差異。當介質表面的電荷積聚到一定數量時,就會形成一個非均勻分布的靜電場來影響電子的運動。所以,模型中靜電場的精確求解,對后續電子的運動具有重要的影響。本文采用數值的方法對介質表面由電荷積聚所產生的靜電場進行求解,設置電荷位于介質表面的r(x,h,0)處,靜電場的求解如式(5)所示:

1637566136687963.png

為了求解模型中靜電場的分布,我們應該先對表面電荷積聚所產生的電勢求解,利用疊加原理,通過增加每個電荷各自的貢獻,就可以得到介質表面的電荷在波導中產生的電勢分布,其表達式如下所示:

1637566198354144.png

其中,G 為單位電荷形成電勢的格林函數,求解出電勢的格林函數就能求解得到電勢。G的展開如式(7)所示:

1637566243404706.png

為了方便計算,在本文中假設點電荷沉積的位置位于介質表面和真空區域的交界處來進行計算。所考慮問題的幾何特征和線性性質使得狄拉克函數的展開可以如式(8)(9)所示:

1637566293178394.png

其中,kxn = nπ/a,kz 是沿著z 方向的傅里葉變換。

1637566359175054.png

其中,image.png是頻域中的格林函數。結合式(7)就可以求解出image.png為:

1637566419351534.png

image.png代入到式(5)中,就能求解出介質表面的電荷在模型中的電勢分布,從而求解出模型中靜電場的空間分布。

為了驗證所求出的靜電場在模型中的分布是否正確,下面在介質表面積聚固定數量的正電荷來進行模擬。令介質表面每個位置都均勻分布了1×107 個正電荷,在模型中產生的靜電場分布如圖3 所示,圖中箭頭的方向代表靜電場的方向,箭頭的長度代表靜電場的強度。從圖中可以看出,靜電場沿著介質表面呈現均勻分布,隨著與介質表面之間距離的增加,靜電場的強度呈現逐漸減弱的趨勢。取x = 5 mm 處的靜電場進行分析,靜電場的強度隨介質表面距離的變化,如圖4 所示,此時,靜電場最強的位置為介質表面達到了Edc = 1.37×105,在金屬表面的靜電場的強度最弱。模擬的結果和理論的結果保持一致,驗證我們理論的可行性。

1637566489232463.png

圖3 靜電場在模型中的分布

1637566516842853.png

圖4 靜電場的強度隨介質表面距離的變化

3   電子的動力學方程

我們通過忽略RF 磁場簡化了模型,可以通過在每個時間步中數值求解洛倫茲方程來計算電子的軌跡。在射頻電場ERF 的作用下,電子在金屬與介質表面之間周期性地運動,當與介質表面碰撞時,在激發新二次電子的同時會積聚相應的電荷,形成靜電場從而影響電子的運動。模型中電子的動力學方程如式(12)~(19)所示:

1637566616511583.png

設置初始條件,初始時刻t=0 時的速度為:

1637566669546021.png

沿著x,y 方向的速度分別為:

1637566712522078.png

沿著x,y 方向的位移分別為:

1637566748648568.png

4   實驗結果分析

平行板波導模型的長度為a = 10 mm,間隙距離為d = 1 mm,介質層的厚度為h = 0.05 mm。因此,金屬與介質表面之間的間距為H = 0.95 mm。設置RF 電壓VRF = 1.9 kV,工作頻率f = 6 GHz,在該參數下為一階的諧振階數。兩種模型下,電子數目隨時間的變化如圖5 所示,通過對比電子數目Ne 隨時間的變化可知,Rice 模型的二次電子倍增速率要大于經典的Vaughan模型, 并且在達到短暫的飽和時電子數目達到了Ne_R = 4×108,比經典的Vaughan 模型達到飽和時的電子數目Ne_V = 6×106 要更多, 兩種模型在達到短暫飽和之后的趨勢也有所不同,Rice 模型在達到飽和之后,Ne_R 仍然以極小的速度繼續增長,但是經典的Vaughan 模型在達到飽和之后,Ne_V 呈現下降的趨勢。與之對應的靜電場Edc 隨時間的變化如圖6 所示,兩種模型在達到飽和之后Edc 都達到了飽和并保持穩定,但Rice 模型達到飽和時的靜電場為Edc_R = 4×105,而經典的Vaughan 模型在達到飽和時的靜電場為Edc_V = 2.3×105。圖7 展示了兩種模型在電子數目達到飽和之后,電子相的空間分布。從圖7(a)中可以看出,在達到飽和時,經典的Vaughan 模型的電子在波導中呈現均勻分布的狀態,而從圖7(b)中可知,在達到飽和時,Rice 模型的電子主要分布在波導的介質表面附近,隨著與介質表面距離的增加,電子數目呈現減少的趨勢。圖8 和圖9 分別展示了兩個模型下平均SEY 和平均電子能量隨時間的變化,從圖8(a)可知,經典的Vaughan 模型在達到飽和之后平均SEY<1,入射電子大于激發的二次電子,電子數目呈下降趨勢,與之對應平均電子能量如圖9(a)所示,穩定在Wave = 21 eV,Rice 模型在達到飽和之后平均SEY > 1,入射電子大于激發的二次電子,電子數目呈現緩慢的增長,如圖8(b)所示,平均電子能量穩定在Wave = 15 eV,如圖9(b)所示。

1637566835195631.png

圖5 電子數目隨時間的變化

1637566866237395.png

圖6 靜電場隨時間的變化

image.png

圖7 電子的相空間分布,(a)為Vaughan模型,(b)為Rice模型

1637566938505134.png

圖8 平均SEY隨時間的變化,(a)為Vaughan模型,(b)為Rice模型

1637566976910993.png

圖9 平均電子能量隨時間的變化,(a)為Vaughan模型,(b)為Rice模型

4   結束語

本文分別對經典的Vaughan 模型和Rice 模型在單面介質加載微波部件中的效應進行了模擬仿真。研究結果表明,Rice 模型在電子的增長速率、達到飽和時電子的數目和靜電場的強度都要高于經典的Vaughan模型,考慮的二次電子發射模型更適合用于研究微波設備的微放電效應。

參考文獻:

[1] 崔萬照,李韻,張洪太,等.航天器微波部件微放電分析及其應用[M].北京:人民郵電出版社,2019.

[2] 樂群,殷希梅,黃學翔,等.微波部件微放電效應研究[C].電工理論學術年會,2007.

[3] VAGUE J, MELGAREJO J C,GUGLIELMI M,et al.Multipactor effect characterization of dielectric materials for space applications[j].ieee transactions on microwave theory and techniques,2018(66):3644-3655.

[4] POPOVIC S,UPADHYAY J,MAMMOSSER J,et al.Resonantfrequency discharge in a multi-cell radio frequency cavity[J/OL].Journal of Applied Physics,2014,116(11).173301.https://www.researchgate.net/publication/268448052_Resonant-frequency_discharge_in_a_multi-cell_radio_frequency_cavity

[5] COVES A,TORREGROSA-PENALVA G,VICENTE C,et al.Multipactor discharges in parallel-plate dielectric-loaded waveguides including space-charge effects[J].Electron Devices,IEEE Transactions on,2008,55(9):2505-2511.

[6] TORREGROSA-PENALVA G,COVES A,GIMENO MARTINEZ B,et al.Multipactor susceptibility charts of a parallel-plate dielectric-loaded waveguide[J].IEEE Transactions on Electron Devices,2010,57(5):1160-1166.

[7] 王新波,白鶴,李韻,等.星載微波部件介質微放電理論研究現狀及發展趨勢[J].空間電子技術,2019(6):1-9.

[8] TORREGROSA G,COVES N,VICENTE C P,et al.Time evolution of an electron discharge in a parallel-plate dielectric-loaded waveguide[J].IEEE Electron Device Letters,2006,27(7):619-621.

[9] BERENGUER,ANDRES,MESA,e t a l.A n a lysis o f mu l t i p a c t o r e f f e c t i n a P a r t i a l l y Di e l e c t r i c - L o a d e d Rectangular Waveguide[J].IEEE Transactions on Plasma Science,2019.

[10] 張雪,王勇,范俊杰,等.金屬壁與介質窗之間次級電子倍增效應的研究[J].物理學報,2014,63(16):397-405.

[11] LI Y,WANG D,YU M,et al.Experimental verification of multipactor discharge dynamics between ferrite dielectric and metal[J].IEEE Transactions on Electron Devices,2018(8):1-8.

[12] 張雪,王滔,倪鑫榮,等.低能二次電子對微波輸能窗擊穿現象的影響[J].強激光與粒子束,2020,32(10).49-54.

(本文來源于《電子產品世界》雜志2021年10月期)



評論


技術專區

關閉
中文天堂最新版在线www-bt天堂网www天堂-电影天堂 欧美在线香蕉在线视频
#耿直真香哥黑化卖惨# 我要我们在一起 外交部回应拜登重申不支持台独 男子写80页PPT拯救爱情却离婚 动物世界 中国共产党第三个历史决议全文发布 林丹世界排名被正式移除 国足最新出线概率0.08% #耿直真香哥黑化卖惨# 国际人士热议中共十九届六中全会 #耿直真香哥黑化卖惨# 周冠宇成为中国首位F1车手 大连一密接者擅自点外卖聚餐被调查 胡锡进谈中美元首会晤 动保组织向上饶信州区申请信息公开 扫黑风暴 北京冬奥火炬宣传片获金花环奖 意大利错失直接晋级世界杯资格 寻梦环游记 我要我们在一起 五个扑水的少年 男子写80页PPT拯救爱情却离婚 许家印为恒大注入超70亿续命资金 中美元首会谈重点内容 胡锡进谈中美元首会晤 大连一密接者擅自点外卖聚餐被调查 中美元首是否达成新共识?中方回应 国足最新出线概率0.08% 入殓师 得知母亲出事男子在地铁痛哭 红色通缉令 24岁救人牺牲消防员获批为烈士 美人鱼 胡锡进谈中美元首会晤 嘉南传 动物世界 苏宁易购回应破产传闻 #耿直真香哥黑化卖惨# 沙丘 许家印为恒大注入超70亿续命资金 大连一密接者擅自点外卖聚餐被调查 中美元首是否达成新共识?中方回应 俄方回应卫星碎片危及国际空间站 24岁救人牺牲消防员获批为烈士 星际穿越 林丹世界排名被正式移除 中美元首是否达成新共识?中方回应 中美元首是否达成新共识?中方回应 男子写80页PPT拯救爱情却离婚 林丹世界排名被正式移除 中美元首会谈重点内容 两个女人 红色通缉令 美人鱼 国足战澳大利亚大名单:4归化在列 罗永浩吐槽苹果文案没文化 24岁救人牺牲消防员获批为烈士 外交部回应拜登重申不支持台独 房价上涨城市创七年新低 拐点来了? 苏宁易购回应破产传闻 大连现超级传播者26人在同一传播链 罗永浩吐槽苹果文案没文化 动物世界 房价上涨城市创七年新低 拐点来了? 动保组织向上饶信州区申请信息公开 国足战澳大利亚大名单:4归化在列 你好李焕英 林丹世界排名被正式移除 得知母亲出事男子在地铁痛哭 动保组织向上饶信州区申请信息公开 意大利错失直接晋级世界杯资格
河源市| 德钦县| 荣昌县| 金门县| 腾冲县| 四川省| 辰溪县| 台北县| 千阳县| 高安市| 克山县| 武定县| 石泉县| 五常市| 杭锦后旗| 喜德县| 化州市| 栖霞市| 淮阳县| 乐至县| 太原市| 汝阳县| 克拉玛依市| 清苑县|