A finite differences numerical model for buried pipe systems is presented, accounting for sensible as well as for latent heat exchanges, so as for fully three dimensional heat diffusion in soil and flexible border conditions. After description of the algorithm, extensive validation against an analytical solution as well as against several long-term monitored real scale installations will be discussed.