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

तकनीकी सारांश: निक्षेपण वेग (Deposition Velocity) का उपयोग करके दो आयामों में संवहन-विसरण समीकरण (Advection-Diffusion Equation) का समाधान

समस्या विवरण
यह अध्ययन अस्थिर वायुमंडलीय स्थितियों के तहत दो-आयामी स्थिर अवस्था संवहन-विसरण समीकरण (ADE) के विश्लेणात्मक समाधान को संबोधित करता है। प्राथमिक उद्देश्य प्रदूषकों के प्रसार का मॉडल तैयार करना है, जो विशेष रूप से सेंटरलाइन (centerline) पर क्रॉसविंड-एकीकृत सामान्यीकृत सांद्रता (crosswind-integrated normalized concentrations) पर केंद्रित है। मॉडल जमीनी स्तर पर एक निक्षेपण वेग सीमा स्थिति (deposition velocity boundary condition) को शामिल करता है और ऊर्ध्वाधर ऊंचाई के फलन के रूप में ऊर्ध्वाधर एडी डिफ्यूज़िविटी (vertical eddy diffusivity) और पवन गति के लिए पावर-लॉ (power-law) फलनों का उपयोग करता है। यह अध्ययन इस विश्लेणात्मक दृष्टिकोण को डेन्मार्क में कोपेनहेगन फील्ड प्रयोग से प्राप्त सल्फर हेक्साफ्लोराइड (SF6) सांद्रता डेटा के विरुद्ध मान्य करने का प्रयास करता है।

कार्यप्रणाली
लेखक ADE को हल करने के लिए पृथक्करण विधि (separation of variables) का उपयोग करते हैं। गणितीय ढांचा निम्नानुसार परिभाषित है:

• शासी समीकरण (Governing Equation): दो-आयामी ADE को सूत्रबद्ध किया गया है जहाँ पवन गति $u(z)$ और ऊर्ध्वाधर एडी डिफ्यूज़िविटी $k_z(z)$ को ऊर्ध्वाधर ऊंचाई के पावर लॉ ($z^p$ और $z^q$ क्रमशः) के रूप में व्यक्त किया गया है।
• सीमा स्थितियाँ (Boundary Conditions):
  1. शीर्ष सीमा (मिश्रण ऊंचाई $h$): शून्य फ्लक्स स्थिति ($k_z \frac{\partial C}{\partial z} = 0$)।
  2. जमीनी सीमा ($z=0$): एक निक्षेपण फ्लक्स स्थिति जहाँ फ्लक्स निक्षेपण वेग ($V_d$) और सांद्रता का गुणनफल होता है ($k_z \frac{\partial C}{\partial z} = -V_d C$)।
  3. स्रोत पद (Source Term): स्टैक की ऊंचाई $h_s$ पर एक निरंतर बिंदु स्रोत जिसे डिराक डेल्टा फंक्शन (Dirac delta function) द्वारा दर्शाया गया है।
• विश्लेणात्मक व्युत्पत्ति (Analytical Derivation):
  • समीकरण को स्थानिक घटकों $X(x)$ और $Z(z)$ में अलग किया जाता है।
  • ऊर्ध्वाधर घटक एक स्टर्म-लिउविल (Sturm-Liouville) प्रकार के अवकल समीकरण की ओर ले जाता है जिसमें परिवर्तनशील गुणांक होते हैं।
  • विशिष्ट घातांकों (जैसे, $q=0.85$) को मानकर और एक परिवर्तनशील रूपांतरण का उपयोग करके, समीकरण को एक सामान्यीकृत बेसेल समीकरण (generalized Bessel equation) में बदल दिया जाता है।
  • प्रथम और द्वितीय प्रकार के बेसेल फलनों का उपयोग करके ऊर्ध्वाधर सांद्रता प्रोफाइल के लिए समाधान व्युत्पन्न किया जाता है, जिन्हें हाफ-ऑर्डर बेसेल फलनों के लिए त्रिकोणमितीय रूपों (साइन और कोसाइन) में सरल बनाया गया है।
  • आइजनवैल्यू (eigenvalues) और आइजनफंक्शंस (eigenfunctions) को सीमा स्थितियों को लागू करके निर्धारित किया जाता है, जिससे सांद्रता $C_y(x, z)$ के लिए एक श्रेणी समाधान (series solution) प्राप्त होता है।
  • गुणांकों की गणना पवन गति प्रोफाइल से प्राप्त ऑर्थोगोनैलिटी (orthogonality) स्थितियों का उपयोग करके की जाती है।

