边界元素法期末报告.ppt

邊界元素法期末報告 Solution for wave propagation through a circular cylinder mounted on different topography ripple-bed profile shoals using DRBEM 在淺水處以雙互換邊界元素法求解波浪傳播透過被設置在不同地形之波動式底床的圓筒 授課老師:陳正宗 老師 報告人:李佩諭 2013.06.13 目錄 前言 研究理論 1.修正型緩坡方程式 2.數值模式 模式驗證 結論 前言(1/2) 緩坡方程式(Mild-Slope Equation, MSE)由Berkhoff提出以線性波理論，並引入依反映地形變化之微小參數，以攝動法展開得到一二階之微分方程式，其求解方式以有限差分法或是有限元素法，但這兩種方法無法處理複雜地形及無窮邊界條件問題。 因此，本文採用修正型緩坡方程式(Modified Mild-Slope Equation, MMSE)做為波浪控制方程式。將方程式化為Helmholtz方程式形式，藉由基本解Hankel函數自動滿足Sommerfeld輻射邊界條件的優點。 前言(2/2) 並以雙互換邊界元素法(Dual Reciprocity Boundary Element, DRBEM)求解，其能仍保有原本傳統的邊界元素法，並對於非齊次項(Non-homogenous term)藉由徑向基函數再次離散至邊界點的處理，則可以避免複雜的領域積分。 二.研究理論(1/3) 修正型緩坡方程式如下: 為波速， 為群波速， 為水平梯度運算子， 為底床斜率項， 為底床曲率項， 二.研究理論(2/3) 數值模式:計算領域可分為兩部分，外部等水  深領域以及內部變水深領域。 由於控制方程式MMSE在外  部等水深區域時會退化成為 Helmholtz方程式 而基本解為 利用Green第二恆等式得 二.研究理論(3/3) 由於海洋是無窮域，而數值計算則為有限域 ，故有限域上的邊界條件卻必須滿足無窮域的物理現象。因此在計算領域邊界 必須滿足 Sommerfeld 提出之輻射邊界條件。 最後可藉以下方程式得變水深區域之速度勢  能 模式驗證 首先對數值模式之準確定性做測試(利用圓錐形島嶼) 結論 由圖可知f1以及f 2變化最劇烈之範圍落在中間性水深條件。因此，可以推論MSE 與MMSE 之間的差異量皆由f1以及f 2所造成。 基本的淺水地形 波動底床的淺水地形 基本的淺水地形 波動底床的淺水地形 波動底床的淺水地形 基本的淺水地形 若地形呈平滑狀態時，其受底床曲率及底床坡度平方項影響較小，而利用MSE可得到其數值結果。 若為一正弦曲線的地形，底床曲率項變化會較明顯，底床坡度平方項變化較小。 若延伸項(底床曲率及底床坡度平方項)不被描述，當地形是有起伏狀態時，則其結果是不準確的。 謝謝大家聆聽 * *