My function will not work to calculate the 'Cumene inlet temperature' section
3 次查看(过去 30 天)
显示 更早的评论
clc
clear
close all
format short g
%% Properties taken from Perry's Chemical Engineering Handbook 8th Edition
T1=414.95;
t1=286.15;t2=379.25;
Ac=61723; Bc=494.81; Cc = 0;%Heat capacity coefficients Cumene
Ab=162940; Bb=-344.94; Cb=0.85562; %Heat capacity coefficients Benzene
B1b=7.5117; B2b=294.68; B3b=-2.794;%viscosity coefficients Benzene
B1c=-24.988; B2c=1807.9; B3c=2.0556;%viscosity coefficients Cumene
C1b=0.23444; C2b=-0.00030572;%Thermal Conductivity coefficients Benzene
C1c=0.1855; C2c=-0.00020895;%Thermal Conductivity coefficients cumene
D1b=1.0259; D2b=0.26666; D3b=562.05; D4b=0.28394; %Density coefficients benzene
D1c=0.58711; D2c=0.25583; D3c=631; D4c=0.28498; %Density coefficients cumene
%% Flow rate calculation
mb=281.4749; mc=281.25;%molar flow rate in kmol/hr
m_b = 281.4749e3; m_c = 281.25e3; %molar flow rate in mol/hr
GFMb=78.11/1000; GFMc=120.19/1000;%GFM of Benzene and Cumene in kg/mol
Mb=(m_b*GFMb)/3600; %flow rate benzene in kg/s
Mc=(m_c*GFMc)/3600;% flow rate cumene in kg/s
%%
%Duty calculation
Cp=@(T)(Ac+(Bc*T));
Hb=integral(Cp,t1,t2);
Qb=(Hb.*Mb)./1000; %Cumene Duty in kW
Qc=Qb.*-1; %Benzene duty in kW
%%Cumene inlet temperature calculation
T=400;
func= @(T) (Mc*Cp)*(T-T1) - (Qb);
T2 = func(T);
0 个评论
采纳的回答
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Thermal Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!