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

Technische Zusammenfassung: Lösung der Advektions-Diffusions-Gleichung in zwei Dimensionen unter Verwendung der Deponierungsgeschwindigkeit

Problemstellung
Diese Studie befasst sich mit der analytischen Lösung der stationären zweidimensionalen Advektions-Diffusions-Gleichung (ADE) unter instabilen atmosphärischen Bedingungen. Das primäre Ziel ist die Modellierung der Schadstoffdispersion, wobei der Schwerpunkt auf den querwindintegrierten normalisierten Konzentrationen an der Mittellinie liegt. Das Modell beinhaltet eine Randbedingung der Deponierungsgeschwindigkeit am Boden und verwendet Potenzfunktionen sowohl für die vertikale Eddy-Diffusivität als auch für die Windgeschwindigkeit als Funktionen der vertikalen Höhe. Die Studie zielt darauf ab, diesen analytischen Ansatz gegen beobachtete Schwefelhexafluorid-Konzentrationsdaten (SF6) aus dem Kopenhagener Feldeperiment in Dänemark zu validieren.

Methodik
Die Autoren verwenden die Methode der Variablentrennung, um die zugrunde liegende ADE zu lösen. Der mathematische Rahmen ist wie folgt definiert:

• Zugrundeliegende Gleichung: Die zweidimensionale ADE wird formuliert, wobei die Windgeschwindigkeit $u(z)$ und die vertikale Eddy-Diffusivität $k_z(z)$ als Potenzgesetze der vertikalen Höhe ($z^p$ bzw. $z^q$) ausgedrückt werden.
• Randbedingungen:
  1. Obere Randbedingung (Mischungsschichthöhe $h$): Nullfluss-Bedingung ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Boden-Randbedingung ($z=0$): Eine Deponierungsfluss-Bedingung, bei der der Fluss dem Produkt aus der Deponierungsgeschwindigkeit ($V_d$) und der Konzentration entspricht ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Quellterm: Eine kontinuierliche Punktquelle in Schornsteinhöhe $h_s$, dargestellt durch eine Dirac-Delta-Funktion.
• Analytische Herleitung:
  • Die Gleichung wird in räumliche Komponenten $X(x)$ und $Z(z)$ zerlegt.
  • Die vertikale Komponente führt zu einer Sturm-Liouville-artigen Differentialgleichung mit variablen Koeffizienten.
  • Durch die Annahme spezifischer Werte für die Exponenten (z. B. $q=0,85$) und die Nutzung einer Variablentransformation wird die Gleichung auf eine verallgemeinerte Bessel-Gleichung reduziert.
  • Die Lösung für das vertikale Konzentrationsprofil wird unter Verwendung von Bessel-Funktionen erster und zweiter Art abgeleitet, welche sich für Bessel-Funktionen halber Ordnung zu trigonometrischen Formen (Sinus und Kosinus) vereinfachen lassen.
  • Eigenwerte und Eigenfunktionen werden durch Anwendung der Randbedingungen bestimmt, was zu einer Reihenlösung für die Konzentration $C_y(x, z)$ führt.
  • Koeffizienten werden mithilfe der Orthogonalitätsbedingungen abgeleitet, die aus dem Windgeschwindigkeitsprofil resultieren.

Wesentliche Beiträge

• Analytischer Rahmen: Die Arbeit präsentiert eine geschlossene analytische Lösung für die 2D-ADE unter instabilen Bedingungen, die explizit eine Deponierungsgeschwindigkeits-Randbedingung einbezieht, was oft ein komplexer Faktor in der Dispersionsmodellierung ist.
• Potenzgesetz-Parametrierung: Das Modell nutzt Potenzgesetz-Darstellungen für Windgeschwindigkeit und Eddy-Diffusivität, was Flexibilität über verschiedene Stabilitätsklassen hinweg (definiert durch die Exponenten $p$ und $q$) ermöglicht.
• Validierung: Das analytische Modell wird rigoros gegen den bekannten Kopenhagener SF6-Datensatz getestet, der verschiedene Stabilitätsklassen (A, B, C, D) und stromabwärts gelegene Distanzen abdeckt.

Ergebnisse
Die Studie vergleicht die vorhergesagten querwindintegrierten normalisierten Konzentrationen mit den beobachteten SF6-Konzentrationen aus dem Kopenhagener Experiment.

• Datenvergleich: Tabelle 2 und Abbildung 1 zeigen einen direkten Vergleich der beobachteten gegenüber den vorhergesagten Konzentrationen über mehrere Durchläufe (Durchlauf 1 bis 9) bei variierenden stromabwärts gelegenen Distanzen.
• Statistische Leistungsfähigkeit: Die Leistung des Modells wird anhand der Standard-Statistikmetriken (Hanna, 1989) bewertet:
  • Normalisierter mittlerer quadratischer Fehler (NMSE): 0,07 (was auf einen geringen Fehler hindeutet).
  • Fraktionaler Bias (FB): 0,003 (was auf einen vernachlässigbaren Bias hindeutet).
  • Korrelationskoeffizient (COR): 0,82.
  • Innerhalb eines Faktors zwei liegend (FAC2): 0,997.
• Visuelle Übereinstimmung: Abbildung 1 illustriert eine hohe Übereinstimmung zwischen den beobachteten und den vorhergesagten Konzentrationsprofilen. Abbildung 2 bestätigt, dass die vorhergesagten Werte innerhalb eines Faktors von zwei zu den beobachteten Daten liegen.

Bedeutung und Behauptungen
Die Autoren behaupten, dass das vorgeschlagene analytische Modell eine „hohe Übereinstimmung“ mit den beobachteten Felddaten erreicht. Insbesondere stellt die Arbeit fest, dass die vorhergesagten Konzentrationen eine 100%ige Übereinstimmung mit den beobachteten Konzentrationen an der Mittellinie erreichten, wie der FAC2-Wert von 0,997 und der nahezu Null befindliche Bias belegen. Die Studie kommt zu dem Schluss, dass die Separationsmethode in Kombination mit Potenzgesetz-Parametrisierungen ein robustes analytisches Werkzeug zur Vorhersage der Schadstoffdispersion unter instabilen atmosphärischen Bedingungen bietet, ohne dass komplexe numerische Simulationen erforderlich sind. Die Arbeit unterstreicht die Nützlichkeit von Transformations- und Separations-basierten analytischen Lösungen zur Erfassung der wesentlichen Physik des Schadstofftransports in instabilen Grenzschichten.