Procesos de Tratamientos de Aguas Residuales y PDF

Proyecto Final de Carrera
Procesos de tratamientos de
aguas residuales y simulacin con
MatLab
Autora: Miriam Mrquez

Ponentes: Moiss Graells
Montserrat Prez
Procesos de tratamientos de aguas residuales y simulacin con MatLab
NDICE
1. RESUMEN .......................................................................................................................................8
1.1 Castellano ......................................................................................................................................8
1.2 Catal.............................................................................................................................................8
1.3 English ...........................................................................................................................................8
2. MOTIVACIN................................................................................................................................9
3. OBJETIVOS ....................................................................................................................................9
4. INTRODUCCIN .........................................................................................................................10
4.1 Espacio Europeo de Educacin Superior (E.E.E.S.)....................................................................10
4.2 Gestin del agua ..........................................................................................................................11
4.3 Simulacin de Procesos ...............................................................................................................12
4.3.1. Descripcin de la ventana inicial ..........................................................................................13
5. DISEO DE LA ASIGNATURA .................................................................................................14

5.1 Perfil del alumno..........................................................................................................................15
5.2 Estructura ....................................................................................................................................15
5.3 Temas ...........................................................................................................................................16
Tema 1 : SOLUBILIDAD ...........................................................................................................17

Prctica 1: SOLUBILIDAD .....................................................................................................18
Horas de trabajo real presencial................................................................................................18
Resumen ...................................................................................................................................18
Objetivos ..................................................................................................................................18
Fsico-Qumicos....................................................................................................................18
MatLab/Informticos ............................................................................................................18
Fundamentos qumicos .............................................................................................................19
Introduccin ..........................................................................................................................19
Solubilidad del oxgeno ........................................................................................................19
Presin parcial del oxgeno igual a 1 atmsfera ................................................................19
Presin parcial del oxgeno diferente a 1 atmsfera ..........................................................19
Programa en MatLab ................................................................................................................21
IDEAS PREVIAS .................................................................................................................21
Conceptos generales ..........................................................................................................21
Clear...............................................................................................................................21
Punto y coma .................................................................................................................21
Ans.................................................................................................................................21
Aadir comentarios........................................................................................................21
Definicin de matrices, vectores y escalares. ....................................................................21
Matrices .........................................................................................................................21
Vectores .........................................................................................................................21
Otras formas de introducir vectores ...............................................................................21
Escalares ........................................................................................................................22
Direccionamiento de matrices y vectores ..........................................................................22
Acceso elemento a elemento..........................................................................................22
Acceso en bloque ...........................................................................................................22
Operaciones de matrices, vectores y escalares ..................................................................22
Elemento-a-elemento .....................................................................................................22
Orientacin del vector o matriz......................................................................................23
Controles de flujo o bucles ................................................................................................23
2
for (desde)......................................................................................................................23
Funciones ..........................................................................................................................23
Funcin lenght ...............................................................................................................23
Funcin abs....................................................................................................................23
Funcin exp ...................................................................................................................23
Graficar..............................................................................................................................23
Opciones ........................................................................................................................23
Estilo de lneas, marcadores y colores ...........................................................................23
Adicin de rejillas y etiquetas........................................................................................24
Ejercicio ................................................................................................................................24
Solucin.............................................................................................................................24
Grfica ...........................................................................................................................25
Tema 2 : AUTODEPURACIN..................................................................................................26
Prctica 2: AUTODEPURACIN ...........................................................................................27
Resumen ...................................................................................................................................27
Objetivos ..................................................................................................................................27
Fsico-Qumicos....................................................................................................................27
Conocimientos bsicos .............................................................................................................28
Biooxidacin .....................................................................................................................28
Biodegradacin..................................................................................................................28
Introduccin ..........................................................................................................................28
Autodepuracin de los ros ...................................................................................................28
Desoxigenacin.....................................................................................................................29
Reoxigenacin.......................................................................................................................30
Dficit de concentracin de oxgeno.....................................................................................30
Determinacin de k1 ..........................................................................................................30
Determinacin de k2 ..........................................................................................................30
Evolucin de la concentracin de oxgeno causada por un vertido.......................................31
IDEAS PREVIAS .................................................................................................................32
Clf ..................................................................................................................................32
Global ............................................................................................................................32
Definicin de matrices, vectores y escalares .....................................................................32
Archivos *.m .....................................................................................................................32
Como crearlos................................................................................................................32
Como ejecutarlos ...........................................................................................................32
Como guiar a MatLab para encontrarlos........................................................................32
Funciones ..........................................................................................................................33
Length ............................................................................................................................33
Optimset.........................................................................................................................33
Fminsearh ......................................................................................................................33
Num2sr ..........................................................................................................................33
Round.............................................................................................................................33
Sum................................................................................................................................33
Graficar..............................................................................................................................33
Funcin Hold .................................................................................................................34
Funcin line ...................................................................................................................34
Ejes a medida .................................................................................................................34
Ejercicio ................................................................................................................................34
Solucin.............................................................................................................................34
Grfica ...........................................................................................................................35
3
Tema 3 : FRMULA EMPRICA...............................................................................................36

Prctica 3: DETERMINACIN DE LA FRMULA EMPRICA DE LOS COMPONENTES
DE LAS AGUAS .........................................................................................................................37
Resumen ...................................................................................................................................37
Objetivos ..................................................................................................................................37
Fsico-Qumicos....................................................................................................................37
Demanda Qumica de Oxigeno (DQO) .............................................................................38
Carbono orgnico total (COT)...........................................................................................38
Introduccin ..........................................................................................................................38
Clculo de la frmula emprica del agua residual .................................................................38
IDEAS PREVIAS .................................................................................................................40
Definicin matrices, vectores y escalares..........................................................................40
Operaciones con matrices, vectores y escalares ................................................................40
Funciones ..........................................................................................................................40
Double ...........................................................................................................................40
Round.............................................................................................................................40
Min.................................................................................................................................40
Resolucin de ecuaciones algebraicas...............................................................................40
Ejercicio ................................................................................................................................41
Solucin.............................................................................................................................41
Tema 4 : DEMANDA BIOQUMICA DE OXGENO ...............................................................42

Prctica 4: PARMETROS CINTICOS DE LA DBO..........................................................43
Resumen ...................................................................................................................................43
Objetivos ..................................................................................................................................43
Fsico-Qumicos....................................................................................................................43
Demana Biologica de Oxgeno..........................................................................................44
Conceptos relacionados con la DBO.....................................................................................44
Cintica de la DBO ...............................................................................................................45
Determinacin de las constantes cinticas de la DBO ..........................................................46
IDEAS PREVIAS .................................................................................................................47
Manipulacin matricial ..................................................................................................47
Funciones ..........................................................................................................................47
Diff.................................................................................................................................47
Log.................................................................................................................................47
Length ............................................................................................................................47
Sum................................................................................................................................47
Exp.................................................................................................................................47
Round.............................................................................................................................47
Num2str .........................................................................................................................47
Graficar..............................................................................................................................47
Subplot...........................................................................................................................47
Ejercicio ................................................................................................................................48
Solucin.............................................................................................................................48
Grfica ...........................................................................................................................50
4
Tema 5 : DEMANDA QUMICA DE OXIGENO ......................................................................51

Prctica 5: DEMANDA QUMICA DE OXGENO................................................................52
Resumen ...................................................................................................................................52
Objetivos ..................................................................................................................................52
Fsico-Qumicos....................................................................................................................52
Blanco................................................................................................................................53
Patrn primario..................................................................................................................53
Medida experimental de la DQO ..........................................................................................53
Definicin de la DQO ...........................................................................................................54
Demanda qumica de oxgeno de sustancias puras ...............................................................54
Patrn para la medida de la DQO .........................................................................................55
Interferencias en la determinacin de la DQO ......................................................................55
Sales ferrosas.....................................................................................................................55
Cloruros.............................................................................................................................55
Procedimiento y reactivos de la medida de la DQO .............................................................56
Para el rango de 150 a 1500 ..............................................................................................56
Preparacin del reactivo de dicromato...........................................................................56
Tcnica operatoria para la determinacin de la DQO ....................................................56
Para el rango de 0 a 150 ....................................................................................................57
Determinacin de la absortividad..........................................................................................58
IDEAS PREVIAS .................................................................................................................60
Operar con matrices, vectores y escalares .........................................................................60
Funciones ..........................................................................................................................60
Sum ...................................................................................................................................60
Input...............................................................................................................................60
Num2str .........................................................................................................................60
Length ............................................................................................................................60
Round.............................................................................................................................60
Graficar..............................................................................................................................60
Ejes a medida .................................................................................................................60
Line................................................................................................................................60
Ejercicio ................................................................................................................................61
Solucin.............................................................................................................................61
Grfica ...........................................................................................................................62
Tema 6 : MODELOS DINMICOS DE TANQUES..................................................................63

Resumen ...................................................................................................................................64
Introduccin..............................................................................................................................64
Prctica 6.1 : ALTURA DE LMINA DE VERTIDO Y VOLUMEN DEL TANQUE .........65
Resumen ...................................................................................................................................65
Objetivos ..................................................................................................................................65
Fsico-Qumicos....................................................................................................................65
IDEAS PREVIAS .................................................................................................................68
Archivos *.m .....................................................................................................................68
5
if (condicional)...............................................................................................................68
Funciones ..........................................................................................................................68
Odeset ............................................................................................................................68
Ode23.............................................................................................................................68
Length ............................................................................................................................69
Zeros ..............................................................................................................................69
Graficar..............................................................................................................................69
Ejercicio ................................................................................................................................69
Solucin.............................................................................................................................70
Grfica ...........................................................................................................................72
Prctica 6.2 : TANQUE DE MEZCLA PERFECTA ...............................................................73
Resumen ...................................................................................................................................73
Objetivos ..................................................................................................................................73
Fsico-Qumicos....................................................................................................................73
IDEAS PREVIAS .................................................................................................................75
Archivos *.m .....................................................................................................................75
if (condicional)...............................................................................................................75
Funciones ..........................................................................................................................75
Odeset ............................................................................................................................75
Ode23.............................................................................................................................75
Length ............................................................................................................................75
Zeros ..............................................................................................................................75
Graficar..............................................................................................................................75
Ejercicio ................................................................................................................................76
Solucin.............................................................................................................................76
Grfica ...........................................................................................................................79
Prctica 6.3 : DISEO DE TANQUES REGULADORES DE CAUDAL..............................80
Resumen ...................................................................................................................................80
Objetivos ..................................................................................................................................80
Fsico-Qumicos....................................................................................................................80
Mediante el grfico Volumen-tiempo ...................................................................................82
Mediante la variacin del caudal frente al tiempo ................................................................83
IDEAS PREVIAS .................................................................................................................84
Archivos *.m .....................................................................................................................84
Funciones ..........................................................................................................................84
Odeset ............................................................................................................................84
Round.............................................................................................................................84
Length ............................................................................................................................84
Num2str .........................................................................................................................84
Spline .............................................................................................................................84
Fzero ..............................................................................................................................84
Graficar..............................................................................................................................84
6
Ejercicio ................................................................................................................................85
Solucin.............................................................................................................................85
Grfica ...........................................................................................................................86
5.4 Conversin a crditos ECTS ........................................................................................................87

5.5 Sistema de evaluacin ..................................................................................................................88
5.6 Logstica.......................................................................................................................................89
5.6.1. Coordinacin.........................................................................................................................89
5.6.2. Sitio de trabajo ......................................................................................................................89
5.7 Bibliografa recomendada ...........................................................................................................90
5.8 Presupuesto..................................................................................................................................91
5.9 Conclusiones ................................................................................................................................93
6. REFERENCIAS.............................................................................................................................94
7. AGRADECIMIENTOS.................................................................................................................95
7
1. RESUMEN
1.1 Castellano
El proyecto presente consiste en el diseo de una asignatura de 4,5 crditos ECTS de la
titulacin de la Ingeniera Tcnica Industrial, especialidad en Qumica Industrial, aplicable para
toda la enseanza Europea.
El objetivo de esta asignatura es combinar los contenidos sobre la caracterizacin de aguas
residuales y sus tratamientos y la metodologa prctica de resolucin de ejercicios mediante
tcnicas de simulacin (MatLab).
1.2 Catal
El projecte present consisteix en el disseny d'una assignatura de 4,5 crdits ECTS de la
titulaci de l'Enginyeria Tcnica Industrial, especialitat en Qumica Industrial, aplicable per a tot
l'ensenyament Europeu.
L'objectiu d'aquesta assignatura s combinar els continguts sobre la caracteritzaci d'aiges
residuals i els seus tractaments i la metodologia prctica de resoluci d'exercicis mitjanant
tcniques de simulaci (MatLab).
1.3 English
The present project consists of the design of a course of 4.5 ECTS credits of the degree of
Industrial Technical Engineering, speciality in Industrial Chemistry, applicable for all European
education.
The objective of this subject is to combine the contents on the waste water characterization and
its treatments and the practical methodology of resolution of exercises by means of simulation
techniques (MatLab).
8
2. MOTIVACIN
La eleccin del diseo de una asignatura de Ingeniera Tcnica Industrial, especialidad en
Qumica industrial que combine los contenidos sobre la caracterizacin de aguas residuales y
sus tratamientos y la metodologa prctica de resolucin de ejercicios mediante tcnicas de
simulacin (MatLab) como proyecto final de carrera nace de las ganas de aprender el
funcionamiento y aplicaciones de nuevos programas informticos, y por la satisfaccin de ser
promotora de este tipo de asignaturas, asignaturas menos tradicionales y actualizadas a las
nuevas tecnologas.
Esta asignatura se debe implantar en la enseanza Europea ya que trata de un tema muy
presente en la sociedad actual, el agua y los tratamientos de depuracin a los que se somete.
Por otra parte, este tema al ser tratado con mtodos informticos, tan desarrollados y utilizados
en la actualidad, estimula a los alumnos al uso de las simulaciones en situaciones reales
ahorrando, as, tiempo, espacio y dinero.
3. OBJETIVOS
El principal objetivo del PFC, Proyecto de Final de Carrera, es que el alumno aprenda, tanto en
el campo que estudia (Ingeniera Tcnica Industrial, especialidad en Qumica Industrial, en mi
caso), como en el organizativo (planificar, estructurar y disear un proyecto).
Por otro lado, los objetivos de este proyecto, disear una herramienta, asignatura, capaz de:
- Desarrollar los conocimientos qumicos y fisico-qumicos que aparecen en los
procesos de determinacin de la contaminacin del agua.
- Progresar en los fundamentos de la Ingenieria Tcnica Industrial, especializada
en Qumica industrial necesarios para el diseo y operacin de los procesos de
tratamientos de agua.
- Estimular la simulacin por ordenador, mediante el programa MatLab.
- Garantizar su aplicacin dentro del Espacio Europeo de Educacin Superior.
9
4. INTRODUCCIN
Para introducir el diseo de la asignatura se ha de mencionar sus tres principales ideas:
1. El contexto en que se encuentra, y por tanto las pautas que ha de seguir la

asigntaura:
El tratado de Bolonia
2. El tema qumico y fisico-qumico a estudiar:

Los procesos de determinacin de la contaminacin del agua.
3. El mtodo utilizado para dicho estudio:

El programa informtico MatLab.
4.1 Espacio Europeo de Educacin Superior (E.E.E.S.)

La sociedad del siglo XXI poco tiene que ver con la de hace 100 aos. Esta diferencia se debe
a los cambios en el sistema de comunicaciones, la adquisicin de la informacin, globalizacin
de nuevas tecnologas, el crecimiento intelectual en la sociedad,
Frente a estos cambios, E.E.E.S. se cuestiona el sistema de Educacin Universitaria,
planteando el abandono de la enseanza tradicional.
En el ao 1998, donde los Ministros de Educacin de Francia, Alemania, Italia y Reino unido
firman la Declaracin de la Sorbona1. Desde entonces, en Europa se ha iniciado un proceso
para promover la convergencia entre los sistemas nacionales de educacin superior. Los
Ministerios de cada pas miembro de la unin han refrendado2, con la firma de la Declaracin
de Bolonia (1999), la importancia de un desarrollo armnico de un Espacio Europeo de
Educacin Superior antes del 2010. Este objetivo est propiciando diversos foros de debate en
todos los pases miembros, para la evaluacin de programas acadmicos convergentes que
aseguren una calidad docente, adoptando el sistema de transferencia de crditos que permitir
un reconocimiento acadmico inmediato de ttulos, una movilidad entre pases y un aprendizaje
a lo largo de toda la vida.
La Declaracin de Bolonia incluye entre sus principales objetivos:
- La adopcin de un sistema fcilmente legible y comparable de titulaciones,
mediante la implantacin, entre otras cosas, de un Suplemento al Diploma.
- La adopcin de un sistema basado, fundamentalmente, en dos ciclos
principales: pregrado y grado.
El ttulo otorgado al terminar el primer ciclo tendr que tener un valor especfico
en el mercado de trabajo Europeo. El segundo ciclo llevar a la obtencin de
un Master y/o Doctorado como ocurre en muchos estados Europeos.
- El establecimiento de un sistema de crditos, como el sistema ECTS.
- La promocin de la cooperacin Europea para asegurar un nivel de calidad
para el desarrollo de criterios y metodologas comparables.
- La promocin de una necesaria dimensin Europea en la educacin superior
con particular nfasis en el desarrollo curricular.
Todos estos objetivos estn siendo discutidos, a distintos niveles, en la mayora de los estados
de la Unin Europea, tratando de llegar a la deseada "armonizacin" preconizada3 desde la
Declaracin de la Sorbona.
1
Declaracin conjunta para la armonizacin del diseo del Sistema de Educacin Superior Europeo.
2
Contenido, dominado o hacer menos violento.
3
Defendida, apoyada.
10
En Espaa, los principales puntos a tener en cuenta, para alcanzar una adaptacin del sistema
actual de docencia y crditos, seran, en una primera fase:
- La implantacin de un sistema de crditos europeos.
- La adopcin de un sistema de calificaciones que permitan una conversin fcil
al sistema de calificaciones y grados ECTS, Sistema europeo de transferencia
de crditos (European Credit Transfer System).
- La implantacin del Suplemento al Diploma. (1)
Por el momento, el objetivo ms avanzado es la implantacin del ECTS. Y por tanto, es la
condicin que aplicaremos al diseo de la asignatura.
El ECTS se ha desarrollado para facilitar el reconocimiento acadmico de los estudios

