本發(fā)明涉及計算流體力學(xué)技術(shù)領(lǐng)域，更具體地，涉及一種加熱爐內(nèi)導(dǎo)熱流體仿真方法。

背景技術(shù)：

目前國內(nèi)外傳統(tǒng)CFD(Computational Fluid Dynamics，計算流體力學(xué)，簡稱CFD)研究方法，例如CFX、Fluent、Phoenics、Star-CD、CFdesign、6SigmaDC等，通常都需要前處理模塊劃分網(wǎng)格，以供后續(xù)程序進(jìn)行計算，而劃分網(wǎng)格一般耗時很長、并且整個操作過程非常復(fù)雜。

技術(shù)實現(xiàn)要素：

本發(fā)明提供一種加熱爐內(nèi)導(dǎo)熱流體仿真方法，用于克服現(xiàn)有技術(shù)的缺陷，旨在簡化操作過程、縮短仿真運算時間。

為實現(xiàn)上述目的，本發(fā)明提供一種加熱爐內(nèi)導(dǎo)熱流體仿真方法，包括：

根據(jù)待模擬的加熱爐及位于所述加熱爐內(nèi)的坯料的空間結(jié)構(gòu)，建立三維仿真模型；

根據(jù)預(yù)設(shè)的參數(shù)對所述三維仿真模型進(jìn)行LBM(Lattice Boltzmann Method，格子玻爾茲曼法，以下簡稱LBM)求解；

根據(jù)LBM的求解結(jié)果進(jìn)行LBM后處理，實現(xiàn)加熱爐內(nèi)部熱傳導(dǎo)流體及溫度的連續(xù)仿真。

所述根據(jù)待模擬的加熱爐及位于所述加熱爐內(nèi)部的坯料的空間結(jié)構(gòu)，建立三維仿真模型包括：

確定待模擬的加熱爐的平焰燒嘴的三維仿真模型；

確定待模擬的加熱爐爐體及位于所述加熱爐爐體內(nèi)的坯料的三維仿真模型。

優(yōu)選地，所述確定模擬的加熱爐的平焰燒嘴的三維仿真模型包括：

在SolidWorks環(huán)境中按照設(shè)定參數(shù)建立平焰燒嘴的三維仿真模型，將平焰燒嘴的三維仿真模型導(dǎo)入Xflow中；

在Xflow環(huán)境中進(jìn)行LBM仿真實驗；

調(diào)整平焰燒嘴的安裝參數(shù)，使得仿真實驗中火焰形狀趨于平焰燒嘴實際工作時的火焰形狀。

優(yōu)選地，所述確定模擬的加熱爐的爐體及坯料的三維仿真模型包括：

根據(jù)所述待模擬的加熱爐爐體的空間結(jié)構(gòu)以及所述平焰燒嘴的安裝參數(shù)，在SolidWorks環(huán)境中建立加熱爐爐體的三維結(jié)構(gòu)；

根據(jù)所述待模擬的坯料的空間結(jié)構(gòu)及所述坯料與所述加熱爐爐體的相對位置關(guān)系，在SolidWorks環(huán)境中建立坯料的三維結(jié)構(gòu)；

對建立的加熱爐爐體及坯料的三維結(jié)構(gòu)表面進(jìn)行簡化處理；

對建立的加熱爐爐體的三維結(jié)構(gòu)中沒有空氣流經(jīng)的部位進(jìn)行封閉處理。

優(yōu)選地，所述根據(jù)預(yù)設(shè)的參數(shù)對所述三維仿真模型進(jìn)行LBM求解包括：

根據(jù)所述三維仿真模型生成格子；

將所述三維仿真模型中平焰燒嘴三維仿真模型和坯料三維仿真模型的格子細(xì)化；

將所述坯料三維仿真模型的格子再次細(xì)化；

根據(jù)坯料的熱物理參數(shù)和設(shè)定的時間參數(shù)，進(jìn)行求解計算；

在根據(jù)定義的求解目標(biāo)值的穩(wěn)定性判斷求解結(jié)果收斂時，導(dǎo)入所述求解結(jié)果。

