UNIVERSIDAD NACIONAL

SANTIAGO ANTUNEZ DE MAYOLO

EPIS –UNASAM
TEMA: Modelo gaussiano (II)
 Elevación de la pluma
• Los gases de escape turbulentos emitidos, se mezclan con el aire
• A esta mezcla en la pluma se le denomina el arrastre, durante el cual la pluma
aumenta su diámetro mientras viaja a sotavento (la parte opuesta a aquella de
donde viene el viento con respecto a un lugar determinado)
• Los gases salen de la fuente con una cierta velocidad  penetran en la
atmósfera con un cierto momento cinético ascendente + en general, salen a
temperatura más alta que la del aire externo  son menos densos que el aire
exterior  flotan en él
• La combinación del momento y la flotabilidad de los gases hace que se eleven
(fenómeno conocido como elevación de la pluma) y permite que los
contaminantes emitidos asciendan a una mayor altura en la atmósfera

• Conforme se elevan a la salida de la chimenea, los gases van perdiendo energía

cinética y además su temperatura se iguala con la del aire ambiental por lo que,
después de un tiempo, ya no flotan en él y son arrastrados por la componente
horizontal del viento que hace que la pluma se incline
•A mayor velocidad del viento, más horizontal será el movimiento de la pluma
•La velocidad del viento aumenta con la distancia al suelo por lo que, medida
que la pluma se eleva, los vientos más fuertes hacen que se incline aún más
•Este proceso persiste, en general, hasta que la pluma parece paralela al suelo
•La distancia donde la pluma parece llana puede encontrarse bastante lejos de
la chimenea (a sotavento)
•La elevación de la pluma debida a su flotabilidad es una función de la
diferencia de temperatura entre la pluma y la atmósfera circundante
• Atmósfera inestable  la flotabilidad de la pluma aumenta al elevarse  la
altura final de la pluma se incrementa
• Atmósfera estable  la flotabilidad de la pluma disminuye a medida que
se eleva
• Atmósfera neutra  constante
La pluma pierde flotabilidad a través del mismo mecanismo que la hace
serpentear, el viento

La mezcla dentro de la pluma arrastra el aire atmosférico

hacia su interior
A mayor velocidad del viento, más rápida será esta mezcla
Formación de "penachos"
Penacho de espiral: condiciones muy
inestables generalmente favorables para la
dispersión
Algunas veces se pueden producir altas
concentraciones momentáneas al nivel del
suelo

Penacho de abanico: condiciones estables

Una inversión impide el movimiento vertical
pero no el horizontal y el penacho se puede
extender varios km a sotavento de la fuente
Ocurren con frecuencia en las primeras
horas de la mañana (inversión por radiación)

Penacho de cono: condiciones neutrales o

ligeramente estables
Mayor probabilidad de producirse entre la
interrupción de una inversión por radiación y
el desarrollo de condiciones diurnas
inestables
Formación de "penachos"

Penacho de flotación: Condiciones

inestables por encima de una inversión

Penacho de fumigación: Se forma justo

debajo de una capa de inversión y
puede producir una grave situación de
contaminación
 Altura de mezcla
En un diagrama adiabático, el punto en el que la porción de aire que se enfría en el
gradiente vertical adiabático seco corta a la "línea" perfil de temperatura ambiental se
conoce como altura de mezcla y es el nivel máximo al que la porción de aire puede
ascender
Cuando no se produce ninguna intersección, la altura de mezcla, H, se puede extender
a mayores alturas en la atmósfera

El aire que se encuentra bajo la altura

de mezcla conforma la capa de
mezclado y mientras más profunda sea
esta capa, mayor será el volumen de
aire disponible para la dispersión de
los contaminantes
En un estrato de inversión térmica (la
temperatura aumenta con la altura) se
da la situación de máxima estabilidad
y equivale a la existencia de una
“barrera” que impide la dispersión de
los contaminantes y determina la
“altura de mezclado” (H)
 En los MD mas completos es necesario indicar el tipo de
atmósfera en que se produce el proceso de dispersión
Tabla 1: Clases de atmósfera según su estabilidad (Turner)

