fluent的一個(gè)實(shí)例(波浪管道的內(nèi)部流動(dòng)模擬).

1、基于FLUENT的波浪管道熱傳遞耦合模擬CFD可以對熱傳遞耦合的流體流動(dòng)進(jìn)行模擬。 CFD模擬可以觀察到管道內(nèi) 部的流動(dòng)行為和熱傳遞，這樣可以改進(jìn)波浪壁面復(fù)雜通道幾何形狀中的熱傳遞。目的：（1）創(chuàng)建由足夠數(shù)量的完整波浪組成的波浪管道，提供充分發(fā)展條件；（2）應(yīng)用周期性邊界條件創(chuàng)建波浪通道的一部分；（3）研究不同湍流模型以及壁面函數(shù)對求解的影響；（4）采用固定表面溫度以及固定表面熱流量條件， 確定雷諾數(shù)與熱特性之間的 關(guān)系。問題的描述：通道由重復(fù)部分構(gòu)成，每一部分由頂部的直面和底部的正弦曲面構(gòu)成， 如圖圖1管道模型 空氣的流動(dòng)特性如下：質(zhì)量流量：m=0.816kg/s;密度：p =1kg/m3;

2、動(dòng)力粘度： 流動(dòng)溫度：卩=0.0001kg/(m s);Tb=300K ；流體其他熱特性選擇默認(rèn)項(xiàng)。流動(dòng)初試條件：x方向的速度=0.816m/s； 湍動(dòng)能=1m2/s2;湍流耗散率=1xi05m/s3。所有湍流模型中均采用增強(qiáng)壁面處理操作過程：-、完整波浪管道模型的數(shù)值模擬(1) 計(jì)算Re=uH/v=0.816 1/ (0.0001/1) =8160-0 2-0 2Cf/2=0.0359Re =0.0359 X8160)=0.00592595 =比號=0.816 匯 0.0059259 = 0.0628+y =uty/vy=0.00159(2) 創(chuàng)建網(wǎng)格入口本例為波浪形管道，管道壁面為我們所感

3、興趣的地方所以要局部細(xì)化。 和出口處的邊界網(wǎng)格設(shè)置如圖。圖2邊網(wǎng)格生成面網(wǎng)格圖3管道網(wǎng)格(3) 運(yùn)用Flue nt進(jìn)行計(jì)算本例涉及熱傳遞耦合，所以在flue nt中啟動(dòng)能量方程，如圖圖4能量方程設(shè)定條件，湍流模型選擇標(biāo)準(zhǔn)k-e模型，近壁面處理選擇增強(qiáng)壁面處理圖5湍流模型設(shè)定材料，密度為1,動(dòng)力粘度改為0.0001如圖圖6材料設(shè)定設(shè)定邊界條件，入口速度為0.816,湍動(dòng)能為1,湍流耗散率為100000。出 口為自由出口，壁面溫度為固定溫度分別為300k, 500k。圖7邊界條件初始化,并計(jì)算。R&siduials continuity XVelocity y-velocity energy kS

