当前位置: X-MOL 学术Adv. Comput. Math. › 论文详情
Our official English website, www.x-mol.net, welcomes your feedback! (Note: you will need to create a separate account there.)
Analysis of a Helmholtz preconditioning problem motivated by uncertainty quantification
Advances in Computational Mathematics ( IF 1.7 ) Pub Date : 2021-09-03 , DOI: 10.1007/s10444-021-09889-0
I. G. Graham 1 , O. R. Pembery 1 , E. A. Spence 1
Affiliation  

This paper analyses the following question: let Aj, j = 1,2, be the Galerkin matrices corresponding to finite-element discretisations of the exterior Dirichlet problem for the heterogeneous Helmholtz equations ∇⋅ (Ajuj) + k2njuj = −f. How small must \(\|A_{1} -A_{2}\|_{L^{q}}\) and \(\|{n_{1}} - {n_{2}}\|_{L^{q}}\) be (in terms of k-dependence) for GMRES applied to either \((\mathbf {A}_1)^{-1}\mathbf {A}_2\) or A2(A1)− 1 to converge in a k-independent number of iterations for arbitrarily large k? (In other words, for A1 to be a good left or right preconditioner for A2?) We prove results answering this question, give theoretical evidence for their sharpness, and give numerical experiments supporting the estimates. Our motivation for tackling this question comes from calculating quantities of interest for the Helmholtz equation with random coefficients A and n. Such a calculation may require the solution of many deterministic Helmholtz problems, each with different A and n, and the answer to the question above dictates to what extent a previously calculated inverse of one of the Galerkin matrices can be used as a preconditioner for other Galerkin matrices.



中文翻译:

由不确定性量化驱动的亥姆霍兹预处理问题的分析

本文分析以下问题:设A j , j = 1,2, 为异构亥姆霍兹方程∇⋅ ( A ju j ) + k 2 n的外狄利克雷问题的有限元离散对应的伽辽金矩阵j u j = − f。必须有多小\(\|A_{1} -A_{2}\|_{L^{q}}\)\(\|{n_{1}} - {n_{2}}\|_{ L ^ {q}} \)是(在以下方面ķ为GMRES -dependence)施加到或者\((\ mathbf {A} _1)^ { - 1} \ mathbf {A} _2 \)21 ) − 1收敛于任意大k 的k无关的迭代次数?(换句话说,对于A 1来说,A 1A 2 的一个很好的左或右预处理器?)我们证明了回答这个问题的结果,为它们的锐度提供了理论证据,并给出了支持估计的数值实验。我们解决这个问题的动机来自于计算具有随机系数An的亥姆霍兹方程的关注量。这样的计算可能需要解决许多确定性亥姆霍兹问题,每个问题都有不同的An,并且上述问题的答案决定了先前计算的 Galerkin 矩阵之一的逆可以在多大程度上用作其他 Galerkin 矩阵的预处理器。

更新日期:2021-09-04
down
wechat
bug