Scribd
Upload
50 views5 pages

Heat Transfer in 2D

A computer code was developed to solve a 2D heat transfer problem using finite difference methods. The code calculates temperatures at interior nodes of a 1.5m x 1m solid body with prescribed boundary temperatures. The temperatures are presented in a table and contour plot showing higher temperatures along boundaries decreasing to 326.9K at the center.

Uploaded by

FaisalAlbader
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
50 views5 pages

Heat Transfer in 2D

A computer code was developed to solve a 2D heat transfer problem using finite difference methods. The code calculates temperatures at interior nodes of a 1.5m x 1m solid body with prescribed boundary temperatures. The temperatures are presented in a table and contour plot showing higher temperatures along boundaries decreasing to 326.9K at the center.

Uploaded by

FaisalAlbader
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

HEAT TRANSFER IN 2D

Problem Statement
Consider a steady two-dimensional heat transfer in a long solid body whose cross section
is given in Figure 1. The measured temperatures at selected points of the outer surfaces
are as shown. The thermal conductivity of the body is k = 1 W/m. oC. The dimension of
the solid body is 1.5m x 1m. Using the finite difference method with a mesh size of
x=y=10.0 cm, write a computer solution code to solve for temperatures of interior
nodes. Determine the temperature of all interior nodes.
Present the temperature values in a table.
Plot a contour plot for the temperature of the domain.

Figure 1 Plane wall specifications


The boundary conditions of the top, left and right edges are prescribed temperature of
500K. In that case, a constant temperature value is assigned to all nodes on the edges
specified above.

For the interior nodes, the energy balance equations:

Tm1, n  Tm1, n  Tm, n1  Tm, n1


Tm, n 
4
The energy balance for the nodes on the bottom edge subjected to convection of 10
W/m2-K and environmental temperature is 300K, it is shown below

y T(m 1,n)  T(m,n) T(m,n 1)  T(m,n)


q1  k q 2  k x
2 x y
y T(m1,n)  T(m,n)
q3  k q 4  x h  T  Tm,n 
2 x
 in
q node  0
y T(m 1,n)  T(m,n) T  T(m,n) y T(m 1,n)  T(m,n)
0k  k x (m,n 1) k  x h  T  Tm,n 
2 x y 2 x
x  y, divide by k and multiply by 2
y T(m 1,n)  T(m,n) T(m,n 1)  T(m,n) y T(m1,n)  T(m,n) 2x h
0 k
2 x
 2 k x
y
k
2 x

k
 T  Tm,n 
2x h
T(m1,n)  2  T(m,n 1)  T(m 1,n)  T
T(m,n)  k
2x h
4
k

A computer code using EES software is used to calculate the temperature matrix and to
plot the temperature contour over the body.
EES Code
!"Input Data"

k_con=1 [W/m-K] "thermal conductivity "


h = 10 [W/m^2-K] "heat transfer coefficient"
T_inf = 300 [K] "enviornment temperature for convection"
T_top=500 [K] "temperature of the top surface"
T_left=500 [K] "temperature of the left surface"
T_right=500 [K] "temperature of the right surface"
DELTAx=0.1[m] "DELTAx"
DELTAy=0.1[m] "DELTAy"

!"Solution Method"
"n steps" k=10
"m steps" l = 15
"vertical left edge, specified temperature"
Duplicate n=0,k
T[0,n] = 500[K]
end
"vertical right edge, specified temperature"
Duplicate n=0,k
T[15,n] = 500[K]
end
" horizontal top edge, specified temperature"
Duplicate m=1,l-1
T[m,k] = 500[K]
end

"energy balance equation for the bottom edge"


Duplicate m=1,l-1
*2T[m,1] + T[m-1,0] + T[m+1,0]+ (2*DELTAx*h/k_con)*T_inf =
(4+2*DELTAx*h/k_con)*T[m,0]
end

"all other nodes"


Duplicate m=1,l-1
Duplicate n=1,k-1
T[m,n] = (T[m-1,n]+T[m+1,n]+T[m,n-1]+T[m,n+1])/4
end
End
The array table obtained from EES for the temperature at each node is shown in Table 2
Table 2 Temperature at the nodes

T[i,10] [K] 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0 500.0
T[i,9] [K] 500.0 497.6 495.4 493.4 491.7 490.4 489.6 489.1 489.1 489.6 490.4 491.7 493.4 495.4 497.6 500.0
T[i,8] [K] 500.0 495.1 490.5 486.4 483.0 480.5 478.7 477.9 477.9 478.7 480.5 483.0 486.4 490.5 495.1 500.0
T[i,7] [K] 500.0 492.3 485.1 478.7 473.5 469.6 467.0 465.7 465.7 467.0 469.6 473.5 478.7 485.1 492.3 500.0
T[i,6] [K] 500.0 489.0 478.7 469.9 462.8 457.5 454.0 452.3 452.3 454.0 457.5 462.8 469.9 478.7 489.0 500.0
T[i,5] [K] 500.0 484.8 471.0 459.4 450.2 443.5 439.2 437.1 437.1 439.2 443.5 450.2 459.4 471.0 484.8 500.0
T[i,4] [K] 500.0 479.4 461.1 446.3 435.1 427.2 422.3 419.9 419.9 422.3 427.2 435.1 446.3 461.1 479.4 500.0
T[i,3] [K] 500.0 471.6 447.8 429.7 416.8 408.1 402.7 400.2 400.2 402.7 408.1 416.8 429.7 447.8 471.6 500.0
T[i,2] [K] 500.0 459.3 428.6 407.8 394.2 385.5 380.4 378.0 378.0 380.4 385.5 394.2 407.8 428.6 459.3 500.0
T[i,1] [K] 500.0 436.8 399.7 378.8 366.6 359.4 355.3 353.4 353.4 355.3 359.4 366.6 378.8 399.7 436.8 500.0
T[i,0] [K] 500.0 388.0 354.8 341.1 334.1 330.1 327.9 326.9 326.9 327.9 330.1 334.1 341.1 354.8 388.0 500.0

The temperature contour inside the body is shown in Figure 2

Figure 2 Temperature contour

You might also like