Untitled

"Andreas Messner"
"ETME 4250"

"Problem statement"



"given"

A=120[m^2]
L=3[m]
k=175[W/m*K]
E_dot_gen=200000[W]
H_r=11[W/m*K]
H_l=11[W/m*K]
E_dot=(E_dot_gen)/(A*L)
T_inf_R=converttemp(C,K,80[C])
T_inf_l=converttemp(C,K,80[C])
T_initial=80 [C]


numcells=10
DELTAx=L/numcells
numnodes=numcells+1
Duplicate i=0,numcells
x[i]=i*DELTAx
End

{
"Find"

"case 1"
H_l=0

T_inf=converttemp(C,K,T_initial)



"BCl"
}

"case 2"
"BCL"
H_l*(T_inf_l-T[numcells])-k*((T[numcells-numcells]-T[numnodes-numcells])/(DELTAx))+(E_dot)*(DELTAx/2)=0


"BCI"
DUPLICATEi=1,numcells-1
(T[i-1]-2*T[i]+T[i+1])/((DELTAx)^2)+(E_dot/k)=0
End

"BCR"

0=H_r*(T_inf_R-T[numcells])-k*((T[numcells-numcells]-T[numnodes-numcells])/(DELTAx))+(E_dot)*(DELTAx/2)


T_s=T[numcells]
t_max=T[0]