fluent爐膛仿真教學(xué)教程文檔

.爐膛仿真過程及其其中的問題一、（Gambit）幾何建模部分1.大體尺寸在本次設(shè)計中，（實際標(biāo)高-5=圖中的標(biāo)高）鍋爐的尺寸為：鍋爐高度為26890mm，精品文檔放心下載寬度為7570mm，深度為7570mm。燃燒器的高度為2.105m，最底層的燃燒器低端距冷灰斗距離為2.1775m。謝謝閱讀采用四角切圓（順時針切圓，假想切圓直徑0.8m）的均等配風(fēng)燃燒方式。其中一次風(fēng)謝謝閱讀2層，二次風(fēng)3層。由低到高燃燒器風(fēng)口布置依次為二、一、二、一、二。燃燒器寬度為謝謝閱讀.0.4m，一次風(fēng)口高度0.2405m，二次風(fēng)口高度0.352/0.315m，風(fēng)口間距為謝謝閱讀0.21/0.12/0.155m。2.簡化處理將水冷壁簡化成一個恒溫平面；將燃燒器簡化成一個平面，各次風(fēng)口為平面中的一個矩形區(qū)域，作為速度入口；感謝閱讀忽略屏式過熱器，將折焰角上方與水平煙道相連結(jié)的平面作為出口（outflow）。感謝閱讀.3.幾何建模過程及網(wǎng)格劃分為了方便鍋爐的網(wǎng)格劃分，我們將整個計算域劃分為5個區(qū)域：冷灰斗下端至燃燒區(qū)感謝閱讀域下端、燃燒區(qū)域、燃燒區(qū)上端至折焰角下端、折焰角區(qū)域、折焰角上端至爐膛出口。感謝閱讀3.1點線面的生成幾何建模的方法通?？梢允亲韵露系?，即先生成體的各個點（通過坐標(biāo)確定位置）；精品文檔放心下載將生成的點依次連接成線；將線圍成體的各個面；最后將面組合成一個實體。謝謝閱讀當(dāng)然建模時也可以通過設(shè)置實體（面）的長寬高（長寬）直接生成。精品文檔放心下載3.2實體分割塊的劃分方法如下：先產(chǎn)生一個面，并將該面平移至該實體要切割的位置，splitvolume選卡中，splitwith謝謝閱讀選擇face（real），然后選中要切割的實體（對應(yīng)splitvolume中的volume）以及用來精品文檔放心下載切割這個體的面（對應(yīng)face欄）（注意：在切割時需要選中Connected，保證切割產(chǎn)生的感謝閱讀兩個體之間的面是公共面，而不是兩個重合的面。因為公共面可以通過物質(zhì)和能量，而重精品文檔放心下載合的面不加定義時是wall），最后點擊APPLY確定。精品文檔放心下載.根據(jù)這種方法，我們可以在Z方向?qū)⑷紵齾^(qū)分為很多層，方便以后設(shè)置一、二次風(fēng)入感謝閱讀口的邊界條件。同時，在xy平面內(nèi)燃燒區(qū)被分為8份，如圖所示：感謝閱讀3.3網(wǎng)格劃分網(wǎng)格劃分的最后記結(jié)果如圖所示：.這種網(wǎng)格的特點是：四個角的地方網(wǎng)格比較密，而中間網(wǎng)格比較稀疏。同時網(wǎng)格線的方感謝閱讀向與流動合速度方向重合度比較高。這樣的網(wǎng)格劃分可以很好的抑制偽擴散的發(fā)生。謝謝閱讀這種網(wǎng)格的劃分步驟如下：在將區(qū)域分塊的基礎(chǔ)上對實體按照線、面、體的順序進行依次劃分。精品文檔放心下載Edge的劃分：為了形成這種對稱的網(wǎng)格，我們需要對edge進行劃分，如圖所示：謝謝閱讀.其中，softlink采用maintain形式，Spacing選用Intervalcount（劃分?jǐn)?shù)目）。在精品文檔放心下載本設(shè)計中，將每條線均分為30份，即ratio為1，intervalcount為30，其他保持不變。謝謝閱讀Face的劃分：由于前面對每條邊進行了劃分，所以對面的網(wǎng)格劃分就只需要設(shè)置網(wǎng)格精品文檔放心下載的形式和類型如圖所示：.其中，Elements采用Quad形式，Type采用Map形式（映射成結(jié)構(gòu)化網(wǎng)格）。此時精品文檔放心下載不需要對Spacing進行設(shè)置了。Volume的劃分：對volume的劃分，我們采用Cooper（制桶）方式。采用這種劃分謝謝閱讀方式時，有一點需要注意，就是上下兩個Face的網(wǎng)格劃分要完全一樣，也就是說組成Face感謝閱讀的Edge的劃分也要一樣。如圖所示：.其中，Element采用Hex/Wedge形式。Sources表示需要選擇制“桶”的上下兩個謝謝閱讀面。Intervalcount表示兩個面之間劃分的數(shù)目。本設(shè)計中，根據(jù)風(fēng)口和墻面的高度進行劃精品文檔放心下載分，每個網(wǎng)格高度在0.1m左右。最后依照上面的方法和步驟對燃燒區(qū)的每一層進行這樣的網(wǎng)格劃分。感謝閱讀對于除了燃燒器區(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è)置成一對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)說過了，下面我們進行一、二次風(fēng)入口的設(shè)置。根據(jù)燃燒器謝謝閱讀的結(jié)構(gòu)確定各次風(fēng)口在模型中的位置，然后將邊界條件的Type設(shè)置成velocity_inlet（速精品文檔放心下載度入口）。名字格式為ofa/pa/sa+兩位數(shù)字，數(shù)字前一位表示在xy平面所處的象限，后一感謝閱讀位表示自高向低同類型風(fēng)口的層數(shù)。如圖所示：二、Fluent仿真過程0.網(wǎng)格導(dǎo)入、Interface設(shè)置以及網(wǎng)格檢查.在完成Gambit中的工作后，需要將生成的.msh文件導(dǎo)入到Fluent中。謝謝閱讀0.1網(wǎng)格導(dǎo)入、檢查以及解法器設(shè)置在General中點擊Check完成網(wǎng)格檢查（網(wǎng)格檢查中不能出現(xiàn)網(wǎng)格體積為負(fù)數(shù)的情況，感謝閱讀否則會出錯，需要重新進行稽核建模）。點擊ReportQuality進行網(wǎng)格質(zhì)量檢查。精品文檔放心下載在解法器中選擇Pressure-Based、Absolute、Steady的情況。勾選Gravity，建立重感謝閱讀力場（z=-9.81m/s2），設(shè)置如圖所示：0.2Interface設(shè)置點擊MeshInterfaces中的Create。在InterfaceZone1中選擇interface11，在精品文檔放心下載InterfaceZone2中點選interface12，MeshInterface名稱為interface1，點擊Create精品文檔放心下載設(shè)置完成。按照同樣的方法設(shè)置interface2，如圖所示：感謝閱讀.1.燃料及邊界條件參數(shù)確定1.1燃料特性及風(fēng)煤計算1.1.1燃料計算工況符號單位3#T-017#T-018#T-01適用標(biāo)準(zhǔn)全水分Mt%GB/T211-2007空氣干燥基水分Mad%2.531.822.17收到基灰分Aar%27.8733.2236.39GB/T212-2008干燥無灰基揮發(fā)分Vdaf%41.0144.0039.92收到基碳Car%51.3647.5045.49收到基氫Har%3.493.433.09DL/T568-1995收到基氮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/T213-2008收到基低位發(fā)熱量Qnet,v,arMJ/kg19.5518.2717.36根據(jù)表格，我們將元素分析數(shù)據(jù)轉(zhuǎn)換成干燥無灰基的揮發(fā)分的元素組成。由于干燥無灰感謝閱讀基無水、無灰，故剩下的成分不受水分和灰分的影響，是表示碳、氫、氧、氮、硫成分百分謝謝閱讀函數(shù)最穩(wěn)定額基準(zhǔn)，所以通常選擇轉(zhuǎn)換為干燥無灰基來計算。感謝閱讀各種煤不同分析基之間的換算公式為X=0×其中，X0，X分別為某成分原基準(zhǔn)與新基準(zhǔn)的質(zhì)量百分?jǐn)?shù)，%；K為換算系數(shù)。精品文檔放心下載收到基與干燥無灰基之間的換算系數(shù)為K=

