This study addresses the challenge of accurately modeling transient heat transfer in complex multilayered cylindrical structures, which are critical in aerospace and industrial applications. The primary objective is to develop and validate a novel numerical approach using the Differential Quadrature Method (DQM) to solve the time-dependent Fourier heat conduction equation in composite cylinders with distinct material properties. While conventional methods like Finite Element Method (FEM) are widely used, they often require computationally expensive meshes. In this work, the governing equations in cylindrical coordinates were discretized using Chebyshev-Gauss-Lobatto grid points for spatial derivatives, coupled with an efficient temporal discretization scheme. This DQM formulation allows for the precise handling of interface conditions between layers without complex mesh generation. The proposed model was validated against ANSYS finite element simulations for a four-layer cylinder under constant boundary temperatures (200 °C and 300 °C). Results demonstrate excellent agreement with maximum relative errors ranging from 0% to 2.4%, proving the method’s superior accuracy and computational efficiency. Parametric investigations revealed that variations in thermal conductivity and heat capacity significantly alter the transient temperature profiles, offering key insights for the thermal design of composite insulation and piping systems.
