摘要：滑坡涌浪鏈生災(zāi)害的演化過(guò)程涉及滑坡體-河流之間復(fù)雜的相互作用，流固耦合方法在顆粒尺度下的應(yīng)用正日益凸顯。為了深入研究這一問(wèn)題，研究采用顆粒尺度的CFD（Computational Fluid Dynamics）-DEM（Discrete Element Method）流固耦合方法。通過(guò)編寫(xiě)UDF，以改進(jìn)計(jì)算流體力學(xué)軟件FLUENT中離散相模型只能進(jìn)行動(dòng)量交換而無(wú)法考慮排水效應(yīng)的缺點(diǎn)。這不僅改善了滑坡體與流體間的體積置換問(wèn)題，而且克服了網(wǎng)格顆粒臨界尺寸比的限制。在顆粒堆積體坍塌-涌浪試驗(yàn)?zāi)M中，數(shù)值模擬計(jì)算得到的最大涌浪高度為8.67 cm，與物理試驗(yàn)結(jié)果較為一致，證明了改進(jìn)后的CFD-DEM流固耦合模型能有效還原顆粒運(yùn)動(dòng)狀態(tài)及涌浪的波動(dòng)變化，能夠滿(mǎn)足滑坡涌浪過(guò)程模擬中滑坡體粗顆粒材料建模和高分辨率真實(shí)地形網(wǎng)格的需求。

關(guān)鍵詞：滑坡涌浪；流固耦合；CFD-DEM耦合方法；數(shù)值模擬

中圖分類(lèi)號(hào)：TV642.4文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1001-9235（2024）07-0111-07

郝悅彤，肖鴻.顆粒離散相模型的改進(jìn)及其在滑坡涌浪數(shù)值模擬應(yīng)用驗(yàn)證［J］.人民珠江，2024，45（7）：111-117.

Improvement of Discrete Phase Model for Particles and Its Validation for Application in Numerical Simulation of Landslide Surges

HAO Yuetong，XIAO Hong*

（State Key Laboratory of Hydraulics and Mountain River Engineering，Sichuan University，Chengdu 610065，China）

Abstract：The evolution of landslide surge chain-generated hazards involves complex landslide-stream interactions，and the application of fluid-solid coupling methods at the granular scale is becoming more and more prominent.In order to study this problem in depth，this study adopts the CFD（Computational Fluid Dynamics）-DEM（Discrete Element Method）fluid-solid coupling method at the granular scale.The UDF is written to improve the disadvantage that the discrete phase model in the Computational Fluid Dynamics software FLUENT can only perform momentum exchange but cannot consider the drainage effect.This not only improves the volume replacement problem between landslides and fluids，but also overcomes the limitation of the critical size ratio of mesh particles.In the simulation of particle accumulation body collapse-surge test，the maximum surge height obtained by numerical simulation is 8.67 cm，which is more consistent with the physical test results，proving that the improved CFD-DEM fluid-solid coupling model proposed in this paper can effectively restore the particle motion state and the fluctuation change of the surge，and it can meet the needs of the coarse-grained material modeling of landslides in the simulation of the landslide surge process and the high-resolution real topography grid in the simulation of landslide surge process.

Keywords：landslide surge;fluid-solid coupling;CFD-DEM coupling method;numerical simulation

庫(kù)岸邊坡由于庫(kù)水位波動(dòng)、暴雨、地震以及人為施工活動(dòng)等影響容易發(fā)生滑坡災(zāi)害［1］，滑坡涌浪作為滑坡體入水后的次生災(zāi)害，所產(chǎn)生的破壞性和影響范圍都超越了滑坡災(zāi)害本身［2］。大規(guī)?；麦w高速入水激起大量水體，水體回落形成涌浪在河道、庫(kù)區(qū)等水域內(nèi)傳播。一方面，涌浪的沖擊會(huì)對(duì)水工建筑物造成破壞甚至產(chǎn)生漫壩或潰壩的風(fēng)險(xiǎn)，對(duì)下游人民生命財(cái)產(chǎn)安全形成威脅；另一方面，滑坡體入水最終堆積在河道形成堰塞體會(huì)進(jìn)一步增加潰壩漫壩的風(fēng)險(xiǎn)［3］。

