×

使用水平集方法和块Krylov解算器反演涡流信号。 (英语) Zbl 1520.35149号

小结:推动这项工作的应用与使用涡流探头识别核电厂蒸汽发生器内部沉积物有关。我们考虑了一个现实的实验过程,该过程依赖于沿着管轴扫掠一个由线圈制成的探针来扫描一个区域,从而扮演源/接收器的角色。使用最小二乘法解决与这些测量相关的反形状问题需要在每个梯度衰减迭代中解决大量右侧的涡流和伴随问题。前向求解器的额外成本来自使用问题的潜在公式,该公式的优点是独立于导电介质的拓扑结构(在迭代过程中可能会发生变化)。我们使用级别集方法来避免重新网格化和处理未知拓扑。我们算法中的关键要素是一种处理大规模问题迭代求解器大量右侧的优化方法。我们首先对各种分块Krylov方法、分块GMRES和分块BGCRODR进行基准测试,以测试它们与标准对应方法(即GMRES和GCRODR)相比的有效性。然后,我们为BGCRODR提出了一种新的实现方法,用于从以前生成的Krylov基中回收信息,这种方法比传统方法具有更好的扩展性。这部分独立于手头的实际逆问题。最后,通过真实的三维合成实例验证了整体反演过程的有效性。

MSC公司:

35Q61问题 麦克斯韦方程组
78A45型 衍射、散射
78A46型 光学和电磁理论中的逆问题(包括逆散射)
2010年第49季度 优化最小曲面以外的形状
65N21型 含偏微分方程边值问题反问题的数值方法
65号55 多重网格方法;含偏微分方程边值问题的区域分解
65K10像素 数值优化和变分技术
65层10 线性系统的迭代数值方法
65F08个 迭代方法的前置条件
35兰特 PDE的反问题
PDF格式BibTeX公司 XML格式引用
全文: 内政部

参考文献:

