固液二相流計算平台開發與應用

Development and application of computation platform for solid-liquid flows

莊益彰(Yi-Chang Chuang)

指導教授 ： 陳國慶

共同指導教授 ： 周逸儒(Yi-Ju Chou)

國立臺灣大學/工學院/應用力學研究所/碩士(2012年)

https://doi.org/10.6342/NTU.2012.10252

摘要

隨著電腦科技的演進，固液二相流模擬日漸受到重視。本研究以電腦數值方法進行固液二相流計算平台開發與實際應用，利用CFD-ACE計算流體力學軟體以有限體積法求解Navier-Stokes方程式計算流體，採用的模型為Flow Module；固體方面，分別採用表面質點近似法與拉格朗日法，模型分別為Macro Particle Module與Spray Module，因此在求解問題時有兩種耦合方式可求解固體與液體間交互作用行為。前者，在顆粒碰撞時壓力場將會發散，在研究中證實利用自行發展的碰撞模型可以避開此問題使數值穩定讓模擬順利運作，使用此種耦合方式計算效率較低但可以微觀角度探討二相流問題；後者，用來實際模擬微米尺度顆粒在微流道中聚集效應，以往認為此效應主因剪力梯度與壁面效應引起的側向力造成顆粒側向位移至平衡位置所造成，壁面效應與剪力梯度都與顆粒體積造成周圍流場不對稱有關，但在Spray Module中並無考慮顆粒體積，也無法解析出顆粒周圍的流場卻能模擬出與實驗相同的趨勢，並根據數值實驗的結果發現，在微米尺度下，邊牆效應以及剪力梯度的效應並不明顯，故相較之下流場才是主導顆粒側向位移的主因，使用此種耦合方式計算效率較高但只能觀察系統整體巨觀行為。

關鍵字

固液二相流模擬 ； Navier-Stokes方程式 ； 表面質點近似法 ； 拉格朗日法 ； 顆粒聚集 ； 剪力梯度 ； 壁面效應

並列摘要

With the advance of the computer technology, simulation of the two-phase for the solid-liquid system has been emphasized. In this research, we develop a computation platform for the solid-liquid system, which utilizes the flow module of the CFD-ACE software to solve the Navier-Stokes equations. There are two coupling methods to solve the interactions between the solid phase and the fluid phase. To solve the motion of solids, we use the surface maker point approach and the Lagrangian particle-tracking method.
In the former method, the pressure field blows up when two particles collide. To resolve this issue, we develop a collision model to calculate the two-phase flow motion that resolves individual micro particles. We use the latter method to simulate particle focusing in micro channels. Previous literatures pointed out that the main factors of lateral forcing to move the micro particles to their equilibrium positions are shear gradient and wall effect. These two forcing factors account for the asymmetry of the flow field surrounding the particle. However, while the spray module does not consider the particle volume and the surrounding flow field, simulation results qualitatively agree with the experimental observations. Moreover, according to a separate numerical experiment, perturbations of velocity and pressure field due to the presence of a single particle are not obvious. We find that the background hydrodynamic forcing is the dominant cause of the lateral displacement. In conclusion, in simulating particle focusing in micro channels, flow perturbation due to the pressure of individual particles may not be important, and the spray module can obtains the relevant approximation. However, the spray module can only be used to observe the macroscopic behavior.

並列關鍵字

simulation of the two-phase for the solid-liquid system ； Navier-Stokes equation ； surface maker point approach method ； Lagrangian particle-tracking method ； particle focusing ； shear gradient ； wall effect

參考文獻

[38] 黃昱傑.(2011). 「顆粒材料二維及三維剪切帶之識別」：國立台灣大學應用力學研究所碩士論文。

[39] 陳尹中.(2011). 「旋轉流體中泰勒渦柱與池盆渦漩交互影響之數值計算分析」：國立台灣大學應用力學研究所碩士論文。

[40] 鄭宇軒.(2011). 「垂直風力發電機之模擬與扭力及效率分析」：國立台灣大學應用力學研究所碩士論文。

[41] 蔡倆俠.(2011). 「以數值方法探討壁面效應對粒子沉降行為之影響」：國立台灣大學應用力學研究所碩士論文。

[42] 林欣瑞.(2011). 「多球體在黏性流的動態運動模擬」：國立台灣大學土木工程研究所碩士論文。