Index exceeds the number of array elements. Index must not exceed 1.
%% Design Project Part 2
% 2.1
% Constants:
kapp = 0.009; % m^3/mol*min
L = 1; % m
Vtot = 0.25; % m^3
FT0 = 20; % mol/ min
v0 = 0.08; % m^3/min
P = 100; % atm
rho = 4.2; % kg/m^3
mu = 0.492*10^-5; % kg/m*s
dp = 0.015; % m
phi = 0.75; % ?
km = 4; % 1/min
Fsg0 = 0.1; % mol/min
% Cases:
vsg0 = 0.1; % m^3/min
%vsg0 = 0.001; % m^3/min
v = @(FA, FB, FR, FP) ((v0/FT0)*(FA+FB+FR+FP));
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
CA = @(FA, FB, FR, FP) FA/v(FA, FB, FR, FP);
CB = @(FA, FB, FR, FP) FB/v(FA, FB, FR, FP);
CP = @(FA, FB, FR, FP) FP/v(FA, FB, FR, FP);
CR = @(FA, FB, FR, FP) FR/v(FA, FB, FR, FP);
CPsg = @(FPsg) FPsg/vsg(FPsg);
dFAdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFBdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFRdV = @(FA, FB, FR, FP) -kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
dFpsgdV = @(FPsg, FP) km(CP-CPsg);
FA = 10; % mol/ min
FB = 10; % mol/ min
FR = 0; % mol/ min
FP = 0; % mol/ min
FPsg = 0.1; % mol/min
dV = 0.01;
Vlim = 0.25;
Vspan = 0:dV:Vlim;
CA_values = zeros(size(Vspan));
CB_values = zeros(size(Vspan));
CP_values = zeros(size(Vspan));
CR_values = zeros(size(Vspan));
CPsg_values = zeros(size(Vspan));
FA_values = zeros(size(Vspan));
FB_values = zeros(size(Vspan));
FR_values = zeros(size(Vspan));
FP_values = zeros(size(Vspan));
FPsg_values = zeros(size(Vspan));
for i=1:length(Vspan)
dFA = dFAdV(FA, FB, FR, FP);
dFB = dFBdV(FA, FB, FR, FP);
dFR = dFRdV(FA, FB, FR, FP);
dFP = dFPdV(FPsg, FA, FB, FP, FR);
dFPsg = dFPsgdV(Fpsg, FP);
CA_values(i) = CA(FA, FB, FR, FP);
CB_values(i) = CB(FA, FB, FR, FP);
CP_values(i) = CP(FA, FB, FR, FP);
CR_values(i) = CR(FA, FB, FR, FP);
CPsg_values(i) = CPsg(Fpsg);
dFA = FA+dFA*dV;
dFB = FB+dFB*dV;
dFR = FR+dFR*dV;
dFP = FP+dFP*dV;
dFPsg = FPsg+dFPsg*dV;
FA_values(i) = FA;
FB_values(i) = FB;
FR_values(i) = FR;
FP_values(i) = FP;
FPsg_values(i) = FPsg;
end
figure()
hold on;
plot(Vspan, FA_values,’LineWidth’, 2);
plot(Vspan, FB_values,’LineWidth’, 2);
plot(Vspan, FP_values,’LineWidth’, 2);
plot(Vspan, FR_values,’LineWidth’, 2);
plot(Vspan, FPsg_values,’LineWidth’, 2);
xlabel(‘PBR volume V [=] m^3’);
ylabel(‘f(V) [=] mol/min’);
set(gca, ‘LineWidth’, 2, ‘FontSize’, 20);
hold off;
Error in Design_Project1>@(FPsg)vsg0((Fsg0+FPsg)/Fsg0) (line 27)
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
Error in Design_Project1>@(FPsg)FPsg/vsg(FPsg) (line 33)
CPsg = @(FPsg) FPsg/vsg(FPsg);
Error in Design_Project1>@(FPsg,FA,FB,FP,FR)(-kapp*CA(FA,FB,FR,FP)*CB(FA,FB,FR,FP))-km(CP(FA,FB,FR,FP)-CPsg(FPsg)) (line 38)
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
Error in Design_Project1 (line 68)
dFP = dFPdV(FPsg, FA, FB, FP, FR);%% Design Project Part 2
% 2.1
% Constants:
kapp = 0.009; % m^3/mol*min
L = 1; % m
Vtot = 0.25; % m^3
FT0 = 20; % mol/ min
v0 = 0.08; % m^3/min
P = 100; % atm
rho = 4.2; % kg/m^3
mu = 0.492*10^-5; % kg/m*s
dp = 0.015; % m
phi = 0.75; % ?
km = 4; % 1/min
Fsg0 = 0.1; % mol/min
% Cases:
vsg0 = 0.1; % m^3/min
%vsg0 = 0.001; % m^3/min
v = @(FA, FB, FR, FP) ((v0/FT0)*(FA+FB+FR+FP));
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
CA = @(FA, FB, FR, FP) FA/v(FA, FB, FR, FP);
CB = @(FA, FB, FR, FP) FB/v(FA, FB, FR, FP);
CP = @(FA, FB, FR, FP) FP/v(FA, FB, FR, FP);
CR = @(FA, FB, FR, FP) FR/v(FA, FB, FR, FP);
CPsg = @(FPsg) FPsg/vsg(FPsg);
dFAdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFBdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFRdV = @(FA, FB, FR, FP) -kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
dFpsgdV = @(FPsg, FP) km(CP-CPsg);
FA = 10; % mol/ min
FB = 10; % mol/ min
FR = 0; % mol/ min
FP = 0; % mol/ min
FPsg = 0.1; % mol/min
dV = 0.01;
Vlim = 0.25;
Vspan = 0:dV:Vlim;
CA_values = zeros(size(Vspan));
CB_values = zeros(size(Vspan));
CP_values = zeros(size(Vspan));
CR_values = zeros(size(Vspan));
CPsg_values = zeros(size(Vspan));
FA_values = zeros(size(Vspan));
FB_values = zeros(size(Vspan));
FR_values = zeros(size(Vspan));
FP_values = zeros(size(Vspan));
FPsg_values = zeros(size(Vspan));
for i=1:length(Vspan)
dFA = dFAdV(FA, FB, FR, FP);
dFB = dFBdV(FA, FB, FR, FP);
dFR = dFRdV(FA, FB, FR, FP);
dFP = dFPdV(FPsg, FA, FB, FP, FR);
dFPsg = dFPsgdV(Fpsg, FP);
CA_values(i) = CA(FA, FB, FR, FP);
CB_values(i) = CB(FA, FB, FR, FP);
CP_values(i) = CP(FA, FB, FR, FP);
CR_values(i) = CR(FA, FB, FR, FP);
CPsg_values(i) = CPsg(Fpsg);
dFA = FA+dFA*dV;
dFB = FB+dFB*dV;
dFR = FR+dFR*dV;
dFP = FP+dFP*dV;
dFPsg = FPsg+dFPsg*dV;
FA_values(i) = FA;
FB_values(i) = FB;
FR_values(i) = FR;
FP_values(i) = FP;
FPsg_values(i) = FPsg;
end
figure()
hold on;
plot(Vspan, FA_values,’LineWidth’, 2);
plot(Vspan, FB_values,’LineWidth’, 2);
plot(Vspan, FP_values,’LineWidth’, 2);
plot(Vspan, FR_values,’LineWidth’, 2);
plot(Vspan, FPsg_values,’LineWidth’, 2);
xlabel(‘PBR volume V [=] m^3’);
ylabel(‘f(V) [=] mol/min’);
set(gca, ‘LineWidth’, 2, ‘FontSize’, 20);
hold off;
Error in Design_Project1>@(FPsg)vsg0((Fsg0+FPsg)/Fsg0) (line 27)
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
Error in Design_Project1>@(FPsg)FPsg/vsg(FPsg) (line 33)
CPsg = @(FPsg) FPsg/vsg(FPsg);
Error in Design_Project1>@(FPsg,FA,FB,FP,FR)(-kapp*CA(FA,FB,FR,FP)*CB(FA,FB,FR,FP))-km(CP(FA,FB,FR,FP)-CPsg(FPsg)) (line 38)
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
Error in Design_Project1 (line 68)
dFP = dFPdV(FPsg, FA, FB, FP, FR); %% Design Project Part 2
% 2.1
% Constants:
kapp = 0.009; % m^3/mol*min
L = 1; % m
Vtot = 0.25; % m^3
FT0 = 20; % mol/ min
v0 = 0.08; % m^3/min
P = 100; % atm
rho = 4.2; % kg/m^3
mu = 0.492*10^-5; % kg/m*s
dp = 0.015; % m
phi = 0.75; % ?
km = 4; % 1/min
Fsg0 = 0.1; % mol/min
% Cases:
vsg0 = 0.1; % m^3/min
%vsg0 = 0.001; % m^3/min
v = @(FA, FB, FR, FP) ((v0/FT0)*(FA+FB+FR+FP));
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
CA = @(FA, FB, FR, FP) FA/v(FA, FB, FR, FP);
CB = @(FA, FB, FR, FP) FB/v(FA, FB, FR, FP);
CP = @(FA, FB, FR, FP) FP/v(FA, FB, FR, FP);
CR = @(FA, FB, FR, FP) FR/v(FA, FB, FR, FP);
CPsg = @(FPsg) FPsg/vsg(FPsg);
dFAdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFBdV = @(FA, FB, FR, FP) kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFRdV = @(FA, FB, FR, FP) -kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP);
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
dFpsgdV = @(FPsg, FP) km(CP-CPsg);
FA = 10; % mol/ min
FB = 10; % mol/ min
FR = 0; % mol/ min
FP = 0; % mol/ min
FPsg = 0.1; % mol/min
dV = 0.01;
Vlim = 0.25;
Vspan = 0:dV:Vlim;
CA_values = zeros(size(Vspan));
CB_values = zeros(size(Vspan));
CP_values = zeros(size(Vspan));
CR_values = zeros(size(Vspan));
CPsg_values = zeros(size(Vspan));
FA_values = zeros(size(Vspan));
FB_values = zeros(size(Vspan));
FR_values = zeros(size(Vspan));
FP_values = zeros(size(Vspan));
FPsg_values = zeros(size(Vspan));
for i=1:length(Vspan)
dFA = dFAdV(FA, FB, FR, FP);
dFB = dFBdV(FA, FB, FR, FP);
dFR = dFRdV(FA, FB, FR, FP);
dFP = dFPdV(FPsg, FA, FB, FP, FR);
dFPsg = dFPsgdV(Fpsg, FP);
CA_values(i) = CA(FA, FB, FR, FP);
CB_values(i) = CB(FA, FB, FR, FP);
CP_values(i) = CP(FA, FB, FR, FP);
CR_values(i) = CR(FA, FB, FR, FP);
CPsg_values(i) = CPsg(Fpsg);
dFA = FA+dFA*dV;
dFB = FB+dFB*dV;
dFR = FR+dFR*dV;
dFP = FP+dFP*dV;
dFPsg = FPsg+dFPsg*dV;
FA_values(i) = FA;
FB_values(i) = FB;
FR_values(i) = FR;
FP_values(i) = FP;
FPsg_values(i) = FPsg;
end
figure()
hold on;
plot(Vspan, FA_values,’LineWidth’, 2);
plot(Vspan, FB_values,’LineWidth’, 2);
plot(Vspan, FP_values,’LineWidth’, 2);
plot(Vspan, FR_values,’LineWidth’, 2);
plot(Vspan, FPsg_values,’LineWidth’, 2);
xlabel(‘PBR volume V [=] m^3’);
ylabel(‘f(V) [=] mol/min’);
set(gca, ‘LineWidth’, 2, ‘FontSize’, 20);
hold off;
Error in Design_Project1>@(FPsg)vsg0((Fsg0+FPsg)/Fsg0) (line 27)
vsg = @(FPsg) vsg0((Fsg0+FPsg)/Fsg0);
Error in Design_Project1>@(FPsg)FPsg/vsg(FPsg) (line 33)
CPsg = @(FPsg) FPsg/vsg(FPsg);
Error in Design_Project1>@(FPsg,FA,FB,FP,FR)(-kapp*CA(FA,FB,FR,FP)*CB(FA,FB,FR,FP))-km(CP(FA,FB,FR,FP)-CPsg(FPsg)) (line 38)
dFPdV = @(FPsg, FA, FB, FP, FR) (-kapp*CA(FA, FB, FR, FP)*CB(FA, FB, FR, FP))-km(CP(FA, FB, FR, FP)-CPsg(FPsg));
Error in Design_Project1 (line 68)
dFP = dFPdV(FPsg, FA, FB, FP, FR); ode MATLAB Answers — New Questions