何霞，楊成宇，王國(guó)榮，鐘林，蔣雯婧

恒定流速下三角型溝槽隔水管減阻性能的仿真分析

何霞a,b，楊成宇a,b，王國(guó)榮a,b，鐘林a,b，蔣雯婧a,b

（西南石油大學(xué) a.機(jī)電工程學(xué)院 b.能源裝備研究院，成都 610500）

為了提高隔水管的使用壽命和保障深水鉆井的安全性，在只考慮同一流速下的海流力、不同溝槽的布置對(duì)隔水管受力影響的工況下，利用織構(gòu)減阻的原理，將三角型溝槽織構(gòu)引入隔水管表面來(lái)起到降低隔水管繞流阻力的作用。采用大渦數(shù)值模擬，建立三角溝槽織構(gòu)隔水管力學(xué)模型，在恒定流速下開(kāi)展三角型溝槽深度和寬度幾何參數(shù)對(duì)隔水管升阻力系數(shù)、速度場(chǎng)分布、壓力系數(shù)分布、渦量分布影響規(guī)律的仿真分析，并與光滑隔水管的速度場(chǎng)分布、壓力系數(shù)分布、速度矢量場(chǎng)分布進(jìn)行對(duì)比。在南海恒定平均流速0.3 m/s的模擬工況下，當(dāng)三角型溝槽織構(gòu)溝槽個(gè)數(shù)為32時(shí)，溝槽織構(gòu)隔水管的升阻力系數(shù)均小于光滑隔水管，起到了減阻抑升的作用，且在深度為0.14 mm、寬度為1.96 mm的條件下，隔水管的阻力系數(shù)降低了65.02%，升力系數(shù)均方根降低了78.3%，減阻效果最佳。建立了一套恒定流速下三角型溝槽隔水管繞流減阻優(yōu)化設(shè)計(jì)與評(píng)價(jià)的新方法，為利用三角型溝槽織構(gòu)改善隔水管的疲勞失效、提升其使用壽命提供了理論參考。

隔水管；三角型溝槽；渦激振動(dòng)；大渦數(shù)值模擬；抑振減阻

隨著海洋油氣開(kāi)發(fā)技術(shù)的提高，油氣資源的勘采進(jìn)入了更深的海域，高流速引起的漩渦泄放會(huì)加劇隔水管振動(dòng)[1]，加劇了隔水管的疲勞損傷。渦激振動(dòng)抑制方法可以從三方面入手，分別是：通過(guò)改變質(zhì)量、隔水管結(jié)構(gòu)設(shè)計(jì)等方式重新設(shè)計(jì)隔水管參數(shù)；增加系統(tǒng)阻尼來(lái)提高隔水管系統(tǒng)的穩(wěn)定性；添加附加抑制裝置以減小振動(dòng)。用于抑制海洋工程管柱渦激振動(dòng)的方法，根據(jù)Zdravkovich[2]所做的研究可分為：管柱上表面突起物，可使得來(lái)流分離，擾亂尾部渦旋的線條裝、螺旋列板、半球體凸面等；附加裝置，如可將來(lái)流分離為許多小漩渦的多孔覆蓋物、控制桿、軸向板條等；近尾流穩(wěn)定裝置，如導(dǎo)流板、分離盤(pán)、飄帶等。

徐志輝[3]研究了國(guó)內(nèi)螺旋列板，此裝置的抑制原理是通過(guò)可改變沿管縱向的波流分離角，使得旋渦的強(qiáng)度被削弱，升力會(huì)降低，從而有效地弱化尾流渦旋，減小立管渦激振動(dòng)響應(yīng)[4]。但是螺旋列板材料的選擇會(huì)因?yàn)槊芏仍黾庸苤亓浚⑶野惭b成本較高，需要對(duì)其進(jìn)行防腐處理[5]。許天羽等[6]運(yùn)用CFD研究了導(dǎo)流板形橫剖面立管對(duì)渦激振動(dòng)的抑制，結(jié)果顯示，導(dǎo)流板兩側(cè)角度選取的不同會(huì)產(chǎn)生不同的渦激抑制效果。趙恩金等[7]研究了導(dǎo)流板對(duì)海底管線渦激振動(dòng)的影響，發(fā)現(xiàn)安裝導(dǎo)流板會(huì)使管線分離點(diǎn)前移，并隨著導(dǎo)流板的增高，周期性作用在管線上的豎直向下力逐漸增大。王曉龍等[8]研究了不同弦厚比、不同流速、不同直徑立管的導(dǎo)流板對(duì)抑制渦激振動(dòng)的影響規(guī)律，說(shuō)明了導(dǎo)流板可以隨來(lái)流方向旋轉(zhuǎn)，在渦激振動(dòng)的抑制方面高達(dá)80%，但是導(dǎo)流板結(jié)構(gòu)設(shè)計(jì)和制作復(fù)雜，正常工作周期比螺旋列板更短，安裝成本高。Roshko[9]在鈍體尾流的研究文獻(xiàn)中提出利用分離盤(pán)裝置穩(wěn)定管柱尾跡、抑制旋渦脫落。許博涵等[10]用風(fēng)洞試驗(yàn)測(cè)試了柔性分離盤(pán)，結(jié)果顯示，加裝分離盤(pán)的長(zhǎng)度為隔水管直徑的0.5～2.5倍時(shí)，最大減阻效率可達(dá)30%。1987年Lesag等[11]通過(guò)在鈍體上游放置控制桿減少其所受的拖曳力和橫向力。Sakamoto等[12]研究了在圓柱系統(tǒng)內(nèi)加控制桿來(lái)對(duì)流體力進(jìn)行控制。Lee等[13]通過(guò)試驗(yàn)添加控制桿對(duì)圓柱所受拖曳力及其流場(chǎng)變化，添加控制桿使得圓柱的拖曳力系數(shù)最大可以減少29%，通過(guò)合理布置附屬管線達(dá)到最大的減阻效果。通過(guò)對(duì)上述傳統(tǒng)的渦激振動(dòng)抑制裝置調(diào)研發(fā)現(xiàn)其有效率的抑制渦激振動(dòng)，針對(duì)不同的裝置有各自的優(yōu)缺點(diǎn)，但是傳統(tǒng)的渦激振動(dòng)抑制裝置普遍具有安裝難度大、結(jié)構(gòu)更加復(fù)雜且成本較高的缺陷，不能適應(yīng)各方向來(lái)流的海洋工況。