4、c-elte4 Rtiduel5Jn 0FLUENT 6.3 (2d. pb圖8殘差殘差中的e和k并沒有減小，沒有達(dá)到10-3一下，并且由于網(wǎng)格很大，計(jì)算 時(shí)間很長。I3.46e-0l2倜曲亦心佃血-2 016-01-3 376-014 74e-D1Hg.11a-01-7 畑 01-1 (Ka*D0-1 16e*D0-1 29fi*00-1 43e*00 1.57ig*00-1 ?0e*00-1 B4a*D0-1 905*00-2.110*00-2 25e*00-2 3*00ContcHTB of Total Pnsssure (pascal)Jan 07, 2012FLUENT 6 3 (2

5、d. pbmi tto)圖9壓力分布圖隨著流體流動(dòng)，管道中壓力分布趨于平穩(wěn)，波浪管道中波谷的壓力最低，在 入口處的壓力較高。圖10速度分布圖從圖中可以看出，在管道 7-11個(gè)波浪處，流動(dòng)已經(jīng)充分發(fā)展，貼近上壁處 速度最大，在波谷出的速度最小，甚至接近于零。1.9024盹畑4 75*02 屯畑礎(chǔ)4 55&+024 4te+(K?4 36*+0?44閔遊4.07e+023.97e*02s.are-os3 7B&+023.6Sa0e3.58a*0fi3.3帥訛3.25fr+0e3d脅礎(chǔ)3.000+02Contaur& of Total Temperature (k)Jan 07. 2012FLUEN

6、T S3 (2d. pbs. ek)圖11溫度分布圖管道后部流貼近波浪壁面出的溫度較高，流動(dòng)充分發(fā)展后，由于換熱作用, 體溫度逐漸升高，在波峰與波谷之間的流體溫度最高，如圖。UOaOO1.53*00 1J5e+00 1.3600 1 2E&+M 1 1!te+001 Ilift+OG&.37e-GlB. 5-017.67e-O1B aie-015.96a-0l5.11e-0l1.?Ge-ai2.56&-011.710018&4a-Oe2.50a-04Jan 07. 2012FLUENT S3 (2d. pbs. ek)圖12速度矢量圖可以看出，波浪壁面出流體出現(xiàn)反流，在波谷出反流的流體最多，速

7、度在波 谷出最小，接近于0,出現(xiàn)滯留區(qū)。若要觀察波峰、波谷處流體流動(dòng)速度，需要在波峰、波谷處創(chuàng)建兩條直線， 觀察直線上的速度。因?yàn)楣艿?7-11節(jié)處流動(dòng)充分發(fā)展，所以在第十一節(jié)波峰、 波谷處建立兩條直線，如圖。一 crest-w圖13波峰、波谷V&l&city MagnitudeJan 07, 20112FLUENT 6.3 (2d. pbns.ske)圖14波峰、波谷的速度波峰貼近兩側(cè)壁面出的速度梯度很大，在管道中間速度隨高度增加而增加， 在0.7m左右達(dá)到最大。波谷處靠近上壁面的速度梯度很大，但是由于有波谷存 在，波谷處的速度梯度不大，在谷中速度先增大再減小，在0處左右達(dá)到最下，隨后逐漸增

8、加，在0.7m高度左右速度最大。高度在0.5m處以上波峰波谷處的速 度基本相等。= =蘭一 三=_=二 二三三= - =-=-_ 一三一ii -=i_= = _= -=i-圖15周期性網(wǎng)格.r.f.-m周期性波浪管道模型的數(shù)值模擬網(wǎng)格密度與完全管道網(wǎng)格相同。在flue nt中輸入以下指令，創(chuàng)建周期性網(wǎng)格。/grid modify-z ones/grid/modify-z on es make-periodicPeriodic zone () inletShadow zone () outletRotati onal periodic? (if no, tran slati on al) yes

9、 noCreate periodic zon es? yes yesAuto detect tran slati on vector? yes yescomputed tran slation deltas: 1.000000 0.000000 all 100 faces matched for zones 6 and 5.zone 5 deletedcreated periodic zon es./grid/modify-z on es邊界條件中可以看到已經(jīng)沒有outlet，inlet也變成了 periodic周期性的。這里要設(shè)置周期性邊界條件。質(zhì)量流量為0.816,其他設(shè)置與之前相同。f

10、PericxJic ConditionsTypeFlow Direction存 Specify Mass FlowSpecify Pressure Gradienlx1Mass Flow Rate |kg/s|9,816lzPressure Gradient pascal/m|Helaxalion Factor|d.sUpstreaim Bulk Temperature |khlumber of Iteratians阿2引OK | Update I Cancel | Help圖16邊界條件Residusiiscontinuity jc-vetacity y-velocityle+oale-03

11、Scaled ResidualsJan 07,2012FLUMT6.$ (2d ftn* 知時(shí)圖17殘差Energy并不收斂，反而隨計(jì)算而發(fā)散。其他參數(shù)都收斂。計(jì)算量很小，計(jì) 算速度明顯提高。2 41*41 Z23tl|1 SMiI 4B* 011 IOa-OI帖切住 3他 4B-7-i-I Mh-QI- Kb-012 Zta-aiMMH-1 Ua411JBQqXII J Tl、4伽亦-ICanliHjn of Total PmuB ijpaualiFLUENTTa. 3l2d.pbrw skij圖18壓力圖彳晦畑燈如02屮嶺*fl?心世4.17b*i02 w*oe3険*強(qiáng) 生內(nèi)小佢 工如竝

12、matted 3.5ta 磁 m&KJ?3 314MI BGfrOq1 jreii丁 qswh2 43311 95#-0l1 47*1glHQQ 1 W5JVflkrciTy Vactera C咼rvd By Vekreity Magn耐圖18壓力圖速度圖溫度圖速度矢量圖管道中間壓力最大并比較均勻，在波浪波谷出壓力最低。溫度分布不太正確圖19波峰波谷速度圖波峰貼近兩側(cè)壁面出的速度梯度很大，與完全模型相似，波峰在管道中間速度隨高度增加先增加在減小，但是變化不大，其他分布趨勢與完全模型相似。IerestpI匚*耳&Velocity MagnitudeJan 07, 2012FLUENT 6.3 (

13、2d. pbns.sto)圖20周期模型與完全模型波峰速度圖比較周期性模型，管道中速度基本相等，沒有呈現(xiàn)完全模型的變化趨勢。圖21周期模型與完全模型波谷速度圖比較兩模型的波谷速度的變化趨勢大致相同，在小于0.3m處，周期性模型速度大于完全模型，0.3m之后，周期性模型速度小于完全管道模型。三、 周期性波浪管道不同湍流模型的數(shù)值模擬應(yīng)用周期模型,如圖。Icrasfptrouah-DVelocityMagnitude fm/s)再選分別應(yīng)用 RNG與Realizable湍流模型計(jì)算。實(shí)驗(yàn)數(shù)據(jù)Jm07, 2012FLUENT 8 3 (2d. pbcis. mgkA)圖21 RNG模型波峰波谷速度圖

14、圖22 Realizable模型波峰波谷速度圖將之前的數(shù)據(jù)放在一起進(jìn)行比較。VelocityMagnitudeVeloc rty MagnrKideJan 07, 2012FLUENT 6.3 (2d.兩事 rhk圖23不同湍流模型波峰速度圖trough-p trough-rea Irauah-fnaVelocityMagnitude(m劇Velocity MogrortudeJaR 07.3012FLUtWT .3 (2d, ptns. rrg帕)cresl-fng cresRea orisibo圖24不同湍流模型波谷速度圖從圖中可以看出realizable和標(biāo)準(zhǔn)k-e模型波峰的速度分布大致

15、相同，RNG模型與完全模型更接近，變化更大。三種湍流模型波谷的速度分布與波峰分布特 點(diǎn)相似。四、整體模型壁面采用固定熱流率的數(shù)值模擬重新選用整體模型，只需改變波浪面即 wall2的邊界條件即可，改為熱流率 1000w/m2，其他條件不變。圖25不同湍流模型波谷速度圖VelocityMagnitude(m/s)Velacrfy Magnitude-Jan 07. 2012FLUENT 6 3 (2d. pbm.圖26固定壁面溫度與固定熱流率波峰速度圖比較圖27固定壁面溫度與固定熱流率波谷速度圖比較從圖中可以看出固定壁面溫度與固定熱流率波峰與波谷的速度分布還是比 較相似的，不同的是在0.6m處固定熱流率波峰波谷速度比固定壁面溫度的速度 稍高，在0.6m之后固定熱流率波峰波谷速度比固定壁面溫度的速度稍低，在接 近兩壁面處基本相等。2倜曲射-2.010-01-3 3333114 75e-D1l卻11X1寸畑01-ft B5e-0l-1 026*00-1 16e*D0-11 30e*00-