Un modele d'evolution thermique d'une etoile a neutrons ayant une structure de densite fixee est developpe. Un nouvel algorithme de resolution des equations d'evolution est presente et compare a d'autres algorithmes. Les emissivites des neutrinos, les capacites de chauffage, les conductivites thermiques et le traitement du superfluide utilises dans le modele sont decrits. Les simulations numeriques sont realisees pour une etoile de masse gravitationnelle M = 1.4 M ○. , avec ou sans champs magnetiques de surface forts, et pour une etoile ayant un mecanisme de refroidissement a acceleration centrale.