跳到主要內容

臺灣博碩士論文加值系統

(3.231.230.177) 您好!臺灣時間:2021/08/04 01:19
字體大小: 字級放大   字級縮小   預設字形  
回查詢結果 :::

詳目顯示

我願授權國圖
: 
twitterline
研究生:陳啟宏
研究生(外文):Chi-Hom Chen
論文名稱:以最小平方有限元素法探討雙橋墩流場變化
論文名稱(外文):Study of the flow fields in a dual bridge-pier by least-squares finite element method
指導教授:鄧志浩鄧志浩引用關係
指導教授(外文):Jyh-Haw Tang
學位類別:碩士
校院名稱:中原大學
系所名稱:土木工程研究所
學門:工程學門
學類:土木工程學類
論文種類:學術論文
論文出版年:2005
畢業學年度:93
語文別:中文
論文頁數:100
中文關鍵詞:雙橋墩最小平方有限元素法淺水波方程式
外文關鍵詞:dual piers of a bridgeleast-squares finite element methodshallow water equations
相關次數:
  • 被引用被引用:4
  • 點閱點閱:119
  • 評分評分:
  • 下載下載:0
  • 收藏至我的研究室書目清單書目收藏:0
以數值模擬方法研究流體經過雙橋墩流場變化。此方法使用最小平方有限元素法求解淺水波方程式。從模擬結果得知,接近橋墩時水深會壅起而過橋墩後方水深會下陷。當雙橋墩間距增加時,其干涉現象會由大逐漸變小。當兩橋墩間距超過4倍的直徑時,干涉現象可以忽視的。對於縱向排列的雙橋墩,干涉現象會較小,而其水深的分布也會類似單一橋墩現象。兩橋墩之間的水深隨著兩橋墩間距的增加水深會降低。本文研究,水深分布可有效地模擬和可能壅起的水體預測。希望結果能為水利工程學界被廣泛設計應用。
A numerical simulation of the flow pass dual piers of a bridge has been presented in this study. The least-squares finite element method is employed for solving the governing shallow water equations. From the simulation results, it can be concluded that the water depth is accumulated proximal to the pier and drops immediately after the pier. The effect of the interference caused by the dual piers is strong and become weaker as the distance between the piers is increased. When the space between two piers is more than four times of the diameter of the pier, the effect of the interference is negligible. For the case of dual piers formed in a line along the flow direction, the effect of interference is small and the formation of the water depth is very similar to the distribution of a single pier. The water depth between two piers is lowered as the distance between piers increasing. In this study, the water depth distribution can be effectively simulated and the possible accumulation of water body is calculated. It is hoped that the results can
be used for the design purpose for the hydraulic engineering.
本 文 目 錄
第一章� 緒論........................................................................01
1-1前言.........................................................................01
1-2研究目的與方法...............................................................02
1-3文獻回顧......................................................................04
第二章� 理論分析....................................................................10
2-1 LSFEM簡介.................................................................10
2-2理論推導...................................................................11
2-3淺水波方程式簡介...............................................................14
2-3-1最小平方有限元素法應用於淺水波方程式............................16
2-4 共軛梯度法介紹(Conjugate Gradient Method) .....................22
2-4-1元素-元素技巧(Element-by-Element Technique)...................23
2-4-2共軛梯度法(conjugate Gradient Method)簡介.......................24
2-5庫倫數(CFL)條件.........................................................................26
第三章� 模式建立及邊界條件設定.........................................28
3-1單一橋墩流場驗證.............................................................28
模式建立.........................................................................28
邊界條件設定.........................................................................30
3-2具自由液面流體受結構物影響之數值驗證..............................30
模式建立.........................................................................30
初始條件設定.........................................................................31
邊界條件設定.........................................................................31
第四章� 結果與分析.........................................................................36
4-1單一橋墩模擬結果與分析.......................................................36
4-2雙橋墩模擬結果與分析.......................................................37
4-3縱向雙橋墩模擬結果與分析...........................................................39
4-4 單一橋墩與雙橋墩間模擬結果分析及比較..............................41
第五章� 結論與建議.........................................................................43
5-1結論.........................................................................43
5-2建議.........................................................................44
參考文獻.........................................................................46

