端木玉，萬德成

( 1. 上海交通大學(xué) 船舶海洋與建筑工程學(xué)院 海洋工程國家重點(diǎn)實(shí)驗(yàn)室 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240; 2. 江蘇海事職業(yè)技術(shù)學(xué)院，江蘇 南京 212003)

雷諾數(shù)為3 900時(shí)三維圓柱繞流的大渦模擬

端木玉1,2，萬德成1

( 1. 上海交通大學(xué) 船舶海洋與建筑工程學(xué)院 海洋工程國家重點(diǎn)實(shí)驗(yàn)室 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240; 2. 江蘇海事職業(yè)技術(shù)學(xué)院，江蘇 南京 212003)

自從人們對層流的圓柱繞流現(xiàn)象有了系列研究及清楚的認(rèn)識后，人們逐漸把目光投向湍流的圓柱繞流，但相關(guān)研究主要關(guān)注于模擬湍流方法的數(shù)值格式和精度問題，而忽略了對高雷諾下圓柱繞流流場本身的認(rèn)識及規(guī)律的總結(jié)?；陂_源代碼OpenFOAM的大渦模擬方法以連續(xù)方程和Navier-Stokes方程為控制方程，選用Smargorinsky模式為亞格子應(yīng)力模型，采用有限體積法和一次預(yù)測兩次修正的PISO算法，對Re=3 900時(shí)三維圓柱繞流問題進(jìn)行了數(shù)值模擬研究，并著重分析了其尾流特征和性質(zhì)。數(shù)值計(jì)算結(jié)果表明：大渦模擬方法可以模擬出細(xì)致的流場結(jié)構(gòu)，該雷諾數(shù)下的圓柱繞流具有很強(qiáng)的三維及湍流效應(yīng)，在圓柱后方約一倍直徑的范圍內(nèi)存在回流區(qū)域，在靠近圓柱壁面的尾流區(qū)域的速度剖面呈“U”型，遠(yuǎn)離壁面的速度剖面呈“V”型。瞬時(shí)速度剖面始終圍繞著時(shí)均速度的周圍脈動(dòng)，且距離圓柱越遠(yuǎn)瞬時(shí)速度場的脈動(dòng)范圍越大。

大渦模擬;圓柱繞流;OpenFOAM;尾流分析

Abstract: After a series of studies of laminar flow around circular cylinders, clear understanding was gained. People then set their sights on turbulent flow problems. However, relevant research mainly focused on the scheme and precision of numerical simulation methods at high Reynolds numbers; less attention was paid to the understanding of turbulent flow around circular cylinders itself. The simulations were carried out using large eddy simulation (LES) solver within the open source code OpenFOAM. The Navier-Stokes equations were discretized by the Finite Volume Method; the Smargorinsky model was chosen as Sub-grid stress model; the PISO algorithm with two prediction steps and one correction step were used to deal with the coupling between velocity and pressure. This paper emphatically analyzes the wake features of the three-dimensional flow around a circular cylinder atRe=3 900. The numerical results show that LES can simulate the detailed flow structures, which exhibit strong three dimensional effects and turbulence effects at the current Reynolds number. The main conclusions of this study are that there exists a recirculation region about one diameter after the cylinder. The mean velocity profile in the region close to the cylindrical wall shows “U” shape, while the velocity profile far from the wall exhibits “V” shape. The instantaneous velocity, which fluctuates around the mean velocity, has a larger amplitude of fluctuation as the distance from the cylinder gets farther.

Keywords: large-eddy simulation; flow past a cylinder; OpenFOAM; wake analysis

圓柱繞流問題一直是流體力學(xué)領(lǐng)域非常經(jīng)典的算例，它能很好地揭示分離流、渦流及脫落過程等復(fù)雜流動(dòng)現(xiàn)象。同時(shí)圓柱繞流問題也有很高的工程應(yīng)用價(jià)值，航行中的飛機(jī)、水中行駛的船舶、海上的石油平臺以及橋墩等，都存在繞流。