本發(fā)明提供的加熱爐內(nèi)導(dǎo)熱流體仿真方法，首先根據(jù)加熱爐及位于加熱爐內(nèi)的坯料的空間結(jié)構(gòu)建立三維仿真模型，然后根據(jù)坯料的屬性等設(shè)定參數(shù)，并根據(jù)預(yù)設(shè)的參數(shù)對三維仿真模型進(jìn)行LBM求解，最后再根據(jù)LBM的求解結(jié)果進(jìn)行LBM后處理，實現(xiàn)加熱爐內(nèi)部熱傳導(dǎo)流體及溫度的連續(xù)仿真；實現(xiàn)對鋼坯等坯料加熱過程的仿真，由于LBM體系在運算過程中不需要網(wǎng)格，因此節(jié)省了大量劃網(wǎng)格的時間，大大縮短了仿真時間，同時物理背景清晰，操作實現(xiàn)簡單，具有巨大的實際應(yīng)用價值。

附圖說明

圖1為本發(fā)明實施例一的提供的加熱爐內(nèi)導(dǎo)熱流體仿真方法的流程圖；

圖2為本發(fā)明實施例一中加熱爐的三維仿真模型的示意圖；

圖3為本發(fā)明實施例一中坯料及加熱爐的三維仿真模型的示意圖；

圖4為圖1中步驟S10的細(xì)化流程圖；

圖5為發(fā)明實施例一中平焰燒嘴的結(jié)構(gòu)示意圖；

圖6為發(fā)明實施例一中平焰燒嘴三維仿真模型的示意圖；

圖7為圖4中步驟S11的細(xì)化流程圖；

圖8為圖4中步驟S12的細(xì)化流程圖；

圖9為圖1中步驟S20的細(xì)化流程圖；

圖10為圖9中步驟S21中加熱爐及坯料三維仿真模型的格子生成示意圖；

圖11為圖10的俯視圖。

具體實施方式

下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例。基于本發(fā)明中的實施例，本領(lǐng)域普通技術(shù)人在沒有作出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護(hù)的范圍。

實施例

如圖1所示，本發(fā)明提供一種加熱爐內(nèi)導(dǎo)熱流體仿真方法，包括以下步驟：

步驟S10，根據(jù)待模擬的加熱爐及位于所述加熱爐內(nèi)的坯料的空間結(jié)構(gòu)，建立三維仿真模型；

首先根據(jù)待模擬的加熱爐的空間結(jié)構(gòu)建立加熱爐三維仿真模型，具體參見圖2，然后根據(jù)坯料的空間結(jié)構(gòu)以及坯料在加熱爐內(nèi)的位置關(guān)系建立坯料的三維仿真模型，最后兩者的綜合就得到了所需的三維仿真模型，坯料A、加熱爐B以及位于加熱爐B頂部的多個平焰燒嘴C的空間結(jié)構(gòu)參見圖3，建模時可以通過三維制圖軟件實現(xiàn)。

步驟S20，根據(jù)預(yù)設(shè)的參數(shù)對所述三維仿真模型進(jìn)行LBM求解；

將步驟S10中建立的三維仿真模型導(dǎo)入LMB系統(tǒng)后，還需要根據(jù)坯料的屬性設(shè)定燒嘴的空氣流量、燃?xì)饬髁康任锢韰?shù)，坯料、壁面的邊界條件的相關(guān)參數(shù)進(jìn)行設(shè)置，以及設(shè)定的總時間步長、頻率等參數(shù)，然后LMB系統(tǒng)根據(jù)上述仿真模型及設(shè)定的參數(shù)進(jìn)行計算求解。

粒子分布函數(shù)f(r,e,t)，drde表示在t時刻，在空間r處，粒子的速度在e到e+de的粒子數(shù)密度。

由統(tǒng)計力學(xué)可知，粒子分布函數(shù)f的演化滿足以下方程：

其中F表示外力，Ω為碰撞項。

根據(jù)現(xiàn)行化的BGK碰撞模型：

