電解水式駐留微氣泡減阻技術及其可行性分析

朱效谷，李勇，李文平（清華大學精密儀器與機械學系摩擦學國家重點實驗室，北京100084）

電解水式駐留微氣泡減阻技術及其可行性分析

朱效谷，李勇，李文平
（清華大學精密儀器與機械學系摩擦學國家重點實驗室，北京100084）

在通常的氣幕減阻技術中，由于氣泡的不斷流失容易產生能耗增大、減阻效果下降等問題。該文提出一種利用駐留微氣泡的減阻技術。微氣泡產生于陣列布置于航行體表面的微凹坑中，并在航行過程中穩定駐留于凹坑處，從而有效降低局部摩阻系數；微氣泡通過電解水產生，并可反過來對電解反應進行自適應控制，從而自適應地控制供氣量和氣泡形狀。通過初步數值模擬和實驗，對影響氣泡駐留的因素進行了定性的分析。設計了一種電解反應自適應控制的凹坑結構，用于氣泡形狀和電解反應的自動控制。通過數值計算，對一特定條件下駐留微氣泡的理論減阻能力進行了討論。

駐留微氣泡；減阻；電解水

1 引言

船只在行進時會受到水流的阻力，采取有效的減阻技術可以節能、增速，因而減阻技術的研究具有重要的經濟和軍事意義。水流阻力主要包括：摩擦阻力、粘壓阻力和興波阻力。其中，摩擦阻力通常是主要的，如水面高速船的摩擦阻力占總阻力的40%~50%，而對于水面中低速船及水下航行體，摩擦阻力占其受到的總阻力的比重更是高達70%以上[1]。因此摩阻減阻技術的研究是船舶減阻的關鍵。

自上世紀30年代至今，已發展出柔順壁減阻、表面微觀結構減阻、疏水表面減阻、聚合物添加劑減阻、氣幕減阻等多種減阻技術[2]。其中，氣幕減阻技術因其突出的減阻效果而廣受關注，是目前研究最多的減阻技術之一。

氣幕減阻的原理是，通過一定的機構向航行體與水接觸的界面以一定速度噴射氣體，變單相流為水—氣泡兩相混合流，通過混合流密度、粘度及流動模式的改變，減小物體表面的摩擦阻力[3]。自上世紀70年代起，俄羅斯、美國等國對該減阻技術從實驗上進行了研究，采用噴氣方法，最高得到了80%的減阻率[4]。在國內，海軍工程大學董文才等最早進行氣幕減阻研究，實驗得到了總阻力減少25%以上的結果[5]。

迄今雖然氣幕減阻技術的研究取得了令人矚目的實驗效果，但氣幕減阻面臨兩個技術難題：一是氣泡難以實現對航行體表面的附著，而仿真表明微氣泡的附著更有利于摩擦阻力的減小[6]；再就是大量的注入氣泡瞬間流失，造成補充氣泡所需功耗較大。因此，本研究室提出一種電解水式駐留微氣泡減阻技術[7]，以期通過實現微小氣泡在航行體表面的附著而不流失，達到降低減阻系統的功耗，增大減阻效果的目的。

本文首先介紹電解水式駐留微氣泡減阻技術的原理，而后著重從微氣泡的駐留因素、電解反應自適應控制方法、以及減阻效果三個方面對這一減阻技術進行可行性分析。

2 電解水式駐留微氣泡減阻技術的原理

根據氣幕減阻對微氣泡附著于壁面的要求，最理想的情況是使微氣泡駐留于最初的位置而不隨水流流失，利用氣體隔離液體環境和固體壁面，從而既獲得更高的減阻率，又減少了氣體的消耗。

