李 霖，張志國，王先洲，馮大奎

(華中科技大學(xué)船舶與海洋工程學(xué)院，湖北武漢430074)

0 引言

一個世紀(jì)以來，圓柱繞流問題一直是經(jīng)典的流體力學(xué)問題之一，也是眾多理論分析、數(shù)值模擬以及實驗研究的對象［1-2］。盡管研究歷經(jīng)時間如此之長，但這種最簡單的流動現(xiàn)象物理本質(zhì)的理解依然不完整。閔強利［3］計算了2種雷諾數(shù)Re=200，800情況下的圓柱繞流，發(fā)現(xiàn)在Re=200時產(chǎn)生的卡門渦街非常規(guī)則，而隨著雷諾數(shù)的增加，Re=800時，渦街出現(xiàn)不規(guī)則現(xiàn)象，并表現(xiàn)出三維效應(yīng)。文獻(xiàn) ［4-5］研究了不同截面形狀在低雷諾數(shù)下的繞流場，并重點研究了阻力、升力的周期性變化特性。

隨著計算機(jī)硬件的不斷發(fā)展，大渦模擬作為一種新的湍流模型相比原來的湍流模型在求解非穩(wěn)態(tài)流動方面具有巨大的潛力。大渦模擬方法的基本思想是用瞬時的N-S方程直接模擬計算湍流中的大尺度渦，而小尺度渦對大渦的影響則通過建立近似的模型來考慮，這種影響模型稱為亞格子尺度模型。王漢青［6］等人介紹了大渦模擬的理論進(jìn)展和發(fā)展趨勢，描述了當(dāng)前大渦模擬在工程中的具體應(yīng)用。指出大渦模擬在模擬計算從層流到湍流轉(zhuǎn)換、非定常湍流和高速湍流方面具有其他湍流模型無可比擬的優(yōu)勢。傅慧萍［7］以潛艇模型SUBOFF為研究對象，采用大渦模擬 (LES)方法求解了流動的非定常解。通過與試驗值以及采用RNGk-ε湍流模型得到的定常結(jié)果比較，驗證了大渦模擬方法的有效性。

本文采用大渦模擬方法求解低雷諾數(shù)下的圓柱繞流場，研究卡門渦街的周期性特性。對一個完整的渦街脫落周期進(jìn)行詳細(xì)討論，揭示了不同時刻旋渦的生成、發(fā)展、脫落和演化過程。通過研究圓柱表面的壓力變化說明產(chǎn)生周期性脫落渦對的原因，以及產(chǎn)生周期性變化的升、阻力的原因。

1 數(shù)值算法

自然界中的流動現(xiàn)象基本可以通過基于連續(xù)性假設(shè)的納維葉-斯托克斯方程 (N-S方程)來描述。在應(yīng)用于大渦模擬方法時，需要對N-S方程進(jìn)行一定過濾處理。

對于亞格子尺度模型［8］，采用標(biāo)準(zhǔn) Smagorinsky［9］模型的方法。

式中，vol是計算單元的體積。

亞格子尺度可表示為

式中，Smagorinsky 常量 Cs取為 0.255［10］。

2 計算模型

本文所采用的計算模型為直徑D=0.04 m的圓柱。來流速度為U=0.01 m/s，計算的雷諾數(shù)Re=400。計算域為方形，沿著來流方向長度為7.5D，尾流去流段長度為40D，兩側(cè)寬度為7.5D。流動在選定的計算域內(nèi)能夠充分發(fā)展，圓柱周圍采用邊界層網(wǎng)格，能夠捕捉圓柱周圍的細(xì)微流動。圖1為計算域示意圖，圖2為圓柱周圍的網(wǎng)格劃分。

3 計算結(jié)果與分析

3.1 計算結(jié)果與實驗結(jié)果對比

表1為計算所得的斯特勞哈爾數(shù)(St)以及阻力系數(shù)(Cd)與文獻(xiàn)中實驗結(jié)果的對比。

表1 實驗結(jié)果與計算結(jié)果Tab.1 Results from experiment and simulation

計算結(jié)果與實驗結(jié)果吻合良好，說明采用本文的數(shù)值計算方法是準(zhǔn)確而可靠的。

3.2 渦街脫落

