張友，王樹博，齊亮，姚克儉

（1浙江工業(yè)大學(xué)化工學(xué)院綠色化學(xué)合成技術(shù)國家重點(diǎn)實(shí)驗(yàn)室培育基地，浙江 杭州 310032；2清華大學(xué)核能與新能源技術(shù)研究院，北京 100084）

引言

全釩液流電池儲能技術(shù)是一種新型的高效電化學(xué)儲能技術(shù)，與其他儲能電池相比，它具有如下優(yōu)勢：能量轉(zhuǎn)換效率高、蓄電容量大、系統(tǒng)設(shè)計靈活、可靠性高、可深度充放電等一系列優(yōu)點(diǎn)，被認(rèn)為是最有前景的儲能技術(shù)[1-7]。

全釩液流電池主要由電極、電解液、質(zhì)子交換膜和集流體等幾個關(guān)鍵部分構(gòu)成，隨著近年來關(guān)鍵材料領(lǐng)域不斷取得技術(shù)突破[8-14]，目前的研究熱點(diǎn)逐漸集中到了電池結(jié)構(gòu)優(yōu)化設(shè)計及對電池內(nèi)部傳遞過程的研究[15-17]。理想情況下，質(zhì)子交換膜既能為正負(fù)電解液提供質(zhì)子傳遞的通道，又能阻隔正負(fù)電極電解液的相互混合[18]。事實(shí)上，在電池運(yùn)行過程中，電解質(zhì)跨膜滲透問題無法完全避免，正負(fù)電解質(zhì)的混合導(dǎo)致的電池自放電現(xiàn)象，是引起的電池容量損失進(jìn)而導(dǎo)致電池性能惡化的重要原因[19-20]，所以無論是在建立電池模型還是研究電池極化以及電池性能壽命等課題上，電池隔膜中的離子傳遞問題都無法回避。關(guān)于電解質(zhì)中釩離子透膜傳質(zhì)問題，前人已經(jīng)有過相關(guān)研究[21-24]。但僅限于考慮質(zhì)子交換膜中單純的濃差擴(kuò)散過程，往往通過直接測得表觀擴(kuò)散系數(shù)來研究釩離子跨膜傳質(zhì)過程，建立的相關(guān)模型也主要是基于 Fick定律離子擴(kuò)散模型推導(dǎo)得到，這并不能完全與電池在運(yùn)行中的實(shí)際狀況符合，難以有實(shí)際應(yīng)用價值。因?yàn)?，釩離子在質(zhì)子交換膜中的傳質(zhì)過程，是處于電場環(huán)境之下，由濃差推動力和電場力協(xié)同作用的結(jié)果，所以除了由于膜兩側(cè)濃度不同而造成的濃差擴(kuò)散外，還應(yīng)該考慮釩離子在場下的電遷移作用。研究這種多物理場耦合作用下離子傳質(zhì)過程對于提高電池性能、優(yōu)化電池結(jié)構(gòu)、提高電池容量和充放電特性以及循環(huán)壽命，具有非常重要的意義，也為掃清全釩液流電池商業(yè)化應(yīng)用的技術(shù)障礙提供了有益的理論基礎(chǔ)。本文正是從協(xié)同場這一角度，通過設(shè)計外加電場作用下釩離子跨膜滲透實(shí)驗(yàn)，系統(tǒng)研究了不同工況下釩離子在質(zhì)子交換膜中的傳質(zhì)過程，并測得了四價釩離子的表觀電遷移率。

1 數(shù)學(xué)模型

對于釩離子在Nafion膜中的傳質(zhì)過程，如果只考慮單純的濃差擴(kuò)散，由Fick定律可知擴(kuò)散傳質(zhì)速率j(d)

在假設(shè)膜左側(cè)濃度c1不變的情況下，通過積分處理可得擴(kuò)散系數(shù)D

在電場作用下，釩離子的電遷移傳質(zhì)速率為j(e)，其數(shù)學(xué)形式為

綜合膜中考慮擴(kuò)散和電遷移這兩種傳質(zhì)形式，總的傳質(zhì)速率公式可以表示為

式中，c1和c2是膜兩側(cè)的濃度；c為質(zhì)子交換膜中釩離子的平均濃度。