在影響圓柱繞流問題的若干影響參數(shù)中，雷諾數(shù)起著決定性作用。因?yàn)橥牧骶哂胁灰?guī)則性、多尺度性及復(fù)雜的非線性，所以一直以來有效地描述湍流的性質(zhì)是流體力學(xué)中的一個(gè)重大難題。自從人們對低雷諾數(shù)下圓柱繞流現(xiàn)象有了系列研究及清楚的認(rèn)識后，便逐漸把目光投向高雷諾數(shù)下的圓柱繞流，如何可以精確描述高雷諾數(shù)下圓柱繞流的湍流尾流特征是今后研究工作的重點(diǎn)。湍流狀態(tài)下的圓柱繞流的尾流表現(xiàn)為高度復(fù)雜的三維非穩(wěn)態(tài)、帶旋轉(zhuǎn)的不規(guī)則的隨機(jī)流動(dòng)，圓柱繞流的尾流是由不同尺度漩渦疊合而成的，大尺度的漩渦破裂后形成小尺度的漩渦，小尺度的漩渦也不斷潰滅，同時(shí)新的漩渦也源源不斷產(chǎn)生，這使得研究高雷諾數(shù)下的圓柱繞流變得非常困難。

對湍流的模擬常用直接模擬(DNS)、大渦模擬(LES)和雷諾時(shí)均(RANS)三種方法。DNS是最精確的方法，但對計(jì)算機(jī)要求很高，網(wǎng)格必須等于或小于最小漩渦尺寸；LES可選用較大的網(wǎng)格，對超出網(wǎng)格尺寸以上的大尺度湍流進(jìn)行直接數(shù)值模擬，而小于網(wǎng)格尺寸的湍流采用網(wǎng)格內(nèi)模型來模擬，LES在復(fù)雜流動(dòng)中可以得到RANS無法獲得的湍流運(yùn)動(dòng)的細(xì)微結(jié)構(gòu)和流動(dòng)現(xiàn)象；RANS抹平了湍流中的微小細(xì)節(jié)，對復(fù)雜精細(xì)的湍流結(jié)構(gòu)的模擬會遇到困難。綜上所述，LES似乎是目前再現(xiàn)圓柱尾流復(fù)雜結(jié)構(gòu)的最佳方法。Breuer[1]采用大渦模擬方法對孤立圓柱繞流問題(Re=140 000)進(jìn)行了三維的數(shù)值模擬，研究了大渦模擬方法對于處理高雷諾數(shù)下流動(dòng)問題的適應(yīng)性并分析了不同亞格子應(yīng)力模式及網(wǎng)格質(zhì)量對結(jié)果的影響。Xu等[2]用大渦模擬方法對變截面半徑圓柱的可壓縮流場進(jìn)行了數(shù)值模擬，分析了波狀圓柱繞流與柱狀圓柱繞流問題的區(qū)別。Lourenco等[3]、Ong等[4]分別對流動(dòng)條件為Re=3 900的圓柱繞流做過模型試驗(yàn)，并提供了尾渦x/D≤3和3≤x/D≤10處的時(shí)均速度場的試驗(yàn)結(jié)果，實(shí)驗(yàn)只能給出一些時(shí)均后的流場參數(shù)信息。Breuer[5]用大渦模擬方法計(jì)算了Re=3 900下的圓柱繞流，但其目的并不在于研究流場物理現(xiàn)象的細(xì)節(jié)，而是為分析影響LES精度的數(shù)值及模型方面。Kravchenko等[6]也研究了Re=3 900時(shí)的圓柱繞流，論證了B樣條法是相比于迎風(fēng)格式和中心差分格式的一種精度更高的數(shù)值格式。趙偉文等[7]用大渦模擬方法模擬了高雷諾數(shù)三維圓柱Spar的渦激運(yùn)動(dòng)問題。王吉飛等[8-9]用大渦模擬方法模擬了高雷諾數(shù)三維方腔驅(qū)動(dòng)流問題。

