4 Scopus citations

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).

Original languageEnglish (US)
Pages (from-to)194-213
Number of pages20
JournalInternational Journal of Numerical Analysis and Modeling
Volume11
Issue number1
StatePublished - Jan 1 2014

Keywords

  • Evolution of probability distribution
  • Maxwell equations
  • Monte Carlo simulation
  • Polynomial chaos
  • Random interface
  • Random media
  • Stochastic partial differential equation

ASJC Scopus subject areas

  • Numerical Analysis

Fingerprint Dive into the research topics of 'Maxwell solutions in media with multiple random interfaces'. Together they form a unique fingerprint.

  • Cite this