劉業(yè)翔, 梁學(xué)民,, 李 劼, 張紅亮, 徐宇杰, 丁鳳其, 鄒 忠

(1. 中南大學(xué) 冶金科學(xué)與工程學(xué)院，長沙 410083；2. 河南中孚實(shí)業(yè)股份有限公司，鞏義 451200)

底部出電型鋁電解槽母線結(jié)構(gòu)與電磁流場(chǎng)仿真優(yōu)化

劉業(yè)翔1, 梁學(xué)民1,2, 李 劼1, 張紅亮1, 徐宇杰1, 丁鳳其1, 鄒 忠1

(1. 中南大學(xué) 冶金科學(xué)與工程學(xué)院，長沙 410083；2. 河南中孚實(shí)業(yè)股份有限公司，鞏義 451200)

從減小槽內(nèi)水平電流和垂直磁感應(yīng)強(qiáng)度進(jìn)而改善磁流體穩(wěn)定性的角度出發(fā)，提出一種底部出電型結(jié)構(gòu)鋁電解槽，該種電解槽采用陰極垂直出電方式代替?zhèn)鹘y(tǒng)水平出電方式。在ANSYS軟件平臺(tái)上，建立400 kA級(jí)該槽型的電磁場(chǎng)模型并進(jìn)行求解，根據(jù)計(jì)算結(jié)果對(duì)母線配置進(jìn)行優(yōu)化，得到一種可使磁場(chǎng)分布最優(yōu)的母線結(jié)構(gòu)，在該母線配置下，垂直磁感應(yīng)強(qiáng)度最大值為1.658 mT，平均值為0.401 mT，遠(yuǎn)低于同規(guī)格普通電解槽，磁流體穩(wěn)定性計(jì)算結(jié)果進(jìn)一步表明該槽能在低極距下穩(wěn)定運(yùn)行，具有較大的節(jié)能潛力。

鋁電解; 母線優(yōu)化; 多物理場(chǎng)

現(xiàn)代鋁電解理論表明，鋁液中水平電流和垂直磁場(chǎng)相互作用引起的電磁力是影響鋁電解槽內(nèi)磁流體穩(wěn)定性的主要因素[1]，而水平電流和垂直磁感應(yīng)強(qiáng)度的大小及分布取決于電解槽及母線的結(jié)構(gòu)設(shè)計(jì)。

在傳統(tǒng)鋁電解槽內(nèi)，電流由陰極底部的水平鋼棒從兩側(cè)導(dǎo)出，再經(jīng)陰極母線匯集至下一槽的立柱母線。母線結(jié)構(gòu)的合理設(shè)計(jì)是實(shí)現(xiàn)槽內(nèi)電磁場(chǎng)優(yōu)化的關(guān)鍵，這一方面國內(nèi)外學(xué)者進(jìn)行了很多研究。早期，TVEDT和NEBELL[2]以及BUIZA[3]分別開發(fā)名為“NEWBUS”和“BUSCAL”的 1D 線單元母線模型，并通過求解電熱耦合方程進(jìn)行母線系統(tǒng)的設(shè)計(jì)與優(yōu)化。KACPRZAK等[4]開發(fā)了電解槽及母線的3D實(shí)體電磁場(chǎng)計(jì)算模型，以調(diào)整和優(yōu)化 100 kA槽陰極軟母線中的電流分布。DUPUIS和BOJAREVICS[5]開發(fā)基于ANSYS軟件的3D實(shí)體母線模型和基于 TECPLOT軟件的 1D槽內(nèi)導(dǎo)體和母線系統(tǒng)模型，進(jìn)行 500 kA槽的母線電流優(yōu)化。同時(shí)，不少學(xué)者從母線優(yōu)化效果的角度開展磁場(chǎng)計(jì)算方法的驗(yàn)證工作[6?8]，國外學(xué)者報(bào)道有關(guān)使用三維結(jié)構(gòu)的母線設(shè)計(jì)或優(yōu)化模型[9?10]，研究 100 kA 和500 kA預(yù)焙槽母線系統(tǒng)的電流均勻分布問題。近些年來，李茂等[11]則應(yīng)用遺傳算法對(duì)母線進(jìn)行優(yōu)化計(jì)算，該方法能得到最優(yōu)化的母線結(jié)構(gòu)，但該方法亦存在優(yōu)化計(jì)算過程時(shí)間較長和母線參數(shù)的設(shè)置復(fù)雜等不足。

