模拟多球体在黏性流的动态运动.pdf

模擬多球體在黏性流的動態運動 1* 1 2 1 陳俊杉 、林欣瑞 、楊馥菱 、謝尚賢 1 2 國立臺灣大學土木工程學系、 國立臺灣大學機械工程學系 *E-mail: dchen@.tw 摘要 以電腦模擬固液二相顆粒流(solid-liquid two-phase granular flow)的運動行為， 有助於我們研究多顆粒流體中固體與液體在微觀尺度顆粒大小( )上的交互作用，以 及此微觀機制與所造成的巨觀現象 ( 多顆粒混合體的相互影響。本研究中整合了許) 多高效率的演算法，包括以晶格波茲曼法（Lattice Boltzmann Method, LBM ）求解 不可壓縮黏性流場，以離散元素法（Discrete Element Method, DEM ）同步求解沉浸 固體顆粒的移動以及接觸行為，並且以沉浸邊界法（ Immersed Boundary Method, IBM ）取代傳統邊界作用力的計算，求解固體跟液體的交互作用行為。藉由這樣的 整合系統，本研究試圖模擬土壤力學中的滲流砂湧 (Sand Boiling)現象，以期能經由 這樣的研究過程，引導出能有效分析及探討二相流相關工程問題與應用及其基礎機 制之物理。 關鍵字：固液二相流模擬、沉浸邊界法、晶格波茲曼法、離散元素法、砂湧 一、前言 近年來，隨著全球氣候的變遷而加劇之風災雨災及地震災害，對地處亞洲於太 平洋沿岸地震帶的台灣，許多可能被地震及暴雨誘發的自然災害，例如土石流、山 崩以及土壤液化等，發生益加頻繁，而且一直嚴重威脅著我們的日常生活。從這次 侵襲日本東北的地震可以看到，除了海嘯造成的洪水及大規模沖刷所造成的損傷外， 1 還有地震發生時，因為劇烈搖晃造成的土壤液化而導致地層下陷樓房橋梁倒塌等災 情，直接侵襲著人類的居住地，對我們的生命財產無疑地造成很大的威脅。 從物理上，我們可以將前述的這些自然災害流（natural hazardous flows ）簡化 為有強烈流固體交互作用的固液二相顆粒流（solid-liquid two-phase granular flow ）， 以方便但有效地探討顆粒和黏滯流體共同運動時，流體浮力及黏滯力等等對顆粒整 體運動的影響。由於兩者的交互作用非穩態且非線性，加上眾多球體的集合的行為 非等向性亦非均勻，因此使得固液二相顆粒流的整體現象相當複雜。在過去，要做 出有工程意義且包含複雜邊界的三維流場數值模擬一直是難以達成的事情。除了數 值演算法的成熟度外，最大的瓶頸為電腦硬體效能，大大地限制了人類以數值模擬 自然與工程問題的尺度。近年來隨著科技的演進，電腦的計算效能大幅提昇，以往 被認為難以達成的三維複雜邊界模擬，已經是變成相當普遍的研究，而具有工程意 義及研究方法挑戰的固液二相顆粒流模擬，已經漸漸盛行於各種顆粒流研究中。 傳統上，求解固液二相流問題，最大的問題是求解包含複雜且快速變化邊界的 流場。統御牛頓流體運動的 Navier-Stokes 方程式，求解時需在每個時距計算流場 速度以及壓力，但這兩個未知數耦合在一起，也因此在求解上本身就有一定的難度； 加上求解 N-S 方程式中非線性的對流加速度項次使得本偏微分方程式在精度及穩 定性兼顧要求下的演算法相當複雜，若再引入形狀複雜的邊界，會讓求解更加困難， 遑論在顆粒流中大量且快速變化的固體邊界。為了有效率求解包含複雜邊界的流場 問題，本研究導入沉浸邊界，此方法最初由 Peskin在 1982 年提出[1] ，為計算血球 在邊界複雜的心臟血管內移動時流場。藉由晶格波茲曼法求解流場並與離散元素法 同步模擬多固體顆粒在黏性流體中的運動，並將這套計算工具應用在砂湧現象的模 擬，試圖展示這套工具的實用性以及計算效率。 除演算法整合之外，與計算效率及穩定度密不可分的議題是時距的選用。純固 體顆粒碰撞計算時