I PDE - Heat Equation - Cylindrical Coordinates.

erobz
Gold Member
Messages
4,442
Reaction score
1,839
Would method of separation of variables lead to a solution to the following PDE?

$$ \frac{1}{r} \frac{ \partial}{\partial r} \left( kr \frac{ \partial T}{ \partial r}\right) = \rho c_p \frac{\partial T }{ \partial t }$$

This would be for the transient conduction of a hollow cylinder, of wall thickness ##r##, and thermal conductivity ##k##, with temperature distribution ##T(r,t)##

How do we tell if it will work out?

I imagine expanding like:

$$ \frac{k}{r} \left( \frac{ \partial T}{\partial r} + r \frac{ \partial^2 T}{\partial r^2} \right) = \rho c_p \frac{ \partial T}{\partial t}$$

Looking for a solution:

$$T(r,t) = f(r)g(t)$$

That would yield the following equation:

$$ \frac{1}{f(r)} \frac{k}{r} \left( f'(r) + r f''(r) \right) = \lambda= \rho c_p \frac{g'(t)}{g(t)} $$

You would get a damped oscillator for ##f(r)##, and exponential function for ##g(t)## from.

$$f''(r) + \frac{1}{r} f'(r) - \frac{ \lambda}{k} f(r) = 0 \tag{2}$$

$$ g'(t) - \frac{ \lambda}{\rho c_p}g(t) = 0 \tag{3} $$

When would I find out whether or not "it makes sense or works"

If at all possible, keep this as low level explanation. I'm on paper thin ice in any type of foundational understanding of PDE's ( or perhaps I have already unknowingly fallen through?)
 
Last edited:
Physics news on Phys.org
erobz said:
Would method of separation of variables lead to a solution to the following PDE?

$$ \frac{1}{r} \frac{ \partial}{\partial r} \left( kr \frac{ \partial T}{ \partial r}\right) = \rho c_p \frac{\partial T }{ \partial t }$$

This would be for the transient conduction of a hollow cylinder, of wall thickness ##r##, and thermal conductivity ##k##, with temperature distribution ##T(r,t)##

What is r? Is is a constant (the wall thickness), in which case the inner radius also needs to be given, or is it the radial coordinate?

How do we tell if it will work out?

I imagine expanding like:

$$ \frac{k}{r} \left( \frac{ \partial T}{\partial r} + r \frac{ \partial^2 T}{\partial r^2} \right) = \rho c_p \frac{ \partial T}{\partial t}$$

Looking for a solution:

$$T(r,t) = f(r)g(t)$$

That would yield the following equation:

$$ \frac{1}{f(r)} \frac{k}{r} \left( f'(r) + r f''(r) \right) = \lambda= \rho c_p \frac{g'(t)}{g(t)} $$

Better to set <br /> \frac{1}{rf(r)}(f&#039;(r) + rf&#039;&#039;(r)) = \lambda = \frac{\rho c_p}{k} \frac{g&#039;(t)}{g(t)} with <br /> \begin{split}<br /> r^2 f&#039;&#039;(r) + rf&#039;(r) - \lambda r^2 f(r) &amp;= 0 \\<br /> g&#039;(t) - \frac{k\lambda}{\rho c_p} g(t) &amp;= 0. \end{split}. Setting \lambda = -c^2 &lt; 0 and x = cr for c &gt; 0 turns the first equation into Bessel's equation with general solution f(r) = AJ_0(cr) + BY_0(cr). J_0 and Y_0 are oscillatory, with amplitutde decaying as r \to \infty. You can use them to construct a series solution \sum_{n= 0}^\infty (a_n J_0(c_nr) + b_n Y_0(c_nr))g_n(t) as you would a Fourier series; the eigenvalues c_n are determined by the boundary conditions, and the eigenfunctions are orthogonal with respect to the inner product <br /> \langle F, G \rangle = \int_a^b rF(r)G(r)\,dr. J_0 is analytic at the origin, but Y_0 blows up, so if your domain includes the origin you can discard it.

If \lambda = c^2 &gt; 0 then we have the modified Bessel's equation with general solution f(r) = AI_0(cr) + BK_0(cr). I_0 and K_0 are monotonic, with I_0 being analytic at the origin and increasing as r \to \infty; K_0 blows up at the origin and decays as r \to \infty. Again, if the domain includes the origin only I_0 is usable. You cannot use these to construct a series solution, just as you cannot use \cosh and \sinh to construct a series solution.

Here, for decaying solutions as t \to \infty, we want \lambda &gt; 0.

We should not ignore the case \lambda = 0, which leads us to the final steady state solution. The radial dependence is given by (rf&#039;)&#039; = 0 with solution f(r) = A + B\ln r.
 
:bugeye:🏳️

Minor In Mathematics here from an undergraduate degree in Mechanical Engineering Technology... I should have probably opened with that!

pasmith said:
What is r? Is is a constant (the wall thickness), in which case the inner radius also needs to be given, or is it the radial coordinate?

Anyhow, ##r## would be the radial coordinate, within the cylindrical wall. I think the wall boundaries would be a part of my set of initial conditions.

I'm sure you answered my question...I'm just not sure where. This is the embarrassment I get from watching a YouTube video on solving PDE's and thinking I could tackle it. It seems like you are saying that it could work, in other words... a solution of this type may be valid depending on the roots?

With that last line you derived the SS solution for the cylinder that is in my heat transfer textbook...So that seems like a good sign.

I wouldn't expect oscillation to show up though, I think the temperature gradient would just decay spatially, and in time, (I could be wrong), does that mean anything for which solution I should consider?

Thank you for your reply.
 
The individual eigenfunctions oscillate; a suitable series of them can be monotonic.
 
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...

Similar threads

Back
Top