Presentación Clase 4

MAESTRÍA EN INGENIERÍA CIVIL
CONTAMINACIÓN EN RÍOS
CINETICA DE LAS REACCIONES
En la figura se esquematizan los procesos que puede
sufrir un contaminante dentro de un sistema de
flujo, que pueden ser de transporte con el flujo de
entrada como las descargas o salidas por
volatilización o sedimentación y la reacción con
otras substancias dentro del sistema.
Un enfoque simple para determinar la evolución del
contaminante es el de suponer una serie de botellas
llenas de agua donde el contenido se mantiene bien
mezclado. Tales recipientes se denominan
comúnmente como reactores en lotes.
Al medir la concentración en cada botella con el
tiempo se determinará la evolución de la
concentración en el tiempo.
FUNDAMENTOS DE LA REACCIÓN
Primero se debe realizar unas definiciones:
Tipo de reacciones:
Reacciones heterogéneas: Involucran más de una

fases (o estados) donde la reacción ocurre
habitualmente en la superficie de interface.
Reacciones Homogéneas: Involucra una sola fase (o

estado). Debido a que es la reacción fundamental en
este curso nos centraremos principalmente en ellas.
Tipo de reacciones:
Reacción reversible: Puede darse en cualquier dirección
y depende de las concentraciones relativas de los
reactivos y de los productos
𝑎𝐴 + 𝑏𝐵 ⇆ 𝑐𝐶 + 𝑑𝐷
Donde las minúsculas corresponden a los coeficientes
estequiometricos y las mayúsculas designan los
compuestos que reaccionan. Estas reacciones tienden a
aproximarse a un estado de equilibrio. Son la base para
la llamada química del equilibrio (pH).
Tipo de reacciones:
Reacciones irreversibles: Estos procesos se dan en una
sola dirección y continúan hasta que se agoten los
reactivos. En estos casos se trata de determinar la tasa
de desaparición de una o más sustancias que
intervienen en la reacción, por ejemplo:
𝑎𝐴 + 𝑏𝐵 → 𝑐𝐶 + 𝑑𝐷
Por ejemplo, se puede estar interesado en la tasa a la
que la sustancia desaparece.
Tipo de reacciones:
El ejemplo común de una reacción irreversible es la
descomposición de la materia orgánica, que en general
se representa como:
𝐵𝑎𝑐𝑡𝑒𝑟𝑖𝑎𝑠
𝐶6 𝐻12 𝑂6 + 6𝑂2 6𝐶𝑂2 + 6𝐻2 𝑂
Donde 𝐶6 𝐻12 𝑂6 es glucosa, que se toma como una
representación simple de la materia orgánica. Cuando el
A.R se descarga en un cuerpo de aguas tiene lugar una
reacción de ese tipo.
Tipo de reacciones:
La materia orgánica es oxidada por bacterias para
producir dióxido de carbono y agua. Aunque la
fotosíntesis es una de las reacciones inversas donde el
crecimiento de las plantas produce materia orgánica y
produce oxigeno, consumiento CO2 y H2O, normalmente
no ocurre en la vecindad de la descomposición. Ambos
procesos son muy lentos y no tienden al equilibrio, por
lo tanto no se consideran como reacciones reversibles.
CINÉTICA DE LA REACCIÓN
La cinética o tasa de la reacción puede expresarse
cuantitativamente por la ley de acción de masas, la cual
estable que la velocidad de la reacción es proporcional
a la concentración de los reactivos y se puede
representar como: dCA
 k f (CA , CB ...) 14
dt
Esta reacción se denomina ley de tasas. En ella se
especifica que la velocidad de reacción depende del
producto de una constante K, dependiente de la
temperatura, y de una función de las concentraciones
de los reactivos .
La función de las concentraciones se determina
habitualmente de forma experimental. Una forma
común es: dCA
 k C  C  15
A B
dt
Las potencias a las que se elevan las concentraciones se

denomina órdenes de la reacción.
En la ecuación la reacción tiene un orden con respecto
al orden del reactante A, con respecto al reaccionante B
(α). El órden global de la reacción es:
El orden global o el de cada uno de los reactantes no
necesariamente son entero. Sin embargo, varias de las
reacciones más importantes en los modelos de calidad
presentan órdenes enteros.
En este parte del curso nos concentramos en un único
reactivo. Para este caso la ecuación de la tasa de
reacción será:
Donde C es la concentración del reactante y n es el

orden de la cinética.
Aunque existe un número infinito de órdenes para la
ecuación, es habitual caracterizarlas como de orden 0, 1
y 2 para las habitualmente empleadas en sistemas de
agua naturales.
Para la cinética de orden cero la ecuación a integrar es:
𝑑𝐶 𝑀
= −𝑘 3 𝑐𝑜𝑛 𝐶 = 𝐶0 𝑒𝑛 𝑡 = 0
𝑑𝑡 𝐿 𝑇
Separando variables e intregrando
0 𝑡
𝑑𝐶 = − 𝑘 𝑑𝑡 ∴ 𝐶 − 𝐶0 = −𝑘𝑡
𝐶0 0
Es decir
𝐶 = 𝐶0 − 𝑘𝑡
En conclusión si se presenta una tasa de decrecimiento constante a partir de una