隨著計(jì)算機(jī)硬件水平的不斷提高，數(shù)值分析方法在研究庫(kù)岸災(zāi)害方面的應(yīng)用變得愈發(fā)重要［4］。河床沖刷、泥石流等諸多巖土工程問(wèn)題往往涉及復(fù)雜的流體-顆粒耦合作用，在對(duì)其機(jī)理進(jìn)行的研究中，顆粒尺度的流固耦合方法正發(fā)揮著越來(lái)越重要的作用［5-6］。利用離散元方法（Discrete Element Method，DEM）構(gòu)建巖土體，結(jié)合常用于流體流動(dòng)建模的計(jì)算流體動(dòng)力學(xué)（Computational Fluid Dynamics，CFD）方法，Tsuji等［7］提出了一種顆粒尺度的流固耦合方法，即CFD-DEM耦合方法，這種方法成功用于模擬流化床問(wèn)題。在CFD-DEM耦合方法中，流體運(yùn)動(dòng)由Euler框架下的連續(xù)方法描述，固體顆粒在考慮粒間碰撞以及流體-顆粒作用的基礎(chǔ)上由Lagrange框架下的離散方法描述其微觀(guān)特性，在空間和時(shí)間尺度上都遠(yuǎn)小于連續(xù)方法描述的宏觀(guān)特性，需要建立宏觀(guān)尺度（連續(xù)介質(zhì)）和微觀(guān)尺度（離散介質(zhì)）的橋梁，也就是流體體積分?jǐn)?shù)［8］。目前，CFD軟件FLUENT中的離散相模型對(duì)顆粒處理尚不完善，僅考慮了動(dòng)量交換而未考慮顆粒與流體之間的體積置換，對(duì)于滑坡涌浪問(wèn)題的模擬計(jì)算會(huì)產(chǎn)生較大誤差。此外，傳統(tǒng)的CFD-DEM耦合方法的使用受限于網(wǎng)格顆粒臨界尺寸比［9-10］，難以滿(mǎn)足滑坡涌浪過(guò)程模擬中滑坡體粗顆粒材料建模和高分辨率真實(shí)地形網(wǎng)格的需求［11］。

本文旨在優(yōu)化FLUENT中CFD-DEM耦合數(shù)值模擬方法，并驗(yàn)證其進(jìn)行滑坡涌浪數(shù)值模擬的可靠性。首先引入虛擬多孔介質(zhì)球孔隙率計(jì)算模型，不僅考慮到各相間體積置換，而且可以解決實(shí)際應(yīng)用過(guò)程中顆粒尺寸大于網(wǎng)格尺寸時(shí)產(chǎn)生計(jì)算不穩(wěn)定的缺點(diǎn)；其次通過(guò)單顆粒沉降模擬、顆粒團(tuán)入水及出水模擬和顆粒堆積體坍塌-涌浪試驗(yàn)?zāi)M，驗(yàn)證了模型的準(zhǔn)確性和有效性。研究成果可以為相關(guān)地區(qū)重大水電工程開(kāi)發(fā)運(yùn)行的安全保障提供參考。

1流固耦合數(shù)值方法

1.1計(jì)算流體運(yùn)動(dòng)學(xué)方法

1.1.1控制方程

對(duì)于流體相，采用以下基本方程：連續(xù)性方程見(jiàn)（1）：

+++=0（1）

動(dòng)量方程見(jiàn)（2）：

式中：u、v、w分別為x、y、z 3個(gè)方向上的速度分量；ρ為密度；t為時(shí)間；μ為動(dòng)力黏度；p為壓強(qiáng)；fx、fy、fz分別為流體單元在3個(gè)方向上所受的體積力，一般考慮為重力；fpfx、fpfy、fpfz分別為流體所受顆粒-流體相互作用在3個(gè)方向上的大小。

1.1.2 VOF方法

使用VOF（Volume of Fluid）方法捕捉自由面，對(duì)涌浪的形成及演變過(guò)程中多相流體交界面進(jìn)行追蹤。通過(guò)針對(duì)多相的體積分?jǐn)?shù)的連續(xù)性方程求解來(lái)完成對(duì)相之間的界面跟蹤，對(duì)于第q相流體，

式中：αq為第q相流體在某個(gè)計(jì)算單元的體積分?jǐn)?shù)值；pq為從第p相到第q相的質(zhì)量輸移；qp為從第q相到第p相的質(zhì)量輸移，默認(rèn)右側(cè)源項(xiàng)sα為0。

