夏建芳，劉 寒，王瀟杰

（中南大學機電工程學院，湖南長沙 410083）

【摘 要】 環(huán)式冷卻機是燒結(jié)礦料的重要冷卻設備，其設計理論均基于 20世紀 80年代日本以臺車高度 1.5 m的實驗研究結(jié)果，導致系列產(chǎn)品所有臺車高度均為1.5 m，且難以預測燒結(jié)礦物性參數(shù)偏離實驗工況下的環(huán)冷機真實冷卻效果。作者在研究燒結(jié)礦風冷熱交換機理和等效仿真模型基礎上，應用 VB.NET語言，基于CFD平臺，開發(fā)了環(huán)式鼓風冷卻機燒結(jié)礦冷卻過程自動仿真軟件。該軟件具備環(huán)冷機幾何參數(shù)、燒結(jié)礦物性參數(shù)、燒結(jié)礦初始溫度、風機風量等友好參數(shù)輸入界面，自動仿真燒結(jié)礦溫度、排出煙氣溫度隨時間的變化曲線并以數(shù)據(jù)文件和圖像顯示。測試結(jié)果與仿真結(jié)果對比表明，燒結(jié)礦冷卻過程仿真結(jié)果與實測結(jié)果誤差在5%以內(nèi)，滿足工程精度要求，可以預測燒結(jié)礦冷卻效果、指導環(huán)式冷卻機的產(chǎn)品設計。

【關鍵詞】 環(huán)式冷卻機；數(shù)值模擬；仿真軟件；CFD軟件；VB.NET

引言

鋼鐵冶煉是一個高能耗、高污染的產(chǎn)業(yè)，據(jù)調(diào)查，在鋼鐵冶煉的過程中，燒結(jié)礦冷卻工序能耗約占企業(yè)總能耗的10%~15%，是僅次于煉鐵的第二大耗能工序^[1-2]。環(huán)式冷卻機（以下簡稱環(huán)冷機）是燒結(jié)礦冷卻工序中的關鍵設備，其結(jié)構圖如下圖 1 所示，隨著我國鋼鐵行業(yè)節(jié)能減排要求的越來越高，如何增大環(huán)冷機產(chǎn)量，提高能源利用率成為了迫在眉睫的問題^[3-4] 。調(diào)研發(fā)現(xiàn)，制約我國環(huán)冷機設備發(fā)展的主要問題是目前國內(nèi)環(huán)冷機的設計，都是參照上世紀八十年代日本基于臺車高度1.5 m的測試數(shù)據(jù)及統(tǒng)計公式。正是因為如此，環(huán)冷機在設計的過程中，很多關鍵參數(shù)無法改變，比如，我國的環(huán)冷機不論尺寸大小，其臺車高度都統(tǒng)一為 1.5m，導致現(xiàn)在只能通過增大環(huán)冷機截面積的方式來增大產(chǎn)量，造成環(huán)冷機安裝場地需求越來越高。不僅如此，環(huán)冷機設計人員在設計完成之后也無法準確的預測環(huán)冷機的冷卻效果。由此可知，傳統(tǒng)的環(huán)冷機設計方法，已經(jīng)無法適應現(xiàn)代的鋼鐵行業(yè)需求。

針對環(huán)冷機，國內(nèi)外學者進行了一系列的研究。在結(jié)構改進方面張忠波等^[5]基于燒結(jié)機結(jié)構開發(fā)新型高效燒結(jié)環(huán)冷機，很大程度上降低設備漏風率，提升能源利用率；在冷卻效率方面夏建芳等^[6]提出冷卻能耗測試與計算方法，評價了多臺套燒結(jié)機、環(huán)冷機的實際能耗。在數(shù)值模擬方面，Caputo^[7]在評估環(huán)冷機燒結(jié)礦冷卻床的換熱時，提出了一維穩(wěn)態(tài)數(shù)學模型；Pelagagge^[8-9]等人提出了環(huán)冷機二維換熱模型，并通過實測燒結(jié)礦的數(shù)據(jù)進行了驗證。

Jang^[10]等提出了環(huán)冷機三維換熱簡易模型，用于模擬環(huán)冷機的換熱及湍流過程；Leong^[11]等基于多孔介質(zhì)-局部熱平衡理論，利用計算流體力學軟件Fluent模擬了環(huán)冷機內(nèi)的流動與換熱過程。以上的研究工作僅研究換熱行為，未能建立相應計算公式或計算軟件，即還無法用于指導環(huán)冷機設計。為了對環(huán)冷機冷卻效果進行有效預測，在設計階段即時調(diào)整環(huán)冷機相關設計參數(shù)，開發(fā)環(huán)冷機冷卻過程仿真軟件迫在眉睫。

