An Algorithm for the Simulation of the Magnetized Neutron Star Cooling

Abstract

A model and algorithm for the cooling of the magnetized neutron stars are presented. The cooling evolution described by a system of parabolic partial differential equations with non-linear coeffcients is solved using the Alternating Direction Implicit method. A difference scheme and the preliminary results of simulations are presented.