其中f^eq是平衡態(tài)時的分布函數(shù)，τ為無量綱弛豫時間，表征流體恢復(fù)到平衡態(tài)時的速度。

其中v為運動粘性，T為溫度。

粒子分布函數(shù)f滿足Lattice Boltzmann方程：

其中，下標(biāo)i表示給定的粒子運動方向，Δt為時間步長。

流體密度ρ和動量ρu可由分布函數(shù)求得：

壓力可由下面的等熵狀態(tài)方程直接計算獲得：

上面的公式考慮密度有微小壓縮性。

按照以下公式限制τ，使得數(shù)值具備穩(wěn)定性：

采用三維模型D3Q19模型進(jìn)行LBM求解：

由以下平衡分布函數(shù)

其中Δx為網(wǎng)格步長，Δt為時間步長；

先由初始時刻t₀的ρ₀和u₀得到f_i^eq，通過遷移和碰撞得到新的粒子分布函數(shù)f_i，在計算出新時刻t₁的ρ₁和u₁，將ρ₁和u₁作為新的初始值進(jìn)行下一個時刻的計算。

步驟S30，根據(jù)LBM的求解結(jié)果進(jìn)行LBM后處理，實現(xiàn)加熱爐內(nèi)部熱傳導(dǎo)流體及溫度的連續(xù)仿真。這里的后處理即可視化操作，也就是將仿真計算結(jié)果(以數(shù)字文檔的形式)以流線圖、溫度分布圖等圖像的方式顯示出來。

通過微元(在LBM中是格子)的計算，傳遞到整個流體域，流體的任意時間的各種物理參數(shù)通過N-S方程(Navier-Stokes equations納維-斯托克斯方程組，描述粘性流體動量守恒的運動方程，反映了粘性流體又稱真實流體流動的基本力學(xué)規(guī)律)與其他的數(shù)學(xué)方程的聯(lián)立求解都可以得到，將獲得溫度、流速等物理參數(shù)隨時間的連續(xù)變化通過圖像呈現(xiàn)出來即實現(xiàn)了加熱爐內(nèi)部熱傳導(dǎo)流體及溫度的連續(xù)仿真。

本發(fā)明提供的加熱爐內(nèi)導(dǎo)熱流體仿真方法，首先根據(jù)加熱爐及位于加熱爐內(nèi)的坯料的空間結(jié)構(gòu)建立三維仿真模型，然后根據(jù)坯料的屬性等設(shè)定參數(shù)，并根據(jù)預(yù)設(shè)的參數(shù)對三維仿真模型進(jìn)行LBM求解，最后再根據(jù)LBM的求解結(jié)果進(jìn)行LBM后處理，實現(xiàn)加熱爐內(nèi)部熱傳導(dǎo)流體及溫度的連續(xù)仿真；實現(xiàn)對鋼坯等坯料加熱過程的仿真，由于LBM體系在運算過程中不需要網(wǎng)格，因此節(jié)省了大量劃網(wǎng)格的時間，大大縮短了仿真時間，同時物理背景清晰，操作實現(xiàn)簡單，具有巨大的實際應(yīng)用價值。

本發(fā)明面向鋼鐵行業(yè)節(jié)能降耗背景下迫切需求，以鋼坯加熱過程為研究坯料，對加熱過程進(jìn)行仿真；本發(fā)明采取CFD的方法對加熱爐進(jìn)行仿真，此方法相當(dāng)于在計算機內(nèi)對步進(jìn)梁式鋼坯加熱爐進(jìn)行“虛擬”的實驗，此方法具有成本低、周期短、環(huán)境寬的優(yōu)點。在加熱爐設(shè)計中引入CFD技術(shù)，可以減少設(shè)計失誤和實驗研究所需的費用，預(yù)測爐內(nèi)氣流、溫度等分布，通過相應(yīng)技術(shù)對各項指標(biāo)進(jìn)行優(yōu)化計算，以達(dá)到高效、低排放的目標(biāo)。

參見圖4，其中步驟S10包括以下步驟：

