魏浩孫鳳舉呼義翔梁天學叢培天邱愛慈1)?

1)(西安交通大學,電力設備電氣絕緣國家重點實驗室,西安 710049)2)(西北核技術研究所,強脈沖輻射環(huán)境模擬與效應國家重點實驗室,西安 710024)(2016年6月30日收到;2016年11月20日收到修改稿)

一種非軸對稱磁絕緣電子鞘層邊界的計算方法?

魏浩1)2)?孫鳳舉2)呼義翔2)梁天學2)叢培天2)邱愛慈1)2)?

1)(西安交通大學,電力設備電氣絕緣國家重點實驗室,西安 710049)2)(西北核技術研究所,強脈沖輻射環(huán)境模擬與效應國家重點實驗室,西安 710024)(2016年6月30日收到;2016年11月20日收到修改稿)

長距離磁絕緣傳輸線內電極偏心、感應腔注入電流非均勻分布引起電子鞘層邊界偏心等非對稱磁絕緣特性.電子鞘層邊界是研究非軸對稱磁絕緣特性的重要參數(shù).本文提出一種計算非軸對稱磁絕緣電子鞘層邊界的方法.通過引入角向非均勻分布的模數(shù),將經典一維軸對稱Creedon穩(wěn)態(tài)磁絕緣理論推廣應用于圓柱坐標系下二維(r,θ)平面.建立了感應電壓疊加器次級非軸對稱磁絕緣的二維Creedon物理模型,給出了非軸對稱磁絕緣電子鞘層邊界的數(shù)值計算方法和計算誤差.當陰極角向磁場(陰極電流)角向分布滿足余弦函數(shù)時,電子鞘層邊界接近高斯分布.陰極電流角向不均勻程度越大,電子鞘層邊界偏心程度越嚴重,計算誤差越大.

磁絕緣感應電壓疊加器,注入電流非均勻分布,非軸對稱磁絕緣,電子鞘層邊界

1 引 言

磁絕緣傳輸線(magnetically insulated transmission lines,MITL)在超高功率脈沖傳輸和匯聚中具有重要應用[1?4].大型脈沖功率裝置常采用數(shù)米至幾十米長的MITL,長距離MITL內電極通常僅在首端與裝置連接固定,這種懸臂式支撐結構易造成MITL內電極偏心、電極間電/磁場非軸對稱分布,從而引起非軸對稱磁絕緣[4?7].此外,在磁絕緣感應電壓疊加器(magnetically induction voltage adders,MIVA)中,由于感應腔初級注入電流非均勻分布,也會導致MIVA次級磁絕緣呈現(xiàn)非軸對稱特性[7?11].

磁絕緣電子鞘層邊界是聯(lián)系MITL微觀特性(空間電子電荷密度、電荷質心、漂移速率)和宏觀特性(線電壓、陰/陽極電流、運行阻抗)的關鍵參數(shù).電子鞘層邊界限定了磁絕緣空間電子的活動區(qū)域,可視作MITL的虛陰極,電子鞘層邊界和陽極構成的傳輸線特征阻抗即是MITL的實際運行阻抗[12].沿電子鞘層邊界對磁場強度進行環(huán)路積分,可獲得MITL陽極電流.在非軸對稱磁絕緣研究中,電子鞘層邊界是評估磁絕緣狀態(tài)和性能的重要物理量.

在通常的一維、軸對稱磁絕緣中,電子鞘層邊界是與MITL陰/陽電極共圓心的同心圓 (rs≡常數(shù)).對于非軸對稱磁絕緣,電子鞘層邊界rs(θ)隨角向位置呈某種分布,即rs(θ)=f(θ).目前已有的磁絕緣模型(1975年Creedon層流理論[13],1979年Mendel[14]任意動量理論,1990年Lawconnell和Neri[15]通用Laminar理論,2006年Ottinger和Schumer[16]基于粒子模擬修正的Mendel理論,以及2006年Stygar等[17]磁絕緣電子彈性碰撞理論)主要用來描述一維、軸對稱磁絕緣問題.由現(xiàn)有模型很難直接獲得非軸對稱磁絕緣rs(θ)的解析解.

