Differential Equations Preview

d(T1)d(t)=W·Cp·(T0T1)+UA·(TsteamT1)M·Cp
d(T2)d(t)=W·Cp·(T1T2)+UA·(TsteamT2)M·Cp
d(T3)d(t)=W·Cp·(T2T3)+UA·(TsteamT3)M·Cp
W=100
Cp=2
T0=20
UA=10
Tsteam=250
M=1000
Boundaries
t(0)=0
T1(0)=20
T2(0)=20
T3(0)=20
t(f)=200

Program validated. Ready to Solve.
# Example 11 - ODE System. Heat Exchange
# Verified Final Values: T1 = 30.9524, T2 = 41.3832 , T3 = 51.3174
# Ref.: Comput. Appl. Eng. Educ. 6: 175, 1998

d(T1)/d(t)=(W*Cp*(T0-T1)+UA*(Tsteam-T1))/(M*Cp)
d(T2)/d(t)=(W*Cp*(T1-T2)+UA*(Tsteam-T2))/(M*Cp)
d(T3)/d(t)=(W*Cp*(T2-T3)+UA*(Tsteam-T3))/(M*Cp)

W=100
Cp=2.0
T0=20
UA=10.
Tsteam=250
M=1000
t(0)=0
T1(0)=20
T2(0)=20
T3(0)=20
t(f)=200

Calculated values

  Variable ▲ Initial Minimal Maximal Final Value
1t00200200
2T1202030.9530.952381
3T2202041.3841.38322
4T3202051.3251.31735
Chart 1
Chart

Data Table (39 records)


#tT1T2T3
10202020
22.37571622.41796722.69325622.714602
34.293166623.97429824.7635324.872118
46.35093425.33026326.83612427.13647
58.55798426.49318728.86099229.483126
610.93404527.47777430.80029931.88762
713.50350628.2993932.62108634.318214
816.29856328.97414434.29668436.738266
919.36362529.51851735.80726239.108425
1022.76395829.9490437.1402841.389442
1126.43283530.26981738.2463543.458527
1230.37435730.50114639.13159645.2695
1334.630130.66375739.82016446.809853
1439.24816530.7746640.33867648.07838
1544.28935230.84770640.7152149.08613
1649.8310430.89388740.9774949.85457
1755.8655130.92134341.14918550.40552
1862.2125130.93644341.25292250.7696
1968.8791830.94446841.31324451.000454
2075.54584530.94845241.3458551.13595
2182.2125130.9504341.36336551.21447
2288.8791830.95141241.37271551.259476
2395.54584530.951941.3776851.285027
24102.2125130.95214341.3803151.299416
25108.8791830.95226341.381751.307457
26115.5458430.95232241.38242351.31192
27122.2125130.95235341.38280551.314384
28128.8791830.95236641.38300351.31574
29135.5458430.95237441.3831151.31648
30142.2125130.95237741.38316351.31688
31148.8791830.9523841.3831951.317097
32155.5458430.9523841.38320551.317215
33162.2125130.95238141.38321351.31728
34168.8791830.95238141.38321751.317314
35175.5458430.95238141.38321751.317333
36182.2125130.95238141.3832251.31734
37188.8791830.95238141.3832251.317345
38195.5458430.95238141.3832251.31735
3920030.95238141.3832251.31735
Raw data presented with 39 records.

Constant Values

ConstantWCpT0UATsteamM
Value100220102501000

Settings and Hints

 Name/SourceValueType
1 Total number of equations 9 Setting
2 Number of differential equations 3 Setting
3 Number of explicit equations 6 Setting
4 Reporting digits 10 Setting
5 Minimal License Free - 2 Setting
6 Elapsed time 0.084 ms Setting
7 Solution method RKF_45 Setting
8 Step size guess. h 6.6666667 Setting
9 Truncation error tolerance. eps 0.000001 Setting
10 Good steps 47 Setting
11 Bad steps 0 Setting
12 Calculated Intermediate data points 50 Setting
13 Version 7.0.73 Setting
14 #@chart_size = 428,300 Default (array2)
15 #@chart_interpolate = 2 Default (integer)
16 #@chart_curve_width = 2 Default (integer)
17 #@report_fix_digits = 8 Default (integer)
18 #@report_show_charts = true Default (boolean)
19 #@report_show_data_points = true Default (boolean)
20 #@report_show_intermediate_data_points = true Default (boolean)
21 #@report_show_source = true Default (boolean)
22 #@report_show_header = false Default (boolean)
23 #@report_show_settings = true Default (boolean)
24 #@chart_all_curves = false Default (boolean)
25 #@chart_curves_stiff_grouping_percent = 50 Default (integer)
26 #@chart_y_curves = Default (array2)
27 #@chart_log_scales = 000 Default (string)

Messages

 SeverityDescription
1Info-Infor- Differential variables (3) ; T1, T2, T3
2Info-Infor- Explicit variables (6) ; W, Cp, T0, UA, Tsteam, M
3Info-Infor- Independent variable = t

Source

# Example 11 - ODE System. Heat Exchange
# Verified Final Values: T1 = 30.9524, T2 = 41.3832 , T3 = 51.3174
# Ref.: Comput. Appl. Eng. Educ. 6: 175, 1998

d(T1)/d(t)=(W*Cp*(T0-T1)+UA*(Tsteam-T1))/(M*Cp)
d(T2)/d(t)=(W*Cp*(T1-T2)+UA*(Tsteam-T2))/(M*Cp)
d(T3)/d(t)=(W*Cp*(T2-T3)+UA*(Tsteam-T3))/(M*Cp)

W=100
Cp=2.0
T0=20
UA=10.
Tsteam=250
M=1000
t(0)=0
T1(0)=20
T2(0)=20
T3(0)=20
t(f)=200