何廣華，莫惟杰，王 威，劉朝綱

(1.哈爾濱工業(yè)大學(xué)(威海) 海洋工程學(xué)院，山東 威海 264209；2.哈爾濱工業(yè)大學(xué) 機(jī)電工程學(xué)院，哈爾濱 150001；3.山東船舶技術(shù)研究院，山東 威海 264209)

天然化石能源依然是當(dāng)今世界能源消耗的主要來源；這些不可再生能源的大量消耗，在為全球經(jīng)濟(jì)高速發(fā)展帶來動(dòng)力的同時(shí)，也給人類的生活環(huán)境質(zhì)量帶來了極大的挑戰(zhàn)[1]。為此，世界各國(guó)都在積極地開發(fā)和利用可再生的環(huán)保能源，其中潮流能依靠太陽和月球的引力產(chǎn)生動(dòng)能，這種可再生的清潔能源具有周期性特點(diǎn)，能量密度是太陽能的30倍、風(fēng)能的4倍[2]，且資源分布豐富。因海水的密度約為空氣的1 000倍，發(fā)電機(jī)輸出相同功率時(shí)，水輪機(jī)的半徑只需風(fēng)機(jī)的1/2[3]，所以在近海處潮流密集地區(qū)，尤其流速高于2 m/s時(shí)，非常適宜布置潮流能發(fā)電裝置。

潮流能發(fā)電的裝置按工作原理的不同主要分為水平軸水輪機(jī)、垂直軸水輪機(jī)、振蕩水翼式發(fā)電裝置。其中水平軸水輪機(jī)的開發(fā)和利用最為完備和成熟；垂直軸水輪機(jī)具有對(duì)于來流方向不敏感的特點(diǎn)，但是能量利用率偏低。振蕩水翼式發(fā)電裝置的能量捕獲效率較高，淺海領(lǐng)域的應(yīng)用優(yōu)勢(shì)明顯[4]，并且加工工藝簡(jiǎn)單，應(yīng)用前景較好。

1981年，Mckinney等[5]進(jìn)行了利用振蕩翼捕獲來流能量的實(shí)驗(yàn)，并就實(shí)驗(yàn)值和理論解進(jìn)行了對(duì)比，并對(duì)理論公式進(jìn)行了修正。Zhu[6]設(shè)計(jì)了被動(dòng)式單翼的獲能機(jī)構(gòu)，依靠彈簧提供的回復(fù)力實(shí)現(xiàn)水翼在線性剪切流作用下的升沉、俯仰運(yùn)動(dòng)，發(fā)現(xiàn)在小剪切率下水翼仍能保持周期性振蕩。王俊芳[7]研究了二維水翼在不同攻角、雷諾數(shù)、擺動(dòng)頻率和幅值條件下的水動(dòng)力性能。

專家學(xué)者針對(duì)雙水翼獲能效率的研究取得了一定的成果。王勇等[8]對(duì)NACA0015型雙水翼的獲能效率進(jìn)行了研究，指出雙水翼耦合振蕩捕獲潮流能裝置在一定條件下獲能效率可達(dá)40%. 劉海賓等[9]研究了二維雙水翼的俯仰振幅存在90°的相位關(guān)系時(shí)的獲能效率，并指出了NACA0018翼型的獲能效率相對(duì)較好。黃鳳躍[10]基于Theodoersen理論分析了大攻角NACA0018振蕩薄翼非定常繞流的水動(dòng)力特性。馬鵬磊等[11]研究了垂直軸振蕩水翼的水動(dòng)力特性，指出水翼在不同的俯仰振幅下均存在“最優(yōu)頻率”使水翼的能量提取效果達(dá)到最佳。謝玉東等[12]研究了振蕩水翼在非均勻流條件下的能量提取性能，指出水翼的攻角在30°附近時(shí)，流速在垂直方向的非均勻分布會(huì)提高水翼的時(shí)均功率系數(shù)。

