基于重置溫度方法的雙參數(shù)介尺度氣固傳熱模型構(gòu)建

劉怡琳，李鈺，余亞雄，黃哲慶，周強(qiáng)，2，3

（1 西安交通大學(xué)化學(xué)工程與技術(shù)學(xué)院，陜西 西安 710049； 2 新能源系統(tǒng)工程與裝備陜西省高校工程研究中心，陜西 西安710049； 3 西安交通大學(xué)動(dòng)力工程多相流國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710049）

引 言

氣固流化床反應(yīng)器由于其良好的混合、傳質(zhì)和傳熱性能，在能源化工領(lǐng)域中得到了廣泛應(yīng)用，掌握氣固兩相流系統(tǒng)的流動(dòng)和傳熱特性對(duì)于提高工業(yè)設(shè)備的性能具有重要意義。決定反應(yīng)器動(dòng)量傳遞及傳熱傳質(zhì)的關(guān)鍵參數(shù)為曳力、傳熱及傳質(zhì)系數(shù)。因此國(guó)內(nèi)外有很多研究者針對(duì)上述系數(shù)開(kāi)展了廣泛的理論、實(shí)驗(yàn)及數(shù)值模擬研究，以獲得上述系數(shù)的計(jì)算關(guān)聯(lián)式。然而，現(xiàn)有研究中一些常用的關(guān)聯(lián)式大都是針對(duì)于顆粒位置隨機(jī)均勻分布的氣固兩相流系統(tǒng)[1]。氣固兩相在非線性相間曳力和固相應(yīng)力的作用下，很容易產(chǎn)生非均勻的介尺度流動(dòng)結(jié)構(gòu)，如顆粒聚團(tuán)等[2-3]。研究表明，計(jì)算中未考慮介尺度結(jié)構(gòu)的影響會(huì)導(dǎo)致高估氣固相間曳力、傳熱和傳質(zhì)[4-6]。

由于氣固系統(tǒng)固有的復(fù)雜性，系統(tǒng)研究介尺度結(jié)構(gòu)中氣固傳熱現(xiàn)象的實(shí)驗(yàn)較少見(jiàn)諸于文獻(xiàn)[7]。在過(guò)去的20 年中，計(jì)算流體動(dòng)力學(xué)（computational fluid dynamics,CFD）已經(jīng)成為模擬氣固系統(tǒng)中流動(dòng)和傳熱的越來(lái)越強(qiáng)大的工具[8]。在模擬中為了解析介尺度結(jié)構(gòu)對(duì)流動(dòng)及傳熱傳質(zhì)的影響，方法之一是采用足夠小的計(jì)算網(wǎng)格（通常為數(shù)倍的顆粒直徑）以滿足近似均勻假設(shè)[9-10]，然而，對(duì)于大型工業(yè)反應(yīng)器而言，細(xì)網(wǎng)格模擬的計(jì)算成本是無(wú)法承受的。因此在實(shí)際應(yīng)用中，研究人員通常采用粗網(wǎng)格雙流體法（two-fluid model, TFM）或多相質(zhì)點(diǎn)網(wǎng)格方法(multi-phase particle-in-cell,MP-PIC)模擬大型流化床，但粗網(wǎng)格模擬無(wú)法解析介尺度結(jié)構(gòu)的影響，即微尺度的近似均勻模型在有介尺度結(jié)構(gòu)存在的粗網(wǎng)格模擬中不再適用，這就需要為適合工業(yè)尺度的粗網(wǎng)格模擬方法提供可以考慮亞網(wǎng)格非均勻性的介尺度模型，以得到合理的預(yù)測(cè)結(jié)果[11]。

