2、網(wǎng)格無(wú)關(guān)性計(jì)算

圖3所示為氣泡直徑為789μm時(shí)不同網(wǎng)格數(shù)量下的氣液界面Marangoni對(duì)流速度變化曲線。在不同網(wǎng)格數(shù)量下，Marangoni對(duì)流速度隨氣泡弧長(zhǎng)增大均呈現(xiàn)先急劇增大、后緩慢減小的趨勢(shì)。當(dāng)氣泡弧長(zhǎng)大于0.1 mm時(shí)，Marangoni對(duì)流速度幾乎相同，因此，主要對(duì)比氣泡接觸點(diǎn)附近的氣泡弧長(zhǎng)區(qū)域的Marangoni對(duì)流速度分布曲線。當(dāng)網(wǎng)格數(shù)量增大時(shí)，Marangoni對(duì)流速度逐漸增大；當(dāng)網(wǎng)格數(shù)量大于313萬(wàn)時(shí)，Marangoni對(duì)流速度曲線幾乎相同。為提高計(jì)算效率和保證計(jì)算精度，后續(xù)相關(guān)數(shù)值模擬均采用313萬(wàn)的網(wǎng)格數(shù)量進(jìn)行網(wǎng)格劃分。

圖3網(wǎng)格無(wú)關(guān)性計(jì)算

3、溫度場(chǎng)和Marangoni對(duì)流結(jié)構(gòu)

微電極表面氣泡生長(zhǎng)過(guò)程中，電極附近的電解液產(chǎn)生焦耳熱，特別是微液層內(nèi)產(chǎn)生的熱量比較大，導(dǎo)致氣泡周圍形成電解液溫度分布特性。如前所述，氫氣泡周圍的物理場(chǎng)參數(shù)呈現(xiàn)以圓形微電極中心的法線為中心軸的對(duì)稱結(jié)構(gòu)，因此，圖4僅顯示單個(gè)氫氣泡垂直中心右邊截面的電解液溫度場(chǎng)分布與Marangoni對(duì)流分布結(jié)果。由圖4(a)可知：氣泡底部微液層內(nèi)溫度較大，氣液界面周圍電解液溫度較小。這是由于氣泡電絕緣和電流線集中扭曲效應(yīng)所導(dǎo)致的，微液層內(nèi)電流密度較大和電導(dǎo)率較小，導(dǎo)致電解液焦耳熱主要產(chǎn)生于微液層內(nèi)。這些熱量隨著對(duì)流擴(kuò)散逐漸影響到更多的氣泡周圍電解液區(qū)域，從氣泡接觸點(diǎn)到氣泡頂點(diǎn)的氣泡弧長(zhǎng)方向，電解液溫度逐漸降低，溫度最大值在氣泡接觸點(diǎn)位置，導(dǎo)致氣液界面附近區(qū)域形成明顯的溫度梯度分布，距離氣液界面較遠(yuǎn)區(qū)域內(nèi)的溫度近似保持為環(huán)境溫度25℃。

圖4單個(gè)氫氣泡周圍溫度場(chǎng)與Marangoni對(duì)流分布

氣液界面溫度梯度分布導(dǎo)致形成表面張力梯度分布，其中氣泡接觸點(diǎn)附近表面張力較小，沿著氣泡弧長(zhǎng)方向表面張力逐漸增大，從而促使電解液從低表面張力向高表面張力方向運(yùn)動(dòng)，形成明顯的不穩(wěn)定的Marangoni對(duì)流旋渦結(jié)構(gòu)，如圖4(b)所示。氣液界面附近區(qū)域的電解液沿著氣泡弧長(zhǎng)方向流動(dòng)，Marangoni對(duì)流速度逐漸減小，最大速度出現(xiàn)在氣泡接觸點(diǎn)附近的右上側(cè)位置。由于微液層內(nèi)不考慮溫度梯度引起的Marangoni對(duì)流結(jié)果，因此，微液層內(nèi)的電解液流動(dòng)速度幾乎為零，近似認(rèn)為微液層邊緣區(qū)域?yàn)楣腆w邊界條件，導(dǎo)致Marangoni對(duì)流到達(dá)接觸點(diǎn)位置受到一定的阻礙，使最大速度出現(xiàn)在氣泡接觸點(diǎn)附近的右上側(cè)位置。Marangoni對(duì)流旋渦的外邊緣區(qū)域速度較大，旋渦內(nèi)部區(qū)域的速度較小，具體Marangoni對(duì)流旋渦結(jié)構(gòu)及演變規(guī)律與整個(gè)氣泡生長(zhǎng)周期不同時(shí)刻的氣泡直徑有關(guān)。