這種駐留微氣泡減阻技術的關鍵在于兩個方面：一是設計利于氣泡駐留的表面微結構；二是設計合理的供氣方式。航行體行進時，表面存在很強的剪切流場，氣泡易受水流沖刷而流失。采用如凹坑形式的微結構，就可以對氣泡形成束縛，使之能抵抗一定強度的剪切流。采用凹坑結構形成無源駐留氣泡的方法已被應用到微流體[8-9]和滑移現象[10]等方面的研究中。然而采用無源氣泡存在如下幾個問題：（1）這是一種不可控的被動機制，不能保證所有微凹陷內都有穩定的氣泡留存；（2）氣泡形狀不可控，實際上由于水下壓力必然高于大氣壓，微結構內留存的氣泡一般都是內凹形狀的，這并不一定是最佳的氣泡形狀，實際上當內凹過深時基本上也就失去了氣泡減阻的功能；（3）當航行體潛深增大時，氣泡會收縮進凹坑底部，或者在航行體行駛過程中由于偶然原因造成氣泡的振蕩、破裂和脫離，不論哪種情況，都難以自動及時地補充恢復到所需尺寸和形狀，因而會失去氣泡減阻的功能。因此，減阻氣泡的產生和補充應采取主動和可控的方法。現有的氣幕減阻研究和應用中多采用經多孔板注入氣體的方式，但是這種方式在控制上很難實現注入氣體量恰好滿足附著和減阻所需的氣泡形態，事實上，即使在低速的水流環境下，大量的注入氣體都會使氣泡脫離微注氣孔，被水流帶走，難以實現駐留。因此，為實現微氣泡駐留，必須從供氣方法上進行改進，實現供氣量的自適應控制。

通過電解反應，直接利用航行體周圍的水環境產生氣體，既提供了氣源，又可望通過產生的氣泡阻斷反應從而實現供氣量的自適應控制。電解水式駐留微氣泡減阻的原理如圖1所示。其中，圖1（a）和（b）示意了在船體表面應用本減阻技術的一種方式，圖1（c）和（d）示意了氣泡駐留的形式和電解產氣的原理。

圖1（d）所示是一個凹坑內駐留氣泡形成的過程。采用合理的凹坑結構設計和電極布置方式，可以使氣泡在恰好生長至合適的形狀時覆蓋住電極，從而阻斷電解反應，避免氣泡的不斷生長。而當氣泡發生收縮或意外脫離時，水流重新進入微凹坑內部并與電解電極接觸，電解水反應自動重啟，重復圖1（d）所示的過程，直至氣泡恢復形狀后電解水反應再次自動關斷。如此的電解反應的停/啟自動控制便可實現供氣量的自適應控制，從而維持氣泡的形狀。再通過對凹坑結構的優化設計，以提高其對氣泡的束縛作用，就有望實現氣泡在隨載體航行的過程中的穩定駐留。

駐留微氣泡的減阻效果來自氣泡的滑移效果，氣泡的存在隔離了固體表面與水的直接接觸，將固-水界面轉變為固—氣—水界面，相比于前者，后者的局部表面摩擦阻力系數將大大減小。最理想的情況是形成大面積的連續氣膜，但連續氣膜的氣液界面易受流場擾動而變形，氣體容易流失而難以保持。駐留微氣泡減阻技術，采用陣列微氣泡，如圖1（c），單元的特征尺寸在10μm~數百μm，利用氣-液界面的表面張力來保持氣泡的形狀和使氣泡駐留。圖1（c）中以簡單矩形柱凹坑和行列對齊的簡單矩陣型陣列布置對駐留氣泡陣列加以說明。由于每一個凹坑中都獨立地進行著電解反應的自適應控制過程，使得每一個氣泡的形狀都維持在設計的范圍內，因而即便氣泡的數量非常龐大，氣泡的形狀也均勻一致，這也是反應自適應控制的優點的一個方面。

圖1 電解水式駐留微氣泡的原理Fig.1 The principle of trapped micro-bubbles by water electrolysis

美國德州大學Stephani等[11]也提出了類似的通過電解水形成駐留氣泡陣列（trapped-bubble array）來獲取減阻效果的方法。然而，他們的研究沒有涉及氣泡駐留問題；同時，他們也沒有對電解反應的自適應控制加以嘗試，而是選擇人工斷電方式。要想真正提高減阻效果并將其應用于實際航行體上，就必須對氣泡駐留的影響因素和電解反應自適應控制的方法加以研究。本文接下來將就這兩個關鍵問題展開分析，并對駐留微氣泡的減阻效果進行初步的估計。

