Heat conduction in a layered domain with imperfect ther-mal contact interfaces is modeled by means of a system of elliptic or parabolic PDEs with suitable boundary and transmission conditions. Well-posedness of this problem is proved and a stability estimate of the solution is given.