# Heat transfer temp. distribution

## Homework Statement

The geometry is attached (an L shaped 2d plate)

taking del_x = del_y = 0.5mm, using finite-volume techniques, construct a matlab program to calculate steady state temps as well as heat transfer for the full 2d geometry.
Tambient = (15+273)K

## Homework Equations

use nodal equations for the 13 nodes
use gauss-seidel iterative method and write a program in matlab to find the temp field
the nodal equations are given below

## The Attempt at a Solution

T(1,1)=[Told(2,1)+Told(1,2)+2*To+q*(del_x/l)]/4; %Node 1

for i=2:64 %Node 2
T(i,1)=[Told(i+1,1)+Told(i,2)+delta_t*Told(i-1,1)+To]/4;
end

T(65,1)=[2*To+Told(65,2)+T(64,1)]/4 %Node 3

for j=2:13 %Node 4
T(65,j)=[To+Told(65,j+1)+Told(64,j)+Told(65,j-1)]/4;
end

T(65,14)=[2*To+Told(64,14)+Told(65,13)]/4; %Node 5

for i=16:64 %Node 6
T(i,14)=[Told(i+1,14)+To+Told(i-1,14)+Told(i,13)]/4;
end

T(15,14)=[0.5*Told(16,14)+0.5*Told(15,15)+Told(14,14)+Told(15,13)]/3; %Node 7

for j=15:33 %Node 8
T(15,j)=[To+Told(15,j+1)+Told(14,j)+Told(15,j-1)]/4;
end

T(15,34)=[2*To+Told(16,34)+Told(15,33)]/4; %Node 9

for i=2:14 %Node 10
T(i,34)=[Told(i+1,34)+To+Told(i-1,34)+Told(i,33)]/4;
end

T(1,34)=[Told(2,34)+2*To+Told(1,33)+q*(del_x/l)]/4; %Node 11

for j=2:33 %Node 12
T(1,j)=[Told(2,j)+Told(1,j+1)+To+Told(1,j-1)+q*(del_x/l)]/4;
end

for i=2:14
for j=14:33
for i=2:64
for j=2:13
T(i,j)=[Told(i+1,j)+Told(i,j+1)+Told(i-1,j)+Told(i,j-1)]/4;
end
end
end
end

can't get any further than this.....need desperate help

#### Attachments

• 17.5 KB Views: 358
Last edited: