楊 實, 黃種亮, 任書慶
(西北核技術(shù)研究所, 西安 710024; 強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國家重點實驗室, 西安 710024)
多針?biāo)_關(guān)是快脈沖低阻抗脈沖功率源脈沖成形開關(guān)采用的主要技術(shù)之一,多用于低能強(qiáng)流電子束源和軔致輻射模擬源[1-4]。作為脈沖成形開關(guān),開關(guān)電感的大小決定輸出脈沖前沿的快慢。因此衡量多針?biāo)_關(guān)性能優(yōu)劣的關(guān)鍵指標(biāo)之一就是開關(guān)電感的大小。
目前,西北核技術(shù)研究所的兩臺低阻抗大型脈沖功率設(shè)備“強(qiáng)光一號”和“閃光二號”加速器[5-6]均采用多針?biāo)_關(guān)作為脈沖成形開關(guān),兩臺設(shè)備的多針?biāo)_關(guān)電感均約60 nH,這種較大的開關(guān)電感導(dǎo)致兩臺加速器的輸出前沿略緩。反觀具有相同技術(shù)途徑和功率水平的美國Blackjack3加速器,其多針?biāo)_關(guān)的電感僅20 nH。升級后的美國Blackjack5加速器電壓為5 MV,阻抗為0.5 Ω,多針?biāo)_關(guān)的電感達(dá)到了10 nH??梢钥闯觯瑢τ诳烀}沖成形開關(guān)而言,電感值的差距非常大。因此,如何降低多針?biāo)_關(guān)的電感,找到多針?biāo)_關(guān)中影響開關(guān)電感的因素是目前需要解決的主要問題之一。
由于多針?biāo)_關(guān)放電時電壓高(MV級)、電流大(MA級)、電磁干擾強(qiáng),因此很難對每只放電通道的電流及導(dǎo)通情況進(jìn)行準(zhǔn)確的實驗測量和分析。本文以“閃光二號”加速器多針?biāo)_關(guān)為模型,通過理論分析和數(shù)值模擬研究了多針?biāo)_關(guān)放電通道電流均勻性對開關(guān)電感的影響,探索了多針?biāo)_關(guān)電感變化的機(jī)理。
“閃光二號”加速器脈沖形成段結(jié)構(gòu)如圖1所示,其中紅線框內(nèi)為多針?biāo)_關(guān)。多針?biāo)_關(guān)針電極由9只φ60 mm的不銹鋼棒組成,均布于φ1 410 mm的圓周,開關(guān)區(qū)間距為255 mm,水線外筒內(nèi)徑為φ2 300 mm。
圖1 “閃光二號”加速器脈沖形成段結(jié)構(gòu)圖Fig.1 The pulse forming part of Flash-II accelerator
根據(jù)多針?biāo)_關(guān)的實際結(jié)構(gòu),針對9個放電通道全部導(dǎo)通的情況,將每個通道與外筒等效成一個回路,則開關(guān)的總電感為9個等效回路相互耦合的總電感。根據(jù)磁場儲能與電流、電感之間的關(guān)系,總電感與每個回路的自感及互感存在式(1)關(guān)系:
(1)
可以得到
(2)
其中,Wtotal為開關(guān)區(qū)磁場總能量;Ltotal和Itotal分別為開關(guān)區(qū)總電感和總電流;Li和Ii分別為第i個回路的自感和電流;Mij為第i個回路與第j個回路之間的互感。對多針?biāo)_關(guān)電極進(jìn)行編號,放電通道排列順序如圖2所示。
圖2 放電通道排列順序分布圖Fig.2 Sequence of the discharge channels
(3)
由于兩個等效回路之間的互感值由回路之間的相對位置關(guān)系決定,所以當(dāng)?shù)刃Щ芈分g的相對位置相同時,則對應(yīng)的Mij相等,且有Mij=Mji,且認(rèn)為當(dāng)?shù)刃Щ芈分睾蠒r,所得互感即為自感,考慮圓周對稱性,則可得到
M0=M11=…=M99
M1=M12=M21=…=M19=M91
M2=M13=M31=…=M18=M81=M29=M92
M3=M14=M41=…=M17=M71=M28=M82=M39=M93
M4=M15=M51=…=M16=M61=M27=M72=M38=M83=M49=M94
即得
互感的大小可由式(4)計算:
(4)
以放電通道1為例,可以認(rèn)為放電通道與加速器外筒構(gòu)成偏心同軸線。采用保角變換方法將開關(guān)橫截面變換到ζ平面,從而將通道1與外筒之間變換成同心圓,再采用無限長同軸電纜的磁場強(qiáng)度計算公式即可求得ζ平面中等效回路1的磁場強(qiáng)度,進(jìn)而求得該回路互感。
基于分式線性變換的保圓性和對稱不變性,對Z平面作分式線性變換到ζ平面上,變換表達(dá)式:
(5)
其中,z為Z平面上的點,z=x+iy;x1和x2分別為z平面上通道1和外筒圓心的鏡像對稱點。則
(6)
求得
(7)
其中,外筒半徑R1為1 150 mm;多針電極半徑R2為30 mm;多針電極軸心與外筒軸心的距離d為705 mm。
在ζ平面上,通道1和外筒同心,變換后通道1和外筒的半徑分別為
(8)
由于R1?R2和d?R2,因此可將放電通道2,3,4看作各自圓心上的一個點。通過變換可求得在ζ平面上通道2,3,4對應(yīng)的半徑分別為ρ2,ρ3,ρ4。因此,互感的表達(dá)式變換為
(9)
針對5種不同的電流分布情況,對多針?biāo)_關(guān)電感進(jìn)行了計算,放電通道排列順序如圖2所示。
1) 只有一個通道導(dǎo)通,即Ik=1,Ii=0(i≠k),則
L1=ITMI=M0=162 nH
根據(jù)開關(guān)電感表達(dá)式及約束條件,構(gòu)造拉格朗日函數(shù),可以得到:
(10)
其中,λ為Lagrange乘子。則條件極值點須滿足下列方程組:
(11)
其中,Bi(i=1,2,…,9)是由矩陣M的第i行組成的行向量。
將方程組(11)中的9個方程相加,即得:
(12)
將λ值代入式(11)可以得到:
上述計算過程證明了當(dāng)開關(guān)通道電流均相等時,開關(guān)電感具有極值。下面證明這個極值為最小值。
取Ltotal在極值點I0的海賽矩陣為
(13)
容易求得K的特征值為121,121,131,131,153,153,190,190,266。所以K為正定矩陣,即W(I0)也為正定矩陣,因此上述所求極值為Ltotal的極小值。
采用電磁場模擬軟件對“閃光二號”加速器多針?biāo)_關(guān)進(jìn)行了等尺寸的數(shù)值模擬研究,圖3為開關(guān)模擬結(jié)構(gòu)剖面圖。開關(guān)模型基本尺寸與“閃光二號”加速器多針?biāo)_關(guān)區(qū)一致, 9只電極、放電通道、電極安裝板及對應(yīng)的外筒區(qū)構(gòu)成模擬開關(guān)電感的回路。
圖3 多針?biāo)_關(guān)模擬結(jié)構(gòu)剖面圖Fig.3 Cross section of the multi-channel water switch used in 3D numerical calculation
根據(jù)磁場能量與電流、電感的關(guān)系式
(14)
可以得出,當(dāng)電流一定時,磁場能量與電感成正比。
根據(jù)式(14),假設(shè)9只電極均導(dǎo)通,即通道數(shù)為9,并且9個通道傳導(dǎo)的電流之和為定值,則可以通過設(shè)定9個通道傳導(dǎo)不同大小的電流,模擬出不同的磁場能量,進(jìn)而求得開關(guān)電感。
這里,假設(shè)9個通道的總電流為270 kA,僅選擇與理論分析部分一致的5種典型情況來說明電流均勻性對開關(guān)電感的影響。通道分布順序如圖2所示。
1) 9個通道中只有通道1導(dǎo)通,導(dǎo)通電流為270 kA,模擬得到該情況下開關(guān)電感為153 nH。
2) 9個通道均導(dǎo)通,電流均為30 kA,模擬得到開關(guān)電感為27 nH。
3) 1,2,3,4,5通道電流為10 kA,其余通道電流為55 kA,模擬得到開關(guān)電感為38 nH。
4) 1,2,4,6,8通道電流為10 kA,其余通道電流為55 kA,模擬得到開關(guān)電感為34 nH。
5) 1,2,4,6,8通道電流為1 kA,其余通道電流為66.25 kA,模擬得到開關(guān)電感為42 nH。
模擬結(jié)果表明,多針?biāo)_關(guān)放電通道電流均勻性越好,開關(guān)電感越小,“閃光二號”多針?biāo)_關(guān)的電感最小值為27 nH,電感最大值為153 nH;各通道電流相同的情況下,分布位置的不同會導(dǎo)致電感的變化;電流分布位置相同的情況下,各通道電流不同也會導(dǎo)致電感的變化。
通過保角變換的方法分析了通道電流不均勻?qū)Χ噌標(biāo)_關(guān)電感的影響,結(jié)果表明,多針?biāo)_關(guān)放電通道電流不均勻性越大,回路電感越大。多針?biāo)_關(guān)電感在各通道電流均勻時,取得極小值。使用電磁場模擬的方法驗證了理論計算的正確性,得到了“閃光二號”加速器多針?biāo)_關(guān)電感的最小值和最大值分別為27 nH和153 nH。
[1]MILLER A R. Sub-ohm coaxial pulse generators, Blackjack3,4 and 5 [C]//3rd IEEE International Pulsed Power Conference, 1981: 200-205.
[2]JOHNSON D L, VANDEVENDER J P, MARTIN T H. High power density water dielectric switching [J]. IEEE Trans Plasma Sci, 1980, 8(3): 204-209.
[3]叢培天, 蒯斌, 邱愛慈, 等. 脈沖電壓下的自擊穿水介質(zhì)開關(guān)擊穿特性和電路參數(shù)實驗研究[J]. 強(qiáng)激光與粒子束, 2005, 17(9): 1 414-1 418.(CONG Pei-tian, KUAI Bin, QIU Ai-ci, et al. Experimental study on circuit parameters and breakdown properties of water switch under pulsed voltage[J]. High Power Laser and Particle Beams, 2005, 17(9): 1 414-1 418.)
[4]MILLER A R, GILBERT C, RAUCH J, et al. Pulsed power design for a compact X-ray simulator[C]//12th IEEE International Pulsed Power Conference, 1999: 271-274.
[5]黃濤, 叢培天, 張國偉, 等. “強(qiáng)光一號”加速器電路模擬與分析[J]. 強(qiáng)激光與粒子束, 2010, 22(4): 897-900. (HUANG Tao, CONG Pei-tian, ZHANG Guo-wei, et al. Circuit simulation and analysis for Qiangguang I accelerator[J]. High Power Laser and Particle Beams, 2010, 22(4): 897-900.)
[6]楊實, 任書慶, 叢培天, 等. “閃光二號”加速器適應(yīng)性改造[J]. 強(qiáng)激光與粒子束, 2016, 28(1): 015103. (YANG Shi, REN Shuqing, CONG Peitian, et al. Transformation of the Flash-Ⅱ accelerator[J]. High Power Laser and Particle Beams, 2016, 28(1): 015103.)