प्रमुख योगदान

• विश्लेणात्मक ढांचा: शोधपत्र अस्थिर स्थितियों के तहत 2D ADE के लिए एक क्लोज्ड-फॉर्म विश्लेणात्मक समाधान प्रस्तुत करता है, जिसमें स्पष्ट रूप से एक निक्षेपण वेग सीमा स्थिति शामिल है, जो प्रसार मॉडलिंग में अक्सर एक जटिल कारक होता है।
• पावर-लॉ पैरामीट्राइजेशन: मॉडल पवन गति और एडी डिफ्यूज़िविटी के लिए पावर-लॉ निरूपण का उपयोग करता है, जो विभिन्न स्थिरता वर्गों (घातांक $p$ और $q$ द्वारा परिभाषित) के लिए लचीलापन प्रदान करता है।
• सत्यापन: विश्लेणात्मक मॉडल का कोपेनहेगन SF6 डेटासेट के विरुद्ध कठोर परीक्षण किया गया है, जिसमें विभिन्न स्थिरता वर्गों (A, B, C, D) और डाउनविंड दूरियों को कवर किया गया है।

परिणाम
अध्ययन कोपेनहेगन प्रयोग से प्राप्त प्रेक्षित SF6 सांद्रता के साथ अनुमानित क्रॉसविंड-एकीकृत सामान्यीकृत सांद्रता की तुलना करता है।

• डेटा तुलना: तालिका 2 और चित्र 1 विभिन्न रन (रन 1 से रन 9 तक) में विभिन्न डाउनविंड दूरियों पर प्रेक्षित बनाम अनुमानित सांद्रता का सीधा तुलनात्मक विवरण प्रस्तुत करते हैं।
• सांख्यिकीय प्रदर्शन: मॉडल के प्रदर्शन का मूल्यांकन मानक सांख्यिकीय मेट्रिक्स (हन्ना, 1989) का उपयोग करके किया गया है:
  • सामान्यीकृत माध्य वर्ग त्रुटि (NMSE): 0.07 (कम त्रुटि का संकेत)।
  • आंशिक पूर्वाग्रह (Fractional Bias - FB): 0.003 (नगण्य पूर्वाग्रह का संकेत)।
  • सहसंबंध गुणांक (Correlation Coefficient - COR): 0.82।
  • दो के कारक के भीतर अंश (Fraction within a Factor of Two - FAC2): 0.997।
• दृश्य सहमति (Visual Agreement): चित्र 1 प्रेक्षित और अनुमानित सांद्रता प्रोफाइल के बीच उच्च स्तर की सहमति को दर्शाता है। चित्र 2 पुष्टि करता है कि अनुमानित मान प्रेक्षित डेटा के दो के कारक के भीतर आते हैं।

महत्व और दावे
लेखक दावा करते हैं कि प्रस्तावित विश्लेणात्मक मॉडल प्रेक्षित क्षेत्र डेटा के साथ "उच्च स्तर की क्लोज सहमति" प्राप्त करता है। विशेष रूप से, शोधपत्र में कहा गया है कि अनुमानित सांद्रता ने सेंटरलाइन पर प्रेक्षित सांद्रता के साथ 100% सहमति प्राप्त की है, जैसा कि FAC2 मान 0.997 और लगभग शून्य पूर्वाग्रह से सिद्ध होता है। अध्ययन निष्कर्ष निकालता है कि पृथक्करण विधि, पवन और डिफ्यूज़िविटी के लिए पावर-लॉ पैरामीट्राइजेशन के साथ मिलकर, जटिल संख्यात्मक सिमुलेशन की आवश्यकता के बिना अस्थिर वायुमंडलीय स्थितियों के तहत प्रदूषक प्रसार की भविष्यवाणी करने के लिए एक मजबूत विश्लेणात्मक उपकरण प्रदान करती है। यह कार्य अस्थिर सीमा परतों में प्रदूषक परिवहन की आवश्यक भौतिकी को पकड़ने के लिए ट्रांसफॉर्म-आधारित और सेपरेशन-आधारित विश्लेणात्मक समाधानों की उपयोगिता को सुदृढ़ करता है।