王爽，盧鵬*，祖永恒，張麗敏，王慶凱，李志軍

(1.大連理工大學(xué) 海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室，遼寧 大連 116024)

1 引言

海冰受到大氣和海洋的動(dòng)力作用產(chǎn)生漂移，其漂移速度與拖曳系數(shù)大小密切相關(guān)[1-3]。拖曳系數(shù)表征氣-冰-海系統(tǒng)的動(dòng)量交換效率。隨著海冰拖曳系數(shù)參數(shù)化趨于完善，其量值從以前的經(jīng)驗(yàn)常數(shù)發(fā)展到隨著海冰形態(tài)變化的變量[4]。海冰的形態(tài)參數(shù)包括密集度、浮冰尺寸、冰脊間距和冰脊厚度等。由于冰脊厚度一般大于海冰的平均厚度，從而影響流速的分布，產(chǎn)生遮掩影響[5]，而冰脊間距又會(huì)影響遮掩作用的強(qiáng)度。因此在海冰拖曳系數(shù)參數(shù)化中，遮掩函數(shù)把冰脊間距和冰脊厚度作為自變量來定量描述相鄰冰脊間下游冰脊受到的拖曳力隨著這兩個(gè)變量發(fā)生的衰減。遮掩函數(shù)直接決定下游冰脊所受拖曳力的大小，對(duì)海冰漂移產(chǎn)生不容忽視的影響。

目前，海冰拖曳系數(shù)參數(shù)化中應(yīng)用的遮掩函數(shù)主要有指數(shù)函數(shù)和冪函數(shù)兩類。Hanssen-Bauer 和Gjessing[6]首次將指數(shù)衰減函數(shù)應(yīng)用到冰-氣界面拖曳系數(shù)參數(shù)化公式中，計(jì)算風(fēng)場中浮冰間的衰減影響。Lüpkes 等[7]將該指數(shù)函數(shù)應(yīng)用到大氣中性分層條件下的冰-氣界面拖曳系數(shù)參數(shù)化方案。Tsamados 等[8]也選用該指數(shù)函數(shù)用于同樣形式的冰-氣和冰-水界面的拖曳系數(shù)參數(shù)化公式中，用來描述冰脊在風(fēng)場和流場中的遮掩效應(yīng)（簡稱TS2014）。此外，Tennekes 和Lumley[9]提出遮掩強(qiáng)度會(huì)隨著流體中結(jié)構(gòu)物間的距離而變化，從而得到遮掩效應(yīng)的冪函數(shù)形式。Steele 等[10]首次將該冪函數(shù)應(yīng)用在冰-水界面形拖曳力參數(shù)化公式中。Lu 等[4]將此衰減函數(shù)應(yīng)用于冰-水界面并提出了系統(tǒng)的冰-水拖曳系數(shù)參數(shù)化方案（簡稱LU2011）。綜上可見，目前使用的冰脊拖曳力的遮掩函數(shù)都是從其他相關(guān)研究中借鑒而來的，并非針對(duì)冰脊的實(shí)際情況進(jìn)行研究的結(jié)果；而且，指數(shù)函數(shù)是從陸地風(fēng)場防護(hù)林研究中得出并直接應(yīng)用于氣-冰-水拖曳系數(shù)參數(shù)化公式中。這種移植和借鑒是否可行，是否適合冰脊間流場遮蔽效應(yīng)的實(shí)際情況仍值得討論。

作為解決上述問題的第一步，本研究關(guān)注冰-水界面上冰脊間的遮蔽效應(yīng)，通過在實(shí)驗(yàn)室物理模型試驗(yàn)中模擬冰下流場，測量前、后冰脊拖曳力的差異，探究冰脊間的尾流遮掩效應(yīng)，初步給出以冰脊的形態(tài)參數(shù)為變量的遮掩函數(shù)公式，用于描述多冰脊間流場衰減情況，為海冰動(dòng)力學(xué)模式的進(jìn)一步完善提供參考。

2 多冰脊拖曳力實(shí)驗(yàn)研究

本次實(shí)驗(yàn)室物理模擬實(shí)驗(yàn)在大連理工大學(xué)海岸和近海工程國家重點(diǎn)實(shí)驗(yàn)室的PIV 水槽完成，該水槽長為22 m，寬為0.45 m，深為0.6 m，有效實(shí)驗(yàn)長度為11.5 m。本次實(shí)驗(yàn)水深為0.45 m，如圖1a 所示。