作者首先分析了環(huán)冷機內(nèi)氣固換熱過程，基于氣體湍流方程、多孔介質(zhì)模型與局部非熱平衡能量方程以及臺車單元環(huán)向移動等效環(huán)冷機模型，應用VB.NET 語言，對環(huán)冷機數(shù)值模擬過程的參數(shù)化建模、自動數(shù)值模擬、結(jié)果后處理等功能進行封裝，隱式生成 Gambit 幾何建模與網(wǎng)格劃分命令流，F(xiàn)luent仿真計算與后處理命令流，生成了燒結(jié)礦冷卻過程自動仿真軟件，實現(xiàn)燒結(jié)礦冷卻過程的自動仿真，并通過與實測數(shù)據(jù)的對比驗證仿真結(jié)果的正確性。

1 仿真等效物理模型與數(shù)學模型的建立

1.1 環(huán)冷機仿真等效物理模型分析

本文以高效燒結(jié)環(huán)冷機為研究對象，結(jié)構展開圖如圖2所示。高效燒結(jié)環(huán)冷機可以分為回風段與鼓風機段兩個部分，主要區(qū)別在于余熱回收段冷卻風由余熱回風管道輸入，溫度約為 120 ℃，煙氣由集熱罩收集進行余熱利用，鼓風機段冷卻風為空氣由風機吹入，煙氣通過百葉窗直接排向大氣。

多功能高效環(huán)冷機臺車回轉(zhuǎn)半徑大，結(jié)構復雜，且臺車繞其回轉(zhuǎn)中心轉(zhuǎn)動動態(tài)冷卻，若要對整個環(huán)冷機進行動網(wǎng)格建模，將對仿真計算機性能要求過高，且仿真過程耗時極長，因此本文探究等效模型的解決方案。等效物理模型以臺車為研究對象，分別建立回風段、鼓風機段臺車等效物理模型，具體模型如圖3、圖4所示。

以圖3為例，為環(huán)冷機回風段臺車模型，由圖可知，臺車由下至上分為3個部分。最下部分為風箱，中間為燒結(jié)礦，上部為集熱罩與煙囪或其它排出煙氣裝置。

1.2 燒結(jié)礦冷卻過程數(shù)學模型

本文數(shù)值分析基于 FLUENT 軟件，對燒結(jié)礦冷 卻過程熱力學分析的重點即分析燒結(jié)礦與冷卻氣體換熱過程所涉及的控制方程。

1.2.1 氣體湍流方程

環(huán)冷機工作過程中冷卻氣體在燒結(jié)礦內(nèi)進行流動，因此氣體的流動為湍流，選用標準 k - ε 模型^[9] 。該模型中，k和ε對應的運輸方程為：

k方程：

式中：ρ代表密度；k代表湍動能；ε是湍流耗散率；u代表速度矢量；i、j代表坐標維數(shù)；S表示平均速度應變率張量；τ_tij代表雷諾應力；C_1ε、C_2ε 、σ_ε、σ_ε為常數(shù)。

1.2.2 多孔介質(zhì)方程

燒結(jié)礦內(nèi)存在很多微小孔隙，并且大多數(shù)孔隙之間能夠互相連通，因此燒結(jié)礦是多孔介質(zhì)模型。

其黏性阻力系數(shù)（1/α）與慣性阻力系數(shù)（C 2 ）可通過Ergun方程計算:

黏性阻力系數(shù)和慣性阻力系數(shù)的定義分別為：

式中：L 代表燒結(jié)礦厚度，D_p是平均粒徑，? 代表燒結(jié)礦孔隙率，Δp 為壓降，u為動力粘度。

1.2.3 局部非熱平衡方程

本文燒結(jié)礦與冷卻氣體間的熱導率與相差較大，燒結(jié)礦與空氣溫度始終存在溫度差，局部熱力學平衡的假設不成立。本文采用局部非熱力學平衡雙能量方程，將燒結(jié)礦溫度和冷卻空氣溫度視為兩個獨立變量，定義燒結(jié)礦的平均溫度為 Ts，冷卻空氣的平均溫度為 Tf ，以平均容積法為基礎建立能量方程：