目前湍流邊界層減阻的研究已經(jīng)引起了國(guó)內(nèi)外學(xué)者的高度重視。NASA蘭利研究所發(fā)現(xiàn)順流向微小的溝槽表面具有降低壁面摩阻的作用。學(xué)者們對(duì)此進(jìn)行了不懈的探究，1999—2005年Yamagishi等[14-16]先后研究了在光滑圓柱表面周向布置不同深度的半圓溝槽、相同數(shù)量的深度半圓溝槽和三角溝槽及不同數(shù)量的三角形溝槽的流場(chǎng)情況，得出當(dāng)溝槽數(shù)目為32時(shí)有較佳的減阻效果，且該三角形溝槽在特定雷諾數(shù)下具有高達(dá)30%的減阻效果。2004年Lim等[17]在光滑圓柱體表面布置了O型圓環(huán)，達(dá)到了控制尾流區(qū)旋渦脫落的目的，使得布置了O型圓環(huán)的圓柱受到的阻力減小。2013年Quintavalla等[18]通過(guò)試驗(yàn)研究了布置不同形狀的溝槽的光滑圓柱的繞流情況，結(jié)果顯示，溝槽使得臨界雷諾數(shù)下降，具有較好的減阻效果。2018年王偉等[19]通過(guò)非定常數(shù)值研究方法對(duì)脊?fàn)罱Y(jié)構(gòu)表面圓柱的流場(chǎng)情況進(jìn)行了研究，結(jié)果顯示，在特定雷諾數(shù)下，阻力系數(shù)最大減少32.56%，升力系數(shù)最大減少65.41%。

傳統(tǒng)的織構(gòu)減阻常以平板為載體，將仿生織構(gòu)運(yùn)用到曲面（隔水管）的研究在國(guó)內(nèi)開(kāi)展得較少。故本文提出把三角型溝槽織構(gòu)引入隔水管表面來(lái)降低繞流阻力的新方法，采用了大渦數(shù)值模擬方法，對(duì)25組布置有不同寬度與深度的三角溝槽的隔水管的流場(chǎng)進(jìn)行分析，并與光滑隔水管流場(chǎng)參數(shù)進(jìn)行對(duì)比，論證三角型溝槽對(duì)隔水管的減阻效果，為利用三角型溝槽織構(gòu)改善隔水管的疲勞失效、提升其使用壽命提供理論參考。

1 數(shù)值計(jì)算方法

1.1 隔水管模型簡(jiǎn)化及三角溝槽建立

本文只考慮海流力、不同溝槽的布置及同一海洋流速對(duì)隔水管受力的影響。采用有限元法方法對(duì)隔水管進(jìn)行分析，并假設(shè)：（1）隔水管單元體材料特性、幾何特性均保持不變；（2）將隔水管上下端面簡(jiǎn)化為鉸支約束的固定連接；（3）隔水管為剛性體，在自重和外力的作用下不發(fā)生形變，且不考慮隔水管內(nèi)部鉆柱及鉆井液對(duì)其影響[20]。建立三角溝槽隔水管簡(jiǎn)化模型，如圖1所示，其中溝槽數(shù)量為32，隔水管直徑為=40 mm，展向長(zhǎng)度=130 mm[21]。

隔水管所受雷諾數(shù)為2×104～8×105，由邊界層公式(1)[22]計(jì)算得到邊界層厚度為0.1～1 mm，式(1)中為邊界層厚度，為隔水管展向長(zhǎng)度，Re為展向長(zhǎng)度為的隔水管所受的雷諾數(shù)，故選取5個(gè)深度（0.14、0.28、0.42、0.56、0.7 mm），并根據(jù)溝槽的間隙選擇5個(gè)寬度（0.98、1.96、2.86、3.38、3.92 mm），如圖2—3所示，并對(duì)以上參數(shù)進(jìn)行正交數(shù)值仿真計(jì)算。試樣代號(hào)如表1所示。