目前，在不考慮傳熱的氣固兩相流研究中，研究者已揭示了構(gòu)建介尺度曳力、介尺度應(yīng)力等模型的必要性[12]，并基于能量最小多尺度方法（energyminimization multi-scale，EMMS）[13-15]、過(guò)濾方法[16-18]等開(kāi)發(fā)了氣固曳力介尺度模型。然而，對(duì)介尺度氣固相間傳熱的研究還比較少。Dong 等[19]采用EMMS方法研究了介尺度結(jié)構(gòu)對(duì)循環(huán)流化床提升管內(nèi)氣固傳質(zhì)的影響，他們發(fā)現(xiàn)，介尺度結(jié)構(gòu)會(huì)極大地影響傳質(zhì)。Hou 等[20]使用EMMS 方法研究了快速流化床中非均勻流動(dòng)結(jié)構(gòu)與傳遞系數(shù)之間的關(guān)系。他們發(fā)現(xiàn)，介尺度結(jié)構(gòu)對(duì)曳力系數(shù)、傳質(zhì)系數(shù)和傳熱系數(shù)有很大影響，他們指出，傳統(tǒng)的均勻模型在有非均勻流動(dòng)結(jié)構(gòu)存在時(shí)不再適用。魯波娜等[21]基于EMMS多尺度動(dòng)量傳遞及傳熱模型對(duì)多產(chǎn)異構(gòu)烷烴催化裂化反應(yīng)器進(jìn)行了模擬，研究表明，相比于傳統(tǒng)的模型，采用多尺度模型可以更為準(zhǔn)確預(yù)測(cè)反應(yīng)器內(nèi)的流動(dòng)結(jié)構(gòu)及溫度分布。Shu 等[22]考慮了顆粒聚團(tuán)結(jié)構(gòu)對(duì)傳熱的影響，對(duì)傳統(tǒng)傳熱系數(shù)模型進(jìn)行了改進(jìn)，采用改進(jìn)的氣固曳力模型及氣固傳熱系數(shù)模型對(duì)下行床反應(yīng)器進(jìn)行了模擬，研究表明采用優(yōu)化模型可以定性捕捉下行床的關(guān)鍵傳熱特性。Guo等[23]基于計(jì)算流體力學(xué)-離散單元法（computational fluid dynamics-discrete element method, CFD-DEM）研究了聚團(tuán)在氣固傳熱中的作用，他們定義了一個(gè)傳熱修正量來(lái)量化由聚團(tuán)引起的局部不均勻性的影響，并提出了一個(gè)基于假定形狀概率密度分布函數(shù)的統(tǒng)計(jì)模型。Lane 等[24-25]通過(guò)過(guò)濾雙流體模型（TFM）的計(jì)算結(jié)果，開(kāi)發(fā)了一個(gè)用于浸沒(méi)水平圓柱體的氣固兩相流傳熱的亞網(wǎng)格模型。在他們的模型中，采用固含率、固相速度、圓柱體幾何形狀（直徑和間距）和Peclet 數(shù)來(lái)封閉Nusselt 數(shù)。Rauchenzauner 等[26]提出了一個(gè)基于漂移溫度的過(guò)濾相間傳熱模型。

Agrawal 等[27]構(gòu)建了基于細(xì)網(wǎng)格TFM 模擬的過(guò)濾相間傳熱模型。他們指出，過(guò)濾后的介尺度相間傳熱系數(shù)比微觀相間傳熱系數(shù)小1～2 個(gè)數(shù)量級(jí)，他們將微觀相間傳熱系數(shù)的修正量（記為Q）建模為過(guò)濾固含率和過(guò)濾器尺寸的函數(shù)，在其模型中，過(guò)濾尺度小于等于16.7dp（dp為顆粒直徑）時(shí)不再對(duì)傳熱系數(shù)進(jìn)行修正（Q=0）。Huang 等[28]通過(guò)過(guò)濾細(xì)網(wǎng)格TFM模擬數(shù)據(jù)，將Q構(gòu)建為過(guò)濾固含率、過(guò)濾器尺寸以及兩相之間的無(wú)量綱過(guò)濾溫差的函數(shù)，在其模型中，過(guò)濾尺度小于等于8.3dp時(shí)不再對(duì)傳熱系數(shù)進(jìn)行修正（Q= 0）。Li 等[29]對(duì)Huang 等的模型進(jìn)行了修正，取消了Q在（0,1）之間的限制。Lei等[30]通過(guò)過(guò)濾CFD-DEM 數(shù)據(jù)，構(gòu)建了介尺度傳熱系數(shù)修正量與過(guò)濾固含率、過(guò)濾器尺寸及氣固相間溫差的函數(shù)關(guān)系。但Huang等[28]、Li等[29]模型中的無(wú)量綱過(guò)濾溫差需要迭代獲得，在實(shí)際的粗網(wǎng)格應(yīng)用中存在困難，Lei 等[30]的模型對(duì)過(guò)濾溫差的無(wú)量綱方法（除以0.0001 K 使得大部分的氣固溫差處于0～10 之間）缺乏普適性，同樣導(dǎo)致模型應(yīng)用困難。