電子鞘層邊界rs(θ)也很難直接實驗測量,目前主要采用粒子模擬(particle in cell,PIC)方法,通過跟蹤大量磁絕緣空間電子的軌跡來確定[16].確定非軸對稱磁絕緣的rs(θ)必須采用三維 PIC模型,這對PIC模擬軟件及代碼提出了較高要求.國外主要采用大型三維粒子模擬軟件LSP(large scale procedure),國內電子科技大學研發(fā)了三維粒子模擬軟件CHIPIC3D,并用于“聚龍一號”裝置四層圓盤錐MITL和1 MV LTD裝置螺旋支撐桿MITL的三維粒子模擬[18?21].但是,大型脈沖功率裝置或器件的三維PIC計算耗時極長,模擬結果因受多種因素和假設條件的影響隨機性很大.因此,本文提出了一種確定MIVA次級非對稱磁絕緣電子鞘層邊界rs(θ)的方法.

2 MIVA次級非軸對稱磁絕緣物理模型

圖1為MIVA次級MITL橫截面(垂直功率傳輸方向)示意圖.rc,ra分別為MITL陰、陽極半徑,rs(θ)為待求解的磁絕緣電子鞘層邊界.MIVA感應腔初級注入電流非均勻分布導致次級磁絕緣陰/陽極電流、電/磁場角向非軸對稱分布[22,23].即使MITL陰/陽電極同軸,仍然會導致rs(θ)偏心[24].

由磁絕緣Creedon層流理論[13],磁絕緣達到穩(wěn)態(tài)時空間電子區(qū)域(rc＜r＜rs(θ))和真空區(qū)域(rs(θ)＜ r ＜ ra) 的相對論因子γi(r,θ)和γo(r,θ)分別滿足數(shù)學物理方程[25?28],

在陰極(r=rc)和陽極 (r=ra)的邊界位置,γi(r,θ)和γo(r,θ)分別滿足如下邊界條件,

其中γ0=1+eV0/(mc2)為陽極的相對論因子,e為電子電荷量,m為電子靜止質量,V0為陰-陽電極間線電壓.

在電子鞘層邊界rs(θ)處,滿足分界面銜接條件,

采取分離變量法[29]求解方程(1),考慮到γi(r,θ)和γo(r,θ)在[0,2π]內關于θ= π 對稱,并代入邊界條件(2),得到相對論因子γi(r,θ)和γo(r,θ)的解析解分別為

其中n為非負整數(shù),表征陰/陽極電流角向非均勻分布的模數(shù),若令n=0,(4)式即為經典的一維、軸對稱磁絕緣Creedon層流模型[13,28];系數(shù)aj,bj(0≤j≤n,j為整數(shù))為待定常數(shù),其取值與陰/陽極電流角向分布有關.

當給定MITL結構參數(shù)(陰/陽極半徑ra,rc)、線電壓V0和陰極電流角向分布Ic(θ)時,聯(lián)合方程(3)和(4),可確定電子鞘層邊界rs(θ).但因方程(4)為超越方程,難以直接獲得rs(θ)的解析解,需要采用數(shù)值方法進行求解.

3 計算電子鞘層邊界rs(θ)的數(shù)值方法

3.1 確定系數(shù)aj(0≤j≤n)

由磁絕緣穩(wěn)態(tài)時空間電子徑向受力平衡[27],得到

(4)和(5)式聯(lián)合,得到在rc＜ r＜ rs(θ)區(qū)域Bθ(r,θ)滿足,

令(6)式中r=rc,得到陰極表面角向磁場Bθ(rc,θ),

采用正交三角函數(shù)系cos(jθ)將Bθ(rc,θ)在[0,2π]區(qū)間內做Fourier級數(shù)展開,得系數(shù)aj,其中κ為常數(shù),κ=?erc/mc.