對于氣相：

式中：k_f,eff為導熱系數(shù)，ρ_f是流體的密度，φ_f為內(nèi)熱源，α_sf為比表面積，k_s,eff是固體的導熱系數(shù)，ρ_s是固體密度，φ_f是固體骨架中內(nèi)熱源。

2 燒結(jié)環(huán)冷機自動仿真系統(tǒng)程序架構

2.1 程序功能設計

本程序基于 VB.NET 高級語言搭建平臺，將編制出集仿真前處理、計算、后處理于一體的，界面友好，操作簡便的，能夠?qū)Νh(huán)冷機的設計與優(yōu)化提供指導的燒結(jié)環(huán)冷機自動仿真系統(tǒng)。因此在功能方面，需要實現(xiàn)用戶從操作界面輸入基本參數(shù)，系統(tǒng)自動對數(shù)據(jù)進行分類與處理、自動調(diào)用 GAMBIT 建模與網(wǎng)格劃分、自動調(diào)用 FLUENT 進行仿真。仿真結(jié)束后，對仿真結(jié)果進行處理，輸出燒結(jié)礦與煙氣曲線圖、風箱流量圖、顯示仿真云圖、自動輸出并保存仿真結(jié)果表格。

2.2 程序架構

燒結(jié)環(huán)冷機自動仿真系統(tǒng)基于 VB.NET 平臺，對 FLUENT、GAMBIT 等仿真軟件進行封裝，包括前處理、仿真計算、后處理三大模塊。為了實現(xiàn)設計目標，達成設計要求，具體的流程圖如圖5。

2.3 交互設計

為了充分實現(xiàn)所需的功能，本系統(tǒng)主要設計了兩類操作界面，分別為參數(shù)輸入界面與仿真界面。

如圖6所示。

在圖6參數(shù)輸入界面中，左側(cè)為參數(shù)輸入欄，右側(cè)為模型示意圖。如圖7所示為燒結(jié)環(huán)冷機自動仿真系統(tǒng)仿真界面，其中下半部分為功能區(qū)，包括建模、仿真功能按鍵，圖形顯示按鍵以及結(jié)果保存按鍵。左上部與右上部分別為結(jié)果的數(shù)值顯示框與圖像顯示框。

3 Gambit 、Fluent 軟件命令語法分析與封裝技術

3.1 Gambit 參數(shù)化建模

本系統(tǒng)的建模與網(wǎng)格劃分環(huán)節(jié)采用了 Gambit軟件，Gambit內(nèi)置了TUI命令語言，要實現(xiàn)參數(shù)化建模，只需利用 TUI 語言將模型參數(shù)與建模步驟編制成jounal文件導入Gambit即可。

3.1.1 Gambit命令語法分析

相較于其他建模軟件的命令流語言，Gambit的TUI 語言相對直觀與嚴謹，不易產(chǎn)生歧義。Gambit的建模以點作為基準，例如創(chuàng)建一個名稱為“in”的長方體的命令流語句：

“volume create "in" width 3.6 depth 2 height

0.5 offset 1.8 1 0.25 brick”

其中 3.6、2、0.5 分別為這個體的長、寬、高。而（1.8，1，0.25）是這個長方體中心點的坐標。因此，在 Gambit 中建立模型只需按照以上語法構建模型建立命令流即可。

網(wǎng)格劃分部分較為簡單，由于環(huán)冷機臺車的形狀是規(guī)整的長方體，因此只需以0.1m為尺寸劃分網(wǎng)格，網(wǎng)格劃分的命令流的語言為：

“volume mesh "in" "porous" "out" cooper

size 0.1”

指的是以 0.1 m 為尺寸，對以上“in”“porous”

“out”三個模型體劃分網(wǎng)格。

3.1.2 參數(shù)化建模及Gambit封裝技術

Gambit 本身無法替換參數(shù)的值，因此，在建模命令語言編輯完成后，需要利用 VB.NET 平臺對其進行參數(shù)化處理，以創(chuàng)建體為例，修改長方體的長度：

Dim width as string ；width=4

Dim jou1 As New IO.StreamWriter("1.jou")

jou1. WriteLine“( volume create "in" width “ +

width + ” depth 2 height 0.5 offset 1.8 1 0.25

brick”)。

這樣將此長方體的長度設置成參數(shù) width，當width的值改變時，長方體的長度也隨之改變。