圖1 實(shí)驗(yàn)示意圖（a）和物理模型實(shí)驗(yàn)裝置圖（b）Fig.1 The sketch of the experimental study (a) and physical model test device diagram (b)

為使本次模擬冰下流場的物理模型實(shí)驗(yàn)接近北極海冰的真實(shí)流況，對(duì)主要參數(shù)采用動(dòng)力學(xué)相似處理。首先問題的依賴參數(shù)有冰脊拖曳力F、冰脊間距L、冰脊入水深度H、水深D、流速V、冰脊底角a、流體密度ρ和流體黏性系數(shù)ν[11]。根據(jù)量綱分析原理選取H、ρ和V為獨(dú)立量綱量，可以得到表1 所示的無量綱量相似依據(jù)。

表1 模型實(shí)驗(yàn)相似依據(jù)Table 1 Similarity basis of model experiment

物理模型實(shí)驗(yàn)的設(shè)計(jì)保證無量綱參數(shù)的范圍和真實(shí)海域的無量綱參數(shù)保持一致[10,12]。北極海冰的漂移速度一般不超過0.3 m/s，冰脊入水深度為2～10 m，冰脊底角范圍為20°～50°[5,12]，雷諾數(shù)Re≥20 000。選用模型和原型的長度比尺λ=1/100，根據(jù)模型相似比計(jì)算得出模型的冰脊入水深度的范圍是0.02～0.10 m。模型形狀為底角為45°的等腰直角三角形，采用有機(jī)玻璃材料制作。與Waters 和Bruno[13]、Pite 等[14]和Zu 等[15]研究類似，該材料的表面粗糙度與真實(shí)海冰不同，但其影響僅限于表面摩阻力，對(duì)于本研究關(guān)注的形拖曳力的影響不超過0.3%，可以忽略[4]。選取的實(shí)驗(yàn)流速V分別為0.1 m/s、0.2 m/s、0.25 m/s、0.3 m/s、0.35 m/s。為使來流更加接近均勻流，本次實(shí)驗(yàn)加設(shè)穩(wěn)流器，將其放置在入水口處。冰脊的入水深度H分別設(shè)定為0.02 m、0.04 m、0.08 m 和0.12 m，則物理模型實(shí)驗(yàn)的雷諾數(shù)Re為2 000～30 000。當(dāng)Re≥10 000時(shí)，我們認(rèn)為物理模型實(shí)驗(yàn)環(huán)境與實(shí)際海洋環(huán)境相似[11]。冰脊的間距L是本次實(shí)驗(yàn)研究的主要變量，由于PIV 水槽長度的限制，4 種入水深度選取了不同的間距，分別為2.5、5、7.5、10、20、30、50、70、100，共9 個(gè)間距。實(shí)驗(yàn)工況共145組。前后冰脊模型在豎直方向上連接量程為10 N 的拉力傳感器，再固定到可移動(dòng)架上，靜置在水槽上方。為了排除造流機(jī)、移動(dòng)架等振動(dòng)對(duì)實(shí)驗(yàn)帶來的干擾影響，本次實(shí)驗(yàn)加設(shè)空載，將直徑為6 mm 的圓桿連接拉力傳感器固定在可移動(dòng)架上，如圖1b 所示。為了減少不同位置流速的差異對(duì)實(shí)驗(yàn)準(zhǔn)確性的影響，將流速儀和前冰脊模型的位置固定不變，通過移動(dòng)后冰脊模型的位置來改變冰脊間距。由此可測量出冰脊模型在水流方向上受到的形拖曳力。

3 結(jié)果

3.1 前冰脊拖曳力

本實(shí)驗(yàn)研究相鄰冰脊間的遮掩影響而產(chǎn)生的冰-水拖曳力衰減，物理模型實(shí)驗(yàn)結(jié)果顯示，前冰脊未受遮掩效應(yīng)的影響，其拖曳力值可以作為后續(xù)遮掩影響分析的基準(zhǔn)值。圖2 給出了前冰脊所受拖曳力F1隨著無量綱參數(shù)間距L/H的變化。

圖2 前冰脊拖曳力F1 隨無量綱間距L/H 的變化Fig.2 Variation of drag force of front ice ridge F1 with dimensionless distance L/H