步驟S11，確定待模擬的加熱爐的平焰燒嘴的三維仿真模型；平焰燒嘴的結(jié)構(gòu)可以通過三維軟件繪制獲得，具體參見圖5。

一臺加熱爐通過多個平焰燒嘴向加熱爐內(nèi)噴火，以對加熱爐內(nèi)的坯料進(jìn)行加熱，平焰燒嘴安裝在加熱爐的頂部，平焰燒嘴的接頭部分位于加熱爐外部用于通入燃料和空氣，平焰燒嘴的風(fēng)盤部分位于加熱爐內(nèi)用于將進(jìn)入空氣以螺旋向前推進(jìn)的方式運動，從而使風(fēng)盤前部噴嘴噴出的燃料燃燒產(chǎn)生螺旋向前的火焰，其中加熱爐的頂部內(nèi)壁上需要預(yù)留或加工用于安裝平焰燒嘴的凹槽，凹槽起到導(dǎo)流作用，與燃料及空氣的配比、壓力、流量等參數(shù)共同改變火焰的形狀，凹槽的形狀尺寸作為平焰燒嘴三維仿真模型的一部分通過LBM仿真實驗來獲得。

步驟S12，確定待模擬的加熱爐爐體及位于加熱爐爐體內(nèi)的坯料的三維仿真模型。

加熱爐爐體的結(jié)構(gòu)可以通過三維軟件繪制獲得，具體參見圖2。

位于加熱爐爐體內(nèi)的坯料的結(jié)構(gòu)可以通過三維軟件繪制獲得，參見圖3。

參見圖7，其中步驟S11包括以下步驟：

步驟S111，在SolidWorks環(huán)境中按照設(shè)定參數(shù)建立平焰燒嘴的三維仿真模型，將平焰燒嘴的三維仿真模型導(dǎo)入Xflow中；

這里設(shè)定參數(shù)指的是按照待模仿的平焰燒嘴的實際結(jié)構(gòu)和尺寸設(shè)定的參數(shù)，在SolidWorks繪圖軟件中繪制平焰燒嘴的三維結(jié)構(gòu)，參見圖5。當(dāng)然也本發(fā)明的其他實施例中也可以在AuctoCAD等其他繪圖軟件中繪制。

步驟S112，在Xflow環(huán)境中進(jìn)行LBM仿真實驗；

步驟S113，調(diào)整平焰燒嘴的安裝參數(shù)，使得仿真實驗中火焰形狀趨于平焰燒嘴實際工作時的火焰形狀。這里的安裝參數(shù)指的是將平焰燒嘴安裝在加熱爐爐體內(nèi)頂壁上時，位于加熱爐爐體內(nèi)頂壁的凹槽的結(jié)構(gòu)參數(shù)。

在LBM仿真實驗過程中觀察平焰燒嘴的火焰形狀，并與實際的步進(jìn)梁式鋼坯加熱爐平焰燒嘴工作時的火焰形狀作對比，并調(diào)整平焰燒嘴的安裝參數(shù)，反復(fù)對比調(diào)整，直到LBM仿真實驗中的火焰形狀與實際的步進(jìn)梁式鋼坯加熱爐平焰燒嘴工作時的火焰形狀相同，確定步進(jìn)梁式鋼坯加熱爐平焰燒嘴的精確模型，參見圖6。

參見圖8，其中步驟S12包括以下步驟：

步驟S121，根據(jù)待模擬的加熱爐爐體的空間結(jié)構(gòu)以及平焰燒嘴的安裝參數(shù)，在SolidWorks環(huán)境中建立加熱爐爐體的三維結(jié)構(gòu)；使用SolidWorks軟件CAD模塊，根據(jù)所需對步進(jìn)梁式鋼坯加熱爐爐體建立三維結(jié)構(gòu)。