因?yàn)槟さ暮穸确浅１。梢哉J(rèn)為膜中從高濃度一側(cè)到低濃度一側(cè)的濃度變化是近似線性的，則式(4)中

模型中，若不考慮由于少量水跨膜傳質(zhì)導(dǎo)致的電解液體積變化，根據(jù)釩離子的物質(zhì)守恒，有

通過以上公式，可以求得在質(zhì)子交換膜中釩離子的電遷移率

為了衡量電場作用對釩離子跨膜傳質(zhì)通量的影響，引入電場因子St的概念，其定義為電場存在下的傳質(zhì)通量與單純濃差擴(kuò)散傳質(zhì)通量之比

2 實(shí) 驗(yàn)

2.1 儀器及藥品

儀器：高精度電池性能測試系統(tǒng)（深圳市新威爾電子有限公司）、紫外可見分光光度計（島津UV-1750）、自動電位滴定儀（上海雷磁）、多功能磁力攪拌器、恒溫水浴鍋、蠕動泵、石墨電極、電熱鼓風(fēng)干燥箱、自制電解槽、Nafion117膜。

藥品：硫酸氧釩（VOSO4·3H2O）(沈陽市海中天精細(xì)化工廠)、硫酸（北京化工廠）、無水硫酸鎂（國藥集團(tuán)化學(xué)試劑有限公司）。

2.2 實(shí)驗(yàn)方案

膜的預(yù)處理：在每一組實(shí)驗(yàn)前將Nafion117質(zhì)子交換膜在質(zhì)量分?jǐn)?shù)為 5%的雙氧水中 80℃煮沸1 h，用去離子水沖洗3次后，放入去離子水中80℃下煮沸1 h，再放入0.5 mol·L-1的 H2SO4溶液中80℃煮沸1 h，最后再在去離子水中80℃煮沸1 h，結(jié)束時用去離子水浸泡備用。

標(biāo)準(zhǔn)曲線的繪制：配制一系列標(biāo)準(zhǔn)濃度的VOSO4的硫酸溶液，應(yīng)用相應(yīng)濃度的 MgSO4溶液作為參比溶液，用紫外分光光度計測出其在波長為762 nm處的吸光度，擬合出標(biāo)準(zhǔn)曲線。

按照圖1所示的裝置進(jìn)行滲透性實(shí)驗(yàn)：電解槽左側(cè)為 VOSO4與 H2SO4的混合溶液，右側(cè)為相應(yīng)離子強(qiáng)度的MgSO4與H2SO4的混合溶液。在不同的實(shí)驗(yàn)工況下，每隔一段時間，用紫外分光光度計檢測右側(cè)電解液（透過側(cè)）中VO2+的濃度，比較VO2+在不同實(shí)驗(yàn)工況下的滲透性，并計算其表觀電遷移率(離子淌度)。

實(shí)驗(yàn)中的參數(shù)見表1。

圖1 釩電池膜滲透實(shí)驗(yàn)裝置Fig.1 Schematic diagram of membrane osmotic apparatus of VFB

表1 實(shí)驗(yàn)中的有關(guān)參數(shù)Table 1 Parameters in experiment

3 結(jié)果與討論

3.1 無電場作用下膜中傳質(zhì)過程

按照圖1的實(shí)驗(yàn)裝置進(jìn)行滲透性實(shí)驗(yàn)，在無外在電場加入時，測試了不同溫度、不同濃度、不同對流狀態(tài)等操作工況下，VO2+的跨膜滲透特性，實(shí)驗(yàn)結(jié)果如圖2所示。從圖中可以看出，對于無電場作用下的滲透傳質(zhì)，各不同操作工況下，透過側(cè)的濃度均隨著時間的增加而增加。在時間為32 h時，電解液溫度為 45℃且有強(qiáng)制對流（250 r·min-1）存在的情況下，透過側(cè)VO2+濃度甚至達(dá)到了0.012 mol·L-1，升高溫度和增加溶液中電解液的湍動，都會加劇電解液的滲透。這是因?yàn)闇囟仍黾?，VO2+運(yùn)動更加劇烈，擴(kuò)散加快。此外，溫度升高，會導(dǎo)致膜體積膨脹，使得Nafion膜致密度下降，從而促進(jìn)了離子的傳質(zhì)過程。當(dāng)電解液被強(qiáng)制湍動時，電解液與膜接觸面的滯流層厚度減小，滯流層內(nèi)的濃度梯度降低，這也加劇了VO2+的跨膜滲透過程。

