×

无射线追踪的一般几何体稠密边界元矩阵的高频渐近压缩。 (英语) Zbl 1422.65430号

摘要:声学中的波传播和散射问题通常用边界元法解决。它们导致了一个典型的密集和大的离散化矩阵:它的大小和条件数随着频率的增加而增长。然而,高频散射问题本质上是局部的,这很好地表现为高度局域的射线在周围反弹。通过使用光线追踪或类似技术从解中明确提取振荡特性,可以使用渐近方法减小线性系统的尺寸,甚至使其与频率无关。然而,在存在具有复杂几何结构的(多)散射障碍物的情况下,光线跟踪变得昂贵甚至难以处理。在本文中,我们从构造完全分解的大稠密矩阵的相同离散化出发,通过显式定位格林函数来实现渐近压缩。这会产生一个大型但稀疏的矩阵,具有更快的相关矩阵-向量乘积,并且如数值实验所示,条件数大大提高。虽然格林函数的适当局部化还取决于一般几何图形所无法获得的渐近信息,但我们可以在从小频率到大频率的扫频中自适应地构造它,从而自动考虑到一般入射波。我们表明,该方法对于非凸、多个甚至近捕获域是鲁棒的,尽管在后一种情况下压缩率明显较低。此外,尽管该方法具有渐近性质,但对于应用中常用的低阶离散化,如分段常数、线性或立方,该方法是稳健的。另一方面,与传统的经典离散化相比,我们没有减少自由度的总数。格林函数的稀疏化修改与其他加速方案(如快速多极子方法)的结合在原则上似乎是可能的,并且是一个未来的研究课题。

理学硕士:

65号38 偏微分方程边值问题的边界元方法
65英尺50英寸 稀疏矩阵的计算方法
45A05型 线性积分方程
2005年4月5日 积分方程解的渐近性
65兰特 积分方程的数值解法
65号35 偏微分方程边值问题的谱、配置及相关方法
65N30型 含偏微分方程边值问题的有限元、Rayleigh-Ritz和Galerkin方法
65层35 矩阵范数、条件、缩放的数值计算
2005年第76季度 水力和气动声学
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] Alouges,F.,Aussal,M.:稀疏基数正弦分解及其在快速数值卷积中的应用。数字。算法70,427-448(2015)·Zbl 1326.65183号 ·doi:10.1007/s11075-014-9953-6
[2] Anand,A.,Boubendir,Y.,Ecevit,F.,Reitich,F.:高频散射问题的多次散射迭代分析。二: 三维标量情况。数字。数学。114, 373-427 (2010) ·Zbl 1196.65180号 ·doi:10.1007/s00211-009-0263-1
[3] Asheim,A.,Huybrechs,D.:在高频散射问题中跨平滑阴影边界提取一致精确的相位函数。SIAM J.应用。数学。74(2), 454-476 (2014) ·Zbl 1295.65118号 ·数字对象标识代码:10.1137/130936075
[4] Babich,V.M.,Buldyrev,V.S.:短波长衍射理论。柏林施普林格(1991)·Zbl 0742.35002号 ·数字对象标识代码:10.1007/978-3-642-83459-2
[5] Bebendorf,M.:边界元的基于层次LU分解的预条件器。计算74225-247(2005)·Zbl 1071.65031号 ·文件编号:10.1007/s00607-004-0099-6
[6] Bleistein,N.,Handelsman,R.A.:积分的渐近展开。多佛出版公司,Mineola(1986)·Zbl 0327.41027号
[7] Brakhage,H.,Werner,P.:Dirichletsche Außenraumproblem für die Helmholtzsche Schwingungsgleichung。架构(architecture)。数学。16, 325-329 (1965) ·Zbl 0132.33601号 ·doi:10.1007/BF01220037
[8] Bruno,O.,Geuzaine,C.,Monro,J.J.,Reitich,F.:任意高频散射问题在固定计算时间内的规定误差容限:凸情况。菲洛斯。事务处理。R.Soc.伦敦。A 362629-645(2004)·Zbl 1073.78004号 ·doi:10.1098/rsta.2003.1338
[9] Beylkin,G.,Kurcz,C.,Monzón,L.:亥姆霍兹-格林函数的快速算法。程序。R.Soc.序列。A 464,3301-3326(2008)·Zbl 1186.65155号 ·doi:10.1098/rspa.2008.0161
[10] Chandler-Wilde,S.、Graham,I.、Langdon,S.和Spence,E.:高频声散射中的数值渐近边界积分方法。数值学报。21, 89-305 (2012) ·Zbl 1257.65070号 ·doi:10.1017/S0962492912000037
[11] Chandler-Wilde,S.N.,Hewett,D.P.,Langdon,S.,Twigger,A.:一类非凸障碍物散射的高频边界元方法。数字。数学。129(4), 647-689 (2015) ·Zbl 1314.65154号 ·doi:10.1007/s00211-014-0648-7
[12] Cheng,H.,Crutchfield,W.Y.,Gimbutas,Z.,Greengard,L.F.,Ethridge,J.F.,Huang,J.,Rokhlin,V.,Yarvin,N.,Zhao,J.:三维亥姆霍兹方程的宽带快速多极方法。J.计算。物理学。216(1), 300-325 (2006) ·Zbl 1093.65117号 ·doi:10.1016/j.jp.2005.12.001
[13] Colton,D.L.,Kress,R.:散射理论中的积分方程方法。威利,纽约(1983年)·Zbl 0522.35001号
[14] Deaño,A.,Huybrechs,D.,Iserles,A.:《计算高振荡积分》,第155卷。费城SIAM(2018)·Zbl 1400.65004号
[15] Domínguez,V.,Graham,I.G.,Smyshlyaev,V.:高频声散射的混合数值渐近边界积分方法。数值。数学。106, 471-510 (2007) ·Zbl 1169.76056号 ·doi:10.1007/s00211-007-0071-4
[16] Ecevit,F.,Reitich,F.:高频散射问题的多次散射迭代分析。一: 二维情况。数字。数学。114, 271-354 (2009) ·Zbl 1196.65184号 ·doi:10.1007/s00211-009-0249-z
[17] Ganesh,M.,Hawkins,S.:三维高频外部声散射的全离散Galerkin方法。J.计算。物理学。230, 104-125 (2011) ·兹比尔1432.76158 ·doi:10.1016/j.jcp.2010.09.014
[18] Ganesh,M.,Langdon,S.,Sloan,I.:高度振荡声散射表面积分的有效评估。J.计算。申请。数学。204, 363-374 (2007) ·Zbl 1122.65025号 ·doi:10.1016/j.cam.2006.03.029
[19] Geuzaine,C.,Bruno,O.,Reitich,F.:关于多重散射问题的O(1)解。IEEE传输。Magn.公司。41(5), 1488-1491 (2005) ·doi:10.1109/TMAG.2005.844567
[20] Geuzaine,C.,Remacle,J.F.:Gmsh:一个带有内置预处理和后处理设施的三维有限元网格生成器。国际期刊编号。方法工程79(11),1309-1331(2009)·Zbl 1176.74181号 ·doi:10.1002/nme.2579
[21] Giladi,E.:高频亥姆霍兹方程的渐近导出边界元。J.计算。申请。数学。198, 52-74 (2007) ·Zbl 1102.65120号 ·doi:10.1016/j.cam.2005.11.024
[22] Greengard,L.,Rokhlin,V.:粒子模拟的快速算法。J.计算。物理学。73(2), 325-348 (1987) ·Zbl 0629.65005号 ·doi:10.1016/0021-9991(87)90140-9
[23] Groth,S.、Huybrechs,D.、Opsomer,P.:单障碍物和多障碍物高频散射射线扩展中的高阶项(2018年)(准备中)
[24] Groth,S.P.,Hewett,D.P.,Langdon,S.:可穿透凸多边形高频散射的混合数值渐近近似。IMA J.应用。数学。80, 324-353 (2015) ·Zbl 1328.35230号 ·doi:10.1093/imamat/hxt040
[25] Harrington,R.F.:时间谐波电磁场。IEEE出版社,Piscatawat(1961)
[26] Huybrechs,D.,Vandewalle,S.:具有高振荡核的积分方程矩阵压缩的二维小波包变换。J.计算。申请。数学。197, 218-232 (2006) ·兹比尔1101.65113 ·doi:10.1016/j.cam.2005.11.001
[27] Huybrechs,D.,Vandewalle,S.:高频散射问题积分方程公式的稀疏离散化。SIAM J.科学。计算。29(6), 2305-2328 (2007) ·Zbl 1154.65376号 ·doi:10.1137/060651525
[28] Huybrechs,D.,Vandewalle,S.:计算成本高昂的格林函数的边界元方法的有效实现。工程分析。已绑定。元素。32(8), 621-632 (2008) ·Zbl 1244.78023号 ·doi:10.1016/j.enganabound.2007.10.018
[29] Khoromskij,B.N.:Rd.J.Constr。约30599-620(2009年)·Zbl 1185.65051号 ·doi:10.1007/s00365-009-9068-9
[30] Khoromskij,B.N.,Veit,A.:使用qtt张量近似高效计算高振荡积分。计算。方法应用。数学。16(1), 145-159 (2016) ·Zbl 1329.65088号 ·doi:10.1515/cmam-2015-0033
[31] Klöckner,A.,Barnett,A.,Greengard,L.,O'Neil,M.:扩展求积:一种评估层势的新方法。J.计算。物理学。252, 332-349 (2013) ·Zbl 1349.65094号 ·doi:10.1016/j.jcp.2013.06.027
[32] Kress,R.,Spassov,W.T.:关于声散射和电磁散射中边界积分算子的条件数。数值。数学。42, 77-95 (1983) ·Zbl 0534.35027号 ·doi:10.1007/BF01400919
[33] Melrose,R.B.,Taylor,M.E.:凸障碍物的近峰值散射和修正的基尔霍夫近似。高级数学。55, 242-315 (1985) ·兹比尔0591.58034 ·doi:10.1016/0001-8708(85)90093-3
[34] Opsomer,P.:发布:渐进压缩版本3。https://github.com/popsomer/bempp.git (2016)
[35] Rokhlin,V.:经典势理论积分方程的快速求解。J.计算。物理学。60, 187-207 (1985) ·Zbl 0629.65122号 ·doi:10.1016/0021-9991(85)90002-6
[36] Śmigaj,W.,Betcke,T.,Arridge,S.,Phillips,J.,Schweiger,M.:用BEM++解决边界积分问题。ACM事务处理。数学。柔和。41(2), 6:1-6:40 (2015) ·Zbl 1371.65127号
[37] Sweldens,W.,Piessens,R.:光滑函数的小波近似的求积公式和渐近误差展开。SIAM J.数字。分析。31(4), 1240-1264 (1994) ·Zbl 0822.65013 ·doi:10.1137/0731065
[38] Wald,I.、Mark,W.R.、Gunther,J.、Boulos,S.、Thiago,I.,Hunt,W.、Parker,S.G.、Shirley,P.:光线追踪动画场景的艺术状态。纽波特Eurograph(2007)
[39] Wong,R.S.:积分的渐近逼近。SIAM,费城(2001)。(1944年重新出版)·Zbl 1078.41001号
[40] Wu,T.:边界元素声学。WIT出版社(2000)。(2005年重印)·Zbl 0987.76500号
[41] Ying,L.:通过局部FFT快速定向计算高频边界积分。SIAM多尺度模型。模拟。13(1), 423-439 (2015) ·Zbl 1317.65244号 ·doi:10.1137/140985123
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。