* La noche se define como el período Categoría Estabilidad

desde una hora después de la puesta de sol, A Extremadamente inestable
hasta una hora antes de la salida del mismo B Moderadamente inestable
** Nubosidad expresada en octavos de C Ligeramente inestable
bóveda celeste cubierta por nubes
D Neutra
Para cielos totalmente cubiertos, tanto
para el día como para la noche, debe asumirse E Ligeramente estable
clase de estabilidad D F Moderadamente estable
 Tipo de atmósfera
Clase (1): Cielos despejados, a una altura solar de 60º sobre el horizonte, típico de
una tarde soleada de verano. Atmósfera muy convectiva
Clase (2): Altura solar entre 35º y 60º - Un día de verano con algunas nubes dispersas
Clase (3): Típico de una tarde de otoño soleada, un día de verano con nubes bajas
dispersas o también de un día con cielos despejados y una altura del sol de 15º a 35º
Clase (4,5): Aunque corresponde a una atmósfera nocturna, se puede usar también
para un día de invierno
 Modelo gaussiano para contaminantes que no reaccionan
Suponiendo constantes la tasa de emisión, Q
(masa de contaminante emitida en la unidad
de tiempo) y las condiciones atmosféricas, se
llega a un estado estacionario, en la cual el
penacho adquiere una forma constante en el
tiempo
La concentración de contaminante es
máxima en el eje del penacho, disminuyendo
hacia los bordes (distribución normal o de
Gauss)
Hipótesis fundamental del modelo gaussiano
La concentración de contaminantes en las
direcciones perpendiculares a la del viento
puede ser descrita utilizando una
distribución normal o de Gauss como la de la
figura (campana de Gauss) cuya forma
depende de los parámetros m y s

- µ (valor medio) indica la posición de la campana

(parámetro de centralización)
- σ es el parámetro de dispersión o desviación
estándar
 Cuanto menor sea σ, más concentrados están los valores alrededor de la
media y cuanto mayor sea σ más "aplastada" será la curva

Gaussianas con diferentes medias Gaussianas con medias iguales (µ= 0) y

( µ= -3, 0, -3) e igual dispersión (σ = 1) diferentes dispersiones (σ = 1, 2, 4)
Un modelo gaussiano parte de unas hipótesis y si las
condiciones reales se alejan mucho de ellas, sus
estimaciones se hacen poco precisas
Es útil para estimar la concentración de un contaminante
para distancias ~ 20 km
No sirve para problemas como la lluvia ácida, que implican
cientos de km
El modelo se basa en la resolución de la ecuación de
difusión atmosférica
Aunque el modelo gaussiano se aplica a una fuente
puntual (chimenea), puede ser usado para considerar
fuentes lineales (carreteras), o fuentes superficiales (que
se modelan como un gran número de fuentes puntuales)
El coeficiente de dispersión se mide en metros e indica cuánto se ha dispersado la
masa inicial cuando la pluma alcanzan una distancia dada desde la fuente de emisión
 Altura efectiva o equivalente de la chimenea
Aunque la pluma tiene su origen a una altura h (la de la chimenea), se eleva una altura
adicional Dh, debido a la capacidad de flotación de los gases que salen a mayor temperatura
que la de su entorno atmosférico y a la cantidad de movimiento cuando salen verticalmente
de la chimenea con una velocidad Vs  la pluma aparece como si se originara en una fuente
puntual a una altura mayor, H, llamada altura efectiva o equivalente de la chimenea
H = h + Dh
El modelo se basa en la difusión de la masa del contaminante en las direcciones y, z cuando un
elemento fluido es arrastrado por el viento en la dirección del eje x con una velocidad u

Hipótesis del modelo gaussiano

 Estado estacionario  C  C(t) y u = cte (en el tiempo y en altura)
 La fuente tiene una emisión constante de un contaminante que es conservador (no se
descompone, reacciona o sedimenta)
 El terreno es relativamente plano y no se producen efectos de absorción u otros
 Fuente puntual sin reflexión en el suelo