3 微氣泡駐留的影響因素

氣泡在固體表面受到多種力共同作用，包括重力、浮力、水流剪切力、相界面的表面力、氣泡生長所產生的微觀流場擾動，以及電解反應中的電場力等。在數十微米的尺度上，影響氣泡形態的主要是水流剪切力和表面力兩個方面。其中，在低水流速度下，表面力起主導作用；在高速水流下，水流剪切作用也不可忽視。下面對這兩方面分別進行初步的分析。

水流速較低情況下，水流對氣泡的剪切作用較微弱。其中，極端情況是外界輸入的剪切率等于零，此時氣泡的分布狀況取決于固壁面的浸潤性和結構。氣泡具有向疏水性較強區域遷移的趨勢，因而氣泡易脫離親水性很強的電極材料并遷移至親水性較弱的區域。凹坑的結構具有約束氣泡的作用，而關鍵的參數便是凹坑的深度。如圖2為由凹坑邊緣抽象而出的臺階結構，當臺階頂面接觸角θt大于底面接觸角θb時，一定高度的臺階可將一定體積以下的氣泡限制于底面上。

從不可壓縮流體的N-S（Navier-Stokes）方程出發，采用VOF（volume of fluid）模型作為運動界面追蹤模型，用CSF（continuum surface force）模型考慮表面張力的影響，使用SIMPLE（semi-implicitmethod for pressure-linked equations）類的有限體積法進行了數值求解。使用了通用流體分析軟件FLUENT。數值計算和分析表明[12]，二維情況下，臺階能約束的最大氣泡對應于圖2中的VC，根據其中的接觸角條件可由幾何關系得出此時半徑為：

圖2 凹坑邊緣氣泡形態的二維模型Fig.2 The 2-Dmodel for a bubble at the boundary of a pit

圖3 靜水情況下不同深度的凹坑內氣泡的駐留效果（a）~（c）1μm深凹坑試片中，氣泡長大后偏移出凹坑；（d）15μm深凹坑試片中，氣泡始終穩定地駐留在凹坑內部，直至填滿凹坑Fig.3 Bubbles in pits of different depth in stationary water(a)~(c)A bubble growing in a 1-μm-deep pitmoves out; (d)Bubbles growing in 15-μm-deep pits are all trapped inside

水流速較高情況下，航行體表面附近形成很強的速度梯度，必須考慮水流對氣泡的剪切作用。同樣采用VOF模型，對二維情況下壁面上一個孤立凹坑中的氣泡受剪切流作用的變形進行了數值模擬，如圖4所示。圖4（a）是計算模型，其中，a為凹坑的寬度，水從左向右流過，設定在凹坑上游5倍于其寬度的位置處的流場為均勻梯度場，來流速度梯度為1.5×105s-1。氣泡初始狀態如圖4（a）中所示，為蓋住整個凹坑的半圓形。如此的二維模型在物理上相當于對如圖4（b）所示的氣柱的行為進行研究，其中水在兩側面擋板上的接觸角為90°。雖然這不完全符合圖1所示的微氣泡形式，但二維下的趨勢和規律仍能為進一步的研究和設計帶來啟發。

數值計算的結果如圖4（c）~（g）所示，層流強剪切流作用下與氣泡駐留有關的因素至少包括：固壁面接觸角、凹坑形狀、氣泡尺寸。定性的結論分別如下：

（1）固壁面的潤濕性的影響

比較圖4（c）~（e），隨著接觸角的依次減小，固壁面的親水性依次增強，氣泡受水流沖刷而引起的變形也逐步加劇。θ=15°時，氣泡被完全剪斷；而θ=150°時，氣泡可以維持在一個穩定的形狀。因此疏水性的固壁面利于氣泡的駐留。

（2）凹坑形狀的影響