cursados en cualquier lugar de Europa.
Este sistema de crditos se ha aplicado y ha mostrado su eficacia, como proyecto piloto, en el
marco del programa ERASMUS. Visto el xito, el sistema se quiere llegar a aplicar para todos
los estudiantes de la Unin Europea, de tal manera que el trabajo desarrollado por un
estudiante sea fcilmente reconocible en cuanto a nivel, calidad y formacin en todos los
estados (reconocimiento acadmico).
Los crditos ECTS representan, en forma de valor numrico (entre 1 y 60). 60 crditos
representan el volumen de trabajo de un ao acadmico, por tanto, 30 crditos equivalen a un
semestre y 20 crditos a un trimestre de estudios.
Dichos crditos se basan en el volumen total de trabajo del estudiante y no se limita
exclusivamente a las horas de asistencia, sino al trabajo real del estudiante (workload).
Un crdito ECTS equivale a 26,7 horas de trabajo. Esta relacin la encontramos estimando que
el alumno trabaja 40h por semana, durante 40 semanas, lo que hacen un total de 1600 horas
por curso acadmico.
Otro requisito a tener en cuenta y seguir las pautas marcadas por el tratado es la perfeccin del
sistema de calificaciones.
Por el momento, la restriccin ms clara es la introduccin de notas numricas con un nico
decimal.
4.2 Gestin del agua

El agua se identifica como motor de desarrollo y fuente de riqueza, y ha constituido uno de los
pilares fundamentales para el progreso del hombre. La ordenacin y gestin de los recursos
hdricos, que ha sido desde siempre un objetivo prioritario para cualquier sociedad, se ha
realizado histricamente bajo directrices orientadas a satisfacer la demanda en cantidades
suficientes, bajo una perspectiva de poltica de oferta.
El incremento de la oferta de agua como herramienta para el impulso econmico, el mayor
nivel de contaminacin, irremisiblemente asociado a un mayor nivel de desarrollo, algunas
caractersticas naturales (sequas prolongadas, inundaciones) y en definitiva una
sobreexplotacin de los recursos hdricos, han conducido a un deterioro importante de los
mismos.
Esto ha hecho necesario un cambio en los planteamientos sobre poltica de aguas, que han
tenido que evolucionar desde una simple satisfaccin en cantidad de las demandas, hacia una
gestin que contempla la calidad del recurso y la proteccin del mismo como garanta de un
abastecimiento futuro y de un desarrollo sostenible. Desde entonces, la poltica de agua sufre
un cambio importante en los conceptos y criterios utilizados en la planificacin hidrolgica e
introducen la calidad de las aguas y la proteccin de los recursos hdricos como puntos
fundamentales para estructurar dicha planificacin. (2)
Para conseguir la calidad adecuada en las aguas, estas se han de someter a procesos de
tratamientos de depuracin.
Estos tratamientos, como es evidente, tienen como objetivo conseguir, a partir de aguas
servidas y mediante diferentes procedimientos fsicos, qumicos y biolgicos, un agua efluente
de mejores caractersticas de calidad, tomando como base ciertos parmetros normalizados,
que se pueden estimar con indicadores como la DBO5 (demanda biolgica de oxgeno a los
11
cinco das) y la DQO (demanda qumica de oxgeno), que son las cantidades de oxgeno que
se necesitan para oxidar la materia orgnica susceptible de ser oxidada bien por va orgnica
(bacterias y microorganismos) o bien por va qumica.
Otro de los parmetros comunes es la cantidad de slidos en suspensin totales (SST) que da
una idea de la cantidad de materia contaminante slida del agua. (3)
Los parmetros caractersticos son:

MENSURABLES:
- Temperatura
- pH
- Slidos en suspensin totales (SST)
- Materia orgnica valorada como DQO y DBO (a veces TOC)
- Nitrgeno total Kjeldahl (NTK)
- Nitrgeno amoniacal y nitratos
- Fsforo total, nitritos, sulfuros, slidos disueltos
NO MENSURABLES:
- Olor
- Sabor
4.3 Simulacin de Procesos

La simulacin virtual de sistemas esta basada en una generalizacin del concepto de
experimentacin del mtodo cientfico, segn el cual en lugar de realizar los experimentos
sobre el sistema real, se realizan sobre un algoritmo matemtico que lo representa, de manera
que si el modelo es una representacin vlida del sistema en estudio, entonces los resultados
de la experimentacin con el modelo pueden transferirse al diseo real del prototipo del
sistema deseado. (4)
El uso de simulacin en computadoras para modelar los sistemas se ha vuelto popular en los
ltimos aos, y uno de los programas destacados en esta categora es el MatLab.
MatLab es el nombre abreviado de MATrix LABoratory. Este programa se utiliza para realizar
clculos numricos con vectores y matrices. Como caso particular puede tambin trabajar con
nmeros escalares (tanto reales como complejos), con cadenas de caracteres y con otras
estructuras de informacin ms complejas.
Una de las capacidades ms atractivas es la de realizar una amplia variedad de grficos en dos
y tres dimensiones.
MatLab dispone de un cdigo bsico y de varias libreras especializadas (toolboxes).
Esta asignatura utiliza la versin 7.0.4 de este programa, aparecida a mediados de 2004.
Figura 1. Icono de MatLab 7.0.4
12
Figura 2. Ventana inicial
Al arrancar el programa se abre una ventana similar a la mostrada en la Figura 2. sta es la

vista que se obtiene eligiendo la opcin Desktop Layout/Default, en el men View. Como esta
configuracin puede ser cambiada fcilmente por el usuario, es posible que en muchos casos
concretos lo que aparezca sea muy diferente.
4.3.1. Descripcin de la ventana inicial

La parte ms importante de la ventana inicial es la Command Window, que aparece en la parte
derecha. En esta sub-ventana es donde se ejecutan los comandos de MatLab, a continuacin
del prompt (aviso) caracterstico (>>), que indica que el programa est preparado para recibir
instrucciones.
En la parte superior izquierda de la pantalla aparecen dos ventanas tambin muy tiles: en la
parte superior aparece la ventana Current Directory, que se puede alternar con Workspace
clicando en la pestaa correspondiente.
La ventana Current Directory muestra los ficheros del directorio activo o actual. Dicho directorio
se puede cambiar desde la propia ventana o desde la barra de herramientas, debajo de la
barra de mens con los mtodos de navegacin de directorios propios de Windows. Clicando
dos veces sobre alguno de los ficheros *.m del directorio activo se abre el editor de ficheros de
MatLab, herramienta fundamental para la programacin.
El Workspace contiene informacin sobre todas las variables que se definen en una sesin y
permite ver y modificar las matrices con las que se trabaje.
En la parte inferior derecha aparece la ventana Command History que muestra los ltimos
comandos ejecutados en la Command Window.
Clicando sobre un comando con el botn secundario, se muestra un men contextual con las
posibilidades disponibles en ese momento.
Para editar uno de estos comandos hay que copiarlo antes a la Command Window.
En la parte inferior izquierda de la pantalla aparece el botn Start, con una funcin semejante a
la del botn Inicio de Windows. (5)
13
5. DISEO DE LA ASIGNATURA
Para disear satisfactoriamente una asignatura se han de analizar los factores que intervienen
en dicho diseo:
1. Perfil del alumno

2. Estructura
3. Temas
4. Conversin a crditos ECTS
5. Sistema de evaluacin
6. Logstica
7. Bibliografia recomendable
8. Presupuesto
9. Conclusiones
14
5.1 Perfil del alumno

Es importante saber a que tipo de alumno va dirigida la asignatura diseada, ya que esto es lo
que limita el contenido.
Para seguir sin problemas esta asignatura, el alumno deber tener conocimientos universitarios
de:
- Qumica Analtica.
- Fisico-qumica.
- Matemticas.
- Informtica.
5.2 Estructura
La asignatura consta de 6 temas relacionados con los procesos de tratamientos de depuracin
de las aguas, y siguen el siguiente orden:
- Tema 1 y 2: Estudian los conceptos de la solubilidad y la autodepuracin.
- Tema 3: Detalla la determinacin de la frmula emprica de los componentes

de las aguas residuales.
- Tema 4 y 5: Tratan los parmetros ms importantes en la determinacin de la

contaminacin de un agua: la DBO y la DQO.
- Tema 6: Introduce las simulaciones de tanques.
Para un mximo rendimiento del alumnado todos los temas son independientes, evitando as
el retraso del alumno que no comprenda uno de los temas tratados.
Y para facilitar la comprensin, cada tema est definido por prcticas que siguen el siguiente
orden:
- Ttulo del tema a tratar.

- Horas de trabajo real presencial.
- Resumen.
- Objetivos:
Fsico-Qumicos.
MatLab/Informticos.
- Conocimientos bsicos*.
- Fundamentos qumicos.
- Programa en MatLab:
Ideas previas.
Ejercicio.
o Solucin (programacin).
Aclaraciones.
*Este apartado no aparece en alguna prctica.
15
5.3 Temas
horas
Se considera que la duracin ptima de esta asignatura es de 1 clase de 2,5 , durante
semana
15 semanas. En ese perodo deben tratarse y asumirse los diferentes temas:
Tema 1 : SOLUBILIDAD...................................................................................................................17
Tema 2 : AUTODEPURACIN ........................................................................................................26
Tema 3 : FRMULA EMPRICA ....................................................................................................36
Tema 4 : DEMANDA BIOQUMICA DE OXGENO ....................................................................42
Tema 5 : DEMANDA QUMICA DE OXIGENO ...........................................................................51
Tema 6 : MODELOS DINMICOS DE TANQUES.......................................................................63
16
Tema 1 :
SOLUBILIDAD
17
Prctica 1: SOLUBILIDAD
Horas de trabajo real presencial

Qumico: 1,5h
Informtico/MatLab: 2,5h
TOTAL 4,0h
Resumen
Se conocer los parmetros que influyen en la solubilidad del oxgeno en el agua, la ecuacin
matemtica que las relaciona y, con la ayuda del MatLab, se representar la cloronidad en
funcin de la concentracin de oxgeno y de la Temperatura, a 1 atm.
Objetivos
La persona que haya realizado correctamente esta actividad de aprendizaje deber ser capaz
de:
Fsico-Qumicos
- Conocer los factores que influyen en la solubilidad del oxgeno en el agua.
- Reconocer la relacin existente entre la temperatura y la solubilidad del
oxgeno disuelto.
MatLab/Informticos
- Definir vectores lineales (linspace).
- Direccionar matrices y vectores.
- Operar con matrices, vectores y escalares.
- Ejecutar las funciones:
lenght
abs
exp
- Crear bucles, controles de flujo (for).
- Graficar (plot).
- Modificar opciones de los graficos:
Color
LineWidth
fontsize
- Aadir etiquetas
18
Fundamentos qumicos
Introduccin
Una de las primeras manifestaciones del efecto de la contaminacin orgnica de las aguas es
la disminucin brusca de la concentracin de oxgeno disuelto, en el lugar de los vertidos,
debido a que el oxgeno es utilizado, y por tanto retirado del medio, por los microorganismos
que proliferan4 en las aguas contaminadas, impidiendo el desarrollo de la vida de las especies
superiores, no slo por el dficit de oxgeno, sino tambin por la insalubridad5 que acarrea la
poblacin microbiana resultante.
Solubilidad del oxgeno

Las aguas prepotables en contacto con el aire se oxigenan, con un nivel de concentracin de
oxgeno disuelto a saturacin, que depende principalmente de:
- La presin parcial del oxgeno.
- La temperatura.
- La salinidad del agua.
Presin parcial del oxgeno igual a 1 atmsfera

Para agua con sales, la concentracin de oxgeno disuelto a saturacin, C*, en un agua en
equilibrio con el oxgeno del aire a 1 atmsfera de presin, est expresada por la siguiente
ecuacin emprica:
1,5757 10 5 6,642310 7 1,24381010 8,6219 1011

log e C* = 139,344 + +
T T2 T3 T4 [1.1]
19,428 3867,3
Cl (3,1929 10 2 + )
T T2
Donde:
T: Temperatura absoluta.
Cl: Clorinidad equivalente6 a KCl, mg/L (salinidad).
Presin parcial del oxgeno diferente a 1 atmsfera

En este caso, la solubilidad del oxgeno disuelto se corrige con la expresin:
(1 Pw / P) (1 P) [1.2]
C p = C * p
(1 Pw) (1 )
Donde:
Cp: Solubilidad del oxgeno a la presin P[atm].
Pw: Presin de vapor del agua, atm.
: 9,75 * 10-4 -1,426 * 10-5 * t + 6,436 * 10-8 * t2, t C.
4
Multiplicar abundantemente.
5 No saludable.
6
La clorinidad equivalente en KCl, es la disolucin de KCl que produce la misma conductividad elctrica
que la disolucin problema.
19
La presin de vapor del agua, Pw, y la presin de vapor a saturacin del agua, Pws, tensin de
vapor, estn relacionados a travs de la humedad relativa, landa, segn la expresin:
Pw
= [1.3]
Pws
La presin de vapor del agua, en mm de Hg, en funcin de la temperatura centgrada, viene

dada por la expresin emprica obtenida:
22, 4077t +1167, 7614
Pws = 0,040032 e t + 245, 9443

[1.4]
20
Programa en MatLab
IDEAS PREVIAS
Conceptos generales
Clear
Formatea el workspace (Tambin utilizado clear all).
Punto y coma
El punto y como se coloca al final de cada orden, dependiendo si queremos o no, que MatLab nos
muestre en el Command Window todo lo procesado.
Ans
Answer aparecer siempre que no se declar a ninguna varialbe el resultado de una operacin, MatLab
utiliza este nombre de variable especial por defecto (ans).
Ejemplo:
Definicin de una matriz A de dimensin (33):
>> 123+45
La respuesta del programa es la siguiente:
ans = 168
Aadir comentarios
Para escribir comentarios en la programacin de MatLab, basta con aadir el smbolo del porcentaje (%)
antes del comentario.
Definicin de matrices, vectores y escalares.

Generalmente se utilizan letras maysculas para matrices y letras minsculas para vectores y escalares.
Matrices
Las matrices se definen mediante corchetes y se introducen por filas.
Los elementos de una misma fila estn separados por espacios o comas (vector fila), mientras que las
filas estn separadas por pulsaciones intro o por caracteres punto y coma(vector columna).
Ejemplo:
Definicin de una matriz A de dimensin (33):
>> A =[1 2 3; 4 5 6; 7 8 9]
A = 1 2 3
4 5 6
7 8 9
Vectores
Para definir vectores solo basta con introducir una fila de numeros, sin punto y coma dentro de los
corchetes.
Ejemplo:
>>a =[1 2 2 3 2 4 3]
Respuesta:
a = 1 2 2 3 2 4 3
Otras formas de introducir vectores

Existen otras dos formas, para vectores donde los elementos individuales estn espaciados linealmente
entre s:
1. Mediante la notacin dos puntos:
(a:b:c): crea un vector que comienza en a, incrementa b y finaliza en c.
2. Mediante la funcin linspace:

linspace(x1,x2,n): genera un vector con n valores igualmente espaciados entre x1 y x2.
21
NOTA: Para espaciado logartmico se utilizari la funcin logspace. (logspace(x1,x2,n)).
Escalares
Igualacin simple. Sin corchetes.
Ejemplo:
>>b =32
Respuesta:
b = 32
Direccionamiento de matrices y vectores
Acceso elemento a elemento

Para acceder a los elementos, se escribe los subndices:
MATRICES: A(i,j)
VECOTRES: a(i)
Siendo,
i: nmero de columnas
j: subndice de las filas
Es importante saber que las matrices se almacenan por columnas (aunque se introduzcan por filas), y
teniendo en cuenta esto puede accederse a cualquier elemento de una matriz con un slo subndice.
Ejemplo:
Si A es una matriz (33) se obtiene el mismo valor escribiendo A(1,2) que escribiendo A(4).
Acceso en bloque
Para acceder a un bloque de elementos a la vez existen dos caminos:
1. La notacin de dos puntos:
a(a:b) nos da los elementos desde la posicin a hasta la posicin b.
Si aadimos un tercer valor, nos maracar el intervalo:
a(a:c:b), obtendriamos el segundo, cuarto y sexto elemento del array.
2. La indicacin de las posiciones de los elementos que queremos obtener:

Ejemplo:
>> a([2 4 1 6])
Respuesta:
ans = 2 3 1 4
o bin,
>>c=a([2 4 1 6])
c = 2 3 1 4
Operaciones de matrices, vectores y escalares

Las operaciones matemticas sencillas entre escalares, matrices y vectores siguen una interpretacin
natural, es decir, se aplica la operacin a todos los elementos.
Elemento-a-elemento
Si tienen la misma longitud y orientacin, la suma, resta, multiplicacin y divisin se aplican sobre la base
de elemento-a-elemento.
Para multiplicar elemento a elemento, escribimos un punto antes del operador deseado.
Ejemplo:
Elevacin al cuadrado elemento-a-elemento de la matriz A:
>> A.^2
ans = 1 4 9
16 25 36
49 64 81
22
Orientacin del vector o matriz

Usando el operador transpuesta ( ) de MatLab, podemos pasar de vector fila a vector columna, y
viceversa.
Ejemplo:
Trasposicin de la matriz A:
>> A'
Respuesta:
ans = 1 4 7
2 5 8
3 6 9
Controles de flujo o bucles

Como en muchos programas de programacin, aparecen los bucles o controles de flujo.
for (desde)
La estructura for es un bloque que ejecuta las rdenes desde i = a hasta i=b, siendo i el subndice del
vectores y/o matrices que intervienen en las rdenes.
for i = a:b
rdenes
end
Funciones
Funcin lenght
Length( ) = Nmero de valores que hay en el vector ( ).
Funcin abs
abs( ) = Valor absoluto de ( ).
Funcin exp
exp( ) = Exponencial de ( ).
Graficar
La orden ms utilizada para graficar es:
plot(X1,Y1,)
Opciones
El grafico puede ser modificado especificando el parmetro, como:
lineWidth: Espesor de la lnea.
color: Color de todas las lneas graficadas.
fontsize: Tamao de letra.
Estilo de lneas, marcadores y colores

Se pueden especificar los colores y estilos de lnea dando un argumento adicional a plot.
SIMBOLO COLOR SIMBOLO ESTILO DE LINEA

y Amarillo . Punto
m Magenta o Crculo
c cien X Marca-x
r Rojo + Ms
g Verde * Estrella
b Azul - Lnea slida
w Blanco : Lnea punteada
K Negro -. Lnea punto-raya
-- Lnea de trazos
23
Adicin de rejillas y etiquetas

- La orden grid on aade una rejilla a la grfica actual en las marcas. Grid off la elimina y grid sin
ningn argumento, las conmuta.
- Los ejes horizontal y vertical se pueden etiquetar, respectivamente, con las rdenes xlabel(texto eje
x) e ylabel(texto eje y).
- La orden title aade una lnea de texto en la parte superior de la grfica.
- Podemos aadir tambin cualquier cadena de texto a cualquier punto en la grfica, especificando las
coordenadas:
text (x,y,texto)
Tambin se puede poner una cadena de texto con el ratn, gtext (texto).
Ejercicio
Representa (plot) la cloronidad en funcin de la concentracin de oxgeno y la temperatura.

Los datos a tener en cuenta sn:
Clorinidad equivalente en KCl de 0 g/L-25g/L, 6 cloronidades linealmente
-
separadas (linspace).
- El rango de temperatura s de 0 a 50 C.
- El agua problema esta en equilibrio con el oxgeno del aire a 1 atmsfera de
presin (Ecuacin [1.1]).
- Eje x: Temperatura, C.
- Eje y: Concentracin de oxgeno disuelto, mg/L.
OPCIONAL: Cuadricula la grfica (grid)
Solucin
%--------------------- SoluOxiTemp.m ----------------------

%
% Solubilidad del Oxgeno disuelto en funcin de
% la temperatura y salinidad en KCl, mg/l.
%
%----------------------------------------------------------
clear all
cl=[linspace(0,25,6)]';
N=length(cl); Ta=-273.15;
T=[linspace(abs(Ta),abs(Ta)+50,150)]';
M=length(T);
for I=1:N
for J=1:M
C(I,J)=-139.344+(1.5757e+5)./T(J)-(6.6423e+7)./...
(T(J).^2)+(1.2438e+10)./(T(J).^3)-...
(8.6219e+11)./(T(J).^4)-cl(I).*...
(3.1929e-2-(19.428/T(J))+3867.3./(T(J).^2));
end
end
C=exp(C);
plot(T+Ta,C,'Color','k','LineWidth',1), grid
xlabel('\fontsize{10}Temperatura, ^oC')
ylabel('\fontsize{10}Conc. oxgeno disuelto, mg/L')
text(16,10.5,'\fontsize{10}Clorinidad = 0 g/L')
text(16,7,'\fontsize{10}Clorinidad = 25 g/L')
%----------------------------------------------------------
24
Grfica
16
Conc. oxgeno disuelto, mg/L 14
12
Clorinidad = 0 g/L
10
Clorinidad = 25 g/L
4
0 5 10 15 20 25 30 35 40 45 50
Temperatura, oC
Figura 1.1.1: Representacin de la cloronidad en funcin de la concentracin de oxgeno y de la
Temperatura
Aclaraciones:
cl=[linspace(0,25,6)] : cl es un vector, de 6 valores, del 0 al 25. Valores equidistantes entre ellos.
cl=cl : El apostrof traspone el vector, y as poder operar con el vector T, tambin traspuesto.
1.5757e+5 : 1.5757 * 10 elevado a +5 (e = 10^5).
: No influyen en el clculo, es totalmente esttico, porque no cabe en una lnea.
plot(T+Ta,C,'Color','k','LineWidth',1) : Grafica temperatura respecto concentracin de oxgeno, en color

negro (k) y con un grosor de 1.
text(6,5,'\fontsize{10}Clorinidad = 25 g/L') : Escribe Clorinidad =25 g/L, a tamao 10 y en la posicin

(6,5).
25
Tema 2 :
AUTODEPURACIN
26
Prctica 2: AUTODEPURACIN

Qumico: 1,5h
TOTAL 4,0h
Resumen
Se introducir el concepto de autodepuracin, sus etapas (desoxigenacin y reoxigenacin) y
se determinarn las constantes (k1 y k2) a partir de:
- Concentracin del oxgeno disuelto de un curso de agua (terico y prctico).
- Concentracin de oxgeno de saturacin.
- DBOu.
- Tiempo.
Se modelizar el oxgeno disuelto respecto el tiempo, y se comparar los resultados prcticos
versus tericos.
Objetivos
de:
Fsico-Qumicos
- Definir el concepto de autodepuracin.
- Identificar las etapas de la autodepuracin (desoxigenacin y reoxigenacin)
- Determinar las constantes de autodepuracin a partir de la ecuacin de
Streeter y Phelps.
- Modelizar un proceso de autodepuracin.
MatLab/Informticos
- Definir vectores y escalares.
- Operar con vectores y escalares.
- Crear y llamar a un archivo *.m.
- Conocer y ejecutar las funciones:
optimset
fminsearch
sum
exp
round
num2str
- Graficar (plot).
hold on
axis
- Aadir etiquetas.
- Dibujar una recta (line).
27
Conocimientos bsicos
Biooxidacin
Proceso natural de disolucin que resulta de la accin de un grupo de bacterias con habilidad de
oxidar minerales sulfurados, permitiendo la liberacin de los valores metlicos contenidos en ellos.
Biodegradacin
Llamada tambin biocorreccin.
Es un proceso en el cual los microorganismos naturales (levaduras, hongos o bacterias)
descomponen o degradan sustancias peligrosas, transformndolas en sustancias menos txicas o
inocuas. Los microorganismos, igual que los seres humanos, comen y digieren sustancias
orgnicas, de las cuales se nutren y obtienen energa. (En trminos qumicos, los compuestos
"orgnicos" son aquellos que contienen tomos de carbono y de hidrgeno).
La biodegradacin puede producirse en presencia de oxgeno (en condiciones aerobias) o sin l

(en condiciones anaerobias).
Los microorganismos descomponen los contaminantes orgnicos en productos inocuos,
principalmente dixido de carbono y agua en el caso de la biodegradacin aerobia. Una vez
degradados los contaminantes, la poblacin de microorganismos disminuye porque ha agotado su
fuente de alimentos. Los microorganismos muertos o una poblacin pequea de microorganismos
sin alimentos no presentan riesgo de contaminacin.
Fundamentos qumicos
Introduccin
El movimiento de las aguas en su discurrir reoxigena las aguas de modo que la carga
contaminante vertida disminuye por biooxidacin y por decantacin, tendiendo la concentracin
de oxgeno disuelto a recuperar la saturacin, a la vez que se degrada la materia microbiana,
denominndose a este proceso autodepuracin.
Autodepuracin de los ros

Los ros tienen una capacidad limitada de absorber y eliminar la contaminacin de los vertidos
que reciben. Esto es debido a la utilizacin de la materia orgnica del vertido por las bacterias
presentes en el ro, de modo que en su actividad mineralizan la materia putrescible7 de los
vertidos, y a su vez las bacterias son fagocitadas8 por los protozoos y stos a su vez por
organismos superiores.
Cuando la capacidad de un ro de eliminar la contaminacin vertida sobre l, autodepuracin,
es rebasada, el ro comienza a presentar signos de contaminacin, como son turbidez, malos
olores
El ndice general de contaminacin utilizado para seguir el proceso de autodepuracin de un ro

es la concentracin de oxgeno disuelto, que disminuir bruscamente en el punto del vertido,
para despus de alcanzar un mnimo gradualmente tender hacia la concentracin
correspondiente a la saturacin.
La evolucin de la concentracin de oxgeno disuelto, a partir de un punto de descarga de un

contaminante, se consigue realizando un balance completo de la concentracin de oxgeno
disuelto.
7
La materia que sufre una descomposicin de protenas animales, especialmente por microorganismos
anaerbicos. El resultado es el de un olor desagradable.
8
Atraeidas por una celula para ser destruidas o digeridas
28
Streeter y Phelps simplificaron el balance a las etapas de desoxigenacin y reoxigenacin.

Constataron que la autodepuracin de un ro se constituye principalmente por dos causas, una
de consumo de oxgeno por los microorganismos, al biooxidar la materia orgnica,
desoxigenacin y otra de aporte de oxgeno del aire al interior del ro, reoxigenacin.
Figura 2.1: Dficit de oxgeno y su concentracin en un ro, despus de un vertido biodegradable.
Desoxigenacin
Al entrar la materia orgnica en el ro, comienza su biodegradacin, proliferando las bacterias
sobre el sustrato orgnico, utilizando el oxgeno disuelto en el agua para su metabolizacin, lo
que da como resultado una disminucin de la concentracin de oxigeno disuelto desde la
saturacin, dficit de oxgeno disuelto:
D = Os O [2.4]
Donde:
D: Dficit de concentracin de oxgeno disuelto.
Os: Concentracin de oxgeno disuelto a saturacin.
O: Concentracin de oxgeno disuelto en el agua.
La velocidad de aumento de dficit de oxgeno es proporcional a la DBO presente:
dD
( ) d = k1 DBO [2.5]
dt
La evolucin de la DBO es una reaccin de primer orden, su evolucin sigue:
DBO(t ) = DBO u e -k1 t [2.6]
Con lo que la variacin de la concentracin de oxgeno disuelto en esta etapa de

desoxigenacin se expresa por:
dD
( ) d = k1 DBOu e -k1 t [2.7]
dt
29
Reoxigenacin
El movimiento del agua en su discurrir provoca que el oxgeno del aire se disuelva en el agua,
influyendo en esta transferencia gas-lquido las condiciones fsicas del curso del agua y su
propia concentracin de oxgeno, de modo que la reoxigenacin, disminucin de dficit de
oxgeno, es contraria al propio dficit de oxgeno:
dD
-( )r = k2 D [2.8]
dt
Dficit de concentracin de oxgeno

Los procesos combinados de reoxigenacin y desoxigenacin conducen a:
dD dD dD
( )= ( )r + ( ) d = k 2 D + k1 DBOu e -k1 t [2.9]
dt dt dt
La ecuacin Diferencial se integra y se obtiene:
k1 DBOu e -k1 t -k1 t k DBOu D (0) (k1 k 2 )

D (t) = - e +( 1 ) e -k 2 t [2.10]
k1 k 2 k1 k 2
Siendo D(0) el dficit de la concentracin del oxgeno disuelto en el momento inicial.
La fiabilidad de la prediccin de la evolucin del oxgeno disuelto a lo largo de un ro, se basa

en la validez de las constantes k1 y k2.
Determinacin de k1
Se obtiene de los datos de la evolucin de la determinacin de la DBO del vertido, o de la
expresin equivalente:
1 DBOa
k1 = log e [2.11]
t ab DBO b
Donde tab es el tiempo del curso del agua entre los puntos a y b, a los que le corresponden
DBOa y DBOb.
Determinacin de k2
Se realiza a travs de medidas experimentales; una de las ecuaciones empleadas es:
v
k2 = C 3
[2.13]
h 2
Donde:
C: Constante.
v: Velocidad media del ro [m/s].
h: Profundidad media del ro [m].
Es importante conocer cundo el dficit de oxgeno ser mximo, y cul es su cuanta, ya que
esta condicin determina dnde estar el tramo de mximo deterioro del ro, por mnima,
concentracin de oxgeno disuelto.
El tiempo en el que el dficit se hace mximo, se obtiene de la derivada y anulacin de la
Ecuacin integrada [2.10]. Con el tiempo del mximo se obtiene el dficit mximo, por su
sustitucin en la misma ecuacin.
30
NOTA: La aplicacin prctica de la ecuacin de Streeter-Phelps, es una aproximacin a la

realidad.
Los factores que provocan desviaciones en su apreciacin son:
- La disminucin de la DBO a consecuencia de la sedimentacin de materias
biodegradables.
- La presencia de algas. Ya que altera la concentracin de oxgeno disuelto,
aumentndolo durante perodos de insolacin, y disminuyndolo en horas de
ausencia de luz.
- La presencia de sedimentos anxicos, que consumen oxgeno disuelto.
- Las constantes k1 y k2 , que varan con la temperatura.
Evolucin de la concentracin de oxgeno causada por un

vertido
La evolucin de la concentracin de oxgeno en el curso de un ro, sometido a un vertido
contaminante, puede seguirse directamente en trminos de la concentracin de oxgeno
disuelto, en lugar del dficit de oxgeno.
La etapa de desoxigenacin se formula como:
dO
( ) d = - k1 DBOu e -k1 t [2.14]
dt
La reoxigenacin:
dO
= k 2 (O s - O) [2.15]
dt r
La ecuacin global de la variacin del oxgeno disuelto es:
dO
= k 2 (O s - O) - k1 DBOu e -k1 t [2.16]
dt
Si integramos:
k1 DBOu -k1 t k (DBOu + O s O(0)) + k 2 (O(0) O s ) -k 2 t

O(t) = O s e +( 1 )e [2.17]
k 2 k1 k 2 k1
Siendo O(0) la concentracin inicial de O.
31
Programa en MatLab
IDEAS PREVIAS
Conceptos generales
Explicado en el Prctica 1.
Clf
Ordena a MatLab que suprima todos los objetos de la figura actual.
Global
Toda variable que est declarada como global, ser compartida por todas las funciones del programa
principal que tambin la declaren global, es decir, cualquier asignacin a esa variable, en cualquier
funcin, est disponible para todas las funciones que la declaran global.
Definicin de matrices, vectores y escalares


Archivos *.m
Se pueden colocar rdenes en un archivo de texto y, a continuacin, decirle a MatLab que lo abra y
evale las rdenes exactamente como si hubiesen sido escritas desde la lnea de orden de MatLab.
Estos archivos, que al fin y al cabo son funciones, se llaman archivos scripts o archivos-M, y deben
finalizar con la extensin m.
Como crearlos
Para crear un archivo-M clicamos New del men File y seleccionamos M-file. Abriendose as, el editor,
dnde se escriben las rdenes que contendr el archivo creada (archivo script).
Como ejecutarlos
Una vez guardado este archivo-M en el disco, MatLab ejecutar las rdenes de dicho archivo cuando el
usuario escriba el nombre del archivo(sin extensin) en la lnea de orden.
Las rdenes dentro del archivo-M tienen acceso a todas las variables en el espacio de trabajo de
MatLab, y todas las variables creadas en el archivo-M se hacen parte del espacio de trabajo.
Como guiar a MatLab para encontrarlos

