programa ar

clc;clear all; clear memory;%[t,y]=ode45('lodoactivadoQr1',[ti tf],[DQO SSV XR O2]);[t,y]=ode45('lodoactivadoQr1',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-b') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-b') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-b') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-b') %Oxigeno disueltohold onformat short g[t(end),y(end,1),y(end,2),y(end,3),y(end,4)]%[t,y(:,1),y(:,2),y(:,3),y(:,4)]% axis([-1 20 0 7.5]) [t,y]=ode45('lodoactivadoQr2',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-r') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-r') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-r') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-r') %Oxigeno disueltohold on[t(end),y(end,1),y(end,2),y(end,3),y(end,4)] [t,y]=ode45('lodoactivadoQr3',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-g') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-g') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-g') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-g') %Oxigeno disueltohold on[t(end),y(end,1),y(end,2),y(end,3),y(end,4)] [t,y]=ode45('lodoactivadoQr4',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-k') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-k') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-k') %XR del sedimentadorhold on

figure (4)plot(t,y(:,4),'-k') %Oxigeno disueltohold on[t(end),y(end,1),y(end,2),y(end,3),y(end,4)]

figure (1); title({'Efecto: Y=0.65';'Variación del caudal de reciclo en degradación de materia organica'})ylabel('Demanda Qimica de Oxígeno - DQO - mg/L'); xlabel('tiempo - d');axis([0 10 0 500]); h = legend('Qr=20','Qr=40','Qr=60','Qr=80',1);figure (2); title({'Efecto: Y=0.65';'Variación del caudal de reciclo en degradación de materia organica'})ylabel('Sólidos Suspendidos Volátiles - SSV - mg/L'); xlabel('tiempo - d');axis([0 15 300 1000]); h = legend('Qr=20','Qr=40','Qr=60','Qr=80',1);figure (3); title({'Efecto: Y=0.65';'Variación del caudal de reciclo en degradación de materia organica'})ylabel('Sólidos Suspendidos Volátiles-Desimentador - Xr - mg/L'); xlabel('tiempo - d');axis([0 15 1000 8000]); h = legend('Qr=20','Qr=40','Qr=60','Qr=80',1);figure (4); title({'Efecto: Y=0.65';'Variación del caudal de reciclo en degradación de materia organica'})ylabel('Concentración de Oxigeno Disuelto - mg/L'); xlabel('tiempo - d');axis([0 8 4.2 8.2]); h = legend('Qr=20','Qr=40','Qr=60','Qr=80',1);

