第13屆計算數學研討會暨台灣工業與應用數學會年會 Providence University, Taichung City, Taiwan, June 1-2, 2013 第13屆計算數學研討會暨台灣工業與應用數學會年會 (海報論文：碩士生組) Yi-Chuan Kao (高怡絹), 河海工程學系, 國立臺灣海洋大學, 台灣 (tsubasa0522@hotmail.com) Advisor: Prof. Jeng-Tzong Chen (陳正宗) 論積分方程法中二維及三維之真假特徵解 摘要 　　本文利用對偶邊界積分方程法搭配退化核函數來求解二維多連通與三維單一扁長型橢球聲場特徵值問題。為了能夠充份地利用圓形、橢圓及扁長橢球外形的幾何特性，在此採用極座標、橢圓座標與扁長橢球座標來求解此特徵值問題。並將閉合型的基本解以退化核的形式分別以極座標、橢圓座標與扁長橢球座標展開，而邊界物理量則使用特徵函數展開。本文分別在Dirichlet與Neumann邊界條件下，觀察真模態發生時對應領域內為非零場而補領域則為零場，當假模態發生時則內外場域亦是本文探討的議題。最後，我們將以兩例題之真假模態與傳統邊界元素法比較，之後再將其延伸至三維問題探討。 問題描述 橢圓座標的退化核 對於二維(圖一及圖二) 及三維(圖三)之聲場特徵值問題，其聲壓的控制方程式為Helmholtz方程式，表示如下： 考慮的邊界條件為 Dirichlet邊界條件： 或Neumann邊界條件： 零場對偶邊界積分方程法 域內點積分方程 零場積分方程 極座標的退化核 圖一 橢圓形含一圓形孔洞(二維) 扁長橢球座標的退化核函數 圖二 圓形含一橢圓形孔洞(二維) 圖三 橢球體問題(三維) 結果與討論 xy null null plane null null BEM null BEM null xz null plane null null Null-BIEM null 圖八 三維問題下零場與非零場之模態圖 (Dirichlet邊界條件) Null-BIEM null null xy null plane null 圖六 BEM及Null-BIEM之零場與非零場之模態圖 (Dirichlet 邊界條件) 圖四 BEM及Null-BIEM之零場與非零場之模態圖 (Dirichlet 邊界條件) null xz null null plane null BEM null BEM 圖九 三維問題下零場與非零場之模態圖 (Neumann邊界條件) 　　圖四至圖七表示在Dirichlet及Neuamm邊界條件下利用零場邊界積分方程法與邊界元素法，畫出聲場的真假模態圖。藉由這四張圖，我們可以發現在真特徵值的情況下，不管聲場形狀為何，以及邊界條件為何，其領域為非零場，而補領域則為零場。但在假特徵值所對應的模態，其情況是完全不一樣。對於Dirichlet邊界條件，其領域為零場，補領域為非零場，而Neuamm邊界條件則是都為非零場。　　 　　另外，而同樣結果也可以延伸至三維，如圖八及圖九，此兩圖引用J. W. Lee and J. T. Chen [3]，所做出的三維結果。 null Null-BIEM null Null-BIEM null 圖七 BEM及Null-BIEM之零場與非零場之模態圖 (Neumann邊界條件) 圖五 BEM及Null-BIEM之零場與非零場之模態圖 (Neumann邊界條件) 註：kT為真根，kS為假根。 結論 1.使用零場邊界積分方程所做出的數值結果比傳統邊界元素法來的更準確。 2.在真特徵值發生時，領域為零場，補領域為非零場。 3.在假特徵值發生時，領域不一定為零場，補領域不一定為非零場。 參考文獻 [1] J. T. Chen, J. W. Lee, I. L. Chen and P. S. Kuo, On the null and nonzero fields for true and spurious eigenvalues of annular and confocal elliptical membranes, Engineering Analysis with Boundary Elements, 37(2013), pp. 42-59. [2] J. T. Chen, J. W. Lee and S. Y. Leu, Analytical investigation for spurious eigensolutions of multiply-connected membranes containing elliptical boundaries using the dual BIEM, International Journal Solids and Structures, 48(2011) , pp. 729-744. [3] J. W. Lee and J. T. Chen, On the null and nonzero fields for true and spurious eigenvalues of a prolate spheroidal cavity, ICOME2012/JASCOME2012, 12-14 December, Kyoto, Japan, 2012. [4]高怡絹，真假特徵值問題之奇異值分解法探討及Mathematica動畫模擬，國立臺灣海洋大學大學生暑期學習實務體驗計畫，2012。(獲得海大暑期計畫工學院組第三名。) [5]陳正宗、郭柏伸、高怡絹，可以聽出鼓的形狀嗎？－對偶邊界元素法分析，數學傳播，審查中。