某些正反散射问题数值算法的研究

发布时间：2018-01-16 23:28

本文关键词：某些正反散射问题数值算法的研究　出处：《吉林大学》2015年博士论文　论文类型：学位论文

【摘要】：声波和电磁波的散射与反散射理论在数学物理领域发挥了重要作用,并且在实际生活中有着很多应用,如医学成像,雷达探测,无损探伤等.本文主要讨论了障碍散射问题和单缝散射问题的数值求解及贝叶斯方法在内腔体反散射问题和单缝反散射问题中的应用.我们所考虑的散射问题模型均为Helmholtz方程.由于障碍散射问题和单缝散射问题都是无界域上的散射问题,进行数值求解时首先要考虑无界域截断问题.在本文中,我们利用Dirichlet-to-Neumann算子(简称DtN算子)将无界域问题转化为有界域问题.与常用的吸收边界条件相比；这种转化是精确的,没有近似.因此,在第一章的绪论中,我们介绍了Hehmholtz方程和 Dirichlet-to-Neumann算子及文中所要探讨的几种散射问题的研究现状,并列出了本文的主要工作.第二章中我们提出了一种处理DtN边界条件的间断有限元方法,以障碍散射问题为模型介绍了算法的具体格式,证明了显波数k的DG模和L2模误差估计,并给出了相应的数值实验结果.第三章中我们利用DtN算子将单缝散射(single slit scattering)问题转化为单缝处的一个算子方程,证明了算子方程解的存在性、唯一性及稳定性.应用Galerkin法对算子方程进行数值求解,证明了数值解的存在性、唯一性和收敛性.最后进行了数值模拟,其结果与实际的物理现象一致,说明了算法的有效性.第四章中介绍了两个反散射问题的贝叶斯方法.以内腔体反散射问题为例,给出了贝叶斯方法的具体步骤及方法适定性分析,数值实验结果说明了方法的可行性.同时,我们还利用贝叶斯方法反演了单缝散射问题中点源的位置,给出了相应的数值算例.本文的最后一章为总结,对全文内容进行了总结并列出了后续研究的几个方向.1.处理DtN边界条件的间断有限元方法考虑障碍散射问题,设区域Ω (?)R2为Lipschitz区域,我们用表示区域Ω的外部.设Ω为包含原点的星形域且边界(?)Ω充分光滑,例如,(?)Ω是C2类的.我们考虑如下的具有Dirichlet边界条件和Sommerfeld辐射条件的Helmholtz方程外问题：给定入射场ui,求散射场us使得对于任意的函数f∈H1/2((?)BR),定义DtN算子如下：借助DtN算子我们将问题(1)转化为如下有界域上的问题其中上的单位外法向量,为了方便,我们将us简写成u.这里,BR1和BR2为包含区域Ω的两个大圆,将问题(1)转化为有界域上的问题(3)后,计算区域为ΩR2.为了处理DtN边界条件我们需要在边界ΓR2附近剖分出一个环形区域AR1R2,在其内部使用环形网格剖分和谱方法基底.在剩余的部分,即ΩR1内采用标准的三角网格剖分和多项式基底.记ΩR2上的整体剖分T=T1+T2,其中T1表示ΩR1中的三角剖分,T2表示AR1R2中的环形剖分.对于任意的单元K∈T1,我们定义其直径若单元我们定义其直径为圆环的宽度相应地,对于剖分T1中的三角单元K的边e,定义he为其长度he：=diam e.在我们的方法中采用拟一致网格剖分,即存在常数C10,C2 0,使得设Fh表示剖分T1的“骨架”,即跃度与平均值的定义与其它标准DG方法相同.定义基于剖分T的分片Sobolev空间：其中那么问题(3)的变分形式如下：求u∈X使得其中ah(·,·)是如下定义的DG双线性形式：且G(v)的定义如下：这里a,b和λ是待定的正参数.从上面的定义可以看出,如果u是问题(3)的解,那么u满足下式定义DG模如下：对于单元K∈T1,无论它是直边三角形还是曲边三角形,在计算中我们都可以将它映射为参考单元K,其顶点为(0,0),(0,1),(1,0).用P1(K)表示单元K上的全体线性多项式的集合.我们定义DG近似空间其中则Vh(?)X是有限维空间.下面我们给出基于变分形式(4)的DG格式：求uh∈Vh使得可以证明DG格式(6)的解是存在且唯一的.定理1.(存在唯一性定理)设则DG格式(6)存在唯一解.为了给出误差估计,我们首先证明了空间Vh的逼近性质及双线性形式定理2.设函数对于(?)φh∈Vh2,如下的迹逆不等式成立其中Gti是与h无关的正的常数,这里“ti”表示“trace inverse",定理3.设ah(·,·)定义式(5)中待定的常数a，，b，λ分别满足则双线性形式bh(·,·)满足下式,|bh(v,v)|≥Ccope,‖v‖DG2,(?)v∈Vh,其中Ccoer是与波数k和剖分h无关的正的常数.定理4.存在与波数k无关的常数C0，使得借助空间Vh的逼近性质及双线性形式bh(·,·)的性质,我们证明了DG模和L2模的误差估计如下：定理5.(DG-模误差估计)设区域Ω为光滑星形域,边界(?)Ω是解析的,定理3的条件成立,当k2h充分小时,我们可以得到如下的误差估计其中C0是与网格剖分h和波数k无关的常数.定理6.(L2模误差估计)设定理5的条件满足,当k2h充分小吐我们可以得到如下的误差估计其中C0是与网格剖分h和波数k无关的常数.最后,我们应用所提出的间断有限元方法对障碍散射问题进行了数值求解.实验结果说明我们的方法对于大波数情形也是有效的,并且其数值误差收敛阶与前面的理论证明相吻合.2.单缝散射问题的Galerkin法设ΓUΓc表示R2中带有单缝的良导体平面,其所在直线设为x轴,其中表示平面上的单缝,}＼r表示良导体平面.令单缝的左端点为原点,单缝长度为L.ui为入射场,ur为平面产生的反射场,即y=0时,ui(x,y)+ur(x,y)=0.若入射场是入射角为θ(关于y轴)的平面波,那么其中当y≥0时,设参考场全场这里us表示散射场.当y0时,全场由于Γc为良导体平面,那么由其物理性质可知在Γc上ut=0.于是全场ut满足：其中波数由于入射场ui和反射场ur都满足Helmhotz方程,且在故散射场us满足：对(10)关于x做Fourier变换可以解得：当y0时,当y0时,其中为零扩张算子,定义如下：定义Γ上的Dirichlet-to-Neumann算子T如下：利用Γ上的DtN悄算子T和全场ut及其法向导数(?)ut/(?)y在Γ上的连续性,我们可以将问题(10)-(12)化成如下r上的算子方程：其中于(16)中解得uΓs后,可由(13),(14)式得到散射场us.为了简单,我们记uΓs为u.于是方程(16)可以改写为：对于实数s,定义Sobolev空间：相应的模为当不发生混淆时,用分别代替借助上面的空间定义,可以证明方程(17)的解存在、唯一并且稳定.其中C0为与g无关的常数.设为H*2(Γ)的有限维子空间,且形成空间VN的一组基底.则解方程(17)的Galerkin法为：求使得可以证明数值解un。存在、唯一并且收敛到真解u.定理8.方程(18)的解存在且唯一定理9.设u是方程(17)的解,un是方程(18)的解,且在中终归稠密,即则当N充分大时,下面的估计式成立,其中C0是与u，N无关的常数.从而在我们所考虑的单缝散射问题中,故所以,我们选取不失一般性,我们假设则算子方程(17)相应的Galerkin方程可化为：解上面的方程可以得到uΓs的数值解通过uΓs的数值解即可得散射场us(x,y)的数值解.最后,我们证明了{VN}N∞=在空间H*1/2(r)中终归稠密.定理10.有限维子空间列{VN}N=1∞在H*1/2(Β)中终归稠密.我们应用Galerkin法对单缝散射问题进行了数值模拟,其结果与实际的物理现象相吻合,说明了算法的有效性.3.两个反散射问题的贝叶斯方法考虑下面的问题：已知变量y∈RJ,求q∈Rn满足如下方程这里我们称y为观测值,q为未知量.这是一个典型的反问题模型,通常情况下观测值y是带有误差的,即我们实际考虑的模型应该是其中δ∈RJ为观测误差.应用贝叶斯方法从统计学的角度来思考以上问题.首先将变量q,y和δ均视为随机变量,那么反问题的解就是给定y,求变量q的概率分布,即求条件分布q|y.这样虽然不知道误差δ的具体值,我们却可以将它的统计性质写入到了模型中.设随机变量q的概率密度函数为π0(q),随机变量δ与q相互独立,且δ的概率密度函数为π(δ).当q已知时,y由(19)式决定,故y|q的概率密度函数为)).于是随机变量y和q的联合概率密度为由贝叶斯公式,我们就可以得到条件分布q|y的概率密度函数.π0(q)通常称为先验密度函数,称为似然函数,πy(q)称为后验密度函数.设脚是具有密度函数πy的后验测度,μ0是具有密度函数π0的先验测度,那么.我们称dμy/dμ0为Radon-Nikodym导数,它关于似然函数成正比,即利用Radon-Nikodym导数的概念可以将上面的分析推广到无穷维的情形,即q∈X,y∈Y,X和Y均为无穷维空间.由上面的分析可以看出,应用贝叶斯方法解反问题主要可以分为以下三步：1.依据关于未知量q的所有先验信息来确定先验密度函数π0(q).2.根据误差函数的性质来确定似然函数π(y-g(q)),实际上似然函数反映了当未知量q已知时观测值y所满足的分布.3.根据前两步的准备和贝叶斯公式,可以确定后验密度函数的表达形式,然后寻找合适的抽样方法将其描述出来即可.3.1.内腔体反散射问题的贝叶斯方法考虑如下的内腔体散射问题,设D(?)R2是单连通区域且具有C2光滑的边界.设点源和观测点都位于区域D内部的曲线C上,则内腔体散射问题可以描述如下：其中k0是波数,ui为入射场,定义如下：这里,是C上的一个固定点,Φ(x,d)是二维Helmhotlz方程的基本解.边界条件(22)中的B代表三类不同的边界条件,即Dirichlet, Neumann 和 Robin边界条件.我们所要考虑的反问题就是已知测量曲线C上的散射场us及边界条件(22),重构腔体形状D.下面为了简化问题,我们设区域D是关于原点的星形域,即其中0r R0.为了后面证明的方便设q=lnr.故问题的模型可以写成如下形式其中为观测点的个数)是对应于散射场方程(21)的有限维观测算子,向量y是具有误差δ的观测值.这里q属于函数空间X,对于X的选择将在后面介绍.假设q满足正态分布(也称高斯分布)其中m0为期望,G0为协方差算子.误差δ满足正态分布N(0,Γ),其中Γ为有界协方差矩阵.由前面的分析可知Radon-Nikodym导数参照文献[99]我们给出q所满足的正态分布.设算子且定义域如下：由于星形域的表示中,r(θ)是2π周期的,所以这里我们在函数空间的符号中用方括号表示函数是2π周期的.假设q"(θ)满足L2[0,2π]上的正态分布N(0,A-1).则由Karhunen-Loeve展开,我们可以得到其中an和bn是独立同分布的随机变量,其满足的分布为Ⅳ(0,1).然后我们通过对q"(θ)积分来得到q(θ).而积分得到的结果并不是唯一的,我们选择其中的一种周期形式来定义q(θ),如文献[99]所示,定义其中q0～N(μ,σ2)是一个正态随机变量.对q"(θ)的展开式逐项积分可得由积分算子的线性和连续性可知这样得到的q(θ)依然是正态分布.根据文献[95]中的引理6.25,可知上面定义的q"(θ)对于α1/2几乎一定落在中.因此上面的逐项积分是可以进行的,且g(θ)几乎一定是C2,a的,因此设可以验证观测算子g满足如下的假设1：假设1(i)对于任意ε0,存在常数使得对于q∈X下式成立,(ii)对于任意t0,存在常数K=K(t)0,使得对于满足t的q1,q2∈X,下式成立,则根据文献[95]中的定理4.1,定理4.2和引理6.31可以得到如下的适定性定理.定理11.(适定性定理)设观测算子g满足假设1且则有(i)后验测度μg关于先验测度μ0绝对连续,且Radon-Nikodym导数(ii)后验测度μg是L2[0,2π]上的一个概率测度.(iii)后验测度μy关于观测值y是Lipschitz连续的,即存在常数使得其中这里，Hellinger距离定义如下:若μ1和μ2都关于μ0绝对连续,则3.2.单缝反散射问题的贝叶斯方法下面考虑点源入射时的单缝散射反问题：假设已知下半平面一条线上的散射场值,目标是反演点源的位置.单缝r的长度仍然设为π,已知y=a上若干点的散射场测量值其中δ为测量误差,由单缝散射问题中的方程(10)-(12)决定,目标是反演点源的位置q=(x,y).假设我们事先已知点源的大致位置在[b,c]×[e,f]的区域中,这样初始分布可以采用均匀分布依然假设误差δ满足正态分布Ⅳ(0,σ2I),则根据贝叶斯公式,当时,后验分布的密度函数正比于当(x,y)不在区域[b,c]×[e,f]中时,后验分布密度函数为0.利用MCMC方法抽样描述上述关于q的后验分布,取平均值后即可做为点源位置的估计.我们对内腔体散射中重构腔体形状的问题及单缝散射中反演点源位置的问题,应用贝叶斯方法分别进行了数值实验,实验说明了贝叶斯方法的有效性.
[Abstract]:......
【学位授予单位】：吉林大学
【学位级别】：博士
【学位授予年份】：2015
【分类号】：TN011

【相似文献】