Abstract
In this work, we introduce a new numerical method, the so-called cone-grid scheme for solving hyperbolic systems of conservation laws. The formulation and the theoretical analysis for this scheme are slightly simpler than other algorithms. Simultaneously, our technique leads to a more efficacious numerical scheme. The Riemann solution is the basic ingredient for this scheme. We also give the main application of this technique, namely the global weak solutions for the one-dimensional model prescribes heat conduction in solids at low temperature, which called phonon-Bose model. This system consists of a conservation equation for the energy density and the heat flux. The scheme also satisfies positivity. The cone-grid scheme was compared with exact solution by three numerical examples, where explicit solutions are known. These numerical results verify the accuracy of the proposed scheme qualitatively and quantitatively.