在圖4（f）的圓弧形凹坑中，氣泡受水流沖刷變形嚴重。與圖4（c）對比可發現，凹坑前緣是造成結果差異的關鍵。對于圓弧形的凹坑，氣泡與凹坑前緣的接觸點無法保持其位置，受水流沖刷后該接觸點后移，接觸點處的氣液界面隨即向后轉動，從而引起氣泡的變形。結合氣泡的變形過程推斷，凹坑前緣為尖角的形狀更利于氣泡的駐留。

（3）氣泡尺寸的影響

圖4（g）中，部分氣泡被剪切流削去并流向下游，與圖4（c）對比后推斷，較小的氣泡在剪切流中更能保持完整。

圖4 對剪切流中二維氣泡的數值模擬Fig.4 Numerical simulation of 2-D bubbles in shear flow

綜合表面力和水流剪切力兩方面的分析來看：固壁面具有較好的疏水性有利于使氣泡在壁面附著時具有抵抗流場剪切力的形狀；從對固壁面的疏水性要求出發，則電極材料具有強親水性的現實就要求采用大深度的凹坑結構；凹坑結構的形狀和尺寸是設計的關鍵，其中尖銳的前緣有利于氣泡的位置保持，而過大的尺寸則不利于氣泡的保持完整。這些定性的結論對于器件的結構、尺寸、物性等參數設計而言具有指導意義，可以為進一步的氣泡駐留實驗提供設計依據。

圖5 反應阻斷與自適應控制的模式（a）模式一：氣泡直接覆蓋住電極；（b）模式二：氣泡切斷電流通路但不與電極接觸Fig.5 Two modes of reaction interruption and adaptive control(a)Mode 1:the bubbles cover the electrodes directly;(b)Mode 2:the bubbles block the electrolyte path but not contact the electrodes

4 電解反應自適應控制的設計

電解水式微氣泡減阻技術還要求對電解水反應能有效控制。如果電解反應不停地進行，則氣泡不斷變大，最終氣泡過分的凸起導致易被水流沖走，從而失去減阻效果，如清華大學陳皓生的實驗[13]中，電解氣泡過大后脫離試片進入液體中，會增大摩阻。美國德州大學Stephani通過人工斷電使氣泡停止生長[11]，但這種簡單斷電的方法存在諸多問題：無法保證斷電時氣泡的形狀利于減阻；同一時間斷電，各微氣泡形狀不一，所得實驗結果的科學性低；氣泡隨時間而慢慢流失后不能及時得到補充，從而失去減阻效果。借助電解反應產生的氣泡來自動地切斷電流通路，從而實現電解反應的自適應控制，則可望大大提高氣泡的一致性，并實現氣泡形狀的高效控制。

反應自適應控制的第一步是電解反應的阻斷。目前鮮有氣泡阻斷反應的研究，一方面因為現有的研究中更多的集中在注入氣泡減阻技術，而對駐留微氣泡減阻技術鮮有關注，另一方面也因為反應的阻斷實現起來存在困難。反應阻斷的困難主要來自材料的浸潤性。實用的電極材料都是強親水性，氣泡無法扁平地覆蓋在電極上，相反，水則易于在電極表面鋪展從而構成電流通路。因而，要實現反應的完全阻斷，需要采用合理的結構設計并選用適當的材料。

采用凹坑結構，除了束縛氣泡使之駐留外，也同時具有阻斷反應的作用。理論上看，假定凹坑側壁是光滑的圓柱面，電極布置在凹坑底部，且一個凹坑內只有一極，而氣液界面與固面的夾角也滿足理想的接觸角條件，則當氣泡體積足夠大后必然可以封住凹坑的截面，便自然實現反應的阻斷。圖5展示了基于凹坑結構的兩種氣泡阻斷電流通路的形式，一種是氣泡直接覆蓋住電極，從而隔斷電極與水的接觸；一種是氣泡將連通陰陽兩極的水從中間某處切斷，而電極依舊保持與水的接觸。