當(dāng)αq=0時(shí)，代表該單元格沒(méi)有該種流體；當(dāng)αq=1，代表該單元格充滿(mǎn)該種流體；當(dāng)0lt;αq＜1，代表該單元格為該流體和其他流體的交界面。根據(jù)單元格上的αq值，分配計(jì)算域中每個(gè)單元格的屬性及其存儲(chǔ)的變量，見(jiàn)式（4）、（5）。

式中：腳標(biāo)對(duì)應(yīng)不同種類(lèi)的流體。

1.2離散元方法

對(duì)于顆粒相，根據(jù)牛頓第二定律，每個(gè)粒子切向的運(yùn)動(dòng)方程為式（6）：

mi=（Fn，k，i+Ft，k，i）+mig+Fi（6）

式中：mi為顆粒i的質(zhì)量；Fn，k，i為相鄰顆粒k產(chǎn)生的法相接觸力；Ft，k，i為相鄰顆粒k產(chǎn)生的切向接觸力；Fi為流體作用在顆粒i上的力。

單個(gè)顆粒的旋轉(zhuǎn)運(yùn)動(dòng)方程見(jiàn)式（7）：

Ii=aεijlniFt，k，l（7）

式中：Ii為顆粒i的轉(zhuǎn)動(dòng)慣量；ωi為顆粒i的轉(zhuǎn)動(dòng)角速度；a為顆粒的半徑；ni為顆粒i表面的單位法向量。

1.3虛擬多孔介質(zhì)球孔隙率模型

計(jì)算網(wǎng)格孔隙率是計(jì)算顆粒相和流體相之間相互作用力的關(guān)鍵［12］。傳統(tǒng)CFD-DEM使用的方法有中心孔隙率計(jì)算方法（Centre Void Fraction Method）和劃分孔隙分?jǐn)?shù)法（Divided Void Fraction Method）。前者誤差較大且顆粒體積大于網(wǎng)格體積時(shí)會(huì)出現(xiàn)流場(chǎng)不連續(xù)的情況，導(dǎo)致整個(gè)CFD-DEM模擬計(jì)算的數(shù)值不穩(wěn)定；后者雖能得到相對(duì)光滑的體積分?jǐn)?shù)場(chǎng)，但在非均勻網(wǎng)格中較難實(shí)現(xiàn)，同時(shí)也無(wú)法解決顆粒體積大于網(wǎng)格體積時(shí)的數(shù)值不穩(wěn)定問(wèn)題。在應(yīng)用過(guò)程中需要滿(mǎn)足網(wǎng)格尺寸與顆粒粒徑尺寸比大于4的限制［8］。

本文引入虛擬多孔介質(zhì)球孔隙率模型，將顆?？紤]為3～5倍大的多孔介質(zhì)球，將顆粒在原本位置體積占有的影響平均到比原有位置處顆粒體積大幾倍的網(wǎng)格上，最后再將所有顆粒對(duì)其周?chē)W(wǎng)格的影響疊加總和，從而完成了顆粒拉格朗日法到流體計(jì)算域歐拉場(chǎng)的映射轉(zhuǎn)換。虛擬多孔介質(zhì)球孔隙率模型對(duì)于實(shí)際顆粒粒徑尺寸接近于網(wǎng)格尺寸時(shí)產(chǎn)生的數(shù)值誤差較小，對(duì)于多孔介質(zhì)球粒徑尺寸為實(shí)際顆粒粒徑3～5倍情況下的計(jì)算穩(wěn)定性更好。

每個(gè)顆粒影響區(qū)域下覆蓋著很多個(gè)流體計(jì)算域網(wǎng)格，每個(gè)顆粒的影響域都會(huì)在其覆蓋流體計(jì)算網(wǎng)格上產(chǎn)生一個(gè)體積占有值，代表著每個(gè)網(wǎng)格被顆粒i所占有的體積比值（圖1）。對(duì)所有顆粒執(zhí)行相同的操作，最終每個(gè)網(wǎng)格上都有被不同顆粒所影響的體積比值，將這些值在每個(gè)流體網(wǎng)格上疊加，就能得到流場(chǎng)中的顆粒所占有體積的信息，見(jiàn)式（8）：

式中：i為第i個(gè)顆粒；di，j為顆粒與周?chē)黧w網(wǎng)格中心的距離；Vp，i為第i個(gè)顆粒的體積；Vcell，j為第j個(gè)流體網(wǎng)格的體積；a為多孔介質(zhì)球與實(shí)際顆粒的半徑尺寸比，這里取值為4。