圖2 三角溝槽深度的布置

圖3 三角溝槽寬度的布置

表1 不同寬度高度的正交計(jì)算

Tab.1 Orthogonal calculation of different width and height

1.2 控制方程與湍流模型

大渦數(shù)值模擬控制方程與Smagorinsky渦黏模型見(jiàn)公式(2)—(7)[23]。

1.3 網(wǎng)格劃分及邊界條件

針對(duì)三角溝槽隔水管計(jì)算域，為保證本文大渦數(shù)值模擬方法的精確可靠性，計(jì)算域采用全局六面體網(wǎng)格，如圖4—5所示，通過(guò)Fluent仿真軟件，設(shè)置模擬時(shí)間為流體質(zhì)點(diǎn)經(jīng)過(guò)3～5倍計(jì)算域長(zhǎng)度所需的時(shí)間。本模型中圓柱壁面邊界層加密設(shè)置為40層，沿壁面法線方向與圓柱壁面相鄰的塊區(qū)域的邊界層生長(zhǎng)率為1.12，使得該網(wǎng)格滿足壁面無(wú)量綱高度因子+=1，展示模型為20.42。

圖4 W2H0.42三角溝槽全局網(wǎng)格

圖5 W2H0.42三角溝槽隔水管近壁面網(wǎng)格

根據(jù)幾何相似和雷諾相似的原理，在雷諾數(shù)相等的情況下，將數(shù)值帶入公式(8)，式中為流速，為密度（1 000 kg/m3），0為特征長(zhǎng)度（此處對(duì)應(yīng)隔水管的外徑0），為動(dòng)力黏度（=1.01×10?3pa·s），0為南海海域平均流速（0=0.3 m/s）[17]，0為附加浮力塊的隔水管外徑（0=0.7 m），通過(guò)計(jì)算得到雷諾數(shù)為2.1×105。本文光滑隔水管模型=40 mm，故在仿真計(jì)算中進(jìn)口速度設(shè)置為=5 m/s。

隔水管是柔性管柱，在海洋工況中承受梯度洋流，近海平面處流速較大且伴隨有風(fēng)浪沖擊，流體沖擊圓柱后無(wú)邊界的自由擴(kuò)散。故為了模擬海洋工況，設(shè)定左側(cè)進(jìn)口邊界面（Inlet）為速度入口，速度方向沿著X軸，并在Z軸方向存在重力，使得流體自然流動(dòng)；右側(cè)出口邊界（Outlet）面為自由出流，為壓力出口且壓力為0；除了左右兩側(cè)邊界的其余四周邊界面（Symmetry）均設(shè)定為對(duì)稱(chēng)壁面邊界，且圓柱壁面（Wall）設(shè)定為固壁無(wú)滑移的邊界條件，如圖6所示，具體邊界條件如表2所示。

表2 邊界條件

Tab.2 Boundary conditions

2 溝槽隔水管的流場(chǎng)分析

2.1 升阻力系數(shù)

升阻力系數(shù)主要反映在該工況下溝槽隔水管瞬時(shí)受力行為特征，其表達(dá)式見(jiàn)式(9)—(10)[24]。

式中：d表示阻力，N；l表示升力，N；為海流密度，kg/m3；為圓柱體外徑，m；為圓柱體展向長(zhǎng)度，m。

如圖7a所示，雷諾數(shù)=2.1×105時(shí)，光滑隔水管阻力系數(shù)為1.155，升力系數(shù)均方根為0.397[25]。在圖7b—f中，所有溝槽隔水管的升阻力系數(shù)均小于光滑隔水管，說(shuō)明三角溝槽織構(gòu)對(duì)隔水管起到了減阻抑升的效果。在同一溝槽深度下，隨著溝槽寬度的減小，阻力系數(shù)均存在先降后升的變化趨勢(shì)；升力系數(shù)均方根在溝槽深度0.28 mm前后，分別呈現(xiàn)出先降后升、一直升高的變化趨勢(shì)。由此可以說(shuō)明，隨著溝槽深度、寬度的增加，減阻效果降低，在溝槽深度為0.14 mm且寬度為0.98～1.96 mm時(shí)存在升阻力最小值，而減阻抑升效果最佳的點(diǎn)就在上述區(qū)間內(nèi)。在25組數(shù)據(jù)中，溝槽30.14組的阻力系數(shù)為0.404、升力系數(shù)均方根為0.086，其減阻效率為65.02%，抑升效率為78.3%，減阻效果最佳；通過(guò)計(jì)算得到該組斯托魯哈爾數(shù)=0.301，符合斯托哈爾數(shù)在粗糙圓柱面上的數(shù)值（0.1～0.47）[26]，說(shuō)明模擬結(jié)果是可靠的。在圖7e和圖7f中，有30.7、40.56、40.7 3個(gè)案例在計(jì)算過(guò)程中產(chǎn)生回流，最終導(dǎo)致計(jì)算終止。

圖7 光滑隔水管升阻力系數(shù)及溝槽隔水管升阻力系數(shù)的變化規(guī)律

2.2 速度與壓力分布