步驟S122，根據(jù)待模擬的坯料的空間結(jié)構(gòu)及坯料與加熱爐爐體的相對位置關(guān)系，在SolidWorks環(huán)境中建立坯料的三維結(jié)構(gòu)；在上述步進(jìn)梁式鋼坯加熱爐爐體三維結(jié)構(gòu)的基礎(chǔ)上，使用SolidWorks軟件CAD模塊，根據(jù)對步進(jìn)梁式鋼坯加熱爐與坯料的位置關(guān)系，建立坯料的三維結(jié)構(gòu)，進(jìn)而形成加熱爐爐體及坯料的三維結(jié)構(gòu)。

步驟S123，對建立的加熱爐爐體及坯料的三維結(jié)構(gòu)表面進(jìn)行簡化處理；具體包括對建立的加熱爐爐體及坯料的三維結(jié)構(gòu)表面的小圓角或小孔(包括工藝孔或吊裝孔等)進(jìn)行去除處理，并對建立的加熱爐爐體及坯料的三維結(jié)構(gòu)表面的小尖角進(jìn)行平滑處理；

步驟S124，對建立的加熱爐爐體的三維結(jié)構(gòu)中沒有空氣流經(jīng)的部位進(jìn)行封閉處理。這里指的是對加熱爐爐體的三維結(jié)構(gòu)中理論上沒有空氣流經(jīng)的部位表面進(jìn)行封閉，減沒有必要的計算，提高計算效率。

參見圖9，其中步驟S20包括以下步驟：

步驟S21，根據(jù)所述三維仿真模型生成格子；這里的三維仿真模型包括加熱爐三維仿真模型和坯料三維仿真模型，其中加熱爐三維仿真模型包括加熱爐爐體三維仿真模型和平焰燒嘴三維仿真模型；在Xflow軟件中按照上述三維仿真模型的結(jié)構(gòu)生成格子。參見圖10、圖11。

步驟S22，將所述三維仿真模型中平焰燒嘴三維仿真模型和坯料三維仿真模型的格子細(xì)化；格子越細(xì)，計算的精度越高。

步驟S23，將所述坯料三維仿真模型的格子再次細(xì)化；根據(jù)需求，有選擇性地對局部三維仿真模型進(jìn)行區(qū)分，能夠在仿真度和計算效率之間找到平衡，保證高度仿真效果的同時，不失效率。

步驟S24，根據(jù)坯料的熱物理參數(shù)和設(shè)定的時間參數(shù)，進(jìn)行求解計算；根據(jù)坯料的材料屬性，設(shè)置坯料(實施例中選用鋼坯)的熱物理參數(shù)(包括密度、熱傳導(dǎo)系數(shù)、比熱容等)，平焰燒嘴的空氣流量、燃?xì)饬髁?，鋼坯、壁面的邊界條件的相關(guān)參數(shù)進(jìn)行設(shè)置，設(shè)定的總時間步長，開始計算。

步驟S25，在根據(jù)定義的求解目標(biāo)值的穩(wěn)定性判斷求解結(jié)果收斂時，導(dǎo)入所述求解結(jié)果。

由于LBM求解的Lattice Boltzmann方程屬于建立離散方程組的方法，只能用數(shù)學(xué)方法求得近似解，而近似解與精確解之間會產(chǎn)生誤差，若誤差通過計算不會擴大，并控制在一定的范圍內(nèi)，則稱求解結(jié)果收斂，反之則為發(fā)散。這個控制的范圍也就是求解目標(biāo)值的穩(wěn)定性判斷標(biāo)準(zhǔn)，如果求解結(jié)果收斂，則認(rèn)為求解結(jié)果可行，并將其導(dǎo)入到后處理程序，通過后處理觀察步進(jìn)梁式鋼坯加熱爐內(nèi)部的空氣流動。

最后應(yīng)說明的是：以上實施例僅用以說明本發(fā)明的技術(shù)方案，而非對其限制；盡管參照前述實施例對本發(fā)明進(jìn)行了詳細(xì)的說明，本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解：其依然可以對前述各實施例所記載的技術(shù)方案進(jìn)行修改，或者對其中部分技術(shù)特征進(jìn)行等同替換；而這些修改或者替換，并不使相應(yīng)技術(shù)方案的本質(zhì)脫離本發(fā)明各實施例技術(shù)方案的精神和范圍。