Spreadsheet Problem Solving for ChE Faculty

Workshop - Ordinary Differential Equations - Countercurrent Tube-in-tube Heat Exchanger

Solution by Euler's Method and Goal Seek (or Solver)
Note: This solution make extensive
use of the IF function to allow for different
L and dz values and the LOOKUP
function to compute the overall enthalpy
changes for the hot and cold fluid
(must be equal for correct solution)

Brice Carnahan
Enter Data into Outlined Cells Only
Note: To insure enthalpy conservation, only hi (not both hi and ho) should be specified.
ho is computed as hi*DiC/DoCo
Inner tube (Cold)
DiC(inches)
DiC
DoC(inches)
DoC
Ai
2
hi(BTU/(ft hr F)
hi
TCin(oC)
vC
rhoC
AflowC
wC
CpC

Step
0
=IF(A34<Nsteps,A34+1," ")
=IF(A35<Nsteps,A35+1," ")
=IF(A36<Nsteps,A36+1," ")
=IF(A37<Nsteps,A37+1," ")
=IF(A38<Nsteps,A38+1," ")
=IF(A39<Nsteps,A39+1," ")
=IF(A40<Nsteps,A40+1," ")
=IF(A41<Nsteps,A41+1," ")
=IF(A42<Nsteps,A42+1," ")
=IF(A43<Nsteps,A43+1," ")
=IF(A44<Nsteps,A44+1," ")
=IF(A45<Nsteps,A45+1," ")
=IF(A46<Nsteps,A46+1," ")
=IF(A47<Nsteps,A47+1," ")
=IF(A48<Nsteps,A48+1," ")
=IF(A49<Nsteps,A49+1," ")
=IF(A50<Nsteps,A50+1," ")
=IF(A51<Nsteps,A51+1," ")
=IF(A52<Nsteps,A52+1," ")
=IF(A53<Nsteps,A53+1," ")
=IF(A54<Nsteps,A54+1," ")
=IF(A55<Nsteps,A55+1," ")
=IF(A56<Nsteps,A56+1," ")
=IF(A57<Nsteps,A57+1," ")
=IF(A58<Nsteps,A58+1," ")
=IF(A59<Nsteps,A59+1," ")
=IF(A60<Nsteps,A60+1," ")
=IF(A61<Nsteps,A61+1," ")
=IF(A62<Nsteps,A62+1," ")
=IF(A63<Nsteps,A63+1," ")
=IF(A64<Nsteps,A64+1," ")
=IF(A65<Nsteps,A65+1," ")
=IF(A66<Nsteps,A66+1," ")
=IF(A67<Nsteps,A67+1," ")
=IF(A68<Nsteps,A68+1," ")
=IF(A69<Nsteps,A69+1," ")
=IF(A70<Nsteps,A70+1," ")
=IF(A71<Nsteps,A71+1," ")
=IF(A72<Nsteps,A72+1," ")

0.76
=B15*0.0254
1
=B17*0.0254
=PI()*DiC
75
=B20*5.673439
10
0.1
1000
=PI()*(DiC^2/4)
=vC*AflowC*rhoC
4200

Distance
0
=IF(B34<L,B34+dz," ")
=IF(B35<L,B35+dz," ")
=IF(B36<L,B36+dz," ")
=IF(B37<L,B37+dz," ")
=IF(B38<L,B38+dz," ")
=IF(B39<L,B39+dz," ")
=IF(B40<L,B40+dz," ")
=IF(B41<L,B41+dz," ")
=IF(B42<L,B42+dz," ")
=IF(B43<L,B43+dz," ")
=IF(B44<L,B44+dz," ")
=IF(B45<L,B45+dz," ")
=IF(B46<L,B46+dz," ")
=IF(B47<L,B47+dz," ")
=IF(B48<L,B48+dz," ")
=IF(B49<L,B49+dz," ")
=IF(B50<L,B50+dz," ")
=IF(B51<L,B51+dz," ")
=IF(B52<L,B52+dz," ")
=IF(B53<L,B53+dz," ")
=IF(B54<L,B54+dz," ")
=IF(B55<L,B55+dz," ")
=IF(B56<L,B56+dz," ")
=IF(B57<L,B57+dz," ")
=IF(B58<L,B58+dz," ")
=IF(B59<L,B59+dz," ")
=IF(B60<L,B60+dz," ")
=IF(B61<L,B61+dz," ")
=IF(B62<L,B62+dz," ")
=IF(B63<L,B63+dz," ")
=IF(B64<L,B64+dz," ")
=IF(B65<L,B65+dz," ")
=IF(B66<L,B66+dz," ")
=IF(B67<L,B67+dz," ")
=IF(B68<L,B68+dz," ")
=IF(B69<L,B69+dz," ")
=IF(B70<L,B70+dz," ")
=IF(B71<L,B71+dz," ")
=IF(B72<L,B72+dz," ")

Goal Seek or Solver is used to make

TH at z = L equal to THin by adjusting
THout.

Outer tube (Hot)

in
m
in
m
2
m per m
2
BTU/(ft hr F)
J/(m2 s C)
C
m/s
kg/m3
m2
kg/s
J/(kg C)
m
m

DiH(inches)
DiH
DoH(inches)
DoH
Ao
ho(BTU/(ft2 hr F)
ho
THin(oC)
vH
rhoH
AflowH
wH
CpH

1.76
=E15*0.0254
2
=E17*0.0254
=PI()*DoC
=B20*B15/B17
=E20*5.673439
50
0.1
1000
=PI()*(DiH^2/4-DoC^2/4)
=vH*AflowH*rhoH
4200

Guess THout 43.2383239474323

TC
=TCin
=IF(B35<=L,C34+dz*E34," ")
=IF(B36<=L,C35+dz*E35," ")
=IF(B37<=L,C36+dz*E36," ")
=IF(B38<=L,C37+dz*E37," ")
=IF(B39<=L,C38+dz*E38," ")
=IF(B40<=L,C39+dz*E39," ")
=IF(B41<=L,C40+dz*E40," ")
=IF(B42<=L,C41+dz*E41," ")
=IF(B43<=L,C42+dz*E42," ")
=IF(B44<=L,C43+dz*E43," ")
=IF(B45<=L,C44+dz*E44," ")
=IF(B46<=L,C45+dz*E45," ")
=IF(B47<=L,C46+dz*E46," ")
=IF(B48<=L,C47+dz*E47," ")
=IF(B49<=L,C48+dz*E48," ")
=IF(B50<=L,C49+dz*E49," ")
=IF(B51<=L,C50+dz*E50," ")
=IF(B52<=L,C51+dz*E51," ")
=IF(B53<=L,C52+dz*E52," ")
=IF(B54<=L,C53+dz*E53," ")
=IF(B55<=L,C54+dz*E54," ")
=IF(B56<=L,C55+dz*E55," ")
=IF(B57<=L,C56+dz*E56," ")
=IF(B58<=L,C57+dz*E57," ")
=IF(B59<=L,C58+dz*E58," ")
=IF(B60<=L,C59+dz*E59," ")
=IF(B61<=L,C60+dz*E60," ")
=IF(B62<=L,C61+dz*E61," ")
=IF(B63<=L,C62+dz*E62," ")
=IF(B64<=L,C63+dz*E63," ")
=IF(B65<=L,C64+dz*E64," ")
=IF(B66<=L,C65+dz*E65," ")
=IF(B67<=L,C66+dz*E66," ")
=IF(B68<=L,C67+dz*E67," ")
=IF(B69<=L,C68+dz*E68," ")
=IF(B70<=L,C69+dz*E69," ")
=IF(B71<=L,C70+dz*E70," ")
=IF(B72<=L,C71+dz*E71," ")
=IF(B73<=L,C72+dz*E72," ")

TH
=Guess
=IF(B35<=L,D34+dz*F34," ")
=IF(B36<=L,D35+dz*F35," ")
=IF(B37<=L,D36+dz*F36," ")
=IF(B38<=L,D37+dz*F37," ")
=IF(B39<=L,D38+dz*F38," ")
=IF(B40<=L,D39+dz*F39," ")
=IF(B41<=L,D40+dz*F40," ")
=IF(B42<=L,D41+dz*F41," ")
=IF(B43<=L,D42+dz*F42," ")
=IF(B44<=L,D43+dz*F43," ")
=IF(B45<=L,D44+dz*F44," ")
=IF(B46<=L,D45+dz*F45," ")
=IF(B47<=L,D46+dz*F46," ")
=IF(B48<=L,D47+dz*F47," ")
=IF(B49<=L,D48+dz*F48," ")
=IF(B50<=L,D49+dz*F49," ")
=IF(B51<=L,D50+dz*F50," ")
=IF(B52<=L,D51+dz*F51," ")
=IF(B53<=L,D52+dz*F52," ")
=IF(B54<=L,D53+dz*F53," ")
=IF(B55<=L,D54+dz*F54," ")
=IF(B56<=L,D55+dz*F55," ")
=IF(B57<=L,D56+dz*F56," ")
=IF(B58<=L,D57+dz*F57," ")
=IF(B59<=L,D58+dz*F58," ")
=IF(B60<=L,D59+dz*F59," ")
=IF(B61<=L,D60+dz*F60," ")
=IF(B62<=L,D61+dz*F61," ")
=IF(B63<=L,D62+dz*F62," ")
=IF(B64<=L,D63+dz*F63," ")
=IF(B65<=L,D64+dz*F64," ")
=IF(B66<=L,D65+dz*F65," ")
=IF(B67<=L,D66+dz*F66," ")
=IF(B68<=L,D67+dz*F67," ")
=IF(B69<=L,D68+dz*F68," ")
=IF(B70<=L,D69+dz*F69," ")
=IF(B71<=L,D70+dz*F70," ")
=IF(B72<=L,D71+dz*F71," ")
=IF(B73<=L,D72+dz*F72," ")

dTC/dz
=hi*Ai/(wC*CpC)*(D34-C34)
=IF(B35<=L,hi*Ai/(wC*CpC)*(D35-C35)," ")
=IF(B36<=L,hi*Ai/(wC*CpC)*(D36-C36)," ")
=IF(B37<=L,hi*Ai/(wC*CpC)*(D37-C37)," ")
=IF(B38<=L,hi*Ai/(wC*CpC)*(D38-C38)," ")
=IF(B39<=L,hi*Ai/(wC*CpC)*(D39-C39)," ")
=IF(B40<=L,hi*Ai/(wC*CpC)*(D40-C40)," ")
=IF(B41<=L,hi*Ai/(wC*CpC)*(D41-C41)," ")
=IF(B42<=L,hi*Ai/(wC*CpC)*(D42-C42)," ")
=IF(B43<=L,hi*Ai/(wC*CpC)*(D43-C43)," ")
=IF(B44<=L,hi*Ai/(wC*CpC)*(D44-C44)," ")
=IF(B45<=L,hi*Ai/(wC*CpC)*(D45-C45)," ")
=IF(B46<=L,hi*Ai/(wC*CpC)*(D46-C46)," ")
=IF(B47<=L,hi*Ai/(wC*CpC)*(D47-C47)," ")
=IF(B48<=L,hi*Ai/(wC*CpC)*(D48-C48)," ")
=IF(B49<=L,hi*Ai/(wC*CpC)*(D49-C49)," ")
=IF(B50<=L,hi*Ai/(wC*CpC)*(D50-C50)," ")
=IF(B51<=L,hi*Ai/(wC*CpC)*(D51-C51)," ")
=IF(B52<=L,hi*Ai/(wC*CpC)*(D52-C52)," ")
=IF(B53<=L,hi*Ai/(wC*CpC)*(D53-C53)," ")
=IF(B54<=L,hi*Ai/(wC*CpC)*(D54-C54)," ")
=IF(B55<=L,hi*Ai/(wC*CpC)*(D55-C55)," ")
=IF(B56<=L,hi*Ai/(wC*CpC)*(D56-C56)," ")
=IF(B57<=L,hi*Ai/(wC*CpC)*(D57-C57)," ")
=IF(B58<=L,hi*Ai/(wC*CpC)*(D58-C58)," ")
=IF(B59<=L,hi*Ai/(wC*CpC)*(D59-C59)," ")
=IF(B60<=L,hi*Ai/(wC*CpC)*(D60-C60)," ")
=IF(B61<=L,hi*Ai/(wC*CpC)*(D61-C61)," ")
=IF(B62<=L,hi*Ai/(wC*CpC)*(D62-C62)," ")
=IF(B63<=L,hi*Ai/(wC*CpC)*(D63-C63)," ")
=IF(B64<=L,hi*Ai/(wC*CpC)*(D64-C64)," ")
=IF(B65<=L,hi*Ai/(wC*CpC)*(D65-C65)," ")
=IF(B66<=L,hi*Ai/(wC*CpC)*(D66-C66)," ")
=IF(B67<=L,hi*Ai/(wC*CpC)*(D67-C67)," ")
=IF(B68<=L,hi*Ai/(wC*CpC)*(D68-C68)," ")
=IF(B69<=L,hi*Ai/(wC*CpC)*(D69-C69)," ")
=IF(B70<=L,hi*Ai/(wC*CpC)*(D70-C70)," ")
=IF(B71<=L,hi*Ai/(wC*CpC)*(D71-C71)," ")
=IF(B72<=L,hi*Ai/(wC*CpC)*(D72-C72)," ")
=IF(B73<=L,hi*Ai/(wC*CpC)*(D73-C73)," ")

in
m
in
m
2
m per m
2
BTU/(ft hr F)
J/(m2 s C)
C
m/s
kg/m3
m2
kg/s
J/(kg C)

Euler Method parameters

dz 0.1
L 5
Nsteps =L/dz

for Goal Seek

dTH/dz
=ho*Ao/(wH*CpH)*(D34-C34)
Computed
=IF(B35<=L,ho*Ao/(wH*CpH)*(D35-C35)," ")
TCout
=IF(B36<=L,ho*Ao/(wH*CpH)*(D36-C36)," ")
THin
=IF(B37<=L,ho*Ao/(wH*CpH)*(D37-C37)," ")
=IF(B38<=L,ho*Ao/(wH*CpH)*(D38-C38)," ")Cold Fluid DH
=IF(B39<=L,ho*Ao/(wH*CpH)*(D39-C39)," ") Hot Fluid DH
=IF(B40<=L,ho*Ao/(wH*CpH)*(D40-C40)," ")
Error
=IF(B41<=L,ho*Ao/(wH*CpH)*(D41-C41)," ")
=IF(B42<=L,ho*Ao/(wH*CpH)*(D42-C42)," ")
=IF(B43<=L,ho*Ao/(wH*CpH)*(D43-C43)," ")
=IF(B44<=L,ho*Ao/(wH*CpH)*(D44-C44)," ")
=IF(B45<=L,ho*Ao/(wH*CpH)*(D45-C45)," ")
=IF(B46<=L,ho*Ao/(wH*CpH)*(D46-C46)," ")
=IF(B47<=L,ho*Ao/(wH*CpH)*(D47-C47)," ")
=IF(B48<=L,ho*Ao/(wH*CpH)*(D48-C48)," ")
=IF(B49<=L,ho*Ao/(wH*CpH)*(D49-C49)," ")
=IF(B50<=L,ho*Ao/(wH*CpH)*(D50-C50)," ")
=IF(B51<=L,ho*Ao/(wH*CpH)*(D51-C51)," ")
=IF(B52<=L,ho*Ao/(wH*CpH)*(D52-C52)," ")
=IF(B53<=L,ho*Ao/(wH*CpH)*(D53-C53)," ")
=IF(B54<=L,ho*Ao/(wH*CpH)*(D54-C54)," ")
=IF(B55<=L,ho*Ao/(wH*CpH)*(D55-C55)," ")
=IF(B56<=L,ho*Ao/(wH*CpH)*(D56-C56)," ")
=IF(B57<=L,ho*Ao/(wH*CpH)*(D57-C57)," ")
=IF(B58<=L,ho*Ao/(wH*CpH)*(D58-C58)," ")
=IF(B59<=L,ho*Ao/(wH*CpH)*(D59-C59)," ")
=IF(B60<=L,ho*Ao/(wH*CpH)*(D60-C60)," ")
=IF(B61<=L,ho*Ao/(wH*CpH)*(D61-C61)," ")
=IF(B62<=L,ho*Ao/(wH*CpH)*(D62-C62)," ")
=IF(B63<=L,ho*Ao/(wH*CpH)*(D63-C63)," ")
=IF(B64<=L,ho*Ao/(wH*CpH)*(D64-C64)," ")
=IF(B65<=L,ho*Ao/(wH*CpH)*(D65-C65)," ")
=IF(B66<=L,ho*Ao/(wH*CpH)*(D66-C66)," ")
=IF(B67<=L,ho*Ao/(wH*CpH)*(D67-C67)," ")
=IF(B68<=L,ho*Ao/(wH*CpH)*(D68-C68)," ")
=IF(B69<=L,ho*Ao/(wH*CpH)*(D69-C69)," ")
=IF(B70<=L,ho*Ao/(wH*CpH)*(D70-C70)," ")
=IF(B71<=L,ho*Ao/(wH*CpH)*(D71-C71)," ")
=IF(B72<=L,ho*Ao/(wH*CpH)*(D72-C72)," ")
=IF(B73<=L,ho*Ao/(wH*CpH)*(D73-C73)," ")

Check Total
Heat fluxes
=LOOKUP(Nsteps,A34:A134,C34:C134)
=LOOKUP(Nsteps,A34:A134,D34:D134)
=wC*CpC*(H35-TCin)
=- wH*CpH*(THin-Guess)
=H38+H39

J/s
J/s
J/s

=IF(A73<Nsteps,A73+1," ") =IF(B73<L,B73+dz," ") =IF(B74<=L,C73+dz*E73," ")

=IF(B74<=L,D73+dz*F73," ")
=IF(B74<=L,hi*Ai/(wC*CpC)*(D74-C74)," ") =IF(B74<=L,ho*Ao/(wH*CpH)*(D74-C74)," ")
=IF(A74<Nsteps,A74+1," ") =IF(B74<L,B74+dz," ") =IF(B75<=L,C74+dz*E74," ")
=IF(B75<=L,D74+dz*F74," ")
=IF(B75<=L,hi*Ai/(wC*CpC)*(D75-C75)," ") =IF(B75<=L,ho*Ao/(wH*CpH)*(D75-C75)," ")
=IF(A75<Nsteps,A75+1," ") =IF(B75<L,B75+dz," ") =IF(B76<=L,C75+dz*E75," ")
=IF(B76<=L,D75+dz*F75," ")
=IF(B76<=L,hi*Ai/(wC*CpC)*(D76-C76)," ") =IF(B76<=L,ho*Ao/(wH*CpH)*(D76-C76)," ")
=IF(A76<Nsteps,A76+1," ") =IF(B76<L,B76+dz," ") =IF(B77<=L,C76+dz*E76," ")
=IF(B77<=L,D76+dz*F76," ")
=IF(B77<=L,hi*Ai/(wC*CpC)*(D77-C77)," ") =IF(B77<=L,ho*Ao/(wH*CpH)*(D77-C77)," ")
=IF(A77<Nsteps,A77+1," ") =IF(B77<L,B77+dz," ") =IF(B78<=L,C77+dz*E77," ")
=IF(B78<=L,D77+dz*F77," ")
=IF(B78<=L,hi*Ai/(wC*CpC)*(D78-C78)," ") =IF(B78<=L,ho*Ao/(wH*CpH)*(D78-C78)," ")
=IF(A78<Nsteps,A78+1," ") =IF(B78<L,B78+dz," ") =IF(B79<=L,C78+dz*E78," ")
=IF(B79<=L,D78+dz*F78," ")
=IF(B79<=L,hi*Ai/(wC*CpC)*(D79-C79)," ") =IF(B79<=L,ho*Ao/(wH*CpH)*(D79-C79)," ")
=IF(A79<Nsteps,A79+1," ") =IF(B79<L,B79+dz," ") =IF(B80<=L,C79+dz*E79," ")
=IF(B80<=L,D79+dz*F79," ")
=IF(B80<=L,hi*Ai/(wC*CpC)*(D80-C80)," ") =IF(B80<=L,ho*Ao/(wH*CpH)*(D80-C80)," ")
=IF(A80<Nsteps,A80+1," ") =IF(B80<L,B80+dz," ") =IF(B81<=L,C80+dz*E80," ")
=IF(B81<=L,D80+dz*F80," ")
=IF(B81<=L,hi*Ai/(wC*CpC)*(D81-C81)," ") =IF(B81<=L,ho*Ao/(wH*CpH)*(D81-C81)," ")
=IF(A81<Nsteps,A81+1," ") =IF(B81<L,B81+dz," ") =IF(B82<=L,C81+dz*E81," ")
=IF(B82<=L,D81+dz*F81," ")
=IF(B82<=L,hi*Ai/(wC*CpC)*(D82-C82)," ") =IF(B82<=L,ho*Ao/(wH*CpH)*(D82-C82)," ")
=IF(A82<Nsteps,A82+1," ") =IF(B82<L,B82+dz," ") =IF(B83<=L,C82+dz*E82," ")
=IF(B83<=L,D82+dz*F82," ")
=IF(B83<=L,hi*Ai/(wC*CpC)*(D83-C83)," ") =IF(B83<=L,ho*Ao/(wH*CpH)*(D83-C83)," ")
=IF(A83<Nsteps,A83+1," ") =IF(B83<L,B83+dz," ") =IF(B84<=L,C83+dz*E83," ")
=IF(B84<=L,D83+dz*F83," ")
=IF(B84<=L,hi*Ai/(wC*CpC)*(D84-C84)," ") =IF(B84<=L,ho*Ao/(wH*CpH)*(D84-C84)," ")
=IF(A84<Nsteps,A84+1," ") =IF(B84<L,B84+dz," ") =IF(B85<=L,C84+dz*E84," ")
=IF(B85<=L,D84+dz*F84," ")
=IF(B85<=L,hi*Ai/(wC*CpC)*(D85-C85)," ") =IF(B85<=L,ho*Ao/(wH*CpH)*(D85-C85)," ")
=IF(A85<Nsteps,A85+1," ") =IF(B85<L,B85+dz," ") =IF(B86<=L,C85+dz*E85," ")
=IF(B86<=L,D85+dz*F85," ")
=IF(B86<=L,hi*Ai/(wC*CpC)*(D86-C86)," ") =IF(B86<=L,ho*Ao/(wH*CpH)*(D86-C86)," ")
=IF(A86<Nsteps,A86+1," ") =IF(B86<L,B86+dz," ") =IF(B87<=L,C86+dz*E86," ")
=IF(B87<=L,D86+dz*F86," ")
=IF(B87<=L,hi*Ai/(wC*CpC)*(D87-C87)," ") =IF(B87<=L,ho*Ao/(wH*CpH)*(D87-C87)," ")
=IF(A87<Nsteps,A87+1," ") =IF(B87<L,B87+dz," ") =IF(B88<=L,C87+dz*E87," ")
=IF(B88<=L,D87+dz*F87," ")
=IF(B88<=L,hi*Ai/(wC*CpC)*(D88-C88)," ") =IF(B88<=L,ho*Ao/(wH*CpH)*(D88-C88)," ")
=IF(A88<Nsteps,A88+1," ") =IF(B88<L,B88+dz," ") =IF(B89<=L,C88+dz*E88," ")
=IF(B89<=L,D88+dz*F88," ")
=IF(B89<=L,hi*Ai/(wC*CpC)*(D89-C89)," ") =IF(B89<=L,ho*Ao/(wH*CpH)*(D89-C89)," ")
=IF(A89<Nsteps,A89+1," ") =IF(B89<L,B89+dz," ") =IF(B90<=L,C89+dz*E89," ")
=IF(B90<=L,D89+dz*F89," ")
=IF(B90<=L,hi*Ai/(wC*CpC)*(D90-C90)," ") =IF(B90<=L,ho*Ao/(wH*CpH)*(D90-C90)," ")
=IF(A90<Nsteps,A90+1," ") =IF(B90<L,B90+dz," ") =IF(B91<=L,C90+dz*E90," ")
=IF(B91<=L,D90+dz*F90," ")
=IF(B91<=L,hi*Ai/(wC*CpC)*(D91-C91)," ") =IF(B91<=L,ho*Ao/(wH*CpH)*(D91-C91)," ")
=IF(A91<Nsteps,A91+1," ") =IF(B91<L,B91+dz," ") =IF(B92<=L,C91+dz*E91," ")
=IF(B92<=L,D91+dz*F91," ")
=IF(B92<=L,hi*Ai/(wC*CpC)*(D92-C92)," ") =IF(B92<=L,ho*Ao/(wH*CpH)*(D92-C92)," ")
=IF(A92<Nsteps,A92+1," ") =IF(B92<L,B92+dz," ") =IF(B93<=L,C92+dz*E92," ")
=IF(B93<=L,D92+dz*F92," ")
=IF(B93<=L,hi*Ai/(wC*CpC)*(D93-C93)," ") =IF(B93<=L,ho*Ao/(wH*CpH)*(D93-C93)," ")
=IF(A93<Nsteps,A93+1," ") =IF(B93<L,B93+dz," ") =IF(B94<=L,C93+dz*E93," ")
=IF(B94<=L,D93+dz*F93," ")
=IF(B94<=L,hi*Ai/(wC*CpC)*(D94-C94)," ") =IF(B94<=L,ho*Ao/(wH*CpH)*(D94-C94)," ")
=IF(A94<Nsteps,A94+1," ") =IF(B94<L,B94+dz," ") =IF(B95<=L,C94+dz*E94," ")
=IF(B95<=L,D94+dz*F94," ")
=IF(B95<=L,hi*Ai/(wC*CpC)*(D95-C95)," ") =IF(B95<=L,ho*Ao/(wH*CpH)*(D95-C95)," ")
=IF(A95<Nsteps,A95+1," ") =IF(B95<L,B95+dz," ") =IF(B96<=L,C95+dz*E95," ")
=IF(B96<=L,D95+dz*F95," ")
=IF(B96<=L,hi*Ai/(wC*CpC)*(D96-C96)," ") =IF(B96<=L,ho*Ao/(wH*CpH)*(D96-C96)," ")
=IF(A96<Nsteps,A96+1," ") =IF(B96<L,B96+dz," ") =IF(B97<=L,C96+dz*E96," ")
=IF(B97<=L,D96+dz*F96," ")
=IF(B97<=L,hi*Ai/(wC*CpC)*(D97-C97)," ") =IF(B97<=L,ho*Ao/(wH*CpH)*(D97-C97)," ")
=IF(A97<Nsteps,A97+1," ") =IF(B97<L,B97+dz," ") =IF(B98<=L,C97+dz*E97," ")
=IF(B98<=L,D97+dz*F97," ")
=IF(B98<=L,hi*Ai/(wC*CpC)*(D98-C98)," ") =IF(B98<=L,ho*Ao/(wH*CpH)*(D98-C98)," ")
=IF(A98<Nsteps,A98+1," ") =IF(B98<L,B98+dz," ") =IF(B99<=L,C98+dz*E98," ")
=IF(B99<=L,D98+dz*F98," ")
=IF(B99<=L,hi*Ai/(wC*CpC)*(D99-C99)," ") =IF(B99<=L,ho*Ao/(wH*CpH)*(D99-C99)," ")
=IF(A99<Nsteps,A99+1," ") =IF(B99<L,B99+dz," ") =IF(B100<=L,C99+dz*E99," ") =IF(B100<=L,D99+dz*F99," ") =IF(B100<=L,hi*Ai/(wC*CpC)*(D100-C100)," =IF(B100<=L,ho*Ao/(wH*CpH)*(D100-C100),"
")
")
=IF(A100<Nsteps,A100+1," ") =IF(B100<L,B100+dz," ") =IF(B101<=L,C100+dz*E100," ")=IF(B101<=L,D100+dz*F100," ")=IF(B101<=L,hi*Ai/(wC*CpC)*(D101-C101)," =IF(B101<=L,ho*Ao/(wH*CpH)*(D101-C101),"
")
")
=IF(A101<Nsteps,A101+1," ") =IF(B101<L,B101+dz," ") =IF(B102<=L,C101+dz*E101," ")=IF(B102<=L,D101+dz*F101," ")=IF(B102<=L,hi*Ai/(wC*CpC)*(D102-C102)," =IF(B102<=L,ho*Ao/(wH*CpH)*(D102-C102),"
")
")
=IF(A102<Nsteps,A102+1," ") =IF(B102<L,B102+dz," ") =IF(B103<=L,C102+dz*E102," ")=IF(B103<=L,D102+dz*F102," ")=IF(B103<=L,hi*Ai/(wC*CpC)*(D103-C103)," =IF(B103<=L,ho*Ao/(wH*CpH)*(D103-C103),"
")
")
=IF(A103<Nsteps,A103+1," ") =IF(B103<L,B103+dz," ") =IF(B104<=L,C103+dz*E103," ")=IF(B104<=L,D103+dz*F103," ")=IF(B104<=L,hi*Ai/(wC*CpC)*(D104-C104)," =IF(B104<=L,ho*Ao/(wH*CpH)*(D104-C104),"
")
")
=IF(A104<Nsteps,A104+1," ") =IF(B104<L,B104+dz," ") =IF(B105<=L,C104+dz*E104," ")=IF(B105<=L,D104+dz*F104," ")=IF(B105<=L,hi*Ai/(wC*CpC)*(D105-C105)," =IF(B105<=L,ho*Ao/(wH*CpH)*(D105-C105),"
")
")
=IF(A105<Nsteps,A105+1," ") =IF(B105<L,B105+dz," ") =IF(B106<=L,C105+dz*E105," ")=IF(B106<=L,D105+dz*F105," ")=IF(B106<=L,hi*Ai/(wC*CpC)*(D106-C106)," =IF(B106<=L,ho*Ao/(wH*CpH)*(D106-C106),"
")
")
=IF(A106<Nsteps,A106+1," ") =IF(B106<L,B106+dz," ") =IF(B107<=L,C106+dz*E106," ")=IF(B107<=L,D106+dz*F106," ")=IF(B107<=L,hi*Ai/(wC*CpC)*(D107-C107)," =IF(B107<=L,ho*Ao/(wH*CpH)*(D107-C107),"
")
")
=IF(A107<Nsteps,A107+1," ") =IF(B107<L,B107+dz," ") =IF(B108<=L,C107+dz*E107," ")=IF(B108<=L,D107+dz*F107," ")=IF(B108<=L,hi*Ai/(wC*CpC)*(D108-C108)," =IF(B108<=L,ho*Ao/(wH*CpH)*(D108-C108),"
")
")
=IF(A108<Nsteps,A108+1," ") =IF(B108<L,B108+dz," ") =IF(B109<=L,C108+dz*E108," ")=IF(B109<=L,D108+dz*F108," ")=IF(B109<=L,hi*Ai/(wC*CpC)*(D109-C109)," =IF(B109<=L,ho*Ao/(wH*CpH)*(D109-C109),"
")
")
=IF(A109<Nsteps,A109+1," ") =IF(B109<L,B109+dz," ") =IF(B110<=L,C109+dz*E109," ")=IF(B110<=L,D109+dz*F109," ")=IF(B110<=L,hi*Ai/(wC*CpC)*(D110-C110)," =IF(B110<=L,ho*Ao/(wH*CpH)*(D110-C110),"
")
")
=IF(A110<Nsteps,A110+1," ") =IF(B110<L,B110+dz," ") =IF(B111<=L,C110+dz*E110," ")=IF(B111<=L,D110+dz*F110," ")=IF(B111<=L,hi*Ai/(wC*CpC)*(D111-C111)," =IF(B111<=L,ho*Ao/(wH*CpH)*(D111-C111),"
")
")
=IF(A111<Nsteps,A111+1," ") =IF(B111<L,B111+dz," ") =IF(B112<=L,C111+dz*E111," ")=IF(B112<=L,D111+dz*F111," ")=IF(B112<=L,hi*Ai/(wC*CpC)*(D112-C112)," =IF(B112<=L,ho*Ao/(wH*CpH)*(D112-C112),"
")
")
=IF(A112<Nsteps,A112+1," ") =IF(B112<L,B112+dz," ") =IF(B113<=L,C112+dz*E112," ")=IF(B113<=L,D112+dz*F112," ")=IF(B113<=L,hi*Ai/(wC*CpC)*(D113-C113)," =IF(B113<=L,ho*Ao/(wH*CpH)*(D113-C113),"
")
")
=IF(A113<Nsteps,A113+1," ") =IF(B113<L,B113+dz," ") =IF(B114<=L,C113+dz*E113," ")=IF(B114<=L,D113+dz*F113," ")=IF(B114<=L,hi*Ai/(wC*CpC)*(D114-C114)," =IF(B114<=L,ho*Ao/(wH*CpH)*(D114-C114),"
")
")
=IF(A114<Nsteps,A114+1," ") =IF(B114<L,B114+dz," ") =IF(B115<=L,C114+dz*E114," ")=IF(B115<=L,D114+dz*F114," ")=IF(B115<=L,hi*Ai/(wC*CpC)*(D115-C115)," =IF(B115<=L,ho*Ao/(wH*CpH)*(D115-C115),"
")
")
=IF(A115<Nsteps,A115+1," ") =IF(B115<L,B115+dz," ") =IF(B116<=L,C115+dz*E115," ")=IF(B116<=L,D115+dz*F115," ")=IF(B116<=L,hi*Ai/(wC*CpC)*(D116-C116)," =IF(B116<=L,ho*Ao/(wH*CpH)*(D116-C116),"
")
")
=IF(A116<Nsteps,A116+1," ") =IF(B116<L,B116+dz," ") =IF(B117<=L,C116+dz*E116," ")=IF(B117<=L,D116+dz*F116," ")=IF(B117<=L,hi*Ai/(wC*CpC)*(D117-C117)," =IF(B117<=L,ho*Ao/(wH*CpH)*(D117-C117),"
")
")
=IF(A117<Nsteps,A117+1," ") =IF(B117<L,B117+dz," ") =IF(B118<=L,C117+dz*E117," ")=IF(B118<=L,D117+dz*F117," ")=IF(B118<=L,hi*Ai/(wC*CpC)*(D118-C118)," =IF(B118<=L,ho*Ao/(wH*CpH)*(D118-C118),"
")
")
=IF(A118<Nsteps,A118+1," ") =IF(B118<L,B118+dz," ") =IF(B119<=L,C118+dz*E118," ")=IF(B119<=L,D118+dz*F118," ")=IF(B119<=L,hi*Ai/(wC*CpC)*(D119-C119)," =IF(B119<=L,ho*Ao/(wH*CpH)*(D119-C119),"
")
")
=IF(A119<Nsteps,A119+1," ") =IF(B119<L,B119+dz," ") =IF(B120<=L,C119+dz*E119," ")=IF(B120<=L,D119+dz*F119," ")=IF(B120<=L,hi*Ai/(wC*CpC)*(D120-C120)," =IF(B120<=L,ho*Ao/(wH*CpH)*(D120-C120),"
")
")
=IF(A120<Nsteps,A120+1," ") =IF(B120<L,B120+dz," ") =IF(B121<=L,C120+dz*E120," ")=IF(B121<=L,D120+dz*F120," ")=IF(B121<=L,hi*Ai/(wC*CpC)*(D121-C121)," =IF(B121<=L,ho*Ao/(wH*CpH)*(D121-C121),"
")
")
=IF(A121<Nsteps,A121+1," ") =IF(B121<L,B121+dz," ") =IF(B122<=L,C121+dz*E121," ")=IF(B122<=L,D121+dz*F121," ")=IF(B122<=L,hi*Ai/(wC*CpC)*(D122-C122)," =IF(B122<=L,ho*Ao/(wH*CpH)*(D122-C122),"
")
")
=IF(A122<Nsteps,A122+1," ") =IF(B122<L,B122+dz," ") =IF(B123<=L,C122+dz*E122," ")=IF(B123<=L,D122+dz*F122," ")=IF(B123<=L,hi*Ai/(wC*CpC)*(D123-C123)," =IF(B123<=L,ho*Ao/(wH*CpH)*(D123-C123),"
")
")
=IF(A123<Nsteps,A123+1," ") =IF(B123<L,B123+dz," ") =IF(B124<=L,C123+dz*E123," ")=IF(B124<=L,D123+dz*F123," ")=IF(B124<=L,hi*Ai/(wC*CpC)*(D124-C124)," =IF(B124<=L,ho*Ao/(wH*CpH)*(D124-C124),"
")
")
=IF(A124<Nsteps,A124+1," ") =IF(B124<L,B124+dz," ") =IF(B125<=L,C124+dz*E124," ")=IF(B125<=L,D124+dz*F124," ")=IF(B125<=L,hi*Ai/(wC*CpC)*(D125-C125)," =IF(B125<=L,ho*Ao/(wH*CpH)*(D125-C125),"
")
")
=IF(A125<Nsteps,A125+1," ") =IF(B125<L,B125+dz," ") =IF(B126<=L,C125+dz*E125," ")=IF(B126<=L,D125+dz*F125," ")=IF(B126<=L,hi*Ai/(wC*CpC)*(D126-C126)," =IF(B126<=L,ho*Ao/(wH*CpH)*(D126-C126),"
")
")
=IF(A126<Nsteps,A126+1," ") =IF(B126<L,B126+dz," ") =IF(B127<=L,C126+dz*E126," ")=IF(B127<=L,D126+dz*F126," ")=IF(B127<=L,hi*Ai/(wC*CpC)*(D127-C127)," =IF(B127<=L,ho*Ao/(wH*CpH)*(D127-C127),"
")
")
=IF(A127<Nsteps,A127+1," ") =IF(B127<L,B127+dz," ") =IF(B128<=L,C127+dz*E127," ")=IF(B128<=L,D127+dz*F127," ")=IF(B128<=L,hi*Ai/(wC*CpC)*(D128-C128)," =IF(B128<=L,ho*Ao/(wH*CpH)*(D128-C128),"
")
")
=IF(A128<Nsteps,A128+1," ") =IF(B128<L,B128+dz," ") =IF(B129<=L,C128+dz*E128," ")=IF(B129<=L,D128+dz*F128," ")=IF(B129<=L,hi*Ai/(wC*CpC)*(D129-C129)," =IF(B129<=L,ho*Ao/(wH*CpH)*(D129-C129),"
")
")
=IF(A129<Nsteps,A129+1," ") =IF(B129<L,B129+dz," ") =IF(B130<=L,C129+dz*E129," ")=IF(B130<=L,D129+dz*F129," ")=IF(B130<=L,hi*Ai/(wC*CpC)*(D130-C130)," =IF(B130<=L,ho*Ao/(wH*CpH)*(D130-C130),"
")
")
=IF(A130<Nsteps,A130+1," ") =IF(B130<L,B130+dz," ") =IF(B131<=L,C130+dz*E130," ")=IF(B131<=L,D130+dz*F130," ")=IF(B131<=L,hi*Ai/(wC*CpC)*(D131-C131)," =IF(B131<=L,ho*Ao/(wH*CpH)*(D131-C131),"
")
")
=IF(A131<Nsteps,A131+1," ") =IF(B131<L,B131+dz," ") =IF(B132<=L,C131+dz*E131," ")=IF(B132<=L,D131+dz*F131," ")=IF(B132<=L,hi*Ai/(wC*CpC)*(D132-C132)," =IF(B132<=L,ho*Ao/(wH*CpH)*(D132-C132),"
")
")
=IF(A132<Nsteps,A132+1," ") =IF(B132<L,B132+dz," ") =IF(B133<=L,C132+dz*E132," ")=IF(B133<=L,D132+dz*F132," ")=IF(B133<=L,hi*Ai/(wC*CpC)*(D133-C133)," =IF(B133<=L,ho*Ao/(wH*CpH)*(D133-C133),"
")
")
=IF(A133<Nsteps,A133+1," ") =IF(B133<L,B133+dz," ") =IF(B134<=L,C133+dz*E133," ")=IF(B134<=L,D133+dz*F133," ")=IF(B134<=L,hi*Ai/(wC*CpC)*(D134-C134)," =IF(B134<=L,ho*Ao/(wH*CpH)*(D134-C134),"
")
")