從目前的研究進(jìn)展來看，雖然已有不少學(xué)者對高雷諾數(shù)下的圓柱繞流問題做過一些試驗(yàn)及數(shù)值研究，但試驗(yàn)的研究很難給出我們所關(guān)注的瞬時(shí)渦量場、速度場、壓力場等信息。而數(shù)值研究方面，大多是側(cè)重于分析不同湍流模型或數(shù)值格式對計(jì)算結(jié)果的影響。而針對高雷諾數(shù)下圓柱繞流問題這一物理現(xiàn)象展開深入研究的并不多見。鑒于前人對Re=3 900下在試驗(yàn)和數(shù)值方面有過諸多貢獻(xiàn)，因此本文選取Re=3 900的三維圓柱繞流進(jìn)行數(shù)值研究，不同的是本文側(cè)重于分析該雷諾數(shù)下的流場細(xì)節(jié)，以便了解湍流圓柱繞流問題的尾流特征，并將數(shù)值結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對比，驗(yàn)證并分析數(shù)值計(jì)算結(jié)果的有效性。文中給出了瞬時(shí)的渦量場、速度場，分析了流場中特征位置處的流場信息，更進(jìn)一步分析了瞬時(shí)場與速度場之間的關(guān)系，使得能深入地理解高雷諾數(shù)下圓柱繞流這一物理現(xiàn)象。

數(shù)值計(jì)算過程采用基于開源程序代碼庫OpenFOAM平臺上的大渦模擬方法，采用有限體積法和一次預(yù)測兩次修正的PISO算法，以三維不可壓縮N-S方程為控制方程。文中選用的OpenFOAM平臺，其全稱是Open Field Operation and Manipulation，采用C++語言編寫的面向?qū)ο蟮挠?jì)算流體力學(xué)開源代碼程序庫。其支持的流動(dòng)模型有層流、基于雷諾時(shí)均的湍流模型、LES及DNS等；支持多種軟件生成的復(fù)雜網(wǎng)格及動(dòng)態(tài)網(wǎng)格技術(shù)；支持SIMPLE、PISO及PIMPLE求解器。

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

1.1 大渦模擬的控制方程

大渦模擬的控制方程是將不可壓縮的N-S方程做過濾，得到：

為亞格子尺度雷諾應(yīng)力(SGS Renolds stress)，則式(1)可寫作

建立亞格子應(yīng)力模式使方程(4)封閉，選用最基本的Smagorinsky模式，其表達(dá)形式：

1.2 初始條件和邊界條件

本文定義的流場的初始條件和邊界條件：

1) 入口，速度入口；

2) 出口，出流邊界；

3) 上下面及前后面，對稱邊界；

4) 圓柱表面，無滑移固壁邊界。

2 數(shù)值計(jì)算與結(jié)果討論

圖1 圓柱繞流的幾何計(jì)算域Fig. 1 Geometry details of flow around a cylinder

采用的計(jì)算模型如圖1所示，計(jì)算域以圓柱底部的中心為坐標(biāo)原點(diǎn)，圓柱直徑D=0.01 m。計(jì)算域長度(x軸的方向)為15D，上游斷面距離圓柱中心為5D，下游斷面距離圓柱中心為15D，方向與來流方向一致；計(jì)算域?qū)挾?z軸方向)為πD，方向沿著圓柱方向；計(jì)算域高度(y軸方向)為10D，上下面距離圓柱為5D，方向平行于圓柱橫截面。入口處為定常來流，流速U0=0.394 m/s，ν=1.01×10-6m2/s，雷諾數(shù)Re=U0D/ν≈3 900。首先從網(wǎng)格收斂性方面進(jìn)行考察，這里選用2套網(wǎng)格，第一套網(wǎng)格總數(shù)為792 320，計(jì)算域xoy平面網(wǎng)格劃分如圖2所示，z向均勻劃分32層網(wǎng)格，圓柱周圍的網(wǎng)格如圖3所示。將第一套網(wǎng)格各個(gè)區(qū)域進(jìn)行加密后得到第二套網(wǎng)格，第二套網(wǎng)格總數(shù)為2 827 500，計(jì)算域xoy平面網(wǎng)格劃分如圖4所示，z向均勻劃分50層網(wǎng)格，圓柱周圍的網(wǎng)格如圖5所示。