100 ar?ar利用上述公式將煤的收到基轉(zhuǎn)換為干燥無灰基的元素組成，如下表所示，這些數(shù)據(jù)也是感謝閱讀Fluent中計算PDF元素分?jǐn)?shù)所需要的。PDF中元素分析元素分?jǐn)?shù)元素 C H O N元素組成 0.820971 0.055766 0.105757 0.017506精品文檔放心下載PDF中工業(yè)分析分析組分組分 V FC A M組分分?jǐn)?shù) 0.2212 0.3329 0.3639 0.082精品文檔放心下載.1.1.2風(fēng)煤計算鍋爐實際燃煤量t/h26.015設(shè)計值一次風(fēng)流速m/s24.3設(shè)計值一次風(fēng)份額%28.28設(shè)計值一次風(fēng)溫K303溫風(fēng)份額%64.09設(shè)計值二次風(fēng)溫K600一次風(fēng)口面積m24*0.1924二次風(fēng)口面積m24*0.4075注：由于不知道乏氣送粉的位置，將乏氣份額歸并到二次風(fēng)中，即二次風(fēng)份額為71.72%。精品文檔放心下載根據(jù)克拉伯龍方程PV=nRT和表格數(shù)據(jù)可知：標(biāo)況下一次風(fēng)速 =. ?；再根據(jù)一、二次風(fēng)的份額和面積可以得到：二次風(fēng)速 = ?二、一次風(fēng)份額之比?一、二次風(fēng)口面積之比=. ?；感謝閱讀根據(jù)克拉伯龍方程可知：實際二次風(fēng)速 =. ?；根據(jù)煤量，可以知道每個一次風(fēng)口煤的質(zhì)量流量：coal= ÷ ÷ ÷ =. ?根據(jù)切圓直徑和爐膛尺寸可以知道風(fēng)煤的入口方向：夾角θ=40.55°；cosθ=0.76；sinθ=0.65.謝謝閱讀1.2邊界條件設(shè)置.現(xiàn)以一次風(fēng)pa11為例，介紹對流場數(shù)據(jù)的設(shè)置。在Fluent中BoundaryConditions菜單下找到pa11項目，如圖所示。然后這個風(fēng)口感謝閱讀進行設(shè)置。點擊Edit進入設(shè)置頁面，如圖所示。.VelocitySpecificationMethod選項中選擇MagnitudeandDirection（速度大小和方向）；在VelocityMagnitude中填入24.3m/s；精品文檔放心下載CoordinateSystem中選擇Cartesian（X,Y,Z）笛卡爾直角坐標(biāo)系，然后在下面依次填入流體流動的方向（X軸為-sin40.55°，Y軸為-cos40.55°）。在第一象限的風(fēng)口方向為（-sin40.55°，-cos40.55°）；在第二象限的為（cos40.55°，-sin40.55°）；第三象限的為（sin40.55°，cos40.55°）；第四象限的為（-cos40.55°，sin40.55°）。謝謝閱讀在設(shè)置湍流參數(shù)時，我們選用IntensityandHydraulicDiameter（湍流強度和水力直徑）方式。湍流強度I我們設(shè)置成10%，為強湍流狀態(tài)，水力直徑D的設(shè)置根據(jù)公式D=4?謝謝閱讀設(shè)置成0.3m（二次風(fēng)口的水力直徑為0.374m/0.352m）。精品文檔放心下載Species選卡中將MeanMixtureFraction設(shè)置為0（氧化劑入口），MixtureFractionVariance設(shè)置為0，如圖所示。感謝閱讀.在DPM選卡中將DPBCType設(shè)置成reflect（反射）精品文檔放心下載其他設(shè)置保持不變。依照上面的方法，可以完成對PA和SA流場參數(shù)設(shè)置。謝謝閱讀2.與流動和燃燒相關(guān)的模型設(shè)置在Fluent中我們打開Models選項。在中意菜單中，我們可以設(shè)置包括流動、傳熱、精品文檔放心下載燃燒等方面的模型。2.1.氣相流動模型本文采用標(biāo)準(zhǔn)k-e雙方程湍流流動模型，同時采用標(biāo)準(zhǔn)壁面函數(shù)處理近壁面的流動問題，精品文檔放心下載其中的流動參數(shù)保持默認(rèn)參數(shù)。其設(shè)置如下：.2.2氣相湍流燃燒模型模擬氣相湍流燃燒過程的關(guān)鍵在于如何模化湍流燃燒反應(yīng)率。針對擴撒火焰的模型有精品文檔放心下載k-ε-g模型，混合數(shù)-概率密度函數(shù)模型。為了減少計算量，采用但混合數(shù)PDF模型。對于煤粉燃燒，我們在Models->Species謝謝閱讀中選用Non-PremixedCombustion（非預(yù)混燃燒）模型，然后再彈出的菜單中進行相關(guān)感謝閱讀.的參數(shù)設(shè)置，如圖所示。在PDFTableCreation欄目中點擊Chemistry選卡。在StateRelation中選擇謝謝閱讀Equilibrium（化學(xué)平衡法），EnergyTreatment中選擇Non-Adiabatic（非絕熱）形式，精品文檔放心下載StreamOptions選擇EmpiricalFuelStream（經(jīng)驗燃料流）。謝謝閱讀點擊CoalCalculator會彈出對燃料特性進行設(shè)置的對話框，根據(jù)燃料特性表的中數(shù)據(jù)，謝謝閱讀可以設(shè)置完燃料工業(yè)分析和元素分析的參數(shù)。其中物料的名稱為coal-particle，HCV為感謝閱讀1.819e+07j/kg，其他的保持默認(rèn)數(shù)值，點擊Apply和OK確認(rèn)，我們可以看到Model感謝閱讀Settings中的數(shù)據(jù)發(fā)生了相應(yīng)的變化。如圖所示。謝謝閱讀.點擊Boundary選卡將燃料溫度設(shè)置為303K，氧化劑的溫度設(shè)置為600K。精品文檔放心下載點擊Table選卡中的CalculatePDFTable進行燃燒的計算。計算完成后，我們可以點感謝閱讀擊DisplayPDFTable查看關(guān)于煤粉燃燒的數(shù)據(jù)，如圖所示。謝謝閱讀.計算結(jié)果顯示燃燒形成的成分有20種，點擊Control選卡可以查看成分名稱，如圖所精品文檔放心下載示。設(shè)置完成后點擊Apply和OK。同時我們可以看出Models中的Species變成了感謝閱讀Non-PremixedCombustion模型。謝謝閱讀.2.3煤粉燃燒模型煤粉燃燒可分為煤粉預(yù)熱、揮發(fā)分析出和燃燒過程、焦炭燃燒等過程。謝謝閱讀在本設(shè)計中，揮發(fā)分析出模型采用單速率析出模型，焦炭燃燒模型選用擴散-動力控制感謝閱讀模型。在Fluent中首先需要對煤粉顆粒的噴射進行相關(guān)的設(shè)置。感謝閱讀在Models菜單中點擊DiscretePhase，彈出對話框，如圖所示。精品文檔放心下載在Interaction中勾選InteractionwithContinuousPhase（對連續(xù)相的影響，但是謝謝閱讀在仿真時候需要先建立無顆粒相的流場，即在仿真開始時不勾該選項）。謝謝閱讀在Trackking選項中，Max.NumberofSteps設(shè)置為3000，勾選SpecifyLength精品文檔放心下載Scale，其中，LengthScale設(shè)置為0.01m；感謝閱讀.在PhysicalModels選項和Numerics選項中的參數(shù)和選項均為默認(rèn)設(shè)置。謝謝閱讀為了讓煤粉能夠噴射入爐膛，我們需要進行對顆粒相的設(shè)置，點擊 DiscretePhase謝謝閱讀Model中的Injections彈出對話框，如圖所示。謝謝閱讀點擊Create對煤粉顆粒進行參數(shù)設(shè)置。以某一次風(fēng)口進入的煤粉為例，如圖所示。感謝閱讀.在InjectionType中選擇surface，然后在ReleaseFromSurface中選擇相應(yīng)的煤粉精品文檔放心下載噴射地點。在ParticleType中勾選Combusting（燃燒）。感謝閱讀在設(shè)置DiameterDistribution時，可以選擇uniform（均勻）形式，當(dāng)然最好選擇感謝閱讀rosin-rammler形式。在PointProperties中設(shè)置好速度方向（與空氣速度矢量一致）、顆粒直徑（1e-06m）、感謝閱讀顆粒溫度（303K）、質(zhì)量流量（0.903kg/s）等參數(shù)；精品文檔放心下載在TurbulentDispersion中的StochasticTracking選擇DiscreteRandomWalk精品文檔放心下載Model（隨機軌道模型），如圖所示。.依照上述方法對其他的煤粉噴射源進行參數(shù)和模型的設(shè)置選擇。謝謝閱讀煤粉噴射設(shè)置完后可以發(fā)現(xiàn)Models->Materials選卡中多了CombustingParticle謝謝閱讀（coal-particle）一項。接下來我們對這個成分進行相關(guān)設(shè)置，如圖所示。感謝閱讀.在彈出的對話框中的Properties中設(shè)置DevolatilizationModel（揮發(fā)分析出模型）精品文檔放心下載為single-rate（單速率）；同時設(shè)置CombustionModel（焦炭燃燒模型）為謝謝閱讀kinetic/diffusion-limited（擴散-動力控制）；謝謝閱讀對于該顆粒相的其他物性參數(shù)保持默認(rèn)參數(shù)。2.3輻射傳熱模型爐內(nèi)的能量主要通過輻射的形式進行傳遞。Fluent軟件提供了多種輻射模型，在本設(shè)感謝閱讀計中，我們選用P1模型，如圖所示。.在Model中點選P1即可。P1模型考慮了輻射散射作用，更適用于光學(xué)厚度較厚以及幾何結(jié)構(gòu)復(fù)雜的燃燒設(shè)備。精品文檔放心下載但是該模型也有缺點，包括對來自內(nèi)部熱源的輻射熱通量有過高估計的趨勢等。精品文檔放心下載在選用P1模型后，Models中Energy選卡自動開啟。感謝閱讀2.4邊界條件的設(shè)置在BoundaryConditions中設(shè)置Supersonic/InitialGaugePressure為0pascal；一感謝閱讀次風(fēng)入口溫度為303K；二次風(fēng)入口溫度為600K；DPBCType為reflect。精品文檔放心下載墻壁為固定無滑移壁面；熱力條件為很穩(wěn)壁面600K，內(nèi)部發(fā)射率為1，壁厚及生熱率感謝閱讀均為0；壁面材料為Al（鋁）；DPM邊界條件數(shù)類型BCType為reflect。其他保持默認(rèn)設(shè)精品文檔放心下載置。出口邊界條件設(shè)置為outflow，保持默認(rèn)設(shè)置。3求解過程（Solution）.3.1SolutionMethods點選Solution->SolutionMethods。謝謝閱讀壓力-速度耦合采用SIMPLE格式；空間離散方法中，梯度采用格林-高斯單元法（Green-GaussCellBased）；感謝閱讀壓力采用Standard方法；其他均采用二階迎風(fēng)格式（SecondOrderUpwind）。精品文檔放心下載3.2SolutionInitialization謝謝閱讀初始化方法采用StandardInitialization方法；感謝閱讀初始溫度為600K；其他保持默認(rèn)參數(shù)和設(shè)置，點擊Initialize完成初始化，如圖所示：感謝閱讀.3.3求解在求解時，需要先建立連續(xù)相的流場，然后再加入顆粒相耦合修正流場，即在開始時精品文檔放心下載Models->DiscretePhase中不勾選InteractionwithContinuousPhase，迭代100步先謝謝閱讀建立連續(xù)相流場；迭代完成后，勾選InteractionwithContinuousPhase，其他保持默認(rèn)設(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)驗表明是順時針切圓；但在仿真精品文檔放心下載中采用順時針切圓，結(jié)果卻是乙側(cè)溫度低于甲側(cè)溫度，標(biāo)高24.2m處溫度情況如圖所示：感謝閱讀標(biāo)高27.5m處溫度情況如圖所示：.4.各標(biāo)高所在平面溫度等值圖4.1下二次風(fēng)4.2下一次風(fēng)4.3中二次風(fēng).4.4上一次風(fēng)4.5上二次風(fēng)4.6標(biāo)高24.2處溫度.X=-1.685處Y坐標(biāo)（m）溫度（K）-3.785 600-3.63361449.03-3.48221511.38-3.33081522.37-3.17941530.99-3.028 1538.02-2.87661543.94.-2.72521549-2.57381553.35-2.42241557.08-2.271 1560.24-2.11961562.91-1.96821565.17-1.81681567.08-1.66541568.69-1.514 1570.07-1.36261571.25-1.21121572.24-1.05981573.07-0.90841573.77-0.757 1574.35-0.60561574.83-0.45421575.22-0.30281575.52-0.15141575.75-7.35393e-18 1575.910.1514 15760.3028 1576.020.4542 150561575.780.757 1575.470.90841575.041.05981574.441.21121573.631.36261572.611.514 1571.331.66541569.761.81681567.851.96821565.542.11961562.782.271 1559.532.42241555.752.57381551.372.72521546.332.87661540.573.028 1533.933.17941526.223.33081517.183.48221506.053.63361449.113.785 600.4.6標(biāo)高27.5處溫度X=-1.467m處Y坐