對(duì)于底部出電型鋁電解槽的研究，國內(nèi)外目前進(jìn)行得很少，鮮有文獻(xiàn)報(bào)道，更多的是一種概念化電解槽，其中，PETERSON等[12]提出了一種“蘑菇狀”陰極導(dǎo)流型鋁電解槽，采用蘑菇狀可潤濕性陰極，但該槽型的陰極材料抗腐蝕和耐沖擊性能難以達(dá)到工業(yè)化要求；GEORGES和NORA[13]在其專利中提出了一種采用惰性陽極的導(dǎo)流槽，采用聚鋁溝溝型陰極，陰極導(dǎo)桿處于槽底陰極塊的幾何中心處，但也存在陰極易于早期破損的問題；此外，還有不少豎式鋁電解槽，部分結(jié)構(gòu)也采用將電流從底部導(dǎo)出[14]。

本文作者提出一種有別于傳統(tǒng)結(jié)構(gòu)的 400 kA級(jí)底部出電型鋁電解槽，在ANSYS有限元平臺(tái)上，根據(jù)由簡到繁的過程，對(duì)新型槽母線配置進(jìn)行系統(tǒng)優(yōu)化研究，得到最優(yōu)化的母線結(jié)構(gòu)，最后通過研究電?磁?流場(chǎng)和磁流體的穩(wěn)定性，驗(yàn)證該電解槽的巨大節(jié)能潛力。

表1 400 kA預(yù)焙鋁電解槽的主要結(jié)構(gòu)參數(shù)Table 1 Main structure parameters of 400 kA reduction cell

1 底部出電型鋁電解槽結(jié)構(gòu)及母線優(yōu)化

1.1 結(jié)構(gòu)簡介

傳統(tǒng)鋁電解槽的陰極鋼棒為水平放置，電流經(jīng)陰極炭塊在陰極鋼棒匯集后，沿水平方向流出電解槽并連接至陰極母線，如此則產(chǎn)生較大的水平電流，造成電解槽內(nèi)鋁液波動(dòng)，從而降低電解槽的電流效率。為此，本文作者提出一種底部出電型鋁電解槽，即將傳統(tǒng)水平鋼棒改變成為水平?垂直組合的形式，電流經(jīng)過垂直鋼棒從電解槽底部導(dǎo)出，從而大大降低槽內(nèi)的水平電流。

底部出電鋁電解槽除陰極炭塊、鋼棒及槽周圍母線外，其他結(jié)構(gòu)與現(xiàn)行普通400 kA鋁電解槽結(jié)構(gòu)基本一致，其主要結(jié)構(gòu)參數(shù)如表1所列，其陰極鋼棒與母線的連接如圖1所示。

圖1 陰極鋼棒結(jié)構(gòu)示意圖Fig.1 Schematic diagram of structure of cathode bar

1.2 母線結(jié)構(gòu)設(shè)計(jì)

與傳統(tǒng)結(jié)構(gòu)槽相同，母線設(shè)計(jì)對(duì)本新型結(jié)構(gòu)槽內(nèi)電磁場(chǎng)的優(yōu)化至關(guān)重要。為此，本文作者遵循由簡至繁的設(shè)計(jì)理念，從最簡單的連接方式，逐漸引入底部母線補(bǔ)償和端部母線補(bǔ)償，再考慮出電母線間的相互抵消作用，直至得到最優(yōu)化的母線設(shè)計(jì)方案，該過程如圖2所示。在整個(gè)過程中，總共設(shè)計(jì)了幾十種母線結(jié)構(gòu)，記為 SG1?SG31，根據(jù)其磁場(chǎng)分布的特點(diǎn)及大小，對(duì)母線配置進(jìn)行評(píng)價(jià)。由于篇幅所限，本研究僅選出其中有代表性的4種方案，如圖3所示，其中，SG1為根據(jù)最初底部出電思路得到的母線最簡單結(jié)構(gòu)，SG7為具備端部磁場(chǎng)補(bǔ)償?shù)哪妇€結(jié)構(gòu)，SG19為同時(shí)具備端部與底部磁場(chǎng)補(bǔ)償?shù)哪妇€結(jié)構(gòu)，SG31為最終經(jīng)過各種優(yōu)化后的母線結(jié)構(gòu)。

