@Article{IJNAM-11-1, author = {Jung, C.-Y. and Kwon, B. and Mahalov, A. and T., Nguyen, B.}, title = {Maxwell Solutions in Media with Multiple Random Interfaces}, journal = {International Journal of Numerical Analysis and Modeling}, year = {2014}, volume = {11}, number = {1}, pages = {193--212}, abstract = {
A hybrid operator splitting method is developed for computations of two-dimensional transverse magnetic Maxwell equations in media with multiple random interfaces. By projecting
the solutions into the random space using the polynomial chaos (PC) projection method, the
deterministic and random parts of the solutions are solved separately.
There are two independent stages in the algorithm: the Yee scheme with domain decomposition
implemented on a staggered grid for the deterministic part and the Monte Carlo sampling in the
post-processing stage. These two stages of the algorithm are subject of computational studies.
A parallel implementation is proposed for which the computational cost grows linearly with the
number of random interfaces. Output statistics of Maxwell solutions are obtained including means,
variance and time evolution of cumulative distribution functions (CDF). The computational results
are presented for several configurations of domains with random interfaces.
The novelty of this article lies in using level set functions to characterize the random interfaces
and, under reasonable assumptions on the random interfaces (see Figure 1), the dimensionality
issue from the PC expansions is resolved (see Sections 1.1.2 and 1.2).