fluent爐膛仿真教程文檔

1、爐膛仿真過程及其其中的問題一、（Gambit）幾何建模部分1.大體尺寸在本次設(shè)計中，（實際標高-5=圖中的標高）鍋爐的尺寸為：鍋爐高度為26890mm，寬度為7570mm，深度為7570mm。燃燒器的高度為2.105m，最底層的燃燒器低端距冷灰斗距離為2.1775m。采用四角切圓（順時針切圓，假想切圓直徑0.8m）的均等配風燃燒方式。其中一次風2層，二次風3層。由低到高燃燒器風口布置依次為二、一、二、一、二。燃燒器寬度為0.4m，一次風口高度0.2405m，二次風口高度0.352/0.315m，風口間距為0.21/0.12/0.155m。2.簡化處理將水冷壁簡化成一個恒溫平面；將燃燒器簡化成一

2、個平面，各次風口為平面中的一個矩形區(qū)域，作為速度入口；忽略屏式過熱器，將折焰角上方與水平煙道相連結(jié)的平面作為出口（outflow）。3.幾何建模過程及網(wǎng)格劃分為了方便鍋爐的網(wǎng)格劃分，我們將整個計算域劃分為5個區(qū)域：冷灰斗下端至燃燒區(qū)域下端、燃燒區(qū)域、燃燒區(qū)上端至折焰角下端、折焰角區(qū)域、折焰角上端至爐膛出口。3.1點線面的生成幾何建模的方法通?？梢允亲韵露系?，即先生成體的各個點（通過坐標確定位置）；將生成的點依次連接成線；將線圍成體的各個面；最后將面組合成一個實體。當然建模時也可以通過設(shè)置實體（面）的長寬高（長寬）直接生成。3.2實體分割塊的劃分方法如下：先產(chǎn)生一個面，并將該面平移至該實體要切

3、割的位置，split volume選卡中，split with選擇face（real），然后選中要切割的實體（對應(yīng)split volume中的volume）以及用來切割這個體的面（對應(yīng)face欄）（注意：在切割時需要選中Connected，保證切割產(chǎn)生的兩個體之間的面是公共面，而不是兩個重合的面。因為公共面可以通過物質(zhì)和能量，而重合的面不加定義時是wall），最后點擊APPLY確定。 根據(jù)這種方法，我們可以在Z方向?qū)⑷紵齾^(qū)分為很多層，方便以后設(shè)置一、二次風入口的邊界條件。同時，在xy平面內(nèi)燃燒區(qū)被分為8份，如圖所示：3.3網(wǎng)格劃分網(wǎng)格劃分的最后記結(jié)果如圖所示：這種網(wǎng)格的特點是：四個角的地方網(wǎng)格

4、比較密，而中間網(wǎng)格比較稀疏。同時網(wǎng)格線的方向與流動合速度方向重合度比較高。這樣的網(wǎng)格劃分可以很好的抑制偽擴散的發(fā)生。這種網(wǎng)格的劃分步驟如下：在將區(qū)域分塊的基礎(chǔ)上對實體按照線、面、體的順序進行依次劃分。Edge的劃分：為了形成這種對稱的網(wǎng)格，我們需要對edge進行劃分，如圖所示：其中，soft link采用maintain形式，Spacing選用Interval count（劃分數(shù)目）。在本設(shè)計中，將每條線均分為30份，即ratio為1，interval count為30，其他保持不變。 Face的劃分：由于前面對每條邊進行了劃分，所以對面的網(wǎng)格劃分就只需要設(shè)置網(wǎng)格的形式和類型如圖所示：其中，E

5、lements采用Quad形式，Type采用Map形式（映射成結(jié)構(gòu)化網(wǎng)格）。此時不需要對Spacing進行設(shè)置了。Volume的劃分：對volume的劃分，我們采用Cooper（制桶）方式。采用這種劃分方式時，有一點需要注意，就是上下兩個Face的網(wǎng)格劃分要完全一樣，也就是說組成Face的Edge的劃分也要一樣。如圖所示：其中，Element采用Hex/Wedge形式。Sources表示需要選擇制“桶”的上下兩個面。Interval count表示兩個面之間劃分的數(shù)目。本設(shè)計中，根據(jù)風口和墻面的高度進行劃分，每個網(wǎng)格高度在0.1m左右。最后依照上面的方法和步驟對燃燒區(qū)的每一層進行這樣的網(wǎng)格劃分