圖 目 錄
圖1-1流體流經圓柱所產生的分離現象及尾流渦( Morkovin 1964).. 05
圖1-2各種雷諾數值下之圓柱尾流型態( Morkovin 1964).............06
圖3-1-1網格示意圖(Yulistiyanto 等, 1998) ........................................29
圖3-1-2 網格示意圖...............................................................................29
圖3-1-3邊界示意圖...............................................................................30
圖3-2-1單一橋墩滑動邊界示意圖 .....................................................32
圖3-2-2單一橋墩非滑動邊界示意圖圓柱上之速度設為零 .............32
圖3-2-3並排雙橋墩間距示意圖 .........................................................33
圖3-2-4縱向雙橋墩間距示意圖 .........................................................34
圖3-2-5單一橋墩網格示意圖 .............................................................34
圖3-2-6並排雙橋墩網格示意圖(以2D間距) ....................................35
圖3-2-7縱向雙橋墩網格示意圖(以2D間距) .......................................35
圖4-1-1速度向量圖(Yulistiyanto 等, 1998) .......................................50
圖4-1-2速度向量圖................................................................................50
圖4-1-3水深變化圖(Yulistiyanto 等, 1998) .........................................51
圖4-1-4水深變化圖................................................................................51
圖4-1-5 滑動邊界與非滑動邊界(單一橋墩)水深變化 T=4
時間4秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖.................52
圖4-1-6 滑動邊界與非滑動邊界(單一橋墩)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 .................53
圖4-1-7 滑動邊界與非滑動邊界(單一橋墩)流線變化 T=12 .............54
圖4-1-8 明顯的看出滑動邊界圓柱下游後方不會有渦度產生 .........54
圖4-1-9非滑動邊界單一橋墩T=42~56秒間的橋墩後方波街變化 .56
圖4-1-10 單一橋墩選一間段時間及P點位置(6,0)與垂直速度的關係
圖 Time step 850~1150(t=42.5~57.5) ...................................57
圖4-1-11 雙橋墩間距1.25倍選一間段時間及P點位置(6,0)與垂直速
度的關係圖Time step 850~1150(t=42.5~57.5) ......................57
圖4-1-12 雙橋墩間距1.5倍選一間段時間及P點位置(6,0)與垂直速度的關係圖Time step 850~1150(t=42.5~57.5) .........................58
圖4-2-1 滑動邊界與非滑動邊界(雙圓柱橋墩1.25倍)水深變化 T=4
時間4秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............59
圖4-2-2 滑動邊界與非滑動邊界(雙圓柱橋墩1.25倍)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............60
圖4-2-3 滑動邊界與非滑動邊界(雙橋墩間距1.5倍)水深變化 T=4
時間4秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............61
圖4-2-4滑動邊界與非滑動邊界(雙橋墩間距1.5倍)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............62
圖4-2-5 滑動邊界與非滑動邊界(雙橋墩間距2倍)水深變化 T=4
時間4秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............63
圖4-2-6滑動邊界與非滑動邊(雙橋墩間距2倍)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............64
圖4-2-7 滑動邊界與非滑動邊界(雙橋墩間距3倍)水深變化 T=4
時間4秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............65
圖4-2-8 滑動邊界與非滑動邊界(雙橋墩間距3倍)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...............66
圖4-2-9 滑動邊界與非滑動邊界(雙橋墩間距4倍)水深變化 T=4
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 .................67
圖4-2-10 滑動邊界與非滑動邊界(雙橋墩間距4倍)水深變化 T=6
時間6秒水深變化滑動邊界(a)2D示意圖(c) 3D示意圖
非滑動邊界(b)2D示意圖(d) 3D示意圖 ...................68
圖4-2-11 非滑動邊界(雙橋墩間距1.25倍)流線變化 ........................70
圖4-2-12 非滑動邊界(雙橋墩間距1.5倍)流線變化 ..........................71
圖4-2-13 非滑動邊界(雙橋墩間距2倍)流線變化 .............................72
圖4-2-14非滑動邊界(雙橋墩間距3倍)流線變化 ..............................73
圖4-2-15非滑動邊界(雙橋墩間距4倍)流線變化 ..............................75
圖4-3-1 縱列雙橋墩1.25倍水深變化4秒(a)2D示意圖
(b)3D示意圖及A、B、C點水深高 ...................................................75
圖4-3-2 縱列雙橋墩1.25倍水深變化5秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E點水深高 ......................................75
圖4-3-3 縱列雙橋墩1.25倍水深變化6秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F點水深高 ...............................76圖4-3-4 縱列雙橋墩1. 5倍水深變化4秒(a)2D示意圖
(b)3D示意圖及A、B、C點水深高 ...................................................76
圖4-3-5 縱列雙橋墩1. 5倍水深變化5秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E點水深高 ....................................77
圖4-3-6 縱列雙橋墩1. 5倍水深變化6秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F點水深高 .................................77
圖4-3-7 縱列雙橋墩2倍水深變化4秒(a)2D示意圖
(b)3D示意圖及A、B、C點水深高 .....................................................78
圖4-3-8 縱列雙圓柱2倍水深變化5秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、G點水深高 ................................78
圖4-3-9 縱列雙橋墩2倍水深變化6秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F、G點水深高 .........................79
圖4-3-10 縱列雙橋墩2倍水深變化7秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F點水深高 .................................79
圖4-3-11 縱列雙橋墩3倍水位變化4秒(a)2D示意圖
(b)3D示意圖及A、B、C點水深高 ...................................................80
圖4-3-12 縱列雙橋墩3倍水位變化5秒(a)2D示意圖
(b)3D示意圖及A、B、C、D點水深高 .............................................80
圖4-3-13 縱列雙橋墩3倍水位變化6秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、G點水深高 ..............................81
圖4-3-14 縱列雙橋墩3倍水位變化7秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F、G點水深高 ........................81
圖4-3-15 縱列雙橋墩4倍水位變化4秒(a)2D示意圖
(b)3D示意圖及A、B、C點水深高....................................................82
圖4-3-16 縱列雙橋墩4倍水位變化5秒(a)2D示意圖
(b)3D示意圖及A、B、C、D點水深高 ............................................82
圖4-3-17 縱列雙橋墩4倍水位變化6秒(a)2D示意圖
(b)3D示意圖及A、B、C、D點水深高 ..............................................83
圖4-3-18 縱列雙橋墩4倍水位變化7秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、G點水深高 ...............................83
圖4-3-19 縱列雙橋墩4倍水位變化8秒(a)2D示意圖
(b)3D示意圖及A、B、C、D、E、F、G點水深高 ............................84
圖4-3-20 縱列雙橋墩1.25倍流線變化15秒 .....................................84
圖4-3-21 縱列雙橋墩1.25倍流線變化25,35,45,60,70,80秒 ..86
圖4-3-22縱列雙橋墩1.5倍流線變化15,25,35,45,60,70,80
秒 ..........................................................................................................87
圖4-3-23 縱列雙橋墩2倍流線變化15,25,35,50,70,80秒 ...89
圖4-3-24 縱列雙橋墩3倍流線變化15,25,40,55,70,80秒 ...90
圖4-3-25 縱列雙橋墩4倍流線變化15,25,40,55,70,80秒 ...92
圖4-4-1橫向雙橋墩各不同間距間時間與水深之比較 ......................92
圖4-4-2單一橋墩在T=40秒之水深變化圖.........................................93
圖4-4-3單一橋墩在T=50秒之水深變化圖........................................93
圖4-4-4單一橋墩在T=40秒之水深變化3D示意圖.........................94
圖4-4-5單一橋墩在T=50秒之水深變化3D示意圖.........................94
圖4-4-6橫向橋墩間距2倍T=30秒以中心線為主(y=0)之水深變化圖
..................................................................................................................94
圖4-4-7橫向橋墩間距2倍T=30秒以橋墩後方(y= )之水深變化圖
..................................................................................................................95
圖4-4-8橫向橋墩間距4倍T=30秒以中心線(y=0)之水深變化圖......95
圖4-4-9橫向橋墩間距4倍T=30秒以橋墩後方(y= )之水深變化圖
..................................................................................................................96
圖4-4-10縱向橋墩間距1.25倍T=30秒以中心線為主(y=0)之水深變化圖..........................................................................................................96
圖4-4-11縱向橋墩間距2倍T=30秒以中心線為主(y=0)之水深變化圖
..................................................................................................................97
圖4-4-12縱向橋墩間距4倍T=30秒以中心線為主(y=0)之水深變化圖
..................................................................................................................97