圖2 無電場作用時滲透側(cè)VO2+的濃度隨時間的變化關(guān)系Fig.2 Relationship between VO2+concentration with time inosmotic side without electric field

通過實(shí)驗(yàn)測定的不同時刻透過側(cè)的 VO2+的濃度，能夠關(guān)聯(lián)出VO2+的擴(kuò)散系數(shù)D，圖3為由透過側(cè)濃度求擴(kuò)散系數(shù)D的關(guān)聯(lián)圖。表2所示為不同工況下濃差擴(kuò)散系數(shù)D的數(shù)值大小。

圖3 無電場作用時由滲透側(cè)VO2+的濃度關(guān)聯(lián)擴(kuò)散系數(shù)圖Fig.3 Diffusion coefficient associated with VO2+concentration in osmotic side without electric field

3.2 電場因素對膜中傳質(zhì)過程的影響

3.2.1 正向電場 全釩液流電池在實(shí)際運(yùn)行過程中，在電場力作用下的電遷移傳質(zhì)過程會對釩離子跨膜滲透產(chǎn)生影響。采用外加恒壓電場的方法，可以很好地等效釩離子真實(shí)的傳質(zhì)過程。

表2 無電場作用時滲透側(cè)VO2+在Nafion117膜中的擴(kuò)散系數(shù)Table 2 Diffusion coefficient of VO2+in Nafion117 membrane without electric field

圖 4表示左側(cè)為 0.5 mol·L-1VOSO4與 2 mol·L-1H2SO4右側(cè)為 0.5 mol·L-1MgSO4溶液的體系，在不同外加電壓下不同時刻透過側(cè)的 VO2+濃度。圖中隨著外加電壓從0 V增加到3 V，透過側(cè)VO2+的濃度隨時間增加的趨勢明顯，這與理論預(yù)期相吻合。

圖 4 正向電場下 0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè)VO2+濃度隨時間的變化關(guān)系Fig.4 Relationship between VO2+concentration with time in osmotic side of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in positive electric field

圖5表示在各不同外加電壓下電場因子的隨時間的變化趨勢。從圖中可以發(fā)現(xiàn)，外加電壓越大，電場因子St越大，電壓為0.5 V時，St為1.2左右，當(dāng)電壓增加到3 V時，St可達(dá)到2.35～2.41，VO2+發(fā)生嚴(yán)重的跨膜滲透傳質(zhì)現(xiàn)象。在不同電壓下，電場因子St均隨著時間的增加而增大，且隨著外電壓的加大，St隨時間增大的速度越來越快。

圖 5 正向電場下 0.5 mol·L-1VO2++2 mol·L-1H2SO4體系電場因子隨時間的變化關(guān)系Fig.5 Relationship between electric factor with time of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in positive electric field

3.2.2 反向電場 當(dāng)外加電場為反向電場時，如圖6所示，電池透過側(cè)的VO2+濃度隨著反向電場數(shù)值的增大而減小。圖 7所示為電場因子在不同反向電壓下隨時間的變化，可以看出：當(dāng)反向電壓從0.5 V增加到1 V時，St從0.78左右降到0.55左右，這與理論預(yù)期相吻合。說明反向電場在抑制離子跨膜滲透防止兩極交叉污染方面有一定的應(yīng)用價值。

圖 6 反向電場下 0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè)VO2+濃度隨時間的變化關(guān)系Fig.6 Relationship between VO2+concentration with time in osmotic side of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system inreversed electric field

3.3 離子濃度對協(xié)同場下膜中傳質(zhì)過程的影響

實(shí)驗(yàn)還考察了 VO2+和 H2SO4濃度對協(xié)同場下離子在膜中傳質(zhì)過程的影響。

圖 7 反向電場下 0.5 mol·L-1VO2++2 mol·L-1H2SO4體系電場因子隨時間的變化關(guān)系Fig.7 Relationship between electric factor with time of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in reversed electric field