4、Marangoni對(duì)流演變規(guī)律

為了深入探明氫氣泡單個(gè)生長(zhǎng)周期內(nèi)氣液界面溫度梯度誘導(dǎo)產(chǎn)生的Marangoni對(duì)流結(jié)構(gòu)的演變規(guī)律，本文對(duì)表2所示的不同氣泡直徑和接觸直徑下的Marangoni對(duì)流效應(yīng)進(jìn)行了數(shù)值模擬。圖5所示為微電極表面電極中心點(diǎn)到氣泡接觸點(diǎn)區(qū)域內(nèi)的電解液溫度差變化規(guī)律，每條溫度差曲線最右側(cè)端點(diǎn)的橫坐標(biāo)為氣泡接觸半徑Rc。由圖5可見(jiàn)：隨著氣泡不斷生長(zhǎng)，氣泡直徑持續(xù)增大，氣泡接觸直徑先增大、再近似保持不變、最后再減小且近似保持不變。從電極中心點(diǎn)到氣泡接觸點(diǎn)位置區(qū)域內(nèi)，不同氣泡生長(zhǎng)時(shí)刻的電解液溫度差均呈現(xiàn)先緩慢減小、然后快速減小的趨勢(shì)。微電極邊緣附近區(qū)域和接觸點(diǎn)區(qū)域內(nèi)的溫度差減小趨勢(shì)非常明顯，這與圖4(a)一致。這是由于靠近氣泡底部的氣液界面的Marangoni對(duì)流速度較大，電解液焦耳熱能夠快速傳輸?shù)礁蠓秶碾娊庖罕倔w區(qū)域內(nèi)。當(dāng)氣泡直徑逐漸增大時(shí)，電解液溫度差逐漸增大，但是溫度差增大的趨勢(shì)逐漸減小。氣泡生長(zhǎng)引起氣泡直徑增大，進(jìn)而導(dǎo)致氣泡電絕緣和電流線集中扭曲效應(yīng)也更加明顯，促使微液層內(nèi)電流密度進(jìn)一步增大，產(chǎn)生更多的焦耳熱量。當(dāng)氣泡直徑從189μm增大到467μm時(shí)，氣泡直徑和接觸直徑增大明顯，上述解釋可進(jìn)一步闡述溫度差增大幅度較大的原因。而當(dāng)氣泡直徑大于578μm時(shí)，氣泡生長(zhǎng)過(guò)程中的接觸直徑變化趨勢(shì)非常小，且氣泡直徑增大的幅度也較小，因此，總體上電解液溫度差的增大幅度較小。

圖5不同氣泡直徑(不同氣泡生長(zhǎng)時(shí)刻)下微液層內(nèi)溫度差變化

圖6所示為不同氣泡直徑下氣液界面電解液溫度差和Marangoni對(duì)流速度的變化規(guī)律。由圖6(a)可知：在任意氣泡直徑下，氣液界面電解液溫度差均沿著氣泡弧長(zhǎng)方向逐漸減小。氣泡直徑越大，氣液界面電解液溫差也較大，這主要與電極中心到氣泡接觸點(diǎn)區(qū)域內(nèi)溫度分布及演變規(guī)律有密切關(guān)聯(lián)。在相同氣泡弧長(zhǎng)下，大直徑氣泡的氣液界面溫度差比較大，意味著氣泡生長(zhǎng)后期的大直徑氣泡周圍的溫度梯度較大。不同氣泡直徑下，明顯溫度梯度均主要體現(xiàn)在氣泡接觸點(diǎn)附近大約0.3 mm距離內(nèi)的氣泡弧長(zhǎng)界面上，距離氣泡接觸點(diǎn)稍遠(yuǎn)的氣液界面上幾乎不存在溫度梯度。因此，當(dāng)表面張力隨溫度變化的系數(shù)保持不變時(shí)，大直徑氣泡界面的表面張力梯度也較大，從而導(dǎo)致氣液界面Marangoni對(duì)流速度近似隨著氣泡直徑增大而增大，但這種變化規(guī)律對(duì)于不同直徑的氣泡來(lái)說(shuō)，體現(xiàn)在不同的氣泡弧長(zhǎng)區(qū)域范圍內(nèi)。