concentración inicial y su representación será una línea recta, la cinética será de orden
cero. Indica un decaimiento de la concentración en el tiempo; por ejemplo
0.3 mg/L día, la concentración decae en 0.3 mg/L cada día.
• Para la cinética de orden uno:

𝑑𝐶
= −𝑘𝐶 donde 𝐶 = 𝐶0 para 𝑡 = 𝑡0
𝑑𝑡
Separando variable e integrando
𝐶 𝑡
𝑑𝐶 𝐶
=− 𝑘 𝑑𝑡 ⇒ ln = −𝑘𝑡 ∴ 𝐶 = 𝐶0 𝑒 −𝑘𝑡
𝐶0 𝐶 0 𝐶0
En esa ecuación se refleja que el agotamiento es exponencial, es decir que la

curva de concentración es asintótica al eje del tiempo y la concentración decae y
tiende rápidamente a cero en el tiempo.
En esta ecuación la tasa de decaimiento tiene unidades de tiempo a la menos

uno, que es más difícil de interpretar, pero se puede realizar dicha interpretación
a través de la aproximación en serie de MacLaurin a una función exponencial.
𝑥 𝑥2 𝑥3
𝑒 −𝑥 =1− + − + …
1! 2! 3!
Si truncamos la serie luego del término de primer orden
𝑒 −𝑥 = 1 − 𝑥
Con base en esto, se puede interpretar que para tasas menores a 0,5 d-1 esta será
la fracción de concentración que pierde en un día. Para tasa mayores a 0.5 la
interpretación se puede obtener realizando un cambio de unidades así por
ejemplo para una tasa K de 6 d-1
 d 
K  6 d 1  6    0,25 hr 1
 24 hr 
Que nos indica que se pierde el 25% cada hora.
Donde se puede observar que la concentración decae desde un valor inicial con una curva cóncava hacia
arriba asintótica al eje del tiempo.
Esta expresión se puede generalizar para valores de n enteros con n ≠ 1
dC
  KC n 25
dt
Separado variable e integrado para condiciones entre las iniciales C= Co para
t= O y C para el tiempo t.
𝐶 𝑡 𝐶
𝑑𝐶 −1
𝑛
=− −𝑘 𝑑𝑡 ⇒ 𝑛−1
= −𝑘 𝑡
𝐶0 𝐶 0 𝑛−1 𝐶 𝐶0
−1 −1
− 𝑛−1 𝑘𝑡 = - ⇒
𝐶 𝑛−1 𝐶0𝑛−1
1 −1
𝑛−1 = + 1−𝑛 𝑘𝑡 ⇒
𝐶 𝐶0𝑛−1
1
𝐶 𝑛−1 = −1 ⇒
+ 1−𝑛 𝑘 𝑡
𝐶𝑛−1
Para la cinética de orden n se tiene:
1
𝐶 𝑛−1 = 𝐶0𝑛−1 ⇒
−1+ 1−𝑛 𝑘𝐶0𝑛−1 𝑡
Finalmente:
1
𝐶 = 𝐶0 1
𝑛−1 𝑛−1
1 + 𝑛 − 1 𝐾𝐶0 𝑡
ANÁLISIS DE LOS DATOS PARA ESTABLECER LAS TASAS DE DECAIMIENTO
Existen una gran variedad de caminos para analizar los datos de los reactores
por lotes. Aquí se partirá de la ecuación general para ilustrar algunas
técnicas.
EL MÉTODO INTEGRAL
Este método consiste básicamente en suponer un valor para n e integral la

ecuación 𝑑𝐶 𝑑𝑡 = 𝑘𝐶 𝑛 para obtener C (t).
Luego se emplean métodos gráficos para establecer si el modelo se ajusta

adecuadamente a los datos.
Los métodos gráficos son la realización de versiones linealizadas de la

funciones subyacentes para los modelos de cinética de orden n enteros
vistos.
EL MÉTODO INTEGRAL
En la tabla se resumen las estrategias a utilizar para aplicar el método
integral.
Unidades Variable Variable Pendien