由溝槽隔水管的升阻力系數(shù)分析可得，溝槽深度為0.14 mm時(shí)，減阻效果最顯著，故對(duì)該深度下的溝槽隔水管進(jìn)行速度場(chǎng)分析。從圖8a—e可以看出，在同一溝槽深度0.14 mm下，隨著溝槽寬度的減小，溝槽隔水管的尾流區(qū)域從溝槽壁面處分離時(shí)，位置逐漸向前移動(dòng)，尾流區(qū)域沿來(lái)流法線方向上的寬度逐漸增大，且沿著來(lái)流方向的中心線上下交替排列，形成了規(guī)則、周期性的渦街排列，當(dāng)寬度減少至0.42 mm后，尾流區(qū)域開(kāi)始形成面積大而紊亂的速度分布區(qū)域。理論上阻力系數(shù)越小，克服阻力所損失的能量越少，對(duì)應(yīng)的速度最大值應(yīng)該增加[27]。從上文中的升阻力系數(shù)變化規(guī)律可知，前4組溝槽的升阻力系數(shù)依次逐漸減少。在圖8中，前4組溝槽中溝槽壁面分離位置處向外的速度變化逐漸緩慢，能量損失逐漸減少，與升阻力系數(shù)的變化規(guī)律吻合，證明同一深度下溝槽隔水管的速度分布規(guī)律滿足該原理。

由于第2組溝槽中，5組數(shù)據(jù)中20.28的升阻力系數(shù)最小，故比較同一寬度2.86 mm下不同溝槽深度的速度場(chǎng)分布情況。如圖9a—e所示，隨著溝槽深度的增加，尾流區(qū)域沿著來(lái)流的法線方向上的寬度逐漸增加，尾流區(qū)域從溝槽壁面分離時(shí)位置同樣是逐漸向前移動(dòng)，當(dāng)深度增加至0.42 mm后，尾流區(qū)域開(kāi)始形成面積大而紊亂的速度分布區(qū)域。由上文中升阻力系數(shù)變化規(guī)律可知，圖9中5組溝槽的升阻力系數(shù)依次增大，在尾流區(qū)域沒(méi)有出現(xiàn)紊亂的前2組溝槽速度分布的變化規(guī)律與升阻力變化規(guī)律吻合，也符合阻力系數(shù)越小，克服阻力所損失的能量越少的原理。

圖8 溝槽深度為0.14 mm時(shí)不同寬度的速度分布云圖

圖9 溝槽寬度為2.86 mm時(shí)不同深度的速度分布云圖

通過(guò)對(duì)不同寬度、深度的溝槽隔水管速度場(chǎng)分布分析可知，在上述溝槽深度、寬度的范圍內(nèi)，隨著溝槽深度的減小、寬度的增加，尾流區(qū)域沿方向的寬度逐漸增加，尾流區(qū)域處的漩渦形成的渦街排列會(huì)更有規(guī)則性、周期性，減阻效果會(huì)更加明顯。由于流體會(huì)在近壁面產(chǎn)生較大的速度梯度，當(dāng)靠近壁面的速度在法線上的梯度為0時(shí)，會(huì)出現(xiàn)流動(dòng)分離，且分離點(diǎn)越靠前，尾流區(qū)域越大[28]，隨著溝槽深度的增加、寬度的減小，溝槽織構(gòu)會(huì)使流動(dòng)分離點(diǎn)逐漸靠前。當(dāng)溝槽深度大于0.42 mm、寬度小于0.98 mm時(shí)，尾流區(qū)域從溝槽壁面分離時(shí)的位置過(guò)于靠前，導(dǎo)致尾流區(qū)域?qū)?huì)出現(xiàn)大面積紊亂的速度分布區(qū)域，對(duì)隔水管的減阻起到負(fù)作用。

溝槽隔水管阻力減小的直接流場(chǎng)因素就是力的分布，壓差產(chǎn)生的阻力是圓柱繞流流場(chǎng)中阻力產(chǎn)生的主要原因[29]。如圖10所示，隨著溝槽寬度的減小，前4組溝槽隔水管后方尾流區(qū)域的壓力（壓力分布云圖上可通過(guò)顏色的深淺判斷壓力的變化規(guī)律）逐漸減小，前后的壓差逐漸增大，阻力隨之逐漸減小，而40.14組溝槽隔水管后方的壓力急增，其前后的壓差也急劇減少，阻力隨之增加，與上文中同一深度0.14 mm、不同寬度的溝槽隔水管阻力系數(shù)的變化規(guī)律一致。

同樣，探究同一寬度2.86 mm下不同溝槽深度的流場(chǎng)分布情況。由圖11的壓力云圖分布可知，隨著溝槽深度的增加，溝槽隔水管后方的尾流區(qū)域的壓力逐漸減小，由最初的中壓區(qū)域逐漸轉(zhuǎn)變?yōu)槲擦鲌?chǎng)的低壓區(qū)壓力分布，隔水管前后的壓差也逐漸增大，阻力隨之減小，與上文中同一寬度2.86 mm、不同深度的溝槽隔水管阻力系數(shù)的變化規(guī)律一致。圖10—11所示壓力分布云圖中的壓力均為靜壓。

圖10 溝槽深度為0.14 mm時(shí)不同寬度的壓力分布云圖