在設(shè)計(jì)過程中，發(fā)現(xiàn)電解槽電場(chǎng)結(jié)果基本類似，因此，本研究中僅列出SG31的電場(chǎng)結(jié)果。此外，在磁場(chǎng)方面，對(duì)槽電壓與電流效率影響最大的為垂直磁場(chǎng)，因此，僅列出各種母線配置下的槽內(nèi)垂直磁場(chǎng)的仿真結(jié)果。

2 電磁流場(chǎng)仿真模型

2.1 電場(chǎng)的計(jì)算原理與模型

由于為電解槽供電的整流電源可以近似認(rèn)為是一個(gè)恒流源，故鋁電解槽內(nèi)電場(chǎng)也可以近似認(rèn)為是一個(gè)穩(wěn)恒電場(chǎng)，電流滿足歐姆定律和守恒方程：

式中：J是電流密度, A/m2；E是電場(chǎng)強(qiáng)度, V/m；σ是電導(dǎo)率, S/m。利用矢量等式?×(φ?)=0，獲得標(biāo)量電勢(shì)φ與電場(chǎng)強(qiáng)度E之間的關(guān)系φ-?=E，即可求解出電場(chǎng)分布[15]。計(jì)算所用有限元網(wǎng)格模型如圖4所示，鋼棒與母線及陰極的連接如圖5所示。

圖2 母線設(shè)計(jì)流程Fig.2 Flowchart for bus bar design

圖3 母線結(jié)構(gòu)的設(shè)計(jì)和演變Fig.3 Design and optimization of bus structure

圖4 3臺(tái)實(shí)體槽及槽周母線電場(chǎng)計(jì)算有限元網(wǎng)格模型Fig.4 Finite element mesh model of three cells and near by bus for calculation of electric field

圖5 陰極鋼棒與陰極炭塊示意圖Fig.5 Schematic diagram of cathode steel bars and cathode carbon blocks

電場(chǎng)邊界條件如下：在電源正極方向上的橫梁母線進(jìn)電位置上施加電流；在電源負(fù)極方向上6個(gè)陽極立柱上施加零電勢(shì)，采用標(biāo)量電位求解電場(chǎng)。

2.2 磁場(chǎng)的計(jì)算原理與模型

鋁電解槽內(nèi)外電流源復(fù)雜，且含有大量的磁介質(zhì)，這使得磁場(chǎng)計(jì)算難度大大增加。電場(chǎng)分布求解合理與否直接影響到磁場(chǎng)的求解精度。本研究采用GP ψ-DP三步標(biāo)量磁位法求解磁場(chǎng)。鋁電解槽的磁場(chǎng)問題滿足穩(wěn)態(tài)麥克斯韋方程組：

式中：H為磁場(chǎng)強(qiáng)度，A/m；J為電流密度，A/m2；B為磁場(chǎng)感應(yīng)強(qiáng)度，T；μ為磁導(dǎo)率，N/A2。

磁場(chǎng)計(jì)算的網(wǎng)格模型如圖6所示。磁場(chǎng)邊界條件如下：鋁電解槽磁場(chǎng)求解屬于開域問題，假設(shè)有限空氣的外表面處于無限遠(yuǎn)處，在外邊界的節(jié)點(diǎn)上施加零磁標(biāo)量位，即Dirichlet邊界條件[16?17]。

2.3 流場(chǎng)及穩(wěn)定性計(jì)算原理與模型

鋁電解槽內(nèi)熔體是不可壓縮粘性流體，其運(yùn)動(dòng)屬于典型的湍流運(yùn)動(dòng)。本研究應(yīng)用非線性淺水模型，通過耦合全槽三維仿真獲得的電磁場(chǎng)分布，對(duì)該槽內(nèi)熔體運(yùn)動(dòng)進(jìn)行動(dòng)力學(xué)分析，并考察該槽不同極距下的磁流體穩(wěn)定性。具體模型及計(jì)算流程見文獻(xiàn)[18]。