對(duì)于每個(gè)網(wǎng)格，被固體顆粒占據(jù)的體積與孔隙率之和為1，孔隙率見(jiàn)式（9）：

1.4流體與顆粒之間的相互作用

流體與顆粒之間的相互作用是流-固耦合模型連接流場(chǎng)信息與顆粒信息的主要內(nèi)容，流體與顆粒之間相互作用力計(jì)算的是否準(zhǔn)確，是流-固耦合成敗的關(guān)鍵［13］。模擬設(shè)置的顆粒密度遠(yuǎn)大于液相密度，流體重力、虛擬質(zhì)量力、Basset力、Saffman力和Magnus力在所有受力中所占量級(jí)較小，研究中將其忽略，主要考慮相間拖曳力的作用。本文考慮目前常用的2種拖曳力模型并進(jìn)行驗(yàn)證。

a）Ergun、Wen和Yu模型。結(jié)合計(jì)算密集顆粒堆積體的Ergun公式以及計(jì)算稀疏顆粒堆積體的Wen和Yu公式，得到Ergun、Wen和Yu模型，見(jiàn)式（10）、（11）［14］：

其中：

式中：fp為流體對(duì)顆粒產(chǎn)生的拖曳力；ε為網(wǎng)格孔隙率；dp為顆粒直徑；ufi為顆粒球心處的流體速度，通常采用顆粒球心所在流體網(wǎng)格中的流速；up為顆粒速度；CP為拖曳力系數(shù)；Rep為顆粒雷諾數(shù)，見(jiàn)式（12）：

Rep=（12）

b）Di Felice模型。Di Felice［15］通過(guò)對(duì)顆粒堆積體內(nèi)滲流實(shí)驗(yàn)數(shù)據(jù)進(jìn)行整理和擬合，提出了Di Felice模型公式，此模型認(rèn)為顆粒堆積體內(nèi)流體作用于顆粒的拖曳力可以根據(jù)單顆粒在流體內(nèi)平穩(wěn)沉降時(shí)受到的拖曳力fp0乘一個(gè)關(guān)于堆積體孔隙率ε的函數(shù)f（ε）來(lái)獲得，見(jiàn)式（13）：

fp=f（ε）fp0（13）

其中fp0見(jiàn)式（14）：

其中，Cp0為拖曳力系數(shù)，見(jiàn)式（15）：

Cp0=（0.63+）2

f（ε）為冪函數(shù)形式，見(jiàn)式（16）：

f（ε）=ε2-Χ

其中，見(jiàn)式（17）：

Χ=3.7-0.65exp-

1.5耦合計(jì)算流程

流體相和顆粒相相互作用的耦合需要得到流體對(duì)顆粒的作用和顆粒對(duì)流體的作用。因此耦合方案由2個(gè)部分構(gòu)成：一是獲得顆粒周邊流體的性質(zhì)，為顆粒模型提供邊界條件；二是獲得并更新相間信息，如體積分?jǐn)?shù)。本文的CFD-DEM耦合計(jì)算步驟依賴(lài)于ANSYS FLUENT的DPM離散相模型的應(yīng)用，其計(jì)算流程見(jiàn)圖2。

2單顆粒沉降模擬

球形顆粒的基本沉降是流體-顆粒相互作用中的重要問(wèn)題，選擇此類(lèi)顆粒沉降問(wèn)題，以驗(yàn)證流體和顆粒的相互作用力是否計(jì)算準(zhǔn)確。選取的計(jì)算域尺寸為4 mm×4 mm×8 mm，初始水面高4 mm，球形顆粒從水面上方1 mm的位置處由靜止?fàn)顟B(tài)自由落體。顆粒直徑0.1 mm，顆粒密度2 500 kg/m3。空氣相密度1 kg/m3，運(yùn)動(dòng)黏度1×10-5 Pa·s；水相密度1 000 kg/m3，運(yùn)動(dòng)黏度0.002 Pa·s。網(wǎng)格尺寸設(shè)置為0.1 mm，網(wǎng)格尺寸和顆粒粒徑尺寸比為1。

對(duì)于低雷諾數(shù)，顆粒的沉降速度規(guī)律符合Stokes定律使用范圍下的理論解，顆粒入水后的速度隨時(shí)間的變化關(guān)系見(jiàn)式（18）：

對(duì)2種拖曳力模型進(jìn)行對(duì)比驗(yàn)證，結(jié)果見(jiàn)圖3。

