二维圆柱外域Helmholtz场虚拟频率产生机制之探讨.pdfVIP

二維圓柱外域Helmholtz場虛擬頻率產生機制之探討 陳義麟 國立台灣海洋大學河海工程研究所博士班研究生 國立高雄海洋技術學院造船工程科講師 摘 要 本文主要目的係以對偶邊界積分方程為基礎，探討虛擬頻率發生的機制，並以對偶級 數模式理論推導虛擬頻率所發生的位置，且採用直接法的對偶邊界元素法作數值實驗來驗 証此理論。研究結果顯示，虛擬頻率發生的位置與使用的數值方法中的核函數有關，而與 該外域問題所給定的邊界條件型式無關，並以二維Helmholtz方程為範例驗證此一結論。 用一個合成式，在含退化邊界問題的應用上，仍然 一.前 言 有其盲點。而對偶邊界元素法使用兩式，可同時解 當一單自由度振動系統承受含自然頻率之外力激 決退化邊界與虛擬頻率的問題，有關對偶積分架構 發時，會產生振動放大的現象，此種現象稱為共 在計算力學所扮演的角色，陳與洪有一回顧性文章 振，此為自然發生的物理現象，在數學與物理上均 之中有詳細討論(陳與洪，1998)。若改以有限元素 可反應此事實。然而，於積分方程解外域荷姆茲 法分析外域聲場問題，雖無所謂的解不唯一的虛擬 (Helmholtz)場問題時，會發生數值不穩定的虛擬頻 頻率問題，但為了滿足遠域輻射條件，需引入無限 率(fictitious frequency) ，又可稱為寄生頻率(parasite 元素的觀念，DtN(Dirichlet to Neumann)法(Givoli et. frequency) 或 稱 之 為 不 規 則 頻 率 (irregular al, 1997)即為一例，Malenica ＆ Chen(1998)利用 frequency) ，此非自然現象，乃是因採取的數值方 改變對應虛擬內域的邊界來確保所有頻率下解的 法過程中所衍生的問題，因為在純解析推導時並不 唯一性。因虛擬頻率係伴隨在我們所使用的積分方 產生此種問題(陳與洪，1992)。為了克服此種問題， 程中，故在陳與洪(1992)的書中稱虛擬頻率為寄生 Schenck (1968) 最 早 使 用 CHIEF (Combined 的數值共振頻率，而非實際的物理問題。此點更可 Helmholtz Integral Equation Formulation)法，以域內 由聲波、水波、彈性波、電磁波、浮體工程、置入 點邊界積分方程式為輔助條件，補足不夠的束制條 物(inclusion)波場、聲波結構互制、流體結構互制、 件。此法所選取的點若不巧為內域問題模態的結點 土壤結構互制(李，1992)與聲場全像術等領域的研 (node)時則將會失效，亦即得到的是無效的方程式 究均會遭遇到，且可看出問題的數學本質是一致 (trivial equation) ，此問題在愈高頻時發生的可能性 的。而其數值觀點則可視為 0/0 的數值不定型。雖 就愈高(Seybert ＆ Rengarajan 1989) 。後來 Burton 然理論上可使用 L’Hospital 法則予以克服，但是數 ＆ Miller (1971)及 Lee ＆ Sclavounos (1989) ，亦 值執行上卻有其困難。 曾提出合成奇異積分方程式與超奇異積分方程式 文獻上對於虛擬頻率所發生位置，Shaw(1979) 的純虛數倍的方法，可在任意波數下均可得唯一解 曾提出：“外域的 Dirichlet 問題所產生的虛擬頻率 (Cunefare ＆ Koopmann 1989) ，此技巧可說是將波 會對應到內域 Neumann 問題的共振頻率，而外域