Scribd
Cargar
766 vistas3 páginas

Ejercicio en Matlab

Este documento presenta un ejercicio de modelado y simulación ambiental en MATLAB. Se analiza un río que recibe un vertido de aguas residuales a 40 km de su desembocadura. Se dan los caudales, temperaturas y concentraciones de materia orgánica y oxígeno disuelto del río y vertido. Se calculan estos valores en el punto de vertido y las constantes de biodegradación a esa temperatura. Luego, se modela la variación de la DBO y oxígeno disuelto a lo largo de los 40 km hasta la de

Cargado por

ERICKA
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
766 vistas3 páginas

Ejercicio en Matlab

Este documento presenta un ejercicio de modelado y simulación ambiental en MATLAB. Se analiza un río que recibe un vertido de aguas residuales a 40 km de su desembocadura. Se dan los caudales, temperaturas y concentraciones de materia orgánica y oxígeno disuelto del río y vertido. Se calculan estos valores en el punto de vertido y las constantes de biodegradación a esa temperatura. Luego, se modela la variación de la DBO y oxígeno disuelto a lo largo de los 40 km hasta la de

Cargado por

ERICKA
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA QUÍMICA

ESCUELA ACADÉMICO PROFESIONAL DE INGENIERÍA

AMBIENTAL

“EJERCICIOS EN MATLAB”

ALUMMNA:

VILLANUEVA SAAVEDRA, ERICKA YESSENIA

DOCENTE:

Mg. MONCADA TORRES ,LUIS DAVID

CURSO:

MODELAMIENTO Y SIMULACIÓN

AMBIENTAL

12 DE NOVIEMBRE DEL 2019

TRUJILLO-PERÚ
UNIVERSIDAD NACIONAL
DE TRUJILLO

Ejercicio 2:
Un río recibe un vertido de aguas residuales 40 km aguas arriba de su desembocadura al mar. Si los
caudales, temperaturas, y concentraciones de materia orgánica y oxígeno disuelto del río antes del vertido
y del vertido son los que se indican en la siguiente tabla, encuentra la concentración de materia orgánica y
oxígeno disuelto en la desembocadura del río.

Rio antes del vertido Vertido


Caudal en m3/s 5.787 0.463
Temperatura en °C 20 28
Oxígeno disuelto en mg/L 7.5 2
Concentración de materia 2 200
orgánica (DBO) en mg/L

𝑘1 (20°𝐶) = 0.5𝑑−1
𝑘2 (20°C) = 1.84 𝑑−1
𝜃1 = 1.047
𝜃2 = 1.024
Velocidad de movimiento = 0.48 m/s

DESARROLLO EN MATLAB
% Datos del río antes del vertido
Q1=5.787; % Caudal (m3/s)
T1=20; % Temperatura (°C)
Od1=7.5; % Oxígeno disuelto (mg/L)
L1=2; % Concentración de materia orgánica (mg/L)
% Datos del vertido
Q2=0.463; % Caudal (m3/s)
T2=28; % Temperatura (°C)
Od2=2; % Oxígeno disuelto (mg/L)
L2=200; % Concentración de materia orgánica (mg/L)
% Datos adicionales
k1=0.5; % Constante de biodegradación para T=20°C (d^-1)
k2=1.84; % Tasa de aireación de oxígeno en el cuerpo de agua para T=20°C (d^-
1)
a1=1.047;
a2=1.024;
d=40; % Distancia del vertido a la desembocadura (km)
vmov=0.48; % Velocidad de movimieto (m/s)
% Cálculo de valores en el punto de vertido
Q0=Q1+Q2;
L0=((L1*Q1)+(L2*Q2))/Q0;
T0=((T1*Q1)+(T2*Q2))/Q0;
Od0=((Od1*Q1)+(Od2*Q2))/Q0;
% Cálculo del Oxígeno disuelto saturado
Te=T0+273; % Convertir la temperatura a grados Kelvin
Ods=2.718282^(-139.34410+((1.575701*(10^5))/Te)-
((6.642308*(10^7))/(Te^2))+((1.2438*(10^10))/(Te^3))-
((8.621949*(10^11))/(Te^4))); % (mg/L)
% Cálculo de las constantes para la nueva temperatura T0
K1=k1*(a1^(T0-T1));
K2=k2*(a2^(T0-T1));
% Cálculo de la DBO y Od en la desembocadura
Página 2|3
UNIVERSIDAD NACIONAL
DE TRUJILLO

tp=((d*1000)/vmov)/60; % Se convierte la distancia a metros, y el tiempo


queda en minutos
f=@(t,x)[(-K1*(x(1))),((K1*(x(1)))-(K2*(Ods-(x(2)))))];
x0=[L0,Od0];
tspan=[0,tp];
[t,f]=ode45(f,tspan,x0);
yyaxis left
plot(t,f(:,1));
yyaxis right
plot(t,f(:,2));
yyaxis left
title('Modelo de DBO y Od')
xlabel('Tiempo (días)|')
ylabel('DBO (mg/l)')
yyaxis right
ylabel('oxígeno disuelto (mg/l)')

Página 3|3

También podría gustarte