近些年，研究改進(jìn)翼型水翼的水動(dòng)力的特性也逐漸受到關(guān)注。馬鵬磊等[13]研究了改進(jìn)翼型做俯仰運(yùn)動(dòng)時(shí)的水動(dòng)力特性，通過改變水翼振蕩的運(yùn)動(dòng)方式使其在潮流中獲取更多的能量。郭春雨等[14]分別通過試驗(yàn)和數(shù)值模擬比較了普通前緣水翼和波浪形前緣水翼的水動(dòng)力性能，發(fā)現(xiàn)波浪形前緣水翼結(jié)構(gòu)可有效減弱表面渦流分離對(duì)水翼的擾動(dòng)，提高水翼的水動(dòng)力性能。

已發(fā)表的數(shù)值研究結(jié)果多是針對(duì)前后串聯(lián)式雙水翼、翼間距較大的并行雙水翼，而關(guān)于地翼效應(yīng)的雙翼研究并不多見。為了實(shí)現(xiàn)小間距雙水翼的模擬，本文采用重疊網(wǎng)格技術(shù)對(duì)振蕩水翼的水動(dòng)力性能進(jìn)行研究，并分析翼地效應(yīng)在振蕩雙翼接近時(shí)的影響。

1 數(shù)值模擬方法及驗(yàn)證

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

本文所用數(shù)學(xué)模型的相關(guān)控制方程如下：

1)連續(xù)性方程。

(1)

式中：ρ為流體密度；t為時(shí)間；ui為流體介質(zhì)在笛卡爾坐標(biāo)軸i方向上的速度分量；xi為i軸方向的坐標(biāo)，i=1,2,3。

2)動(dòng)量方程。

(2)

3)湍流方程。

SSTk-ω模型由Menter[15]提出，它綜合了標(biāo)準(zhǔn)k-ω模型近壁區(qū)計(jì)算和標(biāo)準(zhǔn)k-ε模型遠(yuǎn)場(chǎng)計(jì)算的優(yōu)點(diǎn)，適用于翼型計(jì)算：

(3)

(4)

1.2 水翼的流域邊界條件及獲能效率

本文采用重疊網(wǎng)格技術(shù)處理水翼的運(yùn)動(dòng)問題，計(jì)算域的邊界條件如圖1所示，水翼的邊界為壁面條件。水翼周圍的網(wǎng)格屬于組分網(wǎng)格，它重疊在流域背景網(wǎng)格內(nèi)部[16]，組分網(wǎng)格在背景網(wǎng)格平面內(nèi)運(yùn)動(dòng)。兩部分網(wǎng)格之間通過挖洞和嵌入實(shí)現(xiàn)耦合[17]，即數(shù)值計(jì)算時(shí)先將背景網(wǎng)格中與組分網(wǎng)格對(duì)應(yīng)的區(qū)域刪除，然后通過相鄰網(wǎng)格單元的插值實(shí)現(xiàn)組分網(wǎng)格和背景網(wǎng)格的數(shù)據(jù)傳遞。

圖1 計(jì)算流域邊界條件及重疊網(wǎng)格結(jié)構(gòu)

通過自定義函數(shù)(user define function,UDF)的方式實(shí)現(xiàn)組分網(wǎng)格的運(yùn)動(dòng)，如圖1所示，水翼的運(yùn)動(dòng)方程為:

h=h0cos(2πft)

(5)

θ=θ0sin(2πft)

(6)

式中：h為水翼的垂向位移，h0為升沉振幅；f為水翼的振蕩頻率；θ為水翼繞轉(zhuǎn)軸的擺角(俯仰角)，θ0為擺角幅值。

水翼的阻力系數(shù)CD、升力系數(shù)CL、力矩系數(shù)CM分別為:

(7)

(8)

(9)

式中:FD、FL、M分別為水翼的阻力、升力、力矩；U∞為來流速度；S為水翼的弦長(zhǎng)c與展長(zhǎng)l的乘積，l在二維計(jì)算中取單位長(zhǎng)度。

水翼獲能的瞬時(shí)功率P為升力FL與俯仰力矩M做功之和，即

(10)

(11)

式中T為水翼的振蕩周期。

水翼的瞬時(shí)功率系數(shù)CP為

(12)

(13)

振蕩水翼的總功率P0定義為

(14)

式中H為水翼的掃掠高度，它與展長(zhǎng)l的乘積為水翼的掃掠面積。