實際的固壁面上的接觸角存在滯后性，這造成實際的接觸角與氣-液界面的運動方向有關，因而反應阻斷效果和氣泡形狀會較理想條件下復雜。親水材料的接觸角滯后性更明顯，同時有研究表明親水表面具有特定的紋理時會大大提高液體的鋪展能力[14-17]，從而導致水路的連通，因而應試圖提高所選擇材料的疏水性。圖5（a）中的模式一，凹坑側壁為電極材料，這種反應阻斷模式對電極材料的疏水性能提出了要求，需要開發疏水性好的導電材料；對于圖5（b）中的模式二，凹坑側壁是絕緣的，可以選擇常見的疏水材料，如聚四氟乙烯、硅橡膠等。

除反應的阻斷外，也要求電解反應在氣泡收縮或意外脫離時自動重啟。反應阻斷和重啟的自動控制是駐留微氣泡減阻技術的關鍵問題，其目的是將氣泡的形狀控制在一個最佳的范圍內，以達到最佳的減阻效果。圖5的兩種方式分別對應了不同的自適應控制模式。對第一種，當氣泡縮小至暴露出電極時，反應重啟，氣泡再次生長，至電極重新被覆蓋。這種模式下，由于電極材料與其它區域的絕緣材料的浸潤性不同，需考慮接觸角的變化來確定反應停/啟的不同階段氣泡形狀的變化范圍。對第二種，當氣泡縮小至陰陽兩極重新由水連通時，反應重啟，氣泡再次生長，至水路重新被氣泡切斷。這種模式下，氣泡不與電極材料接觸，則假使各處絕緣材料的物性一致，氣泡的形狀變化范圍就取決于絕緣材料的物性和凹坑的形狀。

采用上述的第二種模式，圖6給出了一種電解反應自適應控制的結構設計。其結構如圖6（a）所示，包括基底、結構層Ⅰ、結構層Ⅲ共三層，結構層Ⅰ和結構層Ⅲ共同構成凹坑結構。凹坑的結構主體上是圓柱形，僅在結構層Ⅰ的側壁一處開有沿軸向貫通的V形槽。基底層中布置有電極，且每一個凹坑內只有一極。這種設計的自適應控制原理如圖6（b）和（c）所示。凹坑內部的電極與其它凹坑中的異性電極之間形成電解電流，產生氣體。當氣泡生長到一定體積后就會封住圓柱形凹坑的截面，而尖銳的V形槽結構中的液體可以繼續保持在槽道中[16]，從而繼續維持電解電流，如圖6（b）。氣泡依靠通過V形槽的電流繼續生長，直至上部的氣-液界面上升至結構層Ⅲ處時，V形槽中的導電水路也被氣體與外界隔開，于是電解反應被阻斷，如圖6（c）。反應阻斷后，氣體停止產生，氣泡由于溶解或逆反應等原因而開始緩慢收縮。上部的氣-液界面逐漸下降至V形槽上端時，導電水路重新連通，此時又回到圖6（b），反應自動重啟。由上述的自適應控制周期過程來看，這種結構設計下的氣泡形狀會在如圖6（c）中所示的δ范圍內變化，δ的具體值由結構尺寸和材料的接觸角決定。

圖6 一種用于電解反應自適應控制的結構設計Fig.6 A design for adaptive control of electrolysis reaction

5 駐留微氣泡的減阻效果

以固—氣—液界面代替固-液界面從而減小摩阻的思想由來已久。如氣墊船、超空泡魚雷[18]等氣膜減阻的技術已被投入實際的應用，并已得到了顯著的減阻效果。相比而言，駐留微氣泡減阻技術則研究甚少，其減阻效果尚處研究階段。

Stephani[11]的實驗，由于測量條件和反應控制技術等方面的限制，尚未能得到明顯的減阻效果。雖然目前還沒有直接的實驗驗證，但邊界滑移效果的研究則間接為駐留微氣泡減阻技術提供了依據。近年來在微流體的研究中越來越多地觀察到滑移現象，這吸引了大量研究者的興趣。不少學者對分布有氣體的表面的滑移特性進行了實驗和理論研究[19-22]，結果顯示出駐留有氣泡的表面具有很強的滑移效果。滑移現象的產生突破了傳統流體理論的限制，為提升減阻效果提供了更大的理論空間。

