Temperature rate of change inside a cylinder

Say you have a solid cylinder which is at temperature T0 inside.

Say you then move the cylinder to a warm location (like an oven) which is at temperature T1.

Is it possible to achieve a rate equation for dT/dt=f(t) where T is the temperature exactly in the centre of the cylinder and t is time?

And therefore work out the time it would take this location to reach a temperature T2 somewhere between T0 and T1.

To keep it as general as possible I would prefer to assume the cylinder is height h and radius r.

I've been trying to use the thermal diffusivity equation:

But my working gets messy quickly, I figured as this is quite a simple shape it should be less complicated than some similar examples I've found.

EDIT:

So using my notes from a differential equations course a few years back I managed to do the following:

But I'm stuck here..

• The top and bottom also has heat flow in, so the solution would only be rotational invariant, i.e. $T(s,\theta,\phi)=T(s,\phi)$ – unsym Dec 17 '13 at 22:56
• You could probably assume that the majority of the heat flow warming the centre (which I'm interested in) will come from the top and bottom as height<<radius. – Edd Dec 17 '13 at 23:00
• Yes. You can solve the heat equation "analytically" for this geometry. Once you have the solution you can find the derivative in time at any point. Using separation of variables with $T=A(r)B(\theta)C(Z)$, since it is cylindrical geometry, you will end up with Bessel functions in $r$ and Fourier series in $\theta$ and $Z$. – SimpleLikeAnEgg Dec 17 '13 at 23:01
• @SimpleLikeAnEgg That rings a bell, I'm going to have to check back through my dusty second year physics notes. – Edd Dec 17 '13 at 23:02

The boundary condition is the tricky part. If you use a constant temperature for the surface, i.e, Dirichlet BC, the problem is considerably easier but would not represent the reality of the situation. In reality, its is convection and radiation that is responsible for the heat transfer in the exposed areas and conduction where the cylinder rest (i would assume the base from the way you have drawn it). If you are using convective boundary condition with known $h$, then you can use these charts and how to use them provided by Cengel (also attached herein).
With radiative heat transfer, its more complicated. You would need to use Stephan-Boltzmann law and somehow estimate the emmisivity of your surface, as well as the view factor. However, you could in come up with an conservative heat transfer coefficient, $h$, that could allow include the heat radiation part as well.