從圖2 可以看出，盡管前后冰脊間距不斷變化，前冰脊拖曳力F1都表現(xiàn)出了良好的穩(wěn)定性，其數(shù)值大小幾乎與L/H無關(guān)，因此可以認(rèn)為前冰脊沒有受到后冰脊的影響。圖2 中只呈現(xiàn)了H=0.04 m 和H=0.08 m兩個(gè)入水深度F1的變化情況，其他入水深度組次結(jié)果的變化趨勢(shì)類似。當(dāng)然，流速較大時(shí)拖曳力值的浮動(dòng)略大，其主要的原因是流速較大時(shí)，造流電機(jī)運(yùn)行不穩(wěn)定。綜上，對(duì)于多冰脊實(shí)驗(yàn)研究，前冰脊拖曳力的變化規(guī)律與單個(gè)冰脊實(shí)驗(yàn)研究中的規(guī)律保持一致。目前已經(jīng)有較完善的研究成果[4,16]，故本文將不再贅述。

3.2 后冰脊拖曳力

3.2.1 流速的影響

圖3 給出了在不同入水深度以及冰脊間距的情況下，后冰脊拖曳力F2隨流速平方的變化規(guī)律。

從圖3 中可以看出，后冰脊拖曳力F2隨著水流速度V的增大而增大，且F2與V2存在良好線性關(guān)系，與經(jīng)典拖曳力公式吻合[17]。值得注意的是，在L/H≤7.5的小間距情況下，F(xiàn)2出現(xiàn)負(fù)值，即后冰脊產(chǎn)生反向拖曳力，其數(shù)值大小隨著流速增大而反向增大，也保持著與V2的線性關(guān)系，與正值拖曳力形成對(duì)稱現(xiàn)象。后冰脊拖曳力F2的大小和反轉(zhuǎn)與冰下流場的形態(tài)密切相關(guān)。前冰脊尾流場的影響區(qū)域隨著冰脊入水深度的增加而增大，而與流速無關(guān)[1]。在尾流后形成的漩渦場會(huì)產(chǎn)生相對(duì)低壓區(qū)，這是產(chǎn)生后冰脊拖曳力的原因。隨著冰脊入水深度增大，漩渦范圍逐漸變大，漩渦中心也逐漸遠(yuǎn)離冰脊[3]，這使得在冰脊間距較近時(shí)，后冰脊落在前冰脊漩渦范圍內(nèi)從而產(chǎn)生反向拖曳力值；并隨著冰脊間距的增加，拖曳力開始由負(fù)值到0 值再到正值。這一過程中，多冰脊遮掩的疊加對(duì)于浮冰所受到的總冰脊形拖曳力會(huì)產(chǎn)生較大改變，對(duì)于冰脊局地拖曳系數(shù)的確定也將產(chǎn)生較大的影響。

另外，圖3 中冰脊拖曳力隨流速平方的增長曲線的斜率實(shí)際上反映了后冰脊拖曳系數(shù)Cr的大小，拖曳系數(shù)Cr明顯被分為兩部分，以L/H=10 為分界。在L/H≤10 時(shí)，后冰脊受到前冰脊遮掩效果明顯，此時(shí)拖曳系數(shù)Cr變化均勻，隨L/H增大逐漸由負(fù)到正。在L/H≥10 時(shí)，拖曳力隨速度平方的增長曲線十分集中，說明此時(shí)后冰脊的拖曳系數(shù)Cr十分接近，基本與L/H無關(guān)。為更清晰地體現(xiàn)這種變化，需要進(jìn)一步分析F2隨L/H的變化規(guī)律。

3.2.2 無量綱間距的影響

圖4 給出了在不同入水深度以及流速下，后冰脊拖曳力F2隨無量綱間距L/H的變化規(guī)律。從圖4 的總體趨勢(shì)可以看出，隨著冰脊間距的增大，后冰脊的拖曳力F2不斷增大直至平穩(wěn)，它的變化分為3 個(gè)階段：當(dāng)L/H=2.5～10 時(shí)，盡管L/H增加較小，F(xiàn)2增長變化十分明顯，稱為迅速增加階段；當(dāng)L/H=10～30時(shí)，F(xiàn)2隨著L/H的增長趨勢(shì)相對(duì)較緩，稱為緩慢增加階段；當(dāng)L/H≥30 時(shí)，隨著L/H的持續(xù)增加，F(xiàn)2相對(duì)平穩(wěn)，變化較小，稱為穩(wěn)定階段。