Concentración de contaminante en un punto de coordenadas (x, y, z) para la emisión de

un foco de altura efectiva H (sin considerar reflexiones en el suelo):

 2
C  x, y, z  =
Q
 2 

exp -
y  
exp -
 z - H  
2 π u σ y σz  2 σ2   2 σ2 
 y  z 

Los coeficiente de dispersión, σy y σz,
σy y σz son los coeficientes de dispersión lateral y vertical forma de la
distribución de concentraciones con la distancia lateral (y) y con la vertical (z)
Se miden en metros e indican cuánto se ha dispersado la masa inicial cuando el
penacho alcanza una distancia dada desde la fuente de emisión
"sin reflexiones“ extraordinariamente importante
La ecuación anterior nos da la concentración en la dirección del viento hasta
llegar a un punto en la dirección x en que la concentración a nivel del suelo (z =
0) sea significativa ya que entonces tendrá lugar una apreciable reflexión del
contaminante gaseoso al difundirse regresivamente a la atmósfera desde el
nivel del suelo
La reflexión es un fenómeno de retrodifusión
de los contaminantes cuando encuentran la
barrera del suelo

Se supone que el suelo

no es un sumidero

es decir, que los contaminantes

no se absorben por lo que
se reflejan volviendo a la
atmósfera

Fuente puntual con reflexión en el suelo

Considerar la reflexión en el suelo es equivalente a considerar dos
fuentes de contaminación, una situada en z = +H y otra situada en
z = –H:
 2  
-  z - H
2 
-  z + H 
2 
C  x, y, z  =
Q
exp  -
y  exp   + exp  
2 π u σ y σz  2σ2    2σ2 
 
 2σ2  
 
 y   z   z 

 A nivel del suelo (z = 0), en la línea central (y =
Las concentraciones a nivel del suelo (z = 0) los receptores reciben los máximos niveles
0) son muy importantes (receptores): de contaminación:
   2   2 
 2 
C  x, y, 0  = exp  -
H 
C  x, 0, 0  = exp  -
Q y Q H 
exp  - 
2    
πuσ σ
y z  2 σy   2 σ2  πuσ σ
y z  2 σ2 
   z  z

σy y σz son función de la posición en la dirección del viento, x, y de la

estabilidad atmosférica (requiere la caracterización del tipo de
atmósfera en una de las categorías de Turner)
Estudiar una fuente elegir la clase de estabilidad atmosférica típica
de la región que conduzca al peor episodio de contaminación posible
A través de numerosas medidas experimentales en la atmósfera, se
ha llegado a obtener la correlación de σy y σz con la distancia y el tipo
de atmósfera
Hay varios métodos para obtener los coeficientes de dispersión σy y
σz: describiremos los métodos de Pasquill (gráfico y analítico) y el
método de Martin (analítico) por ser ampliamente usados en la
bibliografía del tema
 Método Gráfico
Las curvas de Pasquill-
Gifford

Gráficas cuyo
objetivo es la
estimación de los
valores de σy y σz
Los valores de σz
tienen mayor error
que los de σy sobre
todo para distancias
superiores a 1 km en
la dirección del
viento
¡¡ Distancia x km
σy , σz m !!
 Método Gráfico
Curvas de Pasquill - Gifford

Amplitud del penacho, σz,

para una distancia x dada:
σz Máxima

inestabilidad atmosférica
máxima (A)

σz Mínima

atmósfera muy estable (F)

¡¡ Distancia x km
σy , σz m !!
 Determinación de los coeficientes de difusión gaussiana: Métodos analíticos
Debido a la dificultad de leer los valores de σy y σz en las gráficas se han obtenido los
ajustes algebraicos de las mismas (los valores de σ son promedios sobre un intervalo de
10 minutos)
Método de Pasquill:
Corrección de σz por rugosidad del terreno
Pasquill propuso las ecuaciones que se muestran a continuación y en las cuales aparece
una dependencia de un coeficiente de rugosidad del terreno, z0 ,para el cálculo de σz
La rugosidad tiene en cuenta el efecto sobre el coeficiente de dispersión vertical, σz, de
la vegetación exuberante, cultivos, edificios, etc., que cambian la forma vertical del
penacho
El coeficiente de dispersión lateral, σy , no se ve afectado por la rugosidad del terreno
Desviación típica transversal (σy) y vertical (σz) en metros, ajustadas para distancias a la
fuente, x (m), entre 100 m y 10 km, siguen la ley potencial:

