Rapport :Programmes des
RAC
Module :calcul des reacteurs
Realise par : Encadre par :
Asmaa Sadik [Link]
PROFIL DE CONCENTRATION CA EN FONCTION
DU TEMPS DANS UN REACTEUR PARFAITEMENT
AGITE A MARCHE CONTINUE :
v=1
F0=1
Ca0=1
k=0.01;
t=0:0.01:100
for i=1:length(t)
Ca1(i)=Ca0/(1+(k*v/F0))*(1-exp(-
((F0/v)+k)*t(i)));
end
c=0.75*Ca1 (length(t));
temps=(-1/((F0/v)+k)*log(1-
(1+(k*v/F0))*c/Ca0));
to=5*temps;
plot(t,Ca1)
title('figure de Ca1 en fonction de temps' )
xlabel('temps(s)')
ylabel('concentration(mol/l)')
deux réacteurs parfaitement agités à
marche continue en série :
k=0.01
v=5
f0=1
ca0=0.8
t=0:0.01:160;
for i=1:length(t)
ca1(i)=ca0*(1-exp(-
(k+(f0/v))*t(i)))/(1+(k*v/f0));
end
c=0.75*ca1(length(t))
temps=(-1/(k+(f0/v)))*log(1-
(1+(k*v/f0))*c/ca0);
to=5*temps
m=floor(to)/0.01
for i=m+1:length(t)
ca2(i)=ca0*(1-((f0/v)+k)*exp(-
(k+(f0/v)*t(i)*t(i-m)-exp(-(k+(f0/v))*t(i-
m))/(1+(k*v/f0)))))^2
end
plot(t,ca1,t,ca2);
xlabel('temps(s)')
ylabel('concentration(mol/l)');
CALCUL D’UNE CASCADE DEREACTEUR
PARFAITEMENT AGITES :
Methode de Eldridge et Piret :
k1=4.53*10^-2;
teta=12
ca0=0.736;
can=0.278;
x=0:0.01:1;
a=input('donner l''ordre de la réaction=')
x=0:0.1:1;
if a==1
for j=1:length(x)
c(j)=x(j)+(k1*x(j)^a)*teta;
end
yt=x;
plot(x,c,'r',x,yt,'g')
axis([0,1,0,2])
line([ca0 ca0], [0 ca0])
n=0;
while ca0>=can
y=ca0/(teta*k1+1);
line([y y], [y ca0])
line([y ca0],[ca0 ca0])
n=n+1;
fprintf(' Ca%d = %f \n' , n,y)
ca0=y;
end
line([y y],[0 y])
end
if a==2
for j=1:length(x)
c(j)=x(j)+(k1*x(j)^a)*teta;
end
yt=x;
plot(x,c,'r',x,yt,'g')
axis([0,1,0,2])
line([ca0 ca0], [0 ca0])
n=0;
while ca0>=can
y=((-
1/(teta*k1))+sqrt((1/(teta*k1))^2+(4*ca0/(teta*k1))))/2;
line([y y], [y ca0])
line([y ca0],[ca0 ca0])
n=n+1;
fprintf(' Ca%d = %f \n' , n,y)
ca0=y;
end
line([y y],[0 y])
end
fprintf('le nombre d''étages est N= %d \n' , n)
title('détermination de nombre d''étages par la méthode
d''ELDRIGE et PIRET ')
xlabel('concentration Can mol/l')
ylabel('concentration Can-1 mol/l')
la méthode de Jones :
k1=4.53*10^-2;
teta=12
ca0=0.736;
can=0.278;
x=0:0.01:1;
a=input('donner a=');
if a==1
for j=1:length(x)
v(j)=k1*x(j)^a;
end
plot(x,v,'r')
axis([0,1,0,35*10^-3])
n=0;
while ca0>=can
y=ca0/(teta*(k1+1/teta));
k=k1*y;
line([ca0 y],[0 k])
line([y y], [0 k])
n=n+1;
fprintf(' Ca%d = %f \n' , n,y)
ca0=y;
end
end
if a==2
for j=1:length(x)
v(j)=k1*x(j)^a;
end
plot(x,v,'r')
axis([0,1,0,35*10^-3])
n=0;
while ca0>=can
y=((-
1/(teta*k1))+sqrt((1/(teta*k1))^2+(4*ca0/(teta*k1))))/2;
k=k1*y^2;
line([y y], [0 k])
line([ca0 y],[0 k])
n=n+1;
fprintf(' Ca%d = %f \n' , n,y)
ca0=y;
end
end
fprintf('le nombre d''étages est N= %d \n' , n)
title('détermination de nombre d''étages par la méthode
de jones')
xlabel('concentration mol/l')
ylabel('Vitesse de la réaction ')