Abstract
We present an efficient numerical method for simulating the scattering of electromagnetic fields by a multilayered medium with random interfaces. The elements of this algorithm, the Monte Carlo-transformed field expansion method, are (i) an interfacial problem formulation in terms of impedance-impedance operators, (ii) simulation by a high-order perturbation of surfaces approach (the transformed field expansions method), and (iii) efficient computation of the wave field for each random sample by forward and backward substitutions. Our perturbative formulation permits us to solve a sequence of linear problems featuring an operator that is deterministic, and its LU decomposition matrices can be reused, leading to significant savings in computational effort. With an extensive set of numerical examples, we demonstrate not only the robust and high-order accuracy of our scheme for small to moderate interface deformations, but also how Padé summation can be used to address large deviations.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