待所有模型命令進行了參數(shù)化處理后，需要調(diào)用GAMBIT，并且使其執(zhí)行建立好的命令文件，調(diào)用Gambit并執(zhí)行jounal文件“1.jou”命令的語句如下：

Shell("" + gambitpath + " -inp 1.jou -new")

調(diào)用語句使用“shell”函數(shù)，此函數(shù)的使用方法為“shell+調(diào)用路徑+參數(shù)1+參數(shù)2”，在以上語句中，" + gambitpath + "指 Gambit 應用程序路徑，" - inp 1.jou"指的是執(zhí)行根目錄下文件名為1的命令文件，"-new"指建立新的模型文件。

3.2 fluent 命令語法研究與封裝技術

3.2.1 Fluent的GUI語法研究

Fluent是一款功能全面的計算流體力學仿真軟件，滿足燒結(jié)環(huán)冷機的仿真要求。Fluent 的命令語言采用 GUI 語法。因此，實現(xiàn)自動仿真需要使用GUI語言編制Fluent的jounal文件導入Fluent。

GUI語言的部分示例如下：

打開網(wǎng)格文件：“(cx-gui-do cx-set-text-entry

"select File*Text" "zhenduan.msh")”。

其中，“zhenduan.msh”指的是程序根目錄的網(wǎng)格文件，如果網(wǎng)格文件不在程序根目錄，則需要寫入完整的路徑。

設置仿真步長為 5s：“(cx-gui-do cx-set-real

entry - list "Run Calculation * Frame1 * Table1

*Frame6* Table6* RealEntry1 (Time Step Size)" '(

5))”

設置入口溫度為1023k：”(cx-gui-do cx-set-re

al - entry-list "Patch*Frame2*RealEntry1(Value)" '(

1023))”。

從以上這些例子可以看出，F(xiàn)luent 的命令語句看似繁瑣復雜，實際上并不難理解。其結(jié)構通常為：“功能名*框架+編號*表格+編號*輸入框+編號（輸入的值）”構成，只需要理解這一命令結(jié)構就能逐句編出Fluent的命令語句。

3.2.2 Fluent封裝技術

同理，調(diào)用 Fluent 并讀取 jounal 文件自動仿真的語句如下：

Shell("" + fluentpath + " 3d -g -i 2.jou")

調(diào)用Fluent同樣使用“shell”函數(shù)，以上命令中，" + fluentpath + "表示 Fluent 軟件路徑，“3d”表示采用 3d 模型，“-g”表示打開設置菜單，“-i 2.jou ”表示執(zhí)行根目錄下文件名為2的命令文件。

4 工程驗證

為了檢驗環(huán)冷機自動仿真平臺的可靠性，必須將數(shù)值仿真數(shù)據(jù)與實測數(shù)據(jù)進行對比分析。本文以鞍鋼 415㎡高效燒結(jié)環(huán)冷機為研究對象，實地檢測燒結(jié)礦溫度數(shù)據(jù)以及環(huán)冷機自動仿真平臺所需輸入?yún)?shù)。鞍鋼415 ㎡高效燒結(jié)環(huán)冷機中徑44 m，按照進風方式可以劃分為回風段與鼓風段，中間以隔板隔開，周向共布置有 5 臺風機。燒結(jié)礦入口溫度約為750 ℃，要求出口溫度為120 ℃以下。鞍鋼415 ㎡高效燒結(jié)環(huán)冷機關鍵參數(shù)與實測燒結(jié)礦物性參數(shù)如表1。

將表 1 中的參數(shù)輸入環(huán)冷機自動仿真平臺，經(jīng)過數(shù)值模擬過程，得到燒結(jié)礦溫度變化數(shù)據(jù)。415㎡高效燒結(jié)環(huán)冷機運行一周的時間為 4200 s，為了減少篇幅同時保證對比準確度，現(xiàn)選取包括進口與出口在內(nèi)的 11 個時間節(jié)點的燒結(jié)礦溫度數(shù)據(jù)進行對比，燒結(jié)礦溫度變化模擬數(shù)據(jù)與實測數(shù)據(jù)趨勢對比結(jié)果如圖8所示。

如圖 8 可知，燒結(jié)礦溫度實測曲線與模擬溫度曲線變化趨勢一致，而且重合率很高。現(xiàn)定量分析實測溫度與模擬溫度的絕對誤差與相對誤差，結(jié)果如下表2所示。

