鄧露 劉建祥

摘 ? 要：柱體繞流是海洋工程領(lǐng)域的重要問題. 流體繞經(jīng)海洋平臺所產(chǎn)生的周期性渦脫會造成平臺往復運動，這將加劇平臺系泊結(jié)構(gòu)疲勞損傷，降低結(jié)構(gòu)疲勞壽命. 為研究多柱式浮式風機平臺繞流的力學特性，厘清尾流之間的干擾機理，采用Delayed Detached Eddy Simulation（DDES）方法分析了半潛式浮式風機平臺在不同來流角和流速下的力學特性，從相干結(jié)構(gòu)層面研究了尾流干擾機理，并分析了立柱尾流之間的空間相關(guān)性. 結(jié)果表明：下游立柱的阻力系數(shù)平均值隨流速變化而有較大波動;流場相干結(jié)構(gòu)以流向渦和發(fā)夾渦為主;相干結(jié)構(gòu)之間的相互作用是導致尾流干擾的原因;上下游立柱尾流具有不同的相干結(jié)構(gòu)，但兩者仍有較強的空間相關(guān)性. 研究結(jié)果可為深入理解多柱式海洋平臺繞流現(xiàn)象提供理論參考.

關(guān)鍵詞：半潛式平臺;多柱繞流;相干結(jié)構(gòu);尾流干擾;空間相關(guān)性

中圖分類號：O351.2 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

Abstract：Flow past cylinders is an important issue in ocean engineering. Periodical vortex shedding occurring around the offshore platform in current leads to reciprocating motion of the offshore platform， which may aggravate the fatigue damage and thus cause the reduction of fatigue life of the moorings. In order to study the mechanical characteristics of flow past multi-column floating wind turbine platform and understand the mechanism of wake interference clearly， numerical simulations were performed using Delayed Detached Eddy Simulation （DDES）. The effects of headings and current velocities on the mechanical characteristics of the semi-submersible floating wind turbine platform were investigated. The mechanism of wake interference was studied from the level of coherent structures， and spatial correlation between different wakes behind the columns was also investigated. The results demonstrate that： （1） the average drag force coefficient of downstream column has a large fluctuation with the change of current velocity; and streamwise vortex and hairpin vortex are the main coherent structures; （2） the interaction between coherent structures results in wake interference; （3） the upstream wake and the downstream wake have strong spatial correlation， but their coherent structures are different. The research results provide a theoretical reference for understanding the phenomenon of flow past multi-column offshore platform.

Key words：semi-submersible platform;flow past multiple cylinders;coherent structure;wake interference;spatial correlation

日益增長的能源需求促進了海上風電行業(yè)的發(fā)展，半潛式浮式風機平臺適用水深范圍廣，是一種開發(fā)具有儲量豐富、風速大等優(yōu)點的深水風能資源的重要海洋平臺[1]. 長期以來，柱體繞流一直是海洋工程領(lǐng)域的重點問題[2]. 流體流經(jīng)浮式風機平臺的柱體結(jié)構(gòu)時，柱體后方周期性脫落的漩渦會導致平臺大幅往復運動，這將加劇平臺系泊結(jié)構(gòu)的疲勞損傷，甚至導致系泊結(jié)構(gòu)疲勞破壞[3]. 因此開展半潛式浮式風機平臺繞流研究有著重要的工程意義.

在單柱繞流方面，國內(nèi)外學者針對單柱繞流的力學特性、流場特性及其影響因素等已進行了大量的研究，并取得了豐碩的成果[4-7]. 在雙柱和多柱式結(jié)構(gòu)繞流方面，賈曉荷和劉樺[8]采用大渦模擬方法分析了并列或串列排布的雙圓柱的繞流特性，結(jié)果表明串列排布時，上游圓柱對下游圓柱受力有影響，并列排布時，兩柱的阻力系數(shù)均值基本相等. 鄭宇華和顧杰[9]的研究表明，間距比（雙柱中心的間距與柱寬的比值）影響并列雙柱開縫下方的漩渦個數(shù)、形態(tài)等. Kitagawa T和Ohtab H[10]采用大渦方法模擬了串列雙圓柱的繞流特性，結(jié)果表明當間距比小于某一臨界值時，只有下游圓柱有渦脫現(xiàn)象. 李聰洲等[11]研究了串列雙柱繞流的力學特性，結(jié)果表明存在一個臨界間距比使得圓柱的升力和阻力發(fā)生跳躍性變化. 劉為民等[12]的研究表明，0°和45°來流時，陣列四立柱柱群中立柱的渦脫頻率大于單柱的渦脫頻率，渦脫模式表現(xiàn)為2S模式. 劉正浩等[13]的研究表明，來流角對陣列四柱式海洋平臺的立柱的升力和阻力系數(shù)影響較大，不同來流角下平臺底部的流場不同.

