×
M.鲍曼。;Astudillo，R。;邱，Y。;Ang，E.Y.M。;van Gijzen，医学博士。;普莱西斯共和国。

多频率时谐弹性波方程的MSSS预处理矩阵方程方法。 （英语） Zbl 1405.65146号

计算。地质科学。 22，第1号，43-61（2018）.
小结：在这项工作中，我们提出了一个新的数值框架，用于有效求解多频率下的时谐弹性波动方程。我们表明，当离散化问题写成矩阵方程时，可以合并多个频率（和多个右侧）。使用预处理方法可以有效地求解该矩阵方程。我们提出了一种使用MSSS矩阵计算应用单个预条件的高效且稳健的方法。对于3D问题，我们提出了一种节省内存的实现，它利用了一系列2D问题的解决方案。二维和三维空间中的真实示例验证了新算法的性能。

引用于5文件

MSC公司：

65N21型 含偏微分方程边值问题反问题的数值方法
65N30型 偏微分方程边值问题的有限元、Rayleigh-Riz和Galerkin方法

关键词：

时谐弹性波动方程;多个频率;诱导降维法;预处理矩阵方程;多级顺序半可分矩阵

软件：

BiCG选项卡;MUMPS公司;斯帕斯基;算法913
PDF格式BibTeX公司 XML格式引用
\textit{M.Baumann}等人，计算。地质科学。22，第1号，43-61（2018；Zbl 1405.65146）
全文： DOI程序
OA许可证

参考文献：