由(8)式可知,系數(shù)aj取決于陰極表面角向磁場Bθ(rc,θ)角向分布.實際中,Bθ(rc,θ)角向分布可由感應電壓疊加器的電磁模擬或實驗測量獲得.

3.2 確定系數(shù)bj(0≤j≤n)

由方程(3b)可以推導系數(shù)bj(0≤j≤n)滿足如下矩陣方程,

其中a=(a0,a1,...,an)T和b=(b0,b1,...,bn)T為n+1階向量.N(rs),M(rs)為(n+1)×(n+1)階矩陣,其元素分別為

(9)和(10)式的推導過程見附錄A.(9)和(10)式表明,當系數(shù)aj確定后,系數(shù)bj僅與電子鞘層邊界rs(θ)有關.因此,可依據(jù)MIVA次級磁絕緣電子鞘層邊界的物理實際和可能分布,假定rs(θ)滿足某種特定的函數(shù)分布f0(θ)(如高斯分布、余弦分布、二次拋物線分布等),f0(θ)中含有待定系數(shù).

3.3 確定f0(θ)的待定系數(shù)

采用最優(yōu)化問題的全局搜索解法[30],確定函數(shù)f0(θ)中待定系數(shù).最優(yōu)化問題可表述為

其中(11a)式為最優(yōu)化問題的目標函數(shù),表征當rs(θ)為假定函數(shù)分布f0(θ)時, γi(r, θ)和 γo(r, θ)偏離方程(3a)的程度,定義為計算誤差.搜索約束條件(11b)限定了電子鞘層邊界的取值范圍,即rc＜ f0(θ)＜ra.

3.4 確定電子鞘層邊界rs(θ)

重復3.2–3.3節(jié)的步驟,比較多種假定函數(shù)分布(高斯分布、余弦分布、二次拋物線分布等)的計算誤差,計算誤差最小的分布視為電子鞘層邊界rs(θ).

4 計算方法在MIVA次級非軸對稱磁絕緣的應用

給定MIVA次級MITL運行參數(shù):陰極半徑rc=0.1 m,陽極半徑ra=0.2 m,陰-陽極線電壓V0=4 MV,γ0=9,磁絕緣運行在最小電流工作點.

在一維軸對稱情況(陰/陽極電流角向均勻分布)下,由磁絕緣經典Creedon層流理論,計算磁絕緣電子鞘層邊緣的相對論勢γs=2.0,電子鞘層邊界rs=0.119 m,陽極電流Ia=132 kA,陰極電流Ic=65.5 kA,陰極表面角向磁場Bco=0.131 T.

二維、非軸對稱情況下,由于MIVA感應腔初級注入電流角向的非均勻分布,靠近饋入點的角向位置(θ=0)Bθ(rc,θ)明顯偏大,遠離饋入點的角向位置(θ= π)Bθ(rc,θ)偏小.假定角向磁場Bθ(rc,θ)滿足余弦分布,

其中Bco為一維軸對稱(電流角向均勻分布)時陰極角向磁場,不均勻系數(shù)δ用于表征二維非軸對稱情況下Bθ(rc,θ)角向非均勻分布的程度.

假定MIVA次級MITL電子鞘層邊界rs(θ)滿足余弦分布、高斯分布和二次拋物線分布,其分布函數(shù)分別如(13)–(15)式所示,

由本文所述的數(shù)值計算方法,可以確定f1(θ),f2(θ)和f3(θ)中待定系數(shù).三種假定分布函數(shù)下電子鞘層邊界rs(θ)如圖2所示.余弦、高斯和二次拋物線三種分布函數(shù)下的計算誤差分別為2.7084,0.0682和9.6110.高斯分布的計算誤差最小,因此,當陰極角向磁場Bθ(rc,θ)滿足余弦分布時,MITL電子鞘層邊界rs(θ)最接近高斯分布.

改變Bθ(rc,θ)角向不均勻系數(shù)δ,電子鞘層邊界rs(θ)隨之變化.表1為幾種不均勻系數(shù)δ時三種分布函數(shù)的計算誤差.由表1可知,二次拋物線分布的計算誤差最大,高斯分布的計算誤差均最小.注入電流分布越均勻(角向磁場Bθ(rc,θ)不均勻系數(shù)δ越小),計算誤差越小.

