×

准周期散射问题的快速直接求解器。 (英语) Zbl 1349.78044号

摘要:我们考虑了二维均匀背景介质中反射或透射障碍物的无限周期阵列的时谐平面波散射的数值解。边界积分公式是理想的,因为它们将问题减少到障碍物边界上的未知值。然而,对于复杂的几何结构和/或更高的频率,由此产生的密集线性系统变得更大,排除了密集直接方法,并且通常情况下(尽管是第二类)条件不佳,使得快速的基于多极子的迭代方案也效率低下。我们提出了一种基于积分方程的复杂度为O(N)的求解器,该求解器利用“快速”直接线性代数的最新进展来分层反演孤立障碍矩阵,从而处理此类病态条件。这与一种新的周期化方案相结合,该方案对所有入射角(包括Wood’s异常)都具有鲁棒性,基于自由空间格林函数核。当需要多个入射角时,结果解算器非常有效,这在许多应用中都会发生。我们的数值测试包括一个多波长的复杂障碍物,其大小为(N=10^5),解的误差为(10^{-10}),其中每个入射角的解算器比快速多极迭代解快66倍,当选择入射角共享布洛赫相位时,解算器的速度快600倍。

MSC公司:

78A45型 衍射、散射
78M15型 边界元法在光学和电磁理论问题中的应用
65号38 偏微分方程边值问题的边界元方法

软件:

hlib公司
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Alpert,B.K.,混合高斯-三角形求积规则,SIAM J.Sci。计算。,20, 1551-1584 (1999) ·Zbl 0933.41019号
[3] 阿伦斯,T。;Chandler-Wilde,S.N。;DeSanto,J.A.,《衍射光栅散射的积分方程和最小二乘法》,Commun。计算。物理。,1, 1010-1042 (2006) ·Zbl 1137.78346号
[4] 阿特沃特,H.A。;Polman,A.,Plasmonics for improved voltage devices,Nat.Mater.,《改进光伏器件的等离子体学》。,9, 3, 205-213 (2010)
[5] Barnett,A.H。;Greengard,L.,二维准周期散射问题的新积分表示,BIT-Numer。数学。,51, 67-90 (2011) ·Zbl 1214.65061号
[6] Barty,C.P.J.,《激光聚变实验先进射线照相用LLNL高能短脉冲技术概述》,Nucl。融合,44,12,S266(2004)
[7] 拜尔金,G。;科伊夫曼,R。;Rokhlin,V.,《数值分析中的小波》(Wavelet and Their Applications,1992),Jones and Bartlett:Jones and Battlett Boston,MA),181-210·Zbl 0798.65126号
[8] Bonnet-BenDhia,A.-S。;Starling,F.,电磁光栅导波和衍射问题的非唯一性例子,数学。方法应用。科学。,17, 305-338 (1994) ·Zbl 0817.35109号
[9] Börm,S.,《非局部算子的有效数值方法:(H^2)-矩阵压缩、算法和分析》(2010),欧洲数学学会·Zbl 1208.65037号
[10] 布鲁诺,O.P。;Haslam,M.C.,《周期表面散射的高效高阶评估:深光栅、高频和掠入射》,J.Opt。《美国社会学杂志》,第26、3、658-668页(2009年)
[11] Chandrasekaran,S。;Gu,M.,对称块对角加半可分矩阵特征分解的分治算法,Numer。数学。,96, 4, 723-731 (2004) ·Zbl 1047.65023号
[12] 程,H。;克拉奇菲尔德,W.Y。;Gimbutas,Z。;Greengard,L。;黄,J。;Rokhlin,V。;北雅文。;赵,J.,关于实现二维亥姆霍兹方程宽带FMM的评论,(当代数学,第408卷(2006),美国数学学会:美国数学学会普罗维登斯,RI),99-110·Zbl 1107.65107号
[13] 程,H。;Gimbutas,Z。;Martinsson,P。;Rokhlin,V.,《关于低秩矩阵的压缩》,SIAM J.Sci。计算。,26, 4, 1389-1404 (2005) ·Zbl 1083.65042号
[14] Chew,W.C.,《非均匀介质中的波和场》(1999),威利-艾艾出版社
[15] 科尔顿,D。;Kress,R.,散射理论中的积分方程方法(1983),威利·Zbl 0522.35001号
[16] 科尔顿,D。;Kress,R.,《逆声和电磁散射理论》,《应用数学科学》,第93卷(1998年),Springer-Verlag:Springer-Verlag Berlin·Zbl 0893.35138号
[17] Gillman,A。;杨,P。;Martinsson,P.,一维区域积分方程的一个具有O(N)复杂性的直接求解器,Front。数学。中国,7,2,217-247(2012)·Zbl 1262.65198号
[19] Golub,G。;Loan,C.V.,《矩阵计算》,《约翰·霍普金斯数学科学研究》(1996),约翰·霍普金斯大学出版社:约翰·霍普金斯大学出版社,马里兰州巴尔的摩·Zbl 0865.65009号
[20] 顾,M。;Eisenstat,S.C.,《计算强秩揭示QR因式分解的高效算法》,SIAM J.Sci。计算。,17, 4, 848-869 (1996) ·Zbl 0858.65044号
[21] Hackbusch,W.,基于H-矩阵的稀疏矩阵算法;第一部分:H-矩阵导论,计算,6289-108(1999)·Zbl 0927.65063号
[23] Helsing,J.,混合型边界条件椭圆问题的积分方程方法,J.Compute。物理。,228, 8892-8907 (2009) ·Zbl 1177.65176号
[24] Ho,K。;Greengard,L.,《通过递归骨架化快速直接求解结构化线性系统》,SIAM J.Sci。计算34、5、A2507-A2532(2012)·Zbl 1259.65062号
[25] 霍罗申科夫,K.V。;Chandler-Wilde,S.N.,《二维周期和波导声学格林函数的有效计算》,J.Acoust。《美国社会》,1111610-1622(2002)
[26] 卡普尔,S。;Rokhlin,V.,奇异函数的高阶修正梯形求积规则,SIAM J.Numer。分析。,34, 1331-1356 (1997) ·兹伯利0891.65019
[28] Kress,R.,《数值分析》,《数学文本》,第181卷(1998年),斯普林格-Verlag·Zbl 0913.65001号
[29] Linton,C.M.,《亥姆霍兹方程的晶格和》,SIAM Rev.,52,4,603-674(2010)·Zbl 1208.78016号
[30] Martinsson,P。;Rokhlin,V.,《二维边界积分方程的快速直接求解器》,J.Compute。物理。,205, 1, 1-23 (2005) ·Zbl 1078.65112号
[31] Meier,A。;阿伦斯,T。;钱德勒·王尔德,S.N。;Kirsch,A.,一类实线积分方程的Nyström方法及其在衍射光栅和粗糙表面散射中的应用,J.integral Equ。申请。,12, 281-321 (2000) ·Zbl 1170.78365号
[32] Michielssen,E。;博格,A。;Chew,W.C.,《细长物体的散射:(O(N\log^2 N)运算的直接解》,IEE Proc。微型。天线传播。,143, 4, 277-283 (1996)
[33] 型号,R。;Rathsfeld,A。;总高度。;Wurm,M。;Bodermann,B.,《光学掩模计量中的散射测量反问题》,J.Phys.:Conf.序列号。,135, 012071 (2008)
[34] 莫尔斯,P。;Feshbach,H.,《理论物理方法》,第1卷(1953年),McGraw-Hill·Zbl 0051.40603号
[35] Nicholas,M.J.,三维双周期电磁散射问题的高阶数值方法,Commun。数学。科学。,6, 3, 669-694 (2008) ·Zbl 1168.78003号
[36] Peter,文学硕士。;Meylan,M.H。;Linton,C.M.,任意物体周期阵列的水波散射,J.流体力学。,548, 237-256 (2006)
[37] Rokhlin,V.,用第二类积分方程解决声散射问题,波动,5257-272(1983)·Zbl 0522.73022号
[38] 警长,N.P。;阿格拉瓦尔,M。;Peumans,P.,《使用涂层亚波长光栅的高性能太阳能选择性吸收器》,Opt。快递,18,6,5525-5540(2010)
[39] 盛,Z。;Dewilde,P。;Chandrasekaran,S.,《求解分层半可分系统的算法》,(系统理论,舒尔算法和多维分析,系统理论,舒尔算法和多元分析,算子理论:进步和应用,第176卷(2007年),Birkhä用户:Birkhá用户巴塞尔),255-294·Zbl 1123.65020号
[40] Shipman,S.,《开放周期波导的共振散射》,计算物理进展,第1卷(2010年),边沁科学出版社,第7-50页
[41] 斯塔尔,P。;Rokhlin,V.,关于两点边值问题的数值解:II,Commun。纯应用程序。数学。,47, 8, 1117-1159 (1994) ·兹比尔0805.65080
[42] Wu,Y。;Lu,Y.Y.,用边界积分方程Neumann-to-Dirichlet映射法分析衍射光栅,J.Opt。《美国社会学杂志》,第26、11、2444-2451页(2009年)
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。