振蕩水翼在潮流中的獲能效率η為

(15)

1.3 模型驗(yàn)證

運(yùn)用上述數(shù)學(xué)模型對(duì)振蕩水翼做數(shù)值模擬分析，按照Kinsey等[18]給出的運(yùn)動(dòng)參數(shù)，選取NACA0015翼型，取弦長(zhǎng)c=0.25 m，折算頻率[19]f*=fc/U∞=0.14，θ0=60°，h0=0.25 m，水翼轉(zhuǎn)軸位置L1取為c/3弦長(zhǎng)位置處，來流速度U∞=2 m/s，分析流域邊界尺度、網(wǎng)格密度對(duì)計(jì)算收斂性的影響，使得水動(dòng)力計(jì)算結(jié)果對(duì)以上兩者不再敏感。

1.3.1 網(wǎng)格無關(guān)性驗(yàn)證

本文使用了重疊網(wǎng)格技術(shù)，網(wǎng)格密度調(diào)整時(shí)應(yīng)當(dāng)按照一定比例，同時(shí)協(xié)調(diào)組分網(wǎng)格和背景網(wǎng)格的密度，改變翼型表面節(jié)點(diǎn)個(gè)數(shù)的同時(shí)，背景網(wǎng)格的密度以及其各個(gè)加密區(qū)的網(wǎng)格密度也都要一起做相應(yīng)調(diào)整。本文根據(jù)相同的比例設(shè)置不同的網(wǎng)格數(shù)量(見表1)，網(wǎng)格A和B的獲能效率分別和較密的網(wǎng)格C比較，網(wǎng)格B的誤差為1%以內(nèi)，選擇網(wǎng)格B的劃分密度作為本文的網(wǎng)格密度。

表1 網(wǎng)格密度無關(guān)性驗(yàn)證

采用重疊網(wǎng)格技術(shù)中，水翼所在的組分流域網(wǎng)格需在背景流域網(wǎng)格內(nèi)周期運(yùn)動(dòng)，在過小的計(jì)算域中，湍流難以得到充分發(fā)展，容易出現(xiàn)邊壁反射、出口回流等影響計(jì)算精度的問題，因此對(duì)流域的尺寸做無關(guān)性驗(yàn)證。選取網(wǎng)格B中的密度標(biāo)準(zhǔn)，對(duì)不同流域尺度驗(yàn)證結(jié)果見表2?？梢园l(fā)現(xiàn)當(dāng)流域尺度大于70c×70c后，水翼的效率趨于穩(wěn)定，本文選取的流域尺度為70c×70c.

表2 計(jì)算域收斂性驗(yàn)證數(shù)據(jù)表

1.3.2 模型有效性驗(yàn)證

選取流域邊界尺寸為：70c×70c，網(wǎng)格數(shù)目為8.2萬，θ0=60°和θ0=75°，水翼的水動(dòng)力計(jì)算結(jié)果與文獻(xiàn)中的翼型研究結(jié)果對(duì)比見表3、4。從表3、4可以看出，計(jì)算結(jié)果和文獻(xiàn)[18]的結(jié)果相差5%以內(nèi)，充分驗(yàn)證了本文數(shù)值計(jì)算方法的有效性。

表3 振蕩水翼的計(jì)算結(jié)果驗(yàn)證(θ0=60°)

表4 振蕩水翼的計(jì)算結(jié)果驗(yàn)證(θ0=75°)

2 振蕩水翼數(shù)值模擬結(jié)果與分析

2.1 水翼運(yùn)動(dòng)及結(jié)構(gòu)參數(shù)的影響

首先研究單翼在水流中振蕩運(yùn)動(dòng)時(shí)的水動(dòng)力特性，本文針對(duì)NACA0016型水翼進(jìn)行分析，來流速度U∞=2 m/s，弦長(zhǎng)c=0.25 m。

2.1.1 折算頻率

圖2給出了振蕩水翼在不同折算頻率條件下的獲能效率變化曲線，水翼的擺角幅值θ0=60°，升沉振幅h0=1c。

圖2 折算頻率對(duì)獲能效率的影響