function dy=lodoactivadoQr1(t,y)umax=1.97; %(d-1)b=0.259; %(d-1)Y=0.65; %Coef.rendimiento (mg SSV producido/mg DQO consumidoKs=137.3; %(mg/L)- Constante de afinidadYo=0.915; %Yo=a/YKd=0.0601; %(d-1)- Coeficiente de muerteS=210; %(mg/L)- Concentración del aflueunte a procesarCO2=0.3; %(mg/L)- Cocnentración de oxígeno en el alfuenteCsr=7.02; %(Concent.de saturacion de oxig.en el agua residual)kla=420; %(d-1)- Coeficiente de transferencia de oxígenoQ=760; %(m3/d)- caudal de AR a tratarQR=20; %(m3/d)- Caudal de reflujoQw=35; %(m3/d)- caudal de desecho del sedimentadorV=500; %(m3)- volumen de proceso del reactorVs=20; %(m3)-volumen del sedimentadordy=zeros(4,1);dy(1)=Q*S/V - (Q+QR)*y(1)/V - (umax/Y)*y(1)*y(2)/(Ks+y(1));dy(2)=QR*y(3)/V - (Q+QR)*y(2)/V + umax*y(1)*y(2)/(Ks+y(1)) - Kd*y(2);dy(3)=(Q+QR)*y(2)/Vs - (QR+Qw)*y(3)/Vs;dy(4)=Q*CO2/V - (Q+QR)*y(4)/V - (umax/Yo)*y(1)*y(2)/(Ks+y(1))-b*y(1)+kla*(Csr-y(4));


function dy=lodoactivadoQr3(t,y)umax=1.97; %(d-1)b=0.259; %(d-1)

Y=0.65; %Coef.rendimiento (mg SSV producido/mg DQO consumidoKs=137.3; %(mg/L)- Constante de afinidadYo=0.915; %Yo=a/YKd=0.0601; %(d-1)- Coeficiente de muerteS=210; %(mg/L)- Concentración del aflueunte a procesarCO2=0.3; %(mg/L)- Cocnentración de oxígeno en el alfuenteCsr=7.02; %(Concent.de saturacion de oxig.en el agua residual)kla=420; %(d-1)- Coeficiente de transferencia de oxígenoQ=760; %(m3/d)- caudal de AR a tratarQR=60; %(m3/d)- Caudal de reflujoQw=35; %(m3/d)- caudal de desecho del sedimentadorV=500; %(m3)- volumen de proceso del reactorVs=20; %(m3)-volumen del sedimentadordy=zeros(4,1);dy(1)=Q*S/V - (Q+QR)*y(1)/V - (umax/Y)*y(1)*y(2)/(Ks+y(1));dy(2)=QR*y(3)/V - (Q+QR)*y(2)/V + umax*y(1)*y(2)/(Ks+y(1)) - Kd*y(2);dy(3)=(Q+QR)*y(2)/Vs - (QR+Qw)*y(3)/Vs;dy(4)=Q*CO2/V - (Q+QR)*y(4)/V - (umax/Yo)*y(1)*y(2)/(Ks+y(1))-b*y(1)+kla*(Csr-y(4));


clc;clear all; clear memory;%[t,y]=ode45('lodoactivadoQw1',[ti tf],[DQO SSV XR O2]);[t,y]=ode45('lodoactivadoQw1',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-b') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-b') %crecimiento SSVhold onfigure (3)

plot(t,y(:,3),'-b') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-b') %Oxigeno disueltohold onformat short g[t(end),y(end,1),y(end,2),y(end,3),y(end,4)]%[t,y(:,1),y(:,2),y(:,3),y(:,4)]% axis([-1 20 0 7.5]) [t,y]=ode45('lodoactivadoQw2',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-r') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-r') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-r') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-r') %Oxigeno disueltohold onformat short g[t(end),y(end,1),y(end,2),y(end,3),y(end,4)] [t,y]=ode45('lodoactivadoQw3',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-g') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-g') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-g') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-g') %Oxigeno disueltohold onformat short g[t(end),y(end,1),y(end,2),y(end,3),y(end,4)] [t,y]=ode45('lodoactivadoQw4',[0 15],[210 500 300 2]);hold onfigure (1)plot(t,y(:,1),'-k') %disminucion de DQOhold onfigure (2)plot(t,y(:,2),'-k') %crecimiento SSVhold onfigure (3)plot(t,y(:,3),'-k') %XR del sedimentadorhold onfigure (4)plot(t,y(:,4),'-k') %Oxigeno disueltohold onformat short g[t(end),y(end,1),y(end,2),y(end,3),y(end,4)] figure (1); title({'Efecto: Y=0.35';'Variación del caudal de deshecho en degradación de materia organica'})ylabel('Demanda Qimica de Oxígeno - DQO - mg/L'); xlabel('tiempo - d');axis([0 10 0 700]); h = legend('Qw=15','Qw=25','Qw=35','Qw=45',1);figure (2); title({'Efecto: Y=0.35';'Variación del caudal de deshecho en degradación de materia organica'})

