Can somebody run this code for me, It's taking too much time for me
2 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Ajmal R S
am 7 Okt. 2020
Kommentiert: Stephen23
am 9 Okt. 2020
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
0 Kommentare
Akzeptierte Antwort
Ameer Hamza
am 7 Okt. 2020
Use a stiff solve, ode15s()
[t2sol,T2sol] = ode15s(@(t,T2) heat2(t,T2), tspan, T2_0);
1 Kommentar
Stephen23
am 9 Okt. 2020
Simpler without the anonymous function:
[t2sol,T2sol] = ode15s(@heat2, tspan, T2_0);
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Ordinary Differential Equations finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!