當Bθ(rc,θ)不均勻系數(shù)δ分別為25%,10% 和2.5%時,采用高斯函數(shù)分布求解的電子鞘層邊界如圖3所示.電子鞘層邊界rs(θ)的最大值、最小值和電子鞘層不均勻系數(shù)ζ如表2所列,其中電子鞘層角向不均勻系數(shù)ζ定義為

其中max(rs(θ)),min(rs(θ))和mean(rs(θ)) 分別為rs(θ)在θ∈[0,2π]的最大值、最小值和平均值.(16)式定義的ζ物理含義為電子鞘層厚度的極差與平均厚度之比.

圖2 (網刊彩色)三種假定分布函數(shù)時電子鞘層邊界(Bθ(rc,θ)不均勻系數(shù)δ=25%)Fig.2.(color online)The electron sheath profi les under the three assumed dstribu tion functions as the asymmetric coffi cient of Bθ(rc,θ),δ,is abou t 25%.

表1 Bθ(rc,θ)不均勻系數(shù)δ對計算誤差的影響Tab le 1.The in fl uences of the asymmetric coeffi cient of Bθ(rc,θ),δ,on the calcu lation errors.

表2 Bθ(rc,θ)角向不均勻系數(shù)δ不同時電子鞘層邊界rs(θ)的特性Tab le 2.The parameters of the electron sheath profi les function as the asymmetric coeffi cients of Bθ(rc,θ).

由圖3和表2可知,Bθ(rc,θ)不均勻系數(shù)δ越大(注入電流角向不均勻程度越大),電子鞘層邊界rs(θ)偏心程度越嚴重.電子鞘層角向不均勻系數(shù)ζ與Bθ(rc,θ)角向不均勻系數(shù)δ數(shù)值相接近.當Bθ(rc,θ)分布較均勻(δ=0.5%)時,電子鞘層邊界rs(θ)≈0.119 m,此數(shù)值與一維軸對稱Creedon理論計算的結果相同.

圖3 不同Bθ(rc,θ)角向不均勻系數(shù)δ時的磁絕緣電子鞘層邊界 (a)δ=25%;(b)δ=10%;(c)δ=2.5%Fig.3.The electron sheath profi les under th ree different asymmetric coeffi cients of Bθ(rc, θ)for(a)δ=25%,(b)δ=10%,and(c)δ=2.5%.

5 結 論

在穩(wěn)態(tài)磁絕緣經典、一維軸對稱Creedon理論基礎上,引入表征陰/陽極電流角向非均勻分布的參數(shù)(角向非均勻分布模數(shù)n),建立了圓柱坐標下MIVA次級非軸對稱磁絕緣的二維Creedon物理模型,給出了磁絕緣電子鞘層邊界rs(θ)的數(shù)值計算方法和計算誤差.當給定磁絕緣線結構參數(shù)(陰/陽極半徑)、線電壓和陰極電流(陰極角向磁場)角向分布時,可計算磁絕緣電子鞘層邊界rs(θ).本文提出的磁絕緣電子鞘層邊界確定方法具有計算效率高、耗時短等優(yōu)點,該方法已用于MIVA次級非軸對稱磁絕緣電子鞘層邊界的確定.結果表明,當陰極角向磁場Bθ(rc,θ)角向分布滿足余弦函數(shù)時,電子鞘層邊界rs(θ)呈高斯分布.研究了電流均勻性對磁絕緣電子鞘層邊界rs(θ)的影響,結果表明,隨著電流不均勻程度增加,電子鞘層偏心程度加大;電子鞘層rs(θ)的角向不均勻系數(shù)和陰極角向磁場Bθ(rc,θ)的不均勻系數(shù)數(shù)值相近.

感謝西安交通大學數(shù)學學院李東升教授、電氣工程學院馬西奎教授和葛曉宇博士對數(shù)值求解方法給予的良好建議.感謝西北核技術研究所付梅艷、尹佳輝、姜曉峰、曾江濤、張鵬飛、來定國、孫江等的討論.

附錄A 矩陣方程(9)的推導

分界面銜接條件(3)的物理意義是在分界面rs(θ)上,電勢和法向電場連續(xù).

當電子鞘層偏心程度較輕(即(rs?rc)/rc? 1)時,?/?n ～= ?/?r,(3b)式可簡化為

將(4)式中γi(r,θ)和γo(r,θ)表達式代入(A1)式,得到

為了便于分析,將(A2)式等號兩邊均寫成級數(shù)求和形成,即

其中

由于方程(A3)對于[0,2π]內任意θ均成立,因此,將θ離散化處理,在[0,2π]取(n+1)個離散點,

對于任意一個θi均滿足方程(A3),可以得到,

將公式(A6)寫成矩陣方程的形式,即

其中,a=(a0,a1,...,an)T和b=(b0,b1,...,bn)T為n+1階向量;N(rs),M(rs)為(n+1)×(n+1)階矩陣,其元素分別為

[1]Stygar W A,Corcoran P A,Ives HC,Spielman R B,Douglas JW,W hitney BA,MostromMA,W agoner TC,Speas C S,G illiland TL,Allshouse G A,C lark R E,Donovan G L,Hughes TP,Humph reys D R,JaramilloD M,Johnson MF,Kellogg JW,Leeper R J,Long F W,Martin TH,Mu lville TD,Pelock MD,Pey ton BP,Poukey JW,Ramirez J J,Reynolds P G,Seamen J F,Seidel D B,Seth AP,Sharpe AW,ShoupR W,Smith J W,Vandevalde D M,W av rik R W 2009Phys.Rev.Spec.Top.Accel.Beams.12 120401

[2]Stygar W A,Awe TJ,Bailey J E,Bennett N L,E W Breden,Campbell E M,C lark R E,Cooper R A,CuneoME,Ennis J B,Feh l D L,Genoni TC,Gomez MR,G reiser G W,G runer F R,Herrmann MC,Hutsel BT,Jennings C A,Jobe D O,Jones BM,MC Jones,Jones P A,KnappP F,Lash J S,LeChien KR,Leckbee J J,Leeper R J,Lewis SA,Long FW,LuceroD J,Mad rid E A,Martin MR,Matzen MK,Mazarakis MG,McBride R D,McKee G R,Miller C L,Moore J K,MostromC B,Mu lville TD,Peterson KJ,Porter J L,Reisman D B,Rochau G A,Rochau G E,Rose D V,Rovang D C,Savage ME,Sceiford ME,Schmit P F,Schneider R F,Schwarz J,SefkowAB,Sinars D B,Slutz SA,Spielman R B,Stoltzfus BS,Thoma C,Vesey R A,W akeland P E,W elch D R,W isher ML,W oodworth J R 2015Phys.Rev.Spec.Top.Accel.Beams.18 110401

[3]Zou W K,W ang M,Chen L,Zhou L J,GuoF,X ie W P,Deng J J 2013High Power and Laser and Particle Beams25 2487(in Chinese)[鄒文康,王勐,陳林,周良驥,郭帆,謝衛(wèi)平,鄧建軍2013強激光與粒子束25 2487]

[4]Pate R C,Patterson J C,Dowd ican MC,Ramirez J J,Hasti D E,Tolk KM,Poukey JW,Schneider L X,Rosenthal SE,San ford TW,Alexander J A,Heath C E 1987Proceeding of 6th IEEE Pu lsed Power ConferenceArlington,USA,22–26 June,1987 p478

[5]Johnson D,Bailey V,Altes R,Corcoran P,Smith I,Cordova S,Hahn K,Maenchen J,Molina I,PortilloS,Puetz E,Sceiford M,Valde D V,Rose D,Oliver B,W elch D,D roemer D 2005Proceeding of 15h IEEE Pu lsed Power ConferenceMonterey,USA,13–17 June,2005 p314

