On hyperbolic method for diffusion with discontinuous coefficients

It is shown that the hyperbolic method, if formulated with diffusive fluxes as additional solution variables, works directly and naturally for diffusion problems with discontinuous coefficients as well as discontinuous tangential fluxes across a material interface. A suitable formulation is given, and numerical examples are shown to demonstrate that a hyperbolic scheme can produce exact non-differentiable piecewise-linear solutions for problems with discontinuous coefficients and tangential fluxes.