半平面含阻抗边界条件Laplace问题之格林函数推导-NTOU-海洋大学.doc

100年度國科會大專生參與專題研究計畫(半平面含阻抗邊界條件Laplace問題之格林函數推導) PAGE 21 半平面含阻抗邊界條件Laplace問題之格林函數推導 (一) 摘要 本計畫擬推導半平面含阻抗邊界條件Laplace問題之格林函數。阻抗邊界條件是由、兩個參數所組成的方程式()來表示。當或時可將問題的邊界條件各別簡化成Neumann或Dirichlet型邊界條件。在且的阻抗邊界情況下，針對此問題的映像源分佈，則是我們感興趣的議題。本計畫擬將基本解展開成退化核的形式，希望有別於傳統教科書中以微擾法及相似三角形的方式來推求。此外，我們亦將使用基本解法來比對驗證本法的正確性，並使用Matlab程式語言或Mathematica軟體來模擬二維平面問題及三維空間問題之映像源分佈情形與基本解法中源點強度分佈情況作比對。最後，誤差分析和源點最佳分佈位置也將一併探討。  (二) 研究動機與研究問題: 申請人於2010年修習陳正宗老師於海洋大學河海工程研究所開設之偏微分方程課程，故對偏微分方程在數學及工程等相關方面的應用產生濃厚的興趣。因此，希望能夠在此方面有更進一步了解。申請人於2010年代表海洋大學河海工程學系參加中華顧問工程司所舉辦之第一屆大專生工程力學精英研習營競賽。經過一番辛苦的努力與評比後，有幸獲得工程力學精英研習營表現優異獎學金的肯定，如附件一所示。並獲邀擔任第二屆工程力學競賽之助教。能得此殊榮，更激起申請人對於做力學研究的興趣。在修習專題課程時，進入海洋大學河海工程學系力學聲響振動實驗室( NTOU/MSV)研究團隊。看到學長姐在歷年國科會的贊助下，均有成長且在不錯的學術論文期刊可以發表論文，如附表一所示。此一因素下更激發申請人欲到NTOU/MSV團隊學習做研究的方法，並培養自己做研究的能力。 申請人先前修習偏微分方程課程時，對於半無窮長繩，在固定端與自由端邊界條件下之映像法的推導已有些許了解及心得。希望藉由這次專題計畫的聯想，我們可將一維半無窮長繩的映射問題，延伸至二維的平面問題，並可嘗試使用微擾法等方法推回半無窮長繩在含阻尼器及彈簧時的映射。並且在做偏微分方程課程之作業時，曾使用過Mathematica軟體。此軟體具強大的符號運算與繪圖功能，可以用來模擬半無窮長繩在各種邊界條件情況下的波動動畫；此部分研究成果NTOU/MSV學長們已在CAEE之SCI源期刊發表，如附表一所示。學習過程中，申請人體會到許多教科書美中不足的地方，就是僅能以靜態呈現，造成學生學習上的困難。而陳正宗老師所開的偏微分方程課程，是以動畫的方式來表達數學公式所蘊含之物理現象，將它活生生的呈現於學生面前，著實打開了申請人的眼界。其中，海洋大學NTOU/MSV的教學網站(.tw)是很受歡迎的研究與學習並重的網站，所以申請人認為以此模式來呈現研究成果，必定有極深的教育意義。 本計劃擬先推導半平面含阻抗邊界條件Laplace問題之格林函數，嘗試由固定端(Dirichlet)邊界條件，使用微擾法推導至阻抗邊界，希望有別於Greenberg由自由端的推導。並將基本解法展開成退化核的形式來與文獻裡的結果做比對，再利用Matlab程式語言或Mathematica軟體來模擬二維平面問題及三維空間問題之映像源分佈情形與基本解法中源點強度分佈情況。 (三) 文獻回顧 映像法近年來在工程及數學的應用方面愈來愈受到重視。在2006年，陳與吳[1]使用映像法以及零場積分法搭配分離核來推導格林函數與Poisson積分方程，有別於Greenberg的作法。陳等人[2]利用映像法與加法定理，求解同心圓在各種不同邊界條件情況下的格林函數與其映像源的分佈情形。之後2010年陳等人[3]更利用映像法與退化核推導螺旋差排問題的映像源位置與強度。陳等人[4]延伸映像法求解二維含圓形邊界與三維含球形邊界的邊界值問題，與過去[1,2,3]不同的地方在於[1,2,3]是有外力項的作用，而[4]則是沒有外力項的純邊界值問題。陳等人[5]使用映像法求解偏心圓環映射後的映像源分佈，並討論基本解法中源點最佳位置的分佈。陳等人[6]使用雙極座標推導半無窮平面含一圓洞，以及無窮平面含兩圓洞之映像源映射後的分佈情形。但過去曾有學者Keller [7]，在各種不同形狀的邊界條件下推導格林函數，探討集中力的映像源分佈，並且將格林函數延伸推導至阻抗邊界條件。而Moran [8]探討在水面下深度h之位置，有一個源點，此時速度勢為，考慮其在受Froude number的影響下，速度勢函數對自由表面之水面進行映射後的映像源分佈。在求解此問題中，可將速度勢函數對Froude number進行微擾展開，如同阻抗邊界中的參數及之角色，將其代入自由水面邊界條件，使其滿足邊界條件。因此，我們可以求得速度