# Solving for Temperature Variables Using Loops

Posted 9 years ago
11393 Views
|
|
0 Total Likes
|
 I'm currently having trouble creating a code that solves for 9 temperature variables using 60 loops/iterations. Can anyone help me detect what is wrong with my code? (*Known variables*) alpha = 3.91*(10^-6); dt = 10; dx = 0.1; dy = 0.1; k = 120; h1 = 10; h2 = 100; Tinf1 = 15; Tinf2 = 15; c = 3.069*(10^7); a = 0; Ti1 = 800; Ti2 = 800; Ti3 = 800; Ti4 = 800; Ti5 = 800; Ti6 = 800; Ti7 = 800; Ti8 = 800; Ti9 = 800; (*Define functions*) Do [ solutions = NSolve[ {T1 == ((((h1*(dx/2)*(Tinf1 - Ti1)) + (h2*(dy/2)*(Tinf2 - Ti1)) + (k*(dy/2)*((T2i - Ti1)/dx)) + (k*(dx/ 2)*((Ti4 - Ti1)/dy)))*(dt) - a + Ti1)/(c*(dx/2)*(dy/2))), T2 == ((((h1* dx*(Tinf1 - Ti2)) + (k*(dy/2)*((Ti1 - Ti2)/ dx)) + (k*(dy/2)*((Ti3 - Ti2)/dx)) + (k* dx*((Ti5 - Ti2)/dy)))*(dt) - a + Ti2)/(c*dx*(dy/2))), T3 == ((((h1* dx*(Tinf1 - Ti3)) + (k*(dy/2)*((Ti2 - Ti3)/ dx)) + (k*(dy/2)*((Ti2 - Ti3)/dx)) + (k* dx*((Ti6 - Ti3)/dy)))*(dt) - a + Ti3)/(c*dx*(dy/2))), T4 == ((((h2* dy*(Tinf2 - Ti4)) + (k*(dx/2)*((Ti1 - Ti4)/ dy)) + (k*(dx/2)*((Ti7 - Ti4)/dy)) + (k* dy*((Ti5 - Ti4)/dy)))*(dt) - a + Ti4)/(c*(dx/2)*dy)), T5 == ((((k*dx*((Ti2 - Ti5)/dy)) + (k*dy*((Ti4 - Ti5)/dx)) + (k* dy*((Ti6 - Ti5)/dx)) + (k*dx*((Ti8 - Ti5)/dy)))*(dt) - a + Ti5)/(c*dx*dy)), T6 == ((((k*dx*((Ti3 - Ti6)/dy)) + (k*dy*((Ti5 - Ti6)/dx)) + (k* dy*((Ti5 - Ti6)/dx)) + (k*dx*((Ti9 - Ti6)/dy)))*(dt) - a + Ti6)/(c*dx*dy)), T7 == ((((h2*(dy/2)*(Tinf2 - Ti7)) + (k*(dx/2)*((Ti4 - Ti7)/ dy)) + (k*(dy/2)*((Ti8 - Ti7)/dx)) + (h2*(dx/ 2)*(Tinf2 - Ti7)))*(dt) - a + Ti7)/(c*(dx/2)*(dy/2))), T8 == ((((h2*dx*(Tinf2 - Ti8)) + (k* dx*((Ti5 - Ti8)/dy)) + (k*(dy/2)*((Ti7 - Ti8)/ dx)) + (k*(dy/2)*((Ti9 - Ti8)/dx)))*(dt) - a + Ti8)/(c*dx*(dy/2))), T9 == ((((h2* dx*(Tinf2 - Ti9)) + (k*(dy/2)*((Ti8 - Ti9)/ dx)) + (k*(dy/2)*((Ti8 - Ti9)/dx)) + (k* dx*((Ti6 - Ti9)/dy)))*(dt) - a + Ti9)/(c*dx*(dy/2)))}, {T1, T2, T3, T4, T5, T6, T7, T8, T9}]; Ti1 = T1 /. solutions; Ti2 = T2 /. solutions; Ti3 = T3 /. solutions; Ti4 = T4 /. solutions; Ti5 = T5 /. solutions; Ti6 = T6 /. solutions; Ti7 = T7 /. solutions; Ti8 = T8 /. solutions; Ti9 = T9 /. solutions; , {j, 60}] Print[T1 /. solutions]; Print[T2 /. solutions]; Print[T3 /. solutions]; Print[T4 /. solutions]; Print[T5 /. solutions]; Print[T6 /. solutions]; Print[T7 /. solutions]; Print[T8 /. solutions]; Print[T9 /. solutions];