圖3～圖8為圓柱繞流非穩(wěn)態(tài)流動中一個典型的渦街脫落周期中不同時刻圓柱周圍的壓力場分布。反映了一個典型渦街脫落周期內(nèi)漩渦的初生、形成、脫落及在尾流中不斷發(fā)展的整個過程。從圖3可看出，t=0時刻，上一周期脫落形成的渦在圓柱尾流中逐漸發(fā)展，圓柱的右下方約45°的位置形成一片明顯的負(fù)壓區(qū)，旋渦在負(fù)壓區(qū)內(nèi)逐漸生成。圖4中t=1/5T時刻，圓柱右下方的漩渦完全生成，沿著尾流方向移動，并有逐步脫落的趨勢。而圖5中t=2/5T時刻，圓柱右下方形成的漩渦完全脫落并沿著尾流方向發(fā)展。圖6中t=3/5T時刻，之前由圓柱右下方脫落的漩渦在尾流中逐漸發(fā)展，同時在圓柱的右上方約45°的位置形成一片負(fù)壓區(qū)，旋渦在負(fù)壓區(qū)內(nèi)逐漸生成。圖7中t=4/5T時刻圓柱右上方的漩渦完全生成，沿著尾流方向移動，并有逐漸脫落的趨勢。最后圓柱右上方形成的漩渦在t=T時刻完全脫落 (圖8)，脫落的漩渦在尾流中不斷發(fā)展，改變著圓柱周圍的壓力場分布。圓柱右上方區(qū)域和右下方區(qū)域形成的漩渦交替脫落，在尾流中不斷發(fā)展，形成了典型的卡門渦街。

3.3 水動力特性

圓柱在繞流場中受到的主要水動力為阻力和升力。通過無量綱化，可以用阻力系數(shù)Cd及升力系數(shù)CL來表征。

式中:D和L分別為圓柱所受的阻力和升力;U為來流速度;ρ為流體密度;L為特征長度，本文取為圓柱直徑。

圖9為計算求取的穩(wěn)定的非穩(wěn)態(tài)計算時間內(nèi)圓柱的阻力系數(shù)、升力系數(shù)時程曲線。計算時時間步長取為0.05 s。計算所得的結(jié)果是在非穩(wěn)態(tài)計算了較長時間，流動發(fā)展相對較為穩(wěn)定時提取出來的。

圖9 阻力系數(shù)和升力系數(shù)時程曲線Fig.9 Time-history curve of drag and lift coefficient

由圖9可知阻力和升力都呈現(xiàn)周期性的變化，這是由于周期性脫落的卡門渦街引起的。周期性脫落的漩渦造成了圓柱周圍壓力場的不斷變化。升力系數(shù)基本在0左右波動，而阻力系數(shù)基本在1.4左右波動，這說明阻力對圓柱的作用非常強，升力對圓柱的作用一般。文獻(xiàn)［11］指出，升力以斯特勞哈爾渦泄頻率fs變化，阻力以2倍渦泄頻率2fs變化。由圖9可以看出，阻力變化的頻率是升力變化頻率的2倍。進(jìn)一步對升力做功率譜密度分析，可以得到渦泄頻率fs，圖10為升力功率譜密度。圖10中的尖峰對應(yīng)為斯特勞哈爾渦泄頻率fs，其值在0.05 Hz附近，與文中之前計算的斯特勞哈爾數(shù) (St=fs*D/U)是對應(yīng)的。

圖10 升力功率譜密度Fig.10 The power spectral density of lift

為了深入研究圓柱上阻力和升力脈動的變化，可以研究圓柱瞬態(tài)的升力、阻力系數(shù)。非穩(wěn)態(tài)阻力系數(shù)可以定義為

非穩(wěn)態(tài)升力系數(shù)可以定義為

角度θ的定義見圖11;圖12為典型的渦街脫落周期中不同時刻圓柱周向的壓力系數(shù)變化。