[6]Allen R J,CommissoR J,Cooperstein G,Ottinger P F,Schumer JW 2011Proceeding of 18th IEEE Pu lsed Power ConferenceChicago,USA,19–23 June,2011 p1036

[7]Smith ID 2004Phys.Rev.Spec.Top.Accel.Beams.7 064801

[8]Smith ID,Bailey V L,Fock ler J,Gustwiller JS,Johnson D L,Maenchen J E,D roemer D W 2000IEEE Trans.P lasma Sci.28 1653

[9]Thomas K,Beech P,Brown S,Buck J,Bu rscough J,C lough S,Crotch I,Du ff Y J,Goes C,Huckle I,Jones A,King A,Stringer B,Th readgold J,Trenaman S,W heeldon R,W ood roofe M,Carboni V,DaSilva T,Galver B,G lazebrook W,Hanzel K,Pearce J,PhamJ,Pomeroy S,SaundersW,Speits D,Warren T,W hitney B,W ilson J2011Proceeding of18th IEEE Pu lsed Power ConferenceChicago,USA,19–23 June,2011 p1042

[10]W ei H,Sun F J,Liang TX,Y in J H,Dang TF,Zeng J T,Cong P T,Qiu AC 2014Rev.Sci.Instrum.85 083506

[11]W ei H,Sun F J,Liang TX,GuoJ M,Q iu AC,Cong P T,Y in J H,Hu Y X,Jiang X F,W ang Z G,Dang TF 2015IEEE Trans.P lasma Sci.43 3359

[12]Mendel C W,Rosentha SE 1995Phys.Plasmas2 1332

[13]C reedon J M1975J.Appl.Phys.46 2946

[14]Mendel C W 1979J.Appl.Phys.50 3830

[15]Lawconnell R I,Neri J 1990Phys.F luids B2 629

[16]Ottinger P F,Schumer J W 2006Phys.P lasmas13 063109

[17]Stygar W A,W agoner TC,Ives HC,Corcoran P A,CuneoME,Douglas JW,G illiland TL,Mazarakis MG,Ramirez J J,Seamen J F,Seidel D B,Spielman R B2006Phys.Rev.Spec.Top.Accel.Beams.9 090401

[18]Liu L Q,Liu D G,W ang X Q,Zou W K,Yang C 2012Acta Phys.Sin.61 162901(in Chinese)[劉臘群,劉大剛,王學瓊,鄒文康,楊超2012物理學報61 162901]

[19]Liu L Q,Liu D G,W ang X Q,Yang C,X ia MZ,Peng K2012Acta Phys.Sin.61 162902(in Chinese)[劉臘群,劉大剛,王學瓊,楊超,夏蒙重,彭凱2012物理學報61 162902]

[20]DaiW,Zhu X Z,Liu LQ,Zou W K,GuoF,Liu D G 2015Modern E lectron.Tech.38 144(in Chinese)[戴巍,朱旭智,劉臘群,鄒文康,郭帆,劉大剛 2015現(xiàn)代電子技術38 144]

[21]Zou W K,Chen L,Liu D G,Zhang L,Liu L Q,Zhou L J,W ang M,W ei B,GuoF,W ang X Q,Dai Y M,X ie W P,Deng J J 2012Phys.Rev.Spec.Top.Accel.Beams15 110401

[22]Johnson D L,Smith ID,Corcoran P,Bailey V,Maenchen J,Rovang D,Molina I,Hahn K,LuceroR,Kincy M,Kitterman D,Oliver B,W elch D,Rose D,Goldsack TJ,Phillips MA,Sinclair MA,Thomas KJ 2002Proceeding of 14th IEEE High Power Particle Beams ConferenceAlbuquerque,USA,23–28 June,2002 p123

[23]Oliver BV,Genoni TC 2003Proceeding of 14th IEEE Pu lsed Power ConferenceDallas,USA,June 15–18,2003 p395

[24]Bruner N,MostromC,Rose D V,W elch D R,Bailey V,Johnson D L,Oliver BV 2007Proceeding of 16th IEEE Pu lsed Power ConferenceAlbuquerque,USA,June 17–22,2007 p807