由圖6(b)可知：當(dāng)氣泡直徑較小時(shí)，如氣泡直徑小于733μm，上述變化規(guī)律主要體現(xiàn)在除氣泡接觸點(diǎn)附近極小區(qū)域之外的其他絕大部分氣泡弧長(zhǎng)界面內(nèi)。當(dāng)氣泡直徑僅為189μm，氣液界面Marangoni對(duì)流速度反而較大，這可能與氣泡生長(zhǎng)初期極小氣泡的接觸直徑及氣泡周圍的Marangoni對(duì)流旋渦結(jié)構(gòu)有關(guān)。而當(dāng)氣泡直徑較大時(shí)，如氣泡直徑大于844μm，上述變化規(guī)律主要體現(xiàn)在接觸點(diǎn)附近區(qū)域的部分有限氣泡弧長(zhǎng)界面內(nèi)。不同氣泡生長(zhǎng)時(shí)刻下的氣液界面Marangoni對(duì)流速度大于10 mm/s，數(shù)量級(jí)明顯大于實(shí)驗(yàn)獲得的氣泡生長(zhǎng)速度0.01～1.00 mm/s，這再次說(shuō)明本文建立的固定氣泡直徑模型是可靠合理的。

圖6不同氣泡直徑(不同氣泡生長(zhǎng)時(shí)刻)下溫差和Marangoni對(duì)流速度

圖7所示為不同氣泡直徑下氣液界面Marangoni對(duì)流旋渦結(jié)構(gòu)的變化規(guī)律。相對(duì)于圖4(b)的局部顯示結(jié)果，圖7給出了不同直徑氣泡周圍的更加直觀的Marangoni對(duì)流旋渦形態(tài)及演變過(guò)程。當(dāng)氣泡直徑較小時(shí)，旋渦結(jié)構(gòu)中心大概在氣泡的右上端邊緣區(qū)域。隨著氣泡生長(zhǎng)變大，氣泡周圍的Marangoni對(duì)流旋渦形態(tài)發(fā)生演變，特別是旋渦結(jié)構(gòu)中心逐漸上移。對(duì)于任意氣泡直徑下，相對(duì)于氣泡弧長(zhǎng)的上端頂點(diǎn)所在的高度而言，旋渦結(jié)構(gòu)中心位置的高度較小。隨著氣泡直徑增大，氣液界面最大Marangoni對(duì)流速度所在的氣泡弧長(zhǎng)圓周角度φ逐漸減小，但都大于對(duì)應(yīng)的氣泡接觸角θ，如圖8所示。最大電解液溫差出現(xiàn)在氣泡接觸角θ所在的氣泡弧長(zhǎng)位置上，進(jìn)一步說(shuō)明沿著氣泡弧長(zhǎng)方向，最大Marangoni對(duì)流速度的位置稍微偏離于氣泡接觸點(diǎn)位置。

圖7不同氣泡直徑(不同氣泡生長(zhǎng)時(shí)刻)下Marangoni對(duì)流結(jié)構(gòu)演變

圖8不同氣泡生長(zhǎng)時(shí)刻下最大電解液溫差和Marangoni對(duì)流速度的角度位置

5、結(jié)論

2)氣液界面溫度梯度分布導(dǎo)致形成表面張力梯度分布，促使電解液從氣泡底部沿著氣泡弧長(zhǎng)方向流動(dòng)，形成不穩(wěn)定的Marangoni對(duì)流旋渦結(jié)構(gòu)。沿著氣泡弧長(zhǎng)方向的Marangoni對(duì)流速度逐漸減小，最大速度出現(xiàn)在氣泡接觸點(diǎn)附近的右上側(cè)界面位置。

3)從電極中心點(diǎn)到氣泡接觸點(diǎn)位置區(qū)域內(nèi)，電解液溫度差先緩慢減小、后快速減小。當(dāng)氣泡直徑逐漸增大時(shí)，微液層內(nèi)電解液溫度差逐漸增大，但是溫度差增大的趨勢(shì)逐漸減小。

4)在相同氣泡弧長(zhǎng)下，氣泡直徑越大，氣液界面溫度差越大。氣液界面Marangoni對(duì)流速度近似隨著氣泡直徑增大而增大，但這種變化規(guī)律對(duì)于不同直徑的氣泡來(lái)說(shuō)，體現(xiàn)在不同的氣泡弧長(zhǎng)區(qū)域范圍內(nèi)。具體Marangoni對(duì)流旋渦結(jié)構(gòu)及演變規(guī)律與整個(gè)氣泡生長(zhǎng)周期不同時(shí)刻的氣泡尺寸有關(guān)。