Cylinder with heat generation, Separation of variables

Click For Summary
SUMMARY

This discussion focuses on modeling the temperature of a thermistor connected to a constant current source using the heat equation. The setup involves a cylinder with radius ##a## and length ##L##, experiencing internal heat generation defined by ##Q(T) = \frac{I^2 R(T)}{\pi a^2 L}##, where ##R(T)## varies with temperature. The heat equation is expressed as $$\rho C \frac{\partial T}{\partial t} = \kappa \nabla^2 T - \mu T + \beta$$, where ##\mu## and ##\beta## are derived from the resistance-temperature relationship. The primary challenge is solving for ##T(T_s,t)## using separation of variables, complicated by the presence of the ##\beta## term.

PREREQUISITES
  • Understanding of the heat equation and its components (density, heat capacity, thermal conductivity).
  • Familiarity with thermistors and their resistance-temperature characteristics.
  • Knowledge of Sturm-Liouville theory and eigenfunction expansions.
  • Basic skills in solving partial differential equations (PDEs) using separation of variables.
NEXT STEPS
  • Study the derivation and application of the Sturm-Liouville problem in the context of PDEs.
  • Learn about the method of separation of variables for solving heat equations.
  • Explore numerical methods for solving non-linear PDEs, particularly those involving heat generation.
  • Investigate boundary conditions relevant to cylindrical geometries in thermal analysis.
USEFUL FOR

Researchers, engineers, and students in thermal dynamics, particularly those working with thermistors, heat transfer, and partial differential equations.

  • #31
First, I would like to thank everyone's help with this problem. I spent some time thinking about it, and realized that I probably made it more complex than it really was. After discussing the problem with my advisor, he suggested the following solution:

Let ##T## be the temperature of the thermistor, and ##R = R_0 + \alpha T## be the resistance of the thermistor as a function of its temperature. Note that ##\alpha < 0##. Heat generated by the thermistor at constant current is ##I^2 (R_0 + \alpha T)##. Assume heat loss due to the airflow is ##h_c (T - T_0)##, where ##h_c## is the heat transfer coefficient.

Combining everything together:

$$C_p \frac{dT}{dt} = I^2 (R_0 + \alpha T) - h_c (T-T_0)$$

where C_p is the heat capacity of the thermistor. If we define ##m \equiv I^2 \alpha - h_c## and ##b \equiv I^2 R_0 + h_c T_0##, then this can be rewritten as

$$C_p \frac{dT}{dt} = mT + b$$

Note that ##m < 0## and ##b > 0##.

This has the following solution:

$$T(t) = k_1 e^{\frac{m}{C_p}t} - \frac{b}{m}$$

Again, noting that ##m<0## and ##b>0##, the temperature asymptotically approaches a positive value. Furthermore, according to www.engineeringtoolbox.com, ##h_c## can be approximated as ##h_c \approx 10.45 - v + 10v^{1/2}##, so we can also determine ##T## as a function of airflow velocity.

I'm happy with this solution. I was way overthinking it, but this was quite elegant. What do you guys think?
 
Science news on Phys.org
  • #32
16universes said:
$$C_p \frac{dT}{dt} = mT + b$$
...
I'm happy with this solution. I was way overthinking it, but this was quite elegant. What do you guys think?

Orodruin said:
...
You could then just as well formulate the problem as an ODE directly:
$$
\partial_t T = - \kappa T + \beta
$$
for some constants ##\kappa, \beta##.
...
The stationary solution to this equation is given by letting the LHS = 0 yielding
$$
\kappa T = \beta
$$
or solving a 4th order equation if you add the radiation term.
If you linearize the problem around this solution you will find a solution that exponentially approaches the stationary solution.

This is what I think ;)
 
  • #33
Ha I apologize. I was still caught up "seeing" the problem differently. It took some time to change how I was looking at it. But yes, you were correct. Thanks
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K