6、。對于除了燃燒器區(qū)的其他區(qū)域的網(wǎng)格劃分，要求就比較低一些了。對我們依舊采用COOPER的方式對體進行劃分。不過其他地方的Sources是沿y軸方向的兩個面（燃燒器區(qū)域的sources是沿z方向的）。最后的網(wǎng)格為：冷灰斗30*30*30；燃燒器30*30*8*（3+2+3+1+3+1+3+2+3）；燃燒器上端至折焰角：50*50*66；折焰角：50*50*16；折焰角上方：40*50*40。最后網(wǎng)格數(shù)目大概在480000個，其中燃燒器區(qū)域網(wǎng)格為151200個。3.4交接面處的處理在劃分計算域的時候會涉及到interface的設(shè)置。在燃燒區(qū)的上下兩個端面，我們需要分別將這個面與其相重合的那個面設(shè)置

7、成一對interface。因為燃燒器區(qū)與相鄰的兩個實體并不是通過分割而來，是3個獨立的實體，為了能讓物質(zhì)和能量通過該重合的面，需要通過設(shè)置interface來實現(xiàn)，如圖所示：由于燃燒器區(qū)域上端的xy平面被劃分為了8塊，所以需要將這8個面一起設(shè)置為interface11，然后將與燃燒器區(qū)上端重合的面設(shè)置為interface12。對于interface21和22的設(shè)置和上述一樣。4.邊界條件設(shè)置在gambit中需要預(yù)先設(shè)置邊界條件。將折焰角上方與水平煙道相連接的那個面設(shè)置為outflow邊界條件。Interface的設(shè)置上面已經(jīng)說過了，下面我們進行一、二次風入口的設(shè)置。根據(jù)燃燒器的結(jié)構(gòu)確定各次風口

8、在模型中的位置，然后將邊界條件的Type設(shè)置成velocity_inlet（速度入口）。名字格式為ofa/pa/sa+兩位數(shù)字，數(shù)字前一位表示在xy平面所處的象限，后一位表示自高向低同類型風口的層數(shù)。如圖所示：二、Fluent仿真過程0.網(wǎng)格導入、Interface設(shè)置以及網(wǎng)格檢查在完成Gambit中的工作后，需要將生成的.msh文件導入到Fluent中。0.1網(wǎng)格導入、檢查以及解法器設(shè)置在General中點擊Check完成網(wǎng)格檢查（網(wǎng)格檢查中不能出現(xiàn)網(wǎng)格體積為負數(shù)的情況，否則會出錯，需要重新進行稽核建模）。點擊Report Quality進行網(wǎng)格質(zhì)量檢查。在解法器中選擇Pressure-Ba

9、sed、Absolute、Steady的情況。勾選Gravity，建立重力場（z=-9.81m/s2），設(shè)置如圖所示：0.2Interface設(shè)置點擊Mesh Interfaces中的Create。在Interface Zone 1中選擇interface11，在Interface Zone2中點選interface12，Mesh Interface名稱為interface1，點擊Create設(shè)置完成。按照同樣的方法設(shè)置interface2，如圖所示：1.燃料及邊界條件參數(shù)確定1.1燃料特性及風煤計算1.1.1燃料計算工況符號單位3#T-017#T-018#T-01適用標準全水分Mt%9.18

10、.78.2GB/T211-2007空氣干燥基水分Mad%2.531.822.17GB/T212-2008收到基灰分Aar%27.8733.2236.39干燥無灰基揮發(fā)分Vdaf%41.0144.0039.92收到基碳Car%51.3647.5045.49DL/T568-1995收到基氫Har%3.493.433.09收到基氮Nar%0.820.760.72收到基氧Oar%6.676.175.86全硫St,ar%0.690.220.25GB/T214-2007收到基高位發(fā)熱量Qgr,v,arMJ/kg20.4819.1818.19GB/T 213-2008收到基低位發(fā)熱量Qnet,v,arMJ/

