Monte Carlo-transformed field expansion method for simulating electromagnetic wave scattering by multilayered random...

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.