La ventana Current Directory muestra los ficheros del directorio activo o actual, es decir, la direccin
donde se encuentran los archivos con los que trabaja MatLab, en dicha sesin.
El directorio activo se puede cambiar desde el Current Directory o desde la barra de herramientas,
(debajo de la barra de mens) con los mtodos de navegacin de directorios propios de Windows.
(Figura 3)
32
Figura 3. Bsqueda del directorio deseado
2
NOTA: En el caso de que quisieramos crear una funcin de pocas lneas, como es 2x +7, podriamos
hacerlo mediante la Function-Handel.
>> y=@(x) 2*x.^2+7;
Funciones
Length
Optimset
Crear o corregir una estructura de las opciones de la optimizacin.
>>opciones = optimset('parmetro1',valor1,'parmetro2',valor2,...)
Las opciones que se podrian canviar con la funcin optimset son:
TolFun y TolX: tolerancia en la funcin y en el eje x
MaxFunEvals: nmero mximo de evaluaciones
Fminsearh
El archivo fminsearch.m encuentra el mnimo de una funcin escalar de varias variables (fun)
comenzando en una estimacin inicial (x0) y con unas determinadas opciones.
>>x = fminsearch(funcin,x0,opciones)
Num2sr
Funcin que convierte nmeros a sus representaciones de la secuencia.
Round
Por defecto redondea al entero ms cercano, pero puedes introducir cuantos decimales quieres.
n n
round(10 *valor)/ 10 : Redondear el valor con n decimales
Sum
Funcin sumatorio.
Graficar
33
Funcin Hold
Permite modificaciones de un grfico ya definido.
Hold on (aplica las rdenes en el grfico existente)
Hold off (formatea el grfico existente, y muestra las nuevas rdenes en la misma figura).
Funcin line
Crea una linea en el grfico en la posicin (x, y) con las propiedades deseadas:
>>line([X1 X2],[Y1 Y2],'Nombre de la propiedad a canviar',Valor de la
propiedad,...)
Ejes a medida
axis([xmin, xmax, ymin e ymax]) fija los valores mximo y mnimo de los ejes usando los valores
dados en el vector fila.
Ejercicio
Mediante el anlisis de un ro conseguimos los siguientes datos:

- una DBOu de 20mg/L,
- una concentracin de oxgeno disuelto a saturacin de 8,8mg/L,
- la concentracin de oxgeno disuelto en el agua respecto al tiempo:
t 0 0,4 0,8 1,2 2 3 5 10

O 7,8 5,8 5,3 3,8 5,6 5,8 6,8 7,8
- k1 y k2 se estiman cuando la concentracin de oxgeno disuelto es mnima

(fminsearch). Con una tolerancia en la funcin y en el eje x de 10^(-20) y un
nmero mximo de evaluaciones de 10000. Los puntos iniciales de la
estimacin sn 0,4 y 0,9, respectivamente.
A partir de los datos y sabiendo que el sistema estudiado cumple la ecuacin de Streeter-
Phelps, calcular las constantes de autodepuracin ( k1 y k2) y representar la concentracin de
oxgeno disuelto en el agua problema respecto al tiempo.
Solucin
%---------------- Autodepur.m -----------------------------

clear all, clf, global t O DBOu O0 Os
DBOu=20; Os=8.8;
t=[0 .4 .8 1.2 2 3 5 10]';
O=[7.8 5.8 5.3 3.8 5.6 5.8 6.8 7.8]'; O0=O(1);
Opciones=optimset('TolX',1e-20,'TolFun',...
1e-20,'MaxFunEvals',10000)
Est=fminsearch('AutodepurMin',[0.4 0.9],Opciones);
k1=Est(1); k2=Est(2);
tc=linspace(t(1),t(length(t)),100);
Oc=Os-k1*DBOu/(k2-k1)*exp(-k1*tc)+(k1*...
(Os+DBOu-O0)+k2*(O0-Os))/...
(k2-k1)*exp(-k2*tc);
plot(tc,Oc,'k-','LineWidth',1), hold on, plot(t,O,'ko')
axis([0 10 3.5 9])
ylabel('\fontsize{13} Conc. de oxgeno (mg/L)')
xlabel('\fontsize{13}Tiempo (das)')
text(3,8.65,'\fontsize{13}Conc. de oxgeno a saturacin')
34
line([t(1) t(length(t))],[Os Os],'Color','k','LineWidth',1)

Text1=['\fontsize{13}k_1= ',num2str(round(1000*k1)/1000)];
Text2=['\fontsize{13}k_2= ',num2str(round(100*k2)/100)];
text(7,4.5,Text1), text(7,4,Text2)
%----------------------------------------------------------
%------ Funcin llamada por Autodepur.m -------------

function q=AutodepurMin(p);
global t O DBOu O0 Os
k1=p(1); k2=p(2);
Oc=Os-k1*DBOu/(k2-k1)*exp(-k1*t)+(k1*(Os+DBOu-O0)...
+k2*(O0-Os))/(k2-k1)*exp(-k2*t);
q=sum((O-Oc).^2);
%----------------------------------------------------
Grfica
8.5
Conc. de oxgeno a saturacin
8
Conc. de oxgeno (mg/L)
7.5
6.5
5.5
4.5 k1= 0.358

4 k2= 1.04
3.5
0 1 2 3 4 5 6 7 8 9 10
Tiempo (das)
Figura 2.5 Comparacin de la concentracin de oxgeno respecto al tiempo (terica y
experimentalmente).
Aclaraciones:
plot(tc,Oc,'k-','LineWidth',1), hold on, plot(t,O,'ko') : Grfica tc respecto Oc, en color negro, con grosor
de 1, y en la misma grfica (hold on) segunda representacin de t respecto a O,
per con puntos ('ko').
line([t(1) t(length(t))],[Os Os],'Color','k','LineWidth',1) : Dibuja una lnea desde t(1) hasta t(length(t)) en el
eje x, en Os del eje y, en color negro y con grosor de 1.
Text1=['\fontsize{13}k_1= ',num2str(round(1000*k1)/1000)] : El texto 1 est definido por un tamao de

13, y su contenido s k1 = redondeo de k1 con tres
decimales.
35
Tema 3 :
FRMULA EMPRICA
36
Prctica 3: DETERMINACIN DE LA FRMULA

EMPRICA DE LOS COMPONENTES DE LAS
AGUAS

Qumico: 0,5h
TOTAL 2,0h
Resumen
Se calcular la frmula emprica de los microorganismos y sustratos de las aguas,
representadas por frmulas qumicas, (CaHbOcNd), obtenidas de los anlisis qumicos de los
componentes de las aguas.
Dicho calculo se realiza con la finalidad de facilitar el seguimiento y evaluacin de las
reacciones bioqumicas de biodegradacin.
Objetivos
de:
Fsico-Qumicos
- Clcular frmula emprica mediante balances estequiomtricos, teniendo en
cuenta:
la reaccin de oxidacin
la DQO
los SV
el COT
el nitrogeno orgnico
MatLab/Informticos
double
round
min
- Resolver sistema de ecuaciones (solve).
Construir objetos simblicos (syms).
37
Demanda Qumica de Oxigeno (DQO)
Parmetro que mide la cantidad de materia orgnica susceptible de ser oxidada por medios
qumicos que hay en una muestra lquida. Se utiliza para medir el grado de contaminacin y se
expresa en mg O2/litro.
Sin embargo, se ha de saber que puede haber interferencias debido a que haya sustancias
inorgnicas susceptibles de ser oxidadas (sulfuros, sulfitos, yoduros, etc.).
Es un mtodo aplicable en aguas continentales (ros, lagos, acuferos, etc.), aguas residuales o
cualquier agua que pueda contener una cantidad apreciable de materia orgnica. No es aplicable
para las aguas potables debido al valor tan bajo que se obtendra y, en este caso, se utiliza el
mtodo de oxidabilidad9 con permanganato potsico.
Carbono orgnico total (COT)

Se mide mediante la oxidacin de la materia orgnica mediante calor y oxgeno o mediante
oxidantes qumicos y se detecta mediante anlisis de infrarrojo la produccin de CO2.
Fundamentos qumicos
Introduccin
Los microorganismos de un agua que cumple las condiciones de la medida de la DBO,
metabolizan las materias orgnicas biodegradables de las muestras a examinar, de
procedencia urbana o industrial, multiplicndose, a la vez que consumen las sustancias
biodegradables.
Las reacciones bioqumicas de esta actividad metablica se resumen en reacciones de sntesis
de nuevos organismos, reacciones de produccin de energa, para desarrollo de su actividad y
reacciones de degradacin de los propios microorganismos, consumiendo estas tres
reacciones oxgeno, hasta que el sustrato disponible se agota, comenzando entonces la fase
de metabolismo endgeno10, caracterizado por un consumo mnimo de oxgeno.
La DBO mide el consumo de oxgeno en una muestra causado por las reacciones indicadas.
Clculo de la frmula emprica del agua residual

Las reacciones bioqumicas de biodegradacin son ms fciles de seguir y evaluar, si a los
microorganismos y al sustrato de las aguas a tratar se les representan por frmulas qumicas,
obtenidas de los anlisis qumicos de los componentes de las aguas.
El procedimiento de clculo de la formula qumica que representa un agua residual, se obtiene
por balances estequiomtricos de las reacciones implicadas en la determinacin de los
componentes del agua, que a modo de problema seria:
Calcular la frmula emprica del sustrato de un agua residual, compuesto de C, H, O y N, cuya

demanda qumica de oxgeno, DQO, es 555 mg/L, los slidos voltiles, SV, 275 mg/L, el
carbono orgnico total, COT, 165 mg/L, y el contenido de nitrgeno Kjeldhal 19mg/L.
Aclaraciones:
a) El clculo se basa en la reaccin qumica de la oxidacin total de la materia
carbonada de la muestra, pero no de la nitrogenada, tal como se lleva a cabo
en la determinacin de la DQO.
9
Capacidad para que un tomo o un ion pierda electrones.
10
Que se forma o engendra en el interior de algo, como la clula que se forma en el interior de otra.
38
b) El nitrgeno y el carbono se determinan a travs del nitrgeno Kjeldahl y del

carbono orgnico total.
c) El oxgeno utilizado es la DQO.
d) La cantidad de materia orgnica de la muestra est representada por los SV.
a) La reaccin de oxidacin es:
C a H b Oc N d + xO2
aCO2 + yH 2 O + dNH 3 [3.1]
b) Basndose en el COT, se tiene:
COT
a= [3.2]
12
Del mismo modo para el nitrgeno:
N
d= [3.3]
14
c) Para la DQO, se tiene:
DQO
x=
32 [3.4]
d) Para los SV se tiene:
a 12 + b + c 16 + d 14 = SV [3.5]
El balance de los hidrgenos conduce a:
b = 2 y + 3 d [3.6]
Y con el balance de los oxgenos se obtiene:
c + 2 x = 2a + y [3.7]
Con estas ecuaciones el problema queda resuelto, obteniendo que el agua residual responde a
C13,8 H26,11 O4,02 N1,36, frmula que se simplifica dividiendo los subndices de los elementos por
el menor; y posteriormente redondeando, C10 H20 O3 N.
39
Programa en MatLab
IDEAS PREVIAS
Conceptos generales
Explicado en el Prctica 1 y 2.
Definicin matrices, vectores y escalares

Operaciones con matrices, vectores y escalares

Funciones
Double
Transforma fracciones a nmeros reales.
double ( ) : Transforma el valor ( ).
Round
Min
min ( ) : Busca el valor mnimo de ( ).
Resolucin de ecuaciones algebraicas

Para resolver sistemas de ecuaciones se ha de seguir estos pasos:
2. Introducir la orden syms, que nos permitira construir objetos simblicos (las
incognitas).
>>syms a b c n
siendo,
[a b c n] las variables que representan las incognitas del sistema.
3. Llamar a la funcin solve, quin identifica las incognitas del sistema.

>> [a b c .n] = solve (ecuacin a, ecuacin b, ecuacin c,
.ecuacin n)
40
Ejercicio
Calcular la frmula emprica del sustrato de un agua residual, compuesto de C, H, O y N, cuya

demanda qumica de oxgeno, DQO, es 555 mg/L, los slidos voltiles, SV, 275 mg/L, el
carbono orgnico total, COT, 165 mg/L, y el contenido de nitrgeno Kjeldhal 19mg/L.
Solucin
%----------------------EcuFormula.m-----------------------------
%Calculo de la formula qumica de un agua residual, conociendo
%la DQO, los SSV, el COT, y el nitrogeno orgnico.
%---------------------------------------------------------------
DQO=555; SV=275; COT=165; N=19;
A=COT/12; D=N/14; X=DQO/32;
syms b c y
[b,c,y]=solve(12*A+b+16*c+14*D-SV,2*y+3*D-b,c+2*X-2*A-y);
B=double(b); C=double(c);
z=[A B C D];
z=round(z/min(z))
%---------------------------------------------------------------
Aclaraciones:
B=double(b) : Transforma la fraccin b, a un nmero real B.
round(z/min(z)) : Redondea la fraccin al entero ms cercano.
41
Tema 4 :
DEMANDA BIOQUMICA DE
OXGENO
42
Prctica 4: PARMETROS CINTICOS DE LA DBO

Qumico: 1,5h
TOTAL 4,0h
Resumen
Se tratarn los conceptos bsicos relacionados con la DBO y sus parmetros cinticos.
Representando luego la DBO en funcin del tiempo y comparando datos experimentales con
tericos.
Objetivos
de:
Fsico-Qumicos
- Definir y ver la importancia de la DBO.
- Conocer los conceptos relacionados con la DBO (DOU, DON, DBOu y DBO5).
- Reconocer las reacciones bioqumicas que ocurren en las aguas.
- Conocer los parmetros cinticos de la DBO.
- Modelizar la evolucin de la DBO respecto al tiempo.
- Comparar resultados tericos versus experimentales.
MatLab/Informticos
- Definir matrices, vectores y escalares.
length
sum
diff
log
exp
round
num2str
- Crear bucles, controles de flujo (for).
- Graficar (plot)
Ms de un grafico en una sola figura.(subplot)
axis
- Aadir etiquetas.
43
Demana Biologica de Oxgeno
Tambin denominada demanda bioqumica de oxgeno, (DBO).
Es el parmetro que mide la cantidad de materia susceptible de ser consumida u oxidada por
medios biolgicos que contiene una muestra lquida, y se utiliza para determinar su grado de
contaminacin. Normalmente se mide transcurridos 5 das (DBO5) y se expresa en mg O2/litro.
Fundamentos qumicos
Disponiendo de frmulas qumicas, se pueden representar las reacciones bioqumicas que
ocurren en las aguas, como es la sntesis de hetertrofos que se representa con las frmulas
ms sencillas:
8CH 2 O + 3O2 + NH 3 C 5 H 7 NO2 + 3CO2 + 6 H 2 O (3.1)
La reaccin de produccin de energa, consumiendo el sustrato:
8CH 2O + O2 CO2 + H 2O + 40kcal / mol (3.2)
La reaccin de degradacin de la biomasa:
C 5 H 7 NO2 + 5O2 5CO2 + 2 H 2 O + NH 3 (3.3)
Conceptos relacionados con la DBO

Las condiciones de medida, fueron establecidas en 1912 por una UK Royal Commission, en su
octavo informe titulado Standards and Tests for Sewage and Sewage Effluents Discharging into
Rivers and Streams. En este informe se recomend un perodo de incubacin de cinco das,
tiempo medio de los ros ingleses en llegar al mar, a la temperatura de 65F, mxima
temperatura de los ros ingleses en el mes ms clido, que equivale a 18,3C, redondeado a
20C.
El test qued inalterado desde 1912 y se conoce universalmente como DBO5 a 20C, o
simplemente DBO. De la evolucin temporal de la DBO ejercida, oxgeno absorbido frente al
tiempo, se obtiene la DBO5 y parmetros cinticos por ajuste de los datos experimentales.
La evolucin de la DBO con el tiempo, en la medida de la muestra, debidamente neutralizada y

con nutrientes minerales si fuese necesario, sigue la curva continua, tal como se indica en la
Figura 3.1, alcanzndose en unos 20 a 40 das el valor lmite de la Demanda de Oxgeno
ltima, DOU, incluida la Demanda de Oxgeno del Nitrgeno, DON, causado por nitrificacin. El
valor lmite de la Demanda Bioqumica de Oxgeno ltima, DBOu, corresponde a la biooxidacin
total de la materia carboncea.
En la medida de la DBO se excluye la DON por adicin en la muestra de un inhibidor de
crecimiento de las bacterias nitrificantes, nitrosomonas y nitrobacter.
Las aguas industriales, ya desprovistas de microorganismos o con productos xenobiticos,

ajenos a la naturaleza, necesitan un perodo inicial de desarrollo para adaptacin de los
microorganismos, para habituarse, por sucesivas generaciones en contacto a estas sustancias
extraas, siguiendo una trayectoria como la indicada en la Figura 3.1, en la lnea a trazos,
anloga a la lnea continua sometida a un retardo.
44
Figura 3.1: Evolucin de la demanda bioqumica de oxgeno
La reaccin de biodegradacin, para la medida de la DBO, en las condiciones de medida exige:
- pH neturo.
- Inhibidor de nitrificacin.
- Absorcin del CO2 formado en los mtodos manomtricos.
- Inculo inicial mnimo de 104 bacterias/mL.
- Presencia de sales minerales.
- Incubacin en oscuridad.
- Agitacin.
La importancia de la medida de la DBO, radica en que del anlisis de la evolucin de la DBO,

