计算开腔体散射有限元与边界积分对称耦合方法

来源 :吉林大学 | 被引量 : 0次 | 上传用户:zw840909
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
腔体电磁散射与反散射问题的分析在近些年的理论与应用研究中都是很重要的.正散射问题是在已知入射场和腔体的形状信息,预测远离腔体的场的分布;反散射问题则是已知一些人工的测量数据,通过对这些数据的分析与计算,推测出腔体的形状.本文考虑嵌入在一个良导体的无限接地平面的具任意形状且填充的腔体的时谐平面波散射问题.文章给出一种在TE极化和TM极化情形下求解电磁散射问题的有限元和边界积分方程的对称耦合方法.本文基于对称耦合方法讨论了正散射问题和反散射问题,提出了正散射问题的变分公式,研究其弱解的存在性和唯一性,并给出区域导数,对于反散射问题给出了解的唯一性和局部稳定性结果.首先研究TE极化情形的正反散射问题.设入射平面波ui=exp(iαx1iβx2)从上方入射到良电导体表面Γg∪S,其中α=κ0 sinθ,β=k0 Cosθ,θ∈(-π/2,π/2)是关于正x2轴的入射角κ0=ω(?)为自由空间的波数.散射场满足△us+k2us=-(k2-k02)uref于Γg∪S上方. (1)us=-uref于Γg∪S上. (2)此外,散射场必须满足辐射条件对于正散射问题.有引理1正散射问题(1)一(3)至多有一个解.对于R2中的边界为aΩ的有界区域Ω,Tc(?)ΩHs(Ω)和Hs(T)分别为通常的Sobolev空间,其范数分别为‖·‖Hs(Ω)和‖·‖H8(T).定义下列空间:L2(Γ):={u|Γ:u∈L2((?)Ω)}, H1/2(Γ):={u|Γ:u∈H1/2((?)Ω)}, H1/2(Γ):={u∈H1/2(Γ):supp u(?)Γ}.容易证明,在Ω中,下述问题△u+k2u=0于rg∪S上方, (4)u=0于Γg∪S上. (5)▽·(k-2▽u)+u=0于Γg∪S上方, (6)(?)nu=0于Γg∪S上. (7)的等价变分形式为:求u∈Hs1(Ω)={u∈H1(Ω):u=0在S上}使得其中φ是全场u在T上的法向导数,即φ=(?)nu.,<·.·)表示H-1/2(T)与H1/2(T)的对偶运算.为研究边界积分方程,本文引入单层位势算子VTE,超奇异积分算子DTE,双层层位势算子KTE及其共轭算子K*TE,其定义如下全场边界r上的积分方程为:利用引入的算子,上述方程可以写成其中f=uref,g=(?)nuref,I为恒等算子.于是,下面两个方程组成了有限元法与边界积分方法对称耦合的变分公式,可以用于求解腔体正散射问题.利用关于积分算子的结果得到引理2单层位势算子VTEE是由H-1/2(T)映入H~1/2(T)的紧算子;双层位势算子KTE以及其伴随算子K*TE分别由H1/2(T)映入H1/2(T)和由H-1/2(T)映入H-1/2(T),且都是紧算子,超奇异积分算子DTE为由H~1/2(T)映入H-1/2(T)的紧算子.进一步,单层位势算子VTE和超奇异积分算子DTE在H-1/2(T)和H1/2(r)内是强制的,即,存在紧算子VO和DO使得Re[<VTEφ,φ>+<V0φ,φ>]≥C‖φ‖H-1/2(Γ) (?)φ∈H-1/2(Γ). Re[<DTEu,u>+<D0u,u>)≥‖u‖H1/2(Γ) (?)u∈H1/2(Γ).对于任意u=[u,φ]∈VTE,定义VTE=HS1(Ω)×H-1/2(Γ).VTE中的模定义为‖u‖VTE2=‖u‖H12(Ω)+‖φ‖H-12(Γ)利用上述结果,对于正散射问题得到定理1变分问题(14)-(15)存在唯一解[u,φ]∈VTE.反散射问题是从在Γ上测量的全场u来确定腔体壁S.这个问题解的唯一性可用下述定理描述.定理2设[uj,φj]为在Ωj上的解,(?)Ωj=Γ∪Sj,j=1.2如果在Γ上u1=u2,则S1=S2.根据建立的变分形式,得到关于区域导数的下述结果.定理3设[u.φ]为(14)-(15)在Q上的解.n为S上的外法向.则区域导数可以表示为M’TE(S.p)=[u’]Γ.φ’].其中[u’.φ’]∈H1(Ω)×H-1/2(T)为下述边值问题的解△u’+k2u’=0于Ω内. (17)于Γ上. (18)于Γ上. (19)u’=-(p·n)φ于S上. (20)其中φ’=(?)nu’于Γ.在实际应用中,需要讨论实际构造的腔体壁的稳定性,因为它可提供数据在什么程度上可信的信息.定理4如果p∈C2(S,R)且h>0充分小,则d(Ωh,Ω)≤C‖uh-u‖H1/2(Γ). (21)其中C为与h无关的常数.对于TM极化情形的正问题和反问题,由于它与TE极化情形相似,文中阐述了正问题的某些平行结果,并证明了场关于腔体局部稳定性结果,并给出了区域导数.对于正散射问题,考虑与TE极化情形相同的几何问题,设入射平面波ui=exp(iax1-iβx2)从上方入射到完全电导体表面.散射场满足▽·(k-2▽us)+us=-▽·[(k-2-k0-2)▽uref]于Γg∪S上方,(?)nus=-(?)nuref于Γg∪S上另外,要求散射场满足辐射条件类同定理1有下述定理定理5变分问题(9)有一个唯一解[u,φ]∈VTM.而对于反散射问题,本文研究解的唯一性及局部稳定性,并给出区域导数.唯一性的证明与TE形式十分相似,给出简要证明.但是其局部稳定性的证明与区域导数的给出与TE情形有很大的不同,所以本文给出了详细的证明.关于TM情形的主要结果如下:定理6设[uj,φj]为(4.9)在Ωj中的解,(?)Ωj=Γ∪Sj,J=1.2如果在Γ上u1=u2,则有S1=S2.定理7设[u,φ]为(9)在Ω的解,n为S上的外法向.区域导数可表示为M’TM(S.p)=[U’|Γ,φ’],其中[u’,φ’]∈VTM为下述边值问题的弱解▽·(k-2▽u’)+u’=0于Ω内, (22)于Γ上, (23)于Γ上, (24)k-2φ’=▽s·[k-2(p·n)▽Su]+(p·n)u于S上, (25)其中,在Γ上φ’=(?)nu’.定理8若p∈C2(S,R)且h>0充分小,则有d(Ωh,Ω)≤C‖uh-u‖H1/2(Γ)·(26)其中C为一与h无关的常数.Ωh为由Sh和Γ所截的区域,Sh:x+hp(x)n,p∈C2(S,R).
其他文献
本文主要考虑了非标准增长条件下的几类非线性抛物方程解的性质.这类问题来源于具有丰富物理背景的电子流变流体学、非线性弹性力学以及图像处理等实际问题.目前人们主要研究了这类问题稳态解的性质及相关解空间的性质,而本文则讨论发展型问题解的性质.本文分为如下四章:在第一章中,主要介绍了关于非线性抛物方程研究的发展状况,进一步还介绍了变指数Sobolev空间的发展以及它与常指数Sobolev空间的一些差异,最
创立于上世纪五六十年代的KAM理论是迄今为止最伟大的数学成就之一.本文沿用这一经典的思想和证明技巧,研究当代应用数学领域的热点问题.全文共分三个部分,在第一章绪论中,我们介绍了经典KAM理论以及低维环面保持性定理的发展及其应用,第二,三章是本文的主体部分.在第二章中,我们概述了无穷维KAM-型定理的发展过程.在上世纪八九十年代,Kuksin等人首先将经典的KAM理论推广到无穷维哈密顿系统中,证明了
本文利用最优控制理论的思想,研究一类跨共振的二阶常微分方程边值问题.全文共分四章,第一章是绪论,分小节简述常微分方程边值问题、定性理论和最优控制理论的发展,并给出本文所研究问题的背景及国内外发展现状.在第五小节中给出了本文的工作及今后工作的展望.在第二章中给出控制和最优控制问题的一般概念及经典原理,并列出本文需要的一些预备知识.第三章和第四章为论文主体部分.本文致力于探求周期-积分边值问题在跨共振
本文主要研究一类带有真空的燃烧非牛顿流的存在唯一性.在这里燃烧非牛顿流体分成两种情况:第一种情况是燃烧剪切变稀流其中4/32,μ0>0.本文讨论了上述两种情况带有如下条件的初边值问题,克服了非线性,奇异性以及真空出现等实质性困难,得到了局部强解的存在唯一性.
本文研究具时滞非线性扩散方程的一些定性问题,它来源于自然界中许多时滞扩散现象,如传染病的蔓延、物理学中激光模型等.由于时滞非线性扩散模型具有广泛的实际背景及重要的理论价值,近年来,这类问题正吸引着越来越多的数学工作者的注意,并成为国内、外关注的热点问题之一.本文的目的在于研究模型的周期解问题、解的长时间渐近性态及行波解问题.本文第一章是绪论部分.在第二章里,我们首先讨论具多方渗流扩散的Nichol
素环上的导子和三角代数上的映射问题是有着深刻理论意义和丰富研究内容的研究课题.本文主要研究素环上的导子和广义导子以及三角代数上Jordan映射的加性和协交换映射.(1)证明了特征不为2的素环R若具有在某个非零理想上满足恒等式[d(χκ),χκ]∈Z的非零导子d,则R满足4个变量的标准恒等式s4.推广了Carini和Filippis的结果.(2)研究了在素环R的某个非中心Lie理想上满足恒等式χs[
美国怀俄明州始新世硅化木丰富且保存完好。本研究从两个产地(Big Sandy Reservoir的Bridger组和Parnell Draw的绿河组)分析了它们的植物多样性和系统分类关系。第一产地(Big SandyReservoir)多样性较低,包括一单子叶植物Palmoxylon macginitiei Tidwell, Simper et Medlyn (棕榈科);4科4属约4种双子叶植物E
长江刀鱼,学名刀鲚(Coilia ectenes),隶属于鲱形目、鯷科、鲚属,体形狭长侧薄,似剪刀,故名"刀鱼",为"长江三鲜"之首,以鲜、嫩、美而闻名,是长江珍稀名贵的生殖洄游性鱼类。因水域生态环境的恶化、过度捕捞和涉水工程建设等因素影响,
期刊
本文主要研究非线性波动方程的解的构造及相关控制问题.全文共分三章,第一章对最优控制和非线性发展方程行波解的精确求解相关方法做一些简单综述,第二章和第三章是论文的主要部分.第二章研究了一般的非线性波动方程的行波解的显式解的存在性,通过两种方法我们分别得到了同样的行波解,验证了两种方法的一致性.第三章研究了依赖于空间变量的一维变系数非线性波动方程的最优控制问题,所考虑的控制是通过状态的外部函数实现的.
Loop作为蛋白质二级结构的重要组成元件,其柔性和运动性对于蛋白质功能,尤其是酶催化作用具有重要作用,如定位重要的催化残基、参与底物结合以及调控底物的进入和产物的释放。三级结构分析揭示α/β水解酶超家族成员(脂肪酶第五家族、环氧化物水解酶、卤代烷脱卤酶和过水解酶等家族)存在一段loop区(NC-loop)负责连接α/β水解酶催化结构域的N-末端和盖子结构域(Cap),并且NC-loop在序列、长度