由圖3可知，虛擬多孔介質(zhì)球孔隙率模型可以應(yīng)用于計(jì)算網(wǎng)格顆粒尺寸比為1時(shí)的情況，有效消除了傳統(tǒng)模型中網(wǎng)格尺寸與顆粒粒徑比大于4的限制。同時(shí)，在該情形下，采用Di Felice拖曳力模型的結(jié)果與Stokes理論解具有較好的一致性，最終沉降速度為0.004 87 m/s。因此，本文選取Di Felice拖曳力模型。

在0.015～0.020 s模型與理論解相差較大，這是由于Stokes理論解的計(jì)算分為水面上自由落體和水面下顆粒沉降，水面為定值；而數(shù)值模擬時(shí)，水面會(huì)被顆粒擾動(dòng)，水面處液相和顆粒相之間有相互作用。在拖曳力的驗(yàn)證中，主要依據(jù)最終沉降速度來(lái)驗(yàn)證拖曳力計(jì)算的準(zhǔn)確性。本文中顆粒最終達(dá)到了理想的穩(wěn)定沉降速度，可以驗(yàn)證模型的可靠性。

3顆粒入水及出水模擬驗(yàn)證

顆粒團(tuán)入水及出水涉及復(fù)雜的水面運(yùn)動(dòng)情況以及流體-顆粒之間相互作用。顆粒團(tuán)進(jìn)入水中，會(huì)排開(kāi)相應(yīng)體積的水，水面上升；顆粒浮出水面時(shí)，顆粒原本占據(jù)的體積會(huì)被液體填充，自由面高度會(huì)下降。通過(guò)模擬顆粒團(tuán)入水及出水的過(guò)程，驗(yàn)證改進(jìn)后的模型對(duì)于固體顆粒、水和空氣三相體積處理的守恒性。

流體計(jì)算域尺寸為0.05 m×0.05 m×0.15 m，初始時(shí)刻自由水面位置為z=0.05 m處，上方空氣相密度1 kg/m3，運(yùn)動(dòng)黏度1×10-5 Pa·s；下方水相密度1 000 kg/m3，運(yùn)動(dòng)黏度1×10-3 Pa·s。初始時(shí)刻顆粒像柵格一樣在水面上方排列，其x、y、z方向上數(shù)目分別是20×20×30，共12 000個(gè)（顆粒按初始時(shí)刻z方向的位置每10層分為一組，紅色、綠色、藍(lán)色分別表示上、中、下3組顆粒），顆粒粒徑為0.002 m，顆粒間球心間距略大于顆粒粒徑為0.002 01 m，顆粒密度為2 500 kg/m3。

DEM時(shí)間步長(zhǎng)設(shè)置為2.0×10-5 s，CFD時(shí)間步長(zhǎng)設(shè)置為1.0×10-4 s。計(jì)算2 s后顆粒沉降結(jié)束自由面恢復(fù)平靜，此時(shí)更改水體密度為7 500 kg/m3，黏度不變，再計(jì)算2 s顆粒浮出水面并達(dá)到穩(wěn)定。

圖4顯示，未改進(jìn)的模型在T=2 s時(shí)，顆粒完全入水，水面仍為初始液面高度，未考慮顆粒相的體積排水作用。改進(jìn)模型后，T=2 s時(shí)顆粒完全入水，水面上升的高度為顆粒團(tuán)全部固體顆粒的體積，水面達(dá)到理論值0.070 1 m，T=4 s時(shí)顆粒團(tuán)上浮運(yùn)動(dòng)結(jié)束，水面上升的高度為被淹沒(méi)顆粒所排開(kāi)水的體積，水面為理論值0.058 0 m。證明本文模型有效改善了FLUENT中離散相模型無(wú)法考慮顆粒排水效應(yīng)的缺點(diǎn)。

4顆粒堆積體坍塌-涌浪試驗(yàn)?zāi)M

顆粒堆積體坍塌-涌浪試驗(yàn)可以驗(yàn)證CFD-DEM流固耦合模型模擬滑坡-涌浪災(zāi)害的有效性［16］。Robbe-Saule等［17］所做的室內(nèi)水槽試驗(yàn)包含了顆粒堆積體失穩(wěn)、涌浪形成和傳播等類(lèi)似滑坡涌浪災(zāi)害鏈的過(guò)程，被廣泛用作標(biāo)準(zhǔn)檢驗(yàn)?zāi)Ｐ汀?/p>

