應用圖形處理器叢集與晶格波茲曼法模擬二相流問題

Lattice Boltzmann simulation of two-phase flow on multi-GPU cluster

林方婷(Fang-Ting Lin)

指導教授 ： 林昭安

國立清華大學/工學院/動力機械工程學系/碩士(2019年)

摘要

本研究以D3Q19晶格波茲曼法模擬二相流問題，主要比較兩種介面捕捉公式（Cahn-Hilliard equation 與 Allen-Cahn equation）之差異。以靜止液珠的測試為例，兩種介面捕捉公式造成的假性速度相差不大，但Cahn-Hilliard模型模擬的結果會出現氣液組成率超出上限（純液體）的情形，進而導致液珠收縮的現象；單顆液珠震盪與兩顆液珠融合的測試結果顯示兩種介面捕捉公式的模擬結果與解析解一致。不過Allen-Cahn模型無法模擬旋節分解（Spinodal decomposition），只能用Cahn-Hilliard模型表現隨機生成的氣液組成率到最後發展分離成氣液雙相。本研究利用圖形顯示卡叢集搭配一維切割進行平行運算，並使用多個stream將邊界點的傳值時間隱藏在內部點的計算時間中進行加速計算。

關鍵字

晶格波茲曼法 ； Allen-Cahn模型 ； 旋節分解 ； 圖形顯示卡

並列摘要

This thesis presents the simulation of the three-dimensional two-phase lattice Boltzmann model on the graphics processing unit (GPU) cluster. In this thesis, the objective is to develop a program according to Lee's model with the Allen-Cahn equation, then comparing the results to Lee's model with the Cahn-Hilliard equation. Take one stationary drop as an example, the difference of the generated spurious velocity between these two interface capturing equations are not much different. However, the drop shrinkage phenomenon only occurs in the result of the Cahn-Hilliard model. The numerical results of the Allen-Cahn model both in single drop oscillation and two drops merging are identical to the theoretical solution. However, the Allen-Cahn model is unable to simulate the spinodal decomposition phenomenon, so the phase decomposition from the initial random liquid-vapor composition only can be shown in the result of the Cahn-Hilliard model. This thesis uses one-dimensional decomposition to implement parallel programming and uses the overlapping strategy which hiding the boundary exchange time in the computing time to improve computational efficiency.

並列關鍵字

LatticeBoltzmannmethod ； Allen-Cahnequation ； Spinodal decomposition ； GPU

參考文獻

1. T. Y. Hou, J. S. Lowengrub and M. J. Shelley, ``Boundary integral methods for multicomponent fluids and multiphase materials," J. Comput. Phys.169, 302-362, (2001).

2. S. Osher and R. P. Fedkiw, ``Level set method: An overview and some recent results," J. Comput. Phys. 169, 463-502, (2001).

3. R. Scardovelli and S. Zaleski, ``Direct numerical simulation of free-surface and interfacial flow," Annu. Rev. Fluid Mech. 31, 567-603, (1999).

4. U. Frisch, B. Hasslacher and Y. Pomeau, ``Lattice-gas automata for the Navier-Stokes equation," Phys. Rev. Lett. 56, 1505-1508, (1986).

5. S. Wolfram, ``Cellular automaton fluids 1: Basic theory," J. Stat. Phys. 45, 471-526, (1986).