专题二期末报告(蔡育仁) .docVIP

國 立 台 灣 海 洋 大 學 專題二期末報告 題目:半平面含阻抗邊界條件Laplace 問題之格林函數 指導教授:陳正宗 系級:河工4A 姓名:蔡育仁 學號:Green’s function of a half-plane problem subject to impedance boundary Jeng-Tzong Chen1, 2, Yu-Ren Tsai1 NSC-100-2815-C-019 -011-E 1Department of Harbor and River Engineering, National Taiwan Ocean University 2Department of Mechanical and Mechatronic Engineering, National Taiwan Ocean University Keywords: impedance boundary condition, method of fundamental solutions, half plane problem, Green’s function Abstract In this paper, the Green’s function of a half-plane problem subject to impedance boundary condition is discussed. Greenberg has derived an analytic solution by using the perturbation method from Neumann to impedance case in his book. We intend to follow his method and use the perturbation method from Dirichlet to impedance case. Here, we employ the method of fundamental solutions (MFS) to solve half-plane problems subject to Dirichlet, Neumann and impedance boundary conditions. The relationship between the strength distribution of the sources in MFS and the analytical solution derived by Greenberg will be constructed. The reason why the MFS results in the ill-posedness for solving this problem is also addressed in this paper. Finally the analytical solutions are compared with MFS demonstrate the validity of the present method. 摘要 本文探討半平面含阻抗邊界條件Laplace問題之格林函數，Greenberg曾於其書中，使用微擾法由自由端邊界條件推導至阻抗邊界條件。本文將使用微擾法由固定端邊界條件推導至阻抗邊界條件。基本解法求解格林函數在Neumann、DirichletGreenberg所導之解析解的關係做更進一步探討。並針對基本解法在處理半平面含阻抗邊界條件問題時所導致影響係數矩陣產生病態的原因進行了解。本文所得之解析解與基本解法之結果進行比較，以驗證本法的可行性。問題描述 含阻抗邊界條件之半平面問題可將其表示如下圖一 圖一 含阻抗邊界條件之半平面問題 其控制方程 (1) 其中為Dirac-Delta函數，阻抗邊界條件為 (2) Greenberg曾在其書中[1]使用微擾法推導半平面含阻抗邊界條件之格林函數，其解析公式表示如下: (3) 由(3)式得知格林函數包含特解項及補解項，特解項因包含源點，故會使其產生奇異現象，其等高線圖如下圖一所示，而補解項則不包含奇異點，故不會產生奇異現象，其等高線圖如下圖二所示:  圖一 格林函數等高線圖 圖二 補解等高線圖 公式推導 由固定端邊界條件(Dirichlet)推導至阻抗邊界條件 依照Greenberg[9]書中之推導方法，我們可以嘗試由固定端(Dirichlet)邊界條件使用微擾法延伸推導至阻抗係數為的阻抗邊界條件，如圖三所示 圖三 格林函