ylabel('Sólidos Suspendidos Volátiles - SSV - mg/L'); xlabel('tiempo - d');axis([0 15 300 1500]); h = legend('Qw=15','Qw=25','Qw=35','Qw=45',1);figure (3); title({'Efecto: Y=0.35';'Variación del caudal de deshecho en degradación de materia organica'})ylabel('Sólidos Suspendidos Volátiles-Desimentador - Xr - mg/L'); xlabel('tiempo - d');axis([0 15 1000 18000]); h = legend('Qw=15','Qw=25','Qw=35','Qw=45',1);figure (4); title({'Efecto: Y=0.35';'Variación del caudal de deshecho en degradación de materia organica'})ylabel('Concentración de Oxigeno Disuelto - mg/L'); xlabel('tiempo - d');axis([0 8 4.2 6.7]); h = legend('Qw=15','Qw=25','Qw=35','Qw=45',1);

function dy=lodoactivadoQw1(t,y)umax=1.97; %(d-1)b=0.259; %(d-1)Y=0.65; %Coef.rendimiento (mg SSV producido/mg DQO consumidoKs=137.3; %(mg/L)- Constante de afinidadYo=0.915; %Yo=a/YKd=0.0601; %(d-1)- Coeficiente de muerteS=668; %(mg/L)- Concentración del aflueunte a procesarCO2=0.3; %(mg/L)- Cocnentración de oxígeno en el alfuenteCsr=7.02; %(Concent.de saturacion de oxig.en el agua residual)kla=420; %(d-1)- Coeficiente de transferencia de oxígenoQ=760; %(m3/d)- caudal de AR a tratarQR=40; %(m3/d)- Caudal de reflujoQw=15; %(m3/d)- caudal de desecho del sedimentadorV=500; %(m3)- volumen de proceso del reactorVs=20; %(m3)-volumen del sedimentadordy=zeros(4,1);dy(1)=Q*S/V - (Q+QR)*y(1)/V - (umax/Y)*y(1)*y(2)/(Ks+y(1));dy(2)=QR*y(3)/V - (Q+QR)*y(2)/V + umax*y(1)*y(2)/(Ks+y(1)) - Kd*y(2);dy(3)=(Q+QR)*y(2)/Vs - (QR+Qw)*y(3)/Vs;dy(4)=Q*CO2/V - (Q+QR)*y(4)/V - (umax/Yo)*y(1)*y(2)/(Ks+y(1))-b*y(1)+kla*(Csr-y(4));


function dy=lodoactivadoQw3(t,y)umax=1.97; %(d-1)b=0.259; %(d-1)Y=0.65; %Coef.rendimiento (mg SSV producido/mg DQO consumidoKs=137.3; %(mg/L)- Constante de afinidadYo=0.915; %Yo=a/YKd=0.0601; %(d-1)- Coeficiente de muerteS=668; %(mg/L)- Concentración del aflueunte a procesarCO2=0.3; %(mg/L)- Cocnentración de oxígeno en el alfuente

Csr=7.02; %(Concent.de saturacion de oxig.en el agua residual)kla=420; %(d-1)- Coeficiente de transferencia de oxígenoQ=760; %(m3/d)- caudal de AR a tratarQR=40; %(m3/d)- Caudal de reflujoQw=35; %(m3/d)- caudal de desecho del sedimentadorV=500; %(m3)- volumen de proceso del reactorVs=20; %(m3)-volumen del sedimentadordy=zeros(4,1);dy(1)=Q*S/V - (Q+QR)*y(1)/V - (umax/Y)*y(1)*y(2)/(Ks+y(1));dy(2)=QR*y(3)/V - (Q+QR)*y(2)/V + umax*y(1)*y(2)/(Ks+y(1)) - Kd*y(2);dy(3)=(Q+QR)*y(2)/Vs - (QR+Qw)*y(3)/Vs;dy(4)=Q*CO2/V - (Q+QR)*y(4)/V - (umax/Yo)*y(1)*y(2)/(Ks+y(1))-b*y(1)+kla*(Csr-y(4));


programa ar

Documents