基于SPH-DEM耦合方法的液滴沖擊散粒體運(yùn)動(dòng)機(jī)理研究.pdf

1、液滴沖擊離散顆粒是自然界和工程領(lǐng)域中常見(jiàn)的一種現(xiàn)象，也是一種典型而又復(fù)雜的不可壓縮自由表面流動(dòng)和離散顆粒結(jié)構(gòu)耦合作用問(wèn)題。此類(lèi)問(wèn)題的深入研究為復(fù)雜的流體和顆粒間耦合作用的理論和數(shù)值模擬研究提供了新思路，具有重要的理論意義。在工程應(yīng)用方面為泥沙啟動(dòng)、土壤侵蝕、河床演變和地質(zhì)災(zāi)害預(yù)報(bào)等問(wèn)題的研究提供了新的研究手段，具有實(shí)際應(yīng)用價(jià)值。
　　本文將光滑粒子流體動(dòng)力學(xué)（SPH）方法與離散單元法（DEM）方法相結(jié)合，進(jìn)行了液滴沖擊散粒體顆粒問(wèn)

2、題的研究。采用SPH方法模擬了液滴撞擊帶孔壁面的問(wèn)題，提出了隨計(jì)算區(qū)域變化的鏈表搜索法。結(jié)合實(shí)驗(yàn)進(jìn)一步研究了不同物理?xiàng)l件下黏性、重力和內(nèi)部壓應(yīng)力波動(dòng)對(duì)鋪展過(guò)程中液滴在孔口運(yùn)動(dòng)情況的影響，詳細(xì)分析了有限時(shí)間段內(nèi)孔口斷面處的壓強(qiáng)變化。研究結(jié)果表明，內(nèi)部壓應(yīng)力的脈動(dòng)是導(dǎo)致液滴被吸入孔內(nèi)并撞擊另一側(cè)孔壁形成飛濺現(xiàn)象的主要原因，模擬效果和實(shí)驗(yàn)結(jié)果符合良好。
　　在SPH模型中引入范德瓦爾斯?fàn)顟B(tài)方程中的吸引力來(lái)模擬液滴的表面張力。提出了一種改

3、進(jìn)的Quintic核函數(shù)，用于改善模型中排斥力的不穩(wěn)定問(wèn)題。模型中的吸引力采用與排斥力不同的插值核函數(shù)來(lái)處理，用于改善模型的拉力不穩(wěn)定問(wèn)題。通過(guò)模擬真空中液滴的拉伸振蕩過(guò)程，驗(yàn)證了本文計(jì)算模式在 SPH方法中模擬液滴的表面張力中的適用性。并且將模擬拉伸振蕩周期與理論值進(jìn)行對(duì)比，率定了該模式中范德瓦爾斯?fàn)顟B(tài)方程的拉應(yīng)力參數(shù)。并采用該模式模擬了考慮壁面黏性剪切力條件下液滴撞擊超疏水壁面過(guò)程，通過(guò)將模擬數(shù)據(jù)與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比，分析了壁面黏性剪

4、切力和表面張力對(duì)液滴運(yùn)動(dòng)過(guò)程的影響，研究表明：鋪展過(guò)程中液滴受到的壁面黏性剪切力起主要影響，黏性系數(shù)的大小隨著液滴運(yùn)動(dòng)速度增加而增大；回縮階段表面張力提供回縮動(dòng)力，該階段壁面摩擦阻力要小于鋪展階段。
　　提出了一種壁面可濕潤(rùn)性模擬機(jī)制，研究了壁面可濕潤(rùn)性對(duì)液滴靜態(tài)接觸角的影響，模擬了液滴撞擊可濕潤(rùn)壁面的過(guò)程，通過(guò)對(duì)比模擬數(shù)據(jù)和實(shí)驗(yàn)結(jié)果表明：液滴的靜態(tài)接觸角隨著壁面可濕潤(rùn)系數(shù)的增加而減小，壁面可濕潤(rùn)性對(duì)液滴的鋪展過(guò)程影響很??；回縮和

5、回彈階段壁面可濕潤(rùn)性的對(duì)液滴運(yùn)動(dòng)特征的影響較明顯；壁面可濕潤(rùn)系數(shù)會(huì)隨著壁面的飽和程度以及壁面上流體的脈動(dòng)壓強(qiáng)而改變，壁面飽和度越大壁面可濕潤(rùn)系數(shù)越小，流體的脈動(dòng)壓強(qiáng)越大可濕潤(rùn)系數(shù)越小。
　　針對(duì)自由表面流體與離散顆粒相互作用問(wèn)題提出一種新的SPH-DEM模型耦合方法，并且結(jié)合混合編譯原理建立了駕馭式可視化SPH-DEM模擬系統(tǒng)。本文提出的SPH-DEM模型耦合方法中流體粒子和離散顆粒之間的相互作用力包括三部分：排斥力、黏性剪切力和