在迅速增加階段，拖曳力受L/H影響較大，且入水深度H越大F2增長越快，此過程中會(huì)出現(xiàn)反向拖曳力到正向拖曳力的變化情況。如果定義0 值拖曳力所對(duì)應(yīng)的無量綱間距L/H值為臨界間距值，則臨界間距值主要在5～7.5 范圍內(nèi)；且隨著入水深度的增加，臨界間距值呈現(xiàn)出減小趨勢(shì)，但是臨界間距值的變化幾乎與流速無關(guān)。在臨界間距范圍內(nèi)，由于尾流漩渦場的回流作用，冰脊間流場趨近于靜止，因此后冰脊拖曳力F2接近0 值。在緩慢增加階段，小流速情況下，F(xiàn)2隨L/H的增長已經(jīng)十分緩慢；但是流速較大時(shí)，冰脊拖曳力總值較大，拖曳力受間距影響依然較大。在平緩階段，后冰脊拖曳力F2增長緩慢趨于平穩(wěn)，但是遮掩效應(yīng)依然存在，后冰脊拖曳力F2逐漸接近前冰脊的拖曳力F1。當(dāng)L/H超過70 時(shí)，F(xiàn)2的相對(duì)變化減小到不超過5%，可認(rèn)為后冰脊已經(jīng)離開前冰脊的尾流遮掩區(qū)[18]，不再考慮兩個(gè)冰脊間的相互影響。

3.3 遮掩效應(yīng)分析

由上述分析可知，后冰脊拖曳力主要受到流體速度、冰脊入水深度和冰脊間距3 個(gè)因素的影響，而前后冰脊的冰-水拖曳力比值則能直觀地反映出前冰脊尾流遮掩效應(yīng)的強(qiáng)度和范圍。圖5 給出了物理模型實(shí)驗(yàn)中前后拖曳力比值與無量綱間距L/H的關(guān)系。

圖5 中的誤差線表示了流速的影響，可以看到不同流速對(duì)應(yīng)的拖曳力比值變化非常有限，而冰脊入水深度和冰脊間距的影響較大。將前后冰脊的拖曳力比值F2/F1與無量綱間距的變化關(guān)系進(jìn)行擬合，發(fā)現(xiàn)兩者之間有明顯的指數(shù)函數(shù)關(guān)系：

圖5 中還給出了其他公式的計(jì)算結(jié)果作為比較。點(diǎn)劃線表示的是Tsamados 等[8]在冰-氣和冰-水界面拖曳系數(shù)參數(shù)化公式中所使用的指數(shù)衰減函數(shù)，即F2/F1=1-exp(-sL/H)，其中s為無量綱間距L/H的系數(shù)；虛線所表示的是Lu 等[4]在冰-水拖曳系數(shù)參數(shù)化公式中所使用的冪衰減函數(shù)，即F2/F1=[1-(H/L)1/2]2。當(dāng)L/H≤20 時(shí)，指數(shù)衰減函數(shù)TS2014 表示的遮掩作用影響較弱；在L/H＞20 時(shí)，指數(shù)衰減函數(shù)TS2014 已經(jīng)接近1，則遮掩效應(yīng)帶來的衰減不再存在。且指數(shù)衰減函數(shù)TS2014 公式來自陸地防護(hù)林實(shí)驗(yàn)，然而陸地防護(hù)林與海冰的真實(shí)情況存在較大差異，高密度的防護(hù)林也會(huì)允許風(fēng)場通過[6]，因此指數(shù)衰減函數(shù)TS2014的不準(zhǔn)確性主要體現(xiàn)在表示的遮掩影響較弱和遮掩影響的范圍較小。且指數(shù)函數(shù)中s的選取沒有統(tǒng)一確定，Hanssen-Bauer 和Gjessing[6]、Tsamados 等[8]選取s=0.18；而Lüpkes 等[7]在使用時(shí)選取了s=0.5。冪衰減函數(shù)LU2011 表示的遮掩影響相對(duì)較大一些，但是仍然可以發(fā)現(xiàn)與實(shí)驗(yàn)結(jié)果有較大出入。在L/H≤20 時(shí)，LU2011遮掩影響相對(duì)偏小，而在L/H≥20 時(shí)，遮掩影響又相對(duì)偏大。雖然整體遮掩影響的作用范圍相對(duì)指數(shù)衰減函數(shù)TS2014 較大，但遮掩范圍過大，在L/H=500 時(shí)僅趨近于0.91。

圖5 拖曳力比值F2/F1 隨無量綱間距L/H 的變化Fig.5 The drag force ratio F2/F1 changes with the dimensionless distance L/H