圖6 7臺(tái)槽及槽周母線磁場(chǎng)計(jì)算有限元網(wǎng)格模型Fig.6 Finite element mesh model of seven cells and near by bus for calculation of magnetic field

3 仿真計(jì)算結(jié)果與分析

3.1 電場(chǎng)

本研究計(jì)算并分析不同母線配置下400 kA槽底出電槽的電場(chǎng)結(jié)果，發(fā)現(xiàn)在母線設(shè)計(jì)遵循電阻平衡的前提下，電場(chǎng)結(jié)果并未隨母線結(jié)構(gòu)變化而呈現(xiàn)較大變化。故本研究選取最終母線設(shè)計(jì)方案SG31的電場(chǎng)結(jié)果進(jìn)行分析。

以中間槽作為目標(biāo)分析槽，計(jì)算所得的各部分電壓分布列于表2，為進(jìn)行對(duì)比，本研究亦對(duì)傳統(tǒng)400 kA槽進(jìn)行電場(chǎng)計(jì)算，其各部分電壓值如表2所列。

對(duì)比發(fā)現(xiàn)，傳統(tǒng)槽鋁液層壓降有14 mV，而底部出電槽的鋁液層壓降只有4 mV，故相比于傳統(tǒng)槽，該槽陰極壓降也有顯著下降。

3.2 磁場(chǎng)

磁場(chǎng)為母線優(yōu)化的主要判別依據(jù)，故采用三步標(biāo)量磁位法，計(jì)算了母線優(yōu)化(SG1?SG31)過程中所有相應(yīng)母線配置下的磁場(chǎng)分布，在此，僅給出具備代表性的幾種母線配置所對(duì)應(yīng)的磁場(chǎng)計(jì)算結(jié)果，同時(shí)給出傳統(tǒng)結(jié)構(gòu)400 kA槽的磁場(chǎng)計(jì)算結(jié)果。鋁液層空間3個(gè)方向上磁場(chǎng)的最大值、平均值以及對(duì)應(yīng)的母線用鋁量如表3所示。圖7和8所示分別為方案SG31和傳統(tǒng)結(jié)構(gòu)槽中鋁液的磁場(chǎng)分布情況。

磁場(chǎng)計(jì)算結(jié)果表明： SG31的磁場(chǎng)分布十分理想，其中，Bz的最大絕對(duì)值僅為1.658 mT，平均值為0.401 mT，二者都比傳統(tǒng)結(jié)構(gòu)槽的小很多。另外，盡管該方案母線用量比傳統(tǒng)槽增加了約3 t，但相比SG19等方案，用量則大大減少。因此，綜合磁場(chǎng)結(jié)果來權(quán)衡，SG31為一種比較優(yōu)化的方案。

3.3 磁流體穩(wěn)定性

通過母線優(yōu)化得到最優(yōu)化方案SG31，本研究應(yīng)用非線性淺水模型分析該槽在極距分別為5 cm和4 cm情況下的鋁液?電解質(zhì)界面瞬態(tài)波動(dòng)及熔體瞬態(tài)流場(chǎng)。5 cm和4 cm極距下磁流體的穩(wěn)定性計(jì)算結(jié)果分別如圖9和10所示。圖中，ζm/h表示界面平均波動(dòng)量(ζm)與極距(h)的比值，t表示時(shí)間。

從圖9和10可以看出，當(dāng)極距為5 cm時(shí)，鋁液?電解質(zhì)界面波動(dòng)很?。欢?dāng)極距降為4 cm時(shí)，磁流體穩(wěn)定性仍然良好。另外，槽內(nèi)鋁液流場(chǎng)能很快形成穩(wěn)態(tài)，其流動(dòng)形態(tài)為典型的對(duì)稱兩渦流，最大流速和平均流速分別為18 cm/s和4.8 cm/s。從這一計(jì)算結(jié)果看，槽底出電400 kA電解槽能在很低的極距下穩(wěn)定運(yùn)行，能大幅降低槽電壓，節(jié)能潛力巨大。

表2 傳統(tǒng)槽與研究槽的電場(chǎng)結(jié)果對(duì)比Table 2 Comparison of electric field results of conventional and proposed cells

表3 不同母線配置下磁場(chǎng)分量的比較Table 3 Comparison of magnetic results of different bus structures

