同心圆环二维helmholtz方程特征值问题-海洋大学.docVIP

同心圓環二維Helmholtz方程特徵值問題 真假根探討 郭世榮 陳正宗 劉孟龍 國立台灣海洋大學河海工程系，基隆 摘要 本文針對同心圓環之Helmholtz特徵值問題利用實部對偶邊界元素法來求得真假特徵值。在解析解方面，以離散系統模擬連續系統，並利用邊界元素法中影響係數矩陣循環的特性來求得； 一. 前言 以往求解Helmholtz方程的特徵問題時，多半以複數基本解處理，傳統上邊界元素法為避免在複數域下求解，可將Helmholtz方程式中含特徵值項視為Laplace方程式之外力源，但是此外力源會導得一體積分項，因而使得在利用邊界元素法求解時，仍需對內域作離散。多倒易法將此體積分轉成一系列的邊界積分[1,2,3]。後來Kamiya [4]等人與葉、陳等人[5]發現多倒易法和複數邊界元素法只取實部型，並無兩樣。陳等[6,7]已成功地發展實部對偶邊界元素法，然其應用均侷限於求解單連通的問題，如圓形[8]或正方形[9]等。本文將利用實部對偶邊界元素法[10,11]，建立二維同心圓特徵值問題對偶積分推導，此為多連通問題，其假根機制將更為複雜。本文同時進行對偶邊界積分解析推導與數值計算，以及求得二維同心圓環的自然頻率與自然模態，並配合奇 異值分解技巧，以對偶架構來判定真假特徵解。 二. 研究方法 2.1 對偶邊界元素法 對於傳統多倒易法或實部邊界元素法之假根現象，除了利用對偶架構以殘值技巧過濾外，我們也可以把多倒易法或實部邊界元素法中的奇異方程和超奇異方程結合後再使用奇異值分解法，能更有效地將假根檢測出來並濾除之[7,12,13, 14]。本文將以解析方式加以說明多連通問題的假根現象，進而建構一套適用於求解任意二維聲場 (單連通或多連通)之聲頻與聲模的程式。 考慮二維特徵值問題，其控制方程式為: 其中，k為波數，u為場量。 根據邊界元素法U、T、L與M四個核函數的對偶架構，可將場解表示為[15,16,17]: 其中C.P.V.和R.P.V.分別表示柯西主值與黎曼主值，而H.P.V.為阿達馬主值。若僅取核函數的實 部，則上二式可重新表示成： 上二式中，下標R表示取核函數的實部。將上二式做積分方程離散化，並採用常數元素，對 Dirichlet問題而言，可得: 上式(6)、(7)式即成為一般我們所處理的非線性特徵值問題，不過由於少了虛部的束制條件，將導致假根的現象，所以我們會從式子(6)中得到一些特徵值，不過它們並不會滿足式子(7)，同理我們也會從式子(7)得到一些特徵值，他們也不會滿 足式子(6)，此即為殘值技巧的精神。 本文採用奇異值分解法，來判定真假根，我們可以先結合(6)與(7)式，使得: 其中， 對於真根而言，矩陣[C(k)]的秩原本為N降為N(1 以求得非零解。所謂奇異值分解法，說明如下: 一過定矩陣[A]m(n，設m是方程式的數量，n 是未知數的數目： S.V.D.會把[A]m(n分解成： 其中，， 表示共軛轉置，即所謂奇異值，且，當我們有p個奇異值等於零時 (0(p(n)，也就是說矩陣[A]的秩將降為n(p，對真的特徵值而言，若秩將降為n(p時，即表示此根重根p次，可是若為假根的話，矩陣[A]的秩 仍是n，也就是其最小奇異值不為零。 當判斷出是真根之後，若要尋求特徵向量可由V矩陣對應的柱向量決定。利用SVD技巧除了有效率地判斷真假根之外，還有一好處能判定特 徵值的重根數。此方面應用，可參考[6,7]。 2.2 圓形Helmholtz特徵值問題 若將圓形邊界等分後，如圖1所示，則影響 係數矩陣[UR]與[]之計算，說明如後： = = = 圖1 影響係數與配置點關係和積分示意圖 我們進一步發現[]或[](以後均以[A]與[B]表示)，是一對稱的循環矩陣(Circulant Matrix)，可表示如下： 其中， 依據此特性，(11)式可以表示成： 其中，為最基本的循環矩陣，其特性如下： 而此循環矩陣的特徵值為的 根，即為： 而特徵向量可由特徵值表成： 另外[A]矩陣亦可以被對角化，其過程如下： 其中[R]是[A]之特徵向量所組成，而由於[A]和的特徵向量是相同。因此，可以求出[R]矩 陣，可由式所組合而成： 而對角化後對我們在預測真假根上會有所幫助。 2.3 單一圓環真假根預測 首先若僅取Helmholtz方程的基本解中U核函數的實部，可得： 由矩陣理