圓柱繞流的數(shù)值模擬

1、圓柱繞流的數(shù)值模擬張玉靜20070360204過控（2）班化工與能源學院摘要:使用計算流體力學軟件FLUENT,模擬均勻來流繞固定圓柱的流動,模擬雷諾數(shù)為5，20,40,100時的繞流流動,得到流場的流函數(shù)等值線圖和速度矢量圖。計算結果表明：當雷諾數(shù)增加時，流動表現(xiàn)出一系列不同的構造。當Re=5時，流動不發(fā)生分離,其后未形成旋渦，當Re=20,40,100時，流體發(fā)生分離，其后形成旋渦，且旋渦隨著Re的增大而增大。利用計算流體力學軟件FLUENT可以成功地模擬圓柱繞流問題,反映出流動特性。關鍵詞：圓柱繞流；FLUENT；雷諾數(shù)Abstract:Uniformflowaroundamountin

2、gcylinderissimulatedwiththeapplicationofFLUENTsoftwarewhileReynoldsnumberis5,20,40,100.Streamfunctionandvelocityvectordistributionsareindicated.TheresultsshowthataseriesofconstructionappearsasReynoldsnumberincreases.WhenReis5,Flowseparationdoesnotoccur,anditdoesnotformvortex.WhenReis20,40,100,Flowse

3、parationoccurs,anditformsvortex.VortexincreaseswiththeincreaseofRe.UsingcomputationalfluiddynamicssoftwareFLUENTcansuccessfullysimulateflowaroundcylindrical,reflecttheflowcharacteristic.Keywords：Flowaroundacircularcylinder；FLUENT；Reynoldsnumber柱繞流理論分析研究的狀況一個世紀以來,圓柱繞流一直是眾多理論分析、實驗研究及數(shù)值模擬對象。但迄今對該流動現(xiàn)象物理

4、本質(zhì)的理解仍是不完整的。圓柱繞流中,起決定作用的是雷諾數(shù),但還受到許多因素,如阻塞比,來流湍流度,下游邊界條件等的影響。隨著雷諾數(shù)的增加,粘性不可壓縮流體繞圓柱的流動會呈現(xiàn)各種不同的流動狀態(tài),在小雷諾數(shù)時,流動是定常的,隨著雷諾數(shù)的增加,圓柱后會出現(xiàn)一對尾渦。當雷諾數(shù)較大時,尾流首先失穩(wěn),出現(xiàn)周期性的振蕩。而后附著渦交替脫落,瀉入尾流形成Karman渦街,隨著雷諾數(shù)的增加,流動變得越來越復雜,最后發(fā)展為湍流。White認為圓柱渦流具有經(jīng)典性的重要意義。一般認為圓柱繞流有2種定常的流動圖案:雷諾數(shù)為較小時,圓柱后無尾渦；當雷諾數(shù)為較大時,圓柱后有一對對稱的尾渦。關于定常流失穩(wěn)以及出現(xiàn)湍流的臨界雷

5、諾數(shù)主要是通過應用流場顯示技術觀察流動形態(tài)得到的,所以不是準確值。對于分界點雷諾數(shù)就有不同的見解,Kovasznay,Roshko等認為定常流動失穩(wěn)的臨界雷諾數(shù)大約為40。而從周期性尾流到湍流的詳細的轉變過程的實驗研究似乎還是空白。對均勻來流繞固定圓柱的二維平面流動,國內(nèi)外許多學者進行過大量的研究。決定圓柱繞流流態(tài)的是雷諾數(shù)(Re)的值,Re5時，流動不發(fā)生分離,5Re40,在圓柱體后面出現(xiàn)一對位置固定的旋渦；40Re150,旋渦擴大，然后有一個旋渦開始脫落,接著另一個也脫落,在圓柱體后面又生成新的旋渦,這樣逐漸發(fā)展成兩排周期性擺動和交錯的旋渦,即Karman渦街。Re150,渦街是層流，15

6、0Re300,旋渦由層流向湍流轉變。300Re3X105，稱為亞臨界區(qū)。此時,柱體表面上的邊界層為層流,而柱體后面的渦街已完全轉變?yōu)橥牧?并按一定的頻率發(fā)放旋渦。3X105Re3X106，稱為過渡區(qū)。目前,對圓柱繞流的三維數(shù)值研究則不多,然而試驗表明,即使均勻來流垂直流過等截面的圓柱體，當Re足夠大時，也會呈現(xiàn)三維的流動狀態(tài)。圓柱繞流屬于非定常分離流動問題,在工業(yè)工程中的應用非常廣泛,數(shù)值模擬是研究這類問題的有力工具。但是,控制方程的高度非線性以及邊界條件的多樣性嚴重阻礙了這類問題的解決。到目前為止，有限差分、有限元、大渦模擬(LES)、直接數(shù)值模擬(DNS)等方法都先后被應用。但是,這些方法