下面通過數值計算來初步判斷駐留氣泡的減阻能力，而所述內容均以矩陣型布置的圓柱形豎直凹坑陣列這種最簡單的形式展開討論。由于氣體的密度和動力粘度都遠小于液體，因而將氣泡覆蓋的區域近似以完全滑移的邊界條件處理。已有研究表明，將氣泡近似作完全滑移的邊界所得的流場計算結果與實驗結果的規律相符[20]，而由流場結果便可計算出阻力。圖7是駐留氣泡陣列中的一小塊區域的簡化模型，其中，相對于水下航行體表面復雜的邊界層，本模型在計算時采用了三點簡化：

（1）水下航行體表面存在層流邊界層、湍流邊界層以及過渡區，然而即便在湍流邊界層內，貼近壁面的粘性底層內總是處于層流狀態。鑒于微氣泡駐留于航行體表面，且尺寸微小，故而采用層流模型。

（2）水下航行體所受摩阻是由表面剪切流引起的，本模型采用平板Couette流來形成剪切流動，其中的近似包括兩個方面：航行體表面的邊界層是不斷往后發展的，壁面上各點的剪切率都不同，而在小范圍內可近似認為均勻一致；另一方面，邊界層內速度分布雖然隨高度呈非線性，但在很薄的粘性底層內接近線性分布。

（3）微尺度下，氣泡受表面張力的作用尤為明顯，因而忽略其受剪切流的粘性力作用，而近似為剛性的球冠[23]。

圖7 駐留氣泡減阻效果的數值計算模型Fig.7 The numerical analysismodel for drag reduction by trapped bubbles

假定所有的氣泡都駐留于規則分布的凹坑結構內，并且形狀均勻一致，如圖7所示，則穩定時，流場沿流向呈周期性分布。采用周期邊界條件對單個周期范圍內的流場進行計算，即可得出整個流場的分布以及剪切流的阻力。

取圖7中Couette流的上板高度H為100μm，其移動速度u為1m/s，這接近常見水下航行體運行時的層流底層的情況。在此給定的流場條件下，影響阻力的因素應當包括：氣泡覆蓋率[11]，氣泡凸起角度θ[10]和周期λ[21]。當氣泡凸起角度θ為10°，周期λ為500μm，圓柱形凹坑直徑450μm即覆蓋率為63.6%時，計算得出單個周期內的阻力為1.89μN。而對應于同樣條件下的光滑平板Couette流，相同范圍內受到的摩阻理論值易由牛頓內摩擦定律得出為2.51μN。也就是說，以完全滑移邊界條件模擬上述參數條件的駐留氣泡，得出了近25%的減阻效果。值得注意的是，高的減阻率意味著氣泡的存在明顯地改變了邊界層的粘性底層內的流場，使其與Couette流的差異增大，因此使用平板Couette流模型進行計算的誤差也會增大。

上述初步的阻力計算，以及前人有關駐留氣泡所產生的滑移現象的實驗和理論成果，都直接或間接為駐留微氣泡減阻技術提供了支持。但是，確切的減阻效果與多方面的因素有關，除了上文已經提及的，如氣泡覆蓋率、氣泡凸起角度以及周期等因素外，氣泡的形狀和分布方式更是復雜多樣，需要大量的研究來揭示其中的規律以指導設計。另外，要進行精確測定減阻效果的實驗研究，需要以電解反應的自適應控制作為前提手段，來實現氣泡形狀的有效控制。

6 結論

提出一種電解水式駐留微氣泡減阻技術：在航行體外表面布置內有電極的陣列微凹坑，通過電解水反應生成氣泡且穩定駐留在微凹坑中，變固—液界面為固—氣—液界面，有效降低局部摩阻系數；減阻氣泡形成后阻斷電極與水的接觸，電解反應自動中止，避免氣泡繼續生長導致流失。

這種電解水式駐留微氣泡減阻技術與現有氣幕減阻技術相比，氣泡基本附著于航行體表面而不隨水流流失，從而既可以減少氣源系統能耗的浪費，又易于控制氣體分布。同時通過反應的自適應控制，有望獲得優化的氣泡形狀，從而進一步地提高減阻性能。

