Numerical Solution of Hyperbolic Heat Conduction Equation

Hyperbolic heat conduction problem is solved numerically. The explicit and implicit Euler schemes axe constructed and investigated. It is shown that the implicit Euler scheme can be used to solve efficiently paxabolic and hyperbolic heat conduction problems. This scheme is unconditionally stable for both problems. For many integration methods strong numerical oscillations are present, when the initial and boundary conditions axe discontinuous for the hyperbolic problem. In order to regularize the implicit Euler scheme, a simple linear relation between time and space steps is proposed, which automatically introduces sufficient amount of numerical viscosity. Results of numerical experiments are presented.