×

Numerical Methods for Linear Diffusion Equations in the Presence of an Interface. arXiv:1310.8248

Preprint, arXiv:1310.8248 [math.NA] (2013).
Summary: 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.

MSC:

60H10 Stochastic ordinary differential equations (aspects of stochastic analysis)
65U05 Numerical methods in probability and statistics (MSC1991)
65C05 Monte Carlo methods
60J30 Processes with independent increments (MSC1991)
60E07 Infinitely divisible distributions; stable distributions
65R20 Numerical methods for integral equations
arXiv data are taken from the arXiv OAI-PMH API. If you found a mistake, please report it directly to arXiv.