圖2 第一套網(wǎng)格計(jì)算域xoy平面網(wǎng)格的劃分Fig. 2 The first case mesh generation in the xoy plane

圖3 第一套網(wǎng)格圓柱周圍的網(wǎng)格劃分Fig. 3 Mesh in the vicinity of cylinder for the first case

圖4 第二套網(wǎng)格計(jì)算域xoy平面網(wǎng)格的劃分Fig. 4 The second case mesh generation in the xoy plane

圖5 第二套網(wǎng)格圓柱周圍的網(wǎng)格劃分Fig. 5 Mesh in the vicinity of cylinder for the second case

表1 兩套網(wǎng)格的數(shù)值計(jì)算結(jié)果Tab. 1 Flow parameters computed by the meshes of the two cases

圖6 阻力系數(shù)、升力系數(shù)的時(shí)間歷程曲線(Re=3 900) Fig. 6 Time history of the lift coefficient Cl and the drag coefficient Cd

圖6給出了Re=3 900時(shí)升力和阻力系數(shù)的時(shí)間歷程曲線，可以看出升力系數(shù)雖圍繞0值上下振蕩，但幅值是變化的，與層流情況下的周期性振蕩且幅值不變有著顯著區(qū)別。隨著時(shí)間的推移，阻力系數(shù)沒有像層流情況下一樣趨于穩(wěn)定值，而是存在脈動(dòng)變化，以上都是湍流的顯著特征。

圖7、8、9分別是圓柱繞流尾流區(qū)x、y、z方向上的瞬時(shí)速度等值線圖。并與賈曉荷[10]的數(shù)值模擬結(jié)果進(jìn)行對比，該圖中虛線表示u<0，實(shí)線表示u>0。由于這里的上下圖并非描述的是同一時(shí)刻，且速度等值線的劃分也不相同，但從上下圖的對比可見計(jì)算結(jié)果的流場特征能夠很好地吻合。

在這三幅圖中既能看到小尺度的湍流也能看到大尺度的湍流，隨著流向距離的增大，流動(dòng)的結(jié)構(gòu)在尺度上也趨于增大。圖7在靠近圓柱的區(qū)域存在明顯的速度為負(fù)值的區(qū)域，即回流區(qū)，回流區(qū)后x向的速度等值線均為正。圖8可以看出y方向的瞬時(shí)速度有規(guī)律的正負(fù)交替，且在靠近圓柱的瞬時(shí)速度正負(fù)交替明顯，到遠(yuǎn)處就逐漸削弱了，說明在圓柱的尾部形成了一系列的卡門渦街，這些漩渦在圓柱附近形成后脫落，漩渦在運(yùn)動(dòng)的過程中逐漸擴(kuò)散開。圖9中可以看出z向的速度等值線在圓柱后無規(guī)律的正負(fù)交錯(cuò)存在，且距離圓柱越遠(yuǎn)等值線的正負(fù)差異越小，說明圓柱后形成的渦隨著往后運(yùn)動(dòng)而逐漸擴(kuò)散開來。

圖7 Re=3 900圓柱繞流尾流區(qū)x方向的瞬時(shí)速度等值線圖(x-z平面，y=0)Fig. 7 Instantaneous streamwise velocity in the (x-z, y=0) plane of the presented simulation

