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

Résumé Technique : Résolution de l'équation d'advection-diffusion en deux dimensions à l'aide de la vitesse de dépôt

Énoncé du problème
Cette étude traite de la solution analytique de l'équation d'advection-diffusion (ADE) stationnaire en deux dimensions sous des conditions atmosphériques instables. L'objectif principal est de modéliser la dispersion de polluants, en se concentrant spécifiquement sur les concentrations normalisées intégrées sur la largeur du panache au niveau de la ligne centrale. Le modèle intègre une condition aux limites de vitesse de dépôt au niveau du sol et utilise des fonctions de puissance pour la diffusivité turbulente verticale et la vitesse du vent en fonction de la hauteur verticale. L'étude cherche à valider cette approche analytique par rapport aux données de concentration observées de l'hexafluorure de soufre (SF6) issues de l'expérience de terrain de Copenhague, au Danemark.

Méthodologie
Les auteurs emploient la technique de la séparation des variables pour résoudre l'ADE régissante. Le cadre mathématique est défini comme suit :

• Équation régissante : L'ADE en deux dimensions est formulée où la vitesse du vent $u(z)$ et la diffusivité turbulente verticale $k_z(z)$ sont exprimées comme des lois de puissance de la hauteur verticale ($z^p$ et $z^q$, respectivement).
• Conditions aux limites :
  1. Limite supérieure (Hauteur de mélange $h$) : Condition de flux nul ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Limite au sol ($z=0$) : Une condition de flux de dépôt où le flux est égal au produit de la vitesse de dépôt ($V_d$) et de la concentration ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Terme source : Une source ponctuelle continue à la hauteur de la cheminée $h_s$ représentée par une fonction delta de Dirac.
• Dérivation analytique :
  • L'équation est séparée en composantes spatiales $X(x)$ et $Z(z)$.
  • La composante verticale conduit à une équation différentielle de type Sturm-Liouville à coefficients variables.
  • En supposant des valeurs spécifiques pour les exposants (par exemple, $q=0,85$) et en utilisant une transformation de variable, l'équation est réduite à une équation de Bessel généralisée.
  • La solution du profil de concentration verticale est dérivée à l'aide de fonctions de Bessel de première et seconde espèce, qui se simplifient en formes trigonométriques (sinus et cosinus) pour les fonctions de Bessel d'ordre demi-entier.
  • Les valeurs propres et les fonctions propres sont déterminées en appliquant les conditions aux limites, conduisant à une solution en série de la concentration $C_y(x, z)$.
  • Les coefficients sont calculés à l'aide des conditions d'orthogonalité dérivées du profil de vitesse du vent.

Contributions clés

• Cadre analytique : L'article présente une solution analytique sous forme fermée de l'ADE en 2D sous des conditions instables, intégrant explicitement une condition aux limites de vitesse de dépôt, ce qui est souvent un facteur complexe dans la modélisation de la dispersion.
• Paramétrage par loi de puissance : Le modèle utilise des représentations par loi de puissance pour la vitesse du vent et la diffusivité turbulente, permettant une flexibilité à travers différentes classes de stabilité (définies par les exposants $p$ et $q$).
• Validation : Le modèle analytique est rigoureusement testé par rapport au célèbre ensemble de données SF6 de Copenhague, couvrant diverses classes de stabilité (A, B, C, D) et distances en aval.

Résultats
L'étude compare les concentrations normalisées intégrées sur la largeur du panache prédites aux concentrations de SF6 observées.

• Comparaison des données : Le Tableau 2 et la Figure 1 présentent une comparaison directe des concentrations observées et prédites à travers plusieurs séries (Runs 1 à 9) à diverses distances en aval.
• Performance statistique : La performance du modèle est évaluée à l'aide des métriques statistiques standards (Hanna, 1989) :
  • Erreur quadratique moyenne normalisée (NMSE) : 0,07 (indiquant une erreur faible).
  • Biais fractionnaire (FB) : 0,003 (indiquant un biais négligeable).
  • Coefficient de corrélation (COR) : 0,82.
  • Fraction dans un facteur de deux (FAC2) : 0,997.
• Accord visuel : La Figure 1 illustre un haut degré d'accord entre les profils de concentration observés et prédits. La Figure 2 confirme que les valeurs prédites se situent dans un facteur de deux des données observées.

Signification et affirmations
Les auteurs affirment que le modèle analytique proposé atteint un « haut degré d'accord étroit » avec les données de terrain observées. Plus précisément, l'article stipule que les concentrations prédites ont atteint un accord de 100 % avec les concentrations observées à la ligne centrale, comme en témoignent la valeur FAC2 de 0,997 et le biais quasi nul. L'étude conclut que la méthode de séparation, combinée aux paramétrages par loi de puissance pour le vent et la diffusivité, fournit un outil analytique robuste pour prédire la dispersion des polluants sous des conditions atmosphériques instables sans nécessiter de simulations numériques complexes. Ce travail renforce l'utilité des solutions analytiques basées sur les transformations et la séparation pour capturer la physique essentielle du transport de polluants dans les couches limites instables.