2024年05月18日 星期六 登录 EN

学术活动
Direct/iterative hybrid solver for scattering by inhomogeneous media
首页 - 学术活动
报告人:
Ambuj Pandey, Assistant Professor, Indian Institute of Science Education and Research, Bhopal
邀请人:
Tao Yin, Associate Professor
题目:
Direct/iterative hybrid solver for scattering by inhomogeneous media
时间地点:
13:30-14:30 September 15 (Friday), Zoom meeting: 910 9946 0832, Passpord: 882471
摘要:

In this talk, we present a fast high-order method for the solution of two-dimensional problems of scattering by penetrable inhomogeneous media, with application to high-frequency configurations containing (possibly) discontinuous refractivities. The method relies on a hybrid direct/iterative combination of 1) A differential volumetric formulation (which is based on the use of appropriate Chebyshev differentiation matrices enacting the Laplace operator) and, 2) A second-kind boundary integral formulation (which, once again, utilizes Chebyshev discretization, but, in this case, in the boundary-integral context). The approach enjoys low dispersion and high-order accuracy for smooth refractivities, as well as second-order accuracy (while maintaining low dispersion) in the discontinuous refractivity case. The solution approach proceeds by application of Impedance-to-Impedance (ItI) maps to couple the volumetric and boundary discretizations. The volumetric linear algebra solutions are obtained by means of a multifrontal solver, and the coupling with the boundary integral formulation is achieved via an application of the iterative linear-algebra solver GMRES. In particular, the existence and uniqueness theory presented in the present paper provides an affirmative answer to an open question concerning the existence of a uniquely solvable second-kind ItI-based formulation for the overall scattering problem under consideration. Relying on a modestly-demanding scatterer-dependent precomputation stage (requiring in practice a computing cost of the order of O(Nα) operations, with α≈1.07, for an N-point discretization and for the relevant Chebyshev accuracy orders q used), together with fast (O(N)-cost) single-core runs for each incident field considered, the proposed algorithm can effectively solve scattering problems for large and complex objects possibly containing discontinuities and strong refractivity contrasts.