σy = a xp σz = b’ xm → x (metros) y σ(metros)!!!

a y p son coeficientes tabulados que dependen solo de la clase de estabilidad

b’ y m dependen de la clase de estabilidad y del coeficiente de rugosidad del terreno:
b’ = b’(z0) y m= m(z0)
Los valores de z0 dependen del tipo de superficie y suelen tomarse los siguientes
Los coeficientes a, p, b’ y m se dan en la tabla siguiente, en
función de la rugosidad del terreno y del tipo de atmósfera y
son ajustes de valores experimentales con validez estadística
Determinación de los coeficientes de difusión gaussiana: Métodos analíticos

Fórmulas de Martin

También se utilizan otras ecuaciones que no hacen corrección de rugosidad del terreno
como la propuesta por D. O. Martin:
σy = a xb σz = c xd + f → x (kilometros) y σ (metros)!!!
Las constantes a, c, d y f dependen de la categoría de estabilidad de Pasquill y b vale
siempre 0.894
 Elevación del penacho
La elevación del penacho, Dh, se define como la diferencia entre la altura de la línea
central final del penacho y la altura inicial de la fuente y es directamente proporcional al
contenido calorífico y a la velocidad de salida del efluente e inversamente proporcional
a la velocidad del viento
Existen varios métodos para determinar la elevación del penacho y una de las fórmulas
más empleadas para el cálculo de esta elevación es la formula de Holland:

Vs d  k Q h 
Δh =  n +   Q h = Q cp  Ts - Ta 
u  Vs d 

Dh = Elevación del penacho por encima de la fuente emisora (m)

Vs = Velocidad de salida del contaminante (m s-1)
d = Diámetro interior del conducto de emisión (m)
u = Velocidad del viento (m s-1)
Ts, Ta = Temperaturas del contaminante y ambiente respectivamente (K)
n = Constante adimensional = 1.5
k = Constante = 0.0096 m2 kJ-1
Qh = Tasa de emisión de calor de la chimenea (kJ s-1)
Q = Tasa de emisión de gas (kg s-1)
cp = Calor específico del gas emitido (kJ kg-1 K-1)
 Los valores de ∆h obtenidos con la fórmula de Holland deben corregirse
multiplicando por un factor, establecido por Pasquill-Gifford-Turner y que es
función de las condiciones meteorológicas

Categorías de Factor de corrección

estabilidad de Dh
A, B 1.15
C 1.10
D 1.00
E, F 0.85

También se utiliza la ecuación de Carson y Moses para el cálculo de esta elevación:

Vs d Qh
Δh = - 0.029 + 2.62
u u
 Concentración máxima en la dirección x, a nivel del suelo y en la línea
central

Para el caso de condiciones inestables a casi neutras, se cumple que la

relación sy/sz es prácticamente constante e independiente de x
En estas condiciones para y = 0 (línea central del penacho) se cumple la
relación:
sz = 0.707 H
y se calcula la concentración máxima a nivel del suelo, en la línea central
y en la dirección del viento a través de la expresión:

0.117 Q
Cmax, reflex =
u σ y σz
 Perfil de velocidades del viento
Si no se dispone del dato de la velocidad del viento a la altura efectiva de la
chimenea, H, sino que solo se conoce la velocidad, uref, a una altura de referencia
href (las medidas estándar de velocidad de viento son a 2 y 10 m de altura), se utiliza
la expresión de la variación del viento con la altura en la atmósfera
p
 H 
u = uref  
h
 ref 
Los valores del coeficiente p como función de la clase de estabilidad y el entorno en
que se mueve el viento son los siguientes:
GRACIAS