另外，在現(xiàn)有的傳熱介尺度模型研究中，如何保持整個(gè)系統(tǒng)的熱平衡是一個(gè)非常復(fù)雜的問(wèn)題。為了保證氣固兩相間存在持續(xù)的相間溫差，現(xiàn)有文獻(xiàn)中使用了不同的方法。Rauchenzauer 等[26]在每個(gè)時(shí)間步讓氣相溫度保持線性增加狀態(tài)并在固相中設(shè)置熱匯。Lane等[24-25]和Lei等[30]在固相中添加了熱源項(xiàng)，熱源項(xiàng)分別以1 K/s和0.1 K/s的速率加熱氣固系統(tǒng)。他們指出，加熱速率的值是任意的，其大小不會(huì)影響模型結(jié)果。Agrawal 等[27]和Huang 等[28]在固相和氣相中分別添加了熱源和熱匯項(xiàng)，他們也指出，所構(gòu)建模型對(duì)所施加的熱源熱匯項(xiàng)的大小不敏感。以上方法可以實(shí)現(xiàn)氣固相間溫差的維持，然而如何維持氣固相間溫差以獲得更符合真實(shí)物理過(guò)程（氣固兩相自由傳熱）的介尺度傳熱模型仍是一個(gè)問(wèn)題。

針對(duì)以上問(wèn)題，本文提出了一種重置溫度的方法以維持氣固溫差，使得氣固兩相可自由換熱，并基于細(xì)網(wǎng)格CFD-DEM 數(shù)據(jù)過(guò)濾獲得了粗網(wǎng)格方便使用的介尺度氣固相間傳熱系數(shù)修正模型。

1 細(xì)網(wǎng)格CFD-DEM 傳熱計(jì)算