圖11 溝槽寬度為2.86 mm時(shí)不同深度的速度分布云圖

2.3 渦量分布規(guī)律

圖12為第1組溝槽隔水管在=2.1×105時(shí)，=0時(shí)，圓柱中心截面上的渦量分布云圖，從圖中可以看出，溝槽尾流區(qū)域上下兩側(cè)渦量較大，形成交替脫落的漩渦，隨著溝槽寬度的減小，漩渦脫落的頻率增加，渦振漩渦強(qiáng)度增加且隨著流體質(zhì)點(diǎn)的流動(dòng)逐漸減弱。由此可以看出，布置有溝槽結(jié)構(gòu)的隔水管能夠有效改變其尾流的漩渦脫落頻率，改變溝槽隔水管的振蕩頻率，降低渦激振動(dòng)對(duì)隔水管的疲勞損壞。

圖12 深度為0.14 mm時(shí)不同寬度的溝槽隔水管渦量分布云圖

3 三角溝槽減阻效果的評(píng)價(jià)

對(duì)不同參數(shù)的溝槽隔水管的分析可以得出，在此25組的織構(gòu)布置中，30.14是減阻抑升性能最佳的溝槽布置，下面將對(duì)光滑隔水管和30.14的流場(chǎng)各物理量進(jìn)行分析，進(jìn)一步判斷溝槽的布置對(duì)哪些因素造成了影響。

從圖13可以看出，在隔水管上布置了溝槽的結(jié)構(gòu)體速度的最大值更大、壓力值更小，說(shuō)明動(dòng)能向壓力能的轉(zhuǎn)換量減小，尾流近壁面速度存在接近于0的點(diǎn)，速度為0的區(qū)域是由于漩渦的脫落造成大渦破碎成小渦，從而將動(dòng)量傳遞和消耗；光滑隔水管的尾流區(qū)域呈現(xiàn)放射狀的放大趨勢(shì)，而溝槽隔水管的尾流區(qū)域呈現(xiàn)收縮狀，也可以看出溝槽隔水管尾跡區(qū)域的寬度明顯減小，且溝槽隔水管速度回流的區(qū)域沿著流體質(zhì)點(diǎn)運(yùn)動(dòng)方向向后移動(dòng)，速度分離點(diǎn)向后移動(dòng)了約23°，說(shuō)明該溝槽結(jié)構(gòu)能夠控制尾流區(qū)域達(dá)到延緩流動(dòng)分離的目的。

圖13 速度分布對(duì)比圖

圖14可以從顏色數(shù)值大小的分布來(lái)分析，首先光滑隔水管的最小壓力系數(shù)為?2.5，溝槽隔水管的最小壓力系數(shù)為?4，最小壓力值出現(xiàn)在隔水管兩側(cè)的上下端面，說(shuō)明溝槽隔水管兩側(cè)速度更大，保留的動(dòng)能更多；光滑隔水管的較小壓力數(shù)值為?2.2～ ?0.8，大面積地出現(xiàn)在尾流區(qū)域，而溝槽隔水管的尾流區(qū)的較小壓力數(shù)值為?1.6～0。在尾流區(qū)域溝槽隔水管的壓力明顯大于光滑隔水管，使得溝槽隔水管的前后壓差遠(yuǎn)遠(yuǎn)小于光滑隔水管，壓差阻力是造成隔水管橫向振動(dòng)和變形大位移的主要原因，壓差阻力越小，橫向振動(dòng)越小，隔水管變形也越小。

在圖15所示的速度矢量圖中加入速度流線，流線代表流體質(zhì)點(diǎn)在隔水管中心面=0的軌跡，光滑隔水管上紅色矩形框內(nèi)的速度近似為0，光滑隔水管的分離點(diǎn)在90°左右；而溝槽隔水管在90°之后的壁面上發(fā)生了明顯的分離再附現(xiàn)象，漩渦產(chǎn)生的位置明顯向后移動(dòng)，而分離點(diǎn)也向后移動(dòng)，且前者流線的流向是放開(kāi)向外，沿著隔水管壁面呈放大的趨勢(shì)，后者的流線沿著隔水管90°切線方向向=0的中心線偏移，呈現(xiàn)收口的趨勢(shì)，說(shuō)明尾流區(qū)域溝槽隔水管的速度小于光滑隔水管。

綜上可知，在速度分布上溝槽隔水管的最大速度更大，壓力最小值更小，溝槽隔水管尾流區(qū)的再附現(xiàn)象和二次分離點(diǎn)沿著流動(dòng)方向向后延伸，三角溝槽壁面增加了邊界層的擾動(dòng)，使得溝槽凹槽內(nèi)形成了旋轉(zhuǎn)的漩渦。該漩渦方向與流體流動(dòng)方向相反，溝槽凸起處呈現(xiàn)較低的流速，將流動(dòng)分離點(diǎn)延后，使溝槽隔水管的尾流區(qū)域更小，壓差阻力也更小，使得整體的升力阻力得到控制，起到了抑升減阻的效果。

圖14 壓力系數(shù)分布對(duì)比圖

圖15 隔水管速度矢量對(duì)比圖

4 結(jié)論