[1] Allaire,G.、De Gournay,F.、Jouve,F.和Toader,A.-M.,通过水平集方法使用拓扑和形状敏感性进行结构优化,《控制网络》。,34(2005),第59-80页·Zbl 1167.49324号
[2] Amestoy,P.R.、Duff,I.S.、L'Excellent,J.-Y.和Koster,J.,《使用分布式动态调度的完全异步多前沿求解器》,SIAM J.Matrix Anal。申请。,23(2001),第15-41页,doi:10.1137/S0895479899358194·Zbl 0992.65018号
[3] Balay,S.、Abhyankar,S.,Adams,M.F.、Brown,J.、Brune,P.、Buschelman,K.、Dalcin,L.、Dener,A.、Eijkhout,V.、Gropp,W.D.、Karpeyev,D.、Kaushik,D.、Knepley,M.G.、May,D.A.、McInnes,L.C.、Mills,R.T.、Munson,T.、Rupp,K.,Sanan,P.,Smith,B.F.、Zampini,S.和Zhang,H.,PETSc用户手册,3.13版,阿贡国家实验室,2020年。
[4] Balay,S.、Abhyankar,S.,Adams,M.F.、Brown,J.、Brune,P.、Buschelman,K.、Dalcin,L.、Dener,A.、Eijkhout,V.、Gropp,W.D.、Karpeyev,D.、Kaushik,D.、Knepley,M.G.、May,D.A.、McInnes,L.C.、Mills,R.T.、Munson,T.、Rupp,K.,Sanan,P.,Smith,B.F.、Zampini,S.和Zhang,H.,PETSc,网址:http://www.mcs.anl.gov/petsc, 2020.
[5] Bavier,E.、Hoemmen,M.、Rajamanickam,S.和Thornquist,H.、Amesos2和Belos:大型稀疏线性系统的直接和迭代求解器,Sci。程序。,20(2012年),第241-255页。
[6] Bui,C.T.T.、Dapogny,C.和Frey,P.,《求解水平集平流方程的精确各向异性自适应方法》,《国际数值杂志》。《液体方法》,70(2012),第899-922页·Zbl 1412.65133号
[7] Cai,X.C.和Sarkis,M.,一般稀疏线性系统的限制加性Schwarz预条件,SIAM J.Sci。计算。,21(1999),第792-797页,doi:10.1137/S106482759732678X·Zbl 0944.65031号
[8] Calandra,H.、Gratton,S.、Langou,J.、Pinel,X.和Vasseur,X.,块重启GMRES方法的灵活变体及其在地球物理学中的应用,SIAM J.Sci。计算。,34(2012),第A714-A736页,doi:10.1137/10082364X·Zbl 1248.65036号
[9] Costabel,M.和Dauge,M.,《多面体域中电磁场的奇点》,Arch。定额。机械。分析。,151(2000),第221-276页·Zbl 0968.35113号
[10] Dapogny,C.、Frey,P.、Omnès,F.和Privat,Y.,使用freeFem++的流体力学几何形状优化,结构。多磁盘。最佳。,58(2018),第2761-2788页。
[11] Dolean,V.、Jolive,P.和Nataf,F.,《区域分解方法简介:算法、理论和并行实现》,SIAM,2015年,doi:10.1137/1.9781611974065·Zbl 1364.65277号
[12] Dorn,O.和Lesselier,D.,逆散射的水平集方法,逆问题,22(2006),第R67-R131页·Zbl 1191.35272号
[13] Dupuy,T.,Modélisation des transferrets thermiques dans les dépóts d'encrassement des générateurs de vapeur,博士论文,马赛中央大学,2019年。
[14] García-Martín,J.、Gómez-Gil,J.和Vázquez-Sánchez,E.,《基于涡流检测的非破坏性技术》,《传感器》,11(2011),第2525-2565页。
[15] Gutknecht,M.H.,《具有多个右手边的线性系统的Block Krylov空间方法:现实世界系统的现代数学模型、方法和算法简介》,Siddiqui,A.、Duff,I.和Christensen,O.编辑,Anamaya Publishers,2006年,第420-447页。
[16] Gutknecht,M.H.和Schmelzer,T.,更新块三对角矩阵和块Hessenberg矩阵的(QR)分解,应用。数字。数学。,58(2008),第871-883页·Zbl 1143.65032号
[17] Haddar,H.,Jiang,Z.和Riahi,M.K.,同轴涡流测量定量三维形状重建的稳健反演方法,科学杂志。计算。,70(2017年),第29-59页·Zbl 1440.94004号
[18] Hecht,F.,《FreeFem++的新发展》,J.Numer。数学。,20(2012),第251-266页·Zbl 1266.68090号
[19] Hernandez,V.、Roman,J.E.和Vidal,V.,SLEPc:用于解决特征值问题的可扩展且灵活的工具包,ACM Trans。数学软件。,31(2005),第351-362页,https://slepc.upv.es。 ·兹比尔1136.65315
[20] Heroux,M.A.、Bartlett,R.A.、Howle,V.E.、Hoekstra,R.J.、Hu,J.J.,Kolda,T.G.、Lehoucq,R.B.、Long,K.R.、Pawlowski,R.P.、Phipps,E.T.等人,《Trilinos项目概述》,ACM Trans。数学。柔和。(TOMS),31(2005),第397-423页,https://trilinos.github.io。 ·Zbl 1136.65354号
[21] Jiang,Z.、Haddar,H.、Lechleiter,A.和El-Guedri,M.,通过形状优化识别二维轴对称涡流模型中的磁性沉积物,逆问题。科学。《工程师》,24(2016),第1385-1410页·Zbl 1348.74112号
[22] Jolivet,P.,Roman,J.E.,and Zampini,S.,KSPHPDDM and PCHPDDM:使用鲁棒重叠Schwarz预条件和高级Krylov方法扩展PETSc,计算。数学。申请。,84(2021),第277-295页,https://github.com/prj-/jolivet2020petsc。 ·Zbl 1524.65135号
[23] Jolivet,P.和Tournier,P.-H.,Block迭代方法和循环以改进线性解算器的可伸缩性,《2016年高性能计算、网络、存储和分析国际会议论文集》,SC16,IEEE,2016年。
[24] Maurice,L.、Costan,V.、Guillot,E.和Thomas,P.、Eddy使用模拟工具进行NDE性能演示,AIP Conf.Proc。,1511(2013),第464-471页。
[25] Osher,S.和Fedkiw,R.,水平集方法和动态隐式曲面,应用。机械。第57版(2004年),第B15-B15页·Zbl 1026.76001号
[26] Parks,M.L.、de Sturler,E.、Mackey,G.、Johnson,D.D.和Maiti,S.,线性系统序列的Krylov子空间循环,SIAM J.Sci。计算。,28(2006),第1651-1674页,doi:10.1137/040607277·Zbl 1123.65022号
[27] Parks,M.L.、Soodhalter,K.M.和Szyld,D.B.,《块回收GMRES方法与解算器性能方面的研究》,arXiv:1604.017132016。
[28] Pierre,T.和Benjamin,G.,《CND-CF综合体的模拟》,载于《2017年COFREND期刊论文集》,《福柯现代化》,COFREND2017年。
[29] Poulson,J.、Marker,B.、Van de Geijn,R.A.、Hammond,J.R.和Romero,N.A.,《元素:分布式内存密集矩阵计算的新框架》,ACM Trans。数学。软质。,39 (2013), 13. ·Zbl 1295.65137号
[30] Prusek,T.,《模型与模拟》(Modélisation et simulation numérique du colmatageál’échelle du sous-canal dans les geénérateurs de vapeur),博士论文,Aix-Marseille,2012年。
[31] Riahi,M.K.,蒸汽发生器中的快速涡流无损检测有限元解算器,J.耦合系统。多尺度动态。,4(2016),第60-68页。
[32] Rodríguez,A.A.和Valli,A.,《麦克斯韦方程的涡流近似:理论、算法和应用》,第4卷,施普林格科学与商业媒体,2010年·Zbl 1204.78001号
[33] Roux,F.-X.和Barka,A.,应用于三维电磁波散射和辐射的FETI-2LM的Block-Krylov循环算法,IEEE Trans。《天线传播》,65(2017),第1886-1895页·Zbl 1395.78042号
[34] Saad,Y.和Schultz,M.H.,GMRES:求解非对称线性系统的广义最小残差算法,SIAM J.Sci。统计计算。,7(1986),第856-869页,doi:10.1137/0907058·Zbl 0599.65018号
[35] Schweiger,M.、Arridge,S.R.、Dorn,O.、Zacharopoulos,A.和Kolehmainen,V.,《利用水平集技术重建光学层析成像中的吸收和扩散形状轮廓》,Opt。莱特。,31(2006),第471-473页。
[36] Smith,B.F.、Björstad,P.和Gropp,W.D.,《区域分解:椭圆偏微分方程的并行多层方法》,剑桥大学出版社,2004年·兹比尔0857.65126
[37] Soodhalter,K.M.、de Sturler,E.和Kilmer,M.E.,《子空间循环迭代方法综述》,GAMM-Mitt。,43(2020年),e20200016。
[38] Stathopoulos,A.和Wu,K.,具有恒定同步要求的块正交化程序,SIAM J.Sci。计算。,23(2002),第2165-2182页,doi:10.1137/S1064827500370883·Zbl 1018.65050号
[39] Stewart,G.W.,大型特征问题的Krylov-Schur算法,SIAM J.矩阵分析。申请。,23(2001),第601-614页,doi:10.1137/S0895479800371529·Zbl 1003.65045号
[40] Toselli,A.和Widlund,O.B.,《区域分解方法:算法和理论》,Springer,2005年·Zbl 1069.65138号
[41] Tournier,P.-H.,Aliferis,I.,Bonazzoli,M.,de Buhan。,85(2019年),第88-97页。
[42] Wang,M.Y.,Wang,X.和Guo,D.,结构拓扑优化的水平集方法,计算。方法应用。机械。工程,192(2003),第227-246页·Zbl 1083.74573号
[43] Zenzinger,G.、Bamberg,J.、Satzger,W.和Carl,V.,涡流激励热像裂纹检测,无损检测。测试。评估。,22(2007),第101-111页。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。