百葉窗二維穩(wěn)態(tài)流動(dòng)與換熱周期性充分發(fā)展的數(shù)值模擬匯總

1、數(shù)值傳熱學(xué)大作業(yè)百葉窗二維穩(wěn)態(tài)流動(dòng)與換熱周期性充分發(fā)展的數(shù)值模擬院 系：XXXXXXXXXXXX姓 名：XXXXXX 班 級(jí)：XXXXXXXXX 學(xué) 號(hào)：XXXXXXXXXX 授課老師：XXXXXX 百葉窗二維穩(wěn)態(tài)流動(dòng)與換熱周期性充分發(fā)展的數(shù)值模擬（西安交通大學(xué)，710049，西安）摘要：針對(duì)百葉窗的翅片進(jìn)行數(shù)值模擬，在恒溫條件下，采用SIMPLER算法，對(duì)換熱和對(duì)流進(jìn)入周期性發(fā)展階段的翅片進(jìn)行數(shù)值模擬計(jì)算，得到摩擦因數(shù)隨雷諾數(shù)10到500的變化，模擬結(jié)果證明：翅片為25°時(shí)，當(dāng)雷諾數(shù)的變大時(shí)，模擬速度變大，平均Nu數(shù)也隨著變大，但摩擦因數(shù)變小。關(guān)鍵詞：百葉窗、傳熱、穩(wěn)態(tài)流動(dòng)、數(shù)值模

2、擬Abstract: Regarding to the louver fin to carry on the numerical simulation, under the condition of constant temperature, useing the SIMPLER algorithm, exchange heat convection finned into periodic stage of development to carry on the numerical simulation, the change of the friction factor with Reyn

3、olds number 10 to 500, the simulation results prove that: fin is 25 °, when Reynolds number gets bigger, the simulation speed and the average Nu number also with bigger, but smaller friction factor.Keywords: the shutters, heat transfer, steady flow, the numerical simulation表1：本文所用符號(hào)及名稱單位符號(hào)名稱單位摩

4、擦因數(shù)平均Nusselt數(shù)PrPrandtl（普朗特）數(shù)雷諾數(shù)熱擴(kuò)散系數(shù)（定義）表面換熱系數(shù)W/(m2)導(dǎo)熱系數(shù)W/(m)溫度平均溫度翅片壁面溫度方向上的速度m/s翅片長(zhǎng)度mx方向周期長(zhǎng)度my方向周期長(zhǎng)度mP壓力Paq熱流密度W/m2角度o1 引言百葉窗式翅片具有切斷散熱器上氣體邊界層的發(fā)展、減小邊界層厚度、提高散熱器性能的作用 漆 波,李隆鍵,崔文智,陳清華. 百葉窗式翅片換熱器中的耦合傳熱J. 重慶大學(xué)學(xué)報(bào)(自然科學(xué)版). 2005，28（10）：39-41.。國(guó)外不少學(xué)者對(duì)百葉窗式翅片換熱器的流動(dòng)結(jié)構(gòu)、傳熱與阻力特性進(jìn)行了大量的實(shí)驗(yàn)研究 MARLOW E E, THOLE K A. Ex

5、perimentalDesign forFlow-field Studies of Louvered FinsJ.Experimental Thermaland Fluid Science, 1998, (18): 258 -269.百葉窗式翅片的實(shí)驗(yàn)研究包括對(duì)各種百葉窗翅片散熱器的傳熱性能和流動(dòng)特性的整體檢測(cè)以及對(duì)其流動(dòng)的可視化研究。國(guó)內(nèi)外的學(xué)者通過對(duì)其整體傳熱和流動(dòng)性能的檢測(cè)，積累了大量的試驗(yàn)數(shù)據(jù)，推出了相應(yīng)的經(jīng)驗(yàn)公式，對(duì)百葉窗翅片式換熱器在工程中的應(yīng)用和推廣起到了極大的作用 漆波. 百葉窗翅片式散熱器內(nèi)對(duì)流-導(dǎo)熱禍合傳熱的數(shù)值模擬D. 重慶大學(xué).2005年。而對(duì)其內(nèi)部流動(dòng)形態(tài)（主要是氣側(cè)

6、）的分析更有助于了解其傳熱強(qiáng)化的機(jī)理和特性。對(duì)于百葉窗而言，為了加強(qiáng)換熱，采用把空間翅片布置成周期排列并呈一定的角度來達(dá)到這一目的，選區(qū)其中一部分翅片作為研究對(duì)象如圖1方框中部分，針對(duì)翅片的模擬以Fortran語(yǔ)言編程，采用SIMPLER算法模擬計(jì)算10-500的雷諾數(shù)范圍內(nèi)恒溫的百葉窗翅片板換熱特性和流動(dòng)。以模擬計(jì)算推導(dǎo)雷諾數(shù)對(duì)摩擦因數(shù)和平均Nu數(shù)的影響規(guī)律。圖1 模擬部分表2 文中所用數(shù)據(jù)參數(shù)L1TpLpDelta數(shù)值30mm301.525單位mm18.6mmmm°2 SIMPLER算法及計(jì)算步驟1972年SIMPLE算法問世以來在世界各國(guó)計(jì)算傳熱學(xué)及計(jì)算流體力學(xué)界得到了廣泛的應(yīng)