11、kg19.5518.2717.36根據(jù)表格，我們將元素分析數(shù)據(jù)轉(zhuǎn)換成干燥無灰基的揮發(fā)分的元素組成。由于干燥無灰基無水、無灰，故剩下的成分不受水分和灰分的影響，是表示碳、氫、氧、氮、硫成分百分函數(shù)最穩(wěn)定額基準，所以通常選擇轉(zhuǎn)換為干燥無灰基來計算。各種煤不同分析基之間的換算公式為X=X0×K其中，X0，X分別為某成分原基準與新基準的質(zhì)量百分數(shù)，%；K為換算系數(shù)。收到基與干燥無灰基之間的換算系數(shù)為K=100100-Mar-Aar利用上述公式將煤的收到基轉(zhuǎn)換為干燥無灰基的元素組成，如下表所示，這些數(shù)據(jù)也是Fluent中計算PDF元素分數(shù)所需要的。PDF中元素分析元素分數(shù)元素CHON元素組成0

12、.8209710.0557660.1057570.017506PDF中工業(yè)分析分析組分組分VFCAM組分分數(shù)0.22120.33290.36390.0821.1.2風煤計算鍋爐實際燃煤量t/h26.015設(shè)計值一次風流速m/s24.3設(shè)計值一次風份額%28.28設(shè)計值一次風溫K303溫風份額%64.09設(shè)計值二次風溫K600一次風口面積m24*0.1924二次風口面積m24*0.4075注：由于不知道乏氣送粉的位置，將乏氣份額歸并到二次風中，即二次風份額為71.72%。根據(jù)克拉伯龍方程PV=nRT和表格數(shù)據(jù)可知：標況下一次風速V10=21.894ms；再根據(jù)一、二次風的份額和面積可以得到：二次

13、風速V20=V10*二、一次風份額之比*一、二次風口面積之比=26.216ms；根據(jù)克拉伯龍方程可知：實際二次風速V2=57.62ms；根據(jù)煤量，可以知道每個一次風口煤的質(zhì)量流量：Qcoal=26000÷3600÷4÷2=0.903kgs根據(jù)切圓直徑和爐膛尺寸可以知道風煤的入口方向：夾角=40.55°；cos=0.76；sin=0.65.1.2邊界條件設(shè)置現(xiàn)以一次風pa11為例，介紹對流場數(shù)據(jù)的設(shè)置。在Fluent中Boundary Conditions菜單下找到pa11項目，如圖所示。然后這個風口進行設(shè)置。點擊Edit進入設(shè)置頁面，如圖所示。在Velo

14、city Specification Method選項中選擇Magnitude and Direction（速度大小和方向）；在Velocity Magnitude中填入24.3m/s；在Coordinate System中選擇Cartesian（X,Y,Z）笛卡爾直角坐標系，然后在下面依次填入流體流動的方向（X軸為-sin40.55°，Y軸為-cos40.55°）。在第一象限的風口方向為（-sin40.55°，-cos40.55°）；在第二象限的為（cos40.55°，-sin40.55°）；第三象限的為（sin40.55°

15、;，cos40.55°）；第四象限的為（-cos40.55°，sin40.55°）。在設(shè)置湍流參數(shù)時，我們選用Intensity and Hydraulic Diameter（湍流強度和水力直徑）方式。湍流強度I我們設(shè)置成10%，為強湍流狀態(tài)，水力直徑D的設(shè)置根據(jù)公式D=4AL設(shè)置成0.3m（二次風口的水力直徑為0.374m/0.352m）。在Species選卡中將Mean Mixture Fraction設(shè)置為0（氧化劑入口），Mixture Fraction Variance設(shè)置為0，如圖所示。在DPM選卡中將DP BC Type設(shè)置成reflect（反射）其

16、他設(shè)置保持不變。依照上面的方法，可以完成對PA和SA流場參數(shù)設(shè)置。2.與流動和燃燒相關(guān)的模型設(shè)置在Fluent中我們打開Models選項。在中意菜單中，我們可以設(shè)置包括流動、傳熱、燃燒等方面的模型。2.1.氣相流動模型本文采用標準k-e雙方程湍流流動模型，同時采用標準壁面函數(shù)處理近壁面的流動問題，其中的流動參數(shù)保持默認參數(shù)。其設(shè)置如下：2.2氣相湍流燃燒模型模擬氣相湍流燃燒過程的關(guān)鍵在于如何?；牧魅紵磻?yīng)率。針對擴撒火焰的模型有k-g模型，混合數(shù)-概率密度函數(shù)模型。為了減少計算量，采用但混合數(shù)PDF模型。對于煤粉燃燒，我們在Models->Species中選用Non-Premixed