se obtiene informacin de los parmetros cinticos de la biodegradabilidad de la muestra, base
del diseo y operacin de las depuradoras de aguas residuales.
Cintica de la DBO
La variacin temporal de la DBO, en las condiciones de la medida estandarizada, oxgeno
utilizado en las Reacciones (3.1), (3.2) y (3.3), se considera de primer grado, con tendencia a
que la DBO alcance a tiempo infinito la DBOu, que se expresa como:
dDBO(t )
= K (DBOu DBO(t ) ) [3.1]
dt
Siendo K una variable que depende de la temperatura, con funcionalidad tipo Arrhenius.
Como la determinacin se realiza a 20C, K toma un valor fijo. La constante cintica K, indica la
facilidad de la biodegradacin.
Para sustancias fcilmente biodegradables, la DBOu se alcanza en tiempos cortos, de modo
que se pueden caracterizar las aguas residuales por el valor de su K.
La Ecuacin [3.1] integrada conduce a la expresin ms utilizada de la DBO:
(
DBO(t ) = DBOu 1 e K t ) [3.2]
45
Esta expresin general se ajusta bien para la mayora de las aguas residuales.
Determinacin de las constantes cinticas de la DBO

El mtodo ms sencillo para determinar las constantes cinticas de la Ecuacin [3.2], se
obtiene por su transformacin a una ecuacin lineal, y posterior ajuste por mnimos cuadrados
de la recta resultante.
dDBO(t )
= K DBOu e K t [3.3]
dt
Tomando logaritmos neperianos se consigue:
dDBO(t )
log e = log e (K DBOu ) K t [3.4]
dt
Se realizan las siguientes transformaciones de variables:
dDBO(t )
Y = log e
dt
b = log e (K DBOu )
m = K , X =t
Con lo que se obtiene la ecuacin de una recta:
Y = m X +b [3.5]
Las constantes de ajuste de la recta, se obtienen por las ecuaciones:
b=
Y X X X Y
2
[3.6]
N X ( X ) 2 2
N Y X Y X
m= [3.7]
N X 2 ( X )
2
El grado de bondad del ajuste se analiza mediante el coeficiente de correlacin:
N Y X X Y
r=
(N X )( )
[3.8]
( X ) N Y 2 ( Y )
2 2 2
Una vez conocidos m y b, se obtienen los parmetros cinticos de la DBO, deshaciendo los
cambios hechos, y la DBO5:
K = m , DBOu =
eb
K'
, (
DBO5 = DBOu 1 e 5K )
46
Programa en MatLab
IDEAS PREVIAS
Conceptos generales
Explicado en Prctica 1 y 2

Explicado en Prctica 1
Manipulacin matricial
Utilizar el smbolo dos puntos como la designacin de filas o columnas implica, respectivamente,
todas las filas o columnas; por ejemplo, A (:,1) representa todas las filas en la columna uno.
Fijar las filas o columnas de una matriz igual a la matriz vaca [ ] elimina estas filas o columnas.
Funciones
Diff
Clculo diferencial.
diff( ) : Calcula las diferencias entre los elementos que forman el vector ( ).
Log
Calcula el logaritmo del valor deseado.
Length
Explicado en el Prctica 1
Sum
Exp
Round
Num2str

Graficar
Subplot
La orden subplot permite graficar ms de un grfico en una figura.
Subplot (abc) : Divide la figura en a divisiones horizontalmente, b verticalmente y se representa en
la posicin c.
47
Ejercicio
Calcular las constantes cinticas de la DBO (K, DBOu y DBO5) a partir de la variacin de la DBO
frente al tiempo.
Y representar la prueba de bondad de la Ecuacin [3.2] linealizada y los valores experimentales
(con crculos, o), y los calculados (lnea continua).
NOTA: Estas dos graficas, representarlas en una nica figura (subplot).
Las lecturas experimentales de un medidor de la DBO manomtrico, a una temperatura de 20C,

para que K tome un valor fijo se muestran en la siguiente tabla:
Fecha Hora DBO

17 17:30 0
18 9:2 125
18 13:02 150
18 17:30 175
19 9:05 260
19 13:17 275
19 14:15 280
19 18:40 300
20 9:00 340
20 14:00 350
20 17:44 360
21 10:15 380
22 10:55 408
23 9:20 420
23 18:00 424
Solucin
%------------------------- DBOCalc.m --------------------------

clear all, clf
datos=[17 17 30 0
18 09 23 128
18 13 02 150
18 17 30 175
19 09 05 260
19 13 17 275
19 14 15 280
19 18 40 300
20 09 00 340
20 14 00 350
20 17 44 360
21 10 15 380
22 10 55 408
23 09 20 419
23 18 00 424];
d=datos(:,1);
h=datos(:,2);
m=datos(:,3);
DBO=datos(:,4);
t=d+h/24+m/(24*60);
t=t-t(1);
dDBO=diff(DBO);
dt=diff(t);
Y=log(dDBO./dt);
N=length(t);
for i=1:N-1,
X(i)=(t(i)+t(i+1))/2;
end
48
X=X';
N=N-1;
SX=sum(X);
SX2=sum(X.^2);
SY=sum(Y);
SY2=sum(Y.^2);
SXY=sum(X.*Y);
b=(SY*SX2-SX*SXY)/(N*SX2-SX^2);
m=(N*SXY-SX*SY)/(N*SX2-SX^2);
r=(N*SXY-SX*SY)/((N*SX2-SX^2)*(N*SY2-SY^2))^.5;
Xc=linspace(X(1),X(length(X)),200);
Yc=m*Xc+b;
K=-m;
DBOu=exp(b)/K;
DBO5=DBOu*(1-exp(-K*5));
%--------------------------------------------------------------
subplot(211)
plot(X,Y,'o',Xc,Yc,'-')
title('\fontsize{12}Prueba de DBO_5=DBO_u(1-e^{-Kt})')
TEXTO=['r = ',num2str(abs(r))];
text(0.4,2.7,TEXTO),
ylabel('log_e dDBO/dt')
%--------------------------------------------------------------
tc=linspace(t(1),t(length(t)),200);
DBOc=DBOu*(1-exp(-K*tc));
subplot(212)
plot(t,DBO,'o',tc,DBOc,'-')
X=['\fontsize{12}Resultados y comparacin de datos',...
' calculados y experimentales'];
title(X)
Texto1=['K = ',num2str(round(K*100000)/100000)];
text(4.75,150,Texto1)
Texto2=['DBO_u = ',num2str(round(DBOu)),' mg/L'];
Texto3=['DBO_5 = ',num2str(round(DBO5)),' mg/L'];
ylabel('DBO (mg/L)'), xlabel('Tiempo (das)')
axis([t(1) t(length(t)) 0 DBO5+20])
%--------------------------------------------------------------
49
Grfica
Prueba de DBO5=DBOu(1-e-Kt)
6
loge dDBO/dt 5
3
r = 0.97421
2
0 1 2 3 4 5 6
Resultados y comparacin de datos calculados y experimentales

400
300
DBO (mg/L)
200
K= 0.52825
100 DBOu = 454 mg/L
DBO5 = 422 mg/L
0
0 1 2 3 4 5 6
Tiempo (das)
Figura 3.2: Evolucin y resultados del ajuste numrico de la DBO.
Aclaraciones:
d=datos(:,1) : d ser una matriz con todos los valores de la primera columna de la matriz datos.
t=d+h/24+m/(24*60) : Pasa todos los valores referidos a tiempo en das.
t=t-t(1) : Reinicia el tiempo.
SX=sum(X) : Suma de todas las X.
Y.^2 : Eleva miembro a miembro de Y al cuadrado.
Subplot (211) : En una figura hace dos divisiones horizontales y 1 verticlal, y est definiendo la primera
grfica.
num2str : Funcin que convierte en nmeros a sus representaciones de la secuencia.
Texto1=['K = ',num2str(round(K*100000)/100000)] : El texto 1 s K = nmero redondeado de k con 5

decimales.
text(0.4,2.7,TEXTO) : Escribe el TEXTO, anteriormente definido, en la coordenada (0.4, 2.7).
50
Tema 5 :
DEMANDA QUMICA DE
OXIGENO
51
Prctica 5: DEMANDA QUMICA DE OXGENO

Qumico: 2h
Informtico/MatLab: 2h
TOTAL 4,0h
Resumen
Se conocer la definicin de la DQO y su medida, teniendo en cuenta las interferencias. Para
ello aparecer el concepto de absortividad y el procedimiento para determinarla.
Se modelizar la recta de calibracin Absorbancia/Concentracin.
Objetivos
de:
Fsico-Qumicos
- Conocer el concepto de:
DQO
Absorvancia
Patrn primario
- Reconocer las diferentes reacciones, que se producen al querer determinar la
DQO segn el procedimiento seguido.
- Saber determinar la absorbancia (procedimiento).
MatLab/Informticos
- Definir matrices, vectores y escalares.
- Conocer y utilizar las funciones:
lenght
sum
input
round
- Graficar.
- Modificar las opciones del grfico.
52
Blanco
Un blanco es aquella muestra que contiene todos los ingredientes de la muestra problema, excepto
el elemento a analizar sirve para calibrar, es decir, sirve de valor de referencia.
Valor muestra problema Valor blanco = Valor elemento a analizar
Patrn primario
Sustancias qumicamente puras, que nicamente se contaminan con agua o anhdrido carbnico
del aire, y que por calentamiento lo desprenden sin descomponerse, obtenindose as de nuevo el
producto qumicamente puro tal como indique su formula.
Contiene concentraciones conocidas y sirve de modelo para medir y valorar otra substncia.
Fundamentos qumicos
Las limitaciones y debilidad de la oxidacin bioqumica llev a determinar la demanda de
oxgeno por medios qumicos.
La DQO, Demanda Qumica de Oxgeno, es el ndice general de contaminacin ms usado.
(Chemical Oxygen Demand, COD).
Para aguas blancas, alternativamente se utiliza en lugar de la demanda qumica de oxgeno, la
oxidabilidad al permanganato, PV, Permanganate Value.
Para sustancias fcilmente oxidables, la DQO DBOu, no siendo as para sustancias ms

resistentes, pero siempre se cumplir que DQO DBOu.
El anlisis de la DQO usa la oxidacin qumica, para efectuar la misma reaccin que
provocaban los microorganismos con la materia orgnica, pero ahora de manera enrgica. El
carbono es oxidado a anhdrido carbnico, CO2, permaneciendo el nitrgeno amino amoniacal
en su mismo grado de oxidacin y el nitrgeno correspondiente a los nitritos se oxida a nitratos.
Medida experimental de la DQO

El test estndar se realiza en matraces de 250 mL, provistos de un refrigerante, para llevar la
reaccin a reflujo. En dichos matraces se ponen 50 mL del agua problema y el reactivo
oxidante.
Desde hace unos aos, se viene utilizando este anlisis a escala semimicro, tomando muestras
de 2 a 3 mL de agua problema, realizndose la reaccin en tubos de ensayo provistos de tapn
roscado, de modo que la digestin se pueda realizar bajo presin a 150C durante dos horas.
Este mtodo ha sido aprobado por la Agencia de Proteccin Ambiental Americana,
Environmental Protection Agency, EPA o USEPA.
Tambin se han desarrollado mtodos rpidos de determinacin de la DQO, efectundose la
digestin a 165C durante cinco minutos.
En estos mtodos, el oxgeno se mide por medio del consumo del oxidante utilizado, ya por
valoracin directa del dicromato remanente11 o por espectrofotometra, midiendo la extincin de
su color, amarillo naranja, o por la generacin de color verde azulado, debido al Cr+3 formado.
11
Sobrante.
53
Definicin de la DQO
La DQO es el consumo de oxgeno, mg/L, en la oxidacin total por va hmeda de la materia
carboncea, presente en la muestra problema.
Las condiciones de la oxidacin se deben a la accin combinada de
- un oxidante fuerte, dicromato potsico (K2CrO)
- medio sulfrico. Carbonizador de la materia orgnica
- a temperatura elevada
- tiempo suficiente para completar la oxidacin.
Demanda qumica de oxgeno de sustancias puras

La DQO se realiza oxidando con dicromato potsico en medio sulfrico fuertemente cido. En
presencia de un reductor, que utilice el oxgeno liberado por el dicromato.
La reaccin completa del dicromato en medio sulfrico se representa por:
2 K 2 Cr2 O7 + 8H 2 SO4 2Cr2 (SO4 )3 + 2 K 2 SO4 + 8H 2 O + 3O2 (5.1)
En la reaccin podemos leer que por cada 2 moles de dicromato que se reducen se generan 3
moles de oxgeno, para oxidar la materia orgnica.
Calclese la DQO de un representante de azcares, de una disolucin de glucosa de 1g/L.

La reaccin de oxidacin de la glucosa es,
C 6 H 12 O6 + 6O2 6CO2 + 6 H 2 O (5.2)
6 32
Un gramo de glucosa demanda = 1,066 gramos de oxgeno, lo que significa DQO =
180
1066 mg/L..
NOTA: En una disolucin que contenga 1 g/L de glucosa tiene una DQO de 1066 mg/L.
Calclese la DQO de un representante de sustancias proticas, disolucin de cido glutmico

de 1g/L.
La reaccin de la oxidacin del cido glutmico es:
9
C 5 H 9 O4 N + O2 5CO2 + 3H 2 O + NH 3 (5.3)
2
9 16
La DQO de 1g/L de cido glutmico es = 0,98 g/L = 980 mg/L.
147
NOTA: Una disolucin que contenga 1g/L de carbohidratos o protenas (cido glutmico)
equivale a 1g/L de DQO.
54
Patrn para la medida de la DQO

En la seccin anterior se calcula la DQO de una disolucin de 1g/L de glucosa, resultando ser
1066 mg/L. Pero ocurre que la glucosa obtenida con extraordinaria pureza, al ponerse en
contacto con el aire absorbe agua, que por calentamiento lo desprende, a la vez que se
descompone caramelizndose. De modo que un frasco de glucosa abierto varias veces
absorbe agua en cantidad indeterminada, impidiendo conseguir cantidades exactas de este
producto, y algo anlogo ocurre con el cido glutmico.
Las sustancias qumicamente puras, que nicamente se contaminan con agua o anhdrido
carbnico del aire, y que por calentamiento lo desprenden sin descomponerse, obtenindose
as de nuevo el producto qumicamente puro tal como indique su formula, reciben la
denominacin de <<patrones primarios>>, siendo estas sustancias las adecuadas para
preparar disoluciones de concentracin conocida.
El patrn primario para la determinacin de DQO es el hidrgeno-ftalato potsico, KHC8O4H4,

porque calentado a 120C el producto analtico responde a su frmula, estando exento de
impurezas.
La reaccin de oxidacin del hidrgeno-ftalato potsico se expresa como:
2 KHC8 O4 H 4 + 15O2 16CO2 + 5 H 2 O + K 2 O (5.4)
Esta reaccin establece la base para el clculo de disoluciones de DQO patrn.
Calcular la cantidad de hidrgeno-ftalato potsico necesario para preparar, 1 litro de disolucin

que contenga 500 mg/L de DQO.
(peso molecular: hidrgeno-ftalato potsico = 204,2 y oxgeno = 32)
2 204,2
500 = 425,4 mg KHC8 O 4 H 4
15 32 L
Una disolucin de 425,2 mg/L de hidrgeno-ftalato potsico tiene una DQO de 500 mg/L.
Interferencias en la determinacin de la DQO

Las interferencias ms frecuentes son las que provocan las materias no carbonceas oxidables
por el dicromato, como son Fe+2 y los cloruros, ya que consumen dicromato durante la
oxidacin de la muestra.
Sales ferrosas
La reaccin del consumo de dicromato a causa de la presencia de sales ferrosas es:
2
Cr2 O7 + 6 Fe 2+ + 14 H + 2Cr 3+ + 6 Fe 3+ + 7 H 2 O (5.5)
Es posible eliminar la interferencia, por tratamiento previo de la muestra con un oxidante dbil,
tal como el agua oxigenada. Esta interferencia es poco significativa en aguas naturales.
Cloruros
La reaccin de oxidacin de los cloruros es:
2
Cr2 O7 + 6Cl + 14 H + 2Cr 3+ + 3Cl 2 + 7 H 2 O (5.6)
55
Esta interferencia se evita por adicin de sulfato mercrico (HgSO4), segn la reaccin:
Hg 2+ + 2Cl HgCl 2 (5.7)
La reaccin de la DQO exige un catalizador de la destruccin de la materia orgnica, siendo el

sulfato de plata (Ag2SO4) el ms utilizado.
NOTA: Si la adicin de sulfato mercrico fue insuficiente para precipitar los cloruros, stos
precipitarn la plata, que podra ser insuficiente, para favorecer la destruccin de la materia
orgnica:
Ag + + Cl AgCl (5.8)
Procedimiento y reactivos de la medida de la DQO

Las aguas residuales urbanas brutas presentan una DQO que oscila entre 200 y 1000 mg/L,
mientras que las aguas residuales depuradas no deberan superar los 100 mg/L.
Por tanto se procede a adecuar las determinaciones de la DQO en dos rangos de medida:
- 0 - 150mg/L para aguas de ros poco contaminados y aguas residuales
regeneradas
- 150 - 1500 mg/L para aguas residuales.
Para el rango de 150 a 1500
Preparacin del reactivo de dicromato