1）建立了三角溝槽織構(gòu)隔水管大渦數(shù)值模擬模型，整體模型采用全局六面體網(wǎng)格，并進(jìn)行了網(wǎng)格無(wú)關(guān)性驗(yàn)證，使得該網(wǎng)格滿足壁面無(wú)量綱高度因子+=1，保證網(wǎng)格的精度，為利用三角型溝槽織構(gòu)改善隔水管的疲勞失效、提升其使用壽命提供了理論參考。

2）采用大渦數(shù)值模擬方法，開(kāi)展了5種深度和5種寬度溝槽布置的隔水管升阻力系數(shù)數(shù)值計(jì)算，得到了在=2.1×105時(shí)，隨著溝槽深度的增加，阻力系數(shù)逐漸增大。在選取的25組數(shù)據(jù)中，30.14在升阻力系數(shù)的大小上表現(xiàn)最好，并隨著溝槽寬度的減小，升阻力系數(shù)呈現(xiàn)先減小再增大的趨勢(shì)。與光滑隔水管的升阻力系數(shù)數(shù)值對(duì)比，其阻力系數(shù)降低了65.02%，升力系數(shù)均方根值降低了78.3%，這種溝槽的布置大大減小了隔水管上所受的阻力與升力。

3）光滑隔水管與30.14溝槽隔水管的對(duì)比結(jié)果顯示，在30.14溝槽隔水管的布置時(shí)，速度場(chǎng)變大，隔水管尾流區(qū)的壓力顯著增大，前后壓差減小導(dǎo)致壓差阻力減小。對(duì)比光滑隔水管，溝槽隔水管的分離點(diǎn)靠后，尾流區(qū)域的寬度大大減小，進(jìn)一步證實(shí)了溝槽織構(gòu)布置在隔水管上起到了減阻抑升的可行性。

[1] 梁健, 劉成文, 李兆敏. 帶整流罩隔水管流場(chǎng)特性的大渦模擬研究[J]. 力學(xué)季刊, 2020, 41(1): 197-205.

LIANG Jian, LIU Cheng-wen, LI Zhao-min. Large Eddy Simulation Study on Flow Field Characteristics of Risers with Fairings[J]. Chinese Quarterly of Mechanics, 2020, 41(1): 197-205.

[2] ZDRAVKOVICH M M. Review and Classification of Various Aerodynamic and Hydrodynamic Means for Sup-pre-ssing Vortex Shedding[J]. Journal of Wind Engineering and Industrial Aerodynamics, 1981, 7(2): 145-189.

[3] 徐志輝. 螺旋列板在海底管道上的國(guó)產(chǎn)化應(yīng)用[J]. 化工設(shè)計(jì)通訊, 2017, 43(4): 167-168.

XU Zhi-hui. Localization of Spiral Plate in Submarine Pipeline[J]. Chemical Engineering Design Communica-tions, 2017, 43(4): 167-168.

[4] 楊加棟, 張曉靈, 杜寶銀, 等. 螺旋列板——深水立管渦激振動(dòng)抑制裝置[J]. 海洋技術(shù), 2010, 29(4): 88-92, 116.

YANG Jia-dong, ZHANG Xiao-ling, DU Bao-yin, et al. Helical Stakes—Vortex Induced Vibration Suppression Devices of Deepwater Risers[J]. Ocean Technology, 2010, 29(4): 88-92, 116.

[5] 鄭斌, 袁巨龍, 趙萍, 等. 變曲率溝槽精密球研磨加工優(yōu)化實(shí)驗(yàn)研究[J]. 表面技術(shù), 2017, 46(2): 214-219.

ZHENG Bin, YUAN Ju-long, ZHAO Ping, et al. Experi-ment of Optimized Grinding of Precision Ball with Variable-radius Groove[J]. Surface Technology, 2017, 46(2): 214-219.

[6] 許天羽, 鄭婷婷, 胡建新. 基于CFD的導(dǎo)流板形海洋立管渦激抑制研究[J]. 浙江海洋學(xué)院學(xué)報(bào)(自然科學(xué)版), 2011, 30(6): 533-537.

XU Tian-yu, ZHENG Ting-ting, HU Jian-xin. Study on Vortex Suppression of Guide-Vane Style Riser in Ocean Current Based CFD[J]. Journal of Zhejiang Ocean Univer-sity (Natural Science), 2011, 30(6): 533-537.

[7] 趙恩金, 拾兵, 曹坤. 導(dǎo)流板對(duì)海底管線渦激振動(dòng)的影響[J]. 哈爾濱工程大學(xué)學(xué)報(bào), 2016, 37(3): 320-325.

ZHAO En-jin, SHI Bing, CAO Kun. Influence of Reflec-tors on Vortex-Induced Vibration of Subsea Pipelines[J]. Journal of Harbin Engineering University, 2016, 37(3): 320-325.

[8] 王曉龍, 王刻強(qiáng), 胡成. 導(dǎo)流板渦激振動(dòng)抑制模擬分析[J]. 遼寧化工, 2019, 48(11): 1115-1118.