表 目 錄
表4-1 流場在時間T=4所發生的水深變化表....................................98
表4-2 流場在時間T=6所發生的水深變化表....................................99
表4-3 流場在時間T=4所發生的水深變化表....................................99
表4-4 流場在時間T=5所發生的水深變化表...................................100
表4-5 流場在時間T=6所發生的水深變化表...................................100
參考文獻
1. Bradford, Scott F. and Sanders, Brett F. “Finite-Volume Model for Shallow-Water Flooding of Arbitrary Topography,” Journal of
Hydraulic Engineering, March.,2002,pp 289-298.
2. Bramble, J.H. and Shatz, A.H., 1970. ”On The Numerical Solution
Of Elliptic Boundary-Value Problems By Least-Squares Approximation Of The Data.” in B.Hubbad(ed). Numerical Solution of
P.D.E., vol.2,Academic press. New York., pp.107-133.
3. Fennema, R J. and Chaudhry, M.H.“Implicit methods for two-
dimensional unsteady free-surface flow.” J ydr. Res., 7(3),1989,pp.321-
331.
4. Gau, C., Wu, J. M. and Liang, C. Y. "Heat Transfer Enhancement and Vortex Flow Structure over a Heated Cylinder Oscillating in the Cross
flow Dirction", Joural of Heat Transfer, vol.121, 1999, pp.789-795.
5. Guillou, S. and Nguyen, K.D.,” An Improved Technique for Solving Two-Dimensional Shallow Water Problems,” International Journal for
Numerical methods in Fluids, Vol 29,1999: pp 465–483.
6. Galland, J.C., Goutal, N. and Hervouet, J.M, ‘Telemac: A new numerical model for solving shallow water equations’, Ad6. Water
Resources, 14,1991,pp 138–148 .
7. Igarashi, T., "Characteristics of the Flow Around Two Circular Cylinders Arranged in Tandem", Bulletin of the JSME, vol.24, no.188,
1981, pp.532-582.
8. Jimenez, Oscar F. and Chaudhry, M. Hanif, Member , ASCE.”Computa
. ” Journal of Hydraulic Engineering., Vol.114,No4, April , 1988 , pp
377-395.
9. Jiang, Bo-Nan, and Carey, G. F. ,“ The Least-Squares Finite Element Method for Non-Linear Hyperbolic Problem,”Int. J. Num Meth.Engrg.
, Fluid, 8,1988, pp933-942.
10. Jiang, Bo-Nan and Povinelli Louis a “last-squares finite element method for fluid dynamics.” Comp. Meth. Appl. Mech. Engrg.,82,1990
, pp365-383.
11. Jiang, Bo-Nan .“A Least-Square Finite Element Method for Incompressible Navier-Stokes Problem.”Int. J. Num Method Fluids.
, Vol.14,1992,pp843-859.
12. Jiang, Bo-Nan. “The Least-Squares Finite Element Method: A Unified
Method for Numerical Solution of PDE.”,1998.
13. Jiang, Bo-Nan. “The Least-Squares Finite Element Method:”Theory and Applications in Computational Fluid Dynamics and
Electromagnetics “,1998.
14. Lin, J. C., Yang, Y., and Rockwell, D., "Flow past Two Cylinders
in Tandem:Instantaneous and Averaged Flow Structure", Journal of
Fluids and Structures, 2002, pp.1059-1071.
15. Lynn, P.P. and Alani, K., 1974. “Use of The Least Square Criterion
In Finite Element Formulation.”, Int. J. Num. Method Eng.,
vol.6 ,pp75-88.
16. Molls, Thomas and Chaudhry, M.Hanif . “depth-averaged open-
channel model.” Journal of Hydraulic Engineering.,Vol.121,NO.6,June,
1995,pp 453-465.
17. Ng, C. W. and Ko, N. W. M., "Flow Interaction behind Two Circular
Clinders of Equal Diameter-a Numerical Study", Journal of Wind
Engineering and Industrial, Aerodynamics, 1995, pp.277-287.
18. Patnaik, B. S. V., Narayana., P. A. A., and Seetharamu, K. N., "
Finite Element Simulation of Transient Laminar Flow past a Circular
Cylinder and Cylinders in Tandem", International Journal of Numerical
Methods for Heat and Fluid Flow, vol.10, no.6, 2000,pp.560-580.
19. Slaouti, A. and Stansby, P. K., "Flow Around Two Circular Cylinders
by the Random-Vortex Method", Journal of Fluids and Structures, vol.6, 1992, pp.641-70.
20. Stansby, P. K., "A Numerical Study of Vortex Shedding from One and
Two Circular Cylinders",Aeronautical Quarterly, 1981, pp.48-71 .
21.Tang, L.Q. and Tsang, T.T.H.,. “A Least-Squares Finite
Element Method For Time-Dependent Incompressible Flow With
Thermal Convection.”, Int. J. Num. Methods Fluids., vol.17, 1993
pp.271-289.
22. Terzidis, G. and Strelkoff,T. “Computationof open-channel
surges and shocks .”J.Hydr.Div., ASCE,96(12), 1970,pp2581-2610.
23. Tseng, Ming Hseng and Chu Chia R.“Two-dimensional shallow water
flows simulationusing TVD-Mac Cormack Scheme.” Journal of
Hydraulic Research.,Vol.38,NO.2,Octoner, 2000,pp123-130.
24. Von Karman,T., “Uber den Mechanismuss des Windersstandes den
ein bewegter Korper in einen Flussigkeit Erfahart,” Nachrichten
der K. Gesellschaft der Wissenschaften zu Gottingen, 1912,
pp.547-556.
25. Yang, J. Y. and Hsu, C.A., “ Computation of free surface flows part2:
two-dimensional unsteady bore diffraction, ” Journal of
Hydraulic Research, Vol.31, NO.3, December, 1993,pp 403-414.
26. Yulistiyanto, B., Zech, Y. and Graf, W. H., “Flow Around A
Cylinder:Shallow-WaterModeling with Diffusion-Dispersion,” Journal
of Hydraulic Engineering, April, 1998,pp 419-429.
27. Zdravkovich, M. M., “The Effects of Interference between Circular
Cylinders in Cross Flow,” Journal of Fluids and Structures, Vol.1, 1987
, pp. 239-261.
28. Zhao, D. H., Shen, H. W., Lai ,J. S. and Tabios ,G.Q. Approximate Riemann Solners in FVM for 2D Hydraulic Shock Wave Modeling ,
” Journal of Hydraulic Engineering, ASCE, 122(12),1996,pp 692-702.
29. Zhou, J. G. “A lattice Boltzmann model for shallow water
wquations,” Comput. Methods Appl. Mech. Engrg. Vol.191, 2002,
pp3527-3539.
30. Fluent 計算熱流軟體訓練課程,國家高速電腦中心,2003
31. 吳豐盛(2000),“應用最小平方有限元素法模擬流體經過圓球之
流況,” 私立中原大學土木工程學系碩士學位論文.
32. 賴慶賢(2002),“最小平方有限元素法在二維淺水波方程式之數
值模擬,” 私立中原大學土木工程學系碩士學位論文.
33. 賴盈達(2003),”二維淺水波在不規則底床運行之研究”私立中原
大學土木工程學系碩士學位論文.
QRCODE
 
 
 
 
 
                                                                                                                                                                                                                                                                                                                                                                                                               
第一頁 上一頁 下一頁 最後一頁 top