Orden Intercepto
de las tasas Dependiente Y Independiente X te
Cero (n = 0) M/L3 T C t Co ̶ K
Primer (n =1) T-1 ln C t LnCo ̶ K
Segundo (n = 2 L3/MT-1 1/C t 1/Co K
General(n ≠ 1) (𝐿3 𝑀 −1 )𝑛−1 Cn-1 t C0n-1 (n – 1)K

𝑇
EL MÉTODO DIFERENCIAL
Se basa en aplicar logaritmos para transformar la ecuación
En consecuencia si se grafica en un plano logarítmico se obtendrá una recta

donde el intercepto es log K y la pendiente será n condiciones que resultan
ventajas porque el método proporciona una estimulación del orden n. Se
desventaja es que se debe obtener una estimulación de la derivada, para
esto se puede usar diferentes técnicas una puede ser por diferenciación.
EL MÉTODO DIFERENCIAL
Diferenciación numérica
Usando diferencias finitas se puede obtener una aproximación de las
derivadas, por ejemplo; cuando se usan:
Diferencias centrales:
𝑑𝐶 ∆𝐶 𝐶𝑖+1 − 𝐶𝑖−1
≅ =
𝑑𝑡 ∆𝑡 𝑡𝑖+1 − 𝑡𝑖−1
O diferencias hacia adelante:
𝑑𝐶 ∆𝐶 𝐶𝑖+1 − 𝐶𝑖
≅ =
𝑑𝑡 ∆𝑡 𝑡𝑖+1 − 𝑡𝑖
O diferencias regresivas:
𝑑𝐶 ∆𝐶 𝐶𝑖 − 𝐶𝑖−1
≅ =
𝑑𝑡 ∆𝑡 𝑡𝑖 − 𝑡𝑖−1
EL MÉTODO DE LAS TASAS INICIALES
Existen reacciones en las que se presentan complicaciones a lo largo del tiempo.

Por ejemplo se producen reacciones inversas significativas o las reacciones son
muy lentas en el tiempo, requiriendo experimentos con tiempos muy prolongados,
que los hace costosos. Para este tipo de reacción se puede utilizar el método de
las tasas iniciales, para determinar la tasa y el orden del modelo. En este método
se experimenta con varias concentraciones iniciales Co. Para cada experimento
la velocidad inicial de la variación - dCo/dt se halla diferenciado los datos y
extrapolando para el tiempo cero. En cada caso se usa el gráfico:
 dCo 
log   vs log(Co)
 dt 
de la ecuación
En donde la pendiente dará el orden y el intercepto el log k 

n
EL MÉTODO DE LA VIDA MEDIA
La vida media de una reacción es el tiempo que tarda la concentración en caer

hasta a mitad de su valor inicial.
Donde t50 es la vida media. Nuevamente se utiliza, como modelo de tasa, la

ecuación
1
𝐶 = 𝐶0 1
𝑛−1
1 + 𝑛 − 1 𝐾𝐶0 𝑡 𝑛−1
En donde la pendiente dará el orden y el intercepto el valor de log k 

Si C  Co en t  O que puede ser integrado como:
Combinando las ecuaciones resulta que
Al tomar el logarítmo de la ecuación
2𝑛−1 − 1
log 𝑡50 = log + 1 − 𝑛 log 𝐶0
𝑘 𝑛−1
La grafica del logaritmo vida media contra el logaritmo de la concentración inicial

permitirá que con la pendiente se obtendrá n-1 y conocido n, con el intercepto
hallado este permite calcular la tasa k.
Es claro que la selección de la vida media es arbitrario puede haberse seleccionado

cualquier otro porcentaje de la concentración inicial y su respectivo tiempo de
ocurrencia. En esa condición la ecuación será:
n 1
100   
 100   1 1
t  35
K (n  1) COn1
Y la ecuación a graficar:
100 − ∅ 𝑛−1
−1
100
log 𝑡∅ = log + 1 − 𝑛 log 𝐶0
𝑘 𝑛−1
EL MÉTODO DEL EXCESO
Cuando una reacción implica muchos reactivos a menudo es posible añadir

cantidades en exceso de todos los reactivos menos uno de ellos. En estos casos la
reacción dependerá únicamente del reactivo escaso. Ese es el caso de la
descomposición de sustancias tóxicas (degradación e hidrólisis) y se pueden
representar por:
A  B  producto
Donde A es el compuesto tóxico y B es otra cantidad (Bacterias o hidrógeno) que
participe en la reacción. Se acostumbra a emplear la siguiente ecuación para
describir y modelar la reacción.
𝑑𝐶𝑎
= −𝑘𝐶𝑎 𝐶𝑏
𝑑𝑡
Donde Ca y Cb son las concentraciones de los reactivos
EL MÉTODO DEL EXCESO
Si la concentración inicial de B es Cb0 y es mucho mayor que la de A, la reacción

resultante puede tener un efecto medible sobre A , mientras que B será afectada
mínimamente. Por consiguiente la reacción puede retomarse como:
𝑑𝐶𝑎
= − 𝑘𝐶𝑏 0 𝐶𝑎 = −𝑘𝑏2 𝐶𝑎
𝑑𝑡
Donde, KCbo  Kb2 se denomina como coeficiente de pseudo-primer orden.
Las técnicas descritas se pueden usar para evaluar los coeficientes.

EL MÉTODOS NUMÉRICOS
Existen métodos orientados al uso del computador para obtener los valores de las
tasas. El método integral de mínimos cuadrados presenta los beneficios de método
integral y diferencial en un mismo método.
En este método se suponen valores de n y k , la ecuación se resuelve para C(t). En
este caso la solución se resuelve por métodos numéricos. La solución consiste en
elaborar una tabla de concentraciones previstas que corresponden a las
concentraciones medidas. Se pueden hallar los cuadrados de las diferencias entre
las medidas y las calculadas y su suma total. Luego utilizando técnicas de
optimización, que comúnmente tienen las hojas de cálculo, se pueden obtener los
coeficientes k y el orden n de la cinética de la reacción.
EFECTO DE LA TEMPERATURA
Las tasas de la mayoría de las reacciones en aguas naturales aumentan con la

temperatura, una regla general es que la velocidad aumenta aproximadamente el
doble para un aumento de temperatura de 10°C
La ecuación de Arrhenius presenta de una mejor forma la dependencia de las
constantes con la temperatura.
𝐸
−𝑅𝑇
𝑘 𝑇𝑎 = 𝐴𝑒 𝑎
Donde:
A es el factor de frecuencia pre exponencial
E la energía de activación en j/mol
R es la constante 8.314 j/(mol K) de los gases perfectos
Ta es la temperatura absoluta de k
La ecuación anterior se usa para comparar la constante de la velocidad de reacción

a dos temperaturas diferentes, usando una proporción de las tasas así:
𝐸 𝑇𝑞2 −𝑇𝑞1
𝑘 𝑇𝑎2 𝑅𝑇𝑎2 𝑇𝑞1
=𝑒
𝑘 𝑇𝑎1
La ecuación se puede simplificar debido a que:
• La temperatura en los cuerpos de agua varían en un rango estrecho (273 a 315

k) y producto Ta2 y Ta1 es constante.
• La diferencia de temperatura es igual en escala absoluta o relativa
En consecuencia se puede expresar el término siguiente como una constate

𝐸
𝜃 = 𝑒 𝑅𝑇𝑎2𝑇𝑞1
y la ecuación se escribe como:
𝑘 𝑇2 𝑇2 −𝑇1
=𝜃
𝑘 𝑇1
En los modelos de la calidad agua mucha reacción se reporta a 20°C entonces la
ecuación queda:
𝑇−20
𝑘 𝑇 =𝑘 20 𝜃
En la tabla se resumen algunos valores comúnmente utilizados en los modelos de

calidad del agua
𝜽 Q10 Reacción
1,024 1,27 Reaireación del oxígeno
1,047 1,58 DBO degradación
1,066 1,89 Crecimiento del fitoplancton
1,08 2,16 Demanda de oxígeno en sedimentos SOD
Donde Q10 está definida como la relación
𝑘 20 10
𝑄10 = =𝜃
𝑘 10
En la figura se ilustra la dependencia funcional de la temperatura en rangos de

temperatura que presentan de forma común las aguas en la naturaleza.
Ejemplo
Evaluar la dependencia de las reacciones a la temperatura del agua.
Los resultados proporcionados por ensayos de laboratorio para una reacción en el

agua son:
T1 = 4°C K1 = 0.12 d-1
T2 = 16°C K2 = 0,20 d-1
• Evalué 𝜃 para esa reacción

• Determine la tasa de la reacción a T=20°C
Ejemplo
T1 = 4°C K1 = 0.12 d-1
T2 = 16°C K2 = 0,20 d-1
𝑘 𝑇2 𝑇2 −𝑇1
=𝜃
𝑘 𝑇1
log 𝑘 𝑇2 − log 𝑘 𝑇1 = 𝑇2 − 𝑇1 log 𝜃
log 𝑘 𝑇 − log 𝑘 𝑇 log 0.20− log 0.12

2 1
𝜃 = 10 𝑇2 −𝑇1 = 10 16−4 = 1.0435
Ejemplo
T1 = 4°C K1 = 0.12 d-1
T2 = 16°C K2 = 0,20 d-1
𝑘 =𝑘 𝜃 𝑇−20
𝑇 20
𝑘 =𝑘 𝜃 16−20
16 20
0.20
𝑘 20 =
1.0435 16−20
𝑘 20 = 0.2371
COEFICIENTE DE DIFUSIÓN TÉRMICA
Debido a que el efecto de la temperatura es fundamental en las tasas de

decaimiento y en general en el proceso de transporte por advección y en el
proceso de difusión, es necesario modelar los cambios en la temperatura debidos
a las transferencias de calor y a su difusión en el seno del agua. Como se vio en
las ecuaciones diferenciales es necesario establecer el coeficiente de difusión
𝒌
térmica 𝜶 = 𝝆𝑪 donde la conductividad térmica (k), la densidad (ρ) y el calor
𝒑
específica (Cp), propiedades del agua dependen de la temperatura.
En la siguiente tabla y gráficos se presenta el comportamiento de estas

propiedades fundamentales.
T r Cp k  1/ s
°C kg/m3 J/kg.K kW/m.k m2/s s/m2 w/m2
Medido Ajustado Medido Ajustado Medido Ajustado
0 999.98 1000.57127 4217 4217.30514 0.561 0.56601402 0.000560691 1783.51382 315.729433
4 1000 1000.24024 4207 4206.4533 0.569 0.5745202 0.000568863 1757.89146 334.632368
5 999.9 1000.13956 4205 4204.09403 0.571 0.57664202 0.000570977 1751.38275 339.487838
10 999.7 999.528607 4194 4194.19687 0.58 0.58722282 0.000580448 1722.80819 364.562524
15 999.1 998.738402 4185 4187.08358 0.589 0.59775642 0.000590275 1694.12485 391.001231
20 998 997.768947 4182 4182.27351 0.598 0.60824282 0.000600538 1665.17292 418.852553
25 997 996.620242 4180 4179.3355 0.607 0.61868202 0.000610891 1636.95285 448.165932
30 996 995.292287 4178 4177.88785 0.615 0.62907402 0.00062039 1611.888 478.991662
35 994 993.785082 4178 4177.59832 0.623 0.63941882 0.00063068 1585.58968 511.380887
40 992.1 992.098627 4179 4178.18412 0.631 0.64971642 0.00064109 1559.84318 545.385601
45 990.1 990.232922 4180 4179.41195 0.637 0.65996682 0.000649715 1539.13597 581.05865
50 988.1 988.187967 4181 4181.09795 0.644 0.67017002 0.000659546 1516.19337 618.453728
55 985.2 985.963762 4183 4183.10775 0.649 0.68032602 0.000668127 1496.72034 657.625382
60 983.3 983.560307 4185 4185.35642 0.654 0.69043482 0.000676224 1478.79946 698.629008
65 980.4 980.977602 4187 4187.80851 0.659 0.70049642 0.000685209 1459.40886 741.520853
70 977.5 978.215647 4190 4190.47803 0.663 0.71051082 0.000693365 1442.24102 786.358015
75 974.7 975.274442 4193 4193.42844 0.667 0.72047802 0.000701662 1425.1874 833.198441
80 971.8 972.153987 4197 4196.77268 0.67 0.73039802 0.00070919 1410.05857 882.100931
85 968.1 968.854282 4201 4200.67315 0.673 0.74027082 0.000717524 1393.68177 933.125132
90 965.3 965.375327 4206 4205.34171 0.675 0.75009642 0.000724345 1380.55823 986.331545
95 961.5 961.717122 4212 4211.0397 0.677 0.75987482 0.000732136 1365.8656 1041.78152
100 957.9 957.879667 4217 4218.0779 0.679 0.76960602 0.000740012 1351.3298 1099.53726
Medios 998.71 4193.5 0.586 0.594769 0.000588 1702.967 386.427943
Medios
0-25 998.74 4189.2 4189.3967 0.589 0.59770922 0.000590626 1694.0883
Variación de la densidad del agua con la temperatura

1010
1000
990 r = -0.003585T2 - 0.068416T + 1,000.571267

R² = 0.999341
Densidad en kg/m3
980
970
960
950
940
0 20 40 60 80 100 120
Temperatura en °C
Variación del calor específico del agua con la
temperatura
4225 Cp = 0.0000032976T4 - 0.0008057232T3 + 0.0777879090T2 - 3.0114312927T + 4,217.3051440074
R² = 0.9969114032
4220
4215
4210
4205
4200
4195
4190
4185
4180
4175
0 20 40 60 80 100 120
Variación de la conductividad térmica del agua con la
temperatura
0.8
0.7
0.6
0.5
y = -0.00000944x2 + 0.00213032x + 0.56014021
0.4 R² = 0.99981945
0.3
0.2
0.1
0
0 20 40 60 80 100 120
Coeficiente de difusividad térmica del agua

0.0008
0.0007
Coeficiente de difusibidad Térmica α (m2/s)
0.0006
y = 0.0000018090x + 0.0005647840
R² = 0.9971326932
0.0005
0.0004
0.0003
0.0002
0.0001
0
0 20 40 60 80 100 120
Temperatura °C
DETERMINACIÓN COEFICIENTE DE DIFUSIÓN MOLECULAR
Como ya se ha presentado, el transporte de sustancias en un río depende de

procesos advectivos, reactivos y difusivos. Los términos advectivos están
determinados por los parámetros del flujo, los cuales pueden ser obtenidos por
modernos métodos de adquisición de datos (aforos en tiempo real en secciones
calibradas etc.). Los procesos reactivos están determinados por las tasas y la
cinética de las reacciones de las diferentes sustancias reactivas transportadas, que
como ya vimos, estas se determinan usando ensayos de laboratorio y técnicas
matemáticas. La parte difusiva depende de parámetros de difícil cuantificacióny
más aún en tiempo real.
Los coeficientes de difusión molecular o dispersión , que se puede suponer (para
simplificar) que es isotrópico, depende del soluto, del solvente, de la temperatura,
y de aspectos hidrodinámicos como la turbulencia del flujo, está última muy
determinada por aspectos hidrodinámicos como la magnitud del flujo y de las
condiciones de contorno (fuerzas de inercia y de corte).
Para su determinación existen diferentes métodos como:

• Formulas empíricas
• Ajustes matemáticos
• Solución de expresiones analíticas
• Ensayos de trazadores
Las fórmulas empíricas se usan con propósitos predictivos o para dar marco de
referencia.
El método de ajustes matemático requiere de una muy elevada cantidad de
información de campo.
El método de expresiones analíticas fue propuesto por Fischer a partir del análisis
del balance de masa por advección neta vs la difusión transversal para flujo
uniforme; este método requiere de una cantidad detallada de información
hidrodinámica de campo.
El método de los trazadores requiere de un alto volumen de trabajo de campo y
laboratorio, así como un proceso de análisi matemático.
Cuando un contaminante es descargado en un rio lo que sucede puede ser

dividido en tres etapas. En una etapa inicial (A) la cantidad de movimiento y la
flotabilidad de la descarga determinan la tasa de dilución del contaminante.
Cuando los efectos de la cantidad de movimiento y la flotabilidad en los procesos
de mezcla se vuelven despreciables, se entra en la segunda etapa (B), en la cual el
contaminante es mezclado a lo ancho y profundo del rio por la turbulencia propia
del flujo. En ríos de baja profundidad este proceso se realiza en muy poco tiempo.
Finalmente cuando el contaminante está completamente mezclado en la

transversal y la vertical, entra en juego el proceso de dispersión por corte, el cual
tendera a eliminar cualquier variación de concentración longitudinal (etapa C).
Cabe destacar que esta última etapa pertenece a lo que se denomina campo
lejano mientras que las dos primeras pertenecen al campo cercano.
Dado que la dispersión por corte es la que prevalece en grandes longitudes, es la

condición que más se ha estudiado y es la condición para la cual se busca el
coeficiente de difusión molecular.
La primer alternativa es estudiar diferentes situaciones reales y obtener el
coeficiente de dispersión longitudinal, para relacionarlo con los parámetros
globales del flujo (velocidad media, ancho superficial, profundidad), los cuales son
fáciles de determinar. Esta forma empírica de obtener el coeficiente de difusión
molecular busca extrapolar las condiciones estudiadas para aplicarlas a otras
situaciones de flujo en ríos similares. La forma general que presentan estas
fórmulas empíricas es:
𝑏
𝐷 𝑈 𝐵 𝑐
=𝑎
ℎ𝑣∗ 𝑣∗ ℎ
Donde:
D es el coeficiente de difusión molecular, ℎ es la profundidad media del flujo, U es
la velocidad media del flujo, B el ancho superficial y 𝑣∗ es la velocidad de corte. a,b
y c son los coeficientes empíricos de ajuste.
En la tabla se presentan los valores obtenidos para la formula empírica por cada
uno de los autores recopilados por Seo y Baek, 2004
𝑏 𝑐
𝐷 𝑈 𝐵
=𝑎
Formula a b c
Elder (1959) 5.920 0 0
Fischer (1979) 0.011 2 2
Liu (1977) 0.180 0.5 2
Iwasa y Aya (1991) 2 0 1.5
Seo y Cheong (1998) 5.920 1.43 0.62
Koussis y Rodríguez-Mirasol (1998) 0.60 0 2
Para un tramo de río con las siguientes características del flujo, en un tramo en
estudio:
B = 13.8m, U=0.6 m/s, ℎ=0.65m y Sw=0.0032, se tiene:
13.8∗0.65
𝑣∗ = 9.806 ∗ 13.8+2∗.65 ∗ 0.0032=0.1365 m/s
Formula a b c D (m2/s)
Elder (1959) 5.92 0 0 53.1
Fischer (1979) 0.011 2 2 859.3
Liu (1977) 0.18 0.5 2 1525.8
Iwasa y Aya (1991) 2 0 1.5 1755.0
Seo y Cheong (1998) 5.92 1.43 0.62 2933.3
Koussis y Rodríguez-Mirasol (1998) 0.6 0 2 2425.9
𝑏 𝑐
𝐷 𝑈 𝐵
=𝑎
La ecuación por la solución de expresiones analíticas de Fischer es:
𝐵 𝑦 𝑦
1 1
𝐷=− 𝑢´ 𝑦 ℎ 𝑦 𝑢´ 𝑦 ℎ 𝑦 𝑑𝑦𝑑𝑦𝑑𝑦
𝐴 0 0 𝜀𝑡 ℎ 𝑦 𝑜
Donde A es el área de la sección transversal de río, u´(y) es la diferencia entre la

velocidad longitudinal en la progresiva transversal con respecto a la velocidad
media en la sección y εt es el coeficiente de mezcla transversal que se calcula con
la ecuación:
𝜀𝑡 = 0.6𝐻𝑣∗
En la que H es la profundidad media del flujo
h(y) es h para u´(y).
Método de los trazadores.

Consiste en realizar inyecciones de sustancias conservativas en el tramo de río que
se piensa estudiar y monitorear en tiempo real el comportamiento espacial y
temporal de la difusión advectiva y difusiva de la concentración de la sustancia en
el campo de flujo. Es el método más apropiado para determinar el coeficiente de
difusión molecular ya que el trazador se transporta con la hidrodinámica del
sistema, representando fielmente su propagación en el tramo en estudio. El
método requiere solucionar la ecuación de “variación de los momentos”, para la
cual se necesitan distribuciones de concentración-tiempo en al menos dos
secciones del tramo en estudio. El tramo debe estar en una zona del campo lejano
en donde predominantemente la variación en la concentración del trazador este
determinada básicamente por la difusión. Se espera que la masa del trazador al
difundirse en el seno del flujo presente una nube cada vez mas dispersa en el
espacio y su concentración tienda a ser más baja buscando un valor constante y
por lo tanto la campana de difusión presente picos atenuados y una base cada vez
mayor, condición que estará dada por la difusión del trazador.

Con la inyección del trazador se inicia el proceso de difusión donde las
concentraciones altas iniciales disminuyen al dispersarse en el campo de flujo. Se
debe tener presente que la inyección de masa es constante y el principio de
conservación de la masa se debe cumplir.

La masa inyectada se difunde al moverse el flujo en dirección aguas abajo,
ampliando el área donde la masa está presente, disminuyendo la concentración al
incrementarse la masa de fluido con contenido del trazador. Si se toma como un
distribución normalizada, la difusión y en consecuencia el coeficiente de difusión
será función de las varianzas de las distribuciones temporales de C vs t en las
secciones i e i+1 y de la velocidad media del flujo.

El método de solución de expresiones analíticas fue propuesto por Fischer a partir
del análisis del balance de masa por advección neta vs difusión transversal para
flujo uniforme (Fischer et al., 1979) , en donde el coeficiente de difusión se
determina como:
𝑈 2 𝜎𝑡2𝑖 − 𝜎𝑡2𝑖−1
𝐷=
2 𝑡𝑖 − 𝑡𝑖−1
Donde:
U es la velocidad media en el tramo i – i-1
𝑡𝑖 e 𝑡𝑖−1 son los tiempos medios de paso de la nube de trazadores en la sección.
𝜎𝑡2𝑖 y 𝜎𝑡2𝑖−1 son las varianzas de las distribuciones temporales de C vs t en las
respectivas secciones i e i-1.
Una tarea previa fundamental es obtener en el laboratorio la curva de calibración

de la sonda medidora de conductividad para diferentes concentraciones del
trazador. Esta labor se realiza preparando muestras con diferentes concentraciones
de trazador y midiendo la conductividad para una temperatura constante de la
mezcla. Los resultados se grafican en un plano Concentración Vs Conductividad.
Los resultados se ajustan a una función.
Además se debe verificar la calibración del equipo medidor de velocidad. En el

laboratorio se hace circular el sensor a una velocidad constante en un medio fluido
en reposo, se verifica que la velocidad registrada corresponda a la velocidad del
ensayo.
Otra tarea fundamental previa es realizar la modelación hidráulica para estimar la
velocidad media del flujo y el tiempo de viaje de la nube de trazador en el seno del
flujo.
Trabajo de campo
Se establece una velocidad aproximada del flujo con base en la modelación y una
medida de velocidad con el equipo. Con esta información se establece el tiempo
medio de viaje entre los puntos de inyección y los puntos de muestreo, se mide el
ancho superficial, la profundidad media del flujo y la velocidad de corte. Se
determina la longitud de mezcla como:
𝐵2
𝐿𝑚 = 𝑈
𝑒𝑡
𝑇ℎ
𝛾 𝑆0
𝑇+2ℎ
Con 𝑒𝑡 = 0.6 𝑢∗ ℎ , donde 𝑢∗ = 𝜌
Se seleccionan el punto de inyección a una distancia mayor a la longitud de mezcla

determinada anteriormente.
Trabajo de campo
Se prepara el trazador diluyendo la masa seleccionada en el volumen de agua

adecuado teniendo en cuenta la concentración de saturación y el caudal esperado.
Se inyecta el trazador de la forma lo más rápida posible para tratar de considerar la
inyección instantánea. Se activan simultáneamente las sondas de medición de
conductividad y se registran por un periodo de tiempo superior al tiempo de viaje
esperado.
La siguiente tarea consiste en monitorear la concentración de oxígeno disuelto en
los puntos de muestreo para establecer posibles variaciones de esta con el tiempo
y el espacio.
Finalmente se procede a tomar muestras del agua para llevarlas al laboratorio y
establecer las tasas de decaimiento de los parámetros que se desean modelar.
Trabajo de campo
Se afora el caudal en las secciones de estudio procediendo a definir verticales a lo

ancho de la sección y midiendo la velocidad a diferentes profundidades
DETERMINACIÓN COEFICIENTE DE REAIREACIÓN
Se denomina reaireación a la renovación del oxígeno o cualquier otro componente

gaseoso de la atmósfera en el seno del agua que fluye en un río. La reaireación es
un fenómeno también complejo que y los parámetros que lo controlan son muy
diferentes a los de la desoxigenación. Los estudios se iniciaron con Aderney y
Beacherf (1919), quienes demostraron que la tasa de reaireación es proporcional
al déficit de oxigeno del agua.
𝑑𝐷
= −𝑘2 𝐷
𝑑𝑡
Entre otros la tasa o constante de reaireación k2 depende de la temperatura y la
turbiedad del agua. El déficit de oxígeno es la diferencia entre la concentración de
saturación y la concentración de oxígeno. Es importante mencionar que la
transferencia se da en la interfase agua-aire y se difunde por advección y difusión
dentro del cuerpo de fluido líquido.
Para determinar la tasa de reaireación se han desarrollado diferentes ecuaciones

que relacionan la tasa con diferentes parámetros del flujo. A continuación se
presenta una serie de estas ecuaciones.
Ecuación de Texas 1.923 𝑣 0.273

𝑘2 =
𝐻0.894
Donde:
v es la velocidad media del flujo (m/s)
H es la profundidad media del flujo (m)
La ecuación puede usarse en ríos con profundidades someras (0.2m a 1.0m) y con
velocidades bajas (menores a 0.3 m/s.
Ecuación de O´Connor y Dobbins

(𝐷𝑙 𝑣)0.5
𝑘2 = 𝐻 0.5
Donde:

Dl es el coeficiente de difusión molecular del oxígeno en agua, igual 1,76 x 10-4
m2/d a 20°C.
La variación del coeficiente de difusión molecular del oxígeno en el agua por la
ecuación de Arrhenius es:
𝐷𝑙 𝑇 = 1.076𝑥10−4 1.037 𝑇−20
Con T la temperatura en °C.

Ecuación de Churchill
5.03 𝑣 0.969
𝑘2 = 𝐻 1.697
Ecuación de Owens
5.34 𝑣 0.67
𝑘2 = 𝐻 1.85
Donde:

Ecuación de Langbein y Durum

2.2 𝑣
𝑘2 =
𝐻 1.33
Ecuación de Tsivogiou y Wallace

𝑘2(20) = 𝐶 𝑆 𝑣
Donde:
v es la velocidad media del flujo (m/d)
C es el coeficiente de escape (m-1)
S es la pendiente de la línea de energía
El valor del coeficiente de escape es

C = 0.177 m-1 a 20°C para caudales de entre 0.42 y 85 m3/s
C = 0.361 m-1 a 20°C para caudales de entre 0.028 y 0.42 m3/s

Presentación Clase 4

Cargado por

Copyright:

Formatos disponibles

Presentación Clase 4

Cargado por

Información del documento

Descripción 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

Presentación Clase 4

Cargado por

Copyright:

Formatos disponibles

MAESTRÍA EN INGENIERÍA CIVIL

Reacciones heterogéneas: Involucran más de una

Reacciones Homogéneas: Involucra una sola fase (o

Las potencias a las que se elevan las concentraciones se

Donde C es la concentración del reactante y n es el

En conclusión si se presenta una tasa de decrecimiento constante a partir de una

• Para la cinética de orden uno:

Separando variable e integrando

En esa ecuación se refleja que el agotamiento es exponencial, es decir que la

En esta ecuación la tasa de decaimiento tiene unidades de tiempo a la menos

Este método consiste básicamente en suponer un valor para n e integral la

Luego se emplean métodos gráficos para establecer si el modelo se ajusta

Los métodos gráficos son la realización de versiones linealizadas de la

Unidades Variable Variable Pendien

Segundo (n = 2 L3/MT-1 1/C t 1/Co K

General(n ≠ 1) (𝐿3 𝑀 −1 )𝑛−1 Cn-1 t C0n-1 (n – 1)K

Se basa en aplicar logaritmos para transformar la ecuación

En consecuencia si se grafica en un plano logarítmico se obtendrá una recta

Existen reacciones en las que se presentan complicaciones a lo largo del tiempo.

En donde la pendiente dará el orden y el intercepto el log k 

La vida media de una reacción es el tiempo que tarda la concentración en caer

Donde t50 es la vida media. Nuevamente se utiliza, como modelo de tasa, la

En donde la pendiente dará el orden y el intercepto el valor de log k 

Si C  Co en t  O que puede ser integrado como:

Combinando las ecuaciones resulta que

Al tomar el logarítmo de la ecuación

La grafica del logaritmo vida media contra el logaritmo de la concentración inicial

Es claro que la selección de la vida media es arbitrario puede haberse seleccionado

EL MÉTODO DEL EXCESO

Cuando una reacción implica muchos reactivos a menudo es posible añadir

EL MÉTODO DEL EXCESO

Si la concentración inicial de B es Cb0 y es mucho mayor que la de A, la reacción

Las técnicas descritas se pueden usar para evaluar los coeficientes.

Las tasas de la mayoría de las reacciones en aguas naturales aumentan con la

La ecuación anterior se usa para comparar la constante de la velocidad de reacción

La ecuación se puede simplificar debido a que:

• La temperatura en los cuerpos de agua varían en un rango estrecho (273 a 315

En consecuencia se puede expresar el término siguiente como una constate

y la ecuación se escribe como:

En la tabla se resumen algunos valores comúnmente utilizados en los modelos de

Donde Q10 está definida como la relación

En la figura se ilustra la dependencia funcional de la temperatura en rangos de

Evaluar la dependencia de las reacciones a la temperatura del agua.

Los resultados proporcionados por ensayos de laboratorio para una reacción en el

• Evalué 𝜃 para esa reacción

log 𝑘 𝑇2 − log 𝑘 𝑇1 = 𝑇2 − 𝑇1 log 𝜃

log 𝑘 𝑇 − log 𝑘 𝑇 log 0.20− log 0.12

Debido a que el efecto de la temperatura es fundamental en las tasas de

En la siguiente tabla y gráficos se presenta el comportamiento de estas

Medios 998.71 4193.5 0.586 0.594769 0.000588 1702.967 386.427943

Variación de la densidad del agua con la temperatura

990 r = -0.003585T2 - 0.068416T + 1,000.571267

Coeficiente de difusividad térmica del agua

Como ya se ha presentado, el transporte de sustancias en un río depende de

Para su determinación existen diferentes métodos como:

Cuando un contaminante es descargado en un rio lo que sucede puede ser

Finalmente cuando el contaminante está completamente mezclado en la

Dado que la dispersión por corte es la que prevalece en grandes longitudes, es la

La ecuación por la solución de expresiones analíticas de Fischer es: