Analytical solution to transient heat conduction in polar coordinates with multiple layers in radial direction

Closed form analytical double-series solution is presented for the multi-dimensional unsteady heat conduction problem in polar coordinates (2-D cylindrical) with multiple layers in the radial direction. Spatially non-uniform, but time-independent, volumetric heat sources are assumed in each layer. Separation of variables method is used to obtain transient temperature distribution. In contrast to Cartesian or cylindrical (r,z) coordinates, eigenvalues in the direction perpendicular to the layers do not explicitly depend on those in the other direction. The implication of the above statement is that the imaginary eigenvalues are precluded from the solution of the problem. However, radial (transverse) eigenvalues are implicitly dependent on the angular eigenvalues through the order of the Bessel functions which constitute the radial eigenfunctions. Therefore, for each eigenvalue in the angular direction, corresponding radial eigenvalues must be obtained. Solution is valid for any combination of homogenous boundary condition of the first or second kind in the angular direction. However, inhomogeneous boundary conditions of the third kind are applied in the radial direction. Proposed solution is also applicable to multiple layers with zero inner radius. An illustrative example problem for the three-layer semi-circular annular region is solved. Results along with the isotherms are shown graphically and discussed.