黃良雄臺灣大學:土木工程學研究所溫嘉琪Wen, Chia-ChiChia-ChiWen2010-06-302018-07-092010-06-302018-07-092008U0001-1411200815055800http://ntur.lib.ntu.edu.tw//handle/246246/187642中文摘要鍵詞：水平潛沒平板，勢流-渦流分解法，渦旋法，邊界積分法，反射率，高階諧和波研究旨在應用Lin and Huang(2008)利用勢流-渦流分解法所發展之數值方法，進行週期水波通過一水平潛沒平板之數值模擬研究。此方法的特色在於同時應用了邊界積分法(boundary integral method)與渦旋法(vortex method)的優點。此數值方法應用Helmholtz分解法，將流場分為具有渦度之旋性流場與不具渦度之非旋性流場；其中，利用渦旋法(vortex method)來求解渦度場，另外，以邊界積分法求解非旋性流場與自由水面之運動。此方法保留勢流架構下的優點，能以邊界積分法處理自由水面；另一方面，此方法不同於傳統以使用網格的數值方法，毋須對整個計算域進行黏性流求解，而可只針對流場中具有黏滯作用之範圍集中計算。了瞭解水波通過潛沒平板的表現行為，本研究進行一系列水波通過一潛沒薄平板之數值模擬研究。對一薄平板而言，其產生的渦流僅侷限於二端點附近，對於整個流場，此渦流現象不甚明顯；而相較於薄平板本身，二端點之渦度表現範圍則較平板厚度大了許多，顯示出渦流現象並不可被忽略。為能充分比較勢流與黏流分析之差異，本文先將渦流忽略，採用勢流的理論進行分析，並探討反射率與高階諧和波的生成。而後，本文將探討在黏滯效應的作用之下，水波通過平板之渦流生成現象。此外，本文並將針對線性水波之旋性流場與非旋性流場，探討線性水波條件下，以勢流理論分析之適用性。由本研究之數值模擬成果發現，分別於勢流與渦流之假設條件之下，潛沒平板之置放深度對於水波變形之影響極有重要性；且對於薄平板而言，渦流效應不可忽略，其對於反射率值與主頻波振幅的變化皆有一定之影響。而對於水波通過薄平板之問題，即使於線性水波的條件下，薄平板之兩端點仍會產生明顯的小尺度渦流，雖然此渦流散逸現象並不明顯，但其對於平板端點處之有效邊界的影響，將使得勢流理論的分析容易產生較大的誤差。Abstracteywords：Submerged flat plate, Potential-Vorticity decomposition, Boundary integral method, Vortex method, Reflection coefficients, Higher harmonics. n this research, a numerical scheme that employs a potential-vorticity decomposition is used to investigate the interaction of periodic water waves with a submerged flat plate. This method uses Helmholtz decomposition to decompose the flow field into its rotational and irrotational parts. The vortical flow field is solved via a vortex method; the irrotational flow field and the motion of the free surface are solved using a boundary integral technique. The major advantage of this method is the efficiency of the boundary integral method for solving the free surface motion, and the essentially grid-free nature of the vortex method for the vorticity field, which is predominantly confined in compact regions. series of simulations were conducted to study the free surface deformation and the flow pattern. Although the vortical regions are mainly confined near the two sharp edges of the plate, their scales are relatively larger then the thickness of the thin plate. In order to compare the results obtained from a potential-flow approach and the viscous-flow model, at first we present the numerical results in which the vortical part is neglected, and then the generation and evolution of vortices as well as their effects are reported.ccording to the numerical results, the position of the thin flat plate affects the deformation of the free surface significantly. Furthermore, vortex effects are not negligible for the flow near the flat plate, and may affect the reflection of surface waves and the generation of higher harmonics. Even if linear waves cases, vortices around the tips of the thin plate affected the effective boundary. This may lead significant error in the analysis using a potential flow theory.目錄頁碼文摘要 Ⅰbstract Ⅲ錄 Ⅴ表目錄 Ⅷ號表 Ⅹ一章 緒論 1.1 研究背景與動機 1.2 文獻回顧 2.3 研究方法 3.4 論文架構 5二章 數學理論之建立 6 2.1 控制方程式 6 2.2 邊界條件 6 2.2.1 固體表面邊界條件 6 2.2.2 自由水面之運動邊界條件(KFSBC) 7 2.2.3 自由水面之動力邊界條件(DFSBC) 7 2.2.4 側向邊界條件 8 2.3 Helmholtz分解之積分表示式 8.4 渦度邊界條件之黏滯性分離演算(viscous splitting algorithm) 11.5 與 12 2.5.1 的演算 12 2.5.2 的演算 14三章 數值方法 16 3.1 二維流體之渦旋法(vortex method) 16 3.2 數值計算流程 17 3.3 波的衍生與吸收 20 3.4 次網格模擬(sub-grid modeling) 22 3.5 映射法(method of image)於平坦底床邊界條件之應用 22 3.6 無因次化與穩定性條件 23 3.7 薄平板邊界條件之處理 25四章 結果與分析 27 4.1 非旋性(irrotational)流場分析 27 4.1.1 收斂性測試(convergence study) 27 4.1.2 反射率 30 4.1.3 高階諧和波之生成 36 4.2 旋性(rotational)流場分析 54 4.2.1 渦流的生成與消散 54 4.2.2 渦流與勢流之比較 83 4.2.2.1 反射率 83 4.2.2.2 波的變形 84 4.3 線性水波之旋性流場與非旋性流場 86 4.3.1 渦流的生成與消散 86 4.3.2 波的變形與反射行為 95五章 結論與建議 100 5.1 結論 100 5.2 建議 101考文獻 102表目錄頁碼2.1 數值條件示意圖 153.1 映射系統示意圖 243.2 數值模擬步驟圖 253.3求解固體邊界之切線速度示意圖 264.1不同參數設定下之水面波形比較圖 294.2 平板置放深度不同之反射率比較圖 324.3水波通過平板之自由水面波形變化 334.4 m之不同板長的諧和波振幅在空間上的變化 384.5 m之不同板長的諧和波振幅在空間上的變化 434.6 m與 m之不同板長的諧和波振幅在空間上的變化 494.7 渦度場與流場之輸出範圍示意圖 554.8 m, =0.84m之渦度場變化圖 564.9 m, =0.84m之流場與自由水面變化圖 614.10 m, =0.84m之渦度場變化圖 664.11 m, =0.84m之流場與自由水面變化圖 714.12 m, =0.84m, T之平板左端周圍質點速度示意圖 764.13 m, =0.84m, T之平板右端周圍質點速度示意圖 764.14 m, =0.84m之流場與自由水面變化圖 784.15 渦流與勢流情況下之高階諧和波比較 844.16 線性水波之渦度場(含平板左右端點周圍)變化圖 874.17 平板周圍(含平板左右端點周圍)質點速度示意圖 924.18 渦流與勢流情況下之諧和波比較圖(線性水波) 964.19 勢流與渦流情況下之自由水面波形變化(虛線：渦流；實線：勢流) 974.1 非旋性數值收斂性測試之條件參數設定 284.2 各案例之數值條件 354.3 渦流與勢流之反射率比較 834.4 線性水波案例之數值條件 863043075 bytesapplication/pdfen-US勢流-渦流分解法渦旋法邊界積分法反射率高階諧和波Submerged flat platePotential-Vorticity decompositionBoundary integral methodVortex methodReflection coefficientsHigher harmonicsthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/187642/1/ntu-97-R95521314-1.pdf