WANG Xiao-long, WANG Ke-qiang, HU Cheng. Simula-tion Analysis on Vortex-Induced Vibration Suppression of Guide Plate[J]. Liaoning Chemical Industry, 2019, 48(11): 1115-1118.

[9] ROSHKO A. On the Wake and Drag of Bluff Bodies[J]. Journal of the Aeronautical Sciences, 1955, 22(2): 124-132.

[10] 許博涵, 王嘉松, 梁盛平. 柔性分離盤(pán)抑制圓柱渦激振動(dòng)的風(fēng)洞試驗(yàn)[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2017, 32(4): 470-476.

XU Bo-han, WANG Jia-song, LIANG Sheng-ping. Wind Tunnel Experiment of VIV Control on Circular Cylinder with Flexible Splitter Plate[J]. Chinese Journal of Hydro-dynamics, 2017, 32(4): 470-476.

[11] LESAGE F, GARTSHORE I S. A Method of Reducing Drag and Fluctuating Side Force on Bluff Bodies[J]. Journal of Wind Engineering and Industrial Aerodyna-mics, 1987, 25(2): 229-245.

[12] SAKAMOTO H, HANIU H. Optimum Suppression of Fluid Forces Acting on a Circular Cylinder[J]. Journal of Fluids Engineering, 1994, 116(2): 221-227.

[13] LEE S J, LEE S I, PARK C W. Reducing the Drag on a Circular Cylinder by Upstream Installation of a Small Control Rod[J]. Fluid Dynamics Research, 2004, 34(4): 233-250.

[14] YAMAGISHI Y, OKI M. Effect of the Number of Grooves on Flow Characteristics around a Circular Cylinderwith Triangular Grooves[J]. Journal of Visualization, 2005, 8(1): 57-64.

[15] SATO R, KOGANEI N, YAMAGISHI Y, et al. Effect of Groove Shape on Fluid Flow around a Circular Cylinder with Grooves[J]. Journal of the Visualization Society of Japan, 2008, 28-1(2): 1029.

[16] OKI M, AOKI K, NAKAYAMA Y. Effect of Grooves’ Depth for Flow Characteristics around a Circular Cylinder with Grooves[J]. Transactions of the Japan Society of Mechanical Engineers Series B, 1999, 65(631): 870-875.

[17] LIM H C, LEE S J. Flow Control of a Circular Cylinder with O-Rings[J]. Fluid Dynamics Research, 2004, 35(2): 107-122.

[18] QUINTAVALLA S J, ANGILELLA A J, SMITS A J. Drag Reduction on Grooved Cylinders in the Critical Reynolds Number Regime[J]. Experimental Thermal and Fluid Sci-e-nce, 2013, 48: 15-18.

[19] 王偉, 宋保維, 毛昭勇, 等. 脊?fàn)罱Y(jié)構(gòu)表面圓柱繞流的流體特性研究[J]. 華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2018, 46(9): 107-111.

WANG Wei, SONG Bao-wei, MAO Zhao-yong, et al. Study on Hydrodynamic Characteristics of Flow around a Cylinder with Ridged Surface Structure[J]. Journal of Huazhong University of Science and Technology (Natural Science Edition), 2018, 46(9): 107-111.

[20] 劉清友, 周守為, 姜偉, 等. 基于鉆井工況和海洋環(huán)境耦合作用下的隔水管動(dòng)力學(xué)模型[J]. 天然氣工業(yè), 2013, 33(12): 6-12.

LIU Qing-you, ZHOU Shou-wei, JIANG Wei, et al. A Dynamic Model of Marine Risers/Pipes under the Drilling Operation Condition and Sea Environment[J]. Natural Gas Industry, 2013, 33(12): 6-12.

[21] 蔣雯婧, 錢(qián)權(quán), 王國(guó)榮. V型溝槽表面隔水管減阻的大渦數(shù)值模擬研究[C]//第十屆全國(guó)流體力學(xué)學(xué)術(shù)會(huì)議. 杭州: 中國(guó)力學(xué)學(xué)會(huì), 2018: 355.

JIANG Wen-jing, QIAN Quan, WANG Guo-rong. Large Eddy Numerical Simulation Research on Drag Reduction of V-Groove Surface riser[C]//Abstracts of Papers of the Tenth National Symposium on Fluid Mechanics. Hangzhou: Chinese society of mechanics, 2018: 355.

[22] 宮武旗, 李新宏, 黃淑娟. 溝槽壁面減阻機(jī)理實(shí)驗(yàn)研究[J]. 工程熱物理學(xué)報(bào), 2002, 23(5): 579-582.

GONG Wu-qi, LI Xin-hong, HUANG Shu-juan. Experiment Study on the Mechanism of Riblets Drag Reduction[J]. Journal of Engineering Thermophysics, 2002, 23(5): 579-582.

[23] 李啟. 含運(yùn)動(dòng)邊界鈍體繞流風(fēng)場(chǎng)的大渦模擬數(shù)值算法[D]. 北京: 北京交通大學(xué), 2014.

LI Qi. A Numerical Method for Large-Eddy Simulation of the Wind Field around a Bluff Body with Moving Bouda-ries[D]. Beijing: Beijing Jiaotong University, 2014.

