方腔頂蓋驅(qū)動流數(shù)值模擬

1、精選優(yōu)質(zhì)文檔-傾情為你奉上方腔頂蓋驅(qū)動流數(shù)值模擬張鑫（浙江理工大學 動力工程 2013G）摘 要：在計算流體力學的研究中，通常要計算方腔驅(qū)動流問題來檢驗各種N-S數(shù)值方法的有效性。本文利用Fluent軟件對標準計算流體力學測試算例方腔驅(qū)動流問題進行了模擬分析，其計算結果與文獻中的標準解符合的比較好。關鍵字：N-S方程 方腔驅(qū)動流 Fluent數(shù)值求解0引言流體流動的數(shù)值模擬廣泛應用于氣象、航天、機械、采礦等自然研究和工程計算的各個領域。近年來，隨著高性能計算與通信的迅速發(fā)展，針對流體流動的數(shù)值模擬以及求解相應Navier -Stokes方程(簡稱N-S方程)的高級算法研究現(xiàn)已成為目前國內(nèi)外備受

2、關注的熱點和前沿課題。Fluent軟件是用于模擬具有復雜外形的流體流動以及熱傳導的計算機程序，可以有效地模擬方腔驅(qū)動流問題，為計算流體力學的算法理論研究提供仿真參考。高殿榮等學者采用液壓沖擊進行了分析；韓善玲等分析流體在空腔內(nèi)的運動規(guī)律和物理機制，指出微小的凹凸是引起噪聲的原因之一。楊晶用Fluent軟件對方腔驅(qū)動流動進行了模擬分析，研究了不同雷諾數(shù)對計算結果的影響。1模型介紹下圖描述了本文所研究的物理模型，模型為邊長等于0.1m的正方形，上壁面為有一定速度的水，兩側(cè)壁面及地面均固定。流體材料為水，密度為998.2kg/m3，黏度。 a=0.4m2數(shù)值計算2.1、N-S方程本文控制方程采用納維

3、司托克斯方程，納維司托克斯方程是描述粘性不可壓縮流體動量守恒的運動方程。簡稱N-S方程。在直角坐標系中，可表達為如下所示：連續(xù)方程：動量方程：2.2、網(wǎng)格劃分及邊界條件設置在gambit軟件中建立模型劃分網(wǎng)絡，由于模型幾何形狀比較規(guī)則，故全部采用四邊形的的結構化網(wǎng)格，如下圖所示。邊界條件：壁面皆為壁面無滑移條件，其中上頂蓋以一定速度移動。網(wǎng)格總數(shù)為10000。2.3、fluent軟件求解計算導入.mesh文件后，在scale里面同一單位。如下圖所示： 求解器采用基于壓力，定常，層流模型。流體材料為液態(tài)水。邊界條件設置為wall1為平移界面，剩下的為無滑移的固定界面。對殘差進行設置后，初始化運算

4、條件后開始進行數(shù)值計算。3、模擬結果及分析由于流體有粘性，相鄰的流體層之間相互作用，高速流體層帶動低速流體層，低速流體層阻礙高速流體層。由于流體粘性系數(shù)較大，雷諾數(shù)變化較小，但依然從流線圖中可以看出，在Re=100時中心渦在方腔的中間的位置，Re=300時，中心漩渦有向右上方移動的趨勢，隨著雷諾數(shù)的增加，中心漩渦的位置越來越接近運動方向偏上的位置。3.1 方腔內(nèi)流線函數(shù)的分布情況 Re= 100 Re=300 Re=500 Re=700Re=9003.2 tecplot后處理 對Fluent輸出的速度題圖的數(shù)據(jù)進行后處理得到的圖像與上述分析相同的結果。見下圖 Re依次為100、300、500、

5、700、9004、結論方腔驅(qū)動流是數(shù)值計算中比較簡單，具有驗證性的一種流動情況，受到很多研究者的關注。本文通過不同雷諾數(shù)觀察方腔流動，所得結論如下：（1）當邊界以一定的速度移動時，方腔內(nèi)會產(chǎn)生一個較大的漩渦流，漩渦中心靠近運動方向的一側(cè)，雷諾數(shù)相對較小的時候，方腔渦中心偏向于方腔的中間，隨著雷諾數(shù)增加面中心渦，逐漸向右上方移動。（2）隨著雷諾數(shù)的增加，漩渦向方腔右上方移動。（3）方腔壁面上的速度大于其他地方的速度。參考文獻：【1】楊晶.基于Simple 算法的方腔驅(qū)動流問題數(shù)值模擬.電力學報，2010,01【2】E. Ertuk, T.C.Corke and C. Gokcol. Numerical Solutions of 2-D Steady Incompressible Driven Cavity Flow at High Reynolds Numbers, Int.J.Numer.Meth.Fluids,2005, vol.48,pp.747-774 【3】Ronald J. Adrian. Particle-imaging techn