Pasos a seguir:
1. Disolver 8,166 g de dicromato, K2Cr2O7, en 200 mL de agua destilada, en un

vaso de precipitados de 1 L.
2. Disolver 10,0 g de sulfato de plata, Ag2SO4, en 600 mL de cido sulfrico,
H2SO4, concentrado, en un vaso de precipitados de 2 L.
3. Mezclar cuidadosamente las dos disoluciones preparadas anteriormente, en un
vaso de precipitados de 2 L (aadiendo la disolucin cida sobre la disolucin
de dicromato).
Como la dilucin del sulfrico es fuertemente exotrmica, se dejar enfriar la
disolucin entre adiciones, para evitar ebullicin con proyeccin del contenido.
4. Llevar la mezcla obtenida a un matraz aforado de 1 L y enrasar con agua.
En esta operacin se vuelve a generar calor, por lo que se proceder a enfriar
el matraz, y el enrase final se realizar con la disolucin a 25C.
Tcnica operatoria para la determinacin de la DQO

Se lleva a cabo en tubos preparados con los reactivos, a los que se le aade el problema o
patrones, segn el siguiente procedimiento:
56
Poner 3 mL de reactivo de dicromato en tubos especiales para esta

determinacin.
Aadir una dosis de HgSO4 en cada tubo, obtenindose los tubos preparados.
Aadir 2 mL de agua destilada a un tubo preparado. Ser el blanco de la

determinacin.
Aadir 2 mL de disolucin de patrones en algunos tubos preparados, para

calibracin.
Poner 2 mL de las aguas a examinar en otros tubos preparados.
Cerrar y homogeneizar, perfectamente, los tubos, y llevar a un bloque

digestor a 150C durante 2 horas.
Dejar enfriar los tubos en nevera un tiempo de diez minutos.
Analizarlos con un espectrofotmetro, y as poder conocer la absorbancia a

620 nm.
Para el rango de 0 a 150

Se realiza la operatoria anteriormente indicada exactamente igual, variando:
- la cantidad pesada de dicromato. Ahora de 0,8166 g
- la lectura de color, absorbancia, a 420 nm.
57
Determinacin de la absortividad12
La absortividad se determina, para varias muestras de un patrn primario de la DQO, en
concentraciones crecientes, que permitan abarcar el rango de validez del procedimiento
analtico a efectuar.
La absorbancia de una muestra leda en un colormetro o espectrofotmetro, viene dada por la
expresin:
A = a bc [5.1]
Donde:
A: Absorbancia
a: Absortividad
b: Ancho de la cubeta
c: Concentracin
Si en el procedimiento analtico se usan siempre cubetas de medida del mismo paso de luz, y
la concentracin de DQO se representa por sus iniciales, la Ecuacin [5.1], toma la forma de:
DQO = f A [5.2]
En la que f, denominado simplemente por factor, es el inverso del producto de la absortividad

por el ancho de la cubeta.
1
f = [5.3]
a b
El clculo de la f se realiza mediante la medida de las absorbancias de disoluciones
preparadas del patrn, de concentraciones crecientes y procediendo por mnimos cuadrados,
correspondientes al tramo recto de la representacin Absorbancia-Concentracin.
En las lecturas de la absorbancia de patrones de concentracin conocida, siempre habr una
desviacin experimental:
d = f A DQO [5.4]
Elevando las desviaciones al cuadrado, para eliminar el signo y haciendo el sumatorio de ellas,
se obtiene:
d 2
= f 2 A 2 2 f A DQO + DQO 2 [5.6]
La f, que hace mnimo el sumatorio de las desviaciones cuadrticas, se obtiene de las

condiciones del mnimo:
d 2
= 2 f A 2 2 A DQO = 0 [5.7]
f
de donde se obtiene:
f =
A DQO [5.8]
A 2
La absortividad ser por tanto:
1 A
2
a= [5.9]
b A dqo
12
caracterstica de una sustancia que nos dice cunta luz absorbe a una longitud de onda determinada
58
NOTA: En las determinaciones de la DQO se suelen utilizar tubos del mismo paso, tanto para
patrones como para muestras.
Si se utilizasen tubos de diferente paso, la ecuacin a aplicar es:
bp
DQO = f A [5.10]
bm
Donde bp y bm, son el ancho de las cubetas utilizadas de los patrones y de las muestras
problemas.
59
Programa en MatLab
IDEAS PREVIAS
Operar con matrices, vectores y escalares

Funciones
Sum
Input
Esta funcin posibilita la entrada de constantes, sin modificar el cdigo del programa.
X=input('Dime Absorbancia de Problema ') : MatLab pregunta y el usuario introduce el
valor, que ser X.
Num2str
Length
Round
Graficar
Ejes a medida
Line
60
Ejercicio
Para la medida de la DQO se preparan patrones de hidrgeno-ftalato potsico, que se

someten a digestin con reactivo de dicromato durante 2 horas a 150C, generndose un color
verde azulado por la presencia de Cr3+. Se dejan enfriar los patrones y se leen sus
absorbancias:
DQO Absorbancia
0 0,000
100 0,071
200 0,143
300 0,214
400 0,286
500 0,357
Sabiendo que en el procedimiento analtico se utilizan cubetas de medida del mismo paso de
luz, tanto para patrones como para muestras13, graficar la recta de calibracin
Absorbancia/Concentracin DQO (con crculos los datos reales y con lnea continua los
tericos).
Cual ser la DQO de una muestra con absorbancia 0,5? (input) Dibuja la recta que indique
dicho punto (line).
Solucin
%------------------- AbsDQO.m ----------------------------------------

A=[0 0.071 0.143 0.214 0.286 0.357]';
DQO=[0:100:500]';
SAporDQO=sum(A.*DQO);
SAcua=sum(A.*A);
f=SAporDQO/SAcua;
X=input('Dime Absorbancia del Problema ');
DQOp=f*X;
plot(A,DQO,'ko',A,f*A,'k-'), axis([0 A(length(A)) 0 DQO(length(DQO))])
sf=['\fontsize{13}f = ',num2str(round(f*10)/10)];
sx=num2str(X);
sx=['\fontsize{13}A la absorbancia de ',sx];
sfx=[sx,'\fontsize{13} corresponde ',num2str(round(f*X))];
s=[sfx,'\fontsize{13} mg/L de DQO'];
text(0.06,465,sf)
text(0.02,440,s)
line([X X],[0 DQOp],'Color','k','LineStyle',':')
line([X 0],[DQOp DQOp],'Color','k','LineStyle',':')
xlabel('\fontsize{13}Absorbancia'), ylabel('\fontsize{13}DQO mg/L')
%---------------------------------------------------------------------
13
Si no fuse el caso, debera canviar la programacin, haciendo uso de la Ecuacin 5.9 en lugar de la
Ecuacin 5.2.
61
Grfica
500
450
f = 1400.2
A la absorbancia de 0.25 corresponde 350 mg/L de DQO
400
350
300
DQO mg/L
250
200
150
100
50
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
Absorbancia
Figura 5.1: Recta de calibracin Absorbancia/Concentracin.
Aclaraciones:
SAporDQO=sum(A.*DQO) : Primero multiplica el vector A por DQO elemento a elemento, y

luego hace el sumatorio.
plot(A,DQO,'ko',A,f*A,'k-') : Grafica A en funcin de la DQO con circulos (ko) y A en

funcin de f*A con una recta (k-).
line([X X],[0 DQOp],'Color','k','LineStyle',':') : Dibuja una linea desde X hasta X en el eje x, y desde 0
hasta DQOp en el eje y. Esta recta es una lnea punteada de color negra.
axis([0 A(length(A)) 0 DQO(length(DQO))]) : Axis define los lmites del grfico ([xmin, xmax, ymin e
ymax]).
62
Tema 6 :
MODELOS DINMICOS DE
TANQUES
63
Tema 6 : MODELOS DINMICOS DE TANQUES
Resumen
Se desarrollarn los modelos matemticos de diferentes sistemas de alimentacin en un
tanque.
Introduccin
Para la introduccin a la dinmica de reactores, estudiamos operaciones fsicas en tanques,
mezclado, equalizacin de caudales y concentracin.
En este tema, se desarrollan los modelos matemticos y su solucin numrica de las
siguientes operaciones bsicas:
1. Altura de lmina de vertido y volumen de un tanque de agua que se somete a

caudales de entrada variables tipo escaln.
2. Transitorio de mezcla en un tanque de mezcla perfecta.
3. Dinmica de tanques de ecualizacin de caudal y concentracin.
64
Prctica 6.1 : ALTURA DE LMINA DE VERTIDO Y

VOLUMEN DEL TANQUE

Qumico: 1,0h
TOTAL 4,0h
Resumen
Se desarrollarn los modelos matemticos del sistema de alimentacin variable tipo escaln en
un tanque y obtener as todos los parmetros y poder graficarlos.
Objetivos
de:
Fsico-Qumicos
- Esquematizar el sistema de alimentacin de un tanque.
- Adquirir conocimiento de las variables que intervienen en dicho sistema, y la
relacin que existe entre ellas.
MatLab/Informticos
odeset
ode23
lenght
zeros
- Crear bucles, controles de flujo
for
if
- Graficar (plot)
axis
- Aadir etiquetas.
65
Fundamentos qumicos
Cuando un tanque se alimenta con un caudal variable de un lquido hasta que se llena, y
comienza a salir el lquido por rebose, surge una lmina del lquido por encima del borde de
vertido. Esta lmina adquiere un espesor dependiente del caudal de entrada y las dimensiones
del tanque. A su vez, la altura de esa lmina provoca la cuanta del caudal de salida.
El caudal de entrada es una variable impuesta, entrada; en simulacin se le asignarn

diferentes valores para observar la respuesta del sistema.
El volumen indicar el estado del sistema.
La altura de la lmina y el caudal de salida son consecuencia del estado del sistema,
constituyen la salida del sistema.
Siguiendo la nomenclatura normalmente utilizada en el anlisis de sistemas, la representacin
de este sencillo proceso es la siguiente:
Figura 6.1.1: Esquema del sistema dinmico de un tanque sometido a un caudal de entrada
variable.
Donde:
u representa la entrada al sistema, caudal de entrada.
x representa el estado del sistema, volumen del lquido en el tanque.
y representa las salidas del sistema, altura de lmina en el tanque y caudal de salida.
El esquema del sistema considerado se representa de la siguiente manera:
Figura 6.1.2: Esquema de un tanque con rebosadero lateral.
Considerando que el lquido circulante es agua, el caudal que sale por un vertedero
rectangular, [m3/s], de longitud L [m], viene dado por la expresin de Bazin:
Q = L h1,5 2 g [6.1.1]
66
donde es un coeficiente que depende de la altura de lmina h [m], y de la profundidad del

tanque H [m]:
0,003 h2
= 0,405 + 1 + 0,55
(H + h )2
[6.1.2]
h
NOTA: El caudal de salida, calculado con este procedimiento, es utilizable para alturas de
lmina de 2,5 a 80 cm.
La variacin instantnea del volumen de agua en el tanque es la diferencia de los caudales de

entrada al de salida:
dV
= Qi Qe [6.1.3]
dt
La altura de la lmina es el volumen total de agua entre la seccin menos la altura del
rebosadero:
V
h= H [6.1.4]
A
67
Programa en MatLab
IDEAS PREVIAS
Conceptos generales



Archivos *.m

if (condicional)
La estructura condicional, if, es un bloque que ejecuta las rdenes cuando se cumple la orden lgica.
if orden lgica 1
rdenes 1
elseif orden lgica 2
rdenes 2
else
rdenes 3
end
Funciones
Odeset
Crea o altera la estructura de las opciones para la entrada de la ecuacin diferencial.
>>opciones = odeset('nombre1',valor1,'nombre2',valor2,...)
Las opciones que se podrian canviar con la funcin odeset son:
RelT: Marca la tolerancia del error relativo que se aplica a todos los componentes del
vector residual.
AbsTol: Marca la tolerancia absoluta del error que se aplican a los componentes del vector
residual.
Ode23
>>[T,Y] = ode23(nombre de la funcin, t0 , tf, y0, opciones)
donde:
Ode23 integra un conjunto de ecuaciones diferenciales descritas en la funcin desde un instante inicial
t0, a un instante final tf, comenzando con la condicin inicial y0 y teniendo en cuenta las
opciones anteriormente modificadas.
68
Length
Zeros
La funcin zeros crea una matriz que solo contiene valores igual a zero.
zeros(i,j) = Crea una matriz de dimensiones i x j.
Graficar
Explicado en el Prctica 1, 2 y 4.
Ejercicio
Graficar, en una sola figura(subplot), los caudales de entrada, salida, altura de lmina y
volumen de un tanque que responde a las siguientes caractersticas:
- rea de la base de 150 m2
- longitud de 10 m
- altura 4 m
- 590 L iniciales de agua
- caudal medio de entrada escalonada de 100 m3/h, con variaciones de 50%
El tiempo de anlisis para este ejemplo se disminuy considerablemente en relacin a un
ciclo diario, ya que esta dinmica es bastante rpida (suponer 20 horas).
Las opciones para integrar la ecuacin de estado estn marcados por 1*105 de tolerancia
del error relativo y absoluto.
El rango de las ordenadas sn:

- Caudales de entrada: [-5 154]
- Altura de lmina: [-1 14]
- Volumen del tanque: [590 603]
NOTA: La altura de lmina oscila entre 2,5 y 80 cm.

El programa de clculo se estructura en tres partes, un programa y dos funciones. En el
programa principal se definen los datos de partida, desde el que se llama a la primera
funcin, y sta a la segunda, para finalmente realizar los clculos requeridos y presentacin
grfica de resultados. La primera funcin calcula el estado del sistema, volumen de agua en
el tanque X, en funcin de la diferencia de caudales de entrada, impuesto, y del de salida.
La segunda funcin calcula las salidas, Y, altura de lmina y caudal de salida.
69
Solucin
%---------------------- Vertedero.m ---------------

%
clear
clf
global A H L Qe1 Qe2 Qe3 tf
%
% Nomenclatura
%--------------------------------------------------
%
% A Superficie de la base
% H Altura del vertedero
% L Longitud del vertedero
%--------------------------------------------------
% Qe1 Caudal del primer tramo
% Qe2 Caudal del segundo tramo
% Qe3 Caudal del tercer tramo
%--------------------------------------------------
% V Volumen inicial de agua en el tanque
% t0, tf Tiempos inicial y final de anlisis
%--------------------------------------------------
% X Volumen variable del agua
% U Caudal de entrada
%--------------------------------------------------
% Y(:,1) Caudal de entrada
% Y(:,2) Caudal de salida
% Y(:,3) Altura de lmina
%--------------------------------------------------
%
% Parmetros y estado inicial
%--------------------------------------------------
%
A=150; H=4; L=10; tf=20*60;
Qe1=50/3600; Qe2=100/3600; Qe3=150/3600; V=590;
%--------------------------------------------------
%
% Integracin de la ecuacin de estado
%--------------------------------------------------
%
Opciones=odeset('RelT',1e-7,'AbsTol',1e-7);
[t,X]=ode23('VerteX',[0,tf],V,Opciones);
%--------------------------------------------------
%
% Clculo de las salidas
%--------------------------------------------------
%
Y=verteY(t,X);
%--------------------------------------------------
%
% Representacin de resultados
%--------------------------------------------------
%
T=t/60;
tff=T(length(T));
%
%--------------------------------------------------
subplot(311)
plot(T,Y(:,1)*3600,'*-',T,Y(:,2)*3600,'o--')
axis([0 tff -5 154])
ylabel('Cuadales (m^3/h)')
text(.5,59,'Caudal de entrada')
text(10.2,51,'Caudal de salida')
%
%--------------------------------------------------
subplot(312)
plot(T,Y(:,3)*1000,'*-'), grid
70
axis([0 tff -1 14])

ylabel('Altura de lmina (cm)')
%
%--------------------------------------------------
subplot(313)
plot(T,X,'*-'), grid
axis([0 tff 590 603])
ylabel('Volumen del agua (m^3)')
xlabel('Tiempo (min)')
%--------------------------------------------------
%--------------------- VerteX.m --------------------------

% Estado del sistema
% Clculo del volumen del tanque
%---------------------------------------------------------
%
function dX=VerteX(t,X)
%
%
% Clculo del caudal de salida
%---------------------------------------------------------
%
Y=VerteY(t,X); U=Y(1,1);
%
% Variacin de volumen del tanque
%---------------------------------------------------------
%
dX=U-Y(1,2);
%
%---------------------------------------------------------
%-------------------------------- VerteY.m --------------------