CFD-DEM 方法是基于歐拉-拉格朗日坐標(biāo)系的一種離散模擬方法，在歐拉坐標(biāo)系下求解氣相平均方程，在拉格朗日坐標(biāo)系下求解顆粒方程，分別追蹤單個(gè)顆粒的受力、運(yùn)動(dòng)及傳熱。有研究表明，CFD-DEM 采用小至1.75～3倍顆粒直徑（1.75dp～3dp）的網(wǎng)格系統(tǒng)可以獲得與解析到顆粒表面的直接數(shù)值模擬（particle resolved-direct numerical simulation,PR-DNS）近似的計(jì)算結(jié)果[31-32]。因此細(xì)網(wǎng)格（1.75dp～3dp）CFD-DEM 計(jì)算可以用來(lái)為粗網(wǎng)格計(jì)算提供非均勻介尺度封閉模型。本文的CFD-DEM 模擬通過(guò)對(duì)開(kāi)源軟件MFIX（版本19.2.2，https://mfix.netl.doe.gov/)進(jìn)行二次開(kāi)發(fā)實(shí)現(xiàn)，顆粒相和氣相的動(dòng)量、質(zhì)量控制方程及求解過(guò)程可參考Garg 等[33]的研究。顆粒為典型的Geldart A 類顆粒，擬二維周期邊界計(jì)算域垂直方向長(zhǎng)度Ly為960dp，橫向Lx為240dp，展向Lz為6dp，計(jì)算網(wǎng)格數(shù)量為384×96×2，對(duì)應(yīng)x,y方向的網(wǎng)格尺寸為2.5dp。為了平衡周期域上的重力，在垂直方向上設(shè)置了恒定氣體流量邊界條件。整體固含率為0.05，初始時(shí)刻顆粒在計(jì)算域中隨機(jī)均勻分布。為了保證氣固兩相間存在持續(xù)的相間溫差，本文中選取了兩種方法：方法一，參考文獻(xiàn)中的方法給氣相能量方程添加熱源項(xiàng)，該熱源項(xiàng)以0.1 K/s 的速度加熱氣體，該設(shè)定值的大小不影響計(jì)算結(jié)果[27-28]；方法二，每隔固定時(shí)間重置氣體溫度，所間隔時(shí)間為1為顆粒弛豫時(shí)間，單位為s，=18μg，其中，ρp為顆粒密度；dp為顆粒直徑；μg為氣相黏度），氣固兩相在1τstp內(nèi)自由換熱。在計(jì)算初始時(shí)，氣相溫度設(shè)置為1000 K，固相溫度設(shè)置為0 K，氣固兩相間傳熱由能量方程控制，氣相的能量方程如式(1)所示，具體可參考Lei 等[30]和Hou 等[34]的研究。

計(jì)算持續(xù)50τ，每0.2輸出一次結(jié)果，計(jì)算后全場(chǎng)平均滑移速度、平均氣固溫差等統(tǒng)計(jì)數(shù)據(jù)達(dá)到穩(wěn)定狀態(tài)，即系統(tǒng)達(dá)到統(tǒng)計(jì)穩(wěn)態(tài)[36-37]，后30τstp的數(shù)據(jù)用于模型構(gòu)建。模型通過(guò)對(duì)計(jì)算結(jié)果進(jìn)行過(guò)濾、擬合獲得，所采用的空間過(guò)濾方法與Agrawal等[27]相同，具體可參考文獻(xiàn)[9,27]，研究中所采用的過(guò)濾尺度為5dp、10dp、17.5dp、20dp、25dp和40dp。

表1 主要計(jì)算參數(shù)及設(shè)置Table 1 Parameters and settings

2 結(jié)果與討論

與文獻(xiàn)[27-29]中相同，本文針對(duì)介尺度傳熱系數(shù)修正量Q進(jìn)行建模。Q定義為：

2.1 兩種維持氣固溫差方法的討論

為了便于表述，將方法一稱為熱源項(xiàng)方法，方法二稱為重置溫度方法。圖1為兩種溫度設(shè)置方法下的固含率分布，即特定固含率對(duì)應(yīng)的網(wǎng)格數(shù)量占總計(jì)算網(wǎng)格數(shù)量的百分比?？梢钥吹剑诒居?jì)算中，固含率在（0.025,0.03）范圍內(nèi)的網(wǎng)格數(shù)量較多，存在峰值，兩種方法下固含率占比分布基本一致，流動(dòng)結(jié)構(gòu)基本一致，在此前提下開(kāi)展兩種維持氣固溫差的方法的對(duì)比和討論是有效的。

圖1 兩種方法下的固含率分布Fig.1 The distribution of solid volume fraction under two methods

圖2為兩種方法下計(jì)算域中平均氣固溫差隨時(shí)間的變化，由圖可知，兩種方法均可以使得系統(tǒng)中存在一定的氣固溫差。方法一可以使得系統(tǒng)中氣固溫差在計(jì)算達(dá)到穩(wěn)定后處于較為穩(wěn)定的狀態(tài)，波動(dòng)幅度不大；方法二則是在每次重置氣體溫度后，氣固溫差逐漸下降，直至下一次重置氣體溫度。值得注意的是，兩種方法得到的平均氣固溫差的絕對(duì)值不同，但在氣固兩相傳熱穩(wěn)定后該溫差的大小不影響Q的統(tǒng)計(jì)結(jié)果[27-28]。

