REV尺度多孔介質格子Boltzmann方法的數學模型及應用的研究進展.pdf

CHEMICAL INDUSTRY AND ENGINEERING PROGRESS 2016 年第 35 卷第 6 期 1698 化 工 進 展 REV 尺度多孔介質格子 Boltzmann 方法的數學模型及應用的 研究進展 張瀟丹 1，2，雍玉梅2，李文軍3，趙元生4，李媛媛2，楊巧文1，楊超2 （1中國礦業(yè)大學（北京）化學與環(huán)境工程學院，北京 100083；2中國科學院過程工程研究所綠色過程與工程重 點實驗室，北京 100190；3華北科技學院環(huán)境工程學院，河北 廊坊 065201；4中國石油化工研究院渣油加氫實 驗室，北京 102200） 摘要：綜述了多孔介質表征體元尺度（REV）格子 Boltzmann 模型的研究進展，根據對多孔介質處理方式主要分 為部分反彈模型和阻力模型兩類，分析歸納了各類模型的優(yōu)缺點。由于阻力模型中滲流的廣義格子 Boltzmann 方程（GLBE）的作用力是基于 GUO 等的作用力模型，可以準確得到宏觀方程，不存在離散誤差，且模型的平 衡分布函數和作用力項中都包含反應介質特性的孔隙率，因而應用最為廣泛。本文還重點介紹了 REV 尺度多孔 介質 LBE 模型在流動、傳熱、傳質、化學反應及相變等過程中的具體應用，認為 REV 尺度多孔介質內的三傳一 反數學模型中需要加入孔隙尺度因素，在更大工程尺度上應該考慮過程參數的各向異性，展望了 REV 尺度多孔 介質 LBE 模型的發(fā)展和應用前景。 關鍵詞：多孔介質；表征體元尺度；格子 Boltzmann 方法；流動；傳熱；傳質 中圖分類號：TQ021.9 文獻標志碼：A 文章編號：10006613（2016）06169815 DOI：10.16085/j.issn.1000-6613.2016.06.010 Models and application of lattice Boltzmann method at REV-scale in porous media ZHANG Xiaodan1 ，2， YONG Yumei2， LI Wenjun3， ZHAO Yuansheng4， LI Yuanyuan2， YANG Qiaowen1， YANG Chao2 （1School of Chemical & Environmental Engineering，China University of Mining & Technology (Beijing)，Beijing 100083， China； 2Key Laboratory of Green Process and Engineering， Institute of Process Engineering， Chinese Academy of Sciences，Beijing 100190，China；3School of Environmental Engineering，North China Institute of Science and Technology，Langfang 065201，Hebei，China；4Laboratory of Residue Hydrotreating，Research Institute of Petroleum Processing，PetroChina，Beijing 102200，China） Abstract： This paper discusses the lattice Boltzmann model at representative elementary volume (REV) scale for porous media. According to different treatments of porous media， the lattice Boltzmann model at REV-scale for porous media can be classified into two categories，the partially bouncing-back model and the resistance model. The advantages and disadvantages of various models are analyzed. The Generalized lattice Boltzmann equation (GLBE model) in the resistance model is most widely used. Firstly， the force item of the GLBM model is based on the method proposed by Guo et al， which can be 第一作者：張瀟丹（1989） ，女，碩士研究生，主要從事化學工程 數值模擬。聯(lián)系人：雍玉梅，研究員，主要從事化學工程數值模擬 E-mail ymyong 。 特約評述 收稿日期：2016-01-19；修改稿日期：2016-01-31。 基金項目：國家自然科學基金（21276256,21490584） 、國家重點基礎研究發(fā) 展計劃（2012CB224806）及中央高校基本科研業(yè)務費（3142013097）項目。 第 6 期 張瀟丹等：REV 尺度多孔介質格子 Boltzmann 方法的數學模型及應用的研究進展 1699 accurately recovered to the macroscopic equation without discretization error. Secondly，the equilibrium distribution function and the force items involve porosity， which reflects the characteristics of porous media. This review introduces the application of the lattice Boltzmann model at REV-scale in porous media with flow， heat transfer， mass transfer， chemical reaction and phase transition. We should take the pore scale factors into consideration when studying momentum transport， heat transport， mass transport and reaction engineering in porous media at REV-scale. The anisotropic nature of process parameters should be considered at the larger engineering scale. It gives some predictions and perspectives of applications for the lattice Boltzmann model at REV-scale. Key words：porous media；representative elementary volume (REV)；lattice Boltzmann method (LBM)；flow；heat transfer；mass transfer 流體在多孔介質內的流動是許多領域都涉及的 一類復雜流動，在能源化工、生物醫(yī)藥、空氣凈化、 污水處理等領域都有廣泛應用。多孔介質內的流動 是一種典型的多尺度問題，通常涉及 3 個尺度，即 孔隙尺度、表征體元尺度（REV）和宏觀尺度。格 子 Boltzmann 方法（LBM）是介于宏觀和微觀兩種 方法之間的介觀數值模擬方法，它既克服了宏觀方 法需要處理復雜邊界的劣勢，又克服了微觀方法計 算尺度受限的局限性，在模擬多孔介質方面具有獨 特優(yōu)勢1。 由于 LBM 算法簡單， 計算網格容易生成， 可以比較方便處理流體與邊界之間、不同流體組分 之間等復雜的相互作用， 且程序代碼簡單容易并行， 應用逐漸廣泛。格子 Boltzmann 方法模擬多孔介質 滲流主要有孔隙尺度和 REV 尺度?？紫冻叨?LBM 的研究對象是單個或者若干個孔隙內的流體，把介 質的骨架作為流場的邊界，用來探索滲流機理和基 本規(guī)律。REV 尺度 LBM 是以一個控制體為研究對 象， 求解 REV 體積平均的宏觀流動量， 流體與介質 的相互作用通過適當的處理方法來描述。 REV 尺度 LBM 方法相對于孔隙尺度 LBM 方法，優(yōu)點是前者 不需要介質的細節(jié)結構、僅僅依賴于介質的統(tǒng)計參 數，因而計算效率較高，比較適合大區(qū)域的工程滲 流計算。已有文獻總結了 LBM 方法的基本要素2， 介紹了 LBM 的起源、連續(xù) Boltzmann 方程、邊界 條件和數值穩(wěn)定性，討論了 LBM 在流體流動、多 相流和多組分流、反應擴散中的應用。CHEN 等3 則對 LBM 在湍流、多孔介質流、燃燒、多孔介質 內燃燒的應用進行了綜述， 李楠等4分析了 LBM 在 模擬這些復雜流體問題時存在的優(yōu)缺點。SUKOP 等5通過比較兩種 REV 尺度 LBE 模型， 提出 Dardis 和McCloskey的LBE模型不能用于各向異性的多孔 介質的觀點，但沒有詳細論述兩模型的具體差異。 目前還沒有 REV 尺度多孔介質 LBM 方法的綜述， 因而本文總結了 REV 尺度多孔介質 LBE 各類數學 模型，分析了各模型的優(yōu)缺點，綜述了該方法在 REV 尺度多孔介質內流動、傳熱、傳質、化學反應 和相變等領域的應用，分析了未來應用前景，并對 其發(fā)展趨勢提出建議。 1 REV 尺度 LBE 模型 REV 尺度 LBE 模型根據對多孔介質的處理方 式不同主要分為兩類：部分反彈模型6-11，此種 模型對經典的固壁回彈模型進行修改來擬合達西公 式；阻力模型12-21，此模型通過施加外力公式， 得到 Stokes 或 Brinkman 方程。 1.1 部分反彈模型 部分反彈模型由 3 個步構成，如式(1)。前兩個 步是碰撞步和流動步， 與標準 LBE 模型相同。第三 個步是部分反彈步，該過程中流體分布函數根據該 點的固體密度重新定義。部分反彈模型中引入了多 孔介質連續(xù)表示方法，并將反彈邊界的作用耦合到 演化方程中7。 (eq) 1 (,)( , )( , )( , )( , ) ii ttiiii fctftftftt += + xxxxx (1) 式中，i 是速度離散方向；t 是離散時間；f 為 分布函數； 為松弛時間；i (x,t)是流體粒子與固 體骨架之間的碰撞，表示多孔介質對流體影響。 最早的部分反彈模型是 GAO 和 SHARMA6提 出的多孔介質格子氣自動模型，可延伸到 LBM 中 （GS-LBE 模型)。該模型能夠模擬非均質多孔介質 內的流動。該模型中 i (x,t)如式(2)。 s ( , )=( , )( , ) ii i tnftft xxx (2) 式中，ns是固體格點的平均散射密度，取值范 圍為 01，代表多孔介質的固體結構；i 是 i 的相 化 工 進 展 2016 年第 35 卷 1700 反方向。 但是 GS-LBE 模型存在格子氣的統(tǒng)計噪聲。 1998 年， DARDIS 和 MCCLOSKEY7則提出了第一 個部分反彈 LBE 模型（DM-LBE 模型)，是對 Brinkman 方程的間接求解。與 GS-LBE 模型相比， DM-LBE模型用ns代替孔隙尺度多孔介質的結構參 數布爾變量，消除了格子氣的不穩(wěn)定性。滲透率依 賴于固體格點的平均散射密度 ns，并且對應的有效 黏性系數與流體的黏性系數相同。 該模型中為 i (x,t) 如式(3)。 s ( , )(, )( , ) iiti i tnftft=+ xxcx (3) THORNE 和 SUKOP8根據部分反彈形式不同 提出了第二類部分反彈格式（TS-LBE 模型） ，該模 型流出分布函數反映臨近點流入分布函數碰撞過 程，該類反彈格式需計算每個格點的碰撞步來估算 流出分布函數密度。該模型中 i (x,t)如式(4)。 c*c* s ( , )(,)( ,) ii ti i tnf tft =+ xxcx (4) 式中 t*是碰撞前的時間，滿足 tt*t+1，fic 為碰撞后的分布函數。 DM-LBE模型和TS-LBE模型的部分反彈格式 需要根據臨近格點計算得到，計算復雜，不利于并 行計算。WALSH、BURWINKLE 和 SAAR9引入 了新的部分反彈方法（WBS-LBE 模型） ，該方法 在碰撞步前改變流入流體方向， 同時計算碰撞步和 反彈步，且計算中不包含臨近格點的流體分布函 數，使計算更容易并行，從而節(jié)省計算時間。與 DM-LBE模型和TS-LBE模型的部分反彈方法相比 較， 該方法在非均質多孔介質流動過程中滿足質量 守恒，且提高了 LB 方法數值模擬浮力驅動及應力 和濃度梯度引發(fā)的擴散過程的精度。該模型中 i (x,t)如式(5)。 cc* s ( , )( , )( ,) ii i tnftft = xxx (5) WBS-LBE 模型雖然在非均質多孔介質內滿足 質量守恒，更加適用于浮力驅動的流動和擴散過 程， 但是該模型在不同滲透特性的多孔介質界面處 速度不連續(xù)，ns接近 1 時會導致數值不穩(wěn)定。 YEHYA、NAJI 和 SUKOP10修正了該模型，提出 新的模型（YNS-LBE 模型）。該模型也在碰撞前 改變流入流體方向，但是公式中使用 i 方向的分布 函數代替i 方向的分布函數， 使其在界面處速度連 續(xù)，如式(6)。 * s ( , )( , )( ,) c iii tnftft = xxx (6) WBS-LBE 模型在 i 方向上的流體分布函數是 碰撞后的分布函數，而i 方向上的分布函數則是在 碰撞前的分布函數，導致反彈回的流體不僅沿其反 方向流體，也會向其他方向流動，引起數值不穩(wěn)定 性。 ZHU 和 MA11重新分配同一時刻同一格點碰撞 后方向相反的一對流體質點，每個方向上碰撞后的 流體質點分為兩部分，一部分沿該方向流向下一格 點，另一部分由于固體壁面的作用沿相反方向反彈 回去，提出新的部分反彈模型（ZM-LBE 模型） ，該 模型在同一時間的格點僅改變方向相反的一對流體 質點，因此該模型滿足質量守恒，且數值穩(wěn)定性比 WBS-LBE 模型好。該模型的松弛時間接近 1 時， 有效黏度系數與流體黏度系數相等。 該模型中 i (x,t) 如式(7)。 cc* s ( , )( , )( ,) ii i tnftft = xxx (7) 部分反彈模型區(qū)別主要在于如何處理反彈邊 界的作用、有效黏度、平均散射密度及有效滲透 率，總結目前主要部分反彈模型見表 1。表中包 括 DM-LBE 模型、 TS-LBE 模型、 WBS-LBE 模型、 YNS-LBE 模型和 ZM-LBE 模型， 這 5 個模型的滲 透率根據每個格點的平均散射密度 ns計算得到。 前 4 個模型在 ns=0 時，均可簡化為標準 LBE 模 型. 在 ns=0 時，ZM-LBE 模型不能自動滿足反彈 邊界，因為該模型在同一時間格點僅改變方向相 反的一對流體質點。WBS-LBE 模型、YNS-LBE 模型和 ZM-LBE 模型在 ns=1 時能夠滿足全反彈 LBE 模型，其余模型不滿足全反彈 LBE 模型，因 為當 ns=1 時，DM-LBE 部分反彈模型在不透水邊 界處有質量損失，而 TS-LBE 部分反彈格式模型 僅反映一個時間步的流體函數而不是原則上的兩 個時間步。其次，GS-LBE 模型、WBS-LBE 模型、 YNS-LBE 模型和 ZM-LBE 模型能夠滿足質量守 恒，DM-LBE 模型和 TS-LBE 模型在高滲透性和 低滲透性多孔介質邊界之間不對稱而不滿足質量 守恒，需要在格點的邊緣設置 ns來修正模型，只 有在密度恒定或者平均散射密度不變時才滿足質 量守恒。CHEN 等 22比較了 GS-LBE 模型、 DM-LBE 模型和 TS-LBE 模型，分析了在松弛時 間 =1 時， 為了保證流體黏度為正值， 3 個模型中 ns的上限值。GS-LBE 模型中 ns (1/23)0.289 DM-LBE模型中 ns (571/28)0.234，TS-LBE 第 6 期 張瀟丹等：REV 尺度多孔介質格子 Boltzmann 方法的數學模型及應用的研究進展 1701 表 1 REV 尺度多孔介質內部分反彈 LBE 模型（f為孔隙流體黏度）11 模型 ( , ) i tx 有效黏度 ns（=1） 有效滲透率 GS-LBM s ( , )( , ) i i nftft xx () ss s 121(1)(21) 6 12 nn n + 0.289 f s 2n DM-LBM s (, )( , ) i ti i nftft+ xcx () ssss s 121(1)(21)(12) 6 12 nnnn n + + 0.234 f s 2n TS-LBM c*c* s (,)( ,) ti i nftft i + xcx ()()() () () s s s 212611 1 6 121 n n n + + 1 f s 2n WBS-LBM c* s ( , )( ,) i i nftft xx ()() ()() ss ss 21121 6 1211 nn nn + + 1 () sf s 1 2 n n YNS-LBM c* s ( , )( ,) ii nftft xx 1 - ZM-LBM cc* s ( , )( ,) ii nftft xx ()()() ()() 2 ss ss 268112121 6 12112 nn nn + + 0.5 () sf s 12 2 n n 模型和WBS-LBE模型中 ns在01之間。ZM-LBE 模型的 ns上限為0.5， 說明ZM-LBE模型最多有一 半的流體質點可以反彈回相反的方向。綜上， YNS-LBE模型在 ns=0時，均可簡化為標準LBE 模型，在 ns=1時能夠滿足全反彈LBE模型，能夠 滿足質量守恒，數值穩(wěn)定性也較好，因此應用比 較好。 1.2 阻力模型 REV尺度多孔介質的另一種LBE模型是阻力 模型12-21，即多孔介質的影響通過施加外力得到 Stokes或Brinkman方 程 ， 這 個 思 想 是 基 于 NITHIARASU等23提出的求解不同孔隙度滲流的 Navier-Stokes方程，用其他CFD方法（有限差分、 有限體積）也可以方便求解，采用LBM方法是眾 多CFD方法的一種。該類模型的實現(xiàn)途徑有兩類： 第一類是修正速度12-17，即將多孔介質的阻礙作用 通過修正速度來體現(xiàn)；另一種是流體與多孔介質構 造體-固體的相互作用力直接加入分布函數的演化 方程中18-21。 1.2.1 修正速度的阻力模型 修正速度的REV尺度多孔介質LBE阻力模型 是在碰撞步中修正瞬時速度來描述多孔介質對流體 的影響，所以其演化方程與標準BGK演化方程相 同12，如式(8)。 (eq) 1 (,)( , )( , )( , ) iittii ftftftft += xcxxx (8) 1997年SPAID和PHELAN12建立了直接求解 BRINKMAN方程的LBE模型（SP-LBE模型） 。該 模型將流體-固體的相互作用力的影響通過修正平 衡分布函數的速度來實現(xiàn)，這種作用力處理方式對 應的宏觀流動方程中存在一些誤差，忽略誤差后， SP-LBE對應的宏觀方程在穩(wěn)態(tài)時就是Brinkman方 程。SP-LBE模型平衡態(tài)分布函數的計算采用修正 后的平衡態(tài)速度，修正的平衡態(tài)速度包含了多孔介 質的阻力作用，而阻力大小又與當地速度相關，如 式(9)式(11)。 () 2 (eq) (eq)(eq) 2 (eq) 242 sss () ( , )1 22 i i ii ft ccc =+ cu cuu x (9) (eq) F =+uu (10) = Fu (11) 式中， v/K，是滲透勢能參數，是黏度 v 與 滲透率 K 的比值。 與SP-LBE模型類同，F(xiàn)REED13把多孔介質阻 力通過碰撞過程在平均速度中體現(xiàn)Freed-LBE模 型，式(13)，可以用于Mach數較大的流動過程。 與多孔介質流動控制方程相比，F(xiàn)reed-LBE模型對 應的宏觀方程存在人工多余項，需要通過具體問題 來分析。該模型比較適合用于具有裂隙及導電介質 的巖石孔隙中流體及其他工程模擬。Freed-LBE模 型公式如式(12)式(14)。 ()2 2 (eq) 242 sss ( , )1 22 i i ii ft ccc =+ cucuu x (12) 11 = 1) + 22 F =+( uuuu (13) 化 工 進 展 2016 年第 35 卷 1702 Gu =u (14) () =Fuu (15) 式中，u為碰撞后速度， u 碰撞前速度，u 為 平均速度，G是速度張量， 取決于流體的運動黏度、 松弛因子和滲透系數。 Shan-Chen模型14-15最初是應用于多相流和多 組分流動，把體積力 F 在平衡速度中體現(xiàn)，體積力 F 主要包括液液流體間相互作用力、液固相間相互 作用力和外力。LI等16-17基于Shan-Chen模型提出 單相單組分LBM方法（Li-LBE模型） ，其平衡分 布函數可以簡化成式(16)，平衡態(tài)速度為式(17)。該 模型與其他模型相比，不用計算體積力 F，節(jié)省計 算時間，提高計算效率。 () () (eq) eq 2 s ,1 i ii ft c =+ cu x (16) 8 (eq)0 2(12 ) ii i f = =+ c uu (17) 8 1 0.5 1 2 i it i t f K = + = + cG u (18) 式中， K 為滲透率， 為孔隙率， G 為外部體積力。 1.2.2 修正演化方程的阻力模型 修正演化方程的阻力模型是將多孔介質效應直 接加入到演化方程，其演化方程19為式(19)。 (eq) 1 (,)( , )( , )( , ) iittiitii ftftftftF += + xcxxx (19) 式中，F(xiàn)i是 F 的函數，根據其處理方式不同， 可以分成不同的阻力模型。 MARTYS18對SP-LBE模 型 進 行 改 進 （Martys-LBE模型） ， 采用LUO24的作用力處理方 法描述介質阻力，在演化方程中直接增加了一個作 用力項來消除部分離散誤差， 但是不能完全消除， 如式(20)。其平衡態(tài)分布函數與SP-LBE模型相同， 即式(9)，但是作用力不修正平衡態(tài)分布函數的速 度，而是直接修正演化方程。Martys-LBE模型比 SP-LBE模型的數值穩(wěn)定性好，因為SP-LBE模型 修正平衡分布函數的速度，增加了非物理項。 Martys-LBE模型其作用力表達式如式(20)此時 c =Fa， i i f=。 22 ss () iii ii FF cc =+ cucu c (20) GUO等19在2002年提出了求解包含線性阻 力、 黏性項和非線性阻力項的通用滲流模型的REV 尺度滲流的廣義LBE模型（GLBE模型） 。該模型 的平衡分布函數和作用力項中均中包含反映介質特 性的孔隙率 ，如式(21)和式(22)，而且對滲流速度 大小沒有要求。此模型采用GUO等25提出的作用 力模型， 不存在離散誤差， 可得到準確的宏觀方程。 其中，F(xiàn)是包含介質阻力和外部體積力的總力，式 (23)右端第一和第二項是線性（Darcy）和非線性 （Forchheimer）介質阻力。F是結構參數，與孔隙 率有關，對由固體顆粒構成的多孔介質，如 公式(24)。 ()2 2 (eq) 242 sss ( ) ( , )1 22 i i ii ft ccc =+ cucu u x (21) 2 s 24 ss :()1 (1) 2 ii i ii c F cc =+ cFuFc cI (22) F KK +F =uu uG (23) 3 1.75 150 F = (24) 目前還有一種趨勢，將LBM與其他數值方法 結合起來模擬多孔介質內流動和傳遞過程。 ZARGHAMI等20把GLBE模型和有限體積法相結 合，用有限體積方法離散對流擴散項，如式(25)， 并用校正因子來提高數值穩(wěn)定性。該模型通過 Poiseuille流、Couette流和平板驅動流的驗證，用 于模擬均質和非均質多孔介質管道流。 () eq neqneqneq ,1,1 , neqneqneq 1,11,1 neqneqneq 1,11,11,1 1 d 11 48 1 16 ii acbd iii x yxyx y x y iii xyx yxy iii xyxyxy ffA fff A fff fff + + + = + + + (25) 式中，Ax.y為 abcd 的面積，ac、ab、bd、bc 分 別為格子邊界。 ZHANG等21提出了Navier-Stokes方程的顯性格 式解的LBE模型（Zhang-LBE模型） ，該模型為通過 Chapman-Enskog和泰勒展開，省略高階項之后，得 到平均體積Navier-Stokes方程，在演化方程式(26) 中加入了附加項 Pi，如式(28)。該模型只在平衡分布 第6期 張瀟丹等：REV尺度多孔介質格子Boltzmann方法的數學模型及應用的研究進展 1703 函數中加入反應介質特性的孔隙率，如式(27)，但是 與GLBE模型中孔隙率加入方式不同， 作用力項處理 方法采用式(29)，F(xiàn)公式與GLBE模型一致， 見式(23)。 該模型能有效地解決多孔介質流和兩相流。 (eq) 1 (,)( , )( , )( , ) ii ttiiit iti f xctfx tfx tfx tPF += + (26) ()2 2 (eq) 242 sss ( ) ( , )1 22 i i ii u ft ccc =+ cuc u x (27) 24 ss 1 (1) 2 ii iii Pp cc =+ cucu c (28) 24 ss 1 (1) 2 ii iii F cc =+ cucu c F (29) 兩種類型多孔介質REV尺度LBE模型，應 用較為廣泛的是Freed-LBE模型和GLBE模型。 如果考慮不可壓流時，GLBE模型更加精確；對 于穩(wěn)態(tài)不可壓低速多孔介質流動問題，兩者的區(qū) 別在于GLBE模型加入了非線性阻力項，在低速 情況下，該項的作用較弱26。GLBE模型的作用 力是基于GUO等的作用力模型， 可以準確得到宏 觀方程，不存在離散誤差；其次，該模型的平衡 分布函數中包含反應介質特性的孔隙率，而其余 REV尺度LBE模型的平衡分布函數與標準模型 類似，不能反映介質特性；最后，GLBE模型可 以用于高速流動。 2 REV 尺度 LBE 模型應用現(xiàn)狀 2.1 REV尺度LBE模型應用于多孔介質內流動 各REV尺度LBE模型在多孔介質內流動中均 有運用，但是運用最廣泛的是GLBE模型。 WBS-LBE模型、YNS-LBE模型和ZM-LBE模 型均被用于模擬不同滲透特性多孔介質內的流 動10-11，如圖1。模擬結果顯示WBS-LBE模型在 界面處不連續(xù)，YNS-LBE模型和ZM-LBE模型在 界面處速度連續(xù)，而且在 ns較小時的連續(xù)性比較好 （ns1=ns3=0.2，ns2=0.1） ，且多孔介質的孔隙度變化 不會影響中心處的速度。 REV尺度LBE模型不僅被用于分層多孔介質 內速度的連續(xù)性， 還被用于多孔介質內Klinkenberg 效應等基礎研究。孔隙大小在納米到微米范圍內， 克努森數 （Kn 數， 分子平均自由程和流動的宏觀特 征長度之比）相對較大時，氣體不再滿足連續(xù)介質 假設，導致氣體滑移流、瞬變流，甚至真空分子流。 圖1 不同滲透特性多孔介質示意圖及其速度分布圖10-11 由于存在滑移效應，測得表觀滲透率高于固有滲透 率，兩者之差隨著 Kn 數增大而增大。CHEN等27 基于四參數隨機生成方法（QSGS）并運用基本元 素建模模型（EBB模型）重新構造頁巖，采用 BESKOK和KARNIADAKIS28提出的表觀滲透率 公式，在 Kn 數為常數時根據固有滲透率計算表觀 滲 透 率 ， 運 用GLBE模 型 研 究 了 二 維 含 有 Klinkenberg效應的多孔介質內流動。高孔隙率 （=0.8）時，Klinkenberg效應可以忽略，隨著孔隙 率減小，Klinkenberg效應在低滲透多孔介質比較明 顯， 并且隨著 Kn 數增大、Klinkenberg效應越明顯， 化 工 進 展 2016年第35卷 1704 隨著壓力的減小、Klinkenberg效應更加明顯 （圖2） 。 隨后CHEN等29通過Darcy定律預測REV尺度的 整個計算區(qū)域內的有效滲透率， 研究多孔介質成分、 顆粒大小、 壓力、 滑移及吸附作用對REV尺度滲透 率的影響。 在致密儲氣巖層或頁巖層開采天然氣時， 如果不考慮Klinkenberg效應， 用表觀滲透率代替固 有滲透率，不僅會導致定量錯誤而且會導致定性 錯誤。 圖2 含有Klinkenberg效應和不含Klinkenberg效應的 速度分布圖27（1psi=6894.76Pa） 除了用于致密巖石中滲流， 土壤滲流也是REV 尺度多孔介質LBE模型的研究熱點。 樊火30通過繪 制滲流場幾何剖面圖或者讀入圖形文件，格點化流 場、生成LBM格式的邊界數據文件，基于GLBE 型開發(fā)了土壤滲流LBM_Seepage軟件平臺， 生成數 據文件，以直觀可視的等值線或者云圖實現(xiàn)滲流場 速度實時可視化，并對豎向具有不同滲透率的土體 進行了模擬。結果表明，最上層土的相對的滲透能 力在滲透初期對滲流速度的傳播起控制作用，只有 經過一段時間后土體才會表現(xiàn)出等效應的滲透能力 （圖3） 。 申林芳等31進一步運用GLBE模型探討了 土體在壓力作用下孔隙率、滲透率及滲透壓力等影 響因素與滲流速度的相互關系。 Freed-LBE模型被推廣到非均勻網格中32，采 用HE等33提出的插入格式描述非均勻網格，該模 型在流動和碰撞過程之間采用插值法，不僅能保證 數值準確性而且易于實施。該模型對多孔介質、孔 隙和邊界的處理方式一樣，避免了不同尺度采用不 同數學方法而產生的誤差，與其他多尺度LB模型 相比，它消除了內部邊界，而且采用阻力模型，因 此可以用于大尺度。KANG等32在裂隙周圍細化網 格研究裂隙對流動的影響，遠離裂隙的地方采用粗 網格（圖4） ，模擬了復雜裂隙和簡化裂隙的多孔介 質內的流動，通過 u l K p = 預測REV尺度多孔介質 內的有效滲透率（表2） ，并提出可以把該模型用曲 線坐標和復合網格的設想。 REV尺度LBE模型不僅可以模擬含裂隙的均 質巖石中的流動，還能根據通過高分辨率的掃描圖 像確定空腔、黏土、長石、石英不同成分構成的非 均質砂巖結構，如圖5(a)，基于GLBE模型獲得砂 巖樣品內的流體流動速度，然后根據達西定律求得 整體的滲透率34-35。由圖5(b)可以看出，樣品大于 圖3 不同滲透率的的邊界示意圖及滲透率 K 由小到大和 K 由大到小中心線速度圖30 第6期 張瀟丹等：REV尺度多孔介質格子Boltzmann方法的數學模型及應用的研究進展 1705 圖4 裂隙示意圖、非均勻網格示意圖、計算結果32 表2 滲透率值32（格子單位） 項目 上半部分 下半部分 裂隙 整體 傳統(tǒng) LBE 模型 0.64 0.64 35.31 1.94 Freed-LBE 模型 0.61 0.61 37.01 1.99 誤差 4.7% 4.7% 4.8% 2.6% 10mm時整體的滲透率趨于常數，對于砂巖樣品， 大于10mm時可以達到REV尺度的要求。 由圖5(c) 和圖5(d)可以看出，流場主要受無滲透性的石英和 空腔影響，由于各成分的孔隙率不同導致各成分邊 界處的流線不閉合。但是CT掃描大孔時容易忽略 狹小的孔喉結構而造成誤差，因而掃描圖像的分辨 率必須要高36。 REV尺度多孔介質模型除了可以通過滲流速 度預測多孔介質的滲透率外，還可以分析多孔介質 圖5 砂巖樣品 幾何結構特性及外部條件對滲流速度的影響。 GLBE模型被用于在壓力驅動下的非牛頓流體的電 滲 流 過 程 37 ， 在 演 化 方 程 中 加 入 了 基 于 Herschel-Bulkley模型的非牛頓流體力項， 研究外部 電勢驅動流、離子濃度和電場強度、固體顆粒直徑 和孔隙度、冪律指數和屈服應力、德拜長度對滲透 速度的影響。多孔介質內的流速隨著冪律指數和屈 化 工 進 展 2016年第35卷 1706 服應力的減小而增大，壁面處的電勢主要影響壁面 處流速，多孔介質的電勢對通道中心處流速影響較 大。低冪律指數時，壁面處的電勢在對流速的影響 更加敏感，特別是在通道高度與德拜長度的比值較 小時。 HUANG等38對LBM兩相流模型及其應用 做了詳細描述， 目前LBM兩相流模型已經應用于 REV尺度多孔介質內多相流的驅替過程中。 SPAID和PHELAN39對纖維多孔介質的流動進行 了計算，纖維內部使用SP-LBE模型，外部使用 標準LBE模型，模擬圓形纖維內樹脂驅替空氣多 組分流動過程。纖維間的流動比纖維內的流動速 度快，導致纖維內部含有空氣空穴，空穴最初是 長軸垂直于流動方向的橢圓形，隨后由于表面張 力的影響，從橢圓形變?yōu)閳A形（圖6） 。根據達西 定律預測的不飽和滲透率在定性上正確。但是需 要更多研究證明纖維內的空洞隨時間消失模擬結 果是否準確。SCHAAP等 40應用阻力模型中 Shan-Chen格式多相流LBE模型，模擬了充滿玻 璃珠的多孔介質內水-空氣及水-石油精體系驅替 實驗。該模型基于顯微層析構造的多孔介質幾何 結構，忽略重力、慣性力和黏性力，只考慮毛細 作用力，設置合理的表面張力和接觸角，得到壓 力與飽和度的關系。水-空氣體系模擬值與測量值 相差不大，但是水-石油精體系模擬值與測量值差 異很大（圖7） ，原因是模擬條件與實驗條件不一 致，如接觸角。該方法也可以用于微尺度界面現(xiàn) 象的模擬。PORTER等41也運用該LBE模型估算 潤濕-非潤濕相界面的面積，模擬值與實驗值 吻合。 圖6 不同時刻單行圓形纖維樹脂驅替空氣模擬結果圖39 （黑色：空氣；灰色：樹脂） 圖7 壓力與飽和度演化圖 40 阻力模型中的單松弛GLBE模型推廣到多松弛 模型，LIU等42運用二維D2Q8 MRT-GLBE模型模 擬不可壓多孔介質流，在低黏性時的數值穩(wěn)定性比 單松弛模型好。隨后該模型被用于解含有時間導數 和非線性對流項平均體積動量方程，模擬了在多孔 介質圓柱靜止時上下壁面以恒定速度運動和上下壁 面靜止時管道內圓柱多孔介質以恒定速度運動的兩 種情況下的流動43。 結果表明R1和R2的本征相平 均速度吻合，二者的相平均速度在流體區(qū)吻合，但 是在多孔介質區(qū)域內有偏差圖8(b)，并且隨著孔 隙度的減小而增大。多孔介質阻力對本征相平均速 度和相平均速度有類似規(guī)律。只有含有本征相平均 速度的宏觀方程滿足伽利略不變性。多松弛模型與 單松弛模型相比，克服了許多限制條件，比如運動 黏度和體積黏度比不再是定值，但是多松弛模型增 加了計算量。 REV尺度LBM已經應用到工程滲流問題。何 瑩松44將REV尺度GLBE模型應用于煤礦開采過 程中回采工作面瓦斯?jié)B流問題，該研究中認為煤體 的孔隙度受應力重新分布的影響很小，孔隙度基本 不發(fā)生變化，滲透率沿流動方向呈現(xiàn)分段函數的不 均勻分布，得到了不同時刻非均質煤層中瓦斯壓力 分布和速度分布。這些數據能夠生動再現(xiàn)瓦斯流動 第6期 張瀟丹等：REV尺度多孔介質格子Boltzmann方法的數學模型及應用的研究進展 1707 圖8 多孔介質圓柱運動示意圖及=0.7、Re=100通道中心 固有相平均速度、相平均速度圖42 過程，有助于指導煤礦安全生產。 2.2 REV尺度LBE模型應用于多孔介質內流動與 傳熱耦合過程 流動、 傳熱和傳質兩個及兩個以上過程耦合時， LBE模型可以分為3種；即多速模型、雙分布函數 模型（DDF）與差分方法相結合的混合模型。目前 文獻大多采用DDF模型來模擬多過程， 該模型使用 兩個（及以上）分布函數演化方程，分別用于速度 場和溫度場（濃度場） ，格子結構簡單，且具有良好 的數值穩(wěn)定性。 等溫GLBE模型通過DDF模型推廣到多孔介 質熱流動中（DDF-GLBE模型）45，速度場用等溫 GLBE模型模擬，而溫度場使用一個新的分布函數 描述。DDF-GLBE模型被用于多孔介質方腔模擬孔 隙度恒定45-49、孔隙度變化50、含有熱源51的速度 場和溫度場，研究孔隙度、雷諾數（Re 數） 、達西 數（Da 數）等參數對速度場和溫度場的影響。 DDF-GLBE模型還被用于研究邊界充滿多孔介質 和中心充滿多孔介質兩種情況下管道內的流動與傳 熱過程52-53。 REV尺度的GLBE模型可以直接將孔隙尺度模 擬得到的滲透率和有效熱導率帶入REV模型， 避免 了采用經驗公式可能得到的錯誤結論54。模擬結果 表明，REV尺度得到的熱壁上的平均Nusselt數與 孔隙尺度下的趨勢一致， 證明REV尺度在定性上是 正確的。在瑞利數（Ra 數）低于107時，GLBE模 型能定量預測孔隙尺度下自然對流的宏觀規(guī)律，Ra 數大的時候只能定性預測孔隙尺度的宏觀規(guī)律（圖 9） 。REV尺度GLBE模型對滲透率如何影響多孔介 質自然對流宏觀現(xiàn)象刻畫很好，與孔隙尺度模擬結 果一致，但是往往夸大了流固熱導率帶來的影響。 圖9 不同Ra數下熱壁上平均Nusselt數隨方柱個數N的變 化結果對比54 虛線REV；實線孔隙尺度 REV尺度的LBE模型為太陽墻工程應用研究 提供參考， 通過含熱源項的LBE模型對太陽墻系統(tǒng) 中復雜流場及熱傳遞進行模擬55-56，可以預測太陽 墻在不同工作條件下的狀況，并優(yōu)化太陽墻系統(tǒng)的 設計。 針對工程中常見的多孔介質圓管內傳熱問題， 榮 伏 梅 等 57-58 構 造 了REV尺 度 軸 對 稱 熱 DDF-GLBE模型，包含弱可壓軸對稱LBE模型和 不可壓軸對稱LBE模型， 把刻畫多孔介質的參數加 入到分布函數中并構造合適的分布函數，使之能恢 復正確的宏觀方程，具有形式簡單、外力項不含復 雜梯度形式的優(yōu)點， 并應用于非線性滲流數值模擬。 分別研究了覆蓋多孔介質的方柱繞流、充滿多孔介 質內的變黏性Rayleigh-Benard對流、 以及圓管內插 入多孔介質強化對流傳熱非線性滲流問題。由于這 類問題的計算量都很大，單個CPU的計算效率較 低，因此采用GPU（圖形處理器）加速計算，可支 化 工 進 展 2016年第35卷 1708 持復雜的計算任務59。 單松弛DDF-GLBE模型推廣到多松弛模型中， 建立REV尺度二維多松弛（MRT）DDF-GLBE方 程，用D2Q9 MRT-GLBE模型模擬流場，同時用 D2Q5 MRT-LBE模型模擬溫度場， 模擬結果和解 析解與文獻計算結果吻合60，而且在低黏度時， MRT-LBE模型的數值穩(wěn)定性比LBGK模型好，但 是計算量大。該直角坐標系下的MRT-DDF-GLBE 模型延伸到軸對稱模型61。 2.3 REV尺度LBE模型應用于多孔介質內流動、 傳質和化學反應耦合過程 與流動和傳熱過程類同， 仍使用DDF模型描述 流動和傳質兩個過程，即運用兩個分布函數分別模 擬流場和濃度場，并在濃度演化方程中加入反應源 項， 用于模擬化學反應。TIAN等62基于DDF-GLBE 模型模擬了CO2地下填埋地球化學反應過程，研究 CO2注入過程中溶質離子反應及擴散、巖石骨架 CaCO3溶解相互作用， 引起孔隙率和滲透率變化， 反作用于流場。如圖10，開始在一段時間內出口處 的H+、Ca2+和HCO3離子濃度為0，隨著反應的進 行，3種離子濃度逐漸增加。由于反應消耗和對流 擴散的相互作用，在出口處反應物H+離子濃度在 1105時出現(xiàn)變化。Ca2+離子濃度到達最大后保持平 衡，是因為Ca2+離子是可逆反應。 該模型還用于研究有機廢水溶液中化學反應、 燃料電池和質子交換膜燃料電池中56 ，63-66。LIAO 等63和楊艷霞64研究了有機廢水中含擾流光合細 菌包埋顆粒的流動、傳質及其內部的光生化反應過 程，分析光照強度、進口流速、包埋顆粒的滲透率 及孔隙率對流場、濃度場及產氫性能的影響。XU 等65-66把該模型應用于固體燃料電池多組分化學反 應，計算濃差極化電勢，并研究了電池電極厚度、 電解質濃度、孔隙率、燃料組成等結構參數和工作 圖10 出口離子濃度圖62 溫度操作參數對濃差極化、超電勢、電化學效率、 電池電壓和功率密度的影響，LBM模擬結果與其他 模擬結果相比更接近實驗值。該模型也被用于模擬 質子交換膜