[24] 陳進(jìn)娥, 何順利, 劉廣峰. 我國(guó)海洋石油勘探開(kāi)發(fā)裝備現(xiàn)狀及發(fā)展趨勢(shì)[J]. 油氣藏評(píng)價(jià)與開(kāi)發(fā), 2012, 2(6): 67-71.

CHEN Jin-e, HE Shun-li, LIU Guang-feng. Current Situation and Development Trends of China Offshore Oil Exploration and Development Equipments[J]. Reservoir Evaluation and Development, 2012, 2(6): 67-71.

[25] 錢(qián)權(quán). 典型溝槽表面隔水管減阻性能研究[D]. 成都: 西南石油大學(xué), 2018.

QIAN Quan. Research on Drag Reduction Performance of Riser on Typical Grooved Surface[D]. Chengdu: Southwest Petroleum University, 2018.

[26] Strouhal V. Uber Eine Besondere Art Der Tonerre-gung[J]. Annalen Der Physic, 1878, 6: 216-251.

[27] 何霞, 蔣雯婧, 王國(guó)榮, 等. 基于大渦數(shù)值模擬的海洋隔水管繞流分析[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2020, 35(1): 106-112.

HE Xia, JIANG Wen-jing, WANG Guo-rong, et al. Flow around Analysis of Riser Based on Large Eddy Numerical Simulation[J]. Chinese Journal of Hydrodynamics, 2020, 35(1): 106-112.

[28] 田啟龍, 黃熙龍, 王嘉松. 附加附屬管的隔水管渦激振動(dòng)離散渦模擬[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2016, 31(5): 633-640.

TIAN Qi-long, HUANG Xi-long, WANG Jia-song. Vortex- Induced Vibration of the Riser with Affiliated Pipes Based on Discrete Vortex Method[J]. Chinese Journal of Hydro-dy-namics, 2016, 31(5): 633-640.

[29] 鄧欣, 朱賽健, 李保軍, 等. 隔水管尾流誤差補(bǔ)償技術(shù)研究[J]. 國(guó)外電子測(cè)量技術(shù), 2016, 35(8): 42-44, 50.

DENG Xin, ZHU Sai-jian, LI Bao-jun, et al. Research of Error Compensation for Riser Wake[J]. Foreign Elec-tro-nic Measurement Technology, 2016, 35(8): 42-44, 50.

Simulation Analysis of the Drag Reduction Performance of Triangular Groove Riser under Constant Flow Velocity

a,b,a,b,a,b,a,b,a,b

(a. College of Electromechanical Engineering, b. Energy Equipment Research Institute, Southwest Petroleum University, Chengdu 610500, China)

In order to improve the service life of the riser and ensure the safety of deep-water drilling, under the condition that only the marine current force at the same velocity and the influence of different groove arrangement on the stress of the riser are considered, the triangular groove texture is introduced into the riser surface to reduce the flow resistance of the riser by using the principle of texture drag reduction. Using large eddy numerical simulation, the mechanical model of triangular groove textured riser is established. The effects of geometric parameters of triangular groove depth and width on riser lift resistance coefficient, velocity field distribution, pressure coefficient distribution and vorticity distribution are simulated and analyzed at constant velocity. The velocity vector field distribution is compared.The results show that under the simulated condition of constant average velocity of 0.3 m/s in the South China Sea, when the number of triangular groove texture grooves is 32, the rise resistance coefficient of grooved texture riser is less than that of smooth riser, which plays the role of reducing resistance and restraining rise. When the depth is 0.14 mm and the width is 1.96 mm, the resistance coefficient of riser is reduced by 65.02% and the root mean square of lift coefficient is reduced by 78.3%. The drag reduction effect is the best. Based on this, a new method for optimal design and evaluation of flow drag reduction around triangular grooved riser at constant flow velocity is established, which provides a theoretical reference for improving the fatigue failure and service life of riser by using triangular grooved texture.

riser; triangular groove; vortex-induced vibration; large eddy numerical simulation; vibration suppression and drag reduction

TH117

A

1001-3660(2022)04-0157-10

10.16490/j.cnki.issn.1001-3660.2022.04.015

2021-03-23；

2021-10-21

2021-03-23；

2021-10-21

何霞（1976—），女，博士，副教授，主要研究方向?yàn)橛蜌庋b備潤(rùn)滑減磨。

HE Xia (1976—), Female, Doctor, Associate professor, Research focus: lubrication and wear reduction of oil and gas equipment.

王國(guó)榮（1977—），男，博士，教授，主要研究方向?yàn)樘烊粴馑衔镩_(kāi)采工藝及能源裝備。

WANG Guo-rong (1977—), Male, Doctor, Professor, Research focus: natural gas hydrate mining technology and energy equipment.

何霞, 楊成宇, 王國(guó)榮, 等. 恒定流速下三角型溝槽隔水管減阻性能的仿真分析[J]. 表面技術(shù), 2022, 51(4): 157-166.

HE Xia, YANG Cheng-yu, WANG Guo-rong, et al. Simulation Analysis of the Drag Reduction Performance of Triangular Groove Riser under Constant Flow Velocity[J]. Surface Technology, 2022, 51(4): 157-166.

責(zé)任編輯：萬(wàn)長(zhǎng)清