雙柱和多柱式結(jié)構(gòu)繞流時，柱體尾流之間的相互干擾會導致立柱的力學特性和流場結(jié)構(gòu)更復雜[11]. 然而，以上針對雙柱和多柱式結(jié)構(gòu)繞流的研究，主要分析了流速、來流角和間距比等因素對結(jié)構(gòu)的受力、流場形態(tài)和渦脫模式等的影響，而深入分析流場結(jié)構(gòu)以及尾流之間干擾機理的研究極少. 因此，本文采用DDES方法對不同來流角和流速下多柱式半潛式浮式風機平臺的繞流進行了三維數(shù)值模擬，分析了平臺立柱的阻力系數(shù)平均值、升力系數(shù)頻率等力學指標，從流場的相干結(jié)構(gòu)層面研究了尾流的干擾機理，并分析了尾流的空間相關(guān)性. 研究結(jié)果闡釋了多柱式結(jié)構(gòu)繞流的尾流相互干擾機理，為深入理解多柱式結(jié)構(gòu)繞流現(xiàn)象提供了理論依據(jù).

1 ? 數(shù)值建模

1.1 ? DDES方法

DDES方法是一種改進的RANS與LES相結(jié)合的方法，其在模擬流體分離流動時十分有效[16]，已得到廣泛的應用. 該方法具有多個優(yōu)點，如不需要精細的網(wǎng)格也能得到較好的結(jié)果，適合復雜的結(jié)構(gòu)等. 因此，本文采用該方法進行研究.

1.2 ? 平臺模型

選取美國國家可再生能源實驗室設(shè)計的半潛式浮式風機平臺作為分析對象，平臺的立柱（包括邊柱和中柱）、連桿和浮箱底座的結(jié)構(gòu)布置如圖1所示（圖中括號內(nèi)為構(gòu)件編號）. 平臺的結(jié)構(gòu)尺寸，吃水深度等設(shè)計參數(shù)詳見文獻[17].

1.3 ? 計算域與網(wǎng)格劃分

流場計算域及模型平面布置如圖2所示. θ 為來流角，即中柱與邊柱2的連線與x軸負向的夾角. 流域入口邊界距離邊柱2中心為10D，出口邊界距離邊柱2中心為25D，左右邊界距離中柱中心均為10D. D為邊柱直徑.

邊界條件設(shè)定為：入口為速度入口，出口為出流邊界，頂部為自由邊界，其余邊界為對稱邊界，平臺表面采用無滑移固定壁面.

網(wǎng)格劃分示意圖如圖3. 由于平臺幾何結(jié)構(gòu)復雜，因此，將流場分成兩種類型的網(wǎng)格進行劃分，在平臺附近區(qū)域采用四面體網(wǎng)格以適應復雜的幾何結(jié)構(gòu)，在其他區(qū)域采用六面體網(wǎng)格以加快計算速度. 為控制網(wǎng)格密度，將流域分成九個區(qū)域，在平臺附近和尾流區(qū)劃分密集網(wǎng)格，其他區(qū)域則劃分較稀疏的網(wǎng)格. 根據(jù)張新曙等[18]的建議，經(jīng)過驗算，并綜合考慮計算耗時成本，確定網(wǎng)格方案為：加密區(qū)（平臺四周1.4D的區(qū)域）的網(wǎng)格大小為0.12D;尾流區(qū)的初始網(wǎng)格大小為0.25D，且隨著尾流遠離平臺，網(wǎng)格逐漸稀疏.