通過表 2 可以看出，最大的絕對誤差產(chǎn)生在1260 s 為 15.41℃，最大相對誤差產(chǎn)生在出口處為4.9%。由此可以得出結(jié)論，數(shù)值模擬所得到的數(shù)據(jù)誤差值控制在5%以內(nèi)，符合工藝精度要求，因此認為使用文本開發(fā)的環(huán)冷機自動仿真平臺對環(huán)冷機設備進行數(shù)值模擬所得到的結(jié)果是準確、可靠的。

5 總結(jié)

作者主要研究工作及結(jié)論如下：

（1）分析了環(huán)冷機結(jié)構與運行方式建立了環(huán)冷機的等效物理模型，探究了燒結(jié)礦的冷卻機理，選用氣體湍流模型、多孔介質(zhì)模型、局部非熱平衡方程建立了環(huán)冷機冷卻過程仿真模型。

（2）剖析了 Gambit、Fluent 軟件命令流語法，運用 VB.NET程序設計語言開發(fā)了界面友好的集基本參數(shù)輸入、自動前處理、自動仿真計算、自動后處理等功能于一體的環(huán)冷機燒結(jié)礦冷卻過程自動仿真軟件。極大簡化仿真操作過程，設計人員可以通過該軟件即時獲得燒結(jié)礦不同時刻溫度數(shù)據(jù)，也能觀察燒結(jié)礦溫度變化與排出煙氣隨時間溫度變化情況，使不具備計算流體力學知識背景的工程技術人員也能夠輕易的使用該軟件模擬仿真燒結(jié)礦冷卻過程。對指導環(huán)冷機設計具有重要意義。

（3）以鞍鋼 415㎡高效燒結(jié)環(huán)冷機為研究對象，以燒結(jié)礦溫度數(shù)據(jù)為參考指標，對比了實測數(shù)據(jù)與環(huán)冷機自動仿真平臺的模擬數(shù)據(jù)，驗證了本文開發(fā)的環(huán)冷機自動仿真平臺的準確性與易用性。

參 考 文 獻

[1] 王兆鵬,胡曉民. 燒結(jié)余熱回收發(fā)電現(xiàn)狀及發(fā)展趨勢[J]. 燒結(jié)球團,2008(01):31-35.

[2] 倪鯤鵬,劉亮,劉志永,易新建,王魯杰. 燒結(jié)礦顯熱回收發(fā)電技術及系統(tǒng)優(yōu)化[J]. 金屬材料與冶金工程,2012,40(03):18-21.

[3] 張春霞,齊淵洪,嚴定鎏,胡長慶,張旭孝. 中國煉鐵系統(tǒng)的節(jié)能與環(huán)境保護[J]. 鋼鐵,2006(11):1-5.

[4] 孫用軍,董輝,馮軍勝,張琦,王愛華,李磊. 燒結(jié)-冷卻-余熱回收系統(tǒng)熱力學分析[J]. 鋼鐵研究學報,2015,27(01):16-21.

[5] 張忠波,杜武男. 淺談多功能高效燒結(jié)環(huán)冷機設計[J]. 礦業(yè)工程, 2016,14(01):32-35.

[6] 夏建芳,喻向陽,趙先瓊 . 基于環(huán)冷機冷卻能耗最小目標的工藝參數(shù)優(yōu)化[J]. 鋼鐵研究學報,2016,28(01):13-19.

[7] Caputo A C, Cardarelli G, Pelagagge P M. Analysis of heat recovery in gas-solid moving beds using a simulation approach[J].Applied Thermal Engineering, 1996, 16(16):89-99.

[8] Pelagagge P M, Caputo A C, Cardarelli G. Comparing heat recovery schemes in solid bed cooling[J]. Applied Thermal Engineering, 1997, 17(11):1045-1054.

[9] Pelagagge P M, Caputo A C, Cardarelli G. Comparing heat recovery schemes in solid bed cooling[J]. Applied Thermal Engineering, 1997, 17(11):1045-1054.

[10] JANG Jiinyuh, CHIU Yuwei. 3-D transient conjugated heat transfer and fluid flow analysis for the cooling process of sintered bed[J]. Applied Thermal Engineering, 2009, 29(14): 2895 − 2903.

[11] Leong J C, Jin K W, Shiau J S, et al. Effect of sinter layer porosity distribution on flow and temperature fields in a sinter cooler[J]. International Journal of Minerals, 2009, 16(3): 265−272.