圖8 電壓為1 V時不同濃度體系滲透側(cè)VO2+的濃度隨時間的變化關(guān)系Fig.8 Relationship between VO2+concentration with time in osmotic side of different concentration system when voltage is 1 V

如圖8所示，當(dāng)VO2+濃度同樣為0.5 mol·L-1時，在透過側(cè)，加入2 mol·L-1H2SO4的實(shí)驗(yàn)組比不加 H2SO4的實(shí)驗(yàn)組的 VO2+要小，說明 H2SO4的加入有助于抑制電場中VO2+在膜中的傳質(zhì)過程，這是因?yàn)殡娊庖褐?H+能夠比較容易地透過質(zhì)子交換膜，也就占據(jù)了更多的膜內(nèi)通道，客觀上阻礙了釩離子跨膜滲透。

再考察VO2+濃度因素的影響，由于加入側(cè)初始VO2+濃度不同，1.5 mol·L-1VO2++2 mol·L-1H2SO4體系在透過側(cè)的VO2+濃度比0.5 mol·L-1VO2++2 mol·L-1H2SO4體系要大，但是VO2+跨膜傳質(zhì)的程度要輕微得多，這正如圖 9所示，電場對于不含H2SO4的體系跨膜傳質(zhì)過程的強(qiáng)化作用最強(qiáng)，在本實(shí)驗(yàn)組中，1 V電壓的St能夠達(dá)到1.67～1.71，對于較高濃度的VO2+和H2SO4的體系影響最小，1 V電壓的St僅為1.16左右，也即高濃度電解液在真實(shí)運(yùn)行中的電場環(huán)境中，有助于降低正負(fù)極交叉污染現(xiàn)象。

圖9 電壓為1 V時不同濃度體系電場因子隨時間的變化關(guān)系Fig.9 Relationship between electric factor of different concentration system when voltage is 1 V

而且有文獻(xiàn)報道[25]，當(dāng)H2SO4的濃度增加時，也有利于減小膜電阻，這是因?yàn)?，?dāng)溶液中的 H+濃度較大時，由于 Donnan平衡效應(yīng)，離子交換膜中的 H+也相應(yīng)增加，故 H+的遷移速度也會加快，也即膜面電阻減小，膜的電導(dǎo)性增加，這有助于電池性能的提高。

3.4 溫度對協(xié)同場下膜中傳質(zhì)過程的影響

圖10 不同溫度下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè)VO2+濃度隨時間的變化關(guān)系Fig.10 Relationship between VO2+concentration with time in osmotic side of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in different temperature

圖11 不同溫度下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系電場因子隨時間的變化關(guān)系Fig.11 Relationship between electric factor with time of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in different temperature

在釩電池運(yùn)行過程中，溫度的影響至關(guān)重要，溫度升高，離子電遷移率增加，電遷移速度加快，此外溫度升高會使離子運(yùn)動加劇，并且會導(dǎo)致膜的致密性降低，這均會使得釩離子在膜兩側(cè)相互滲透的現(xiàn)象嚴(yán)重化；而且在電池真實(shí)運(yùn)行狀況中，溫度過高，會導(dǎo)致正極的五價釩離子VO2+發(fā)生聚集，致使電解液性能下降，這會嚴(yán)重影響釩電池的正常運(yùn)行。圖10顯示了靜態(tài)下電池左側(cè)為0.5 mol·L-1VOSO4與 2 mol·L-1H2SO4溶液，右側(cè)為0.5 mol·L-1MgSO4與2 mol·L-1H2SO4溶液的體系，在低溫25℃和高溫45℃下，加入1 V電場后透過側(cè)VO2+濃度隨時間的變化情況。溫度對電場影響VO2+在膜中傳質(zhì)的情況也可以從圖11中看出來，在協(xié)同場中，25℃時電場影響因子St為1.44～1.46，當(dāng)溫度升高為45℃時，St為1.54～1.58，并且隨著時間的增加，這種電場加劇透膜傳質(zhì)的程度越發(fā)嚴(yán)重，溫度越高，這一趨勢越明顯。