圖8 Re=3 900圓柱繞流尾流區(qū)y方向的瞬時(shí)速度等值線(x-z平面，y=0)Fig. 8 Instantaneous cross-flow velocity in the (x-z, y=0) plane of the presented simulation

圖9 Re=3 900圓柱繞流尾流區(qū)z方向瞬時(shí)速度等值線(x-z平面，y=0)Fig. 9 Instantaneous spanwise velocity in the (x-z, y=0) plane of the presented simulation

圖10給出了xoy(z=πD/2)平面處圓柱繞流的渦量等值線圖；圖11給出了圓柱尾部的瞬時(shí)流場渦量，圖中的渦量等勢面用Hunt等人[15]建議的Q準(zhǔn)則(Q-criterion)表示。Q的定義如下：

圖10與圖8的敘述也是吻合的，可以更加直觀的看出漩渦在圓柱附近形成后脫落，并且漩渦在運(yùn)動(dòng)的過程中逐漸擴(kuò)散開。但是Re=3 900時(shí)形成的渦街與低雷諾數(shù)的情況存在明顯的差異，低雷諾數(shù)時(shí)通常觀察到物體尾流左右兩側(cè)產(chǎn)生成對的、交替排列的、旋轉(zhuǎn)方向相反的反對稱渦旋，漩渦的脫落過程具有非常的規(guī)律性，循環(huán)往復(fù)。Re=3 900的工況下，每個(gè)瞬時(shí)在圓柱尾部的一側(cè)會有不同大小的多個(gè)漩渦脫落，隨后向中心線移動(dòng)，移動(dòng)的過程中漩渦的尺度逐漸增大并直至脫落，雖然整個(gè)過程都是在圓柱兩側(cè)交替生成漩渦并逐漸脫落的過程，但每一次漩渦生成及脫落的過程都不是完全相同的，存在一定的隨機(jī)性。圖11展示了Re=3 900下圓柱繞流的尾流呈現(xiàn)明顯的湍流及三維特性，渦面上的深色及淺色分布代表了相反的旋渦旋轉(zhuǎn)方向。

圖10 xoy (z=πD/2)平面處渦量等值線圖Fig. 10 Contours of instantaneous vorticity magnitude in the (x-y,z=πD/2) cross-sectional plane

圖11 瞬時(shí)的三維渦量等值面圖(Q=100) Fig. 11 Isosurface of instantaneous vorticity magnitude (Q-criterion, Q=100)

表2不同長細(xì)比圓柱繞流的流場參數(shù)與前人工作的對比

Tab.2Theparametersoftheflowaroundthecylinderwithdifferent-aspectratioandcomparedtoexistingresultsintheliterature

算例長細(xì)比St數(shù)基礎(chǔ)吸力系數(shù)平均阻力系數(shù)回流區(qū)長度-umin/U￥CaseA112π0.21510.84991.14761.7540.243CaseA2π0.21540.92631.13561.4170.325CaseA32π0.21910.93761.11491.3600.355Lysenko[12]:LESwithSMAGπ0.190.81.180.90.26Lysenko[12]:LESwithTKEπ0.2090.910.971.670.27Norberg[13]150-0.875---Norberg[13]26-0.78---LourencoandShih[3]20.50.22-0.991.220.247KravchenkoandMoin[6]2π0.210.941.04-0.35

圖12分別給出了同一時(shí)刻，相同渦量等勢面的瞬時(shí)三維渦量圖(俯視圖)。從圖中可見，展向的長度越長，圓柱尾流處的漩渦相對越加紊亂，漩渦尺度也相對較小一些，但這種差別并不是十分明顯。另外一個(gè)特征是對于展向長度較長的圓柱，漩渦要早于展向長度短的圓柱脫落。

(b) 圓柱展向長度為πD

(c) 圓柱展向長度為πD圖12 不同圓柱展向長度下的瞬時(shí)三維渦量等值面圖(Q=100)Fig. 12 Isosurface of instantaneous vorticity magnitude in the wake of a cylinder with axial length (Q-criterion, Q=100)

