New Monte Carlo schemes for simulating diffusions in discontinuous media

We introduce new Monte Carlo simulation schemes for diffusions in a discontinuous media divided in subdomains with piecewise constant diffusivity. These schemes are higher order extensions of the usual schemes and take into account the two dimensional aspects of the diffusion at the interface between subdomains. This is achieved using either stochastic process techniques or an approach based on finite differences. Numerical tests on elliptic, parabolic and eigenvalue problems involving an operator in divergence form show the efficiency of these new schemes.