從圖2可以看出：水翼獲能效率η隨著折算頻率f*的增加呈現(xiàn)先增大后減小的趨勢(shì)。f*在0.085～0.095之間的獲能效率較高，最高可達(dá)29%左右；f*<0.085時(shí)，因水翼折算頻率的增加，水翼的垂向運(yùn)動(dòng)速度和擺動(dòng)角速度逐漸增大，水翼的獲能效率逐步增加；f*>0.095以后，水翼的獲能效率下降。為分析獲能效率下降的原因，f*分別取0.095和0.143時(shí)，水翼的3自由度水動(dòng)力時(shí)歷曲線如圖3所示。

圖3 振蕩水翼兩個(gè)周期內(nèi)的受力時(shí)歷曲線

從圖3可以看出，隨著擺動(dòng)頻率的增大，水翼的阻力幅值略有提高，但文中水翼并無水平方向的運(yùn)動(dòng)，所以阻力對(duì)獲能效率的影響較小；升力幅值對(duì)應(yīng)的時(shí)間點(diǎn)向前移動(dòng)，但幅值變化不大，對(duì)獲能效率的影響不大；而力矩系數(shù)因水翼振蕩頻率的減小發(fā)生了較大變化，頻率較小的時(shí)候，力矩系數(shù)的峰值處會(huì)出現(xiàn)反方向的“塌陷”，這種力矩變化造成了水翼獲能效率的差異。因力矩的做功和角速度有關(guān)，文中定義量綱一的角速度Θ為

(16)

力矩和角速度為同一方向時(shí)，水翼是捕獲能量的，取力矩和角速度的時(shí)歷曲線如圖4所示。圖4中以量綱一的時(shí)間t/T=1.5左右的范圍為例，折算頻率為0.095條件下，力矩系數(shù)的峰值處的“塌陷”使得力矩和水翼的轉(zhuǎn)動(dòng)同向，力矩對(duì)水翼做正功，水翼處于獲能狀態(tài)；而折算頻率為0.143條件下，力矩系數(shù)和水翼的轉(zhuǎn)動(dòng)反向，力矩對(duì)水翼做負(fù)功，水翼的獲能效率下降。

圖4 振蕩水翼的力矩系數(shù)和量綱一的角速度的時(shí)歷曲線

分析小頻率條件下水翼力矩系數(shù)峰值“塌陷”的原因，流場(chǎng)漩渦結(jié)構(gòu)的變化是影響水動(dòng)力的直接因素[13]，取t/T=1.35～1.50的范圍為例，水翼的漩渦結(jié)構(gòu)如圖5所示。該時(shí)刻水翼逐漸擺動(dòng)到最下方的水平位置，因迎/背流面的壓差作用，水翼前緣表面因流動(dòng)分離而產(chǎn)生漩渦，漩渦附近的壓力較低，漩渦沿水翼表面向下游運(yùn)動(dòng)時(shí)，會(huì)影響水翼所受的力矩方向，從而出現(xiàn)了圖3(c)中的力矩系數(shù)曲線的“塌陷”變化。

圖5 振蕩水翼運(yùn)動(dòng)過程中的渦脫落

2.1.2 擺角幅值

由式(15)知，水翼的掃掠高度H對(duì)水翼的獲能效率η有重要影響。從圖1可知水翼轉(zhuǎn)軸距離前緣長(zhǎng)度為L(zhǎng)1，距離后緣長(zhǎng)度為L(zhǎng)2，那么前緣和后緣的掃掠高度H1和H2分別為:

H1=2(h-L1sinθ0)

(17)

H2=2(h+L2sinθ0)

(18)

那么振蕩水翼的實(shí)際掃掠高度H為

H=max(H1,H2)

(19)

可見振蕩水翼的擺角幅值θ0對(duì)獲能效率η存在影響。設(shè)定水翼的轉(zhuǎn)軸位置L1=c/3，研究不同的擺角幅值對(duì)水翼獲能效率的影響，如圖6所示。

圖6 擺角幅值對(duì)水翼獲能效率的影響