7、用，這種算法提出不久很快就成為計(jì)算不可壓流場(chǎng)的主要方法，后來SIMPLE算法以及其后的各種改進(jìn)方案成功的推廣到可壓縮流場(chǎng)計(jì)算中，已成為一種可以計(jì)算任何流速的流動(dòng)的數(shù)值方法?；炯僭O(shè)：速度場(chǎng)的假定與壓力場(chǎng)的假定各自獨(dú)立進(jìn)行，二者無任何聯(lián)系。對(duì)假定壓力場(chǎng)的修正通過已求解的速度場(chǎng)的質(zhì)量守恒條件得到。中間速度通過求解當(dāng)前壓力得到，如果求解速度不能滿足質(zhì)量守恒條件，對(duì)壓力添加一個(gè)修正量修正，速度場(chǎng)也隨之得以修正。第二假設(shè)：在做速度修正時(shí)，忽略不同位置的速度修正量之間的影響。采用SIMPLER算法實(shí)施速度分量和壓力代數(shù)方程的分離式求解時(shí)，計(jì)算步驟如下：a、假定一個(gè)速度分布，記為u0,v0,w0，以此計(jì)算動(dòng)

8、量離散方程中的系數(shù)及常數(shù)項(xiàng)；b、假設(shè)一個(gè)壓力場(chǎng)p0 ；c、依次求解動(dòng)量方程，得 u1,v1,w1；d、對(duì)壓力加以修正，得p1 ；e、根據(jù)p1改進(jìn)速度值；f、利用改進(jìn)后的速度場(chǎng)求解那些通過源項(xiàng)物性等與速度場(chǎng)耦合的變量，如果變量并不影響流場(chǎng)，則應(yīng)在速度場(chǎng)收斂后再求解；g、利用改進(jìn)后的速度場(chǎng)重新計(jì)算動(dòng)量離散方程的系數(shù)，并利用改進(jìn)后的壓力場(chǎng)作為下一層次迭代計(jì)算的初值。重復(fù)上述步驟，直到獲得收斂的解。3 劃分網(wǎng)格及數(shù)值模擬在一般情況下，為了化簡(jiǎn)模型，模擬都要采用特定條件下的假設(shè)，本文對(duì)百葉窗翅片的模擬同樣也采用一部分假設(shè)：a 假設(shè)流體中沒有體積力和粘性耗散；b 常物性；c 流動(dòng)換熱處在穩(wěn)態(tài)、層流、周期循

9、環(huán)充分發(fā)展的階段。截取百葉窗的一個(gè)翅片為研究模型，對(duì)翅片進(jìn)行數(shù)值模擬，然后再擴(kuò)展到整個(gè)百葉窗模型，研究所選翅片如圖1所示的方框區(qū)域，網(wǎng)格數(shù)取50×50。圖2 網(wǎng)格劃分模擬計(jì)算的控制方程 陶文銓編著，數(shù)值傳熱學(xué)，第二版，西安，西安交通大學(xué)出版社，2012年，9-10.：動(dòng)能方程：Tx+Ty=a2TX2+2TY2 ()連續(xù)方程：x+y=0 ()動(dòng)量方程：x+y=-1px+2x2+2y2 () x+y=-1py+2x2+2y2 () x,y=Tx,y-TwTbx-Tw ()邊界條件： Tx,Tp=Tx,0 () x,Tp=x,0 () x,Tp=x,0 () 0,y=Lp,y () 0,y

10、=Lp,y () 0,y=Lp,y 對(duì)于邊界上的節(jié)點(diǎn)值采用線性插值的方法：對(duì)于Y方向有： i，1=i，M1=*i,M22+*i,2 在公式（2）中*表示上一次迭代，M1、M2表示y向最后一行和從后面數(shù)第二行的網(wǎng)格。通過前面y向第二行和從后面數(shù)第二行的迭代數(shù)值來替代模擬區(qū)域上面和下面的邊界條件再進(jìn)行下一次的迭代計(jì)算。在X方向有： =12Tx2-2,y-TwTbx2-2-Tw+Tx1-1,y-TwTbx1-1-Tw 輸入和輸出的溫度的更新值以下面的公式來計(jì)算：對(duì)于x方向上的A-B和C-D所代表的邊界區(qū)域圖2所示。其中ABCD所選模擬計(jì)算區(qū)域。平均速度由公式（6）決定： T0,y=Tw+Tb0-Tw

