進階搜尋


下載電子全文  
系統識別號 U0026-3107201911035900
論文名稱(中文) 一圓盤飄落之流場與其周圍壓力震盪之研究
論文名稱(英文) Analysis of flow field and Pressure fluctuation for a falling disk.
校院名稱 成功大學
系所名稱(中) 航空太空工程學系
系所名稱(英) Department of Aeronautics & Astronautics
學年度 107
學期 2
出版年 108
研究生(中文) 湯發民
研究生(英文) Fa-Ming Tang
學號 P46061403
學位類別 碩士
語文別 中文
論文頁數 83頁
口試委員 指導教授-林三益
口試委員-陳文立
口試委員-吳村木
口試委員-胡舉軍
中文關鍵字 六自由度  重疊網格  快速複利葉轉換  動態網格  DNS  自由落體  三維不可壓縮奈維爾-史托克方程式 
英文關鍵字 Computation Fluid dynamics  six-degree of freedom  overset mesh  Direct Numerical Simulation 
學科別分類
中文摘要 本研究旨在以計算流體力學模擬一剛體圓盤在特定之無因次參數下模擬自由落體的六自由度運動行為與流場分析,並取周圍流場監測點之壓力,觀察其振盪情形並分析其週期性與運動行為關係。
流場計算使用商業計算流體力學軟體Ansys Fluent,用二階準確度直接求解三維非穩態不可壓縮奈維爾-史托克方程式。網格用Ansys ICEM 做為網格劃分之工具,分別產生圓盤周圍之結構性主網格與運動場域之結構動態網格,在分別匯入Ansys Fluent中做重疊網格的處理。
在物理模型上以穩態與非穩態的中低雷諾數自由流流經一固定圓盤,分析其所受之升阻力與流場做為靜態上的驗證;動態上模擬一顆圓球墜落於不同密度與黏滯係數之流體,與文獻比對其位移時間圖與速度對時間圖做為動態上的驗證。
本研究探討的無因次物理參數有無因次轉動慣量、阿基米德數及平均雷諾數,在三組不同的物理參數下,圓盤分別以固定週期與位移來回擺盪的簡協飄移(fluttering)、不停的項同一水平方向做翻滾(tumbling)與繞一軸心做打轉的螺旋(spiral)的運動行態飄落。本研究觀察其在三種不同運動行態下的運動與受力情形,並取樣周圍流場若干點作為壓力監測點,以快速傅立葉轉換分析壓力震盪頻率。
英文摘要 The vortex structure, pressure fluctuation, and three kinds of behavior of motion for a free falling disk are simulated in this study. The three kinds of behaviors of motion are fluttering, tumbling, and spiral. There are three dimensionless parameters are chosen to produce the kinds of behavior of motion, which are dimensionless moment of inertia I^*, Archimedes Number Ar, and Reynolds Number Re, will determine the behavior of motion in this study
For each case, we analysis the periodicities about the disk’s motions, the forces acted on the disk, and the pressure fluctuation monitored at some points in the middle field. By analyzing the FFT spectrogram, the results show that vertical disk motion has same period as the pressure fluctuation. We also found that the vortex dominates the pressure distribution in the flow field as well as on the disk.
論文目次 摘要 I
誌謝 IV
符號說明 IX
第一章 緒論 1
1-1. 內容大綱 1
1-2 研究目的與動機 2
1-3文獻回顧 3
第二章 統御方程式與數值方法 5
2-1 無因次化物裡參數 5
2-2 統御方程式 6
2-3 六自由度求解 7
2-4 壓力之頻譜 8
2-5 網格生成與重疊網格 9
第三章 物理模行驗證與網格獨立性分析 11
3-1 低雷諾數均勻流流經零攻角之三維圓盤 11
3-2 於雷諾數500下流經40度攻角之三維圓盤 12
3-3 圓球掉落之動態驗證 13
3-4 穩態低雷諾數網格獨立性分析 14
3-5 動態之網格獨立性分析 15
第四章 結果與討論 16
4-1簡協飄移(Fluttering) 16
4-1-1 運動情形與壓力震盪週期分析 17
4-1-2 流場分析 18
4-2 翻滾(Tumbling) 19
4-2-1運動情形與壓力震盪分析 20
4-2-2 流場分析 21
4-3 螺旋(Spiral) 22
4-3-1運動情形與壓力震盪分析 23
4-3-2 流場分析 24
第五章 結論與建議 26
5-1結論 26
5-2 未來工作與建議 28
參考文獻 29
表[1]升阻力與文獻[6]比對表 32
表[2]粗、中、細網格設定 32
表[3]各飄落行為相對應的無因次參數 32
圖[2-1]ϕ、θ、ψ所代表之相對應角度 33
圖[2-2](a)重疊處理前 (b)重疊處理後 33
圖[3-1]驗證3-1之圓柱流場 34
圖[3-2]3-1網格切面1 34
圖[3-3] 3-1網格切面2 35
圖[3-4]左右分別為雷諾數10、100,上排圖片為模擬之渦流剖面圖,中排為流線圖,與下排文獻[4]的圖片比對流場現象。 35
圖[3-5]雷諾數為10之渦流 36
圖[3-6]雷諾數140時之渦流情形,上為模擬結果,下為文獻[4]結果 36
圖[3-7]雷諾數180時之渦流情形,上為模擬結果,下為文獻[4]結果 37
圖[3-8](a)雷諾數300時(b)雷諾數2250之渦流 37
圖[3-9]阻力係數比較圖 38
圖[3-10]2-1-2計算之圓柱流場場域 38
圖[3-11]3-2網格切面1 39
圖[3-12]上圖為文獻[4]渦流結構,下圖為本研究模擬比對圖 39
圖[3-13]壓力係數與文獻[6]比對 40
圖[3-14]3-3計算域圖 40
圖[3-15]重疊後的網格剖面 41
圖[3-16]動力黏度為0.212 kg/m‧s時的位置與時間圖 41
圖[3-17]動力黏度為0.212 kg/m‧s時的速度與時間圖 42
圖[3-18]動力黏度為0.113 kg/m‧s時的位置與時間圖 42
圖[3-19]動力黏度為0.113 kg/m‧s時的速度與時間圖 43
圖[3-21] 3-5計算域示圖 44
圖[3-22] X與Y平面上軌跡圖 44
圖[3-23]定向角ψ對時間變化圖 45
圖[4-1]X-Y平面運動軌跡 45
圖[4-2]每25時間步所記錄之圓盤圖 46
圖[4-3]4-1計算場域 ,l2、l3、b2、b3、r2、r3為欲觀察之監測點 46
圖[4-4]定向角ψ定義示圖 47
圖[4-5] ψ與速度隨時間變化圖 47
圖[4-7] ψ、X方向受力與速度時間變化圖 48
圖[4-8] ψ、Y方向受力與速度時間變化圖 48
圖[4-9]X方向受力與Z方向力矩隨時間變化圖 49
圖[4-10]速度與受力頻譜圖 49
圖[4-11]b2點壓力與圓盤受力隨時間變化圖 50
圖[4-12]監測點l2壓力與圓盤受力隨時間變化圖 50
圖[4-13]r2點壓力與圓盤受力隨時間變化圖 51
圖[4-14] b3點壓力與圓盤受力隨時間變化圖 51
圖[4-15]l3點壓力與圓盤受力隨時間變化圖 52
圖[4-16]r3點壓力與圓盤受力隨時間變化圖 52
圖[4-17]b2、l2、r2與受力頻譜比對圖 53
圖[4-18] b3、l3、r3與受力頻譜比對圖 53
圖[4-19]路徑上脫離之渦流與壓力分佈(XY視角) 54
圖[4-20] 路徑上脫離之渦流與壓力分佈(ZY視角) 54
圖[4-21]上為壓力分部,下為渦流結構 55
圖[4-22]上為壓力分部,下為渦流結構 55
圖[4-23] 上為壓力分部,下為渦流結構 56
圖[4-24] 上為壓力分部,下為渦流結構 56
圖[4-25]上為壓力分部,下為渦流結構 57
圖[4-26] 上為壓力分部,下為渦流結構 57
圖[4-27]上為壓力分部,下為渦流結構 58
圖[4-28]上為壓力分部,下為渦流結構 58
圖[4-29] 渦流表面流線與圓盤表面壓力圖 59
圖[4-30]翻滾於X-Y平面上軌跡圖 59
圖[4-31] 每25時間步所記錄之圓盤圖 60
圖[4-32]翻滾之運動場域,l1,f1,r1,l2,f2,r2為監測點 60
圖[4-33] ψ隨時間變化圖 61
圖[4-34]速度與受力隨時間變化圖 61
圖[4-35]X-Y方向受力頻譜圖 62
圖[4-36]X-Y方向受力與Z方向受力矩圖 62
圖[4-37] f1點壓力與圓盤受力隨時間變化圖 63
圖[4-38] l1點壓力與圓盤受力隨時間變化圖 63
圖[4-39] r1點壓力與圓盤受力隨時間變化圖 64
圖[4-40] f2點壓力與圓盤受力隨時間變化圖 64
圖[4-41] l2點壓力與圓盤受力隨時間變化圖 65
圖[4-42] r2點壓力與圓盤受力隨時間變化圖 65
圖[4-43] 監測點壓力f1、r1、l1與圓盤Y方向受力頻譜圖 66
圖[4-44] 監測點壓力f2、r2、l2與圓盤Y方向受力頻譜圖 66
圖[4-45]翻滾的路徑上所脫離之渦流與壓力分部圖 67
圖[4-46] 上為壓力分部,下為渦流結構 67
圖[4-47] 上為壓力分部,下為渦流結構 68
圖[4-48] 上為壓力分部,下為渦流結構 68
圖[4-49] 上為壓力分部,下為渦流結構 69
圖[4-50] 上為壓力分部,下為渦流結構 69
圖[4-51] 上為壓力分部,下為渦流結構 70
圖[4-52]圓盤渦流表面流線與圓盤表面壓力分佈圖1 70
圖[4-53] 渦流表面流線與圓盤表面壓力分佈圖2 71
圖[4-54]螺旋運動場域,其中l2、r2、l3、r3為壓力監測點 71
圖[4-55]Z-X平面軌跡圖 72
圖[4-55] X:Y:Z為2:1:2的每25時間步之三維視圖 72
圖[4-56]X與Z方向受力隨時間變化圖 73
圖[4-57]各轉角與隨時間變化圖 73
圖[4-58]角速度隨時間變化圖 74
圖[4-59]X、Z方向之轉速與受力頻譜比較圖 74
圖[4-60] Y方向之轉速與受力頻譜比較圖 75
圖[4-61]各方向受力與X、Z方向速度隨時間變化圖 75
圖[4-62] 各方向受力與Y方向速度隨時間變化圖 76
圖[4-63] l2點壓力與圓盤受力隨時間變化圖 76
圖[4-64] r2點壓力與圓盤受力隨時間變化圖 77
圖[4-65] l3點壓力與圓盤受力隨時間變化圖 77
圖[4-66] r3點壓力與受力隨時間變化圖 78
圖[4-67] 監測點壓力l2、r2與圓盤Y方向受力頻譜圖 78
圖[4-68] 監測點壓力l3、r3與圓盤Y方向受力頻譜圖 79
圖[4-69] 上為壓力分部,下為渦流結構 79
圖[4-70] 上為壓力分部,下為渦流結構 80
圖[4-71] 上為壓力分部,下為渦流結構 80
圖[4-72] 上為壓力分部,下為渦流結構 81
圖[4-73] 渦流表面流線與圓盤表面壓力圖 81
圖[4-74]渦流表面流線與圓盤表面壓力分佈圖 82
圖[4-75] 路徑上脫離之渦流 82
圖[4-76] 路徑上之壓力分佈情形 83
參考文獻 [1] D.Marshell, T.E.Stanton, and F.R.S., On the Eddy System in the Wake of Flat Circular Plates in Three Dimension Flow. The Royal Society, 1930.
[2] F.W.Roos, and W.W.Willmarth, Some Experimental Results on Sphere and Disk Drag, AIAA Journal, Vol 9, No.2 1971.
[3] A.R.Ashenoy,amd C.Kleinstreuer, Flow Over a Thin Circular Disk from Low to Moderate Reynolds number. Journal of Fluid Mechanics, Vol.605, pp. 253-262, 2008.
[4] J.R.Calvert, Experiment on the Flow Past an Inclined Disk, Journal of Fluid Mechanics, Vol.29, part 4, pp.691-703, 1967.
[5]S.Gao, L.B.Tao, X.L.Tian, and J.M.Yang, Journal of Fluid Mechanics. Vol.851, pp. 687-714, 2018.
[6] X.L.Tian, Z.H.Hu, and J.M.Yang. Direct Numerical Simulation of Flow Past an Inclined Disk Journal of Fluid and Structures 72, pp.152-168, 2017.
[7] W.W.Willmarth, N.E.hawk, and R.L.Harvey. Steady and Unsteady Motions and Wakes of Free Falling Disk, The Physics of Fluids 7, 1964.
[8] P.C.Fernandes, F.risso, P.Ern, and J.Magnaudet. Oscillatory Motion and Wake Instability of Freely Rising Axisymmetric Bodies, Journal of Fluid Mechanics, Vol.573, pp.479-502, 2007.
[9] H.J.Zhong, S.Y.Chen, and C.B.Lee, Experimental of Free Falling disks: Transition from Planar Zigzag to Spiral. The Physics of Fluids 23, 011702, 2011.
[10] F.Auguste, J.Magnaudet, and D.Fabre, Falling Style of Disks, Journal of Fluid Mechanics, Vol.719, pp.388-405, 2013.
[11] P.C.Fernandes, P.Ern, Risso, J.Magnaudet. On the zigzag Dynamics of Freely Moving Axisymmetric Bodies. Physics of Fluids 17, 098107, 2005.
[12] S.Scoot M.Murman, M.J.Aftosmis, and M.J.Berger. Implicit Approaches for Moving Boundaries, 41st Aerospace Sciences Meeting and Exhibit, 2003.
[13] S.Y.Lin, Y.H.Chin, J.J.Hu, and Y.C.Chen, A Pressure Correction Method for Fluid-Particle Interaction Flow:Direct-Forcing Method and Sedimentation Flow, International Journal for Numerical Methods in Fluid, Vol.67, pp.1771-1798, 20113
[14] A.T.Cate, C.H.Nieuwstad, J.J.Derksen, and H.E.A.Van Den Akker, Particle Imaging Velocimetry Experiments and Lattice-Boltzmann Simulations on a Sigle Sphere SettlingUnder Gravity, Physics of Fluid 14, 4012, 2002.
[15] H.J.Zhong, C.B. Lee, Z.Su, S.Y.Chen, M.G.Zhou, and J.Z.Wu. Experimental Investigation of Freely Falling Thin Disks, Part 1. The Flow Structures and Reynolds Number Effects on the Zigzag Motion. Journal of Fluid Mechanics, Vol.716, pp.228-250, 2013.
[16] M.Chrust, G.Bouchet, J.Dusek. Numerical simulation of the dynamcis of freely falling discs, Physics of Fluids 25, 044102, 2013.
[17] S.B.Field, M.Klaus, M.G.Moose, and F.Nori. Chaotic dynamics of falling disks, Nature 388, 252-254, 1997.
[18] C.C.Chueh, P.K.Wang, and T.P.Hashino. Numerical Study of Motion of Falling Conical Graupel. Atmospheric Reseach 199(2018)82-92.
[19] K.O.L.F.Jayweera and B.J.Mason. The Behaviour of Freely Falling Cylinders and Cones in a Viscous Fluid, Journal of Fluid Mechanics, Vol.22, part 4, pp.709-720, 1965.
[20] A.Andersen, U.Pesavento, and Z.J.Wang. Unsteady Aerodynamics of Fluttering and Tumbling Plates, Journal of Fluid Mechanics, Vol. 541, pp.65-90, 2005.
[21] P.Ern, F.Risso, D.Fabre, and J.Magnaudet. Wake-Induced Oscillatory Paths of Bodies Freely Rising or Falling in Fluids, Annual Reviews Fluid Mechanics, 44:97-121,2012.
[23]P.Meliga, J.M.Chomaz, and D.Sipp. Global Mode Interaction and Pattern Selection In the Wake of a Disk : a Weakly Nonlinear Expansion, Journal of Fluid Mechanics, Vol.633, pp.159-189. 2009.
[24]P.C.Fernandes, F.Risso, P.Ern and Jacques Magnaudet, Osillatory Motion and Wake Instability of Freely Rising Axisymmetric Bodies, Journal of Fluid Mechanics, Vol.573, pp.479-502, 2007.
[25] R.J.Wu, and S.Y.Lin, The Flow Of a Falling Ellipse: Numerical Method and Classification, Jounary of Fluid Mechanics, Vol.31, No.6, 2015.
[26] L.Vicent, W.S.Shambaugh and E.Kanso, Holes Stablize Freely Falling Coins, Journal of Fluid Mechanics, Vol.801, pp.250-259.2016.
[27] Haervig, Jakob, Jensen, A.Lyhne; Pedersen, Marie Cecilie; and Sorensen Henrik, General Observations of the Time-Dependent Flow Field Around Plates in Free Fall, Fluids Engineering Conference [12708] (pp.978-983), 2015.
[28] A.BelMonte, H.Eisenberg, and E.Moses, From Flutter to Tumble: Inertial Drag and Froude Similarity in Falling Paper, Physical Review Letters, Vols.81, Number 2. 13 July 1998.
[29]戴宏源,重疊式網格流場數值方法之發展,成功大學航太所九十一學年度碩士論文,2002。
論文全文使用權限
  • 同意授權校內瀏覽/列印電子全文服務,於2019-08-02起公開。
  • 同意授權校外瀏覽/列印電子全文服務,於2019-08-02起公開。


  • 如您有疑問,請聯絡圖書館
    聯絡電話:(06)2757575#65773
    聯絡E-mail:etds@email.ncku.edu.tw