Cylinder with heat generation, Separation of variables

Click For Summary
The discussion centers on modeling the temperature of a thermistor connected to a constant current source, with a focus on deriving the temperature as a function of time. The setup involves a cylinder with internal heat generation dependent on temperature, leading to a heat equation that incorporates a non-linear term. Participants debate the correct approach to solve the equation, particularly the challenge posed by the constant term that complicates variable separation. Suggestions include solving the homogeneous equation first and finding a particular solution, while also emphasizing the need to establish appropriate boundary conditions. The conversation highlights the complexities of accurately setting up the problem and the importance of understanding the underlying physics for effective modeling.
  • #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 2 ·
Replies
2
Views
2K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 4 ·
Replies
4
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 35 ·
2
Replies
35
Views
5K
  • · Replies 3 ·
Replies
3
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K