1.4 ? 計算設(shè)置

流域采用Ansys Fluent進行求解. 時間離散采用二階隱式格式，空間離散采用二階迎風格式，壓力與速度耦合方程采用PISO方法（Pressure-Implicit with Splitting of Operators）. 時間步長取0.05 s，每個時間步長內(nèi)最大迭代次數(shù)為20，模擬時長為6 000 ～10 000 s，以保證獲得穩(wěn)定的結(jié)果. 數(shù)據(jù)處理時，選取最后約0.6倍模擬時長內(nèi)較穩(wěn)定的計算結(jié)果進行處理.

1.5 ? 計算工況

本文選取3種典型來流角：θ = 0°、30°、60°. 根據(jù)中國南海的實際情況[19]，選取12種來流速度：u = 0.32 m/s、0.42 m/s、0.53 m/s、0.63 m/s、0.74 m/s、0.85 m/s、0.95 m/s、1.06 m/s、1.16 m/s、1.27 m/s、1.37 m/s、1.48 m/s，來流方向沿x軸正向.

2 ? 計算結(jié)果及分析

2.1 ? 力學特性

2.1.1 ? 阻力系數(shù)平均值

圖4 ～ 圖7分別為中柱、邊柱1、邊柱2、邊柱3的阻力系數(shù)平均值隨來流角和流速變化的關(guān)系圖. 從這四幅圖中可知，在0°來流下，中柱的阻力系數(shù)平均值隨流速的變化而小幅波動，其余立柱的阻力系數(shù)平均值隨流速的改變而基本不變. 在30°來流下，邊柱2和邊柱3的阻力系數(shù)平均值隨流速的改變而基本不變，但中柱和邊柱1的阻力系數(shù)平均值隨流速的改變而變化較大，離散性較強. 這是因為流體流經(jīng)上游立柱后發(fā)展成為湍流，湍流本身具有較強的隨機波動性，所以，受上游尾流的影響，下游中柱和邊柱1的阻力系數(shù)平均值隨流速的變化而波動較大[12]. 在60°來流下，中柱、邊柱2和邊柱3的阻力系數(shù)平均值隨流速的改變而基本不變，而邊柱1的阻力系數(shù)平均值出現(xiàn)了跳躍波動的現(xiàn)象. 這也是因為邊柱1處于下游位置，其阻力系數(shù)平均值受上游尾流的影響較大.

由圖4可知，中柱的阻力系數(shù)平均值隨來流角的增大而減小. 阻力是漩渦脫落引起的立柱前后表面的壓力差和立柱表面的摩擦力共同作用所形成的. 隨著來流角的改變，立柱和立柱、立柱和底座在垂直于來流方向的相對距離發(fā)生了改變. 由圖2的平面布置圖可知，在0°來流時，中柱后方尾流距離立柱和底座較遠，下游柱體的阻擋作用小，上游中柱的尾流漩渦能夠充分地發(fā)展，所以阻力系數(shù)最大. 在30°來流時，邊柱1和底座1靠近中柱尾流，下游柱體的阻擋作用增大，下游中柱后方的漩渦脫落不能充分發(fā)展，故阻力系數(shù)減小. 在60°來流時，邊柱1和底座1位于中柱的正后方，完全阻擋了中柱尾流漩渦的發(fā)展，故阻力系數(shù)最小. 由圖5可知，邊柱1的阻力系數(shù)平均值隨來流角的增大呈現(xiàn)先減小后增大的趨勢. 由圖6和圖7可知，邊柱2和邊柱3的阻力系數(shù)平均值均隨來流角的增大而增大. 這些也均是因為隨著來流角的改變，柱體之間在垂直于來流方向的相對距離發(fā)生了改變，柱體與柱體之間的相互影響所導致的.

2.1.2 ? 升力系數(shù)時程曲線及頻譜分析