從圖6可以看出：隨著擺角幅值的增大，總體上水翼獲能效率會(huì)得到提高；因?yàn)樾》鹊卦黾訑[角幅值會(huì)提高水翼的有效攻角，會(huì)提高獲能效率。但是過大的擺角幅值會(huì)使振蕩水翼的有效攻角接近失速角，降低振蕩翼的水動(dòng)力性能，并且擴(kuò)大了水翼對(duì)來流的掃掠高度H，由式(15)知水翼的獲能效率η也會(huì)減小。

從圖6可以發(fā)現(xiàn)：每個(gè)擺角幅值曲線中最大獲能效率所對(duì)應(yīng)的折算頻率f*也會(huì)隨擺角幅值的增大而增大；擺角幅值在80°附近時(shí)，振蕩水翼最大獲能效率η達(dá)到42.89%，說明振蕩水翼的擺角幅值處于75°～80°時(shí)，水翼的獲能性能較好，這與文獻(xiàn)[9]的結(jié)論一致。雖然較大擺角幅值的水翼在高頻條件下的獲能效率要比較小擺角幅值的水翼在低頻條件下的獲能效率高，但是較高的擺動(dòng)頻率對(duì)潮流的流動(dòng)環(huán)境及振蕩水翼發(fā)電裝置的穩(wěn)定性也提出了較高的要求，所以追求提升振蕩水翼獲能效率的時(shí)候，要兼顧最優(yōu)獲能頻率和擺角幅值的綜合影響。

2.1.3 升沉振幅

圖7 升沉振幅對(duì)水翼獲能效率的影響

2.1.4 轉(zhuǎn)軸位置

為綜合研究振蕩水翼的水動(dòng)力性能，進(jìn)一步分析了水翼的轉(zhuǎn)軸位置對(duì)獲能效率的影響，轉(zhuǎn)軸距離水翼前緣的距離為L(zhǎng)1，不同轉(zhuǎn)軸位置水翼的獲能效率變化如圖8所示。

從圖8可以看出：隨著轉(zhuǎn)軸位置L1的逐漸增大，水翼的獲能效率η曲線呈現(xiàn)先增大后減小的趨勢(shì)；當(dāng)水翼的轉(zhuǎn)軸位于距離前緣的1/3弦長(zhǎng)附近時(shí)，水翼的獲能效率最佳。后文的研究中取1/3弦長(zhǎng)的位置作為振蕩水翼的轉(zhuǎn)軸。

圖8 轉(zhuǎn)軸位置對(duì)水翼獲能效率的影響

2.2 翼地效應(yīng)

翼地效應(yīng)是指運(yùn)動(dòng)物體貼近地面飛行時(shí)，地面對(duì)物體產(chǎn)生的空氣動(dòng)力干擾，影響到空氣對(duì)飛行器的繞流特性，使飛行器上下壓力差增大、升力增加，獲得比空中飛行更高升阻比的流體力學(xué)效應(yīng)。對(duì)于平行對(duì)稱布置的雙翼，隨著上下兩個(gè)水翼的不斷接近，雙翼之間會(huì)形成類似翼地效應(yīng)的高壓，使雙翼分別獲得更大的升力，提高雙翼的能量捕獲性能。

2.2.1 雙翼間距

根據(jù)上文的研究結(jié)果，選取水翼的擺角幅值θ0=80°，折算頻率f*=0.16，轉(zhuǎn)軸位置L1=c/3，雙翼之間的最小距離設(shè)為1～6倍弦長(zhǎng)時(shí)的能量捕獲效率見表5。

表5 振蕩水翼間距對(duì)能量捕獲效率的影響

從表5可以看出對(duì)稱雙翼的上下間距擴(kuò)大過程中，水翼的能量捕獲效率逐漸減小并趨于平穩(wěn)，但是變化范圍并不大；說明水翼間距較遠(yuǎn)時(shí)雙翼彼此的影響會(huì)逐漸減小，并與各翼單獨(dú)運(yùn)動(dòng)的情況趨于相同，所以研究對(duì)稱雙翼間的翼地效應(yīng)對(duì)獲能效率影響主要側(cè)重于雙翼間距小于1倍弦長(zhǎng)的情況。

2.2.2 翼地效應(yīng)分析

對(duì)稱雙翼之間的運(yùn)動(dòng)規(guī)律是對(duì)稱的，它們的相位差為π. 為避免雙翼發(fā)生碰撞，選取雙翼的擺角幅值θ0=60°；根據(jù)圖6的結(jié)果，選擇折算頻率f*=0.08，振蕩水翼的工況見表6。

表6 翼地效應(yīng)的工況

在水翼擺動(dòng)數(shù)個(gè)周期后回到平衡位置時(shí)，提取表6中流場(chǎng)的壓力分布見表7。比較發(fā)現(xiàn)，相較于單振蕩水翼的工況1，雙振蕩水翼的工況2系統(tǒng)中，兩水翼的距離接近時(shí)，翼面之間形成了高壓區(qū)。這是因?yàn)殡S著雙翼距離的不斷靠近，兩個(gè)水翼相互充當(dāng)“地面”，依靠翼地效應(yīng)在雙水翼之間形成了高壓區(qū)，并使水翼兩側(cè)壓差增大，提高了水翼的升力，使水翼的能量捕獲效率η從工況1中單翼情況的24.86%達(dá)到工況2中雙翼情況的25.27%，而工況3的能量捕獲效率達(dá)到25.41%。

