A stable high-order perturbation of surfaces method for numerical simulation of diffraction problems in triply layered media

JOURNAL OF COMPUTATIONAL PHYSICS

卷:330

Article

Hong, Youngjoon¹  Nicholls, David P.¹ 

[1] Univ Illinois, Dept Math Stat & Comp Sci, Chicago, IL 60607 USA

关键词: High-order spectral methods;    Time-harmonic linear wave scattering;    Periodic layered media;    High-order perturbation of surfaces methods;   

DOI  :  10.1016/j.jcp.2016.10.057

来源: Elsevier

【 摘 要 】

The accurate numerical simulation of linear waves interacting with periodic layered media is a crucial capability in engineering applications. In this contribution we study the stable and high-order accurate numerical simulation of the interaction of linear, time-harmonic waves with a periodic, triply layered medium with irregular interfaces. In contrast with volumetric approaches, High-Order Perturbation of Surfaces (HOPS) algorithms are inexpensive interfacial methods which rapidly and recursively estimate scattering returns by perturbation of the interface shape. In comparison with Boundary Integral/Element Methods, the stable HOPS algorithm we describe here does not require specialized quadrature rules, periodization strategies, or the solution of dense non symmetric positive definite linear systems. In addition, the algorithm is provably stable as opposed to other classical HOPS approaches. With numerical experiments we show the remarkable efficiency, fidelity, and accuracy one can achieve with an implementation of this algorithm. (C) 2016 Elsevier Inc. All rights reserved.

【 授权许可】

Free