對電解水式駐留微氣泡減阻技術的可行性進行了初步的分析。從微氣泡的駐留因素、反應自適應控制的方法、駐留微氣泡的減阻效果三個關鍵方面分析了駐留微氣泡減阻技術的關鍵問題。分析得出，固壁面材料的浸潤性、微凹坑結構的形狀和尺寸均是影響氣泡駐留的重要因素。設計了一種側壁具有溝槽的凹坑結構，有望實現電解反應的自適應控制。將邊界層的粘性底層近似為couette流，并將氣-液界面處理成完全滑移的剛性球冠，對駐留氣泡的減阻效果進行了數值計算。得出當駐留氣泡陣列的周期為500μm，氣泡凸起角度為10°，表面的氣體覆蓋率為63.6%時，對于相距100μm且相對速度為1m/s的平板couette流，可獲得接近25%的減阻率。

[1]王晉軍.溝槽面湍流減阻研究綜述[J].北京航空航天大學學報,1998,24(1):31-34. Wang JJ.Reviews and prospects in turbulent drag reduction over riblets surface[J].Journal of Beijing University of Aeronautics and Astronautics,1998,24(1):31-34.(In Chinese)

[2]柯貴喜,潘光,黃橋高,胡海豹,劉占一.水下減阻技術研究綜述[J].力學進展,2009,39(5):546-554. Ke G X,Pan G,Huang Q G,Hu H B,Liu Z Y.Reviews of underwater drag reduction technology[J].Advances in Mechanics,2009,39(5):546-554.(In Chinese)

[3]董文才,郭日修.氣幕減阻研究進展[J].船舶力學,1998,2(5):73-78. Dong W C,Guo R X.Research progress of air injection drag reduction[J].Journal of Ship Mechanics,1998,2(5):73-78. (In Chinese)

[4]Madavan N K,Deutsch S,Merkle C L.Reduction of turbulent skin friction bymicrobubbles[J].Physics of Fluids,1984, 27(2):356-363.

[5]董文才,郭日修,陳小玲,呂巖松.滑行艇氣層減阻試驗[J].中國造船,2002,43(4):13-18. Dong W C,Guo R X,Chen X L,Lu Y S.Experimental study on resistance reduction of planing craft by air injection[J]. Ship Building of China,2002,43(4):13-18.(In Chinese)

[6]吳乘勝,何術龍.微氣泡流的數值模擬及減阻機理分析[J].船舶力學,2005,9(5):30-37. Wu C S,He SL.Numerical simulation ofmicrobubble flow and analysis of themechanism of drag reduction[J].Journal of Ship Mechanics,2005,9(5):30-37.(In Chinese)

[7]陳旭鵬,李勇,朱效谷,黃偉峰.柔性MEMS減阻蒙皮及其制造方法[P].中國:ZL200910079713.0,2009-7-22. Chen X P,Li Y,Zhu X G,Huang W F.Flexible MEMSanti-drag skin and its fabrication[P].China:ZL200910079713.0, 2009-7-22.

[8]Ahmed D,Mao X,Juluri B K,Huang T J.A fastmicrofluidic mixer based on acoustically driven sidewall-trapped microbubbles[J].Microfluid Nanofluid,2009,7(5):727-731.

[9]Kang E,Lee D H,Kim C-B,Yoo S J,Lee S-H.A hemisphericalmicrofluidic channel for the trapping and passive dissipation ofmicrobubbles[J].Journal ofMicromechanics and Microengineering,2010,20(4):045009.

[10]Steinberger A,Cottin-Bizonne C,Kleinmann P,Charlaix E.High friction on a bubblemattress[J].Nature Materials,2007, 6(9):665-668.

[11]Stephani K A,Goldstein D B.An examination of trapped bubbles for viscous drag reduction on submerged surfaces[J]. Journal of Fluids Engineering-Transactions of the ASME,2010,132(4):041303.