限于篇幅，僅給出了θ = 0°、u = 1.06 m/s時各個立柱的升力系數(shù)時程曲線，如圖8所示. 從圖8中可以看出，邊柱2、邊柱3和中柱的升力系數(shù)時程變化表現(xiàn)出較好的諧波特性，而邊柱1的升力系數(shù)變化波動較大. 邊柱2和邊柱3的升力系數(shù)變化相反，中柱的升力系數(shù)變化滯后于邊柱2. 當邊柱2的升力系數(shù)達到波峰時，中柱的升力系數(shù)稍有增大，且曲線變化趨勢發(fā)生了改變. 以圖8中t1-t2時間段為例，此時間段內(nèi)中柱的升力系數(shù)未按圖中粗劃虛線趨勢變化，而是按實線趨勢變化. 這可能是因為中柱處于邊柱2的尾流中，當上游邊柱2的升力系數(shù)達到波峰時，其后方的漩渦脫落，隨后撞擊到中柱上，導致中柱的升力系數(shù)稍微增大而改變了曲線的變化趨勢.

圖9 ～ 圖11為u = 1.06 m/s時，三種來流角下平臺各立柱的升力系數(shù)頻譜圖. 立柱后方漩渦周期性地脫落是導致升力系數(shù)周期性變化的原因. 由圖9可知，0°來流角下，中柱的頻譜除主峰外，還存在多個高頻的小峰，小峰的幅度遠小于主峰的幅度，說明中柱后方漩渦脫落以主峰頻率為主，并伴隨多個高頻成分. 邊柱1的頻譜圖中存在兩個較大的峰值，低頻對應的幅度顯著大于高頻對應的幅度，說明邊柱1后方低頻的渦脫引起了較大的升力. 邊柱2和邊柱3的頻譜圖中均只有一個峰值，說明這兩立柱后方的的漩渦為單頻脫落. 此外，4個立柱的頻譜圖中均有一峰值對應的頻率為f = 0.018 5 Hz. 綜上可知，平臺立柱后方漩渦脫落的主頻率為f= 0.018 5 Hz，部分立柱還伴有高頻或低頻成分.

由圖10可知，30°來流角下，中柱、邊柱2和邊柱3的升力系數(shù)頻譜雖然有多個峰值，但對應的頻率極低. 受上游尾流的干擾，下游邊柱1的升力系數(shù)頻譜沒有明顯的峰值. 可見在此來流角下，平臺立柱后方漩渦脫落不規(guī)則.

由圖11可知，60°來流角下，中柱頻譜沒有峰值，這是因為流體繞過中柱后，撞擊在下游的底座1和邊柱1上，沒有形成周期性的渦脫. 其余立柱的頻譜均有多個峰值，但沒有相同頻率的成分. 由此可見，平臺立柱后方的漩渦呈多頻脫落模式. 其他流速下，各柱的升力系數(shù)頻譜也具有類似的性質(zhì)，限于篇幅，不再贅述.

2.2 ? 流場分析

2.2.1 ? 流場相干結(jié)構(gòu)與尾流干擾

湍流場中有組織性的漩渦結(jié)構(gòu)稱之為相干結(jié)構(gòu)，其在湍流的發(fā)展演化中起著重要的作用[20-21]. u = 1.06 m/s時，0°、30°來流角下的渦量等值面圖（按x方向速度u著色）分別如圖12（a）、（b），渦量等值面采用Hunt等[22]建議的Q判據(jù)（Q-criterion）表示. 從圖12中可以看出，流體流經(jīng)平臺后，形成了多種形式和尺度不一的漩渦. 圖12（a）和（b）中的弧狀渦包成月牙形分布在浮箱底座的底面和頂面邊緣，反映出流體流經(jīng)浮箱底座時，在底面和頂面邊緣發(fā)生流動分離，并在隨后的區(qū)域形成了馬蹄渦. 圖12（a）中，邊柱1、邊柱2和邊柱3后方有成對的、錯位排布的管狀渦包，說明這些立柱后方形成了交替脫落、錯位排布的尾渦，尾渦流經(jīng)中柱后，逐漸拉長而發(fā)展成為流向渦. 從圖12（b）中可以看出，邊柱3后方形成了斜向下的錐形渦包，說明該處漩渦結(jié)構(gòu)的半徑隨著水深的增加而逐漸縮小. 浮箱底座3后方呈螺旋狀上升的螺旋渦包，反映出其后方形成了發(fā)夾渦[21]，該漩渦結(jié)構(gòu)由渦頭和兩條渦腿組成，渦頭受到渦腿旋轉(zhuǎn)向上的抬升力而上揚. 靠近浮箱底座的發(fā)夾渦的尺度明顯小于遠離底座的發(fā)夾渦的尺度，說明隨著湍流的發(fā)展，發(fā)夾渦逐漸發(fā)展壯大.

