解析及数值实验探讨积分方程外域问题CHIEF方法.docVIP

解析及數值實驗探討積分方程外域問題CHIEF方法 陳義麟 朱皖山 陳正宗 梁明德 郭世榮 摘要 積分方程用於解外域聲場Helmholtz問題由來已久。如果僅使用奇異積分方程（式）會產生虛擬頻率的非物理現象。CHIEF法因易於使用而受到歡迎，且可以克服外域虛擬頻率的問題；但是如果所選取內域點配置不當則會導致失效。本文主要目的係以CHIEF方法為基礎，以圓形循環矩陣（circulant）的特性及可分離核函數（degenerate kernel function）配合SVD的技巧，用解析的方法詳細探討虛擬頻率發生的機制及過濾虛擬頻率時可能的失敗的位置，並且以圓柱外域問題為例，數值實驗結果與解析的結果亦非常吻合。 前言 於積分方程求解外域之散射或輻射問題時，當入射波頻率或輻射頻率與對應之內域相關頻率相同時，則其數值結果將不唯一並呈現不穩定而偏離正解，我們可稱此為〝數值共振〞（陳與洪1992）。此非自然現象，而是因採取的數值方法過程中所衍生的問題，此頻率寄生於所使用的方法中，可稱為寄生頻率，一般文獻稱為虛擬頻率（fictitious frequency），此乃因純解析推導不會有此問題的。在以頻率域求解聲波及電磁波或彈性波問題時，亦有類似現象。 由於在實際解析推導過程中並無虛擬頻率的問題，此因數值處理才產生的虛擬頻率，其實是等於在處理一不定型的極限問題，因未定係數僅含分子部分，無法反映出此存在的有限值，且因為要造成分母為零的頻率為無理數，於電腦運算時，雖不會有數值無效情形，但於虛擬頻率附近其數值敏感度很大，易造成偏離正解的誤差。Burton Miller（1971）及Lee Sclavounos（1989），提出合成奇異積分方程式與超奇異積分方程式的純虛數倍的方法，此方法有對偶邊界元素法的架構，是以奇異積分方程加上超強奇異積分方程乘上一個虛數常數的組合式。陳等（1999）提出對偶邊界元素法，此法是將奇異積分方程得到的值與超強奇異積分方程得到的值予以比較，可有效的過濾虛擬頻率，對偶邊界元素法更可有效的解含退化邊界的問題。Schenck（1968），Benthien Schenck（1997）則使用CHIEF（Combined Helmholtz Integral Equation Formulation）法，以域內點邊界積分方程式為輔助條件，補足不夠的束制條件。此法所選取的點若不巧為內域問題模態的節點（node）時則將會失效，亦即得到的是無效的方程式，此問題在愈高頻時發生的可能性就愈高（Seybert ＆ Rengarajan（1989））。而近年來SVD已被發展成為線性系統之最重要的工具，Francis（1989）應用此法來解電磁場的共振問題，陳等（1999）、陳等（1999）及葉等（1998、1999）亦應用此法有效的過濾內域問題的假特徵值(spurious eigenvalues)。Juhl（1994）與Polin（1997）將CHIEF配合SVD技巧更可有效的過濾虛擬頻率。 本文主要目的係以CHIEF方法為基礎，以圓形循環矩陣（circulant）的特性及可分離核函數（degenerate kernel function）配合SVD的技巧，用解析的方法詳細探討虛擬頻率發生的機制及過濾虛擬頻率時可能失敗的位置，並配合數值實驗探討此問題。 二.理論推導 線性化聲波方程式的控制方程式為 （1） 其中為速度勢位，為問題的定義域，為聲速，為聲源項。若無聲源項，則在頻率域的控制方程式（1）可改寫成 （2） 其中為的Fourier轉換，k為波數表為，其中為角頻率。 定義基本解滿足下式 （3） 其中為Dirac delta 函數。引入格林第三定理，可得到域內點（外域問題）積分式如下： （4） 其中為一常數（和的維度有關），B為邊界，為外域的領域，，。 x推到平滑邊界時，則可得邊界點的積分方程如下： （5） 其中R.P.V.代表黎曼主值（Riemann Principal Value），C.P.V.代表柯西主值（Cauchy Principal Value），而U及T為積分方程的二個核函數，其分別為 （6） （7） 其中，表示第一類n階Hankel函數，為場點與源點之間的距離，為源點的法向量之第i個分量，而，為場點的法向量之第i個分量。將式（5）、邊界離散化可得： （8） [U]及[T]矩陣分別為U和T核函數所產生之對應影響係數矩陣，細部推導可在