边界元素法和保角映射于二维laplace问题的退化尺度之数值与解析研究.pdf

邊界元素法和保角映射於二維 Laplace問題的退化尺度 之數值與解析研究 (一)摘要 本計劃利用邊界元素法和複變函數中的保角映射法來研究 Laplace問 題中不同幾何外形的退化尺度。此退化尺度可由兩個方向來做說明：一為 邊界元素法中 log 核所引起的解不唯一；另一為複變保角轉換中的保角半 徑對應的 logarithmic capacity 。因此，我們將分別利用保角半徑和對偶邊 界元素法來得到退化尺度的解析和數值的結果。藉由使用保角映射法可得 到含 Gamma 函數的解析公式之 Gamma輪廓。此外，我們將用複變函數 解析導得邊界密度函數與使用邊界元素法求得數值離散資料，來驗證兩者 的結果是否一致。並在正常尺度與退化尺度下各別解析推導其內場與外場 的勢能分佈。在對偶邊界元素法架構下，針對弱奇異 [U]和強奇異 [T積分] 算子在奇異值分解法下的左酉向量裡埋藏假訊息是退化尺度與強奇異 [T] 和超強奇異 [M積分算子的右酉向量隱含剛體運動的訊息作一檢驗。為了] 處理非唯一解問題，邊界勢能通量方程式、剛體模態法、 CHEEF 法和超 強奇異邊界積分方程法，將用來提升影響係數矩陣的秩數。最後，幾何外 形的解析和數值的結果將一併探討。 - 1 - (二研究動機與研究問題) 申請人成績優異且具學習熱忱，在 工程數學與流體力學方面表現出強 烈的興趣與上進的態度，且 2011 年學期成績班排名第一。 平常在課堂上所接觸到的工程數學都是著重於數值計算的部分，比較 少去了解數學公式背後所隱含的物理意義以及工程實務上的應用。所以， 希望藉此次參與國科會大專生專題計畫來深入了解邊界元素法的精神以 及使用時所會遭遇的問題 (退化尺度 ) 。 在勢能問題與平面彈性力學中，當欲解問題的幾何尺寸為一特殊尺度 時，利用邊界元素法求解會發生解不唯一的情況 ，此時影響係數矩陣會產 生秩降行為，此即所謂的退化尺度問題。 本計畫擬探討問題之發生機制與解決之道。此外，本計畫也將探 討退化尺度、 Gamma輪廓與 logarithmic capacity 三者之間的關係， 以及驗證文獻中 14種特殊幾何外形 (如表一所示 )退化尺度的解析公 式，並作一數值探討。 - 2 - (三)文獻回顧與探討 在二維 Laplace或彈力問題中，當問題的幾何尺寸為某一特殊尺度 時，使用邊界元素法來求解會發生解不唯一現象 ，我們稱此一特殊尺度為 退化尺度。而退化尺度的發生是由於數值方法所導致的，其一為邊界元素 法中因 log 核所導致的解不唯一行為；另一則是複變函數中保角轉換之保 角半徑產生 logarithmic capacity [1, 2] 為 1時所對應的尺度。 關於退化尺度的研究，在 Laplace問題中， Jaswon 與 Symm [3] 將發 生此退化尺度所對應的輪廓稱為做 Gamma輪廓。然而為何稱為 Gamma 輪廓並沒有多作說明。後來我們發現 Landkof [4] 早於 1966 年已提出 14 種幾何外形發生退化尺度時，所對應的退化尺度解析公式，然而於其書 中，並未見其詳細推導過程。 Constanda [1] 於 1993 年探討二維 Laplace 方程Dirichlet邊界條件的問題。陳等人 [5] 於 2001 年使用循環矩陣來解 析研究圓的退化尺度，隔年 [6]則延伸探討同心圓環的退化尺度問題。 Hansen 與 Christiansen [7] 於 1985 年研究圓、橢圓與線段的退化尺度問 題，並用奇異值分解法 (SVD)來找出退化尺度發生的位置。陳等人 [8] 於 2005 年時則探討關於正三角形、矩形等不同幾何外形的退化尺度，而在 文中也使用奇異值分解法(SVD)來找出退化尺度的發生位置，並針對秩數 不足的問題，提出了三種補救的方法，剛體模態法、CHEEF (Combin