流體流經(jīng)平臺時，平臺的立柱、連桿和浮箱底座的尾流之間會相互干擾. 如圖12（a）中所示，邊柱2后方尺度較大的漩渦結(jié)構(gòu)與下游斜桿2后方尺度較小的漩渦結(jié)構(gòu)相互融合，最終匯合成為一個整體，相當于大漩渦“吞并”了小漩渦. 如圖12（b）中所示，浮箱底座3后方形成的發(fā)夾渦在上揚的過程中，由于渦腿的旋轉(zhuǎn)效應，誘導了邊柱3的尾渦向斜下方運動，最終與浮箱底座3的尾渦匯合在一起，并導致了次生流向渦的形成. 此外，由于該發(fā)夾渦渦腿的旋轉(zhuǎn)效應，也誘導了中柱后方的流向渦發(fā)生偏轉(zhuǎn)彎折，從而改變了中柱尾流的流動方向. 由以上分析可知，尾流干擾的機理正是相干結(jié)構(gòu)之間的相互作用.

2.2.2 ? 空間相關(guān)性分析

流場內(nèi)i、j兩點的脈動量Xi、Xj在空間上的相關(guān)性可用相關(guān)系數(shù)R來表示，其定義[23]表示為：

當i = j時，表示該點的脈動量自相關(guān)，當i ≠ j時，表示兩點的脈動量互相關(guān). 相關(guān)系數(shù)的大小表征了脈動量在流場不同位置的關(guān)聯(lián)程度，R = ±1表示完全相關(guān)，R = 0表示不相關(guān).

圖13為0°來流角下y = 0平面結(jié)構(gòu)的剖面圖. 為考察上下游立柱尾流間的空間相關(guān)性，選擇x方向速度u作為脈動量，計算了y = 0平面，z = 1.17D高度處邊柱2（bz2）尾流點與中柱（zz）尾流點的相關(guān)系數(shù)，計算結(jié)果如圖14所示，其中，圖14（b）為圖14（a）的水平投影圖.

由圖14可知，邊柱2尾流與中柱尾流的x方向速度在一定區(qū)間內(nèi)具有較強的相關(guān)性. 經(jīng)計算，x坐標在[0.85D，1.1D]區(qū)段的邊柱2尾流與x坐標在[3.1D，3.4D]區(qū)段的中柱尾流的x方向速度的相關(guān)系數(shù)為0.4 < R < 0.55，表現(xiàn)出較強的正相關(guān)性，而x坐標在[1.7D，2D]區(qū)段的邊柱2尾流與x坐標在[3.1D，3.4D]區(qū)段的中柱尾流的x方向速度的相關(guān)系數(shù)為-0.55 < R < -0.4，表現(xiàn)出較強的負相關(guān)性. 這說明雖然邊柱2后方交替脫落的尾渦流經(jīng)中柱后，變成了流向渦，但兩立柱的尾流之間仍有著較強的空間相關(guān)性.

3 ? 結(jié) ? 論

本文通過DDES方法研究了0°、30°和60°來流下半潛式浮式風機平臺繞流的力學特性，并對平臺的尾流進行了深入分析，得出以下結(jié)論：

1）由于尾流的干擾，下游立柱的阻力系數(shù)平均值隨流速變化而波動，離散性大，而上游立柱的阻力系數(shù)平均值變化不大. 在0°來流時，平臺立柱后方的漩渦脫落以某一主頻率為主，并伴有其他高頻或低頻成分;60°來流時，平臺立柱后方的漩渦脫落呈多頻模式.

2）流場相干結(jié)構(gòu)具有多尺度性，且類型多樣，包括發(fā)夾渦、流向渦、次生流向渦、馬蹄渦等. 相干結(jié)構(gòu)之間的相互作用是尾流干擾的內(nèi)在原因. 尾流干擾存在兩種模式：處于上下游位置的立柱與連桿之間的尾渦相互交融，導致它們的尾流成為一個整體;浮箱后的發(fā)夾渦誘導立柱漩渦改變運動方向，導致立柱尾流匯入浮箱尾流或者運動發(fā)生偏轉(zhuǎn).

