Scribd
Upload
67 views

CFD Problem

The document summarizes solving two 1D heat conduction problems using the TDMA (Tri-Diagonal Matrix Algorithm) method. For the first problem, the temperatures at four points along a bar are calculated. The equations are set up and coefficients are determined. The system is solved in MATLAB to get the temperatures. The second problem considers heat conduction in an insulated rod with constant temperature boundaries, and the temperature distribution is calculated using n=5 control volumes. Fortran code and results are included in appendices.

Uploaded by

denizumit89
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
67 views

CFD Problem

The document summarizes solving two 1D heat conduction problems using the TDMA (Tri-Diagonal Matrix Algorithm) method. For the first problem, the temperatures at four points along a bar are calculated. The equations are set up and coefficients are determined. The system is solved in MATLAB to get the temperatures. The second problem considers heat conduction in an insulated rod with constant temperature boundaries, and the temperature distribution is calculated using n=5 control volumes. Fortran code and results are included in appendices.

Uploaded by

denizumit89
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

CFD HOMEWORK-2

Write a program that solves 1D Heat Conductivity problem with constant


temperature boundary conditions using TDMA method.
Validation Problem-1:
Consider a steady conduction in 1d bar with known temperature at x=0 (T1 = 0
K) and x=L (T6= 16 K ).
The lenght of the bar is 8.0 m, conductivity is constant at 1.5 W/K.m, A=1 m^2
and there is no heat sources generated.We would like to find the temperature at
the center of 4 equal control volumes as shown.

Solution:

d
dT
k
+S= 0
dx dx

( )

is our general equation.

In order to solve our problem lets write our basic equation;

apTp=aWTw+aETE+Su
Finite volume equation can be written as, 2 i n+1,

aiTi=biTi+1+ciTi-1+di ;
bi=k*A/(dx)e , ci= k*A/(dx)w ;
ai= bi+ ci and di=Si+dx+A; dx=length/n (n=4)
We implement the formula;

b2=1.5*1/2=1.5/2

c2=1.5*1/1=1.5

a2=4.5/2

b3=1.5*1/2=1.5/2

c3=1.5*1/2=1.5/2

a3=3.0/2

b4=1.5*1/2=1.5/2

c4=1.5*1/2=1.5/2

a4=3.0/2

b5=1.5*1/1=1.5

c5=1.5*1/2=1.5/2

a5=4.5/2

di=0 for all terms. (There is no source term.)


1

With these coefficients, our equations becomes, after multiply by 2,

4.5T2=1.5T3+3T1+0
3 T3=1.5T4+1.5T2+0
3 T4=1.5T5+1.5T3+0
4.5T5=3 T6+1.5T4+0
Thus;
0
0 T2
4.5 1.5
0
1.5
0
3
1.5
0 T3


0
1.5
3
1.5 T4
0



0
1.5 4.5 T5
0
48
We solve this equation in Matlab and our results are;

T1=0C, T2=2C, T3=6C, T4=10C, T5=14C, T6=16C.


Validation Problem-2:

Consider the problem of source-free


heat conduction in an insulated rod
whose ends are maintained at
constant temperatues of 100C and
500C respectivley. The one dimensial
problem is governed by
d
dT
k
=0 . Calculate the steady
dx dx

( )

state temperature distribution in the rood. Thermal conductivity


k=1000W/mK, cross-sectional area A is 10x10-3 .

n=5

Graphic-1 : Validation Problem-1 Temperature distribtion

Graphic-2 : Validation Problem-2 Temperature distribtion

Appendix A Fortran code for problem 1 and results


Appendix B: Fortran code for problem 2 and results

Referances
1. H.K VERSTEEG & W. MALALASEKERA, An introduction to Computational
Fluid Dynamics The Finite Volume Method,Longman Scientific &
Technical,1995
2. http://www.cae.tntech.edu/~shan/me4730%3b5730(spring2015)/

You might also like