圖7 SG31的鋁液層磁感應(yīng)強(qiáng)度分布Fig.7 Magnetic flux density distribution of aluminium layer of SG31: (a) Bx; (b) By; (c) Bz

圖8 傳統(tǒng)400 kA槽的鋁液層磁感應(yīng)強(qiáng)度分布Fig.8 Magnetic flux density distribution of aluminium layer of conventional 400 kA cell: (a) Bx; (b) By; (c) Bz

圖9 極距為5 cm時(shí)鋁液?電解質(zhì)界面波動(dòng)情況Fig.9 Aluminum-electrolyte interface wave curve at anode cathode distance of 5 cm

圖10 極距為4 cm時(shí)鋁液?電解質(zhì)界面波動(dòng)情況Fig.10 Aluminum-electrolyte interface wave curve at anode cathode distance of 4 cm

4 結(jié)論

1) 從減小槽內(nèi)水平電流和垂直磁場(chǎng)進(jìn)而改善磁流體穩(wěn)定性的角度出發(fā)，提出一種底部出電型結(jié)構(gòu)鋁電解槽，該種電解槽采用陰極垂直出電的方式代替?zhèn)鹘y(tǒng)的水平出電方式。

2) 對(duì)底部出電型鋁電解槽進(jìn)行了母線設(shè)計(jì)，通過物理場(chǎng)優(yōu)化，得到一種最佳的母線配置設(shè)計(jì)方案。

3) 在最佳的母線配置方案下，垂直磁場(chǎng) Bz絕對(duì)值的最大值為1.658 mT，平均值為0.401 mT，該最大值和平均值比水平出電槽的小很多。

4) 磁流體穩(wěn)定性計(jì)算也顯示，該電解槽在低極距下還能維持較高的穩(wěn)定性，具備較大的節(jié)能空間。

REFERENCES

[1] LINDSAY R I, DAVIDSON P A. Stability of interfacial waves in aluminum reduction cell[J]. J Fluid Mech, 1998, 362: 327?331.

[2] TVEDT T, NEBELL H G. Newbus, a simulation program for calculation of the current distribution in the bus bar system of alumina reduction cells[C]//BOXALL L G. Light Metals.Phoenix, AZ: Metallurgical Soc of AIME, 1988: 567?573.

[3] BUIZA J I. Electromagnetic optimization of the V-350 cell[C]//CAMPBELL P G. Light Metals. Las Vegas, NV:Metallurgical Soc of AIME, 1989: 211?214.

[4] KACPRZAK D, GUSTAFSSON M, LI L, TAYLOR M.Numerical analysis of the collector bar current distribution of a reduction cell[C]//GALLOWAY T J. Light Metals. San Antonio,TX: TMS, 2006: 367?369.

[5] DUPUIS M, BOJAREVICS V. Weakly coupled thermo-electric and MHD mathematical models of an aluminium electrolysis cell[C]//KVANDE H. Light Metals. San Francisco, CA: TMS,2005: 449?454.

[6] BLANC J M, ENTNER P. Application of computer calculations to improve electromagnetic behaviour of pots[C]//MCMINN C J.Light Metals. Las Vegas, NV: Metallurgical Soc of AIME, 1980:285?295.

[7] SEVERO D S, SCHNEIDER A F, PINTO E C V, GUSBERTI V,POTOCNIK V. Modeling magnetohydrodynamics of aluminum electrolysis cells with ANSYS and CFX[C]//KVANDE H. Light Metals. San Francisco, CA: TMS, 2005: 475?480.

[8] 周 萍. 鋁電解槽內(nèi)電磁流動(dòng)模型及鋁液流動(dòng)數(shù)值仿真的研究[D]. 長沙: 中南大學(xué), 2002: 76?85.ZHOU Ping. A research on mathematical models of electromagnetic hydrodynamics and numerical simulation of metal pad flow in aluminum reduction cells[D]. Changsha:Central South University, 2002: 76?85.

[9] 吳建康, 黃 珉, 黃 俊, 姚世煥. 鋁電解槽電解質(zhì)?鋁液流動(dòng)及鋁液表面變形計(jì)算[J]. 中國有色金屬學(xué)報(bào), 2003, 13(1):241?244.WU Jian-kang, HUANG Min, HUANG Jun, YAO Shi-huan.Computation of flow field of electrolyte-aluminium liquid and surface distortion of aluminum liquid in reduction cell[J]. The Chinese Journal of Nonferrous Metals, 2003, 13(1): 241?244.

