We implement a numerical model to analyze the thermodynamics that can be used for navigation and orbit determination of future space missions to small bodies. Unlike past models that use a spherical homogeneous model, our model includes the real non-spherical shape of asteroid 433 Eros. The surface temperature map is expressed as a function of latitude and longitude directions for different initial isothermal temperatures. The lack of spherical symmetry is modeled along with the long and short term temperature variations as a consequence of the orbital and spin axis rotations of the body.