多体散射问题的高精度迭代算法研究
发布时间:2021-01-29 12:04
声波的多体散射问题在工程与医学上有广泛的应用。本文提出了当散射体在均匀或者局部非均匀介质中时求解多体散射问题的一种高效迭代算法。该算法的核心思想是利用人工边界分别将每个散射体包围,对于局部非均匀介质我们要求人工边界要将非均匀介质一同包围,然后利用波的叠加原理,将散射波分解为单纯向外传播的场的叠加。原始多体散射问题将被分解为一系列有限多个单体散射问题,其中每个单体散射问题与其他散射问题的相互作用仅发生在散射体边界或者人工边界上。因此每个问题相互独立,可以在每一步迭代时采用并行计算。可将迭代法和不同的单体散射问题求解器进行耦合,这样可以应对各种不同的情形,因此该算法具有高度的灵活性。对于单体散射问题,本文采用的是高精度的谱元算法,在人工边界上加上无反射人工边界条件(NRBC),并且在计算特殊函数的展开系数时采用了精确的半解析格式进行计算。特别地,当散射体为狭长形时,本文使用椭圆形人工边界,并提出一种半解析格式用以计算椭圆形DtN边界条件中的积分。这种半解析格式是基于在谱元离散时选取的适当的谱单元映射,用来计算无反射人工边界条件中的Mathieu函数展开系数,从而得到高精度的谱元算法。该格式...
【文章来源】:湖南师范大学湖南省 211工程院校
【文章页数】:157 页
【学位级别】:博士
【部分图文】:
图2:左:椭圓坐标系,右:椭圆人工边界与圆形人工边界对比示意图??类似宁圓的情形,我们首先介绍椭圖坐标系⑷办笛卡尔坐??
?多体散射问题的高精度迭代算法研究???其中U表示对函数v取共轭,在i:述内积意义下,,L2(叫为Hilbert??空间。??给定整数m?>?0,可以给出标准Sobolev空间iF1识)的定??义:??Hm(n)?:=?G?L2(Q),?dsu?e?L2(0)5i?|s|?=?1,?2y?.??.?ym}.??其中s为多重指标。其上的范数可定义为??\\u\\m,n=?(?[?|5s'u|2d^|?.??\kKm?)??给定n_?1维闭合流形r?C股'?我们定义如下的边界迹积分:??('?,?v)y?=?^?uvds.??2.2.2圆形人工边界的谱元离散格式??本节我们介绍用圆形人工边界截断的边值问题(2.30)的谱元??离散格式。??H?广?r.-1--in??■1編漏眷??(a)曲边网格单元?(b)参考单元£的LGL点(c)物理单元Jfe上的LGL点??图3:圆形人工边界曲边网格季元示意:图,在参考单元上的LGL点以及通过??Gordon-Hall变换映射到某个物理单元上的:LGL点示意图??16??
?多体散射问题的高精度迭代算法研究???/flt?丨::El?f—网??i僵遲暑??響?I::::::::!:::::::::::::;??'SW?M?s?I? ̄?d?)??(a)曲边网格单元?(b)参考单元f上的LGL貞?(c)物理单元护的’?LGL貞??图4:椭圆人工边輿?曲边网格单元示意图,在参考单元上的LGL点以及通过??Gordon-Hall变换映射到某个物理单元上的LGL点示意图.??给定网格剖分T,只有参考单元兌到物理单元Xe的映射??尺和圆的不一样,其他的设置都是一样的。故椭圓人工边界的谱??元离散可表示为:寻找#?e%(-wine),使得??A{:u^w)?=?0,?VwGFp(O),?(2.42)??其中??乂(m,,?w)?=?-(Vm,,?十?k2(<,十〈办#,?w)〉rf,K,?(2.43)??我们同样使用截断的椭歸DtN映射,即??補??n=o?Mcn、R'q)?(2?44)??^Td^Ms^\R-q)?.??+?}?Ml?^?sen(?7;g),??其中是散射场在椭圖人工边界上的Mathiei!展开系??数。??从下一节开始将介绍本章的重点s即如何利用谱元网格推导??出一种“半解析格式”,从而实现对全局DtN算子中的积分(展开??19??
【参考文献】:
期刊论文
[1]An Efficient and Accurate Spectral Method for Acoustic Scattering in Elliptic Domains[J]. Qirong Fang~1,Jie Shen~(1,*) and Li-Lian Wang~2 1 Department of Mathematics,Purdue University,West Lafayette,IN 47907,USA. 2 Division of Mathematical Sciences,School of Physical and Mathematical Sciences,Nanyang Technological University,637616,Singapore.. Numerical Mathematics:Theory,Methods and Applications. 2009(03)
本文编号:3006855
【文章来源】:湖南师范大学湖南省 211工程院校
【文章页数】:157 页
【学位级别】:博士
【部分图文】:
图2:左:椭圓坐标系,右:椭圆人工边界与圆形人工边界对比示意图??类似宁圓的情形,我们首先介绍椭圖坐标系⑷办笛卡尔坐??
?多体散射问题的高精度迭代算法研究???其中U表示对函数v取共轭,在i:述内积意义下,,L2(叫为Hilbert??空间。??给定整数m?>?0,可以给出标准Sobolev空间iF1识)的定??义:??Hm(n)?:=?G?L2(Q),?dsu?e?L2(0)5i?|s|?=?1,?2y?.??.?ym}.??其中s为多重指标。其上的范数可定义为??\\u\\m,n=?(?[?|5s'u|2d^|?.??\kKm?)??给定n_?1维闭合流形r?C股'?我们定义如下的边界迹积分:??('?,?v)y?=?^?uvds.??2.2.2圆形人工边界的谱元离散格式??本节我们介绍用圆形人工边界截断的边值问题(2.30)的谱元??离散格式。??H?广?r.-1--in??■1編漏眷??(a)曲边网格单元?(b)参考单元£的LGL点(c)物理单元Jfe上的LGL点??图3:圆形人工边界曲边网格季元示意:图,在参考单元上的LGL点以及通过??Gordon-Hall变换映射到某个物理单元上的:LGL点示意图??16??
?多体散射问题的高精度迭代算法研究???/flt?丨::El?f—网??i僵遲暑??響?I::::::::!:::::::::::::;??'SW?M?s?I? ̄?d?)??(a)曲边网格单元?(b)参考单元f上的LGL貞?(c)物理单元护的’?LGL貞??图4:椭圆人工边輿?曲边网格单元示意图,在参考单元上的LGL点以及通过??Gordon-Hall变换映射到某个物理单元上的LGL点示意图.??给定网格剖分T,只有参考单元兌到物理单元Xe的映射??尺和圆的不一样,其他的设置都是一样的。故椭圓人工边界的谱??元离散可表示为:寻找#?e%(-wine),使得??A{:u^w)?=?0,?VwGFp(O),?(2.42)??其中??乂(m,,?w)?=?-(Vm,,?十?k2(<,十〈办#,?w)〉rf,K,?(2.43)??我们同样使用截断的椭歸DtN映射,即??補??n=o?Mcn、R'q)?(2?44)??^Td^Ms^\R-q)?.??+?}?Ml?^?sen(?7;g),??其中是散射场在椭圖人工边界上的Mathiei!展开系??数。??从下一节开始将介绍本章的重点s即如何利用谱元网格推导??出一种“半解析格式”,从而实现对全局DtN算子中的积分(展开??19??
【参考文献】:
期刊论文
[1]An Efficient and Accurate Spectral Method for Acoustic Scattering in Elliptic Domains[J]. Qirong Fang~1,Jie Shen~(1,*) and Li-Lian Wang~2 1 Department of Mathematics,Purdue University,West Lafayette,IN 47907,USA. 2 Division of Mathematical Sciences,School of Physical and Mathematical Sciences,Nanyang Technological University,637616,Singapore.. Numerical Mathematics:Theory,Methods and Applications. 2009(03)
本文编号:3006855
本文链接:https://www.wllwen.com/shoufeilunwen/jckxbs/3006855.html
教材专著