A multilevel approach towards unbiased sampling of random elliptic partial differential equations

Partial differential equation is a powerful tool to characterize various physics systems. In practice, measurement errors are often present and probability models are employed to account for such uncertainties. In this paper, we present a Monte Carlo scheme that yields unbiased estimators for expectations of random elliptic partial differential equations. This algorithm combines multilevel Monte Carlo [Giles, 2008] and a randomization scheme proposed by [Rhee and Glynn, 2012, Rhee and Glynn, 2013]. Furthermore, to obtain an estimator with both finite variance and finite expected computational cost, we employ higher order approximations.