3）流體流經(jīng)成串列排布的立柱時，即使流場的相干結(jié)構(gòu)發(fā)生了轉(zhuǎn)變，上下游立柱尾流之間仍有著較強的空間相關(guān)性.

參考文獻

[1] ? ?鄧露，黃民希，肖志穎，等. 考慮氣動阻尼的浮式風機頻域響應分析[J]. 湖南大學學報（自然科學版），2017，44（1）：1—8.

DENG L，HUANG M X，XIAO Z Y，et al. Analysis on frequency response of floating wind turbine considering the influence of aerodynamic damping[J]. Journal of Hunan University（Natural Sciences），2017，44（1）：1—8. （In Chinese）

[2] ? ?周強，廖海黎，曹曙陽. 高雷諾數(shù)下方柱繞流特性的數(shù)值模擬[J]. 西南交通大學學報，2018，53（3）：533—539.

ZHOU Q，LIAO H L，CAO S Y. Numerical study of flow characteristics around square cylinder at high reynolds number[J]. Journal of Southwest Jiaotong University，2018，53（3）：533—539. （In Chinese）

[3] ? ?DUAN F，HU Z Q，WANG J. Investigation of the VIMs of a spar-type FOWT using a model test method[J]. Journal of Renewable and Sustainable Energy，2016，8（6）：063301.

[4] ? ?KUMAR A，CHANG H S，GOWDA L，et al. Influence of comer radius on the near wake structure of a transversely oscillating square cylinder[J]. Journal of Mechanical Science and Technology，2009，23（9）：2390—2416.

[5] ? ?ZHANG W，SAMTANEY R. Low-Re flow past an isolated cylinder with rounded corners[J]. Computers and Fluids，2016，136：384—401.

[6] ? ?端木玉，萬德成. 雷諾數(shù)為3900時三維圓柱繞流的大渦模擬[J]. 海洋工程，2016，34（6）：11—20.

DUAN M Y，WAN ?D C. Large-eddy simulation of the flow past a cylinder with Re=3900[J]. The Ocean Engineering，2016，34（6）：11—20. （In Chinese）

[7] ? ?槐文信，柳夢陽，楊中華. 基于LES的圓柱尾流漩渦特性分析[J]. 華中科技大學學報（自然科學版），2018，46（5）：95—99.

HUAI W X，LIU M Y，YANG Z H.Analysis of vortex characteristics in wake of circular cylinder based on LES[J]. Journal of Huazhou University of Science and Technology （Natural Science Edition），2018，46（5）：95—99. （In Chinese）

[8] ? ?賈曉荷，劉樺. 雙圓柱繞流的大渦模擬[J]. 水動力學研究與進展A輯，2008（6）：625—632.

JIA X H，LIU Y. Large eddy simulation of flow around two circular cylinders[J]. Chinese Journal of Hydrodynamics，2008（6）：625—632. （In Chinese）

[9] ? ?鄭宇華，顧杰. 兩并列矩形柱繞流的PIV試驗研究[J]. 應用力學學報，2018，35（3）：465—471+682.

ZHENG Y H，GU J. PIV experimental study on flow around two parallelized rectangular columns[J]. Chinese Journal of Applied Mechanics，2018，35（3）：465—471+682. （In Chinese）

[10] KITAGAWA T，OHTAB H. Numerical investigation on flow around circular cylinders in tandem arrangement at a subcritical reynolds number[J]. Journal of Fluids Structure，2008，24（5）：680—699.

[11] ?李聰洲，張新曙，胡曉峰，等. 高雷諾數(shù)下多柱繞流特性研究[J]. 力學學報，2018，50（2）：233—243.

LI C Z，ZHANG X S，HU X F，et al. The study of flow past multiple cylinders at high reynolds number[J]. Chinese Journal of Theoretical and Applied Mechanics，2018，50（2）：233—243. （In Chinese）