% Salidas del sistema
% Clculo de caudal de salida y altura de lmina
%--------------------------------------------------------------
%
function Y=VerteY(t,X)
%
%
N=length(t);
Y=zeros(N,3);
for i=1:N
% Caudal de entrada
if t(i)<tf/3
Y(i,1)=Qe1;
elseif t(i)<2*tf/3
Y(i,1)=Qe2;
else
Y(i,1)=Qe3;
end
h=X(i)/A-H;
% Altura de lmina
if h<=0
Y(i,2)=0; Y(i,3)=0;
else
mu=0.405+0.003/h*(1+ 0.55*h^2/(h+H)^2);
% Caudal de salida
Y(i,2)=mu*L*h*(2*9.8*h)^(1/2);
Y(i,3)=h;
end
end
%---------------------------------------------------------
71
Grfica
Cuadales (m3/h) 150
100
50 Caudal de entrada Caudal de salida
0
0 2 4 6 8 10 12 14 16 18 20
Altura de lmina (cm)
10
0
0 2 4 6 8 10 12 14 16 18 20
Volumen del agua (m3)
600
595
590
0 2 4 6 8 10 12 14 16 18 20
Tiempo (min)
Figura 6.1.3: Caudales de entrada, salida, altura de lmina y volumen de un tanque con caudal de
entrada escalonada.
Observacin acerca de los resultados
En la Figura 6.1.3, se aprecia que en muy pocos minutos, la altura de lmina y el caudal de
salida alcanzan los valores de estado estacionario, por ello se suele tomar el caudal de salida,
igual al de entrada en condiciones de rebose.
Aclaraciones:
Zeros(N,3) : Matriz de zeros, (N x 3).
Subplot (311) : Divide la figura en 3 divisiones horizontales y 1 verticlal, y est definiendo la primera grfica.
plot(T,Y(:,1)*3600,'*-',T,Y(:,2)*3600,'o--') : Grafica T en funcin de todas las filas de la primera

columna de la matriz Y multiplicadas por 3600 (horas) con *
para resaltar los puntos y para completar la recta. Una
segunda representacin con la segunda columna de Y con o
para resaltar los puntos y - - para completar la recta con lnea a
trazos.
axis([0 tff -5 154))]) : Axis define los lmites del grfico ([xmin, xmax, ymin e ymax]).
72
Prctica 6.2 : TANQUE DE MEZCLA PERFECTA

Qumico: 1,0h
TOTAL 4,0h
Resumen
Desarrollo de los modelos matemticos del sistema de alimentacin en un tanque con el
objetivo de calcular todos los parmetros y posteriormente hacer su representacin.
Dicho tanque es alimentado por un caudal y composicin variable en una sustancia que
consigue mezcla instantnea con el contenido del tanque.
Objetivos
de:
Fsico-Qumicos
- Esquematizar el sistema de alimentacin de un tanque con caudal y sustrato
de entradas escalonadas.
MatLab/Informticos
odeset
ode23
lenght
zeros
for
if
- Graficar (plot)
axis
- Aadir etiquetas.
73
Fundamentos qumicos
Para el anlisis de la mezcla en un tanque, se considera la instalacin de la prctica anterior
con un diseo que logre conseguir mezcla homognea de un sustrato influente, Si, con el
contenido del tanque.
De modo que el supuesto de funcionamiento es el siguiente:
Una alimentacin de caudal y composicin variable en una sustancia se introduce en un tanque
con el diseo representado en la Figura 6.2.1, para conseguir mezcla instantnea con el
contenido del tanque, se facilita por una tubera perforada y agitacin si fuese necesario; caso
semejante al que ocurre en un reactor biolgico, que es agitado en toda su extensin mediante
difusores de aire, siendo innecesaria la tubera de distribucin perforada para conseguir la
mezcla homognea.
La composicin del sustrato de entrada, variable en el tiempo, Si, en el supuesto analizado,
pasa instantneamente de la composicin de entrada en la tubera a la composicin
homognea del tanque, S, cuyo flujo sale por rebose, siendo la composicin de salida igual a la
del contenido del tanque.
Figura 6.2.1: Caudales y concentraciones de entrada, salida, altura de lmina y volumen de un

tanque con caudal y sustrato de entradas escalonadas.
El clculo de la composicin variable del contenido del tanque se realiza considerando que, en
un instante, la variacin de la cantidad de sustrato del interior del tanque es la diferencia de la
cantidad que entra a la que sale en el instante considerado, segn esto:
d (V S ) = Qi S i dt Qe S dt [6.2.1]
Reordenando y despejando la variacin del sustrato con el tiempo se obtiene:
dS Qi Q S dV
= Si e S [6.2.2]
dt V V V dt
La Ecuacin [6.2.2] se podria simplificar a:
dS Q
= (S i S ) [6.2.3]
dt V
- Si el objetivo de estudio son las concentraciones, y por tanto la dinmica del caudal de
salida puede despreciarse y tomarse como de respuesta instantnea.
- Si incrementramos el caudal de 50 a 150 m3/h, en un tanque de 600 m3,el aumento de

volumen del tanque seria de entre 2 y 3 m3, por lo que puede muy bien despreciarse a
efectos de anlisis de las concentraciones de mezclado.
74
Programa en MatLab
IDEAS PREVIAS
Conceptos generales



Archivos *.m

if (condicional)
Explicado en el Prctica 6.1.
Funciones
Odeset
Ode23
Length
Zeros
Graficar
75
Ejercicio
Representar en una misma figura, los caudales de entrada, salida, altura de lmina y volumen
de un tanque con caudal y sustrato de entradas escalonadas. Sabiendo que:
- La base del tanque tiene un rea de 150 m2, una longitud de 10 m y una altura
de 4 m
- 550 L iniciales de agua en el vertedero (tanque).
- El caudal de entrada escalonada es de 50, 150 y 50 mg/L.
- La concentracin de sustrato de entrada es de 100, 300 y 100 mg/L, primero,
segundo y tercer tramo.
- Se supone una mezcla instantnea entre el sustrato de entrada y el contenido
del tanque. Considerando que, en un instante, la variacin de la cantidad de
sustrato del interior del tanque es la diferencia de la cantidad que entra a la
que sale en el instante considerado.
- El tiempo de anlisis es de 48 horas.
Las opciones para integrar la ecuacin de estado estn marcados por 1*105 de tolerancia
del error relativo y absoluto.
El rango de las ordenadas sn:

- Caudales de entrada: [-5 154]
- Concentracin de sustrato: [95 305]
NOTA: Recuerda las ecuaciones de la prctica 6.1-
Solucin
%------------------------------- VertederoM.m ------------

%
clear, clf
global A H L Qe1 Qe2 Qe3 tf Se1 Se2 Se3
%
% Nomenclatura
%---------------------------------------------------------
%
% A Superficie de la base
% H Altura del vertedero
% L Longitud del vertedero
% Qe1 Caudal de entrada del primer tramo
% Qe2 Caudal de entrada del segundo tramo
% Qe3 Caudal de entrada del tercer tramo
% Se1 Conce. sustrato de entrada, primer tramo
% Se2 Conce. sustrato de entrada, segundo tramo
% Se3 Conce. sustrato de entrada, tercer tramo
% V Volumen inicial de agua en el tanque
% t0, tf Tiempos inicial y final de anlisis
%
% X Volumen variable del agua
% U(:,1) Caudal de entrada
% U(:,2) Concentracin de sustrato de entrada
% Y(:,1) Caudal de entrada
% Y(:,2) Caudal de salida
% Y(:,3) Altura de lmina
% Y(:,4) Concentracin de sustrato de entrada
%
76
% Parmetros
%---------------------------------------------------------
%
A=150; H=4; L=10; tf=48*60*60;
Qe1=50/3600; Qe2=150/3600; Qe3=50/3600;
Se1=100; Se2=300; Se3=100; %---- mg/l
%
% Estado inicial
%---------------------------------------------------------
%
V=550; t0=0; X0=V; S0=100;
%
% Integracin de la ecuacin de estado
%---------------------------------------------------------
%
Opciones=odeset('RelTol',1e-5,'AbsTol',1e-5);
[t,X]=ode23s('VerteXM',[t0,tf],[X0,S0],Opciones);
%
%
% Clculo de las salidas
%---------------------------------------------------------
%
Y=verteYM(t,X); Y3=Y(:,3)*1000;
%
% Representacin de resultados
%---------------------------------------------------------
%
T=t/(60*60); tff=T(length(T));
subplot(211)
plot(T,Y(:,1)*3600,'*-',T,Y(:,2)*3600,'o--', T,Y3,'d-')
axis([t0 tff -5 154]),xlabel('Tiempo (horas)')
ylabel('Caudales (m^3/h)')
text(3,58,'Caudal de entrada')
text(1.3,25,'Caudal de salida')
text(17,22,'Altura de lmina (mm)')
subplot(212); plot(T,Y(:,4),'*-',T,X(:,2),'*--')
axis([t0 tff 95 305]), ylabel('Conc. sustrato (mg/L)')
text(33,288,'Conc. sustrato de salida')
text(32.3,114,'Conc. sustrato de entrada')
%
%---------------------------------------------------------
%--------------------- VerteXM.m --------------------------

%
% Clculo del volumen del tanque
%---------------------------------------------------------
%
function dX=VerteXM(t,X)
%
%
% Clculo del caudal de salida
%---------------------------------------------------------
%
Y=VerteYM(t,X); U(1)=Y(1,1); U(2)=Y(1,4);
%
77
% Variacin de volumen del tanque

%---------------------------------------------------------
%
dX(1)=U(1,1)-Y(1,2);
dX(2)=U(1,1)/X(1)*U(2)-Y(1,2)/X(1)*X(2)-X(2)/X(1)*dX(1);
%
%---------------------------------------------------------
%------------------------------- VerteYM.m --------------------

%
% Clculo de Caudal de salida y altura de lmina
%--------------------------------------------------------------
%
function Y=VerteY(t,X)
%
%
N=length(t);
Y=zeros(N,3);
for i=1:N
% Caudal de Entrada
if t(i)<tf/3
Y(i,1)=Qe1; Y(i,4)=Se1;
elseif t(i)<2*tf/3
Y(i,1)=Qe2; Y(i,4)=Se2;
else
Y(i,1)=Qe3; Y(i,4)=Se3;
end
h=X(i)/A-H;
if h<=0
Y(i,2)=0; Y(i,3)=0;
else
mu=0.405+0.003/h*(1+ 0.55*h^2/(h+H)^2);
% Caudal de salida
Y(i,2)=mu*L*h*(2*9.8*h)^(1/2);
Y(i,3)=h;
end
end
%-------------------------------------------------------------
78
Grfica
150
Caudales (m3/h)
100
Caudal de entrada
50
Caudal de salida Altura de lmina (mm)
0
0 5 10 15 20 25 30 35 40 45
Tiempo (horas)
300
Conc. sustrato de salida
Conc. sustrato (mg/L)
250
200
150
Conc. sustrato de entrada
100
0 5 10 15 20 25 30 35 40 45
Figura 6.2.2: Evolucin de caudales y concentraciones de entrada, salida y altura de lmina de un

tanque con caudal y sustrato de entradas escalonadas.
Aclaraciones:
Zeros(N,3) : Matriz de zeros, (N x 3).
Subplot (211) : Divide la figura en 2 divisiones horizontales y 1 verticlal, y definie la primera grfica.
axis([0 tff -5 154))]) : Axis define los lmites del grfico ([xmin, xmax, ymin e ymax]).
Comentarios
La dinmica del mezclado continuo en un tanque, contrasta en lentitud con la de altura de
lmina o del caudal de salida, analizada en el ejemplo anterior, puesto en evidencia al
observar el tiempo en alcanzar el estado estacionario de estas magnitudes consideradas.
En la grfica obtenida del programa, da la impresin de que las dinmicas de caudal de

salida y altura de lmina son instantneas.
79
Prctica 6.3 : DISEO DE TANQUES REGULADORES

DE CAUDAL

Qumico: 1,0h
TOTAL 4,0h
Resumen
Se calcular el valor del volumen de un tanque de ecualizacin14.
Objetivos
de:
Fsico-Qumicos
- Comprender la importancia de un ecualizador.
- Esquematizar el sistema de un tanque con ecualizacin.
MatLab/Informticos
lenght
spline
max
round
fzero
for
if
- Graficar (plot)
axis
- Aadir etiquetas.
14
Encargado de ajustar la frecuencia de una seal.
80
Fundamentos qumicos
Los vertidos urbanos e industriales generan corrientes de aguas residuales de fuerte
variabilidad en caudal y concentracin, que originan serios problemas en la operacin y control
de las depuradoras.
Las ventajas de regular el caudal, mediante la instalacin de tanques de ecualizacin, conlleva

la suavizacin de la concentracin, presentando las siguientes ventajas con respecto a las
etapas de depuracin posteriores:
- Se evita que los procesos de depuracin sean perturbados por aumentos o
disminuciones de caudal inesperados, pudiendo causar ya desabastecimiento
o lavado del reactor biolgico.
- Se logra que en el tanque de retencin las descargas puntuales de
contaminantes se diluyan en el contenido global de este tanque.
- Se consigue controlar mejor la depuradora, ya que la alimentacin es
constante.
Un tanque de ecualizacin previo a una depuradora, permite alimentar la depuradora con un

caudal constante diario, como valor de la media mvil de las ltimas 24 horas, reteniendo el
agua durante caudales que superan su valor medio, para devolverlo en horas valle15.
Se plantea disear un tanque de ecualizacin de caudal que acumule la alimentacin en

exceso del valor del caudal medio, para devolverlo en las horas de alimentacin bajas, a la vez
que determinar la fraccin de tiempo de aportaciones mximas, para calcular el valor medio del
caudal en el tiempo de mximos, y lo mismo durante las aportaciones inferiores al valor medio
del caudal de entrada.
Existen dos metodologas para el diseo, segn los datos de partida:

- Volumen.
- Caudal.
15
Horas donde la alimentacin es menor.
81
Mediante el grfico Volumen-tiempo

El clculo del volumen del ecualizador se determina grficamente representando el volumen
aportado en un ciclo completo.
Ejemplo con datos imaginarios:
Figura 6.3.1: Clculo grfico del volumen del ecualizador.
El caudal medio es el cociente del volumen total vertido en el ciclo, dividido por el tiempo de
duracin del ciclo.
Qmedio =
Vvertido [6.3.1]
T
siendo T el tiempo de un ciclo productivo.
El volumen del tanque es la distancia mxima de las tangentes a la curva de volumen frente al
tiempo, paralelas a la bisectriz de la representacin.
82
Mediante la variacin del caudal frente al tiempo

Cuando se tiene la variacin del caudal con el tiempo en un ciclo productivo, el clculo del
caudal de salida se obtiene por el siguiente procedimiento:
Calcular el valor medio del caudal, Qmedio.
Trazar la lnea del caudal medio, en la grfica del caudal, obtenindose las
regiones de mximos y mnimos..
Obtener los valores medios de las zonas de mximos, Qmax, y de la de

mnimos, Qmin.
Calcular el tiempo, t, durante las aportaciones mximas.
T Qmedio = t Qmax + (T t ) Qmin [5.4.1]
Obtener el volumen del tanque, definido por el exceso del volumen aportado
en el tiempo de mximos:
V = t (Qmax Qmedio ) [5.4.2]
83
Programa en MatLab
IDEAS PREVIAS
Conceptos generales



Archivos *.m

Explicado en el Prctica 1 y 6.1.
Funciones
Odeset
Round
Length
Num2str
Spline
Interpolacin cbica de los datos de los vectores introducidos
Fzero
Encuentra el valor zero de la funcin, desde el valor x0 y teniendo en cuenta las opciones
definidas.
>> x = fzero(funcin,x0,opciones)
Graficar
84
Ejercicio
Sabiendo que:
- Se alimentar un tanque con ecualizador cada 2 horas, desde un tiempo

inicial 0 hasta 22 (lo que hace un total de 12 tiempos).
- Los valores de los caudales de entrada respecto al tiempo sn:
Q=[211 168 155 162 220 330 460 533 475 370 290 250]
- La funcin cuadal y tiempo se cruza 2 veces con el caudal medio. Una primera
vez despues del tiempo 9.6 y una segunda despues del tiempo 19,4 (fzero).
Calcular:
a) el volumen del ecualizador
b) representar el caudal respecto al tiempo, identificando(subplot):
a. el caudal medio mnimo
b. el caudal medio mximo
c. el caudal medio
d. el tiempo exacto donde se produce el cambio del caudal medio
mnimo al caudal medio mximo
NOTA: El diseo de este ecualizador se estudia bajo la condicin de no prever de:

