How to Solve 1-D Transient Heat Transfer for Constant Temperature Surfaces?

Click For Summary
SUMMARY

This discussion focuses on solving the 1-D Transient Heat Transfer equation for constant temperature surfaces using finite difference methods. The primary equation discussed is the heat conduction PDE, represented as ∂T/∂t = α ∂²T/∂x², with the initial condition T(x,0) = 200x. Two finite difference schemes are evaluated: an explicit scheme that requires adherence to the CFL condition for stability, and the Crank-Nicholson method, which is implicit and necessitates solving a linear system Ax = b. The reference book "Heat Conduction" by M. Necati Özışık is recommended for further reading.

PREREQUISITES
  • Understanding of Partial Differential Equations (PDEs)
  • Familiarity with finite difference methods for numerical analysis
  • Knowledge of stability conditions, particularly the CFL condition
  • Basic concepts of heat conduction and boundary conditions
NEXT STEPS
  • Study the separation of variables technique for solving PDEs
  • Learn about the Crank-Nicholson method for numerical solutions
  • Explore the CFL condition and its implications for stability in numerical methods
  • Read "Heat Conduction" by M. Necati Özışık for in-depth understanding
USEFUL FOR

Students and professionals in mechanical engineering, applied mathematics, and computational physics who are involved in heat transfer analysis and numerical methods for PDEs.

Link-
Messages
99
Reaction score
0
I am a little confuse here. I am trying to solve the 1-D Transient Heat Transfer for a Constant Temp. Surfaces "wall".
The PDE is
\frac{\partial T}{\partial t}=\alpha \frac{\partial ^2 T}{\partial x^2}
T(to)=200*x
I am not so sure about how I'm suppose to get the Temperature profile

any help?
-link
 
Physics news on Phys.org
I don't understand what you've written, what do you mean by "T(to)=200*x"? Is that the temperature profile at time t=0?

To solve the PDE, have you attempted the separation of variables technique?
 
This is a very simple problem in the area of heat conduction.
There are several well written books in this area. However I have followed one of those best in the field. I will suggest
Özışık (Ozisik for non-UTF viewers), M. Necati, "Heat Conduction", John-Wiley, 1993

and inquirer meant the initial condition by "T(to)=200*x" here. Btw, boundary conditions were not supplied. Most probably solution is sought for an identical Dirichlet type boundary conditions.

I'm sure the reference will help.
 
siddharth said:
I don't understand what you've written, what do you mean by "T(to)=200*x"? Is that the temperature profile at time t=0?

To solve the PDE, have you attempted the separation of variables technique?
Yes, T=T(x,t) T(x,0)=200*x for 0<=x<=.5L (symetric at 0.5L) T(0,t)=0 T(L,t)=0
No, not separation of variables..
The problem was that I was using the finite difference method to find the temperature profile in the wall at different times and I had two equations for the problem but I didn't knew which one was useful.

Here are the equations
T_m^{p+1}= Fo(T_{m+1}^p+T_{m-1}^p)+(1-2Fo)T_m^p
Subscript are for spatial nodes.
Superscript are for the time nodes.

T_m^{n+1}-T_m^n=\frac{Fo}{2}(T_{m-1}^{n+1}-2T_m^{n+1}+T_{m+1}^{n+1}+T_{m-1}^{n}-2T_m^{n}+T_{m+1}^{n})

I found that both equations are good for solving the pde using finite difference.

Thanks siddharth

PS This how should look, (axis are not labeled...:P)
 

Attachments

  • Graph2.png
    Graph2.png
    5.6 KB · Views: 513
OK. Both are different FD schemes of the problem.
The first one is explicit in time* so you shoul be aware of the CFL condition that ensures stability.

The other is the Crank-Nicholson type discretization which is stable but needs system solver. In this one you have to collect unknown values (i.e. the values for the new time step which are indexed as n+1 here) to one side and solve the resulting Ax=b type linear system. This is an implicit form(*)


* Explicit in time. Values for the new time step can be obtained by direct substitution of previously known values. Spatial derivatives are obtained from the n-th time step.

* Implicit form. Spatial derivatives involves unknown values of the n+1st time step.
 
bilgealp said:
OK. Both are different FD schemes of the problem.
The first one is explicit in time* so you shoul be aware of the CFL condition that ensures stability.

The other is the Crank-Nicholson type discretization which is stable but needs system solver. In this one you have to collect unknown values (i.e. the values for the new time step which are indexed as n+1 here) to one side and solve the resulting Ax=b type linear system. This is an implicit form(*)


* Explicit in time. Values for the new time step can be obtained by direct substitution of previously known values. Spatial derivatives are obtained from the n-th time step.

* Implicit form. Spatial derivatives involves unknown values of the n+1st time step.

Thanks a lot for the help and info.
 
My pleasure.
 

Similar threads

Replies
13
Views
3K
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
30
Views
4K
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 52 ·
2
Replies
52
Views
7K