[12] ?劉為民，谷家揚，盧燕祥. 陣列四柱繞流的數(shù)值模擬[J]. 江蘇科技大學學報（自然科學版），2012，26（6）：534—540.

LIU W M，GU J Y，LU Y X. Study on the characteristics of flow around a four square cylinder array[J]. Journal of Jiangsu University of Science and Technology（Natural Science Edition），2012，26（6）：534—540. （In Chinese）

[13] ?劉正浩，趙偉文，萬德成. 四柱式半潛平臺繞流的三維數(shù)值模擬[C]//. 中國力學大會-2017暨慶祝中國力學學會成立60周年大會論文集（A）. 北京：中國力學學會，2017：1—12.

LIU Z H，ZHAO W W，WAN D C. Numerical investigation of characteristics of flow past semi-submersible platform[C]//. Proceedings of the Chinese Congress of Theoretical and Applied Mechanics （CCTAM-2017） and the celebration of the 60th anniversary of Chinese Society of Theoretical and Applied Mechanics（A）. Beijing：Chinese Society of Theoretical and Applied Mechanics，2017：1—12. （In Chinese）

[14] ?SPALART P R，JOU W H，STRELETS M，et al. Comments on the Feasibility of LES for wings，and on a hybrid RANS/LES approach[C]// Proceedings of First AFOSR International Conference on DNS/LES Approach. Ruston：Greyden Press，1997：137—147.

[15] ?SPALART P R，DECK S，SHUR M L，et al. A new version of detached-eddy simulation，resistant to ambiguous grid densities[J]. Theoretical Computational Fluid Dynamics，2006，20（30）：181—195.

[16] ?杜若凡，閻超，韓政，等. DDES延遲函數(shù)在超聲速底部流動中的性能分析[J]. 北京航空航天大學學報，2017，43（8）：1585—1593.

DU R F，YAN C，HAN Z，et al. Performance of delayed functions in DDES for supersonic base flow[J]. Journal of Beijing University of Aeronautics and Astronautics，2017，43（8）：1585—1593. （In Chinese）

[17] ?ROBERTSON A，JONKMAN J，MASCIOLA M，et al. Definition of the semi-submersible floating system for phase II of OC4[R]. Golden：National Renewable Energy Laboratory，2014：1—38.

[18] ?張新曙，胡曉峰，尤云祥，等. 深海多立柱浮式平臺渦激運動特性研究[J]. 力學學報，2016，48（3）：593—598.

ZHANG X S，HU X F，YOU Y X，et al. Investigation on the characterisitics of vortex induced motion of a deep sea muti-column floating platform[J]. Chinese Journal of Theoretical and Applied Mechanics，2016，48（3）：593—598. （In Chinese）

[19] ?SHEIKH R，BROWN A. Extreme vertical deep water current profiles in the South China Sea，Offshore Borneo[C]// Proceedings of the 29th International Conference on Ocean，Offshore and Arctic Engineering. Shanghai：American Society of Mechanical Engineers，2010：585—595.

[20] ?PEREIRA F S，ECA L，VAZ G，et al. Challenges in Scale-Resolving Simulations of turbulent wake flows with coherent structures[J]. Journal of Computational Physics，2018，363：98—115.

[21] ?陳槐. 封閉槽道紊流相干結(jié)構(gòu)研究[D]. 北京：清華大學土木水利學院，2015：3—5.

CHEN H. Coherent structures in turbulent channel flows[D]. Beijing：School of Civil Engineering，Tsinghua University，2015：3—5. （In Chinese）

[22] ?HUNT J C R，WRAY A A，MOIN P. Eddies，streams and convergence zones in turbulent flows[C]// Proceedings of the 1988 Summer program：studying turbulence using numerical simulation databases，2. Palo Alto：Stanford University，1988：193—208.

[23] ?胡昊，李新凱，戴麗萍，等. 渦發(fā)生器控制平板邊界層分離的大渦模擬[J]. 航空動力學報，2016，31（2）：493—502.

HU H，LI X K，DAI L P，et al. Vortex generators control flat boundary layer separation of large eddy simulation[J]. Journal of Aerospace Power，2016，31（2）：493—502. （In Chinese）