- Un volumen de reserva mnimo.
- Una reserva de volumen para perturbaciones extraordinarias.
Solucin
%---------------RegCaudal.m--------------------------------
clear, clf
global TN QN QM
T=[0:2:22]';
Q=[211 168 155 162 220 330 460 533 475 370 290 250]';
TN=[T(1):.01:T(length(T))]';
QN=spline(T,Q,TN);
plot(TN,QN,'k-')
axis([T(1) T(length(T)) 0 max(1.05*Q)])
xlabel('Tiempo (horas)')
ylabel('Caudal (m^3/h)')
hold on
plot(T,Q,'k*')
QM=sum(QN)/length(QN);
line([T(1) T(length(T))],[QM QM],'Color','k')
t1=fzero('RegCaudF',9.6);
t2=fzero('RegCaudF',19.4);
n=length(TN);
n1=round(n*t1/TN(length(TN))); % N datos hasta t1
n2=round(n*t2/TN(length(TN))); % N datos hasta t2
n3=n-n2; % N datos desde t2 hasta n
Qmin=(sum(QN(1:n1))+sum(QN(n2:n)))/(n1+(n-n2));
Qmax=sum(QN(n1:n2))/(n2-n1);
line([TN(1) TN(n1)],[Qmin Qmin],'Color','k')
line([TN(n1) TN(n2)],[Qmax Qmax],'Color','k')
line([TN(n2) TN(n)],[Qmin Qmin],'Color','k')
line([TN(n1) TN(n1)],[Qmin Qmax],'Color','k')
line([TN(n2) TN(n2)],[Qmax Qmin],'Color','k')
V=(Qmax-QM)*(t2-t1);
Texto1=['Qmin = ',num2str(round(Qmin*10)/10),' m^3/h'];
text(2,218,Texto1)
Texto2=['Qmax = ',num2str(round(10*Qmax)/10),' m^3/h'];
85
text(12,447,Texto2)
Texto3=['Qmed = ',num2str(round(10*QM)/10),' m^3/h'];
text(2,318,Texto3)
Texto4=['Volumen del ecualizador = ',num2str(round(10*V)/10),' m^3'];
text(8,80,Texto4)
Texto5=['t_1 = ',num2str(round(100*t1)/100),' horas'];
Texto6=['t_2 = ',num2str(round(100*t2)/100),' horas'];
%--------------------------------------------------------
%---------RegCaudF.m-------
function D=RegCaudF(TT)
global TN QN QM
D=spline(TN,QN,TT)-QM;
%--------------------------
Grfica
500
Qmax = 436.5 m3/h
400
Caudal (m3/h)
Qmed = 308.1 m3/h

300
Qmin = 205.6 m3/h

200 t1 = 9.65 horas t2 = 19.45 horas
100
Volumen del ecualizador = 1258.2 m3
0
0 2 4 6 8 10 12 14 16 18 20 22
Tiempo (horas)
Figura 6.3.2: Clculo numrico del volumen del ecualizador
Aclaraciones:
Qn(1:n1) : Valores de Qn. Desde Qn(1) hasta Qn(n1)
NOTA: El volumen del ecualizador calculado, por el mtodo grfico o numrico, ha de

sobredimensionarse en un factor de prudencia que asegure que el tanque ecualizador no se
quede vaco por un desabastecimiento, ni que puedan ocurrir rebosamientos por aportaciones
extraordinarias, lo que lleva a dimensionar el tanque, pensando que haya un volumen de
reserva mnimo, y que en el punto de volumen mximo previsible an reserve suficiente
capacidad para perturbaciones extraordinarias.
86
5.4 Conversin a crditos ECTS

De acuerdo con las directrices marcadas por el Tratado de Bolonia, se establece a
continuacin el nmero de crditos ECTS de cada prctica y totales, estimando que el
alumno trabaja 40h/semana, durante 40 semanas,
1 curso acadmico = 1600 h de trabajo real
y teniendo en cuenta que,

1 curso acadmico = 60 crditos ECTS
entonces,
1 crdito ECTS = 26,7 h de trabajo real
Trabajo real del alumno (h)

Clase Crditos ECTS
Presencial No presencial
Introduccin 2 3,0 0,2
Prctica 1 4 3,0 0,3
Prctica 2 4 4,0 0,3
Prctica 3 2 2,5 0,2
Prctica 4 4 4,0 0,3
Prctica 5 4 4,0 0,3
Prctica 6.1 4 5,0 0,3
Prctica 6.2 4 5,0 0,3
Prctica 6.3 4 5,0 0,3
Tutora 2 0,5 0,1
Proyecto 4 45,0 1,8
TOTAL 38 81,0 4,5
87
5.5 Sistema de evaluacin

Para la evaluacin de esta asignatura, se contabilizan los siguientes conceptos:
1. Los objetivos MatLab/Informticos de cada prctica, asumidos correctamente.

2. La realizacin de un proyecto con 10 objetivos MatLab/Informticos, estudiados
en el transcurso de la asignatura y 2 ms que no hayan aparecido (se acepta
otro porcentaje de objetivos nuevos respecto a los estudiados, siempre que se
cumpla un mnimo de 2 nuevos).
3. La defensa del proyecto, donde se evaluarn 3 objetivos:
o Presentacin.
o Contenido.
o Explicacin.
4. La asistencia de las 15 clases.
Los conceptos 1, 2, 3 y 4 se evalan con los siguientes porcentajes:
Porcentaje por
Concepto Objetivos asumidos Porcentaje total
objetivo
Prctica 1 10 0,5% 5,0%
Prctica 2 14 0,4% 5,0%
Prctica 3 6 0,6% 3,5%
Prctica 4 14 0,4% 6,0%
Prctica 5 8 0,5% 4,0%
Prctica 6.1 14 0,5% 7,0%
Prctica 6.2 14 0,5% 7,0%
Prctica 6.3 14 0,5% 7,5%
Proyecto 12 1,7% 20,0%
Defensa 3 3,3% 10,0%
TOTAL 100,0%
El cuarto concepto se evala aplicando el siguiente porcentaje:
Concepto Porcentaje por clase Porcentaje total

Asistencia 1,7% 25,0%
La suma de los porcentajes, como es lgico, dar el 100% de la nota total de la asignatura.
Dicha nota se expresar con un nico decimal, tal y como dictamina el Tratado de Bolonia.
88
5.6 Logstica
Es importante considerar una serie de aspectos relacionados con el apoyo logstico y la
organizacin que debe tener una asignatura de esta naturaleza. Por ello, se han incluido
algunas recomendaciones que, aunque obvias, son importantes de tener en cuenta, pues de
ello depende en gran medida el xito que se pueda lograr en la implantacin de la asignatura.
5.6.1. Coordinacin
La entidad encargada de la asignatura debe nombrar un coordinador, quien ser el
responsable de la ejecucin de las actividades tericas y prcticas del curso, as como de
todos los aspectos relacionados con los materiales y equipos requeridos durante la actividad.
5.6.2. Sitio de trabajo

El ambiente y comodidad que ofrezca el aula donde se implante la asignatura es de suma
importancia pues, adems de facilitar la labor del profesor, permite que los asistentes
aumenten su rendimiento.
En la seleccin de las aulas para la imparticin de dicha asignatura, es necesario tener en

cuenta la comodidad de los participantes, as como propiciar un ambiente de trabajo favorable.
Las instalaciones deben ser confortables, funcionales y con espacios suficientes.
El aula ha de disponer durante todo el tiempo de imparticin de la asignatura de los siguientes

equipos:
o Ordenadores (profesor y alumnos) con capacidad suficiente como para ejecutar el
programa MatLab sin problemas.
o Proyector
o Pantalla
o Pizarra
Todo ello supervisado por el plan de mejor de las condiciones ergonmicas en trabajos con
pantallas de visualizacin de datos (6).
89
5.7 Bibliografa recomendada

5.7.1. Gestin de aguas
Manuel Gil Rodrguez. Procesos de descontaminacin de aguas. Clculos Avanzados

Informatizados. Editorial Thomson.
http://www.ina.gov.ar/cartillas_edu/cartilla_6.htm
http://www1.ceit.es/Asignaturas/Ecologia/Hipertexto/11CAgu/180Depur.htm
http://es.wikipedia.org/wiki/Tanque_de_agua
5.7.2. MatLab
http://www.tayuda.com/ayudainf/aprendainf/varios.htm
http://www.sisoft.ucm.es/Manuales/Matlab.pdf
http://nereida.deioc.ull.es/~pcgull/ihiu01/cdrom/matlab/contenido/node1.html
90
5.8 Presupuesto
Los costes para realizar este proyecto se dividen en:

a) Costes de material.
b) Costes de energia.
c) Costes destinados a recursos humanos.
a) Costes de material.
En la siguiente tabla se muestra el material necesario para la realizacin del proyecto, as como
su coste econmico y su amortizacin.
Nmero Coste
Vida til Coste Amortizacin
Unidades de Unitario
[aos] total [] []
unidades [/unidad]
Licencia programa MatLab 2 0* --------- 0,00
Libro Procesos de
2 17,9 --------- 35,80
descontaminacin de aguas
Ordenador 2 900 6 1800,00 300,00
TOTAL 1835,80 300,00
NOTA: La amortizacin se obtiene dividiendo el coste total entre la vida til de la unidad.
*Subvencionado por la universidad.
b) Costes de energa.
Sabiendo que,
- La potencia media consumida de un ordenador es de 0,154 kW.
- Las horas se encontuentran, aproximadamente, mediante las horas de trabajo
de la ingeniera restadas a las horas de reunin con los ponentes. Esta
igualacin es correcta, ya que se desprecia las horas de energia necesarios
para los ponentes suponiendo que son equivalentes a las horas de trabajo de
la ingeniera que no hace uso del ordenador.
- La empresa subministradora es Endesa y por tanto:
Se le aplica un 16 % de IVA.
El precio esta establecido en 0,08742 ./kWh
Obtenemos:
Energia Kw Horas Precio Unitario [/kWh] IVA [] Coste []

Electricidad 0,154 200 0,0874 0,43 3,12
c) Costes de recursos humanos.

El coste del personal se divide en dos categoras, Ingeniera Tcnica Industrial, especializada
en qumica y ponente.
Para obtener el coste del personal se ha de tener en cuenta esta ecuacin:
SBA + SS
Coste = RH
horas
ao
91
donde,
RH = Horas * Persona
SBA = Salario Bruto Anual
SBAIngeniera tcnica = 12500 / trabajador
SBAPonente = 24000 / trabajador
SS = Cuota de la empresa de la Seguridad Social (0,32*SBA)
Horas/ao 1600 (1 curso acadmico)
Aplicando dicha ecuacin, obtenemos:
Horas de
RH Nde SS SBA
Categora trabajo Coste []
[horas*persona] Personas [/unidad] []
anual
Ingeniera tcnica industrial,
200 1 4000 12500 1600 2062,50
especializada en qumica
Ponente 15 2 7680 24000 1600 594,00
TOTAL 2656,50
Por lo tanto:
Tipo de Coste Coste Total []

Coste de material 1.535,80
Coste de energa 3,12
Coste de recursos humanos 2.656,50
TOTAL 4.195,42
El diseo de esta asignatura ha producido un gasto total de 4.195,42 .
92
5.9 Conclusiones
Una vez finalizado el diseo de una asignatura de la titulacin de Ingeniera Tcnica Industrial,
especialidad en Qumica Industrial, conveniente para la sociedad actual, podemos concluir el
proyecto diciendo que:
- La asignatura responde como herramienta que ejercita:

El desarrollo a los conocimientos qumicos y fisico-qumicos que
aparecen en los procesos de determinacin de la contaminacin del
agua como:
o La solubilidad.
o La autodepuracin, sus etapas (Desoxigenacin y
reoxigenacin) y constantes correspondientes (k1 y k2)
o Formulacin emprica.
o DBO, conceptos relacionados (DBOU, DBO5, DOU,
DON) y su cintica.
o DQO y su medida experimental.
El progreso en los fundamentos de la Ingenieria Qumica necesarios

para el diseo y operacin de los procesos de tratamientos de agua
como:
o Los sistemas de alimentacin de un tanque.
La estimulacin de la simulacin por ordenador, con el MatLab,

haciendo uso de multiples rdenes como:
o Definir, operar y direccionar matrices, vectores y
escalares.
o Crear y utilizar funciones. (archivos *.m)
o Programar controles de flujo o bucles.
o Graficar.
- La asignatura pude ser aplicada dentro del Espacio Europeo de Educacin

Superior, ya que cumple con las condiciones ms conocidas y aceptadas,
hasta el momento, del tratado de Bolonia, como:
El uso de crditos ECTS.
El sistema de evaluacin.
- Dicho diseo ha provocado un coste total de 4.195,42 .
93
6. REFERENCIAS
(1) http://www.crue.org/apadsisuniv.htm
(2) http://es.wikibooks.org/wiki/Ingenier%C3%ADa_de_aguas_residuales/La_calidad_del_a
gua._Marco_jur%C3%ADdico
(3) http://es.wikipedia.org/wiki/Planta_de_depuraci%C3%B3n
(4) http://fc.uni.edu.pe/labosin/
(5) http://www.tayuda.com/ayudainf/aprendainf/varios.htm
(6) http://sites.upc.edu/~w-preven/AvalErgonomica/pla/index.html
(7) Manuel Gil Rodrguez. Procesos de descontaminacin de aguas. Clculos Avanzados

Informatizados. Editorial Thomson.
(8) Octavo informe de UK Royal Commission. Standards and Tests for Sewage and
Sewage Effluents Discharging into Rivers and Streams.
(9) Agencia de Proteccin Ambiental Americana, Environmental Protection Agency, EPA o

USEPA (http://www.epa.gov/).
(10) International Association on Water (http://www.iwa.org/)
(11) http://es.wikipedia.org/wiki/Tanque_de_agua
(12) http://www.sisoft.ucm.es/Manuales/Matlab.pdf
(13) http://www.au.uva.es/web/publicaciones/convergenciaeuropea
(14)http://www.criecv.org/es/proyectos/pag_agua/agua.html
(15) David Robson. Simulation saves time.

Scientific Computing World. December 2006/January 2007. ISSU 91. (Page 18)
(16) http://www.aguamarket.com/diccionario/index.asp
(17)http://www.quimicaviva.qb.fcen.uba.ar/contratapa/aprendiendo/capitulo2.htm
(18)http://www.unavarra.es/genmic/curso%20microbiologia%20general/40-
tratamiento%20aguas%20residuales.htm
(19) http://www.proyectosfindecarrera.com/
(20) Factura de Endesa.
94
7. AGRADECIMIENTOS
El proyecto se ha finalizado con xito gracias a la colaboracin de los ponentes Moiss Graells
y Montserrat Prez. A los que agradezco su motivacin e implicacin en el proyecto.
Miriam Mrquez Boada
95

Procesos de Tratamientos de Aguas Residuales y PDF

Cargado por

Copyright:

Formatos disponibles

Procesos de Tratamientos de Aguas Residuales y PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Procesos de Tratamientos de Aguas Residuales y PDF

Cargado por

Copyright:

Formatos disponibles

Proyecto Final de Carrera

Autora: Miriam Mrquez

5. DISEO DE LA ASIGNATURA .................................................................................................14

Tema 1 : SOLUBILIDAD ...........................................................................................................17

Tema 3 : FRMULA EMPRICA...............................................................................................36

Tema 4 : DEMANDA BIOQUMICA DE OXGENO ...............................................................42

Tema 5 : DEMANDA QUMICA DE OXIGENO ......................................................................51

Tema 6 : MODELOS DINMICOS DE TANQUES..................................................................63

5.4 Conversin a crditos ECTS ........................................................................................................87

1. El contexto en que se encuentra, y por tanto las pautas que ha de seguir la

2. El tema qumico y fisico-qumico a estudiar:

3. El mtodo utilizado para dicho estudio:

4.1 Espacio Europeo de Educacin Superior (E.E.E.S.)

El ECTS se ha desarrollado para facilitar el reconocimiento acadmico de los estudios

4.2 Gestin del agua

Los parmetros caractersticos son:

4.3 Simulacin de Procesos

MatLab dispone de un cdigo bsico y de varias libreras especializadas (toolboxes).

Figura 1. Icono de MatLab 7.0.4

Figura 2. Ventana inicial

Al arrancar el programa se abre una ventana similar a la mostrada en la Figura 2. sta es la

4.3.1. Descripcin de la ventana inicial

1. Perfil del alumno

5.1 Perfil del alumno

- Tema 1 y 2: Estudian los conceptos de la solubilidad y la autodepuracin.

- Tema 3: Detalla la determinacin de la frmula emprica de los componentes

- Tema 4 y 5: Tratan los parmetros ms importantes en la determinacin de la

- Tema 6: Introduce las simulaciones de tanques.

- Ttulo del tema a tratar.

*Este apartado no aparece en alguna prctica.

Tema 2 : AUTODEPURACIN ........................................................................................................26

Tema 3 : FRMULA EMPRICA ....................................................................................................36

Tema 4 : DEMANDA BIOQUMICA DE OXGENO ....................................................................42

Tema 5 : DEMANDA QUMICA DE OXIGENO ...........................................................................51

Tema 6 : MODELOS DINMICOS DE TANQUES.......................................................................63

Horas de trabajo real presencial

Solubilidad del oxgeno

Presin parcial del oxgeno igual a 1 atmsfera

1,5757 10 5 6,642310 7 1,24381010 8,6219 1011

Presin parcial del oxgeno diferente a 1 atmsfera

La presin de vapor del agua, en mm de Hg, en funcin de la temperatura centgrada, viene

22, 4077t +1167, 7614

Pws = 0,040032 e t + 245, 9443

Definicin de matrices, vectores y escalares.

Otras formas de introducir vectores

2. Mediante la funcin linspace:

NOTA: Para espaciado logartmico se utilizari la funcin logspace. (logspace(x1,x2,n)).

Direccionamiento de matrices y vectores

Acceso elemento a elemento

2. La indicacin de las posiciones de los elementos que queremos obtener:

Operaciones de matrices, vectores y escalares

Orientacin del vector o matriz

Controles de flujo o bucles

Estilo de lneas, marcadores y colores

SIMBOLO COLOR SIMBOLO ESTILO DE LINEA

Adicin de rejillas y etiquetas

Representa (plot) la cloronidad en funcin de la concentracin de oxgeno y la temperatura.

%--------------------- SoluOxiTemp.m ----------------------

Conc. oxgeno disuelto, mg/L 14