3.5 強(qiáng)制對流對協(xié)同場下膜中傳質(zhì)過程的影響

水力學(xué)條件是通過改變邊界層厚度從而影響膜與溶液的接觸面的離子的濃度和組成，進(jìn)而影響離子跨膜滲透傳質(zhì)過程。若增加溶液的擾動，例如本實(shí)驗(yàn)中對釩電池溶液進(jìn)行機(jī)械攪拌，則這種強(qiáng)制對流條件下形成的膜與溶液界面的滯流層相對于無攪拌的靜止?fàn)顟B(tài)厚度降低，更加有利于離子的質(zhì)量傳遞過程。圖12為常溫25℃下不同對流狀態(tài)下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè) VO2+的濃度隨時間的變化關(guān)系，可以發(fā)現(xiàn)對流狀態(tài)下透膜傳質(zhì)現(xiàn)象更加嚴(yán)重。這也可以用St隨時間變化圖定量表征，如圖13所示。協(xié)同場中，對流狀態(tài)下的St在實(shí)驗(yàn)時間下由1.465增加到1.495，而靜止?fàn)顟B(tài)下，St的變化范圍是 1.445～1.462。雖然對流加劇了協(xié)同場下釩離子的跨膜傳質(zhì)過程，但是，對流作用降低了膜與溶液界面的滯流邊界層，這有助于H+的傳遞過程，抑制了電池濃差極化，有助于提高電池的能量效率。所以，釩電池在實(shí)際應(yīng)用中，往往刻意設(shè)計外循環(huán)從而造成一種電解液對流環(huán)境以提高電池性能。

圖12 不同對流狀態(tài)下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè)VO2+濃度隨時間的變化關(guān)系Fig.12 Relationship between VO2+concentration with time in osmotic side of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in different convection condition

圖13 不同對流狀態(tài)下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系電場因子隨時間的變化關(guān)系Fig.13 Relationship between electric factor with time of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in different convection condition

3.6 對流和溫度綜合因素對協(xié)同場下膜中傳質(zhì)過程的影響

本實(shí)驗(yàn)還通過研究電池在高溫和攪拌情況下VO2+跨膜滲透過程，考察了對流因素和高溫因素對協(xié)同場中VO2+在Nafion膜中傳質(zhì)過程的綜合影響。

圖14表示了左側(cè)為0.5 mol·L-1VOSO4與2 mol·L-1H2SO4溶液體系，右側(cè)為 0.5 mol·L-1MgSO4與2 mol·L-1H2SO4溶液體系，在加入電場前后滲透側(cè)VO2+濃度隨時間的變化趨勢。從圖中可以看出，在溫度為45℃攪拌轉(zhuǎn)速為250 r·min-1的對流狀態(tài)時，滲透側(cè) VO2+濃度相比于低溫 25℃且電解液體系靜止的情形均要大得多。

圖14 不同溫度和對流狀態(tài)下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系滲透側(cè)VO2+濃度隨時間的變化關(guān)系Fig.14 Relationship between VO2+concentration with time in osmotic side of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system indifferent temperature and convection condition

圖15 不同溫度和對流狀態(tài)下0.5 mol·L-1VO2++2 mol·L-1H2SO4體系電場因子隨時間的變化關(guān)系Fig.15 Relationship between electric factor with time of 0.5 mol·L-1VO2++2 mol·L-1H2SO4system in different temperature and convection condition

表3 各工況下VO2+電遷移率計算值Table 3 Calculated value of electric mobility of VO2+in different condition

圖15直觀反映了在同樣施加大小為1 V的電壓，高溫對流狀態(tài)時電場因子St相比于低溫靜止時要大得多，在45℃并且有對流存在時，St值在實(shí)驗(yàn)時間范圍內(nèi)從1.625增加到1.655，而25℃靜止?fàn)顟B(tài)的St的范圍為1.442～1.474，這表明：在協(xié)同場下高溫對流能明顯加劇釩離子的跨膜傳質(zhì)過程。

3.7 各工況下 VO2+在 Nafion117膜中表觀電遷移率計算值

由實(shí)驗(yàn)所得數(shù)據(jù)，根據(jù)式（7）擬合出VO2+在Nafion117膜中的表觀電遷移率，見表3。