[25]Q ian BL,Liu Y G,Li C L,Liu C B2000IEEE Trans.P lasma Sci.28 760

[26]Fuks MI,Schamiloglu E 2014Phys.Plasmas21 053102

[27]BluhmH(translated by Jiang W H,Zhang C)2008Pu lsed Power Systems:Principles and Application(Beijing:Tsinghua University Press)pp128–135(in Chinese)[布魯姆H著(江偉華,張弛譯)2008脈沖功率系統(tǒng)的原理與應用 (北京:清華大學出版社)第128–135頁]

[28]W ang MY 1978Appl.Phys.Lett.33 284

[29]Ma X K2000Theory and Application of Electromagnetic Field(X i’an:X ian Jiaotong University Press)pp206–216(in Chinese)[馬西奎 2000電磁場理論及應用(西安:西安交通大學出版社)第206–216頁]

[30]Xue D Y Chen Y Q 2011Technology and Application of SystemSimu lation Based on MATLABand Simu link(Beijing:Tsinghua University P ress)pp108–110(in Chinese)[薛定宇,陳陽泉 2011基于MATLAB/Simu link的系統(tǒng)仿真技術與應用(北京:清華大學出版社)第108–110頁]

PACS:84.70.+p,07.85.Fv,41.20.–q,52.35.–gDOI:10.7498/aps.66.038402

Amethod tocalcu late the electron sheath profi le of the nonaxisymmetricalmagnetic insu lation?

Wei Hao1)2)?Sun Feng-Ju2)Hu Yi-Xiang2)Liang Tian-Xue2)Cong Pei-Tian2)Qiu Ai-Ci1)2)?

1)(State Key Laboratory of E lectrical Insu lation and Power Equipment,X i’an Jiaotong University,X i’an 710049,China)2)(State Key Laboratory of Intense Pu lsed Rad iation Simu lation and Effect,Northwest Institute of Nuclear Technology,Xi’an 710024,China)(Received 30 June 2016;revised manuscript received 20 November 2016)

The nonaxisymmetricalmagnetic insulation would occur due tothe disalignment of inner electrodes in longmagnetically insulated transmission lines,or the nonuniformdistributionsof injected currents in induction cavities ofmagnetically insulated induction voltage adders(MIVA).The electron sheath profi le is a very important parameter tocharacterize the nonaxisymmetricalmagnetic insulation.In the past,the three-dimensional particle in cell simu lation was usually used todetermine the electron sheath profi le,which is extremely time-consuming and ineffi cient.In this paper,a fast and effi cient calcu lation method is proposed.The classical one-dimensional Creedon theory of the magnetic insu lation equilibriumis generalized toa two-dimensional plane of(r,θ)via introducing a parameter defined as the azimuthalmode number.Two-dimensional Creedon is developed tomodel the asymmetric magnetic insulation of the MIVA.Provided the azimuthal distributions ofmagnetic fl ux density on the cathode,which is in proportion tothe cathode current,the two-dimensional Creedon model is numerically solved.Anumerical solution method tocalculate the electron sheath profi le is proposed,and then the calculation error is alsogiven.As the azimuthal distribution ofmagnetic flux density on the cathodemeets a cosine function,the profi le of the electron sheath is approximate tothe Gauss function.As the nonuniformportion of cathode current increases,the electron sheath becomesmore eccentric,and the calculation error is alsomuch larger.

magnetically insulated induction voltage adders,asymmetrical distribution of injected currents,nonaxisymmetricalmagnetic insulation,electron sheath profi le

10.7498/aps.66.038402

?國家自然科學基金(批準號:11505138,51577156)資助的課題.

?通信作者.E-mail:weihaoyy@nint.ac.cn

?通信作者.E-mail:qiuac@cae.cn

*Project supported by the National Natural Science Foundation of China(G rant Nos.11505138,51577156).

?Corresponding author.E-mail:weihaoyy@nint.ac.cn

?Corresponding au thor.E-mail:qiuac@cae.cn