為便于與實(shí)驗(yàn)結(jié)果進(jìn)行對比以顯示圓柱周圍近尾流區(qū)及遠(yuǎn)尾流區(qū)的流場特征，現(xiàn)將圓柱附近的特征位置描述如圖13所示。所取的截面是xoy(z=πD/2)處的平面，y/D=0是該平面上從圓柱中心沿x軸正方向到x=15D處的直線，x/D=0.58,1.06,1.54,2.02是該平面上與y軸平行的直線，直線的范圍是-3.0D～3.0D。

圖13 圓柱周圍的特征位置Fig. 13 The characteristic locations near the cylinder

圖14是圓柱中心線y/D=0上x方向的平均流速分布。將數(shù)值計(jì)算結(jié)果與Lourenco等[3]及Ong等[4]的實(shí)驗(yàn)結(jié)果及Kravchenko等[6]的數(shù)值結(jié)果進(jìn)行對比，模擬出的趨勢和實(shí)驗(yàn)結(jié)果完全吻合，數(shù)值模擬的精度也較Kravchenko等[6]更接近于實(shí)驗(yàn)值。從圖中發(fā)現(xiàn)在圓柱后面約1個(gè)直徑的范圍內(nèi)平均速度為負(fù)值，該范圍是回流區(qū)域，這與前面圖7的描述是一致的?；亓鲄^(qū)的范圍與實(shí)驗(yàn)值很接近，但是平均速度值較實(shí)驗(yàn)值偏大一些，這可能是由于實(shí)驗(yàn)存在外界擾動(dòng)，從而導(dǎo)致近圓柱處尾流速度的提前過渡。

圖14 圓柱中心線y/D=0上x方向的時(shí)均流速分布Fig. 14 Streamwise mean velocity on the center line in the wake of a cylinder

圖15分別給出了x/D=0.58,1.06,1.54,2.02處x方向一系列瞬時(shí)的速度剖面及時(shí)均速度。這里共列出了37個(gè)瞬時(shí)的x速度剖面，時(shí)間間隔ΔT=0.02s，用淺色的線表示，深色的線是這37個(gè)速度剖面的平均值。從圖中可見瞬時(shí)速度剖面始終圍繞著時(shí)均速度的周圍脈動(dòng)，且距離圓柱越遠(yuǎn)脈動(dòng)的范圍越大，因?yàn)榫嚯x圓柱越遠(yuǎn)的地方物體對流體的控制越弱。時(shí)均速度剖面在所取的四個(gè)特征位置處的形狀也有其特征，平均速度的谷值一定在y/D=0的中心線上，但是各瞬時(shí)的谷值圍繞中心線振蕩。x/D=0.58處的時(shí)均速度呈“U”型，其余呈“V”型，隨著距離的增大，“V”的形狀也趨于平緩。在x/D=0.58處及x/D=1.54處的時(shí)均速度的谷值都接近于0，在x/D=1.06處的時(shí)均速度谷值在負(fù)值最大處，而在x/D=2.02處的時(shí)均速度谷值為正，以上這些特征都和圖14相符。

(a) x/D=0.58

(b) x/D=1.06

(c) x/D=1.54

(d) x/D=2.02

圖16 圓柱尾流在4個(gè)特征位置處x方向的時(shí)均速度分布 Fig. 16 Profiles of streamline mean velocity at four locations downstream of a cylinder

將上述四個(gè)特征位置處x方向的時(shí)均速度值與實(shí)驗(yàn)值進(jìn)行對比，如圖16所示。數(shù)值計(jì)算結(jié)果的趨勢與Lourenco & Shih[3]的實(shí)驗(yàn)結(jié)果完全吻合，距離圓柱越近的位置與實(shí)驗(yàn)結(jié)果越相近，距離圓柱越遠(yuǎn)的地方則有較大誤差。