7、都需要增加網(wǎng)格生成的附加工作,尤其是在遇到復雜幾何區(qū)域時更需要巨大的網(wǎng)格數(shù)量,大大影響計算效率。本文利用N-S方程,對固定圓柱繞流進行了三維數(shù)值模擬,利用計算流體力學軟件FLUENT,選取Re=20,40,100,模擬了圓柱周圍的流場，得到流場的渦量等值線和速度矢量圖。2數(shù)學模型與數(shù)值方法對不可壓縮粘性流體,在直角坐標系下,其運動規(guī)律可用N-S方程來描述,連續(xù)性方程和動量方程分別為：(1) 連續(xù)性方程dudvdwn+=0dxdydz式中，u、v、w表示速度在x、y、z方向上的分量。(2) 動量方程1)4)0(uu)1dpd2uij=f一+Vi_dxipdxdxxjiij式中，xi(i=l,2,

8、3)為坐標系坐標，Xj(j=l,2,3)為坐標系坐標，ui(i=1,2,3)為沿i方向的速度分量，為沿i方向的質(zhì)量力，P是靜壓，P為空氣密度。3模型建立與網(wǎng)格劃分一個無窮長直徑為2.0m的圓柱體，放置在無窮遠來流速度為l.Om/s不受干擾的均勻橫流中,計算區(qū)域如圖1所示。網(wǎng)格劃分采用四邊形網(wǎng)格，網(wǎng)格數(shù)為16000。Gambit中網(wǎng)格劃分結果如圖2所示。圖1計算區(qū)域圖2網(wǎng)格劃分4邊界條件與求解設置進口為Velocity-inlet，出口為Outflow，其他邊界默認為Wall，計算區(qū)圖4指定流體區(qū)域域設為fluid,在gambit中設置界面如圖3所示。圖3邊界條件置界面求解設置如下：1）控制方程

9、采用連續(xù)方程和動量方程，而不考慮能量方程。（2）求解器選擇COUPLEDIMPLICITSOLVER，選擇Steady狀態(tài)。3）為了提高計算精度，差分格式采用二階迎風格式。（4）速度場和壓力場的求解基于Simple算法，松弛因子默認。5）模型選擇層流流動。6）殘差設為e-6。在fluent中的操作具體如下：（1）讀入從gambit中導出的msh文件，檢查網(wǎng)格質(zhì)量，并確定模型單位如圖5所示。2）設置求解器，如圖6所示。3）選擇層流模型，如圖7所示。圖5確定單位圖6選擇求解器圖7選擇層流模型4）改變流體介質(zhì)，如圖8所示。5）設置入口速度，如圖9所示。6）設置殘差，如圖10所示。7）初如化流場，進行

10、迭代求解。圖8改變流體介質(zhì)ZoneName|uelocityinlet.16MomentumThermalRadiationSpeciesDPMMultiphaseUDSVelocitySpecificationMethodReferenceFrameMagnitudeNormaltoBoundaryAbsoluteVelocityMagnitude(m/s)孑constantOK|Cancel|Help|圖9設置入口速度PlottingIterationsNormalizePScaleonitorsStorageIterations1000NormalizationWindowConverg

11、enceCriterion圖10設置殘差5數(shù)值計算與結果討論雷諾數(shù)由圓柱體直徑和自由來流速度確定。通過改變粘度,得到不同的雷諾值(5,20,40,100),Re=pUD/p。邊界條件:進口和側面邊界始終為自由來流條件U=1.0m/s,流動出口為零法向梯度出口邊界。流體特性假定為常數(shù),結果見表1。表1流體特性Re52040100密度p/(kg/m3)1111動力粘度p/(Pa.s)0.40.10.050.02Re=20時的殘差曲線如圖11所示。(a)Re=5時的速度分布b)Re=20時的速度分布(c)Re=40時的速度分布(d)Re=100時的速度分布圖12不同Re下的速度分布a)Re=5時的流

12、線分布(b)Re=20時的速度分布(c)Re=40時的流線分布(d)Re=100時的流線分布圖13不同Re下的流線分布圖12為不同Re下的速度分布，圖13為不同Re下的流線分布，由圖12、圖13可以看出流體流過圓柱體時，受圓柱體影響，圓柱后面會形成一低速區(qū)。當Re=5時，流動不發(fā)生分離,其后未形成旋渦，當Re=20,40,100時，流體發(fā)生分離，其后形成旋渦，且旋渦隨著Re的增大而增大。6結論本文使用計算流體力學軟件FLUENT,模擬均勻來流繞固定圓柱的流動，模擬雷諾數(shù)為5，20,40,100時的繞流流動,得到流場的速度分布云圖和流線圖。計算結果表明:(1)當雷諾數(shù)增加時，流動表現(xiàn)出一系列不同的構造。當Re=5時，流動不發(fā)生分離,其后未形成旋渦，當Re=20,40,100時，流體發(fā)生分離，其后形成旋渦，且旋渦隨著Re的增大而增大。(2)利用計算流體力學軟件FLUENT可以成功地模擬圓柱繞流問題,反映出流動特性。參考文獻:酆慶增.圓柱繞流的非線性動力學J.力學進展,1994,24(4):525544.2蘇銘德,康欽軍.亞臨界雷諾數(shù)下圓柱繞流的大渦模擬J.力學學報,1999,31(1):100105.葉春明，吳文權.數(shù)