圖2 兩種方法下計(jì)算域中平均氣固溫差隨時(shí)間的變化Fig.2 The average gas-solid temperature difference in the calculation domain versus time under the two methods

為了對(duì)比兩種方法對(duì)計(jì)算域內(nèi)氣固兩相傳熱的影響，本文統(tǒng)計(jì)了不同局部固含率的傳熱量占計(jì)算域總傳熱量的百分比，以及局部單位體積傳熱量（某網(wǎng)格內(nèi)的傳熱量除以網(wǎng)格體積）與總單位體積傳熱量（總傳熱量除以總體積）之比，如圖3、圖4 所示[為了凸顯兩種方法在大固含率下的區(qū)別，對(duì)縱坐標(biāo)取對(duì)數(shù)，如圖3(b)和圖4(b)所示]。由圖3 可知，兩種方法下傳熱量占比均在0.025≤?s≤0.03 時(shí)存在峰值，即兩種方法下占總網(wǎng)格數(shù)11.0%（圖1）的固含率處于0.025≤?s≤0.03 的網(wǎng)格分別貢獻(xiàn)了25.5%（重置溫度方法）和18.65%（熱源項(xiàng)方法）的傳熱量。由圖4 可知，局部單位體積傳熱量與總單位體積傳熱量之比同樣在0.025≤?s≤0.03時(shí)存在峰值，重置溫度方法中該固含率范圍內(nèi)局部單位體積傳熱量是總單位體積傳熱量的2.32 倍，熱源項(xiàng)方法中該固含率范圍內(nèi)局部單位體積傳熱量是總單位體積傳熱量的1.71 倍。重置溫度方法在固含率較小時(shí)（?s< 0.03）局部單位體積傳熱量與總單位體積傳熱量之比大于熱源項(xiàng)方法，而在固含率?s> 0.03 的范圍內(nèi)該比值小于熱源項(xiàng)方法，二者差異在固含率大于0.1 時(shí)尤為明顯。

圖3 計(jì)算域內(nèi)總傳熱量在不同局部固含率下的分布Fig.3 The distribution of total heat transfer in the calculation domain under different local solid volume fraction

圖4 局部單位體積傳熱量與總單位體積傳熱量之比Fig.4 The ratio of local gas-solid heat transfer per unit volume to the total gas-solid heat transfer per unit volume

為了分析0.025≤?s≤0.03 的特殊性，圖5 展示了某時(shí)刻的顆粒位置圖和固含率云圖，可以看到0.025≤?s≤0.03（圖5右側(cè)圖片綠色部分）的網(wǎng)格大多處于顆粒聚團(tuán)邊界位置，由此，下文將固含率范圍為0.025≤?s≤0.03、?s<0.025和?s>0.03的區(qū)域分別稱為界面、稀相和濃相區(qū)域。事實(shí)上，界面所對(duì)應(yīng)固含率范圍受整體固含率影響較大，本文整體固含率為0.05，所以界面所對(duì)應(yīng)的固含率值較小。結(jié)合圖3～圖5 可以發(fā)現(xiàn)，兩種方法均表明聚團(tuán)界面位置的局部單位體積氣固傳熱量與總單位體積傳熱量之比最大，重置溫度方法在稀相和界面位置的局部單位體積傳熱量與總單位體積傳熱量之比大于熱源項(xiàng)方法，而在濃相位置的局部單位體積傳熱量與總單位體積傳熱量之比小于熱源項(xiàng)方法。為了更直觀地了解本文所研究的系統(tǒng)中傳熱的分布情況，圖6 統(tǒng)計(jì)了稀相、界面和濃相的網(wǎng)格數(shù)量占比（除去固含率為0的網(wǎng)格）及傳熱量占比，可以看到兩種方法的固含率占比一致，而傳熱量占比卻有所區(qū)別。兩種方法的結(jié)果中，濃相、界面和稀相的網(wǎng)格占比分別為55.3%、11%和33.7%，而其對(duì)應(yīng)的傳熱量占比卻差別顯著。在重置溫度方法下，濃相、界面和稀相分別貢獻(xiàn)了36.7%、25.5%和37.8%的傳熱量；熱源項(xiàng)方法下，濃相、界面和稀相分別貢獻(xiàn)了55.4%、18.6%和26.0%的傳熱量。總地來(lái)說(shuō)，重置溫度方法下稀相和界面的傳熱量占比高于熱源項(xiàng)方法，而濃相的傳熱量占比則低于熱源項(xiàng)方法。究其原因，熱源項(xiàng)方法在氣相中添加的熱源（Sg）由式（3）計(jì)算：