4 結(jié) 論

通過研究協(xié)同場下VO2+在Nafion117膜中的傳遞過程，可知：外加電壓越大，電場因子St越大。在各不同電壓下，電場因子St均隨著時間的增加而增大，且隨著外電壓的加大，St隨時間增大的速度越來越快；反向電場在抑制離子跨膜滲透防止兩極交叉污染方面有一定的應(yīng)用價值；高濃度電解液在真實(shí)運(yùn)行中的電場環(huán)境中，St明顯偏小，有助于降低正負(fù)極交叉污染現(xiàn)象；在協(xié)同場中溫度的提高使得St增加，加劇了釩離子在膜中的傳質(zhì)過程，并且隨著時間的增加，這種電場加劇透膜傳質(zhì)的程度越發(fā)嚴(yán)重，溫度越高，這一趨勢越明顯；電解液在對流狀態(tài)下相比于靜止?fàn)顟B(tài)，St值增加，電場影響傳質(zhì)過程更加明顯。

符號說明

A——膜面積，m2

c——濃度，mol·L-1

D——擴(kuò)散系數(shù)，cm2·s-1

E——電場強(qiáng)度，V·m-1

F——法拉第常數(shù)，C·mol-1

i——電流密度，mA·cm-2

j——傳質(zhì)通量，mol·m-3·s-1

L——正負(fù)極平行電極板間距，m

l——膜厚度，μm

St——電場因子

t——時間，s

U——電壓，V

u——電遷移率，m2·s-1·V-1

V——電解液體積，L

z——電荷數(shù)

Γ——遷移數(shù)

下角標(biāo)

d——擴(kuò)散傳質(zhì)性質(zhì)

e——電遷移性質(zhì)

0 ——初始性質(zhì)

1 ——釩液側(cè)（左側(cè)）性質(zhì)

2 ——透過側(cè)（右側(cè)）性質(zhì)

[1] Ludwig J, Juergen G, Fabjan C,et al. Possible use of vanadium redox flow batteries for energy storage in small grids and stand-alone photovoltaic systems [J].J. Power Sources, 2004, 127(1/2)：98-104.

[2] Baker J. New technology and possible advances in energy storage [J].Energy Policy, 2008, 36(12)：4368-4373.

[3] Jia Zhijun (賈志軍), Song Shiqiang (宋世強(qiáng)), Wang Baoguo (王保國). A critical review on redox flow batteries for electrical energy storage applications [J].Energy Storage Science and Technology(儲能科學(xué)與技術(shù)), 2012, 1(1)：50-57.

[4] Rychcik M, Skyllas-Kazacos M. Characteristics of a new all-vanadium redox flow battery [J].Power Sources, 1988, 22(1)：59-67.

[5] Huang K L, Li X G, Liu Su Q,et al. Research progress of vanadium redox flow battery for energy storage in China [J].Renewable Energy,2008, 33(2)：186-192.

[6] Bae C H, Roberts E P L, Dryfer R A W. Chromium redox couples for application to redox flow batteries [J].Electrochim Acta, 2002, 48(3)：279-287.

[7] Vijayakumar M, Sarah D Burton, Cheng Huang,et al. Nuclear magnetic resonance studies on vanadium(Ⅳ) electrolyte solutions for vanadium redox flow battery [J].Journal of Power Sources, 2010,195：7709-7717.

[8] Li Mengnan (李夢楠), Xie Xiaofeng (謝曉峰) , Yang Chun (楊春),et al. Characteristics of BMIMBF4as anode electrolyte additive for vanadium redox flow battery [J].CIESC Journal(化工學(xué)報), 2011,62(S2)：135-139.

[9] Yang Chun (楊春), Wang Shubo (王樹博), Xie Xiaofeng (謝曉峰) ,et al. Performance influence of hydroxyl radical on graphite felt electrode used in allvanadium redox flow battery [J].CIESC Journal(化工學(xué)報), 2012, 63(S1)：188-193.

[10] Liu Ran (劉然), Liao Xiaoyan (廖孝艷), Yang Chun (楊春),et al.Different treatments of graphite electrode materials for vanadium redox flow battery [J].Chemical Industry and Engineering Progress(化工進(jìn)展), 2011, 30：762-766.

