Abstract This paper considers the one-dimensional model of heat conduction in solids at low temperature, the so called phonon-Bose model. The nonlinear model consists of a conservation equation for the energy density e and the heat flux Q with ∣Q∣ < e. We present a simple and accurate class of finite volume schemes for numerical simulation of heat flow in arteries. This scheme consists of predictor and corrector steps, the predictor step contains a parameter of control of the numerical diffusion of the scheme, which modulate by using limiter theory and Riemann invariant, the corrector step recovers the balance conservation equation, the scheme can compute the numerical flux corresponding the real state of solution without relying on Riemann problem solvers and it can thus be turned to order 1 in the regions where the flow has a strong variation and to order 2 in the regions where the flow is regular. The numerical test cases demonstrate high resolution of the proposed finite volume scheme (modified Rusanov) and confirm its capability to provide accurate simulations for heat flow under flow regimes with strong shocks.