Longitudinal dispersion in fluvial environments: new empirical formulations and applications to pollutant transport modeling

Resumo Técnico: Dispersão Longitudinal em Ambientes Fluviais

Definição do Problema
A avaliação confiável do transporte de poluentes e dos processos de mistura é crítica para a gestão da qualidade da água e a conservação ambiental em rios. O coeficiente de dispersão longitudinal ($K$) é um parâmetro fundamental na equação de advecção-dispersão unidimensional, regendo a precisão das simulações de transporte de poluentes. No entanto, as formulações empíricas existentes para $K$ frequentemente dependem de conjuntos de dados que combinam canais estreitos e largos, potencialmente negligenciando os distintos mecanismos de dispersão inerentes a diferentes geometrias de canais. Além disso, nenhuma formulação existente foi demonstrada como capaz de desempenhar consistentemente bem em ambos os conjuntos de dados experimentais e de campo para diferentes cenários de liberação (instantânea vs. contínua). Este estudo aborda a necessidade de formulações específicas por regime e uma análise detalhada de como a geometria do canal influencia a dinâmica de transporte de poluentes.

Metodologia
A pesquisa emprega uma abordagem de duas frentes: o desenvolvimento de novas formulações empíricas e a simulação numérica do transporte de poluentes.

1. Desenvolvimento de Formulação Empírica:

   • O estudo utiliza análise dimensional para expressar o coeficiente de dispersão adimensional ($\tilde{K} = K/uh$) em termos da razão de velocidade de cisalhamento ($\Xi = u_*/u$) e da razão de aspecto do canal ($\beta = w/h$).
   • Uma abordagem baseada em regimes é adotada, classificando os canais como "estreitos" ($\beta \leq 10$) e "largos" ($\beta > 10$) com base em tendências observadas em conjuntos de dados derivados de campo.
   • A análise de regressão foi realizada em 105 gráficos de dados de campo de várias fontes da literatura (ex: Sooky, Godfrey e Frederick, Yotsukura et al.) para determinar os coeficientes e expoentes para a relação de lei de potência $\tilde{K} = \mathcal{C} \Xi^a \beta^b$.
   • Parâmetros distintos foram derivados para canais estreitos e largos para capturar seus comportamentos hidráulicos únicos.

2. Simulação Numérica:

   • A equação de advecção-dispersão unidimensional foi resolvida numericamente utilizando um método de diferenças finitas.
   • O termo de advecção foi discretizado usando o esquema Lax-Wendroff (segunda ordem de precisão), enquanto o termo de difusão utilizou diferença central com passo de tempo Euler progressivo (primeira ordem de precisão).
   • A estabilidade foi garantida via análise de von Neumann, aderindo aos critérios $\lambda_1 + 2\lambda_2 \leq 1$.
   • As simulações foram conduzidas para dois cenários:
     • Liberação Instantânea: Modelada como uma função Gaussiana representando derramamentos acidentais.
     • Liberação Contínua: Modelada com uma concentração de fonte sustentada ao longo de um intervalo de tempo específico.
   • A análise de sensibilidade de grade determinou um passo de tamanho espacial ($\Delta x$) de 1 m para independência de tamanho de grade.

Principais Contribuições

• Novas Formulações Específicas de Regime: O artigo propõe novas equações empíricas para o coeficiente de dispersão longitudinal especificamente adaptadas para canais estreitos e largos, em vez de uma equação unificada.
• Identificação de Parâmetros Governantes: A análise revela que a dispersão longitudinal em canais estreitos é governada principalmente pela razão de aspecto do canal ($\beta$), enquanto em canais largos, é governada pela razão de velocidade de cisalhamento ($\Xi$).
• Definição de Limiar: O estudo estabelece a razão de aspecto do canal de 10 como o limiar para demarcar os regimes de canais estreitos e largos, um valor escolhido para alinhar-se melhor com as tendências de dados de campo do que limiares anteriores (ex: 28 ou 30,6).
• Dinâmica de Transporte Abrangente: O estudo resolve numericamente a equação de advecção-dispersão para examinar a evolução espaço-temporal das concentrações de poluentes para ambas as liberações instantâneas e contínuas através de diferentes geometrias de canais.

Resultos

• Comportamento do Coeficiente de Dispersão:
  • Em canais estreitos, o coeficiente de dispersão adimensional ($\tilde{K}$) aumenta conforme a razão de velocidade de cisalhamento ($\Xi$) aumenta, impulsionado pela dispersão induzida pelo cisalhamento e pelo fenômeno de queda de velocidade.
  • Em canais largos, $\tilde{K}$ diminui conforme $\Xi$ aumenta. Isso é atribuído à dominância da mistura turbulenta transversal, que promove a homogeneização rápida da concentração do poluente através da seção transversal, reduzindo assim o espalhamento longitudinal.
  • Constatou-se que as formulações existentes falham em capturar esses comportamentos contrastantes, frequentemente superestimando ou subestimando $\tilde{K}$ dependendo do regime do canal.
• Dinâmica de Transporte de Poluentes:
  • Canais Estreitos (Instantânea): Os perfis de concentração exibem um pico distinto que atenua conforme a pluma viaja rio abaixo.
  • Canais Estreitos (Contínua): A concentração permanece em magnitude máxima por um intervalo temporal e espacial finito, formando um platô amplo sem decaimento significativo, indicando dominância de advecção sobre a dispersão.
  • Canais Largos (Instantânea): A dispersão aumentada leva a uma atenuação intensificada do pico de concentração em comparação com canais estreitos. Os perfis mostram uma curva de ascensão rápida e uma curva de queda mais longa.
  • Cais Largos (Contínua): Diferente dos canais estreitos, os canais largos exibem um decaimento significativo da concentração de pico na direção de jusante devido à dispersão aumentada.
• Validação:
  • As formulações propostas mostraram forte concordância com observações de campo, com a maioria dos pontos de dados dentro de bandas de erro de $\pm 20\%$.
  • As previsões numéricas para liberação instantânea alinharam-se bem com soluções analíticas clássicas, embora picos ligeiramente menores tenham sido observados devido à difusão numérica.
  • Comparações com dados de traçadores de campo do Rio Severn (Atkinson e Devis, 2000) mostraram concordância geral, embora as observações de campo tenham exibido picos ligeiramente menores e transporte mais rápido, provavelmente devido a fatores não modelados como armazenamento transiente e influxos laterais.

Significância e Alegações
Os autores alegam que este estudo oferece um arcabouço mais confiável para estimar o coeficiente de dispersão longitudinal ao reconhecer que uma única equação universal é insuficiente para todas as geometrias de rios. Ao distinguir entre regimes estreitos e largos, as formulações propostas fornecem uma representação mais precisa da física da dispersão. Além disso, o estudo elucida as distintas dinâmicas de transporte de poluentes sob diferentes cenários de liberação e geometrias de canais, demonstrando que a interação entre advecção e dispersão varia significamente entre canais estreitos e largos. O trabalho serve como uma ferramenta robusta para simular o transporte de poluentes, embora os autores notem limitações relativas à exclusão de curvatura de rios, troca com planícies de inundação, resistência da vegetação e processos de transporte reativo, sugerindo estes como áreas para investigação futura.