重疊網(wǎng)格技術(shù)下振蕩水翼的水動力特性

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

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

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

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

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

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

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

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

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

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

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

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

(1)

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

2)動量方程。

(2)

3)湍流方程。

SSTk-ω模型由Menter[15]提出，它綜合了標準k-ω模型近壁區(qū)計算和標準k-ε模型遠場計算的優(yōu)點，適用于翼型計算：

(3)

(4)

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

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

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

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

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為水翼的弦長c與展長l的乘積，l在二維計算中取單位長度。

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

(10)

(11)

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

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

(12)

(13)

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

(14)

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

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

(15)

1.3 模型驗證

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

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

本文使用了重疊網(wǎng)格技術(shù)，網(wǎng)格密度調(diào)整時應(yīng)當按照一定比例，同時協(xié)調(diào)組分網(wǎng)格和背景網(wǎng)格的密度，改變翼型表面節(jié)點個數(shù)的同時，背景網(wǎng)格的密度以及其各個加密區(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)性驗證

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

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

1.3.2 模型有效性驗證

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

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

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

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

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

首先研究單翼在水流中振蕩運動時的水動力特性，本文針對NACA0016型水翼進行分析，來流速度U∞=2 m/s，弦長c=0.25 m。

2.1.1 折算頻率

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

圖2 折算頻率對獲能效率的影響

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

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

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

(16)

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

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

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

圖5 振蕩水翼運動過程中的渦脫落

2.1.2 擺角幅值

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

H1=2(h-L1sinθ0)

(17)

H2=2(h+L2sinθ0)

(18)

那么振蕩水翼的實際掃掠高度H為

H=max(H1,H2)

(19)

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

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

從圖6可以看出：隨著擺角幅值的增大，總體上水翼獲能效率會得到提高；因為小幅度地增加擺角幅值會提高水翼的有效攻角，會提高獲能效率。但是過大的擺角幅值會使振蕩水翼的有效攻角接近失速角，降低振蕩翼的水動力性能，并且擴大了水翼對來流的掃掠高度H，由式(15)知水翼的獲能效率η也會減小。

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

2.1.3 升沉振幅

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

2.1.4 轉(zhuǎn)軸位置

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

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

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

2.2 翼地效應(yīng)

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

2.2.1 雙翼間距

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

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

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

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

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

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

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

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

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

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

K=CL/CD

(20)

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

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

3 結(jié) 論

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

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

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

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