11、 TL,y=Tw+TbL-Tw 圖3 Tbx=TpTx,yx,ydyTpx,ydy 模擬計(jì)算中采用SIMPLER算法。控制方程使用有限容積法(Finite Volume Method,DVM)離散，以乘方來求解離散方程中擴(kuò)散項(xiàng)和對(duì)流的系數(shù)。以搜尋方式邊界的節(jié)點(diǎn)編號(hào)號(hào)，因?qū)λx翅片部分的處理，在進(jìn)行迭代計(jì)算時(shí)，該部分中的速度和溫度分別取為0、Tw，擴(kuò)散系數(shù)在所選部分中的值取為/Pr，然后進(jìn)行一萬次迭代。又因百葉窗的翅片與流體間的熱交換與平均Nu數(shù)有關(guān)，所以平均Nu數(shù)通過垂直壁面數(shù)以予確定如下： N=hmL1 hm=qmTbL-Tb0lnTw-Tb0-lnTw-TbL qm=RmGpT出-T入 f

12、=pm0-pmLpm22 其中：Rm是流入質(zhì)量流量，Gp是流體比熱容模擬數(shù)據(jù)結(jié)果與討論：相關(guān)雷諾數(shù)下的熱交換數(shù)據(jù)與溫度場(chǎng)和速度場(chǎng)（每隔5°C取一條等溫線）。模擬計(jì)算的溫度場(chǎng)： 圖4 雷諾數(shù)為30 圖5雷諾數(shù)為40 圖6雷諾數(shù)為100 圖7雷諾數(shù)為200 圖8雷諾數(shù)為300 圖9雷諾數(shù)為500模擬計(jì)算的速度場(chǎng)： 圖10雷諾數(shù)為10 圖11雷諾數(shù)為15 圖12雷諾數(shù)為20 圖13雷諾數(shù)為100 圖14雷諾數(shù)為200 圖15雷諾數(shù)為500 從上面不相同的雷諾數(shù)在二維周期百葉窗翅片流道內(nèi)溫度場(chǎng)和速度場(chǎng)模擬結(jié)果來看，在雷諾數(shù)為20附近時(shí)，在百葉窗翅片部的右邊開始有漩渦產(chǎn)生，當(dāng)雷諾數(shù)逐漸增大是，

13、所出現(xiàn)的漩渦逐漸向左移動(dòng)，隨著雷諾數(shù)的值增加到200時(shí)，模擬計(jì)算的部分產(chǎn)生很大的渦流漩渦，所以模擬結(jié)果證明，當(dāng)雷諾數(shù)逐漸增大時(shí)流體流動(dòng)速度逐漸增大并且流體的擾動(dòng)也隨之增加。由圖中可以看出，在雷諾數(shù)增大的同時(shí)，百葉窗翅片附近的等溫線逐漸加密，故溫度梯度增強(qiáng)，模擬結(jié)果證明，當(dāng)速度逐漸增大時(shí)擾動(dòng)和熱交換提高。在不同雷諾數(shù)下的平均Nu數(shù)和摩擦系數(shù)的比較模擬結(jié)果如圖16所示：a在不同雷諾數(shù)下的平均Nu數(shù) b 在不同雷諾數(shù)下的摩擦系數(shù)值圖16 由圖16 a可以看出，百葉窗翅片平均Nu數(shù)隨著雷諾數(shù)的增大而增加。在常物性條件下雷諾數(shù)與流體的速度成正比，即速度隨雷諾數(shù)增大而提高，所以熱交換就越強(qiáng)，故此平均Nu數(shù)也會(huì)增大，所以較小的雷諾數(shù)時(shí)，如雷諾數(shù)為15時(shí)，出現(xiàn)較小的渦流漩渦，隨著雷諾數(shù)增加到100的時(shí)候，所選的模擬計(jì)算部分出現(xiàn)較大的渦流漩渦，平均Nu數(shù)的增長(zhǎng)率也提高不少，故此換熱率也相對(duì)提高。又如圖16 b可以看出，在雷諾數(shù)增大的同時(shí)，摩擦差系數(shù)卻在相對(duì)減小，其減小速率也在降低，當(dāng)雷諾數(shù)大于200時(shí)，摩擦系數(shù)減小速率已經(jīng)變得非常緩慢，引起此現(xiàn)象的因素為，當(dāng)流體渦流產(chǎn)生壓力降低，摩擦系數(shù)減小，使雷諾數(shù)的變化減緩。4 本文所得結(jié)論在百葉窗翅片為25°時(shí)，對(duì)其在穩(wěn)態(tài)、層流、周期性充分發(fā)展的條件下進(jìn)行模擬計(jì)算，以雷諾數(shù)為10，15，20，30，100，200，300，500下的