表7 翼地效應(yīng)下壓力云圖

工況3中流場(chǎng)的壓力分布與工況2相似，工況4主要起到對(duì)比分析的作用，給出了單個(gè)振蕩水翼與固定平板組成的能量捕獲系統(tǒng)的流場(chǎng)壓力分布情況。從表7工況4中可以看出當(dāng)水翼與平板的距離較近時(shí)，兩者之間的確形成了高壓區(qū)域，水翼會(huì)獲得更大的升力，此時(shí)的獲能效率為27.51%。而水翼與平板或者雙水翼的距離較大時(shí)，流場(chǎng)的壓力分布狀況與單水翼相似。通過以上分析可知，雙翼彼此之間在系統(tǒng)中互相充當(dāng)了“地面”，形成了翼地效應(yīng)。

進(jìn)一步分析水翼的升阻比性能，升阻比K定義為水翼升力系數(shù)和阻力系數(shù)的比值如下:

K=CL/CD

(20)

在水翼逐漸遠(yuǎn)離運(yùn)動(dòng)軌跡平衡位置(或水翼從平板附近遠(yuǎn)離)的時(shí)候，分析水翼的升阻比變化特性。因形成翼地效應(yīng)的時(shí)間較短，所以只分析前T/64的時(shí)間段即可，各工況水翼的升阻比K的變化如圖9所示?？梢园l(fā)現(xiàn)工況1中單水翼的升阻比相對(duì)較小，其最大值僅為38.4；工況2中雙水翼的升阻比最大值為51.8；工況3中雙水翼的升阻比最大值為87.4；工況4單水翼和平板系統(tǒng)中，升阻比高達(dá)101.9；說明雙翼的距離越小，翼地效應(yīng)對(duì)振蕩水翼升阻比的影響越明顯。

圖9 翼地效應(yīng)時(shí)各工況水翼的升阻比曲線圖

3 結(jié) 論

1)振蕩水翼的能量捕獲效率η隨著折算頻率f*的升高而先增加后減少，折算頻率處于0.085～0.095區(qū)間時(shí)，水翼的能量捕獲效率較高，而渦脫落對(duì)獲能效率存在重要影響。

2)振蕩水翼的獲能效率η總體上隨著擺角幅值θ0的增加而增加，但受有效攻角αe和掃掠高度H的影響，獲能效率的提升要兼顧最優(yōu)獲能頻率和擺角幅值的綜合影響，當(dāng)折算頻率處于獲能最優(yōu)區(qū)間時(shí)，擺角幅值在75°～80°的范圍內(nèi)，振蕩水翼的獲能性能較好。

3)振蕩水翼的升沉振幅h0在1倍弦長(zhǎng)附近，轉(zhuǎn)軸距離前緣1/3弦長(zhǎng)位置時(shí)，水翼的能量捕獲效率較好。

4)雙振蕩水翼互相靠近時(shí)，兩翼之間相互作用形成了翼地效應(yīng)的現(xiàn)象，水翼之間的流場(chǎng)形成高壓區(qū)域，水翼的升阻比大幅度提高，提升了水翼的能量捕獲性能。