[1] Airaksinen，T.，Pennanen，A.，Toivanen，J.：流体和弹性材料中时间谐波方程的阻尼预处理器。J.计算。物理学。228(5), 1466-1479 (2009) ·Zbl 1157.76046号 ·doi:10.1016/j.jcp.2008.10.036
[2] Amestoy，P.，Ashcraft，C.，Boiteau，O.，Buttari，A.，L'Excellent，J.Y.，Weisbecker，C.：通过块低秩表示改进多面方法。SIAM J.科学。计算。37，A1451-A1474（2015）·Zbl 1314.05111号 ·数字对象标识代码：10.1137/120903476
[3] Amestoy，P.、Brossier，R.、Buttari，A.、L'Excellent，J.Y.、Mary，T.、Métiver，L.、Miniussi，A.和Operto，S.：采用并行块低阶多面直接解算器的快速三维频域全波形反演：应用于北海的OBC数据。地球物理学81（6），R363-R383（2016）·doi:10.1190/geo2016-0052.1
[4] Astudillo，R.，van Gijzen，M.B.：求解线性矩阵方程的诱导降维方法。《Procedia Computer Science》80，222-232（2016）
[5] Baumann，M.：二维和三维时间谐波弹性波动方程的两个基准问题。数字对象标识：http://dx.doi.org/https://github.com/ManuelMBaumann/fractice_benchmarks网址（2016年9月）。doi:10.5281/zenodo.154700
[6] Baumann，M.，van Gijzen，M.B.：移位线性系统的嵌套Krylov方法。SIAM J.科学。计算。37（5），S90-S112（2015）·Zbl 06502306号
[7] Baumann，M.，van Gijzen，M.B.：用于多频波传播问题的有效两级预处理器。技术代表，DIAM报告17-03代尔夫特理工大学（2017）·Zbl 1406.65104号
[8] Chandrasekaran，S.，Dewilde，P.，Gu，M.，Pals，T.，Sun，X.，van der Veen，A.，White，D.：序列半可分表示的一些快速算法。SIAM J.矩阵分析。申请。27(2), 341-364 (2005) ·Zbl 1091.65063号
[9] Chandrasekaran，S.，Dewilde，P.，Gu，M.，Pals，T.，van der Veen，A.J.：序列半可解线性方程组的快速稳定解算器。劳伦斯·利弗莫尔国家实验室技术代表（2003年）·Zbl 1032.65508号
[10] Cottrell，J.A.、Hughes，T.J.R.、Bazilevs，Y.：等几何分析。走向CAD和FEA的集成。John Wiley&Son有限公司（2009年）·Zbl 1378.65009号
[11] De Basabe，J.：地震波传播的高阶有限元方法。博士论文德克萨斯大学奥斯汀分校（2009）
[12] Dewilde，P.，Van der Veen，A.：时间变量系统和计算。Kluwer学术出版社，波士顿（1998）·Zbl 0937.93002号
[13] Eidelman，Y.，Gohberg，I.：关于拟可分有限块矩阵的生成元。Calcolo 42（3），187-214（2005）·Zbl 1146.15300号 ·doi:10.1007/s10092-005-0102-4
[14] Elman，H.，Silvester，D.，Wathen，A.：有限元和快速迭代求解器：在不可压缩流体动力学中的应用。数值数学与科学计算。牛津大学出版社（2014）·Zbl 1304.76002号
[15] Etienne，V.，Chaljub，E.，Virieux，J.，Glinsky，N.：三维弹性波建模的hp自适应间断Galerkin有限元方法。地球物理学。《国际期刊》183（2），941-962（2010）·doi:10.1111/j.1365-246X.2010.04764.x
[16] van Gijzen，M.B.，Sonneveld，P.：算法913：高效利用双正交性的优雅IDR变体。ACM事务处理。数学。软件38（1），5:1-5:19（2011）·Zbl 1365.65089号
[17] Jbilou，K.，Messaoudi，A.，Sadok，H.：矩阵方程的全局FOM和GMRES算法。申请。数字。数学。31, 49-63 (1999) ·Zbl 0935.65024号 ·doi:10.1016/S0168-9274（98）00094-4
[18] Kavcic，A.，Moura，J.：带逆矩阵：高斯-马尔科夫过程的反演算法和因式分解。IEEE传输。通知。理论46（4），1495-1509（2000）·Zbl 0998.65037号 ·doi:10.1109/18.954748
[19] Knibbe，H.，Vuik，C.，Oosterlee，C.W.：通过图形处理单元减少基于亥姆霍兹方程的最小二乘偏移的计算时间。计算。地质科学。20(2), 297-315 (2016) ·Zbl 1390.86008号 ·doi:10.1007/s10596-015-9546-z
[20] Liesen，J.，Strakos，Z.：Krylov子空间方法：原理与分析。牛津大学数学与科学计算学院（2013）·Zbl 1263.65034号
[21] 马丁，GS；马福特，KJ；Larsen，S.，Marmousi-2：结构复杂地区AVO调查的更新模型，1979-1982（2002）
[22] Mulder，W.A.，Plessix，R.E.：如何在频域有限差分偏移中选择频率子集。地球物理学。《国际期刊》158，801-812（2004）·文件编号：10.1111/j.1365-246X.2004.02336.x
[23] Petrov，P.V.，Newman，G.A.：傅里叶域阻尼弹性波传播的三维逆模拟。地球物理学。《国际期刊》1981599-1617（2014）·doi:10.1093/gji/ggu222
[24] Plessix，R.E.：三维地震成像问题的亥姆霍兹迭代求解器。地球物理学72（5），SM185-SM194（2007）·doi:10.1190/1.2738849
[25] Plessix，R.E.：用迭代求解器进行三维频域全波形反演。地球物理学74，WCC149-WCC157（2009）·doi:10.1190/1.3211198
[26] Plessix，R.E.，Mulder，W.A.：作为迭代亥姆霍兹解算器的预处理器的变量分离。申请。数字。数学。44, 385-400 (2004) ·Zbl 1013.65117号 ·doi:10.1016/S0168-9274（02）00165-4
[27] Plessix，R.E.，Pérez Solano，C.A.：低频广角活动陆地地震数据弹性波形反演的修正表面边界条件。地球物理学。《国际期刊》第201卷，第1324-1334页（2015年）·doi:10.1093/gji/ggv087
[28] Pratt，R.：频域地震波形反演，第1部分：物理尺度模式下的理论和验证。地球物理学64（3），888-901（1999）·数字对象标识代码：10.1190/1.144597
[29] 邱，Y.，van Gijzen，M.B.，van Wingerden，J.W.，Verhaegen，M.，Vuik，C.：多层序贯半可分矩阵结构PDE约束优化问题的有效前置条件。电子。事务处理。数字。分析。44, 367-400 (2015) ·Zbl 1330.65168号
[30] 邱，Y.，van Gijzen，M.B.，van Wingerden，J.W.，Verhaegen，M.，Vuik，C.：使用不可压缩流和迭代求解器软件评估CFD基准问题上的多级顺序半可分预条件。数学。方法应用。科学。38 (2015) ·Zbl 1386.65104号
[31] Rice，J.：分布式控制的有效算法：结构化矩阵方法。代尔夫特理工大学博士论文（2010年）
[32] Rice，J.，Verhaegen，M.：分布式控制：空间异质线性系统的顺序半分离方法。IEEE传输。自动垫。控制54（6），1270-1283（2009）·Zbl 1367.93186号 ·doi:10.1109/TAC.2009.2019802
[33] Riyanti，C.D.，Erlangga，Y.A.，Plessix，R.E.，Mulder，W.A.，Vuik，C.，Osterlee，C.：时谐波动方程的新迭代求解器。地球物理学71，E57-E63（2006）·数字对象标识代码：10.1190/1.2231109
[34] Rizzuti，G.，Mulder，W.：基于多重网格的时谐弹性波方程的“位移-拉普拉斯”预处理。J.计算。物理学。317, 47-65 (2016) ·Zbl 1349.74364号 ·doi:10.1016/j.jcp.2016.04.049
[35] Saad，Y.：SPARSEKIT:稀疏矩阵计算的基本工具包。明尼阿波利斯明尼苏达大学技术代表（1994年）
[36] Saad，Y.：稀疏线性系统的迭代方法：第二版。工业和应用数学学会（2003）·Zbl 1031.65046号
[37] Saad，Y.，Schultz，M.:GMRES：求解非对称线性系统的广义最小残差算法。SIAM J.科学。统计计算。7, 856-869 (1986) ·Zbl 0599.65018号 ·doi:10.1137/0907058
[38] Saibaba，A.，Bakhos，T.，Kitanidis，P.：一种适用于位移系统的灵活Krylov解算器，并应用于振荡液压层析成像。SIAM J.科学。计算。35, 3001-3023 (2013) ·Zbl 1288.65041号 ·doi:10.137/120902690
[39] Sleijpen，G.L.G.，Sonneveld，P.，van Gijzen，M.B.：BiCGStab作为一种诱导降维方法。申请。数字。数学。60, 1100-1114 (2010) ·Zbl 1200.65024号
[40] Sleijpen，G.L.G.，van der Vorst，H.A.：在有限精度算法中保持BiCGstab方法的收敛特性。数字。算法10203-223（1995）·Zbl 0837.65030号
[41] Sonneveld，P.，van Gijzen，M.B.：IDR（S）：求解大型非对称线性系统的一系列简单快速算法。SIAM J.科学。计算。31(2), 1035-1062 (2008) ·Zbl 1190.65053号
[42] Tsuji，P.，Poulson，J.，Engquist，B.，Ying，L.：用谱元方法对弹性波传播进行扫描预处理。ESAIM数学。模型。数字。分析。48(2), 433-447 (2014) ·Zbl 1423.74950号 ·doi:10.1051/m2安/2013114
[43] Vandebril，R.，Van Barel，M.，Mastronardi，N.：矩阵计算和半可分矩阵：线性系统。约翰·霍普金斯大学出版社，巴尔的摩（2007）·Zbl 1141.65019号
[44] Virieux，J.，Operto，S.：勘探地球物理学中的全波形反演概述。地球物理学73（6），VE135-VE144（2009）
[45] van der Vorst，H.A.：BiCGStab:非对称线性系统解的bi-CG的一种快速平滑收敛变体。SIAM J.科学。统计计算。13(2), 631-644 (1992) ·Zbl 0761.65023号
[46] Wang，S.，de Hoop，M.V.，Xia，J.，Li，X.：三维各向异性介质中时谐弹性波的大规模并行结构多面解算器。地球物理学。《国际期刊》191（1），346-366（2012）
[47] Xia，J.：一般大型稀疏矩阵的高效结构化多前沿因式分解。SIAM J.科学。计算。35（2），A832-A860（2013）·Zbl 1266.15022号 ·数字对象标识代码：10.1137/120867032
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配，可能包含数据转换错误。在某些情况下，zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献，而不要求完整或完全匹配。