圖5 某時(shí)刻顆粒位置及固含率分布Fig.5 Particle position and distribution of solid volume fraction at a certain time

圖6 稀相、界面和濃相的網(wǎng)格數(shù)量占比及傳熱量占比Fig.6 The percentage of grid number and heat transfer of dilute phase,interface and dense phase

式中，Π?為氣相的能量輸入速率，該設(shè)定值的大小不影響計(jì)算結(jié)果；?sd為計(jì)算域整體固含率；?gd為計(jì)算域整體氣含率。由式（3）可知，氣相所加熱源項(xiàng)的大小與所在網(wǎng)格的氣含率（?g）有關(guān)，不論網(wǎng)格內(nèi)真實(shí)傳熱過(guò)程如何，均有熱源加入促使其傳熱。當(dāng)體系中存在化學(xué)反應(yīng)時(shí)，熱源項(xiàng)方法可能與實(shí)際較為吻合；當(dāng)體系中不存在化學(xué)反應(yīng)時(shí)，在濃相區(qū)域中氣固相間傳熱一般會(huì)較快完成致使氣固相間溫差較小進(jìn)而傳熱量減少，但熱源項(xiàng)方法會(huì)強(qiáng)迫濃相區(qū)域保持一定的溫差進(jìn)而促進(jìn)濃相區(qū)域的氣固相間傳熱。重置溫度方法保持了氣固相間自由傳熱特性，僅通過(guò)反復(fù)重置溫度獲得足夠的統(tǒng)計(jì)數(shù)據(jù)來(lái)完成建模。在無(wú)化學(xué)反應(yīng)的氣固傳熱系統(tǒng)中，氣固兩相間由于溫差的存在會(huì)在能量方程控制下自發(fā)傳熱，該過(guò)程與本文所提出的重置溫度方法所描述的物理過(guò)程類似，因此采用重置溫度方法所構(gòu)建的介尺度傳熱模型可能會(huì)具有較高的準(zhǔn)確度。而在有化學(xué)反應(yīng)的氣固傳熱系統(tǒng)中，由于反應(yīng)和流動(dòng)、傳熱之間存在復(fù)雜的耦合關(guān)系，熱源項(xiàng)方法所構(gòu)建模型是否能較好地近似系統(tǒng)內(nèi)的傳熱過(guò)程尚需更多的模擬或?qū)嶒?yàn)來(lái)進(jìn)行分析和驗(yàn)證。

2.2 介尺度氣固傳熱模型構(gòu)建

圖8 為兩種方法在不同尺度下過(guò)濾得到的Q的平均值隨過(guò)濾固含率的變化，兩種方法下獲得的Q有所不同，在固含率較大時(shí)，二者差異更為明顯，整體上重置溫度方法的Q大于熱源項(xiàng)方法獲得的Q。

基于圖8 的數(shù)據(jù)，可將Q構(gòu)建為過(guò)濾固含率和過(guò)濾尺度的函數(shù)，重置溫度方法獲得的模型擬合效果如圖9所示，模型為：

圖7 重置溫度后不同時(shí)刻的過(guò)濾Q隨過(guò)濾固含率的變化Fig.7 The change of filtered Q with the filtered solid volume fraction at different time after resetting the temperature

