指導教授：楊照彥臺灣大學：應用力學研究所蔡承翰Tsai, Cheng-HanCheng-HanTsai2014-11-302018-06-292014-11-302018-06-292014http://ntur.lib.ntu.edu.tw//handle/246246/264084在連體模型中，主要的統御方程式是透過歐拉或納維爾方程式來求解。然而在氣體稀薄程度提升時，將不再適用，此時方程式為波茲曼方程式。在相關的研究中，對波茲曼方程式作空間及速度離散，即可推導出格子波茲曼法 (Lattice Boltzmann Method, LBM)。 在本文中，根據Uehling-Uhlenbeck Boltzmann-BGK方程（Uehling-Uhlenbeck Boltzmann Bhatnagar-Gross-Krook Equation）及橢圓統計BGK方程（Ellipsoidal Statistical BGK equation, ESBGK）所發展出的半古典橢圓統計格子波茲曼法，透過這個方法，我們可以有效利用其主導的分布函數來計算出所需之宏觀量。利用此法以及D2Q9速度格子模型為基礎，模擬在三種不同粒子統計下，流體流過使用沉浸邊界速度修正法(Immersed Boundary Velocity Correction Method, IBVCM)的圓柱所產生的渦漩、渦度、壓力以及阻力係數之比較。由模擬結果顯示，ESBGK模型之間三種不同的粒子統計差異(Bose-Einstein統計、Fermi-Dirac統計和Maxwell-Boltzmann統計)，並且隨著雷諾數的增加，流線出現了馮卡門漩渦現象。除此之外，還可運用半古典格子波茲曼可還原成古典格子波茲曼法的特性(MB統計)，與其相關研究作驗證。當圓柱在流道中持續的移動，此時模擬比較而觀察出尾流渦漩在三種統計中有明顯不同的差異結果。Modeling gases in the continuum level is traditionally achieved with macroscopic level by using Euler or Navier-Stokes equations. However, as the degree of rarefaction of a gas increses, the governing equation becomes Boltzmann Equation. The Lattice Boltzmann method is derived by discretizing Boltzmann equation in physical and velocity space. In the study, the development of a semiclassical lattice Boltzmann–Ellipsoidal Statistical method is based on the Uehling-Uhlenbeck Boltzmann-BGK equation. According to the method, we can effectively link its dominant distribution function to calculate the quantities of macroscopic properties. Here, we present simulations of the flow over cylinder for several Reynolds numbers based on D2Q9 lattice model and the semiclassical lattice Boltzmann–Ellipsoidal Statistical method. In this work, the Immerse Boundary Velocity Correction method (IBVCM) has been used to model the boundary of the cylinder. We compare the results of vortices, pressure, drag coefficient for different particle statistics : Bose-Einstein, Fermi-Dirac and Maxwell-Boltzmann statistics. By studying the streamlines, we observed von Karman vortex street phenomenon as the Reynolds number increases. In addition, the movement of the cylinder boundary is not stationary in the flow channel by taking advantage of the boundary using IBVCM. We observed the wake and found differences in the results of the three statistics.致謝……………………………………………………………………………………I 中文摘要…………………………………………………………………………II Abstract……………………………………………………………………III 目錄…………………………………………………………………………………IV 圖目錄………………………………………………………………………………VI 符號表………………………………………………………………………………IX 第一章 緒論 1 1-1 計算流體力學 1 1-2 格子Boltzmann法簡介 2 1-3 量子理論簡介 3 1-4 文獻回顧 4 1-5 本文目的 6 第二章 波茲曼傳輸方程式 7 2-1 氣體運動理論 7 2-2 分布函數 9 2-3 波茲曼傳輸方程式 11 2-4 Boltzmann H定理 15 2-5 馬克斯威爾分布 18 2-6 Boltzmann BGK 方程 21 2-7 格子Boltzmann方程 22 2-8 單鬆弛時間離散速度模型 24 2-9 Hermite展開平衡態分布函數 26 第三章 半古典格子波茲曼方程式 30 3-1 量子統計 30 3-2 半古典格子Boltzmann-BGK方程式 32 3-3 單鬆弛時間Chapman-Enskog分析 38 第四章 半古典橢圓統計格子波茲曼方程式 42 4-1 半古典格子Boltzmann-ESBGK方程式 42 4-2 宏觀物理量之求法 47 4-3 格子Boltzmann-ESBGK方程之Chapman-Enskog分析 47 第五章 基本模型與邊界處理方法 52 5-1 格子Boltzmann-ESBGK方程式 52 5-2 邊界條件 53 5-2-1 沉浸邊界速度修正法 53 5-2-2 標準反彈邊界 55 5-3 模擬參數定義及收斂條件 56 第六章 模擬結果與討論 58 6-1 圓柱繞流(Flow Over Cylinder) 58 6-2 模擬問題描述 59 6-3 模擬結果分析與探討 61 第七章 結論與未來展望 99 7-1 結論 99 7-2 未來展望 100 文獻回顧 10115148023 bytesapplication/pdf論文公開時間：2018/08/31論文使用權限：同意有償授權(權利金給回饋本人)半古典格子波茲曼方法半古典橢圓統計格子波茲曼法D2Q9模型圓柱繞流沉浸邊界速度修正法thesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/264084/1/ntu-103-R01543013-1.pdf