圖12可看出，圓柱壓力系數(shù)在1個周期內(nèi)變化基本是對稱的。圓柱駐點處(θ=0°和θ=360°)壓力系數(shù)Cp有最大值，在1附近;隨著來流向圓柱兩側(cè)擴(kuò)展，壓力系數(shù)迅速減小。在θ=80°和θ=280°附近，流動產(chǎn)生分離，圓柱駐點后方的壓力在不同時刻呈現(xiàn)出較為規(guī)律的周期性變化。對比t=1/5T和t=4/5T時刻可以看出，在t=1/5T時刻圓柱的右下方區(qū)域 (180°～270°)的壓力大于右上方區(qū)域(90°～180°)，此時圓柱具有最大的負(fù)升力;在t=4/5T時刻圓柱的右下方區(qū)域 (90°～180°)的壓力大于右下方區(qū)域 (180°～270°)，此時圓柱具有最大的正升力。這樣升力就完成了1個渦街脫落周期中由最大正升力向最大負(fù)升力的轉(zhuǎn)變。上述過程不斷進(jìn)行，就形成了升力的周期性變化。

4 結(jié)語

采用大渦模擬的方法研究了雷諾數(shù)Re=400的圓柱繞流場，準(zhǔn)確捕捉了卡門渦街這一經(jīng)典流動現(xiàn)象，同時通過對水動力特性的研究，得出以下結(jié)論:

1)采用大渦模擬的方法進(jìn)行低雷諾數(shù)圓柱繞流場的研究是準(zhǔn)確可靠的，大渦模擬作為一種新的湍流模型是值得深入研究和普遍運用的;

2)圓柱的渦泄頻率是和圓柱的升力功率譜相互對應(yīng);

3)圓柱表面壓力的周期性變化是誘發(fā)漩渦周期性生成和脫落的原因。

本文僅用大渦模擬對二維圓柱進(jìn)行了繞流場的分析，但實際的流動都是具有三維效應(yīng)的，故今后將運用大渦模擬對圓柱的三維繞流場繼續(xù)研究。

［1］MAIR W A.Bluff bodies and vortex shedding［J］.Journal of Fluid Mechanics.，1971，45:209-224.

［2］MUJUMDAR A S.，DOUGLAS W J M.Vortex shedding from slender cylinders of various cross sections［J］.Journal of Fluids Engineer，1973，95(3):474-476.

［3］閔強利.低雷諾數(shù)卡門渦街?jǐn)?shù)值模擬［J］.四川兵工學(xué)報，2009，30(11):81-84.MIN Qiang-li.Simulation of Karman Vortex Street at low Reynolds number［J］.Journal of Sichuan Ordnance，2009，30(11):81-84.

［4］LEI Yue.Investigation of the flow and heat transfer around cylinders at low Re［J］.Advanced Materials Research，2011，156-157:1630-1634.

［5］BOGRAD D D，GARRISON D H.Moment and heat transfer from cylinders in laminar cross-flow at 10-4≤Re≤200［J］.International Journal of Heat and Mass Transfer，1998，41(22):3409-3430.

［6］王漢青，王志勇，寇廣孝.大渦模擬理論進(jìn)展及其在工程中的應(yīng)用［J］.流體機(jī)械，2004，32(7):23-27.WANG Han-qing，WANG Zhi-yong，KOU Guang-xiao.Theoretical development and application of Large Eddy Simulation in engineering［J］.Fluid Machinery，2004，32(7):23-27.

［7］傅慧萍，繆國平，大渦模擬在噪聲計算中的應(yīng)用［J］.上海交通大學(xué)學(xué)報，2009，43(8):1307-1312.FU Hui-ping，MIAO Guo-ping.Application of LES method to flow noise calculation［J］.Journal of Shang hai Jiaotong University，2009，43(8):1307-1312.

［8］FENG W，F(xiàn)ERNANDO，STOLL R.Evaluation of dynamic subgrid-scale models in large-eddy simulations of neutral turbulent flow over a two-dimensional sinusoidal hill［J］.Atmospheric Environment，2007，41(13):2719-2728.

［9］BOGEY C，BAILLY C.Large eddy simulations of round free jets using explicit filtering with/without dynamic Smagorinsky model［J］.International Journal of Heat and Fluid Flow，2006，27:603-610.

［10］MOIN P，KIM J.Numerical investigation of turbulent channel flow［J］.Journal of Fluid Mechanics，1982，118:341-377.

［11］SCHLICHTING H.邊界層理論［M］.北京:科學(xué)出版社，1979.SCHLICHTING H.Theory of boundary layer［M］.Beijing:Science Press，1979.