圖8 兩種方法在不同尺度下獲得的Q的平均值隨過(guò)濾固含率的變化Fig.8 The average value of Q obtained by the two methods varied with the filtered solid volume fraction at different filter sizes

圖9 重置溫度方法下介尺度模型擬合Q和真實(shí)Q隨過(guò)濾固含率的變化Fig.9 The variation of fitting Q obtained by the mesoscale model and real Q with the filtered solid volume fraction under the resetting temperature method

值得注意的是，由圖9 可知，當(dāng)過(guò)濾尺度較大時(shí)，如為25dp及40dp時(shí)，過(guò)濾固含率的最大值小于小過(guò)濾尺度的最大值。造成此現(xiàn)象的原因?yàn)橛?jì)算區(qū)域有限（240dp× 960dp），過(guò)濾尺度為40dp時(shí)已占據(jù)了橫向尺寸的1/6，很難出現(xiàn)粗網(wǎng)格上的高固含率。加之全場(chǎng)整體固含率為0.05，隨著過(guò)濾尺度的增大，粗網(wǎng)格內(nèi)的固含率會(huì)逐漸接近整體固含率。

2.3 介尺度氣固傳熱模型先驗(yàn)分析

本節(jié)將對(duì)上述模型進(jìn)行先驗(yàn)分析，評(píng)價(jià)指標(biāo)采用Pearson相關(guān)系數(shù)[36]，其計(jì)算式為：

式中，x、y分別為模型預(yù)測(cè)值（Qpredict）的數(shù)據(jù)集和粗網(wǎng)格過(guò)濾值（Qexact）的數(shù)據(jù)集；xˉ、yˉ分別為數(shù)據(jù)集的平均值。采用本文所構(gòu)建的模型[式(4)]對(duì)不同過(guò)濾尺度和固含率下的Q進(jìn)行預(yù)測(cè)，預(yù)測(cè)值（Qpredict）與基于細(xì)網(wǎng)格CFD-DEM 模擬的粗網(wǎng)格過(guò)濾值（Qexact）的Pearson 相關(guān)系數(shù)如圖10 所示，本文所構(gòu)建模型的表現(xiàn)明顯優(yōu)于Agrawal 等[27]的模型，而且隨著過(guò)濾尺度的增大模型表現(xiàn)在持續(xù)變好。這是因?yàn)楸疚牡哪Ｐ褪腔诩?xì)網(wǎng)格CFD-DEM（網(wǎng)格尺寸為2.5dp）的模擬數(shù)據(jù)，而Agrawal 等[27]的細(xì)網(wǎng)格模擬采用的網(wǎng)格尺寸為16.67dp，所構(gòu)建的模型是基于至少2 倍的基礎(chǔ)網(wǎng)格即33.34dp的數(shù)據(jù)，因此理論上模型適用范圍為過(guò)濾尺度大于等于33.34dp，因此導(dǎo)致其模型在小過(guò)濾尺度時(shí)表現(xiàn)不如本文所構(gòu)建的模型。更重要的是，采用網(wǎng)格大小為16.67dp的模擬數(shù)據(jù)作為建模數(shù)據(jù)，本身已完全忽略了比16.67dp尺度更小的非均勻結(jié)構(gòu)對(duì)更大尺度上流動(dòng)和傳熱的影響，因此所構(gòu)建模型準(zhǔn)確度較低。事實(shí)上，大量研究表明，對(duì)氣固兩相流動(dòng)進(jìn)行完全解析需要在3dp左右的網(wǎng)格尺寸[38-39]。另外，與文獻(xiàn)中所采用的熱源項(xiàng)方法不同，本文所構(gòu)建的模型基于重置溫度方法，數(shù)據(jù)源的不同也是導(dǎo)致模型表現(xiàn)不同的原因之一。最后，相較于TFM，CFDDEM 可以追蹤顆粒的運(yùn)動(dòng)，計(jì)算得到的介尺度結(jié)構(gòu)更為準(zhǔn)確。圖11 為所構(gòu)建模型在不同過(guò)濾尺度下模型預(yù)測(cè)值（Qpredict）與粗網(wǎng)格過(guò)濾值（Qexact）相對(duì)誤差的概率密度分布。可以看到隨著過(guò)濾尺度的增大，相對(duì)誤差逐漸向0 集中，即在模型應(yīng)用中，在本文所研究的網(wǎng)格尺寸范圍內(nèi)，所采用的粗網(wǎng)格尺寸越大，模型的表現(xiàn)越好。本文所構(gòu)建的模型可適用于過(guò)濾尺度小于等于40dp時(shí)的粗網(wǎng)格氣固傳熱計(jì)算。另外需要指出的是，本文所構(gòu)建的模型是基于周期性邊界的模擬數(shù)據(jù)，未考慮壁面效應(yīng)，在實(shí)際應(yīng)用過(guò)程中，若壁面剪切效應(yīng)較強(qiáng)，則可能會(huì)存在誤差。

