1. Not finding help here? Sign up for a free 30min tutor trial with Chegg Tutors
    Dismiss Notice
Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Heat transfer temp. distribution

  1. Apr 10, 2008 #1
    1. The problem statement, all variables and given/known data
    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


    2. Relevant 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

    3. 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
     

    Attached Files:

    • fin.JPG
      fin.JPG
      File size:
      17.5 KB
      Views:
      64
    Last edited: Apr 10, 2008
  2. jcsd
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Can you help with the solution or looking for help too?



Similar Discussions: Heat transfer temp. distribution
Loading...