Solution of the Advection-Diffusion Equation in Two Dimensions using the Deposition Velocity

Resumo Técnico: Solução da Equação de Advecção-Difusão em Duas Dimensões utilizando a Velocidade de Deposição

Definição do Problema
Este estudo aborda a solução analítica da equação de advecção-difusão (ADE) estacionária bidimensional sob condições atmosféricas instáveis. O objetivo principal é modelar a dispersão de poluentes, focando especificamente em concentrações normalizadas integradas na direção transversal ao vento na linha de centro. O modelo incorpora uma condição de contorno de velocidade de deposição no nível do solo e utiliza funções de potência tanto para a difusividade de turbulência vertical quanto para a velocidade do vento como funções da altura vertical. O estudo busca validar esta abordagem analítica contra dados observados de concentração de hexafluoreto de enxofre (SF6) do experimento de campo de Copenhague, na Dinamarca.

Metodologia
Os autores empregam a técnica de separação de variáveis para resolver a ADE governante. O arcabouço matemático é definido como segue:

• Equação Governante: A ADE bidimensional é formulada onde a velocidade do vento $u(z)$ e a difusividade de turbulência vertical $k_z(z)$ são expressas como leis de potência da altura vertical ($z^p$ e $z^q$, respectivamente).
• Condições de Contorno:
  1. Limite Superior (Altura de Mistura $h$): Condição de fluxo zero ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Limite do Solo ($z=0$): Uma condição de fluxo de deposição onde o fluxo é igual ao produto da velocidade de deposição ($V_d$) e a concentração ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Termo de Fonte: Uma fonte pontual contínua na altura da chaminé $h_s$ representada por uma função delta de Dirac.
• Derivação Analítica:
  • A equação é separada em componentes espaciais $X(x)$ e $Z(z)$.
  • O componente vertical leva a uma equação diferencial do tipo Sturm-Liouville com coeficientes variáveis.
  • Ao assumir valores específicos para os expoentes (por exemplo, $q=0,85$) e utilizar uma transformação de variável, a equação é reduzida a uma equação de Bessel generalizada.
  • A solução para o perfil de concentração vertical é derivada usando funções de Bessel de primeira e segunda espécie, que são simplificadas para formas trigonométricas (seno e cosseno) para funções de Bessel de ordem seminteira.
  • Autovalores e autofunções são determinados pela aplicação das condições de contorno, levando a uma solução em série para a concentração $C_y(x, z)$.
  • Os coeficientes são calculados usando condições de ortogonalidade derivadas do perfil de velocidade do vento.

Principais Contribuições

• Arcabouço Analítico: O artigo apresenta uma solução analítica de forma fechada para a ADE 2D sob condições instáveis, incorporando explicitamente uma condição de contorno de velocidade de deposição, que é frequentemente um fator complexo na modelagem de dispersão.
• Parametrização por Lei de Potência: O modelo utiliza representações de lei de potência para a velocidade do vento e a difusividade de turbulência, permitindo flexibilidade entre diferentes classes de estabilidade (definidas pelos expoentes $p$ e $q$).
• Validação: O modelo analítico é rigorosamente testado contra o conhecido conjunto de dados de SF6 de Copenhague, cobrindo várias classes de estabilidade (A, B, C, D) e distâncias a jusante do vento.

Resultos
O estudo compara as concentrações normalizadas integradas na direção transversal previstas com as concentrações observadas de SF6 do experimento de Copenhague.

• Comparação de Dados: A Tabela 2 e a Figura 1 apresentam uma comparação direta entre as concentrações observadas e previstas em múltiplas execuções (Execução 1 a Execução 9) em várias distâncias a jusante.
• Desempenho Estatístico: O desempenho do modelo é avaliado utilizando métricas estatísticas padrão (Hanna, 1989):
  • Erro Quadrático Médio Normalizado (NMSE): 0,07 (indicando baixo erro).
  • Viés Fracionário (FB): 0,003 (indicando viés desprezível).
  • Coeficiente de Correlação (COR): 0,82.
  • Fração dentro de um Fator de Dois (FAC2): 0,997.
• Concordância Visual: A Figura 1 ilustra um alto grau de concordância entre os perfis de concentração observados e previstos. A Figura 2 confirma que os valores previstos estão dentro de um fator de dois dos dados observados.

Significância e Alegações
Os autores afirmam que o modelo analítico proposto alcança um "alto grau de concordância fechada" com os dados de campo observados. Especificamente, o artigo afirma que as concentrações previstas atingiram 100% de concordância com as concentrações observadas na linha de centro, conforme evidenciado pelo valor de FAC2 de 0,997 e o viés próximo de zero. O estudo conclui que o método de separação, combinado com parametrizações de lei de potência para vento e difusividade, fornece uma ferramenta analítica robusta para prever a dispersão de poluentes sob condições atmosféricas instáveis sem a necessidade de simulações numéricas complexas. O trabalho reforça a utilidade de soluções analíticas baseadas em transformadas e separação para capturar a física essencial do transporte de poluentes em camadas limite instáveis.