相比前兩種衰減函數(shù)，本文結(jié)果顯示的衰減函數(shù)可以更好地描述多冰脊之間的遮掩效應(yīng)。尤其在小間距情況下，遮掩效應(yīng)最為明顯，且在L/H≤5 時(shí)出現(xiàn)反向拖曳力，即F2/F1＜0，而TS2014 和LU2011 的衰減函數(shù)均不能描述出這一特征，因?yàn)樗鼈儧]有考慮到由于旋渦回流而產(chǎn)生的反向拖曳力會(huì)對(duì)冰脊總拖曳力產(chǎn)生的影響。而在密集度較高且浮冰較小的情況下，尾流的遮掩作用最強(qiáng)[10]，所以小間距條件在實(shí)際的海冰情況中顯得尤為重要。且整體來看，式（1）所表示的遮掩影響更加符合實(shí)際情況。根據(jù)后冰脊拖曳力隨無量綱間距變化（圖4）可知，當(dāng)L/H＞70 時(shí)，可以認(rèn)為后冰脊已經(jīng)離開尾流區(qū)，因此遮掩式（1）的適用范圍為L/H=0～70。需要注意的是，隨著間距的增大，式（1）的拖曳力比值F2/F1在超出尾流區(qū)后穩(wěn)定在0.86，并未達(dá)到理想值1，與LU2011 冪衰減函數(shù)類似。這是由于在進(jìn)行實(shí)驗(yàn)室物理模型實(shí)驗(yàn)時(shí)，冰脊的入水深度與實(shí)驗(yàn)水深在同量級(jí)，流場的演變會(huì)受到水深的明顯影響；而且實(shí)驗(yàn)水槽水頭會(huì)因水流的影響而產(chǎn)生損失從而導(dǎo)致流速衰減，影響拖曳力值的測量值[19]。該問題可以通過開展大水深條件下的物理模擬或者數(shù)值模擬研究來解決，最終得到無關(guān)水深條件的遮掩規(guī)律。

4 結(jié)論

本文以探究冰-水界面的多冰脊遮掩效應(yīng)為研究目標(biāo)，開展了測量前冰脊遮掩對(duì)后冰脊拖曳力影響的物理模型實(shí)驗(yàn)，得到了前后冰脊冰-水拖曳力在尾流遮掩情況下的變化規(guī)律，并總結(jié)出以下結(jié)論。

前冰脊拖曳力不受遮掩效應(yīng)的影響，其拖曳力與流速和冰脊入水深度等變量的關(guān)系與單冰脊情況相同。后冰脊受到遮掩影響，拖曳力與水流速度的平方呈線性關(guān)系。當(dāng)L/H≤7.5 時(shí)，后冰脊處在前冰脊尾流漩渦場中，從而出現(xiàn)反向拖曳力。后冰脊拖曳力由負(fù)值到0 值再到正值的過程中，拖曳力與流速平方的線性關(guān)系保持不變，而拖曳系數(shù)Cr呈現(xiàn)出先減小后增大至不變的趨勢(shì)。后冰脊拖曳力F2隨著L/H的增大而增大，0 值拖曳力對(duì)應(yīng)的臨界間距值主要集中在5～7.5 范圍內(nèi)，且隨著入水深度的增加呈現(xiàn)減小趨勢(shì)。在平緩階段L/H＞70 時(shí)，可認(rèn)為不再需要考慮尾流的衰減影響。

多冰脊間的遮掩效應(yīng)主要受到冰脊入水深度和冰脊間距的影響，與流速無關(guān)。本實(shí)驗(yàn)發(fā)現(xiàn)，前后冰脊的拖曳力比值與無量綱間距L/H之間有明顯的指數(shù)函數(shù)關(guān)系，并得到了可以更加準(zhǔn)確描述冰下流場由于冰脊尾流所產(chǎn)生的遮掩影響的衰減函數(shù)，該函數(shù)直接從冰脊拖曳的物理模型實(shí)驗(yàn)中得到，將冰脊形態(tài)和冰脊間距對(duì)流速衰減函數(shù)的影響進(jìn)行參數(shù)化表示，與現(xiàn)應(yīng)用在海冰模式中的兩種衰減函數(shù)相比，更適用于描述冰脊密集區(qū)域冰-水界面上的流速衰減過程。

作為冰脊間遮蔽作用研究的第一步，本文給出了傾角為45°的冰脊在有限水深下的實(shí)驗(yàn)結(jié)果，其意義在于給出初步規(guī)律以及驗(yàn)證研究方法的可行性。目前式（1）還不能直接應(yīng)用于模式計(jì)算，下一步將通過更多的物理模型實(shí)驗(yàn)與數(shù)值模擬的計(jì)算，得到不同類型冰脊在無限水深情況下的結(jié)果，為海冰拖曳系數(shù)參數(shù)化的完善提供支持。