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

Technische Samenvatting: Oplossing van de Advectie-Diffusievergelijking in Twee Dimensies met behulp van de Depositiesnelheid

Probleemstelling
Deze studie behandelt de analytische oplossing van de tweedimensionale stationaire advectie-diffusievergelijking (ADE) onder instabiele atmosferische condities. Het primaire doel is het modelleren van de dispersie van verontreinigingen, waarbij specifiek wordt gefocust op dwarswind-geïntegreerde genormaliseerde concentraties op de centerline. Het model incorporeert een depositiesnelheid randvoorwaarde op grondniveau en maakt gebruik van machtswetfuncties voor zowel de verticale eddy-diffusiviteit als de windsnelheid als functie van de verticale hoogte. De studie streeft ernaar deze analytische benadering te valideren tegen geobserveerde zwavelhexafluoride (SF6) concentratiegegevens van het Copenhagen veldeperiment in Denemarken.

Methodologie
De auteurs maken gebruik van de scheiding van variabelen techniek om de leidende ADE op te lossen. Het wiskundige kader wordt als volgt gedefinieerd:

• Leidende Vergelijking: De tweedimensionale ADE wordt geformuleerd waarbij de windsnelheid $u(z)$ en de verticale eddy-diffusiviteit $k_z(z)$ worden uitgedrukt als machtswetten van de verticale hoogte ($z^p$ en $z^q$, respectievelijk).
• Randvoorwaarden:
  1. Bovenste Rand (Menghoogte $h$): Voorwaarde van nul flux ($k_z \frac{\partial C}{\partial z} = 0$).
  2. Grondrand ($z=0$): Een depositieflux voorwaarde waarbij de flux gelijk is aan het product van de depositiesnelheid ($V_d$) en de concentratie ($k_z \frac{\partial C}{\partial z} = -V_d C$).
  3. Bronterm: Een continue puntbron op schoorsteenhoogte $h_s$, gerepresenteerd door een Dirac-delta functie.
• Analytische Afleiding:
  • De vergelijking wordt gescheiden in ruimtelijke componenten $X(x)$ en $Z(z)$.
  • De verticale component leidt tot een Sturm-Liouville type differentiaalvergelijking met variabele coëfficiënten.
  • Door specifieke waarden voor de exponenten aan te nemen (bijv. $q=0.85$) en gebruik te maken van een variabele transformatie, wordt de vergelijking gereduceerd tot een gegeneraliseerde Bessel-vergelijking.
  • De oplossing voor het verticale concentratieprofiel wordt afgeleid met behulp van Bessel-functies van de eerste en tweede soort, die vereenvoudigd worden naar trigonometrische vormen (sinus en cosinus) voor Bessel-functies van de halve orde.
  • Eigenwaarden en eigenfuncties worden bepaald door de randvoorwaarden toe te passen, wat leidt tot een reeksoplossing voor de concentratie $C_y(x, z)$.
  • Coëfficiënten worden berekend met behulp van orthogonaliteitscondities afgeleid van het windsnelheidsprofiel.

Belangrijkste Bijdragen

• Analytisch Kader: Het artikel presenteert een gesloten analytische oplossing voor de 2D ADE onder instabiele condities, waarbij expliciet een depositiesnelheid randvoorwaarde wordt opgenomen, wat vaak een complexe factor is in dispersiemodellering.
• Machtswet-Parametrisering: Het model maakt gebruik van machtswet-representaties voor windsnelheid en eddy-diffusiviteit, wat flexibiliteit biedt over verschillende stabiliteitsklassen (gedefinieerd door exponenten $p$ en $q$).
• Validatie: Het analytische model wordt rigoureus getest tegen de bekende Copenhagen SF6-dataset, die diverse stabiliteitsklassen (A, B, C, D) en downwind afstanden beslaat.

Resultaten
De studie vergelijkt de voorspelde dwarswind-geïntegreerde genormaliseerde concentraties met de geobserveerde SF6-concentraties uit het Copenhagen-experiment.

• Gegevensvergelijking: Tabel 2 en Figuur 1 presenteren een directe vergelijking van geobserveerde versus voorspelde concentraties over meerdere runs (Run 1 tot en met Run 9) bij variërende downwind afstanden.
• Statistische Prestaties: De prestaties van het model worden geëvalueerd met standaard statistische metrieken (Hanna, 1989):
  • Genormaliseerde Gemiddelde Kwadratische Fout (NMSE): 0,07 (wat wijst op een lage fout).
  • Fractionele Bias (FB): 0,003 (wat wijst op een verwaarloosbare bias).
  • Correlatiecoëfficiënt (COR): 0,82.
  • Fractie binnen een factor twee (FAC2): 0,997.
• Visuele Overeenkomst: Figuur 1 illustreert een hoge mate van overeenstemming tussen de geobserveerde en voorspelde concentratieprofielen. Figuur 2 bevestigt dat de voorspelde waarden binnen een factor twee van de geobserveerde gegevens vallen.

Betekenis en Claims
De auteurs beweren dat het voorgestelde analytische model een "hoge mate van gesloten overeenstemming" bereikt met de geobserveerde veldgegevens. Specifiek stelt het artikel dat de voorspelde concentraties een overeenstemming van 100% bereikten met de geobserveerde concentraties op de centerline, zoals aangetoond door de FAC2-waarde van 0,997 en de bijna nul-bias. De studie concludeert dat de scheidingsmethode, gecombineerd met machtswet-parametriseringen voor wind en diffusiviteit, een robuust analytisch instrument biedt voor het voorspellen van de dispersie van verontreinigingen onder instabiele atmosferische condities zonder de noodzaak voor complexe numerieke simulaties. Het werk bevestigt het nut van transformatie- en scheidingsgebaseerde analytische oplossingen voor het vatten van de essentiële fysica van vervuilingstransport in instabiele grenslagen.