Difusão de Poluentes na Atmosfera Aplicações

Aplicação da modelagem de dispersão de poluentes:

™ Calcular se a população está ou não exposta ou se estará exposta a material perigoso Exposição

Depende basicamente de 3 fatores:

™ Quão misturada a pluma está ™ Horizontal ou vertical x tempo

™ Quão próxima a pessoa está do centro da pluma ™ Relação com o levantamento da pluma

™ Quanto tempo ela está dentro da pluma ™ Velocidade do vento, duração da emissão, direção,... Dispersão e Difusão

™ Termos usados indistintamente => processos turbulentos na atmosfera

™ É um processo aleatório (i.e. Monte Carlo)

™ Não isotrópico (e depende do vento)

™ Cria uma distribuição gaussiana (H/V)

™ Sempre age no sentido de diluir as concentrações Dispersão

Escoamento do ar: ™ Quanto mais ar, mais diluída

Turbulência na atmosfera:

™ Mecânica Dispersão

Escoamento do ar: ™ Quanto mais ar, mais diluída

Turbulência na atmosfera:

™ Mecânica

™ Cisalhamento Dispersão

Escoamento do ar: ™ Quanto mais ar, mais diluída

Turbulência na atmosfera:

™ Mecânica ™ Cisalhamento

™ Empuxo Estabilidade

™ O quanto a atmosfera suporta, tolera ou suprime movimentos turbulentos na vertical

Critérios:

™ Instável se dT/dz < Γd = -9.8 K / Km

™ Neutro dT/dz = Γd

™ Estável se dT/dz > Γd Estimativas

™ σθ – desvio padrão da direção do vento horizontal

™ σw – desvio padrão da velocidade vertical Estimativas

™ Calculando a velocidade de fricção u*

Vertical moment flux (a measure of turbulent diffusion)

(1) cisalhamento do vento

Perfil logaritmo do vento na PBL (2) Comprimento de mistura Tipos de pluma -1, 2 Tipos de pluma – 3, 4 Tipos de pluma - 5 “Plume rise” e vento “Plume rise” e vento

Os ventos vão entortar a pluma e mistura-la com o ar ambiente

A pluma continua a subir enquanto Tp > Tamb

O Efeito da velocidade de ejeção é curto Exemplos especiais

™ Containers pressurizados ™ Velocidade de ejeção pode não ser vertical

™ Incêndios ™ Velocidade zero, mas muito empuxo (calor)

™ Gases densos ™ Se comporta como um liquido (segue a topografia) Plume rise Holland’s Eq. “Plume rise” e vento

Durante a manhã

Estabiliza a temperatura mais Previsão de rapidamente ventos fortes a tarde Tipos de modelos de dispersão

Assume umaExemplos de alguns ™ Dispersão gaussiana distribuição normal modelos de dispersão: Difusão proporcional ™ Transporte dos gradientes ao gradienteADMS3 e a cte de difusividadeAERMOD CALPUFF Baseados em ™ Estatísticos DISPERSION21 Monte-Carlo ISC3 MERCURE Monin-Obuhkov ™ Similaridade PUFF-PLUME Tuneis de vento SIRANE HYSPLIT ™ Mistura uniforme STILT “Caixa” instantânea FLEXPART ... Box model Pluma Gaussiana

™ É uma aproximação simples para o processo de

dispersão (válido para t>>t0)

Ex.: 99% da massa Pluma Gaussiana

™ É preciso conhecer a capacidade de dispersão da atmosfera!

™ Método 1: medir σθ e σw Pasquill turbulence types Pasquill turbulence types

Divergem dependendo da instabilidade

logaritmicas Pluma Gaussiana

™ Além disso, a estabilidade suprime os movimentos verticais, mas não horizontais! Geometria da pluma

Distribuições gaussianas em y e z Dispersão vento-abaixo na direção x Ex.: fonte pontual contínua

™ Predição da concentração no centro da pluma Geometria da pluma Ex.: fonte pontual contínua

™ Q = 100 g/s de SO2 ™ U = 10 m/s e insolação moderada (Pasquill=>D)

™ 1km vento abaixo, teremos (no centro da pluma):

™ σy=76m

™ σz=38m Ex.: fonte pontual contínua

™ Fora do centro da pluma, temos a equações completa: Ex.: software Screen3

Pluma Gaussiana

Hipóteses:

™ É uma solução idealizada

™ Mais estabilidade => maiores concentrações

™ Menor estabilidade => menores concentrações

™ Concentrações ~1/vento

Falha quando:

™ Dispersão em longas distâncias Transporte dos Gradientes

™ Equação de conservação:

∂N ! 2 + ∇⋅(NV)− D∇ N = 0 ∂t

Advecção Difusão molecular

™ Decompondo:

χ = χ + χʹ e V = V + Vʹ Transporte dos Gradientes

™ Chegamos em: Difusão turbulenta

∂N ! 2 + ∇⋅(NV)− D∇ N = −∇⋅ NʹVʹ ∂t

™ Em primeira ordem, a turbulência causa um transporte parecido com difusão molecular. Pela lei de Fick: ∇⋅ NʹVʹ = K∇N

™ Um modelo numérico então resolve a eq completa. Turbulência seca

™ Este método é utilizado tanto em modelos numéricos de dispersão de pluma: ™ Hysplit ™ Flexpart ™ STILT

™ Quanto em modelos de previsão de tempo e clima! Modelo de dispersão precisa: • Condições meteorológicas • Circulação, turbulência e convecção, etc.. • Parâmetros de emissão • Parâmetros sobre a superfície • Topografia, rugosidade, etc.. • Reações químicas Plume going forward is diluted! T0-zf2

Plume going T0-embp backwards is concentrated! Average IOP1 @ 10km

Log10 (Normalized Concentration) Longitude= T2

Longitude=T3 Também importante:

™ Transporte por convecção ™ Remoção úmida (partículas) por precipitação ™ Remoção seca por deposição na superfície ™ Terreno