發(fā)布時間：2018-02-16 13:13

  本文關鍵詞： 聲波傳輸問題 DtN算子 有限元方法 傳輸特征值問題 C~0內罰間斷伽遼金法　出處：《重慶大學》2016年博士論文　論文類型：學位論文

【摘要】：時諧聲波傳輸問題和傳輸特征值問題在實際科學和工程領域都有廣泛的應用。傳輸特征值能用來估計散射體材料的性質,并且在逆散射理論中對于證明解的唯一性和重構有著重要的作用。本論文主要研究聲波傳輸問題和傳輸特征值問題的數值方法。本文主要由兩部分組成:第一部分關于聲波傳輸問題的數值方法;第二部分關于傳輸特征值問題的數值方法。對于聲波傳輸問題,最常用的處理無界區(qū)域的方法是邊界積分方程方法及其與有限元方法的耦合,比如Hsiao和Xu用邊界積分方程方法處理無界區(qū)域的計算(Hsiao and Xu,2011)。為了應用有限元方法求解此類問題,最常用的方法是通過引入一個包含障礙物的人工邊界將無界區(qū)域分解為一個有界區(qū)域和一個無界區(qū)域。我們根據人工邊界外聲波滿足的散射問題,分別利用傅里葉級數和邊界積分算子定義了兩個不同的Dirichlet-to-Neumann(DtN)算子。接著,通過在人工邊界上引入DtN算子,可以將原無界傳輸問題轉化為等價的非局部邊值問題。在適當的索伯列夫空間下,證明了相應的變分問題解的存在唯一性。對于聲波傳輸特征值問題,為了避免直接計算非自伴特征值問題,我們將原傳輸特征值問題轉化為等價的一系列自伴的四階特征值問題。接著,提出了用拉格朗日元的C~0內部懲罰間斷伽遼金(C~0 IPG)方法去研究四階特征值問題。由于低階項對離散算子范數收斂性的影響,不能直接應用Babu?ka-Osborn理論證明收斂性。為了克服這個困難,根據Descloux等人在(Descloux et al,1978a)中發(fā)展起來的收斂性理論以及Antonietti等人用DG方法求Laplace特征值問題的思想(Antonietti et al,2006),我們首先證明了C~0 IPG方法是譜正確的,接著證明了其最優(yōu)收斂性。而對于原傳輸特征值問題轉化的一個等價的非自伴的四階特征值問題,我們給出了離散四階傳輸特征值問題的C0 IPG方法,并證明了C~0 IPG方法最優(yōu)收斂階。對于高階橢圓問題,相比于經典的協(xié)調有限元方法,C0 IPG方法的數值實現更簡單,因為它的基函數簡單并且有更少的自由度。對于每一種提出的數值方法,我們給出一些數值例子來驗證方法的有效性和準確性。
[Abstract]:The time-harmonic acoustic wave propagation problem and the propagation eigenvalue problem are widely used in practical scientific and engineering fields. The transmission eigenvalues can be used to estimate the properties of scattering materials. And in the inverse scattering theory, it plays an important role in proving the uniqueness and reconstruction of the solution. In this paper, we mainly study the numerical methods of acoustic wave transmission problem and transmission eigenvalue problem. This paper mainly consists of two parts: the first part. The numerical method of acoustic wave transmission; The second part deals with the numerical method of transmitting eigenvalue problem. For acoustic wave transmission problem, the most commonly used method to deal with the unbounded region is the boundary integral equation method and its coupling with the finite element method. For example, Hsiao and Xu use the boundary integral equation method to deal with the computation of unbounded region. The most commonly used method is to decompose the unbounded region into a bounded region and an unbounded region by introducing an artificial boundary containing an obstacle. Two different Dirichlet-to-Neumann-DtN operators are defined by using Fourier series and boundary integral operators, respectively. Then, by introducing DtN operators on artificial boundaries, The original unbounded transmission problem can be transformed into an equivalent nonlocal boundary value problem. The existence and uniqueness of the solution of the corresponding variational problem are proved in the appropriate Soberlev space. In order to avoid directly computing the non-self-adjoint eigenvalue problem, we transform the original transmission eigenvalue problem into an equivalent series of self-adjoint fourth-order eigenvalue problems. In this paper, the fourth order eigenvalue problem is studied by using the Lagrangian C0 internal penalty discontinuous Galerkin C0 IPG method. Because of the influence of low order terms on the convergence of discrete operator norm, Babu? Ka-Osborn theory proves convergence. In order to overcome this difficulty, according to the convergence theory developed by Descloux et al in Descloux et al1978a and the idea of Antonietti et al. to solve Laplace eigenvalue problem by DG method, we first prove that the C0 IPG method is spectral correct. Then we prove its optimal convergence. For an equivalent non-self-adjoint fourth-order eigenvalue problem transformed from the original transmission eigenvalue problem, we give a C0 IPG method for discrete fourth-order transmission eigenvalue problem. It is proved that the optimal convergence order of the C0 IPG method is simpler than that of the classical concordant finite element method (C0 IPG) for higher order elliptic problems. Because its basis function is simple and has less degree of freedom, we give some numerical examples to verify the validity and accuracy of the proposed numerical method.
【學位授予單位】：重慶大學
【學位級別】：博士
【學位授予年份】：2016
【分類號】：O241.82