17、Combustion（非預(yù)混燃燒）模型，然后再彈出的菜單中進行相關(guān)的參數(shù)設(shè)置，如圖所示。在PDF Table Creation欄目中點擊Chemistry選卡。在State Relation中選擇Equilibrium（化學平衡法），Energy Treatment中選擇Non-Adiabatic（非絕熱）形式，Stream Options選擇Empirical Fuel Stream（經(jīng)驗燃料流）。點擊Coal Calculator會彈出對燃料特性進行設(shè)置的對話框，根據(jù)燃料特性表的中數(shù)據(jù)，可以設(shè)置完燃料工業(yè)分析和元素分析的參數(shù)。其中物料的名稱為coal-particle，HCV為1.819e

18、+07j/kg，其他的保持默認數(shù)值，點擊Apply和OK確認，我們可以看到Model Settings中的數(shù)據(jù)發(fā)生了相應(yīng)的變化。如圖所示。點擊Boundary選卡將燃料溫度設(shè)置為303K，氧化劑的溫度設(shè)置為600K。點擊Table選卡中的Calculate PDF Table進行燃燒的計算。計算完成后，我們可以點擊Display PDF Table查看關(guān)于煤粉燃燒的數(shù)據(jù)，如圖所示。計算結(jié)果顯示燃燒形成的成分有20種，點擊Control選卡可以查看成分名稱，如圖所示。設(shè)置完成后點擊Apply和OK。同時我們可以看出Models中的Species變成了Non-Premixed Combustion

19、模型。2.3煤粉燃燒模型煤粉燃燒可分為煤粉預(yù)熱、揮發(fā)分析出和燃燒過程、焦炭燃燒等過程。在本設(shè)計中，揮發(fā)分析出模型采用單速率析出模型，焦炭燃燒模型選用擴散-動力控制模型。在Fluent中首先需要對煤粉顆粒的噴射進行相關(guān)的設(shè)置。在Models菜單中點擊Discrete Phase，彈出對話框，如圖所示。在Interaction中勾選Interaction with Continuous Phase（對連續(xù)相的影響，但是在仿真時候需要先建立無顆粒相的流場，即在仿真開始時不勾該選項）。在Trackking選項中，Max. Number of Steps設(shè)置為3000，勾選Specify Length

20、Scale，其中，Length Scale設(shè)置為0.01m；在Physical Models選項和Numerics選項中的參數(shù)和選項均為默認設(shè)置。為了讓煤粉能夠噴射入爐膛，我們需要進行對顆粒相的設(shè)置，點擊Discrete Phase Model中的Injections彈出對話框，如圖所示。點擊Create對煤粉顆粒進行參數(shù)設(shè)置。以某一次風口進入的煤粉為例，如圖所示。在Injection Type中選擇surface，然后在Release From Surface中選擇相應(yīng)的煤粉噴射地點。在Particle Type中勾選Combusting（燃燒）。在設(shè)置Diameter Distributi

21、on時，可以選擇uniform（均勻）形式，當然最好選擇rosin-rammler形式。在Point Properties中設(shè)置好速度方向（與空氣速度矢量一致）、顆粒直徑（1e-06m）、顆粒溫度（303K）、質(zhì)量流量（0.903kg/s）等參數(shù)；在Turbulent Dispersion中的Stochastic Tracking選擇Discrete Random Walk Model（隨機軌道模型），如圖所示。依照上述方法對其他的煤粉噴射源進行參數(shù)和模型的設(shè)置選擇。煤粉噴射設(shè)置完后可以發(fā)現(xiàn)Models->Materials選卡中多了Combusting Particle（coal-pa

