A non-isothermal two-dimensional general rate model is formulated and analytically solved to analyze the effects of temperature changes inside liquid chromatographic columns of cylindrical geometry. The model equations form a system of convection-diffusion partial differential equations. The finite-Hankel transformation, the Laplace transformation, the eigen-decomposition technique and a conventional solution technique of ordinary differential equations are used to solve the equations of the model. The coupling between concentration and temperature fronts is demonstrated and important parameters that affect the performance of the column are evaluated. To find the ranges of validity of our analytical results, a semi-discrete high resolution finite volume method is applied to solve the same system of equations for both linear and nonlinear isotherms. The results of this contribution can be helpful to optimize non-isothermal liquid chromatographic processes in which both radial and axial gradients occur.