圖17分別給出了x/D=1.06,1.54,2.02處y方向一系列瞬時(shí)的速度剖面及時(shí)均速度。同上一樣共列出了37個(gè)瞬時(shí)的x速度剖面，時(shí)間間隔ΔT=0.02s，用淺色的線表示，深色的線是這37個(gè)速度剖面的平均值。

圖17 不同x/D條件下y方向的瞬時(shí)速度和時(shí)均速度Fig. 17 Profiles of cross-flow velocity at different x/D position

將上述三個(gè)特征位置處y方向的時(shí)均速度值與實(shí)驗(yàn)值進(jìn)行對比，如圖18所示。從圖中可見所有y向平均速度剖面的形狀均關(guān)于y/D=0反對稱，與實(shí)驗(yàn)的趨勢基本吻合，數(shù)值上存在一定誤差。

圖18 圓柱尾流在3個(gè)特征位置處y方向的時(shí)均速度分布Fig. 18 Profiles of cross-flow mean velocity at three locations downstream of a cylinder

3 結(jié) 語

基于OpenFOAM的大渦模擬方法研究了Re=3 900時(shí)圓柱繞流問題，并對數(shù)值模擬結(jié)果進(jìn)行了分析，得到以下結(jié)論：

1) 在雷諾數(shù)較高的情況下必須采用三維的計(jì)算模型，二維簡化模型會對計(jì)算結(jié)果產(chǎn)生較大誤差，且在較高雷諾數(shù)下，圓柱繞流的漩渦脫落已經(jīng)具有了三維效應(yīng)；

2) 從升力系數(shù)和阻力系數(shù)的時(shí)歷曲線中可以明顯看出存在明顯的脈動(dòng)，無特征規(guī)律可循，說明在Re=3 900時(shí)已經(jīng)存在明顯的湍流效應(yīng)，另漩渦脫落的過程也說明了湍流效應(yīng)；

3)Re=3 900的工況下，每個(gè)瞬時(shí)在圓柱尾部的一側(cè)會有尺度不同的多個(gè)漩渦脫落，隨后向中心線移動(dòng)，移動(dòng)的過程中漩渦的尺度逐漸增大并直至脫落，且每一次漩渦生成及脫落的過程都不是完全相同的，存在隨機(jī)性；

4) 數(shù)值模擬出的St=0.222，換算得漩渦脫落的頻率fst=8.75 Hz；

5) 在靠近圓柱尾部約1D的范圍內(nèi)是回流區(qū)，之后沿x方向速度逐漸增大，并趨于平衡；

6) 對圓柱后x/D=0.58,1.06,1.54,2.02的四個(gè)特征位置處的流場進(jìn)行了分析，所有x向時(shí)均速度剖面的形狀均關(guān)于y/D=0對稱，x/D=0.58處的時(shí)均速度呈“U”型，其余呈“V”型，且隨著距離的增大，“V”的形狀也趨于平緩。瞬時(shí)速度剖面始終圍繞著時(shí)均速度的周圍脈動(dòng)，且距離圓柱越遠(yuǎn)瞬時(shí)速度場的脈動(dòng)范圍越大，這是因?yàn)榫嚯x圓柱越遠(yuǎn)的地方物體對流體的控制越弱。

[1] BREUER M. A challenging test case for large eddy simulation: high Reynolds number circular cylinder flow [J]. Heat and Fluids Flow, 2000, 21: 648-654.

[2] XU Changyue, CHEN Liwei, LU Xiyun. Large-eddy simulation of the compressible flow past a wavy cylinder[J]. Journal of Fluid Mechanics, 2010, 665: 238-273.

[3] LOURENCO L M, SHIH C. Characteristics of the plan turbulent near wake of a circular cylinder: A partical image velocimetry study[R]. Private Communication, 1993.

[4] ONG L, WALLACE J. The velocity field of the turbulent very near wake of a circular cylinder[J]. Experiments in Fluids, 1996, 20: 441-453.

