三維多鬆弛時間半古典格子波茲曼法之流場模擬

A Three-Dimensional Semiclassical Lattice Boltzmann Method Using Multiple Relaxation Times For Quantum Hydrodynamic Flow Simulations

黃柏堯(Bo-Yao Huang)

指導教授 ： 楊照彥

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

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

摘要

發展基於Uehling-Uhlenbeck Boltzmann-BGK(Uehling-Uhlenbeck Boltzmann Bhatnagar-Gross-Krook Equation)方程的三維多鬆弛時間半古典格子Boltzmann法。此種方法藉由Hermite展開法得到離散半古典平衡態分布函數的Hermite展開形式，並且引入多鬆弛時間的概念，推導出三維半古典平衡態分布函數的矩空間形式，最後藉由此半古典矩空間平衡態分布函數加以計算便可得到各種巨觀物理量。
以D3Q19速度模型模擬方腔流，使用單一Re數，比較三維多鬆弛時間半古典格子Boltzmann法和三維單鬆弛時間半古典格子Boltzmann法的結果後，可以發現多鬆弛時間模擬的流場結構更為完整；並且為了更加了解多鬆弛時間半古典格子Boltzmann法的可行性，以二維D2Q9速度模型模擬長寬比例為1:2與2:1長方腔流流場問題，經由幾種Re數與三種遵循Bose-Einstein統計、Fermi-Dirac統計與Maxwell-Boltzmann統計的粒子模擬下，由結果分析可以明確展示下游次渦流在不同量子效應下產生不同的變化形態；最後經由各種邊界條件的計算後確定了此種方法的穩定性及準確性。

關鍵字

多鬆弛時間 ； D2Q9 ； D3Q19速度模型 ； 方腔流 ； 格子波茲曼法 ； 半古典格子波茲曼法

並列摘要

According to two dimensional multiple-relaxation-time semiclassical lattice Boltzmann method, a three dimensional multiple-relaxation-time semiclassical lattice Boltzmann based on the Uehling-Uhlenbeck Boltzmann-BGK equation has been successfully developed by using D3Q19 lattice model. The method is derived by expanding the equilibrium distribution function in term of Hermite polynomials. Then, we combine the concept of multiple-relaxation-time with lattice Boltzmann method to obtain the moment matrix and various physical quantities.
Simulations with SRT and MRT for the lid-driven cavity flows based on D3Q19 lattice model have been carried out for three quantum particles including Bose-Einstein , Fermi-Dirac and Maxwell-Boltzmann statistics. Simulations show the streamline and pressure contour of 3D lid-driven cubic cavity flows. The result can be found in the multi-relaxation time simulation of the flow field structure is more complete. Moreover, comparing with different length and width of cavity flows such as 1:2 and 2:1, it is easy to notice that the shape of downstream eddy of rectangle flows are different for three quantum particles. Finally, after testing various boundary conditions, the stability and the accuracy of this method is illustrated.

並列關鍵字

Mutiple-relaxation-time ； D3Q19 lattice model ； Cavity flows ； Semiclassical lattice Boltzmann method ； Lattice Boltzmann method

參考文獻

[15] 洪立昕 半古典晶格波滋曼方法，國立台灣大學工學院應用力學所博士論文，台北（2011）。

[2] Yang, J. Y, & Hung, L.H., “Lattice Uehling-Uhlenbeck Boltzmann Bhatnagar- Gross-Krook Hydrodynamics of Quantum Gases.” Phys. Rev. E, vol. 79, p.056708, (2009).

Phil. Trans. R. Soc., vol. 360, pp.437–451, (2002).

[4] 郭照立、鄭楚光 格子Boltzmann方法的原理及應用（Theory and Applications of Lattice Boltzmann Method），科學出版社（2009）。

[5] Chen, S.Y., Martinez, D. & Mei, R.W. “On Boundary Conditions in Lattice Boltzmann Methods.” Physics of Fluids, vol. 8, pp.2257–2536, (1996).