[11] Ding C, Zhang H, Li X,et al. Morphology and electrochemical properties of perfluoro sulfonic acid ionomers for vanadium flow battery applications：effect of side-chain length [J].Chem Sus Chem,2013, 6：1262-1269.

[12] Yang Chun (楊春), Wang Jinhai (王金海), Xie Xiaofeng (謝曉峰),et al. AC impedance of influence of glycerin on all vanadium redox flow battery anodic electrolyte [J].CIESC Journal(化工學(xué)報), 2011,62(S1)：163-167.

[13] Zhao Chengming (趙成明), Xie Xiaofeng (謝曉峰). Calculation of solvation free energies of VO V/VO2+ions in aqueous solution by using cluster-continuum model [J].CIESC Journal(化工學(xué)報), 2012,63(S2)：132-135.

[14] Wang Wenpin (王文嬪), Wang Shubo (王樹博), Xie Xiaofeng (謝曉峰),et al. Nonfluorinated anion exchange membranes based on imidazolium salts for all-vanadium redox flow battery [J].CIESC Journal(化工學(xué)報), 2013, 64(S1)：203-208.

[15] Dong J Y, Zhang H M, Chen J. A simple model for the vanadium redox battery [J].Electrochimical Acta, 2009, 54：6827-6836.

[16] Al-Fetlawi H, Shah A A, Walsh F C. Non-isothermal modelling of the all-vanadium redox flow battery [J].Electrochimical Acta, 2009, 55：78-89.

[17] Zheng Q, Zhang H M, Xing F, Ma X K, Li X F, Ning G L. A three-dimensional model for thermal analysis in a vanadium flow battery [J].Applied Energy, 2014, 113：1675-1685.

[18] Wang Baoguo(王保國). Review of the development of proton exchange membranes in the renewable energy technology [J].Membrane Science and Technology(膜科學(xué)與技術(shù)), 2010, 30(1)：1-8.

[19] Yin Haitao(尹海濤), Wang Baoguo(王保國). Influence of membrane diffusivity on the performance of a single vanadium flow battery [J].Battery Bimonthly(電池), 2006, 36(1)：60-61.

[20] Lü Zhengzhong(呂正中), Hu Songlin(胡嵩麟) ,Luo Xuanli(羅絢麗) ,et al. Influence of proton exchange membrane on the performance of vanadium redox flow battery [J].Chemical Journal of Chinese Universities(高等學(xué)?；瘜W(xué)學(xué)報), 2007, 28(1)：145-148.

[21] Chen Jinqing(陳金慶), Lü Hongling(呂宏凌), Wang Baoguo(王保國) ,Yin Haitao(尹海濤). Study on the process of vanadium ions across membrane based on adsorption-diffusion mechanism [J].Contemporary Chemical Industry(當(dāng)代化工), 2011, 40(9)：893-895.

[22] Sun Hong (孫紅), Luan Lihua(欒麗華), Wu Tiejun(吳鐵軍), Tang Yulan（唐玉蘭）, Wang Xun（王遜）. Mass transfer in proton exchange membrane [J].Journal of Engineering Thermophysics(工程熱物理學(xué)報)，2012, 33(2)：255-258.

[23] Kim Soowhan, Yan Jingling, Schwenzer Birgit, Zhang Jianlu, Li Liyu,Liu Jun, Yang Zhenguo (Gary), Hickner M A. Investigation of sulfonated poly(phenylsulfone) membrane for vanadium redox flow batteries[J].Electrochem. Comm., 2010 , 12(11) ：1650-1653.

[24] Li Wenyue, Liu Jianguo, Yan Chuanwei. Multi-walled carbon nanotubes used as an electrode reaction catalyst for VVO2+for a vanadium redox flow battery [J].Carbon, 2011, 49 (11)：3463-3470.

[25] Tan Ning(譚寧), Huang Kelong（黃可龍）,Liu Suqin(劉素琴). Studies on ion exchange membrane performance in vanadium solution for all vanadium flow battery [J].Chinese Journal of Power(電源技術(shù))，2004, 28(12)：775-778.