Paper detail

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface

We consider numerical methods for linear parabolic equations in one spatial dimension having piecewise constant diffusion coefficients defined by a one parameter family of interface conditions at the discontinuity. We construct immersed interface finite element methods for an alternative formulation of the original deterministic diffusion problem in which the interface condition is recast as a natural condition on the interfacial flux for which the given operator is self adjoint. An Euler-Maruyama method is developed for the stochastic differential equation corresponding to the alternative divergence formulation of the equation having a discontinuous coefficient and a one-parameter family of interface conditions. We then prove convergence estimates for the Euler scheme. The main goal is to develop numerical schemes that can accommodate specification of any one of the possible interface conditions, and to illustrate the implementation for each of the deterministic and stochastic formulations, respectively. The issues pertaining to speed-ups of the numerical schemes are left to future work.