[12]Zhu X G,Huang W F,Li Y.Role of step-shape structures and critical failure bubble volume in microbubble constraint [J].Chinese Science Bulletin,2012,57(22):2941-2946.

[13]Chen H,Li J,Chen D.Study of drag forces on a designed surface in bubbly water lubrication using electrolysis[J].Journal of Fluids Engineering-Transactions of the ASME,2006,128(6):1383-1389.

[14]Seemann R,Brinkmann M,Kramer E J,Lange F F,Lipowsky R.Wettingmorphologies atmicrostructured surfaces[J]. Proceedings of the National Academy of Sciences of the United States of America,2005,102(6):1848-1852.

[15]Brinkmann M,Blossey R.Blobs,channels and cigars-Morphologies of liquids at a step[J].European Physical Journal E, 2004,14(1):79-89.

[16]Khare K,Herminghaus S,Baret J-C,Law B M,Brinkmann M,Seemann R.Switching liquid morphologies on linear grooves[J].Langmuir,2007,23(26):12997-13006.

[17]Yuan Q,Zhao Y-P.Topology-dominated dynamic wetting of the precursor chain in a hydrophilic interior corner[J].Proceedings of the Royal Society A,2012,468(2138):310-322.

[18]曹偉,魏英杰,王聰,鄒振祝,黃文虎.超空泡技術現狀、問題與應用[J].力學進展,2006,36(4):571-579. CaoW,Wei Y J,Wang C,Zou Z Z,HuangW H.Current status,problems and applications of supercavitation technology [J].Advances in Mechanics,2006,36(4):571-579.(In Chinese)

[19]Ou J,Perot B,Rothstein JP.Laminar drag reduction inmicrochannels using ultrahydrophobic surfaces[J].Physics of Fluids,2004,16(12):4635-4643.

[20]Ou J,Rothstein JP.Direct velocitymeasurements of the flow past drag-reducing ultrahydrophobic surfaces[J].Physics of Fluids,2005,17(10):103606.

[21]Joseph P,Cottin-Bizonne C,Benoit J-M,Ybert C,Journet C,Tabeling P,Bocquet L.Slippage of water past superhydrophobic carbon nanotube forests inmicrochannels[J].Physical Review Letters,2006,97(15):156104.

[22]Sbragaglia M,Prosperetti A.Effective velocity boundary condition at amixed slip surface[J].Journal of Fluid Mechanics, 2007,578:435-451.

[23]Hyvaluoma J,Harting J.Slip flow over structured surfaceswith entrapped microbubbles[J].Physical Review Letters,2008, 100(24):246001.

Feasibility analysis of drag reduction using trapped m icro-bubbles by water electrolysis

ZHU Xiao-gu,LIYong,LIWen-ping
(State Key Laboratory of Tribology,Departmentof Precision Instrumentsand Mechanology, Tsinghua University,Beijing 100084,China)

In air injection drag reduction,bubbleswill be inevitably lostwith water flow,and thismay cause energy waste and goes against drag reduction.This paper provides a new method of drag reduction by trapped micro-bubbles.Micro-bubbles are generated in arrayed pits on ships by electrolysis reaction,and each bubble stably stays there,so as to reduce the local friction.Further,this electrolysis reaction in every pit can be adaptively controlled,so each bubble can automatically keep its volume and even its shape.To show the feasibility of thismethod,influencing factors on bubbles'stability are qualitatively discussed by numerical analysis and experiments.A special structurewith a V-groove is given for adaptive control of electrolysis reaction.And the drag reduction ability of trapped micro-bubbles under a certain condition is predicted through numerical calculation.

trapped micro-bubbles;drag reduction;water electrolysis

U661.3

A

10.3969/j.issn.1007-7294.2014.10.002

1007-7294（2014）10-1165-10

2014-04-25

高等學校博士學科點專項科研基金資助項目（20110002110077）

朱效谷（1985-），男，清華大學博士研究生，E-mail:sxg03@mails.tsinghua.edu.cn；

李勇（1962-），男，研究員，博士生導師，E-mail:liyong@mail.tsinghua.edu.cn。