試驗(yàn)裝置左側(cè)為顆粒堆積區(qū)，高0.41 m，寬0.165 m，總體積0.102 m3；右側(cè)為一個(gè)長(zhǎng)2.00 m、高0.30 m、寬0.15 m的玻璃水槽，水槽中充滿(mǎn)深度為0.05 m的水。數(shù)值模擬計(jì)算參數(shù)見(jiàn)表1，DEM碰撞參數(shù)見(jiàn)表2。

圖5為顆粒堆積體坍塌過(guò)程不同時(shí)刻水槽試驗(yàn)結(jié)果（每組狀態(tài)左圖）與數(shù)值模擬結(jié)果（每組狀態(tài)右圖）的對(duì)比。通過(guò)對(duì)比圖中顆粒堆積體坍塌過(guò)程及水面涌浪運(yùn)動(dòng)狀態(tài)可以看出，在整個(gè)演進(jìn)過(guò)程中，數(shù)值模擬結(jié)果與水槽試驗(yàn)結(jié)果相差較小。當(dāng)顆粒與水面相撞時(shí)會(huì)形成初始涌浪，隨著更多的顆粒侵入水體，涌浪高度逐漸增加；隨后，涌浪在傳播過(guò)程中逐漸變陡，波峰在重力影響下濺落到水體表面。對(duì)比水槽試驗(yàn)結(jié)果和數(shù)值模擬結(jié)果可知，CFD-DEM流固耦合模型能夠較為準(zhǔn)確地捕捉這一現(xiàn)象，初步驗(yàn)證了該模型解決滑坡涌浪問(wèn)題的可靠性。

圖6顯示，未改進(jìn)的模型由于未考慮顆粒的排水效應(yīng)會(huì)低估涌浪高度，改進(jìn)模型后可縮減涌浪高度模擬的誤差。水槽試驗(yàn)中測(cè)量的最大涌浪高度為8.30 cm，數(shù)值模擬計(jì)算得到的最大涌浪高度為8.67 cm，結(jié)果較為一致。數(shù)值模擬形成涌浪達(dá)到最大涌浪高度的時(shí)間早于水槽試驗(yàn)結(jié)果約0.1 s，可能是由于物理試驗(yàn)提拉擋板的時(shí)間誤差及其對(duì)顆粒的擾動(dòng)造成的。同時(shí)數(shù)值模擬將不規(guī)則顆粒簡(jiǎn)化為單一粒徑球形顆粒，低估了顆粒間以及顆粒與流體間的作用力，造成結(jié)果偏大。該研究結(jié)果與肖華波等［16］數(shù)值模擬結(jié)果類(lèi)似。通過(guò)對(duì)顆粒坍塌運(yùn)動(dòng)過(guò)程和涌浪高度演化過(guò)程的對(duì)比，發(fā)現(xiàn)顆粒運(yùn)動(dòng)及涌浪傳播模擬結(jié)果與物理試驗(yàn)結(jié)果較為吻合，證明了本文所提出改進(jìn)后的CFD-DEM流固耦合模型能較好還原顆粒運(yùn)動(dòng)狀態(tài)及涌浪的波動(dòng)變化。

5結(jié)論

本文對(duì)FLUENT中CFD-DEM流固耦合模型進(jìn)行改進(jìn)，并進(jìn)行多角度驗(yàn)證，得到以下主要結(jié)論。

a）通過(guò)引入虛擬多孔介質(zhì)球孔隙率模型，改進(jìn)了FLUENT離散相模型無(wú)法做到考慮顆粒體積分?jǐn)?shù)的缺點(diǎn)，克服了傳統(tǒng)CFD-DEM模型網(wǎng)格顆粒臨界尺寸限制在實(shí)際應(yīng)用方面的困難。

b）通過(guò)單顆粒沉降模擬，驗(yàn)證了網(wǎng)格顆粒尺寸比為1時(shí)，采用Di Felice拖曳力模型的結(jié)果與理論解具有較好的一致性，最終沉降速度為0.004 87 m/s。

c）通過(guò)多顆粒入水及出水模擬，顆粒完全入水后，水面上升高度達(dá)到理論值0.070 1 m，驗(yàn)證了虛擬多孔介質(zhì)球孔隙率模型可以有效改進(jìn)FLUENT離散相模型，考慮顆粒體積的排水效應(yīng)，從而實(shí)現(xiàn)強(qiáng)耦合。

