![二維圓柱繞流的數(shù)值模擬_第1頁(yè)](http://file4.renrendoc.com/view/1759b203b6620354f7413e93bc1fccfe/1759b203b6620354f7413e93bc1fccfe1.gif)
![二維圓柱繞流的數(shù)值模擬_第2頁(yè)](http://file4.renrendoc.com/view/1759b203b6620354f7413e93bc1fccfe/1759b203b6620354f7413e93bc1fccfe2.gif)
![二維圓柱繞流的數(shù)值模擬_第3頁(yè)](http://file4.renrendoc.com/view/1759b203b6620354f7413e93bc1fccfe/1759b203b6620354f7413e93bc1fccfe3.gif)
![二維圓柱繞流的數(shù)值模擬_第4頁(yè)](http://file4.renrendoc.com/view/1759b203b6620354f7413e93bc1fccfe/1759b203b6620354f7413e93bc1fccfe4.gif)
下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
二維圓柱繞流的數(shù)值模擬
0大試現(xiàn)象的應(yīng)用一個(gè)世紀(jì)以來(lái),圓曲流問(wèn)題一直是經(jīng)典水流問(wèn)題之一,也是許多理論分析、數(shù)值模擬和實(shí)驗(yàn)研究的對(duì)象。盡管研究歷經(jīng)時(shí)間如此之長(zhǎng),但這種最簡(jiǎn)單的流動(dòng)現(xiàn)象物理本質(zhì)的理解依然不完整。閔強(qiáng)利計(jì)算了2種雷諾數(shù)Re=200,800情況下的圓柱繞流,發(fā)現(xiàn)在Re=200時(shí)產(chǎn)生的卡門(mén)渦街非常規(guī)則,而隨著雷諾數(shù)的增加,Re=800時(shí),渦街出現(xiàn)不規(guī)則現(xiàn)象,并表現(xiàn)出三維效應(yīng)。文獻(xiàn)研究了不同截面形狀在低雷諾數(shù)下的繞流場(chǎng),并重點(diǎn)研究了阻力、升力的周期性變化特性。隨著計(jì)算機(jī)硬件的不斷發(fā)展,大渦模擬作為一種新的湍流模型相比原來(lái)的湍流模型在求解非穩(wěn)態(tài)流動(dòng)方面具有巨大的潛力。大渦模擬方法的基本思想是用瞬時(shí)的N-S方程直接模擬計(jì)算湍流中的大尺度渦,而小尺度渦對(duì)大渦的影響則通過(guò)建立近似的模型來(lái)考慮,這種影響模型稱(chēng)為亞格子尺度模型。王漢青等人介紹了大渦模擬的理論進(jìn)展和發(fā)展趨勢(shì),描述了當(dāng)前大渦模擬在工程中的具體應(yīng)用。指出大渦模擬在模擬計(jì)算從層流到湍流轉(zhuǎn)換、非定常湍流和高速湍流方面具有其他湍流模型無(wú)可比擬的優(yōu)勢(shì)。傅慧萍以潛艇模型SUBOFF為研究對(duì)象,采用大渦模擬(LES)方法求解了流動(dòng)的非定常解。通過(guò)與試驗(yàn)值以及采用RNGk-ε湍流模型得到的定常結(jié)果比較,驗(yàn)證了大渦模擬方法的有效性。本文采用大渦模擬方法求解低雷諾數(shù)下的圓柱繞流場(chǎng),研究卡門(mén)渦街的周期性特性。對(duì)一個(gè)完整的渦街脫落周期進(jìn)行詳細(xì)討論,揭示了不同時(shí)刻旋渦的生成、發(fā)展、脫落和演化過(guò)程。通過(guò)研究圓柱表面的壓力變化說(shuō)明產(chǎn)生周期性脫落渦對(duì)的原因,以及產(chǎn)生周期性變化的升、阻力的原因。1亞品格尺度模型自然界中的流動(dòng)現(xiàn)象基本可以通過(guò)基于連續(xù)性假設(shè)的納維葉-斯托克斯方程(N-S方程)來(lái)描述。在應(yīng)用于大渦模擬方法時(shí),需要對(duì)N-S方程進(jìn)行一定過(guò)濾處理。?ˉui?t+?ˉuiˉuj?xj=-1ρ?ˉp?xi+v?2ˉui?xj?xj-?τij?xj?uˉi?t+?uˉiuˉj?xj=?1ρ?pˉ?xi+v?2uˉi?xj?xj??τij?xj,(1)?ˉui?xi=0?i?{1,2,3}?uˉi?xi=0?i?{1,2,3}。(2)對(duì)于亞格子尺度模型,采用標(biāo)準(zhǔn)Smagorinsky模型的方法。τij-13τkkδij=-2vtˉSij=-2(CSˉΔ)2∥S∥ˉSijτij?13τkkδij=?2vtSˉˉij=?2(CSΔˉˉˉ)2∥S∥Sˉˉij。(3)式中:S為應(yīng)變速率張量,S=12(?uj?xi+?ui?xj)S=12(?uj?xi+?ui?xj),‖S‖為進(jìn)行濾過(guò)后的應(yīng)變速率張量,∥S∥=√2ˉSijˉSij∥S∥=2SˉˉijSˉˉij???????√;vt為渦粘度,vt=l2‖S‖。本文計(jì)算所采用的網(wǎng)格均為四邊形結(jié)構(gòu)網(wǎng)格,濾過(guò)范圍可取為ˉΔ=2vol13Δˉˉˉ=2vol13。(4)式中,vol是計(jì)算單元的體積。亞格子尺度可表示為l=CsˉΔ(1-exp(-y+?25)3)1/2l=CsΔˉˉˉ(1?exp(?y+/25)3)1/2。(5)式中,Smagorinsky常量Cs取為0.255。2計(jì)算雷諾數(shù)和流動(dòng)本文所采用的計(jì)算模型為直徑D=0.04m的圓柱。來(lái)流速度為U=0.01m/s,計(jì)算的雷諾數(shù)Re=400。計(jì)算域?yàn)榉叫?沿著來(lái)流方向長(zhǎng)度為7.5D,尾流去流段長(zhǎng)度為40D,兩側(cè)寬度為7.5D。流動(dòng)在選定的計(jì)算域內(nèi)能夠充分發(fā)展,圓柱周?chē)捎眠吔鐚泳W(wǎng)格,能夠捕捉圓柱周?chē)募?xì)微流動(dòng)。圖1為計(jì)算域示意圖,圖2為圓柱周?chē)木W(wǎng)格劃分。3計(jì)算與分析3.1實(shí)驗(yàn)結(jié)果與實(shí)驗(yàn)結(jié)果的比較表1為計(jì)算所得的斯特勞哈爾數(shù)(St)以及阻力系數(shù)(Cd)與文獻(xiàn)中實(shí)驗(yàn)結(jié)果的對(duì)比。計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合良好,說(shuō)明采用本文的數(shù)值計(jì)算方法是準(zhǔn)確而可靠的。3.2u3000學(xué)習(xí)前后建設(shè)內(nèi)自由振動(dòng)的分布圖3~圖8為圓柱繞流非穩(wěn)態(tài)流動(dòng)中一個(gè)典型的渦街脫落周期中不同時(shí)刻圓柱周?chē)膲毫?chǎng)分布。反映了一個(gè)典型渦街脫落周期內(nèi)漩渦的初生、形成、脫落及在尾流中不斷發(fā)展的整個(gè)過(guò)程。從圖3可看出,t=0時(shí)刻,上一周期脫落形成的渦在圓柱尾流中逐漸發(fā)展,圓柱的右下方約45°的位置形成一片明顯的負(fù)壓區(qū),旋渦在負(fù)壓區(qū)內(nèi)逐漸生成。圖4中t=1/5T時(shí)刻,圓柱右下方的漩渦完全生成,沿著尾流方向移動(dòng),并有逐步脫落的趨勢(shì)。而圖5中t=2/5T時(shí)刻,圓柱右下方形成的漩渦完全脫落并沿著尾流方向發(fā)展。圖6中t=3/5T時(shí)刻,之前由圓柱右下方脫落的漩渦在尾流中逐漸發(fā)展,同時(shí)在圓柱的右上方約45°的位置形成一片負(fù)壓區(qū),旋渦在負(fù)壓區(qū)內(nèi)逐漸生成。圖7中t=4/5T時(shí)刻圓柱右上方的漩渦完全生成,沿著尾流方向移動(dòng),并有逐漸脫落的趨勢(shì)。最后圓柱右上方形成的漩渦在t=T時(shí)刻完全脫落(圖8),脫落的漩渦在尾流中不斷發(fā)展,改變著圓柱周?chē)膲毫?chǎng)分布。圓柱右上方區(qū)域和右下方區(qū)域形成的漩渦交替脫落,在尾流中不斷發(fā)展,形成了典型的卡門(mén)渦街。3.3圓柱升力、阻力的變化圓柱在繞流場(chǎng)中受到的主要水動(dòng)力為阻力和升力。通過(guò)無(wú)量綱化,可以用阻力系數(shù)Cd及升力系數(shù)CL來(lái)表征。Cd=D12ρU2LCd=D12ρU2L,(6)CL=L12ρU2L。(7)式中:D和L分別為圓柱所受的阻力和升力;U為來(lái)流速度;ρ為流體密度;L為特征長(zhǎng)度,本文取為圓柱直徑。圖9為計(jì)算求取的穩(wěn)定的非穩(wěn)態(tài)計(jì)算時(shí)間內(nèi)圓柱的阻力系數(shù)、升力系數(shù)時(shí)程曲線。計(jì)算時(shí)時(shí)間步長(zhǎng)取為0.05s。計(jì)算所得的結(jié)果是在非穩(wěn)態(tài)計(jì)算了較長(zhǎng)時(shí)間,流動(dòng)發(fā)展相對(duì)較為穩(wěn)定時(shí)提取出來(lái)的。由圖9可知阻力和升力都呈現(xiàn)周期性的變化,這是由于周期性脫落的卡門(mén)渦街引起的。周期性脫落的漩渦造成了圓柱周?chē)鷫毫?chǎng)的不斷變化。升力系數(shù)基本在0左右波動(dòng),而阻力系數(shù)基本在1.4左右波動(dòng),這說(shuō)明阻力對(duì)圓柱的作用非常強(qiáng),升力對(duì)圓柱的作用一般。文獻(xiàn)指出,升力以斯特勞哈爾渦泄頻率fs變化,阻力以2倍渦泄頻率2fs變化。由圖9可以看出,阻力變化的頻率是升力變化頻率的2倍。進(jìn)一步對(duì)升力做功率譜密度分析,可以得到渦泄頻率fs,圖10為升力功率譜密度。圖10中的尖峰對(duì)應(yīng)為斯特勞哈爾渦泄頻率fs,其值在0.05Hz附近,與文中之前計(jì)算的斯特勞哈爾數(shù)(St=fs*D/U)是對(duì)應(yīng)的。為了深入研究圓柱上阻力和升力脈動(dòng)的變化,可以研究圓柱瞬態(tài)的升力、阻力系數(shù)。非穩(wěn)態(tài)阻力系數(shù)可以定義為C′d(t)=12∫2π0C′p(θ,t)cosθdθ;(8)非穩(wěn)態(tài)升力系數(shù)可以定義為C′l(t)=12∫2π0C′p(θ,t)sinθdθ。(9)角度θ的定義見(jiàn)圖11;圖12為典型的渦街脫落周期中不同時(shí)刻圓柱周向的壓力系數(shù)變化。圖12可看出,圓柱壓力系數(shù)在1個(gè)周期內(nèi)變化基本是對(duì)稱(chēng)的。圓柱駐點(diǎn)處(θ=0°和θ=360°)壓力系數(shù)Cp有最大值,在1附近;隨著來(lái)流向圓柱兩側(cè)擴(kuò)展,壓力系數(shù)迅速減小。在θ=80°和θ=280°附近,流動(dòng)產(chǎn)生分離,圓柱駐點(diǎn)后方的壓力在不同時(shí)刻呈現(xiàn)出較為規(guī)律的周期性變化。對(duì)比t=1/5T和t=4/5T時(shí)刻可以看出,在t=1/5T時(shí)刻圓柱的右下方區(qū)域(180°~270°)的壓力大于右上方區(qū)域(90°~180°),此時(shí)圓柱具有最大的負(fù)升力;在t=4/5T時(shí)刻圓柱的右下方區(qū)域(90°~180°)的壓力大于右下方區(qū)域(180°~270°),此時(shí)圓柱具有最大的正升力。這樣升力就完成了1個(gè)渦街脫落周期中由最大正升力向最大負(fù)升力的轉(zhuǎn)變。上述過(guò)程不斷進(jìn)行,就形成了升力的周期性變化。4大渦流場(chǎng)模擬分析采用大渦模擬的方法研究了雷諾數(shù)Re=400的圓柱繞流場(chǎng),準(zhǔn)確捕捉了卡門(mén)渦街這一經(jīng)典流動(dòng)現(xiàn)象,同時(shí)通過(guò)對(duì)水動(dòng)力特性的研究,得出以下結(jié)論:1)采用大渦模擬的方
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年抗?jié)B碳磚項(xiàng)目可行性研究報(bào)告
- 2025年多槽超聲波清洗機(jī)項(xiàng)目可行性研究報(bào)告
- 2025年印染用助劑項(xiàng)目可行性研究報(bào)告
- 2025年低溫空調(diào)項(xiàng)目可行性研究報(bào)告
- 2025至2030年中國(guó)高密度聚乙烯管材數(shù)據(jù)監(jiān)測(cè)研究報(bào)告
- 2025至2030年鍍前處理酸性除油劑項(xiàng)目投資價(jià)值分析報(bào)告
- 2025至2030年警燈三輪童車(chē)項(xiàng)目投資價(jià)值分析報(bào)告
- 2025至2030年福美胂可濕性粉劑項(xiàng)目投資價(jià)值分析報(bào)告
- 2025至2030年火炮半軸模鍛件項(xiàng)目投資價(jià)值分析報(bào)告
- 2025至2030年桿套項(xiàng)目投資價(jià)值分析報(bào)告
- 四年級(jí)語(yǔ)文下冊(cè)第六單元【集體備課】(教材解讀+教學(xué)設(shè)計(jì))
- 2024版義務(wù)教育小學(xué)科學(xué)課程標(biāo)準(zhǔn)
- 蘇教版小學(xué)信息技術(shù)五年級(jí)下冊(cè)五年級(jí)下冊(cè)教案全集
- 學(xué)校托管工作方案
- 腎性高血壓的護(hù)理查房
- 蘇教版八年級(jí)數(shù)學(xué)上冊(cè)期末試卷及答案【完美版】
- 法院拍賣(mài)議價(jià)協(xié)議書(shū)
- 2021年人教版八年級(jí)物理上冊(cè)期末考試卷(完美版)
- TB 10009-2016 鐵路電力牽引供電設(shè)計(jì)規(guī)范
- 2024年?yáng)|南亞雞蛋分級(jí)包裝設(shè)備市場(chǎng)深度研究及預(yù)測(cè)報(bào)告
- 2023高考數(shù)學(xué)藝考生一輪復(fù)習(xí)基礎(chǔ)講義(學(xué)生版)
評(píng)論
0/150
提交評(píng)論