(海軍工程大學電氣工程學院 武漢 430033)
凈輻射系數法是由Lowke在1974年提出的,用來計算電弧等離子體高溫區域輻射能量損失的方法,其假定沒有考慮電弧自身對輻射熱的吸收問題[1~6]。很多學者針對不同的氣體介質進行了計算研究:J Menart[7]研究了從3000K~25000K 氬氣等離子體中存在不同比例(0、0.1%、1%、10%、100%)的鐵元素的特定電弧等離子體的凈輻射系數隨電弧溫度與電弧長度變化關系,其結論是在鐵元素占電弧比例相同、相同氣壓及相同電弧長度下,凈輻射系數隨溫度的升高而增大;在鐵元素占電弧比例相同、相同氣壓和相同溫度下,凈輻射系數隨電弧長度的增大而減?。辉谙嗤瑲鈮?、相同溫度及相同電弧長度下,凈輻射系數隨鐵元素占電弧比例的增加而增大。Y Naghizadeh-Kashani[8]研究了在假設熱力學平衡和等離子體等溫的情況,且在大氣壓下溫度范圍為3000K~40000K 的空氣等離子體的凈輻射系數;結果顯示在溫度高于6000K時,空氣等離子體的輻射特性主要由大量的氮離子發揮作用的;溫度低于6000K 時,空氣等離子體的輻射特性主要由大量的氧氣和一氧化氮發揮作用的。X Zhang[9]研究了在中壓電氣裝置中由于故障電弧引起的溫度的計算,為了確定故障電弧散發的熱量,必須求得電弧的凈輻射系數。Vladimir Aubrecht[10]研究了空氣和SF6等溫等離子體在外界各種不同氣壓、電弧半徑在0.01cm~10cm 以及電弧溫度從5000K~30000K 的凈輻射系數。
由于國內外這方面的研究還較少,且開關觸頭材料的不同,導致產生的空氣電弧中含有的金屬元素也不一樣,為了研究不同觸頭材料的高速混合型限流熔斷器中開斷器的開斷特性,需要分析電弧的熱平衡關系,這就使得有必要研究一些含金屬元素的空氣電弧的凈輻射系數,本文就實驗室的需要和現有條件進行了銅絲引弧的空氣電弧平板實驗,計算得到含銅元素的空氣電弧的凈輻射系數。
實驗系統由實驗電路、銅絲引弧空氣電弧平板實驗裝置、高速攝像儀、示波器及測溫儀器五個部分組成。其中空氣電弧平板實驗裝置中兩極板材料為T1 銅,純度99.95%,厚度2.98mm,銅絲為純銅,直徑50μm;高速攝像儀是每10μs拍攝一張圖片。設計了銅觸頭銅絲引弧空氣電弧平板實驗電路,如圖1所示;該實驗電路主要由電流回路與自然換流回路組成。電路參數:電路主回路電容C為50mF,電阻R=1Ω,F1與F2均為晶閘管,試品為空氣電弧平板實驗裝置,如圖2所示;利用該電路進行了電容充電100V 或200V 和極板間距為1mm、2mm、3mm 等一系列的銅絲引弧空氣電弧平板實驗。

圖1 空氣電弧平板實驗電路

圖2 空氣電弧平板實驗裝置
在電容充電100V,極板間距為1mm 的實驗條件下,使晶閘管F1導通,電容放電;延時600μs,晶閘管F2導通,回路進行自然換流,直至試品中電弧熄滅;實驗結束后,通過示波器得到電弧的電流和電壓隨時間變化的波形圖,如圖3所示;高速攝像機記錄一組銅絲引弧空氣電弧平板實驗從起弧到熄弧的高速攝像圖片,如圖4所示。

圖3 電弧的電流和電壓實驗波形

圖4 銅絲引弧空氣電弧平板實驗高速攝像
1)電弧直徑及熱邊界層直徑的量取
任取高速攝像儀記錄的空氣電弧圖片一張,將其導入Autocad軟件中,標注兩極板間的長度、電弧長度、電弧直徑及熱邊界層的直徑分別得到四個數據,如圖5所示,其中認為圖中最亮部分的寬度為電弧直徑、電弧周圍稍暗的部分為電弧熱邊界層,則圖中4.05為兩極板間長度與電弧長度、7.08為電弧直徑、11.67為電弧熱邊界層直徑,又由于兩極板間的實際距離為1mm,則可計算得到電弧長度的實際值為1mm。
同理還可計算得到電弧半徑及電弧熱邊界層半徑的實際值分別為0.88mm、1.44mm。
通過一系列的量取計算可以得到電弧電弧半徑與電弧電流的關系曲線,如圖6所示。在t=300μs左右時,燃弧趨于穩定,此時電弧半徑為恒定電流時的電弧半徑0.73mm。

圖6 電弧半徑隨電流的變化曲線
2)不同時刻電弧的電流及電壓
由于高速攝像儀是每10μs拍攝一幀,則對應每一張空氣電弧圖片的電流電壓只需在實驗示波器得出的電流電壓實驗波形中讀出,將示波器所得空氣電弧電流電壓實驗波形導入Matlab軟件中,在Matlab軟件中直接讀取每一幀的電流電壓,例圖5所示照片是在260μs時刻的空氣電弧平板實驗的高速攝像,其電流、電壓分別為77.28A、27.2V。
工業化發展階段的變化,意味著經濟發展驅動因素的改變,即從工業化中期的主要依靠資本驅動,向主要依靠技術進步要素驅動轉變。據此,本文對青島市未來的經濟發展提出以下建議:
1)在所研究的空間內氣體處于局部熱力學平衡狀態;
2)本文空氣電弧平板實驗燃弧時間僅僅600μs,且電弧基本靜止不動,故忽略電弧對流散熱,電弧只是以輻射與傳導的方式散熱:

其中U為電弧電壓,Ⅰ為電弧電流,u為輻射熱功率,v為傳導熱功率;
又利用凈輻射系數法求輻射熱,則

其中ε是凈輻射系數,單位為W/(m3*sr-1);r為電弧半徑;L為電弧長度;
另外,熱傳導分為兩個部分:一部分是電弧與觸頭極板間的傳導;由于有上下兩塊極板,則其傳導功率計算式:

其中λ1為觸頭材料熱導率;由于純銅的熔點為1356K,環境溫度為300K,則觸頭極板的溫差ⅠT=1056K,d為觸頭極板厚度;
另一部分是區域1 電弧圓柱體區通過區域2熱邊界層區域向區域3 周圍空氣區域的傳導[11],為了簡化計算,假設氣體的熱導率λ是常數,則電弧由傳導散發的功率可用下式計算:

其中λ2為熱邊界層氣體熱導率,T為電弧溫度,T0為環境溫度,R為電弧熱邊界層半徑。則電弧由傳導所散發的熱功率:

3)極板觸頭材料為T1 銅,純度99.95%,觸頭熱導率近似為常數,λ1=373W/(m*K);
4)電弧近似直圓柱體,設定為恒溫;且在電弧外還有一層厚度為a的熱邊界層,從文獻[6],可以得出熱邊界層的熱導率λ2近似為1W/(m*K);
對于高速攝像的每一幀圖片均滿足熱力學平衡狀態,取圖5為例來計算這一幀圖片的銅絲引弧的空氣電弧凈輻射系數。
根據2.3節的數據處理所得數值帶入到相應的公式中可得:
電流功率:

第一部分的傳導熱功率:

第二部分的傳導熱功率:

聯立式(1)~(5)可得到:

將上述計算式帶入式(6)中解得銅絲引弧空氣電弧在這一時刻的凈輻射系數:

同理可以計算不同時刻銅絲引弧的空氣電弧凈輻射系數。同時對于不同實驗條件下的電流趨于穩定時的凈輻射系數與電流、極板間距的關系,根據表1數據可以推廣得到:1)在電弧電流一定時,極板間距越大,電弧凈輻射系數越?。?)在極板間距一定時,電弧電流越大時,電弧凈輻射系數越小。

表1 穩定電弧凈輻射系數隨電流的變化
然而當電弧半徑為0.88mm 時且電弧溫度未知時,可以從式(6)推導出銅絲引弧的空氣電弧凈輻射系數隨溫度的變化公式:

其中凈輻射系數ε單位為W*cm-3*sr-1。
同理可以推出在同一時刻即電弧半徑相同時,電弧凈輻射系數隨溫度變化的變化關系曲線。
本文利用銅絲引弧的空氣電弧平板實驗可以較為簡單、準確的求解電弧凈輻射系數,在實驗滿足局部熱力學平衡下,利用實驗中得到的電弧相關數據,求解得到電弧半徑為0.88mm 時,銅絲引弧的空氣電弧凈輻射系數為4.4*104W*cm-3*sr-1;在電弧半徑一定,電弧溫度未知時,銅絲引弧的空氣電弧凈輻射系數隨溫度的變化曲線;并推廣得到穩定電弧的凈輻射系數隨電流的變化滿足以下規律:1)在電弧電流一定時,極板間距越大,電弧凈輻射系數越小;2)在極板間距一定時,電弧電流越大時,電弧凈輻射系數越小。
[1]R.W.Libermann,J.J.Lowke.Radiation emission for sulfur hexafluoride arc plasmas[J].Journal of quantitative spectroscopy and radiative transfer,1976,16:23-264.
[2]J.J.Lowke.Predictions of arc temperature using approximate emission coefficients for radiation losses[J].Journal of quantitative spectroscopy and radiative transfer,1974,14:111-122.
[3]P.J.Shayler,M.T.C.Fang.Radiation transport in wall-stabilised nitrogen arcs[J].Journal of Physics D:Appl.Physics.,1978,11:1743-1756.
[4]A.Gleizes,B.Rahmani,J.J.Gonzalez,et al.Calculation of net emission coefficient in N2,SF6and SF6-N2arc plasmas[J].Journal of Physics D:Appl.Physics.,1991,24:1300-1309.
[5]Naghizadeh Kashani,Y.Cressault.Net emission coefficient of air thermal plasmas[J].Journal of Physics D:Appl.Physics.,2002,35(22):2925-2934.
[6]吳翊.低壓空氣電弧多場耦合過程的仿真及實驗研究[D].西安:西安交通大學,2006.
[7]J Menart,S Malik.Net emission coefficients forargoniron thermal plasmas[J].Institute of Physics Publishing,2002:867-874.
[8]Y Naghizadeh-Kashani,Y Cressault,A Gleizes.Net emission coefficient of air thermal plasmas[J].Institute of Physics Publishing,2002:2925-2934.
[9]X Zhang,J Zhang,E Gockenbach.Calculation of pressure and temperature In medium-voltage electrical installations due to fault arcs[J].Journal of Physics D:Applied Physics,2008.
[10]Vladimir Aubrecht,Milada Bartlova.Net Emission Coefficients of Radiation in Air and SF6Thermal Plasmas[J].Plasma Chem Plasma Process,2009(29):131-147.
[11]李新福.低壓電器電弧仿真研究[D].天津:河北工業大學,2004.