圖10 不同過(guò)濾尺度下模型預(yù)測(cè)值（Qpredict）與粗網(wǎng)格過(guò)濾值（Qexact）的Pearson相關(guān)系數(shù)Fig.10 Pearson correlation coefficients between Qpredict and Qexact under different filter sizes

圖11 不同過(guò)濾尺度下模型預(yù)測(cè)值（Qpredict）與粗網(wǎng)格過(guò)濾值（Qexact）相對(duì)誤差的概率密度分布Fig.11 The probability density distribution（PDF）of relative errors between Qpredict and Qexact at different filter sizes

3 結(jié) 論

本文采用CFD-DEM 研究了氣固兩相流相間傳熱問(wèn)題，比較了兩種維持氣固相間傳熱溫差的方法的優(yōu)缺點(diǎn)。方法一中給氣相能量方程添加了熱源項(xiàng)；方法二中每間隔一段時(shí)間進(jìn)行氣相溫度重置，重置后氣固兩相進(jìn)行自由傳熱?；谥刂脺囟确椒ǖ倪^(guò)濾數(shù)據(jù)構(gòu)建了適用于粗網(wǎng)格計(jì)算的氣固相間介尺度傳熱模型，先驗(yàn)分析表明所構(gòu)建模型具有優(yōu)越性。具體結(jié)論如下。

（1）聚團(tuán)界面位置的局部單位體積氣固傳熱量占比最大，重置溫度方法在稀相和界面位置的局部單位體積傳熱量與總單位體積傳熱量之比大于熱源項(xiàng)方法，而在濃相位置該比值小于熱源項(xiàng)方法。分析表明重置溫度方法更適合計(jì)算氣固相間自由傳熱的情況。

（2）基于過(guò)濾固含率和過(guò)濾尺度構(gòu)建了介尺度氣固相間傳熱系數(shù)修正因子的封閉模型，所構(gòu)建模型的Pearson相關(guān)系數(shù)隨過(guò)濾尺度的增大而增大，所構(gòu)建模型預(yù)測(cè)得到的Q與真實(shí)Q的Pearson 相關(guān)系數(shù)高于文獻(xiàn)中已有的雙參數(shù)模型，模型適用于粗網(wǎng)格尺度小于等于40 倍顆粒直徑的氣固兩相流傳熱計(jì)算。

由于介尺度傳熱模型的后驗(yàn)十分復(fù)雜，涉及介尺度曳力模型的選擇與驗(yàn)證，目前文獻(xiàn)中尚未有介尺度傳熱模型的后驗(yàn)報(bào)道。后續(xù)將繼續(xù)探索介尺度傳熱模型的后驗(yàn)方法，開(kāi)展相關(guān)的后驗(yàn)研究，以進(jìn)一步評(píng)估所構(gòu)建介尺度傳熱模型在實(shí)際應(yīng)用中的表現(xiàn)。