合作客戶/
拜耳公司 |
同濟(jì)大學(xué) |
聯(lián)合大學(xué) |
美國保潔 |
美國強(qiáng)生 |
瑞士羅氏 |
相關(guān)新聞Info
-
> 溴化鋰及離子液體水溶液密度、黏度和表面張力測定與計(jì)算
> 定性分析聚合物界面張力與系統(tǒng)黏度比之間的關(guān)系——數(shù)值模擬、實(shí)驗(yàn)研究結(jié)果
> 以大豆為原料合成的N-椰子油酰基復(fù)合氨基酸表面活性劑表面張力、乳化起泡潤濕性能測定（二）
> 全自動(dòng)張力測定儀揭示子細(xì)胞表面張力對胞質(zhì)分裂結(jié)局的主導(dǎo)作用（二）
> 陽-非離子復(fù)合表面活性劑體系表面張力測定及基礎(chǔ)性能評價(jià)（三）
> 液滴爆炸現(xiàn)象：酒精蒸發(fā)引起的馬蘭戈尼流動(dòng)現(xiàn)象影響參數(shù)（一）
> 不同類型的堿、pH值對孤東油田原油界面張力的影響（下）
> 納米顆粒間相互作用對界面張力的影響
> 氫氧化鈉用量對重烷基苯磺酸鈉水溶液/油體系界面張力的影響（二）
> 基于深度神經(jīng)網(wǎng)絡(luò)模型分析明膠溶液荷電量與表面張力之間的關(guān)系（一）
推薦新聞Info
-
> 面向高效環(huán)保滅火劑的界面張力最小化：短鏈氟碳復(fù)配體系的設(shè)計(jì)與解析（三）
> 面向高效環(huán)保滅火劑的界面張力最小化：短鏈氟碳復(fù)配體系的設(shè)計(jì)與解析（二）
> 面向高效環(huán)保滅火劑的界面張力最小化：短鏈氟碳復(fù)配體系的設(shè)計(jì)與解析（一）
> 基于連續(xù)表面張力模型微觀層面研究凝結(jié)顆粒動(dòng)力學(xué)變化規(guī)律及能量耗散的影響（四）
> 基于連續(xù)表面張力模型微觀層面研究凝結(jié)顆粒動(dòng)力學(xué)變化規(guī)律及能量耗散的影響（三）
> 基于連續(xù)表面張力模型微觀層面研究凝結(jié)顆粒動(dòng)力學(xué)變化規(guī)律及能量耗散的影響（二）
> 基于連續(xù)表面張力模型微觀層面研究凝結(jié)顆粒動(dòng)力學(xué)變化規(guī)律及能量耗散的影響（一）
> 基于表面張力、Walker沉降等研究3種表面活性劑對低階煤(LRC)潤濕作用的差異（三）
> 基于表面張力、Walker沉降等研究3種表面活性劑對低階煤(LRC)潤濕作用的差異（二）
> 基于表面張力、Walker沉降等研究3種表面活性劑對低階煤(LRC)潤濕作用的差異（一）
基于連續(xù)表面張力模型微觀層面研究凝結(jié)顆粒動(dòng)力學(xué)變化規(guī)律及能量耗散的影響（二）
來源：化工進(jìn)展 瀏覽 30 次 發(fā)布時(shí)間:2026-03-02
1 數(shù)值計(jì)算方法
濕顆粒碰撞過程的直接數(shù)值模擬通過流體體積模型（volume of fluid，VOF）模型追蹤氣液界面，連續(xù)表面張力模型（continuum surface force，CSF）模型計(jì)算界面的表面張力，兩者共同作用，控制液膜的形變過程。顆粒運(yùn)動(dòng)則采用六自由度模型（six degrees of freedom，6 DOF）結(jié)合重疊網(wǎng)格技術(shù)進(jìn)行描述，為了模擬顆粒的反彈運(yùn)動(dòng)并提取分析數(shù)據(jù)，本文開發(fā)并應(yīng)用了用戶自定義函數(shù)。
1.1 VOF模型
采用VOF模型來模擬氣液兩相流體的流動(dòng)問題，通過求解體積分?jǐn)?shù)的對流方程實(shí)現(xiàn)相界面的瞬態(tài)追蹤。在VOF模型中，通過跟蹤不同相的體積分?jǐn)?shù)來描述不同相的位置和界面形狀[式(1)]。
?α?t+??(αu)=0(1)
式中，α為液體的體積分?jǐn)?shù)，%，α=1處的網(wǎng)格充滿液體，α=0處的網(wǎng)格充滿氣體，0<α<1的網(wǎng)格區(qū)域含有自由面；t為時(shí)間，s；u為流體的速度矢量，m/s；?為梯度算子。
假設(shè)流體為不可壓縮流體，且氣液兩相的密度保持不變，涉及質(zhì)量與動(dòng)量守恒。
連續(xù)性方程為式(2)。
??u=0(2)
動(dòng)量方程為式(3)。
?(ρu)?t+u??(ρu)=-?p+??[μ?(?u+?uT)]+g+FSV(3)
式中，ρ為流體的局部平均密度，kg/m3；p為流體的壓力，Pa；μ為動(dòng)力黏度，Pa·s；g為流體自身重力，N；FSV為表面張力的等價(jià)體積力，N。
局部平均密度ρ和動(dòng)力黏度μ表達(dá)式為式(4)和式(5)。
ρ=αρl+(1-α)ρg(4)
μ=αμl+(1-α)μg(5)
式中，下角標(biāo)l為液體相；下角標(biāo)g為氣體相。
1.2 連續(xù)表面張力模型
連續(xù)表面張力模型能夠?qū)⒁后w界面上的力轉(zhuǎn)換為周圍的體積內(nèi)連續(xù)分布的力，被廣泛應(yīng)用于液滴流動(dòng)模擬中。界面處表面張力可以表示為式(6)。
Fsf=σκ?φ(6)
式中，F(xiàn)sf為表面張力，N；σ為表面張力系數(shù)；?φ為界面切向的梯度算子；κ為表面曲率。
界面處的單位法向量計(jì)算如式(7)、式(8)。
n=-?φ|?φ|(7)
n=nwcosθ+twcosθ(8)
式中，θ為氣-液界面在壁面處的接觸角，(°)；nw和tw分別為與壁面法向和切向的單位向量。得到界面的曲率為式(9)。
κ=???φ|?φ|(9)
1.3 顆粒運(yùn)動(dòng)模型
控制顆粒運(yùn)動(dòng)位移的方程為式(10)。
dupdt=1mpFp+Fv+Fsf+Fg(10)
式中，F(xiàn)p為壓差阻力，N；Fv為黏性阻力，N；Fg為重力，N；mp為顆粒質(zhì)量，kg；up為顆粒速度，m/s。
為進(jìn)一步分析濕顆粒碰撞過程中的能量損失，可以將其分為以下四部分：壓差阻力引起的能量損失比Ep、黏性阻力引起的能量損失比Ev、由濕顆粒間非完全彈性碰撞引起的能量損失比Ec、由表面張力引起的動(dòng)能耗損比Esf。并將顆粒1剛接觸液膜時(shí)刻的初始動(dòng)能定義為E0。計(jì)算見式(11)~式(13)。
Ep=(∫h0htFpdh+∫hthmFpdh)/E0Fp=-∫Spiz?ndS(11)
Ev=(∫h0htFvdh+∫hthmFvdh)/E0Fv=μ∫S?uz?ndS(12)
Ec=(12mut2-12mus2)/E0(13)
式中，h0為顆粒1接觸液膜時(shí)質(zhì)心高度，μm；ht為顆粒1與顆粒2碰撞時(shí)質(zhì)心高度，μm；hm為顆粒1碰撞后反彈到最高點(diǎn)的質(zhì)心高度，μm；ut為顆粒1碰撞前的速度，m/s；us為顆粒1碰撞后的速度，m/s。
1.4 重疊網(wǎng)格
本文使用的重疊網(wǎng)格方法廣泛應(yīng)用于模擬移動(dòng)網(wǎng)格的邊界運(yùn)動(dòng)問題。重疊網(wǎng)格技術(shù)的工作原理包括挖孔、重疊最小化和供體搜索。首先，生成互不約束的背景網(wǎng)格和重疊網(wǎng)格，滿足N-S求解器的三階離散化。在求解的迭代過程中，流場數(shù)據(jù)通過內(nèi)嵌邊界進(jìn)行傳輸。其次，確定執(zhí)行各自職責(zé)的單元，包括挖孔單元、重疊邊界單元（接收單元）和計(jì)算單元（供體單元）。挖孔單元是指通過挖掘操作丟棄的一些煩瑣單元。重疊邊界單元主要圍繞挖孔單元，并被視為每個(gè)重疊網(wǎng)格中的邊界。重疊邊界單元與計(jì)算單元之間的關(guān)系通過插值確定。具體實(shí)現(xiàn)流程如圖1所示。
圖1 重疊網(wǎng)格計(jì)算求解流程