Exponential Krylov time integration for modeling multi-frequency optical response with monochromatic sources

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:340

Article

Botchev, M. A.^1,2  Hanse, A. M.^3,4,5  Uppu, R.^5,6 

[1] Russian Acad Sci, Keldysh Inst Appl Math, Miusskaya Sq 4, Moscow 125047, Russia

[2] Skolkovo Inst Sci & Technol, Skolkovo Innovat Ctr, Bldg 3, Moscow 143026, Russia

[3] Int Sch Twente, Stedelijk Lyceum, Tiemeister 20, NL-7541 WG Enschede, Netherlands

[4] Univ Twente, Dept Appl Math, POB 217, NL-7500 AE Enschede, Netherlands

[5] Univ Twente, MESA Inst Nanotechnol, POB 217, NL-7500 AE Enschede, Netherlands

[6] Univ Twente, Complex Photon Syst Grp COPS, POB 217, NL-7500 AE Enschede, Netherlands

关键词: Light scattering;    Mesoscopic optics;    Disordered media;    Finite difference time domain (FDTD) methods;    Krylov subspace methods;    Exponential time integration;   

DOI  :  10.1016/j.cam.2017.12.014

来源: Elsevier

【 摘 要 】

Light incident on a layer of scattering material such as a piece of sugar or white paper forms a characteristic speckle pattern in transmission and reflection. The information hidden in the correlations of the speckle pattern with varying frequency, polarization and angle of the incident light can be exploited for applications such as biomedical imaging and high-resolution microscopy. Conventional computational models for multi-frequency optical response involve multiple solution runs of Maxwell's equations with monochromatic sources. Exponential Krylov subspace time solvers are promising candidates for improving efficiency of such models, as single monochromatic solution can be reused for the other frequencies without performing full time-domain computations at each frequency. However, we show that the straightforward implementation appears to have serious limitations. We further propose alternative ways for efficient solution through Krylov subspace methods. Our methods are based on two different splittings of the unknown solution into different parts, each of which can be computed efficiently. Experiments demonstrate a significant gain in computation time with respect to the standard solvers. (C) 2018 Elsevier B.V. All rights reserved.

【 授权许可】

Free