We develop an algorithm for the numerical treatment of nonlinear heat conduction problems, which is adapted to the architecture of high-performance computing systems. The technique is based on a hyperbolic heat conduction model and explicit difference scheme. The new scheme provides the second-order temporal resolution of nonlinearity with an acceptable time step. Its application to plasma dynamic simulations is discussed.