[5] BREUER M. Large eddy simulation of the subcritical flow past a circular cylinder: numerical and modeling aspects[J]. International Journal For Numerical Methods in Fluids, 1998, 28: 1281-1302.

[6] KRAVCHENKO A G, MOIN P. Numerical studies of flow over a circular cylinder atReD=3 900[J]. Physics of Fluids, 2000, 12(2):403-417.

[7] 趙偉文，萬德成. 用大渦模擬方法數(shù)值模擬Spar平臺渦激運(yùn)動(dòng)問題[J]. 水動(dòng)力學(xué)研究與進(jìn)展，2015, 30(1): 40-46. (ZHAO Weiwen, WAN Decheng. Numerical investigation of vortex-induced motions of Spar platform based on large eddy simulation [J]. Chinese Journal of Hydrodynamics, 2015, 30(1): 40-46. (in Chinese))

[8] 王吉飛，萬德成.三維頂板斜向驅(qū)動(dòng)方腔流的有限元并行計(jì)算[J]. 海洋工程，2015, 33(2): 1-12. (WANG Jifei, WAN Decheng. Parallel simulation of 3D lid-driven cubic cavity flows at yaw by finite element method [J]. The Ocean Engineering, 2015, 33(2): 1-12. (in Chinese))

[9] WANG Jifei, WAN Decheng. Parallel finite element method for 3D lid-driven cubic cavity flows[J]. International Journal of Offshore and Polar Engineering, 2014, 24(2): 106-113.

[10] 賈曉荷,單圓柱及雙圓柱繞流的大渦模擬[D]. 上海：上海交通大學(xué)，2008. (JIA Xiaohe. Large eddy simulation of flow around one and two circular cylinders [D]. Shanghai: Shanghai Jiao Tong University, 2008. (in Chinese))

[11] WELLER H G, TABOR G, JASAK H, et al. A tensorial approach to computational continuum mechanics using object-oriented techniques [J]. Computers in Physics, 1998, 12(6): 620-631.

[12] LYSENKO D A, ERTESV?G I S, RIAN K E. Large-eddy simulation of the flow over a circular cylinder at reynolds number 3 900 using the OpenFOAM toolbox [J]. Flow, Turbulence and Combustion, 2012, 89(4): 491-518.

[13] NORBERG C. Effects of Reynolds number and low-intensity free stream turbulence on the flow around a circular[R]. Gothenburg: Chalmer University of Technolgy, 1987.

[14] WORNOM S, OUVRARD H, SALVTTI M V, et al. Variational multiscale largeeddy simulations of the flow past a circular cylinder: Reynolds number effects [J]. Comput. Fluids, 2011, 47(1), 44-50.

[15] HUNT J C R, WRAY A A, MOIN P. Eddies, streams and convergence zones in turbulent flows [R]. Center for Turbulence Research Report CTR-S88, 1988: 193-208.

Large-eddy simulation of the flow past a cylinder with Re=3 900

DUAN Muyu1, 2, WAN Decheng1

(1. State Key Laboratory of Ocean Engineering, School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai 200240, China; 2. Jiangsu Maritime Institute, Nanjing 212003, China)

O351.2

A

10.16483/j.issn.1005-9865.2016.06.002

1005-9865(2016)06-0011-10

2015-10-28

國家自然科學(xué)基金項(xiàng)目(51379125，51490675，11432009，51411130131)；長江學(xué)者獎(jiǎng)勵(lì)計(jì)劃(2014099)；上海高校特聘教授(東方學(xué)者)崗位跟蹤計(jì)劃(2013022)；工信部高技術(shù)船舶數(shù)值水池創(chuàng)新專項(xiàng)課題(2016-23/09)

端木玉(1981-),博士研究生，主要從事海洋工程水動(dòng)力學(xué)方面研究。

萬德成。E-mail：dcwan@sjtu.edu.cn