[10] LI Jie, LIU Wei, LAI Yan-qing, WANG Zhi-gang, LIU Ye-xiang.Analysis of cathode voltage drop in aluminum electrolysis cells with an electric contact model[C]//S?RLIE M. Light Metals.Orlando, FL: TMS, 2007: 465?470.

[11] LI Mao, ZHOU Jie-min. Modeling and optimization of busbar configuration in aluminum electrolysis cells with genetic algorithm[C]//HAGNI A M. Light Metals. Seattle, WA: TMS,2010: 489?492.

[12] PETERSON R D, RICHARDS N E, TABEREAUX A T. Results of 100-hour electrolysis test of a cermet anode: operational results and industry perspective [C]//BLCKERT C M. Light Metals. Warrendale, LA: TMS, 1990: 385?393.

[13] GEORGES B, NORA V D. Aluminum production cell and cathode. US, 6358393[P]. 2002?03?19.

[14] NORA V D. Cell for aluminum electrowinning. US 6093304[P].2000?07?25.

[15] 劉 偉, 李 劼, 賴延清. 鋁電解槽電磁流場(chǎng)數(shù)學(xué)建模及應(yīng)用[J]. 中國有色金屬學(xué)報(bào), 2008, 18(5): 909?916.LIU Wei, LI Jie, LAI Yan-qing. Development and application of electro-magneto-flow mathematic model of aluminum reduction cells[J]. The Chinese Journal of Nonferrous Metals, 2008, 18(5):909?916.

[16] LI Jie, LIU Wei, LAI Yan-qing, LIU Ye-xiang. An improved finite element model for electro-magnetic analysis in aluminum cells[J]. JOM, 2008, 60(2): 58?61.

[17] LI Jie, LIU Wei, LAI Yan-qing, LI Qin-yu, LIU Ye-xiang.Coupled simulation of 3D electro-magneto-flow field in hall-heroult cells using finite element method[J]. Acta Metallurgica Sinica: English Letters, 2006, 19(2): 105?116.

[18] 徐宇杰. 鋁電解槽內(nèi)熔體運(yùn)動(dòng)數(shù)學(xué)建模及應(yīng)用研究[D]. 長沙:中南大學(xué), 2010.XU Yu-jie. A study of mathematical modeling and its application for the melts motion in aluminium reduction cells[D]. Changsha:Central South University, 2010.

Simulation and optimization of bus structure and
electro-magneto-flow field of aluminum reduction cells with vertical bottom bars

LIU Ye-xiang1, LIANG Xue-min1,2, LI Jie1, ZHANG Hong-liang1, XU Yu-jie1, DING Feng-qi1, ZHOU Zhong1
(1. School of Metallurgical Science and Engineering, Central South University, Changsha 410083, China;2. Henan Zhongfu Industrial Co., Ltd., Gongyi 451200, China)

Based on the target of reducing the horizontal current and vertical magnetic field strength, thereby improving the magnetohydrodynamic stability, a kind of aluminium reduction cell with vertical bottom was improved, in which the current flowed out of cathode vertically instead of horizontally. An electromagnetic simulated model of a 400 kA cell with this structure was developed by ANSYS software, and a bus bar scheme that can make the best electromagnetic field distribution was obtained after being optimized according to the computing result. Adopting this bus bar scheme, the maximum and average values of vertical magnetic flux density are 1.658 mT and 0.401 mT, respectively, which are far lower than those in traditional cells, and the MHD computation further proves that the above cell can operate stably under a small anode-cathode distance and has great potential in energy saving.

aluminium electrolysis; bus bar optimization; multi-physical field

TF821

A

1004-0609(2011)07-1688-08

國家“十一五”科技支撐計(jì)劃資助項(xiàng)目(2009BAE85B00)；國家自然科學(xué)基金資助項(xiàng)目(50874020)

2011-01-12；

2011-05-20

李 劼，教授，博士；電話：0731-88876454；E-mail：csulightmatels@126.com

(編輯 李艷紅)