22、rticle）一項。接下來我們對這個成分進行相關(guān)設(shè)置，如圖所示。在彈出的對話框中的Properties中設(shè)置Devolatilization Model（揮發(fā)分析出模型）為single-rate（單速率）；同時設(shè)置Combustion Model（焦炭燃燒模型）為kinetic/diffusion-limited（擴散-動力控制）；對于該顆粒相的其他物性參數(shù)保持默認參數(shù)。2.3輻射傳熱模型爐內(nèi)的能量主要通過輻射的形式進行傳遞。Fluent軟件提供了多種輻射模型，在本設(shè)計中，我們選用P1模型，如圖所示。在Model中點選P1即可。P1模型考慮了輻射散射作用，更適用于光學厚度較厚以及幾何結(jié)構(gòu)復雜的

23、燃燒設(shè)備。但是該模型也有缺點，包括對來自內(nèi)部熱源的輻射熱通量有過高估計的趨勢等。在選用P1模型后，Models中Energy選卡自動開啟。2.4邊界條件的設(shè)置在Boundary Conditions中設(shè)置Supersonic/Initial Gauge Pressure為0pascal；一次風入口溫度為303K；二次風入口溫度為600K；DP BC Type為reflect。墻壁為固定無滑移壁面；熱力條件為很穩(wěn)壁面600K，內(nèi)部發(fā)射率為1，壁厚及生熱率均為0；壁面材料為Al（鋁）；DPM 邊界條件數(shù)類型BC Type為reflect。其他保持默認設(shè)置。出口邊界條件設(shè)置為outflow，保持默認

24、設(shè)置。3求解過程（Solution）3.1Solution Methods點選Solution->Solution Methods。壓力-速度耦合采用SIMPLE格式；空間離散方法中，梯度采用格林-高斯單元法（Green-Gauss Cell Based）；壓力采用Standard方法；其他均采用二階迎風格式（Second Order Upwind）。3.2Solution Initialization初始化方法采用Standard Initialization方法；初始溫度為600K；其他保持默認參數(shù)和設(shè)置，點擊Initialize完成初始化，如圖所示：3.3求解在求解時，需要先建立連

25、續(xù)相的流場，然后再加入顆粒相耦合修正流場，即在開始時Models->Discrete Phase中不勾選Interaction with Continuous Phase，迭代100步先建立連續(xù)相流場；迭代完成后，勾選Interaction with Continuous Phase，其他保持默認設(shè)置，然后進行迭代，如圖所示：三、仿真結(jié)果中出現(xiàn)的問題在迭代步數(shù)為2670步基礎(chǔ)上的仿真結(jié)果1.爐膛內(nèi)切圓直徑與假想切圓直徑相差不大，說明有可能迭代步數(shù)不夠，如圖所示：2.在燃燒器上端至出口的區(qū)域內(nèi)，溫度變化不大，所散點圖所示：3.現(xiàn)場測試數(shù)據(jù)顯示乙側(cè)溫度都高于甲側(cè)溫度，根據(jù)經(jīng)驗表明是順時針切圓

26、；但在仿真中采用順時針切圓，結(jié)果卻是乙側(cè)溫度低于甲側(cè)溫度，標高24.2m處溫度情況如圖所示：標高27.5m處溫度情況如圖所示：4.各標高所在平面溫度等值圖4.1下二次風4.2下一次風4.3中二次風4.4上一次風4.5上二次風4.6標高24.2處溫度X=-1.685處Y坐標（m） 溫度（K）-3.785600-3.63361449.03-3.48221511.38-3.33081522.37-3.17941530.99-3.0281538.02-2.87661543.94-2.72521549-2.57381553.35-2.42241557.08-2.2711560.24-2.11961562

27、.91-1.96821565.17-1.81681567.08-1.66541568.69-1.5141570.07-1.36261571.25-1.21121572.24-1.05981573.07-0.90841573.77-0.7571574.35-0.60561574.83-0.45421575.22-0.30281575.52-0.15141575.75-7.35393e-181575.910.151415760.30281576.020.45421575.960.60561575.780.7571575.470.90841575.041.05981574.441.21121573.631.36261572.611.5141571.331.66541569.761.81681567.851.96821565.542.11961562.782.2711559.532.42241555.752.57381551.372.72521546.332.87661540.573.0281533.933.17941526.223.33081517.183.48221506.053.63361449.113.7856004.6標高27.5處溫度X=-1.467m處Y坐標（m） 溫度（K）-3.785600-3.63361436