請用此 Handle URI 來引用此文件:
http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/62218
完整後設資料紀錄
DC 欄位 | 值 | 語言 |
---|---|---|
dc.contributor.advisor | 楊照彥(Jaw-Yen Yang) | |
dc.contributor.author | Yun-Da Tsai | en |
dc.contributor.author | 蔡昀達 | zh_TW |
dc.date.accessioned | 2021-06-16T13:34:31Z | - |
dc.date.available | 2013-08-14 | |
dc.date.copyright | 2013-08-14 | |
dc.date.issued | 2013 | |
dc.date.submitted | 2013-07-18 | |
dc.identifier.citation | [1]P. Degond, S. Jin, L. Mieussens, A smooth transition model between kinetic and hydrodynamic equation, Journal of Computational Physics, 209 (2005) 665-694.
[2]P. Degond, G. Dimarco, L. Mieussens, A moving interface method for dynamic kinetic-fluid coupling, Journal of Computational Physics, 227 (2007) 1176-1208. [3]J.-F. Bourgat, P. Le Tallec, MD. Tidriri, Coupling Boltzmann and Navier-Stokes equations by friction, Journal of Computational Physics, 127 (1996) 227-245. [4]Klar, H. Neunzert, J. Struckmeier, Transition from kinetic theory to macroscopic fluid equations: a problem for domain decomposition and a source for new algorithms, Transport Theory Stat. Phys. 29 (1–2) (2000) 93–106. [5]J. Schneider, Direct coupling of fluid and kinetic equations, Transport Theory Stat. Phys. 25 (6) (1996) 681–698. [6]J.-F. Bourgat, P. LeTallec, B. Perthame, Y. Qiu, Coupling Boltzmann and Euler equations without overlapping, in: Domain decomposition methods in science and engineering, Contemporary Mathematics, vol. 157, AMS, Providence, RI, 1994, pp. 377–398. [7]P. Degond, J.G. Liu, L. Mieussens, Macroscopic fluid model with localized kinetic upscaling effects, MMS 5 (3) (2006) 940–979. [8]R.E. Caflisch, Monte Carlo and quasi-Monte Carlo Methods, Acta Numer. (1998) 1–49. [9]Muljadi, B. P., Yang, J. Y., A direct Boltzmann-BGK equation solver for arbitrary statistics using the conservation element/solution element and discrete ordinate method, Kuzmin, A. (ed.) Computational Fluid Dynamics 2010, pp. 637–642. [10]Muljadi, B. P., Yang, J. Y., Simulation of Shock Wave Diffraction over 90° Sharp Corner in Gases of Arbitrary Statistics, Journal of Statistical Physics, 145, pp. 1674-1688. [11]Shi, Y. H., Huang, J. C. & Yang, J. Y., High resolution kinetic beam schemes in generalized coordinates for ideal quantum gas dynamics, Journal of Computational Physics, Vol. 222, 2007 573–591. (SCI) [12]Shu, C. W. & Osher, S., Efficient Implementation of Non-oscillatory Shock Capturing Schemes II, Journal of Computational Physics, 83, pp. 32-78. [13]Uehling, E. A. & Uhlenbeck, G. E., Transport phenomena in einstein-bose and fermi-dirac gases, I, Physical Review, 43, pp. 553-561. [14]Van Leer, B., Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov’s method, Journal of Computational Physics, 32, pp. 101–136 . [15]Yang, J. Y., Hsieh, T. Y. & Shi, Y. H., Kinetic Flux Vector Splitting Schemes for Ideal Quantum Gas Dynamics, SIAM Journal on Scientific Computing, 29, pp. 221-244. [16]Yang, J. Y., Hsieh, T. Y., Shi, Y. H. & Xu, K., High Order Kinetic Flux Vector Splitting Schemes in General Coordinates for Ideal Quantum Gas Dynamics, Journal of Computational Physics, 227, pp. 967-982. [17]J. C. Huang, T. Y. Hsieh, J. Y. Yang and K. Takayama, A Numerical Study of Oblique Shock Wave Reflections over Wedges in a Quantum Gas, Shock Wave Journal, Vol. 18, No. 3, Aug. 2008 193–207. (SCI) [18]謝澤揚 (2007) 聲子熱傳輸與理想量子氣體動力學之高解析算則,國立台灣大學工學院應用力學所博士論文,台北。 [19]Yang, J. Y., Huang, J. C, Rarefied Flow Computations Using Nonlinear Model Boltzmann Equations, Journal of Computational Physics, Vol.120, Mar. 1995, pp.323-339. (SCI) [20]M.A. Katsoulakis, A.J. Majda, A. Sopasakis, Multiscale couplings in prototype hybrid deterministic/stochastic systems: part i, deterministic closures, Comm. Math. Sci. 2 (2) (2004) 255–294. [21]L. Pareschi, Hybrid multiscale methods for hyperbolic and kinetic problems, Esaim Proceedings, vol. 15, EDP Sciences, SMAI, 2005,2005. [22]W.-L. Wang, I.D. Boyd, Predicting continuum breakdown in hypersonic viscous flows, Physics of Fluids, 15 (2003) 91–100. [23]S. Tiwari, Coupling of the Boltzmann and Euler equations with automatic domain decomposition, Journal of Computational Physics, 144 (1998) 710–726. [24]M.A. Katsoulakis, A.J. Majda, A. Sopasakis, Multiscale couplings in prototype hybrid deterministic/stochastic systems: part i, deterministic closures, Comm. Math. Sci. 2 (2) (2004) 255–294. [25]B. Weinan E, Engquist, Multiscale Modeling and computation, Notices AMS 50 (9) (2003) 1062–1070. [26]B. Weinan E, Engquist, The heterogeneous multiscale methods Comm. Math. Sci. 1 (2003) 87–133. [27]C. Cercignani, R. Illner, M. Pulvirenti, The Mathematical Theory of Dilute Gases, Springer-Verlag, New York, 1995. [28]D. Levermore, W.J. Morokoff, B.T. Nadiga, Moment realizability and the validity of the Navier Stokes equations for rarefied gas dynamics Physics of Fluids 10 (12) (1998). [29]V.I. Kolobov, R.R. Arslanbekov, V.V. Aristov, A.A. Frolova, S.A. Zabelok, Unified solver for rarefied and continuum flows with adaptive mesh and algorithm refinement, J Journal of Computational Physics, 223 (2) (2007) 589–60 | |
dc.identifier.uri | http://tdr.lib.ntu.edu.tw/jspui/handle/123456789/62218 | - |
dc.description.abstract | 針對稀薄流問題求解時,通常選擇波茲曼BGK模型來求解。但當BGK的鬆弛時間很小時,也就是當碰撞頻率變成很大時將會造成BGK方程的計算成本非常巨大。因此本文主要在研究當流場遇到此類型問題時,能透過適當的機制將巨觀Euler與微觀的BGK做轉換以達到節省計算成本的目的。
本文利用直接解數值方法來求解半古典波茲曼方程與Euler方程組成的耦合系統以模擬稀薄流。數值方法方面使用離散座標法(D.O.M.)與迎風法(Upwind),並加入高解析算則如全變量消逝法(Total Variation Diminishing, TVD)。 在傳統耦合模型求解稀薄流時。常使用介面邊界條件來轉換。但此方法必須添加額外的邊界條件,在維度拓展到二維時會非常複雜。因此耦合方法選擇加入緩衝區的設定,透過此緩衝區的設定能平滑的將巨觀與微觀尺度做連結。尺度轉換的關鍵在於緩衝區的設計,因此本文也針對三種不同類型的緩衝區做尺度轉換的比較。三種緩衝區分別是線性緩衝區、餘弦緩衝區、雙曲正切函數的緩衝區。 做耦合模型時需要求兩種分布函數。一種為BGK系統下的分布函數,另一種則為巨觀尺度下所求得的分布函數-平衡態分布函數。透過此兩種分布函數相加可得到整體系統的分布函數,進而得到系統的巨觀物理量。使用耦合尺度來求解震波管問題,可成功透過緩衝區捕捉震波、膨脹波。然而,計 顯示緩衝區大小若低於無因次長度0.25時模擬效果較容易出現不守恆量,而當緩衝設定位置遠離震波時三種型態的緩衝區就無明顯分別。 | zh_TW |
dc.description.abstract | For targeting flows problem with different degrees or regimes of rarefaction, we can rely on Boltzmann-BGK model to find tractable method to solve this kind of problems [15] [17] [19]. But it is also well known that when collision source term becomes ‘large’ , i.e. the collisions become more frequent, however, its computational cost to use BGK model will increase. Hence, the primary goal of this work is to introduce an appropriate strategy to combine an Euler solver and a BGK solver and achieve cost reduction. In this work we follow and look for ways to improve upon ideas of previous approaches [1] [2]. We also extend the method to semiclassical Boltzmann-BGK equation. The numerical method is based on discrete ordinate method to approach the recovery of macroscopic quantities in velocity space, and upwind method, high-resolution methods such as TVD, to evolve information in physical space.
Traditional coupled methods [3], rely on boundary conditions for the coupling two different hydrodynamic solvers, This strategy becomes more complicated when space domain expands to 2-dimensional problems. In this work, the idea is to use a transition region called buffer zone, in which, both semiclassical hydrodynamic and kinetic equations will be solved simultaneously. A buffer zone using three types of cut functions was designed and tested, namely a linear transition, trigonometric transition and hyper-tangent transition function were used. For our coupling method, two distribution functions are employed. Through the sum of both two distribution functions, we can recover total distribution function for the whole system and therefore recovered physical quantities. Applying the present couple method to solve sod’s problem, it is shown that shock wave and expansion wave can all be captured successfully. However, our results also show that there is a minimum dimensionless length of 0.25 for the buffer zone to work properly. If the buffer zone is far away from the shock region, results from the three transition functions show no differences. | en |
dc.description.provenance | Made available in DSpace on 2021-06-16T13:34:31Z (GMT). No. of bitstreams: 1 ntu-102-R00543077-1.pdf: 8261260 bytes, checksum: 4abd2cba3accc75134f72d69fedc5bb6 (MD5) Previous issue date: 2013 | en |
dc.description.tableofcontents | 口試委員會審定書 #
誌謝 1 中文摘要 2 ABSTRACT 3 大綱 5 圖附錄 7 第一章 緒論 11 1.1 引言 11 1.2 Riemann Problem簡介 12 1.3 文獻回顧 13 1.4 本文架構 14 第二章 波茲曼方程 15 2.1 氣體動力論 15 2.2 Liouville方程 16 2.3 Boltzmann方程 17 2.4 鬆弛時間近似 19 2.5 連續體模型方程 20 2.6 緩衝區與截斷函數 22 2.7 尺度耦合法 23 2.7.1 Kinetic/Kinetic Coupling 25 2.7.2 Kinetic/Hydrodynamic Coupling 26 第三章 半古典Boltzmann方程 28 3.1 三種統計 28 3.2 半古典Boltzmann-BGK方程 28 3.3 半古典粒子運動模型 31 第四章 數值方法 34 4.1 離散座標法 34 4.2 空間離散 34 4.2.1 迎風算則(Upwind Scheme) 36 4.2.2 高解析算則(High Resolution Schemes) 37 4.3 時間離散 38 4.4 初始和邊界條件 39 4.5 無因次化控制方程與離散方程 41 4.5.1 Kinetic/kinetic coupling 41 4.5.2 Kinetic/Hydrodynamic Coupling 42 4.5.3 Roe-Euler解法 43 第五章 數值模擬結果與討論 47 第六章 結論與展望 76 6.1 結論 76 6.2 展望 77 參考文獻 78 | |
dc.language.iso | zh-TW | |
dc.title | 半古典波茲曼方程與尤拉方程耦合之稀薄流模擬 | zh_TW |
dc.title | A Semi-Classical Kinetic and Hydrodynamic Coupling Method for Rarefied Flow Computation | en |
dc.type | Thesis | |
dc.date.schoolyear | 101-2 | |
dc.description.degree | 碩士 | |
dc.contributor.oralexamcommittee | 黃俊誠,陳旻宏,謝澤揚,楊世昌 | |
dc.subject.keyword | 尤拉方程,半古典波茲曼BGK,緩衝區,耦合法,離散座標法,全變量消逝法, | zh_TW |
dc.subject.keyword | Euler,Semi-Classical Boltzmann BGK,Buffer Zone,Couple Method,Discrete Ordinate Method,Total Variation Diminishing, | en |
dc.relation.page | 80 | |
dc.rights.note | 有償授權 | |
dc.date.accepted | 2013-07-18 | |
dc.contributor.author-college | 工學院 | zh_TW |
dc.contributor.author-dept | 應用力學研究所 | zh_TW |
顯示於系所單位: | 應用力學研究所 |
文件中的檔案:
檔案 | 大小 | 格式 | |
---|---|---|---|
ntu-102-1.pdf 目前未授權公開取用 | 8.07 MB | Adobe PDF |
系統中的文件,除了特別指名其著作權條款之外,均受到著作權保護,並且保留所有的權利。