在粉體生產(chǎn)行業(yè),陰極板清洗常采用的技術(shù)手段有超聲波清洗、水射流清洗等方式[1,2]。在鎳粉電解車間中,依然主要依靠工人雙手持板刷進(jìn)行刷粉作業(yè)[3],人力成本高,勞動(dòng)強(qiáng)度大,且危害工人健康,用機(jī)械刷粉方式取代人工刷粉就成為一種必然的發(fā)展趨勢(shì)。然而,當(dāng)陰極板鎳粉自動(dòng)清刷機(jī)(以下簡稱“刷粉機(jī)”)的刷板進(jìn)入電解液時(shí),勢(shì)必會(huì)引起電解液中電流和電壓的變化[4],增加電耗的同時(shí)還會(huì)影響鎳粉的產(chǎn)量和品質(zhì)。因此,為縮短刷粉機(jī)刷板的研發(fā)周期,就需要探究刷粉過程中刷板結(jié)構(gòu)對(duì)電解槽內(nèi)電流密度分布的影響規(guī)律。
本文從電化學(xué)理論的角度出發(fā),用電流密度分布理論建立了刷粉機(jī)運(yùn)行條件下的1/4電解槽數(shù)學(xué)模型,然后根據(jù)電解生產(chǎn)工藝參數(shù)設(shè)置邊界條件[5],使用多物理場仿真軟件COMSOL的電化學(xué)模塊對(duì)電解過程中電解槽內(nèi)電流密度分布及變化情況進(jìn)行了數(shù)值仿真,得到了A、B和C三種型號(hào)刷板工作過程中陽極表面無量綱電流密度分布云圖、電解液中電位分布等值面圖和電流變化曲線。最后,通過與無刷板情況下的結(jié)果進(jìn)行比較,得到了對(duì)電解液電流場影響最小的刷板結(jié)構(gòu)設(shè)計(jì)方案并進(jìn)行了實(shí)驗(yàn)驗(yàn)證。
如圖1所示,為刷粉機(jī)運(yùn)行下的1/4電解槽示意圖。電解槽主要由陰極板、陽極棒以及充斥于陰極與陽極之間的電解液(不斷循環(huán)補(bǔ)充中)組成,刷板與陰極板接觸部分為刷板的刷絲區(qū)。刷粉機(jī)兩刷板的運(yùn)動(dòng)方式為沿著陰極板與陽極之間的間隙上下同時(shí)運(yùn)動(dòng),通過刷絲區(qū)與陰極板左右兩側(cè)面的摩擦使附著于陰極板表面的鎳粉脫落并沉降在電解槽底部。
圖1 刷粉機(jī)運(yùn)行下的1/4電解槽示意圖
對(duì)于刷板結(jié)構(gòu),在保證其強(qiáng)度和剛度的前提下,且考慮到刷板絕緣性要求,根據(jù)設(shè)計(jì)先簡單后復(fù)雜的一般規(guī)律[5],首先設(shè)計(jì)的是非金屬材料板材加工制成的A型結(jié)構(gòu)形式刷板;為減弱刷板對(duì)陰陽極間電流場的阻礙作用,在A型刷板的基礎(chǔ)上在刷板中間開槽設(shè)計(jì)了B型刷板;最后,考慮到陽極的形狀及其位置分布規(guī)律,在刷板上開設(shè)與陽極形狀及位置分布規(guī)律一致的回形槽,設(shè)計(jì)了C型刷板。A、B和C三種刷板的結(jié)構(gòu)形式如圖2所示。
圖2 三種刷板的結(jié)構(gòu)形式
在本文研究的問題中,在刷板上下運(yùn)行時(shí)電解槽中陽極表面的電流密度分布受電極的尺寸、刷板與陰陽極的相對(duì)位置、電解液的電導(dǎo)率、質(zhì)量傳遞速率以及離子濃度等因素的影響,均包含于電流密度分布理論的研究范疇。
電解質(zhì)的凈電流密度用所有離子的通量總和來表示:
且:
其中,為電解質(zhì)的電流密度矢量(A/m2);F為法拉第常數(shù)(C/mol);Ni為物質(zhì)i的通量(mol/(m2·s));zi為離子電荷數(shù)。電解液的離子通量可以用Nernst-Planck方程描述:
該方程利用三個(gè)添加項(xiàng)(擴(kuò)散、遷移和對(duì)流)解釋了電解液中溶質(zhì)物質(zhì)的通量,ci為離子i的濃度(mol/m3);Di為擴(kuò)散系數(shù)(m2/s);um,i為遷移率(s·mol/kg);為電解質(zhì)電位(V);u為速度矢量(m/s)。
在電解生產(chǎn)中,電解質(zhì)溶液處于不斷地補(bǔ)充與濃度調(diào)節(jié)中,可以認(rèn)為電解液不存在濃度梯度,可以忽略式(2)中的擴(kuò)散項(xiàng);且電極-電解質(zhì)界面處的電化學(xué)反應(yīng)進(jìn)行得非???,由于電極極化引起的極化電阻與電解液電阻比起來非常小,即暫且可以忽略電極動(dòng)力學(xué)的影響,認(rèn)為電極-電解質(zhì)界面處的電位差不會(huì)偏離平衡值,即假設(shè)電解質(zhì)為電中性,所以可以忽略式(2)中的對(duì)流項(xiàng)。因此,將式(2)代入式(1)并整理得:
且:
隨后,將zi、um,i、ci和法拉第常量F等用一個(gè)物理量—電導(dǎo)率(1/(W?m2))等效代替得:
即電解質(zhì)中的電流密度矢量:
這個(gè)方程就是電解槽內(nèi)電解質(zhì)的歐姆定律,其中電導(dǎo)率=14.6[S/m]。因此,可以用一次電流密度分布理論來近似描述鎳電解槽中電壓與電流的關(guān)系[6]。為了更加準(zhǔn)確地描述電解槽中電流密度分布情況,使用線性Butler-Volmer方程[7,8](考慮電解液歐姆電阻的同時(shí)也考慮了電極-電解質(zhì)界面處的極化電阻)來描述電解液中的電流分布規(guī)律。鑒于電流在電極表面上的不均勻分布,與電流在通過陰-陽極間電解液時(shí)所遇到的阻力不同有關(guān),這首先是取決于幾何因素[9]。于是,刷板在電解槽中運(yùn)行時(shí)的陽極表面電流密度:
陰極表達(dá)式與式(6)同。其中,i0為交換電流密度(A/m2);S0為無刷板時(shí)陰陽極之間的法向通流面積(m2);Sj(x)-j型刷板進(jìn)入電解液x(m)深度時(shí)陰陽極之間的法向通流面積(m2),j為刷板類型。不同結(jié)構(gòu)形式刷板的表達(dá)式Sj(x)有所區(qū)別,以B型刷板為例:
1)0≤x≤h時(shí):
其中,h為刷板刷絲區(qū)高度(m);D為刷板寬度(m);d為回形槽寬度(m);hc為陰極板沒入電解液深度(m);H為刷板進(jìn)入電解液的最大深度(m)。
另外,R為氣體常數(shù)(J/(K?mol));T為熱力學(xué)溫度(K);αa、αc為陰陽極對(duì)稱因子;η為超電勢(shì)(V),陰極和陽極的超電勢(shì):
φl,c、φl,a為與陰、陽極相鄰電解質(zhì)電位;槽電壓Eslot=12.6[V];陽極平衡電位Eeq,a=0.24[V];陰極平衡電位Eeq,c=0[V]。
刷粉機(jī)刷板結(jié)構(gòu)對(duì)電解槽中電流密度分布規(guī)律影響的研究是刷粉機(jī)研發(fā)過程中的核心技術(shù)問題之一。如圖3所示,為使用COMSOL建立的刷板運(yùn)用條件和無刷板條件下的1/4電解槽三維模型。圖3中深色部分為刷板浸入電解液的表面,本文以刷板進(jìn)入電解液深度200mm位置時(shí)刻的模型為例。
圖3 刷板運(yùn)用條件及無刷板條件下的三維模型
仿真過程中,固定陰陽極間距(30mm),反應(yīng)溫度T=368K,忽略電解過程中電解液溫差變化及其引起的能量耗散[5]。陰陽極的交換電流密度均為i0=300A/m2,陰陽極對(duì)稱因子均為0.5,電解槽電壓12.6V,從而建立了陰陽極的邊界條件。為滿足電極-電解質(zhì)界面處的特殊計(jì)算需求,在陰陽極表面處設(shè)置邊界層網(wǎng)格,以增加其網(wǎng)格單元密度。為模型中的陽極表面添加“積分”運(yùn)算和“平均”運(yùn)算耦合關(guān)系,通過計(jì)算陽極表面電流密度的絕對(duì)值和平均值的比值,可以得到陽極表面的無量綱電流密度分布情況。如圖4所示,為無刷板情況下與在A型、B型和C型刷板進(jìn)入電解液200mm位置時(shí)的陽極表面無量綱電流密度分布仿真結(jié)果。
不均勻的電流密度分布,會(huì)使電極表面的局部反應(yīng)速率處于失控狀態(tài),由此產(chǎn)生的副反應(yīng),會(huì)降低電流效率和產(chǎn)品的質(zhì)量[10]。由圖4可知,在靠近陰極板一側(cè)且無刷板遮擋處陽極表面的電流密度較高。一般而言,無量綱電流密度的最大值與最小值的差值越小,說明電流密度分布相對(duì)越均勻,陽極的溶解也會(huì)越均勻,對(duì)電解反應(yīng)的影響也越小。
圖4 陽極表面無量綱電流密度分布情況
表1 陽極表面無量綱電流密度分布差值對(duì)比
由表1可知,三種刷板中,相較于A型和B型刷板,C型刷板下陽極表面的無量綱電流密度分布差值與無刷板情況下的結(jié)果差別最小。
三種型號(hào)刷板處于電解液中同一位置時(shí)刻下的電位分布情況如圖5所示。
圖5 電解液中電位分布等值面圖
在各種工業(yè)電化學(xué)過程中,電位分布具有十分重要的意義,除了電解加工,都要求具有均勻的電位分布[10]。由圖5可知,刷板所處位置隔斷了陰陽極之間的電位分布。相對(duì)而言,A型刷板對(duì)電解液中電位分布的阻隔作用最嚴(yán)重,B型和C型刷板對(duì)電解液中電位分布的影響較小,其中C型刷板的電位分布情況最接近無刷板條件下的電位分布。
在刷粉機(jī)運(yùn)行過程中,不同結(jié)構(gòu)形式的刷板所引起的槽電壓波動(dòng)規(guī)律有所區(qū)別。為方便研究,把刷板向上和向下運(yùn)行各一次作為一個(gè)刷粉周期,并假設(shè)刷粉機(jī)刷板上行和下行均是勻速狀態(tài),則在電解生產(chǎn)穩(wěn)定進(jìn)行且刷粉機(jī)平穩(wěn)運(yùn)行狀態(tài)下,同一結(jié)構(gòu)形式刷板在電解液中的各個(gè)時(shí)間節(jié)點(diǎn)對(duì)應(yīng)的刷板位置和槽電壓一定。通過實(shí)驗(yàn)測得了三種型號(hào)刷板在刷粉過程中槽電壓隨時(shí)間的變化規(guī)律作為仿真計(jì)算的依據(jù)。為便于計(jì)算,對(duì)于同一刷板,將刷粉過程中刷板進(jìn)入電解液的深度和時(shí)間均分為16等份。于是,在刷板進(jìn)入電解液中的深度和對(duì)應(yīng)電壓值已知的條件下可以計(jì)算出對(duì)應(yīng)的電流值,將所有數(shù)據(jù)整理之后就可以得到每種刷板在一個(gè)刷粉周期中電流的變化規(guī)律。實(shí)驗(yàn)測得的A、B和C三種型號(hào)刷板在一個(gè)刷粉周期中的槽電壓波動(dòng)規(guī)律如圖6所示。
圖6 三種型號(hào)刷板下的槽電壓波動(dòng)規(guī)律
由圖6可知,在刷粉機(jī)運(yùn)行過程中,A型刷板的槽電壓波動(dòng)幅度最大,為2.8V;B型和C型刷板的槽電壓波動(dòng)幅度相對(duì)較小且接近,均為1.3V。
A型、B型和C型刷板運(yùn)用條件下的1/4電解槽的電流波動(dòng)規(guī)律仿真結(jié)果如圖7所示。刷粉過程中的陰陽極間法向通流面積波動(dòng)規(guī)律如圖8所示。
由圖7可見,對(duì)于A型刷板,電流波動(dòng)幅度為1370A;對(duì)于B型刷板,在8.8s電流的下降速度減緩,9.8s降到了829.9A,9.8s~10.2s電流有所回升,10.2s~10.8s電流維持在平穩(wěn)值876A左右,10.8s~11.2s電流開始下降,11.2s時(shí)電流降低至最小值826.8A,11.2s~12.2s電流回升,12.2s時(shí)電流回升速度加快,直至刷板完全脫離電解液,電流回升至正常水平1824.2A左右,電流波動(dòng)幅度為997A;對(duì)于C型刷板,刷板運(yùn)動(dòng)過程中電流波動(dòng)規(guī)律與B型刷板情況下的波動(dòng)規(guī)律類似,只是C型刷板電流波動(dòng)幅度為699A。由此可見,B型和C型刷板運(yùn)用條件下的電流波動(dòng)幅度明顯低于A型刷板,說明刷板開槽可以大幅降低對(duì)電流場的影響。
圖7 三種型號(hào)刷板刷粉過程中電流波動(dòng)規(guī)律仿真結(jié)果
圖8 刷粉過程中陰陽極間法向通流面積波動(dòng)規(guī)律
另外,由圖6和圖8可知,陰陽極間法向通流面積變化幅度決定了刷粉過程中槽電壓變化幅度。而且,由圖7和圖8可知,B型板和C型板陰陽極間法向通流面積變化幅度接近,但C型刷板情況下電流波動(dòng)幅值為B型刷板的70%。由此可見,在陰陽極間法向通流面積接近的情況下,不同結(jié)構(gòu)形式的刷板對(duì)電流場的影響程度不同。
綜合仿真分析結(jié)果可知,開槽方式與陽極形狀及位置分布規(guī)律一致的刷板結(jié)構(gòu)方案(C型板)在同等條件下對(duì)電解槽中電流場的影響最小。
在本次實(shí)驗(yàn)中用到的主要設(shè)備有Multi System 8050數(shù)據(jù)采集器(精度±0.1%)、CSF2-2kA電流傳感器(精度0.2%,測量范圍0~3000A)、磁屏蔽線和三種類型刷板等。其中,電流傳感器將檢測到的大電流轉(zhuǎn)換成0~20mA電流信號(hào)(即1mA代表實(shí)際的100A),再將其傳輸至Multi System 8050數(shù)據(jù)采集器的10號(hào)通道并實(shí)時(shí)記錄;黑(負(fù))、紅(正)兩電壓引線將陰陽極間的電壓信號(hào)傳輸至數(shù)據(jù)采集器的12號(hào)通道并實(shí)時(shí)記錄。圖9為本實(shí)驗(yàn)中的某新型刷粉機(jī)實(shí)驗(yàn)臺(tái),電解槽位于刷粉機(jī)的正下方,同一電解槽中共有四塊陰極板,四條支路為并聯(lián)關(guān)系。圖10為實(shí)驗(yàn)中的主要設(shè)備,電流傳感器安裝于第二塊陰極板所連接的支路母排上。
圖9 刷粉機(jī)實(shí)驗(yàn)臺(tái)
圖10 主要實(shí)驗(yàn)設(shè)備
在保證生產(chǎn)工藝參數(shù)與仿真參數(shù)一致的條件下,使用Multi System 8050數(shù)據(jù)采集器分別記錄無刷板條件下和三種類型刷板刷粉條件下1/4電解槽的槽電壓和電流變化規(guī)律如圖11所示。其中,A型刷板下的電壓電流變化規(guī)律是對(duì)安裝電流傳感器的第二塊陰極板重復(fù)刷4次測得的結(jié)果,B型板和C型板下的變化規(guī)律是對(duì)同一電解槽中四塊陰極板依次進(jìn)行自動(dòng)刷粉作業(yè)(一個(gè)刷粉工作循環(huán))時(shí)測得的結(jié)果。由圖11可知,在刷粉機(jī)對(duì)電解槽的四塊陰極板依次進(jìn)行刷粉作業(yè)時(shí),刷第二塊陰極板時(shí)電流降低的幅值與刷另外三塊陰極板時(shí)電流上升幅值的和幾乎相等??梢哉J(rèn)為在一個(gè)工作循環(huán)(刷粉機(jī)刷四塊陰極板)中,同一塊陰極板電流上升和下降同一幅值的時(shí)間相等,即可等效為刷粉過程中電流不變,只有電壓升高。對(duì)于一個(gè)工作循環(huán),一共有4個(gè)刷粉周期的時(shí)間槽電壓處于高于正常值狀態(tài),刷粉機(jī)在刷粉時(shí)會(huì)產(chǎn)生附加電耗。
圖11 四種情況下槽電壓和電流變化規(guī)律
將實(shí)驗(yàn)數(shù)據(jù)和仿真數(shù)據(jù)整理之后,如圖12所示。
圖12 實(shí)驗(yàn)和仿真數(shù)據(jù)的對(duì)比結(jié)果
由圖12可知,A、B和C三種類型刷板下的仿真數(shù)據(jù)相對(duì)于實(shí)驗(yàn)數(shù)據(jù)的平均偏差分別為4.97%、3.80%和2.13%。在研究過程中,忽略了一些非主要因素,存在模型誤差;測量儀器本身也有測量誤差,實(shí)驗(yàn)過程中電解設(shè)備及刷粉機(jī)實(shí)驗(yàn)臺(tái)運(yùn)行中存在不穩(wěn)定性,都有可能是誤差的來源。實(shí)驗(yàn)結(jié)果表明,刷粉過程中A型、B型和C型刷板下的電流變化規(guī)律與仿真結(jié)果一致,且C型刷板在同等條件下對(duì)電解槽中電流場的影響最小,可以認(rèn)為仿真結(jié)果與實(shí)驗(yàn)結(jié)果吻合,由此驗(yàn)證了數(shù)學(xué)模型的正確性和仿真過程的可靠性。
本文針對(duì)刷粉機(jī)核心部件刷板及其應(yīng)用場所電解槽,將電化學(xué)理論方法應(yīng)用于機(jī)械結(jié)構(gòu)設(shè)計(jì)中,建立了包含刷板的電解槽電流密度分布數(shù)學(xué)模型,使用COMSOL分別對(duì)三種類型刷板情況下的電解槽中電流密度分布情況進(jìn)行了數(shù)值仿真,并對(duì)仿真結(jié)果進(jìn)行了試驗(yàn)驗(yàn)證。主要結(jié)論如下:
1)刷粉過程中,B型刷板和C型刷板對(duì)電解液中電流場的影響明顯弱于A型刷板,即在刷板上開設(shè)回形槽有利于降低刷板對(duì)電解液中電流場的影響。
2)刷板運(yùn)行時(shí)陰陽極之間的法向通流面積變化幅度決定了槽電壓波動(dòng)幅度,且幅值與刷板回形槽尺寸有關(guān),可以通過優(yōu)化回形槽尺寸的方式降低電壓波動(dòng)幅值,降低附加電耗。
3)在保證強(qiáng)度和剛度的條件下,在刷板上采取與陽極形狀及其位置分布規(guī)律相近的開槽方式對(duì)電解槽中電流密度分布規(guī)律的影響最小。
采用數(shù)學(xué)建模和COMSOL仿真的方法可以減少實(shí)驗(yàn)次數(shù)降低實(shí)驗(yàn)成本,為后續(xù)設(shè)備刷板結(jié)構(gòu)的優(yōu)化提供可靠的設(shè)計(jì)依據(jù)和理論參考,對(duì)于降低因刷粉機(jī)而產(chǎn)生的附加電耗和保證鎳粉品質(zhì)具有重要意義。
參考文獻(xiàn):
[1]吳敏,張園,劉海賓.立式電解金屬錳極板清洗機(jī)設(shè)計(jì)[J].中國錳業(yè),2011,29(04):36-38.
[2]吳連杰,楊玉龍,康永欣,潘國華,王琢.基于超聲波的自動(dòng)清洗機(jī)設(shè)計(jì)[J].林業(yè)機(jī)械與木工設(shè)備,2013,41(05):43-45.
[3]孫磊,李軍義,郭順,羅文,楊國啟.超細(xì)鎳粉工藝現(xiàn)狀與發(fā)展趨勢(shì)[J].湖南有色金屬,2014,30(05):57-60.
[4]黨小慶,楊春方,王迪,袁勝利,馬廣大.電除塵器收塵極板表面電流密度分布實(shí)驗(yàn)研究[J].重型機(jī)械,2005,44(02):32-35.
[5]王穎,康敏.基于COMSOL的電噴鍍陽極噴嘴設(shè)計(jì)研究[J].中國機(jī)械工程,2014,25(09):1180-1185.
[6]劉曉坤,王建軍,邵樹淵,等.電解池中電位及電流分布的邊界元分析[A].第三屆航空航天專業(yè)學(xué)術(shù)年會(huì)論文集[C].1988:11.
[7]Eduardo Laborda. A comparison of Marcus-Hush vs. Butler-Volmer electrode kinetics using potential pulse voltammetric techniques[J].Journal of Electroanalytical Chemistry,2011,660(1):169-177.
[8]Yoshinao Okajima.A phase-field model for electrode reactions with Butler-Volmer kinetics[J].Computational Materials Science,2010,50(1):118-124.
[9]黃慧民,曾振歐.應(yīng)用電化學(xué)[M].廣州:華南理工大學(xué)出版社,1995:249-253.
[10]李紹芬.反應(yīng)工程[M].北京:化學(xué)工業(yè)出版社,2015.