Analytical solutions for one-dimensional time dependent multilayer heat conduction problems were developed several decades ago. Mathematical theory for such problems in more than one dimensions was also developed during that time. Several of these methods were based on separation of variable and finite integral transform. However, the application of these methods was hindered by the fact that the eigenvalue problems, which are essential for this methodology are difficult to solve. Moreover, in two and three dimensional Cartesian coordinates these eigenvalues were imaginary rendering their solutions even more difficult. It has been recently shown that similar problems in two dimensional cylindrical and spherical coordinates do not have imaginary eigenvalues. It is also helpful that the softwares which are capable of analytical manipulations are now ubiquitous. This paper discusses the methodology as well as possible application in nuclear reactors of analytical solutions of two-dimensional multilayer heat conduction in spherical and cylindrical coordinates.