Can somebody run this code for me, It's taking too much time for me

Question

Ajmal R S 2020-10-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/607156-can-somebody-run-this-code-for-me-it-s-taking-too-much-time-for-me

评论： Stephen23 2020-10-9

clc;
clear all;
T2_0 = 320.1;
tspan = [0 50];
[t2sol,T2sol] = ode45(@(t,T2) heat2(t,T2), tspan, T2_0);
plot(t2sol,T2sol)
grid on
xlabel('t')
ylabel('T2')
title('Solutions of Heat Balance on Tank-2')

This is a code that I modelled to study about Stirred Tank Heater Performance

function dT2 = heat2(t,T2) 
% Inputs
A2  = 7.854e-05;  %Cross-Sectional Area Tank 2
Tc  = 303;        %Cooling Water Temp
T1  = 325.4;      %Tank 1 Temp
Ta  = 298;        %Atmospheric Temp
h2  = 0.41;       %Level h2
R2  = 0.05;       %Radius Tank 2
rho_w = 1000;     %Density in kg/m3
Cp  = 4182;       %Cp value in J/KgK
U   =  235.1043;   %Heat Transfer Coeff
u1  = 0.45;       %Flow F1
F1_A =  8.0368e-11;
F1_B =  456.85e-11;
F1_C =  42379e-11;
F1 = F1_A*(u1^3)-F1_B*(u1^2)+F1_C*(u1);
u2  = 0.45;       %Flow F2
F2_A =  1.2943e-11;
F2_B =  190.64e-11;
F2_C =  8796.8e-11;
F2_D =  196620e-11;
F2 = -F2_A*(u2^4)+F2_B*(u2^3)-F2_C*(u2^2)+F2_D*(u2);
u3  = 0.45;       %Flow F3
Fr_A =  (1/1800)*1e-03;
Fr = Fr_A*u3;
u5  = 0.45;       %Heat Input Q2
Q2_A = 14.4;
Q2_B = 0.96;
Q2_C = 0.008;
Q2 = Q2_A*(u5)-Q2_B*(u5^2)-Q2_C*(u5^3);
dT2 = (F1*(T1-T2)+F2*(Tc-T2)-Fr*(T2-T1)+(Q2-(2*pi*R2*h2*U*(T2-Ta)))/rho_w*Cp)/(A2*h2);
end

And this is the function "heat2" shown in code. Can Anybdody what output this shows. Since it's taking a really long time for and my laptop is heating up and switching off before it completes

Thank you