d）通過(guò)顆粒堆積體坍塌-涌浪試驗(yàn)?zāi)M，驗(yàn)證了改進(jìn)后的CFD-DEM流固耦合模型能有效還原顆粒運(yùn)動(dòng)狀態(tài)及涌浪的波動(dòng)變化，數(shù)值模擬計(jì)算得到的最大涌浪高度為8.67 cm，與物理試驗(yàn)結(jié)果較為一致，驗(yàn)證了該模型解決滑坡涌浪問(wèn)題的可靠性。

參考文獻(xiàn)：

［1］周家文，陳明亮，瞿靖昆，等.水庫(kù)滑坡災(zāi)害致災(zāi)機(jī)理及防控技術(shù)研究與展望［J］.工程科學(xué)與技術(shù)，2023，55（1）：110-128.

［2］魏云鵬.基于CFD-DEM方法的散粒體滑坡涌浪研究［D］.昆明：昆明理工大學(xué)，2022.

［3］王雷，解明禮，黃會(huì)寶，等.不同規(guī)模滑坡入水誘發(fā)涌浪災(zāi)害特征差異性分析［J］.人民珠江，2024，45（2）：18-28.

［4］鄧成進(jìn)，黨發(fā)寧，陳興周，等.庫(kù)區(qū)滑坡涌浪三維數(shù)值模擬分析［J］.水利水運(yùn)工程學(xué)報(bào)，2020（6）：64-71.

［5］WANG T，ZHANG F S，F(xiàn)URTNEY J，et al.A review of methods，applications and limitations for incorporating fluid flow in the discrete element method［J］.Journal of Rock Mechanics and Geotechnical Engineering，2022，14：1005-1024.

［6］ZHANG P Y，MU L L，HUANG M S.A coupled CFD-DEM investigation into hydro-mechanical behaviour of gap-graded soil experiencing seepage erosion considering cyclic hydraulic loading［J］.Journal of Hydrology，2023，624.DOI：10.1016/j.jhydrol.2023.129908.

［7］TSUJI Y，KAWAGUCHI T，TANAKA T.Discrete particle simulation of two-dimensional fluidized bed［J］.Powder Technology，1993，77（1）：79-87.

［8］劉德天，傅旭東，王光謙.CFD-DEM耦合計(jì)算中的體積分?jǐn)?shù)算法［J］.清華大學(xué)學(xué)報(bào)（自然科學(xué)版），2017，57（7）：720-727.

［9］ZHAO T，HOULSBY G T，UTILI S.Investigation of granular batch sedimentation via DEM-CFD coupling［J］.Granular Matter，2014，16：921-932.

［10］姚鵬，王遠(yuǎn)，冒劉鵬，等.CFD-DEM模型最優(yōu)網(wǎng)格尺寸的理論推證［J］.泥沙研究，2023，48（5）：1-7，34.

［11］李東陽(yáng)，年廷凱，吳昊，等.滑坡-堵江-涌浪災(zāi)害鏈模擬的DEM-CFD耦合分析方法及其應(yīng)用［J］.工程科學(xué)與技術(shù)，2023，55（1）：141-149.

［12］吳亮.含固體顆粒的兩相流界面變化的數(shù)值研究：DEM-VOF方法的實(shí)現(xiàn)［D］.天津：天津大學(xué)，2018.

［13］彭愷然，劉紅帥，平新雨，等.CFD-DEM耦合模擬中拖曳力模型精度［J］.吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2021，51（5）：1400-1407.

［14］ERGUN S.Fluid Flow Through Packed Columns［J］.Chemical Engineering Progress，1952，48：89-94.

［15］DI FELICE R.The Voidage Function for Fluid-Particle Interaction Systems［J］.International Journal of Multiphase Flow，1994，20（1）：153-159.

［16］肖華波，王澤皓，石偉明，等.猴子巖水庫(kù)色玉滑坡涌浪災(zāi)害鏈CFD-DEM耦合數(shù)值模擬［J］.工程科學(xué)與技術(shù)，2023，55（1）：161-170.

［17］ROBBE-SAULE M，MORIZE C，HENAFF R，et al.Experimental in-vestigation of tsunami waves generated by granular col-lapse into water［J］.Journal of Fluid Mechanics，2021，907.DOI：/10.1017/jfm.2020.807.

（責(zé)任編輯：程茜）