Cylinder with heat generation, Separation of variables

Click For Summary

Discussion Overview

The discussion revolves around modeling the temperature of a thermistor connected to a constant current source using a cylindrical geometry with internal heat generation. Participants explore the application of separation of variables to solve the heat equation, considering the complexities introduced by the temperature-dependent resistance of the thermistor.

Discussion Character

  • Exploratory
  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant sets up the problem using the heat equation and proposes a form for the internal heat generation based on the thermistor's resistance.
  • Another participant suggests finding the eigenfunctions of the spatial operator and notes that the spatial operator is Sturm-Liouville, which leads to an orthogonal basis for the functional space.
  • Some participants express concern that the presence of the constant term ##\beta## complicates the separation of variables, making it difficult to isolate spatial and temporal components.
  • There is a suggestion to solve the homogeneous equation first and then find a trivial solution for the inhomogeneous equation, emphasizing that the general solution is a combination of both.
  • Participants discuss the necessity of separating variables before finding eigenfunctions and the implications of doing so on the overall solution process.
  • A later reply provides an illustrative example using a one-dimensional PDE to demonstrate the principles of Sturm-Liouville theory and the expansion of solutions in terms of eigenfunctions.

Areas of Agreement / Disagreement

Participants do not reach a consensus on the best approach to handle the ##\beta## term in the equation or the proper sequence of steps to solve the problem. Multiple competing views on how to proceed with the separation of variables and the Sturm-Liouville framework remain evident throughout the discussion.

Contextual Notes

Participants note limitations regarding the separation of variables due to the inhomogeneous term and the need for particular solutions, which complicates the overall solution process. There is also uncertainty about how to derive the Sturm-Liouville equation from the original PDE.

  • #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
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 5 ·
Replies
5
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · 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
4K