1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143
|
function T_max = hassane_25_essai(Par)
clf
global k0_perh k_perh Ea_perh H_perh Keq_perh_30 Keq_perh Rperh k0_decomp_1 k_decomp_1 Ea_decomp_1 H_decomp_1 R_decomp_1 k0_decomp_2 k_decomp_2 Ea_decomp_2 H_decomp_2 R_decomp_2 k0_Ep_1 k_Ep_1 Ea_Ep_1 R_Ep_1 k0_Ep_2 Ea_Ep_2 k_Ep_2 R_Ep_2 k0_Ep_3 k_Ep_3 R_Ep_3 Ea_Ep_3 H_Ep H_RO k0_RO_W Ea_RO_W k_RO_W R_RO_W k0_RO_FA Ea_RO_FA k_RO_FA R_RO_FA k0_RO_PFA Ea_RO_PFA k_RO_PFA R_RO_PFA Temp_ref mj Cpj Qj UA R Cpr mr rho Tj0 Tj_degree Temp_degree qrx qexh Temp Vaq0 Vorg0 Vtot0 Taj FA_aj CpFA FA PFA W HP OA LA Int Ep RO Vol F Vaq tadd KFA K_FA W_initial CHP_initial
% LES SIX VARIABLES DENTREE
F= Par(1);
tadd=Par(2);
UA=Par(3);
Taj=Par(4);
Tj0=Par(5);
CHP_initial=Par(6);
k0_perh=0.15;
Ea_perh=150000;
H_perh=-5580;
Keq_perh_30=0.96;
k0_decomp_1=0.001;
Ea_decomp_1=20000;
H_decomp_1=-359000;
k0_decomp_2=0.0009;
Ea_decomp_2=20200;
H_decomp_2=-163000;
k0_Ep_1=0.00576;
Ea_Ep_1=61600;
k0_Ep_2=0.00437;
Ea_Ep_2=15000;
k0_Ep_3=0.004;
Ea_Ep_3=15000;
H_Ep=-230000;
H_RO=-90000;
k0_RO_W=0.000237;
Ea_RO_W=150000;
k0_RO_FA=0.00339;
Ea_RO_FA=42100;
k0_RO_PFA=0.0592;
Ea_RO_PFA=70000;
Temp_ref=66.85+273.15;
R=8.314;
mj=0.3;
Qj=15/60;
Cpj=4000;
Cpr= 2000;
mr= 0.4;
rho=1;
Vaq0=0.14;
Vorg0=0.12;
Vtot0=Vaq0+Vorg0;
FA_aj=24;
CpFA=20;
KFA=2.9E-6;
K_FA=9;
W_initial=((1-(Vaq0*CHP_initial)*(34/(Vaq0*1000)))*1000)/18;
tspan=[0:10:10000];
y0=[0 CHP_initial 0 W_initial 0.5 1.70 0.00 0.00 0.00 Tj0 Vtot0];
[t, y]=ode23s(@semibatch,tspan,y0);
T_max = max(y(:,10)); %la variable sortie recherchée.
FA=y(:,1);
HP=y(:,2);
PFA=y(:,3);
W=y(:,4);
OA=y(:,5);
LA=y(:,6);
Int=y(:,7);
Ep=y(:,8);
RO=y(:,9);
Temp=y(:,10);
Vol=y(:,11);
Vaq=Vol-Vorg0;
Tj_degree=Tj0-273.15;
Temp_degree=Temp-273.15;
k_perh=k0_perh*exp(-(Ea_perh/R)*((1./Temp)-(1/Temp_ref)));
k_decomp_1=k0_decomp_1*exp(-(Ea_decomp_1/R)*((1./Temp)-(1/Temp_ref)));
k_decomp_2=k0_decomp_2*exp(-(Ea_decomp_2/R)*((1./Temp)-(1/Temp_ref)));
k_Ep_1=k0_Ep_1*exp(-(Ea_Ep_1/R)*((1./Temp)-(1/Temp_ref)));
k_Ep_2=k0_Ep_2*exp(-(Ea_Ep_2/R)*((1./Temp)-(1/Temp_ref))); k_Ep_3=k0_Ep_3*exp(-(Ea_Ep_3/R)*((1./Temp)-(1/Temp_ref))); k_RO_W=k0_RO_W*exp(-(Ea_RO_W/R)*((1./Temp)-(1/Temp_ref))); k_RO_FA=k0_RO_FA*exp(-(Ea_RO_FA/R)*((1./Temp)-(1/Temp_ref))); k_RO_PFA=k0_RO_PFA*exp(-(Ea_RO_PFA/R)*((1./Temp)-(1/Temp_ref))); Keq_perh=Keq_perh_30*exp((-H_perh/R)*((1./Temp)-(1/303.15)));
Rperh=k_perh*(KFA^0.5).*((FA./W).^0.5).*(FA.*HP-(1./Keq_perh).*PFA.*W);
R_decomp_1=k_decomp_1.*PFA;
R_decomp_2=k_decomp_2.*PFA;
R_Ep_1=k_Ep_1.*PFA.*OA;
R_Ep_2=k_Ep_2.*PFA.*LA;
R_Ep_3=k_Ep_3.*PFA.*Int;
R_RO_W=k_RO_W*(KFA^0.5).*Ep.*((FA.*W).^0.5);
R_RO_FA=k_RO_FA*(KFA^0.5).*Ep.*FA.*((FA./W).^0.5);
R_RO_PFA=k_RO_PFA*(KFA^0.5).*Ep.*PFA.*((FA./W).^0.5);
qrx=-Vaq.*(Rperh*H_perh+R_decomp_1*H_decomp_1+R_decomp_2*H_decomp_2)-Vorg0*((R_Ep_1+R_Ep_2+R_Ep_3)*H_Ep+(R_RO_W+R_RO_FA+R_RO_PFA)*H_RO);
qexh=UA*(Tj0-Temp);
function dydt=semibatch(t,y)
global k0_perh k_perh Ea_perh H_perh Keq_perh_30 Keq_perh Rperh k0_decomp_1 k_decomp_1 Ea_decomp_1 H_decomp_1 R_decomp_1 k0_decomp_2 k_decomp_2 Ea_decomp_2 H_decomp_2 R_decomp_2 k0_Ep_1 k_Ep_1 Ea_Ep_1 R_Ep_1 k0_Ep_2 Ea_Ep_2 k_Ep_2 R_Ep_2 k0_Ep_3 k_Ep_3 R_Ep_3 Ea_Ep_3 H_Ep H_RO k0_RO_W Ea_RO_W k_RO_W R_RO_W k0_RO_FA Ea_RO_FA k_RO_FA R_RO_FA k0_RO_PFA Ea_RO_PFA k_RO_PFA R_RO_PFA Temp_ref mj Cpj Qj UA R Cpr mr rho Tj0 Tj_degree Temp_degree qrx qexh Temp Vaq0 Vorg0 Vtot0 Taj FA_aj CpFA FA PFA W HP OA LA Int Ep RO Vol F Vaq tadd KFA Vorg Vtot a K_FA mcp
k_perh=k0_perh*exp(-(Ea_perh/R)*((1/y(10))-(1/Temp_ref)));
k_decomp_1=k0_decomp_1*exp(-(Ea_decomp_1/R)*((1/y(10))-(1/Temp_ref))); k_decomp_2=k0_decomp_2*exp(-(Ea_decomp_2/R)*((1/y(10))-(1/Temp_ref)));
k_Ep_1=k0_Ep_1*exp(-(Ea_Ep_1/R)*((1/y(10))-(1/Temp_ref)));
k_Ep_2=k0_Ep_2*exp(-(Ea_Ep_2/R)*((1/y(10))-(1/Temp_ref))); k_Ep_3=k0_Ep_3*exp(-(Ea_Ep_3/R)*((1/y(10))-(1/Temp_ref))); k_RO_W=k0_RO_W*exp(-(Ea_RO_W/R)*((1/y(10))-(1/Temp_ref))); k_RO_FA=k0_RO_FA*exp(-(Ea_RO_FA/R)*((1/y(10))-(1/Temp_ref))); k_RO_PFA=k0_RO_PFA*exp(-(Ea_RO_PFA/R)*((1/y(10))-(1/Temp_ref))); Keq_perh=Keq_perh_30*exp((-H_perh/R)*((1/y(10))-(1/303.15)));
Rperh=k_perh*(KFA^0.5)*((y(1)/y(4))^0.5)*(y(1)*y(2)-(1/Keq_perh)*y(3)*y(4));
R_decomp_1=k_decomp_1*y(3);
R_decomp_2=k_decomp_2*y(3);
R_Ep_1=k_Ep_1*y(3)*y(5);
R_Ep_2=k_Ep_2*y(3)*y(6);
R_Ep_3=k_Ep_3*y(3)*y(7);
R_RO_W=k_RO_W*(KFA^0.5)*y(8)*((y(1)*y(4))^0.5);
R_RO_FA=k_RO_FA*(KFA^0.5)*y(8)*y(1)*((y(1)/y(4))^0.5);
R_RO_PFA=k_RO_PFA*(KFA^0.5)*y(8)*y(3)*((y(1)/y(4))^0.5);
if t > tadd
F=0;
end
Vaq=y(11)-Vorg0;
a=Vaq/y(11);
mcp=(Vaq*1.00+Vorg0*0.93)*Cpr;
dydt=[(1/(1+((1-a)/(a*K_FA))))*((F/Vaq)*(FA_aj-y(1))-Rperh+R_decomp_2+(1-a)*(R_Ep_1+R_Ep_2+R_Ep_3-R_RO_FA)/a);...
-Rperh-(F*y(2)/Vaq); ((-F*y(3)/Vaq)+Rperh-R_decomp_2-R_decomp_1-(1-a)*(R_Ep_1+R_Ep_2+R_Ep_3+R_RO_PFA)/a);...
Rperh+R_decomp_1-(F*y(4)/Vaq)-(1-a)*R_RO_W/a; -R_Ep_1; -R_Ep_2; R_Ep_2-R_Ep_3; (R_Ep_1+R_Ep_2+R_Ep_3)-R_RO_W-R_RO_FA-R_RO_PFA;...
R_RO_W+R_RO_FA+R_RO_PFA; (1/mcp)*((-Vaq*(Rperh*H_perh+R_decomp_1*H_decomp_1+R_decomp_2*H_decomp_2)-Vorg0*((R_Ep_1+R_Ep_2+R_Ep_3)*H_Ep+(R_RO_W+R_RO_FA+R_RO_PFA)*H_RO))+UA*(Tj0-y(10))+FA_aj*F*CpFA*(Taj-y(10))); F]; |
Partager