【相似文獻】

相關期刊論文 前10條

1 解惠青;戴華;;多參數二次特征值問題重特征值的靈敏度分析[J];計算數學;2006年01期

2 王琳;;用Fourier Galerkin方法求解間斷特征值問題[J];數學理論與應用;2006年02期

3 高瑞華;孫惠娟;;一類非線性特征值的存在性[J];河南科技;2010年16期

4 欒世霞;孫欽福;;一類二階半正特征值問題的正解[J];應用數學學報;2011年01期

5 李國發(fā);;一類二階兩點邊值特征值問題的非平凡解[J];曲靖師范學院學報;2011年06期

6 嚴深海;;雙特征值約束下的兩類逆二次特征值問題[J];江西理工大學學報;2012年05期

7 程藝;一類離散特征值問題的反問題[J];中國科學技術大學學報;1985年04期

8 顧祝全;含譜參數高次冪的Zakharov-Shabat特征值問題的跡公式[J];數學物理學報;1991年03期

9 喬志軍;WKI族的換位表示[J];科學通報;1992年08期

10 和福生;唐民英;;參數矩陣虛特征值的判別及應用[J];云南師范大學學報(自然科學版);1993年02期

相關會議論文 前5條

1 陳正宗;;使用基本解法求解特征值問題[A];中國力學學會學術大會'2005論文摘要集（下）[C];2005年

2 袁明武;陳璞;孫樹立;孟昭龍;;有限元快速稀疏靜力及特征值問題求解器[A];第十一屆全國結構工程學術會議論文集第Ⅰ卷[C];2002年

3 李紅云;沈為平;;大型結構特征值問題的混合粒度并行算法[A];西部大開發(fā) 科教先行與可持續(xù)發(fā)展——中國科協(xié)2000年學術年會文集[C];2000年

4 吳志剛;鐘萬勰;;最優(yōu)H_∞范數及Hamiltonian微分方程邊值問題特征值計算(Ⅱ)[A];第二十屆中國控制會議論文集（上）[C];2001年

5 黃斌;巫文君;張海洋;;隨機結構特征值遞推求解方法的改進[A];隨機振動理論與應用新進展——第六屆全國隨機振動理論與應用學術會議論文摘要集[C];2008年

相關博士學位論文 前10條

1 曾令輝;二次特征值問題向后穩(wěn)定算法[D];復旦大學;2014年

2 何軍;張量的譜理論和數值代數幾個問題的迭代解法[D];電子科技大學;2015年

3 來鑫;一類趨化模型和一類經理期權模型的定性性質[D];哈爾濱工業(yè)大學;2015年

4 耿宏瑞;聲波傳輸問題和傳輸特征值問題的數值方法研究[D];重慶大學;2016年

5 黃鑫;低秩阻尼非線性特征值問題的理論和算法[D];復旦大學;2011年

6 李慶春;矩陣相合型特征值[D];吉林大學;2007年

7 張磊;高階系統(tǒng)狀態(tài)反饋控制部分特征值配置[D];哈爾濱工業(yè)大學;2014年

8 王順緒;特征值問題的并行計算[D];南京航空航天大學;2008年

9 吳春紅;幾類矩陣的逆特征值問題[D];廈門大學;2009年

10 魏瑩;若干代數特征值反問題[D];南京航空航天大學;2015年

相關碩士學位論文 前10條

1 閆聚品;一類能量依賴于速度的特征值問題及其Bargmann約束下的可積性[D];河北工業(yè)大學;2015年

2 楊樂樂;一個三階特征值問題及其對應的Bargmann完全可積系[D];河北工業(yè)大學;2015年

3 崔鳳娟;帶脈沖條件的不定S-L問題的非實特征值[D];山東大學;2016年

4 董昭婷;能量依賴于速度的三階特征值問題的可積系[D];石家莊鐵道大學;2016年

5 郭景艷;能量依賴于速度的二階特征值問題及其可積系統(tǒng)[D];石家莊鐵道大學;2016年

6 劉貴梅;張量特征值的幾個新包含集[D];云南大學;2016年

7 游園藝;張量的秩和特征值的若干問題[D];哈爾濱工業(yè)大學;2016年

8 孫付;微分算子不定譜問題在不同條件下非實特征值的估計[D];曲阜師范大學;2016年

9 王靜;二階半線性微分方程周期特征值問題的譜[D];西北師范大學;2016年

10 高婷;四階線性不定權邊值問題非實特征值的研究[D];西北師范大學;2016年

，

本文編號：1515599