Geometric Realism Without Angular Resolution Structural Classification of Multilayer Kubelka-Munk Theory within Radiative Transport

यहाँ क्लाउड ज़ेलर के शोध पत्र "Geometric Realism Without Angular Resolution: Structural Classification of Multilayer Kubelka–Munk Theory within Radiative Transport" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

कुबेल्का-मुंक (KM) सिद्धांत परतों वाले, प्रकीर्णन (scattering) और अवशोषण (absorption) करने वाले माध्यमों (जैसे पेंट, कागज, वस्त्र) में प्रकाश परिवहन को मॉडल करने के लिए उद्योग मानक है। अपनी अनुभवजन्य सफलता के बावजूद, इसे ऐतिहासिक रूप से एक "प्रपंचशास्त्रीय" (phenomenological) मॉडल के रूप में देखा जाता है जिसका पूर्ण रेडिएटिव ट्रांसफर इक्वेशन (RTE) के साथ अस्पष्ट सैद्धांतिक संबंध है।

• मुख्य मुद्दा: KM सिद्धांत पूर्ण कोणीय वितरण (angular distribution) को दो स्केलर फ्लक्स (फॉरवर्ड और बैकवर्ड) में संकुचित कर देता है। यह आइसोट्रोपिक प्रकीर्णन और अत्यधिक एनिसोट्रोपिक (फॉरवर्ड-पीक्ड) प्रकीर्णन के बीच अंतर नहीं कर सकता है।
• अंतराल (The Gap): हालांकि पिछले कार्यों ने हेमिस्फेरिक एवरेजिंग के माध्यम से RTE से KM गुणांकों को व्युत्पन्न किया था, लेकिन वे इस न्यूनीकरण (reduction) की गणितीय संरचना को पहचानने में विफल रहे। विशेष रूप से, उन्होंने यह स्पष्ट नहीं किया कि बहु-परत स्टैकिंग (composition) अभ्यास में इतनी अच्छी तरह से क्यों काम करती है, जबकि कोणीय जानकारी का नुकसान होता है, और न ही मॉडल की सटीकता की सीमाओं को कठोरता से परिभाषित किया।

2. कार्यप्रणाली (Methodology)

लेखक KM सिद्धांत को एक अनौपचारिक सन्निकटन (ad hoc approximation) के रूप में नहीं, बल्कि RTE के एक कठोर रैंक-2 ऑर्थोगोनल गैलरकिन प्रोजेक्शन (rank-2 orthogonal Galerkin projection) के रूप में पुनर्गठित करते हैं।

• गणितीय ढांचा:
  • RTE को वर्ग-समावेशी फलनों (square-integrable functions) $L^2([-1, 1])$ के स्थान पर कार्य करने वाले एक अनंत-आयामी ऑपरेटर के रूप में माना जाता है।
  • एक प्रोजेक्शन ऑपरेटर $P$ को हेमिस्फेरिक इंडिकेटर फलनों $\chi_+$ (फॉरवर्ड हेमिस्फीयर, $\mu > 0$) और $\chi_-$ (बैकवर्ड हेमिस्फीयर, $\mu < 0$) द्वारा स्पैन किए गए दो-आयामी उप-स्थान (subspace) पर परिभाषित किया गया है।
  • KM ट्रांसपोर्ट ऑपरेटर को ठीक उसी तरह व्युत्पन्न किया गया है जैसे कि "सैंडविच" ऑपरेशन: $T_{KM} = P T P$, जहाँ $T$ पूर्ण RTE ऑपरेटर है।
• व्युत्पत्ति के चरण:
  1. स्ट्रीमिंग टर्म (Streaming Term): स्ट्रीमिंग टर्म को प्रोजेक्ट करने से ट्रांसपोर्ट गुणांकों के बीच क्लासिक 'कारक-2' संबंध प्राप्त होता है।
  2. प्रकीर्णन टर्म (Scattering Term): प्रकीर्णन इंटीग्रल को हेमिस्फीयरों पर औसत निकाला जाता है। फेज फंक्शन $p(\mu, \mu')$ को उसके हेमिस्फेरिक मोमेंट्स, विशेष रूप से बैकस्कैटर अंश $\bar{p}_{-+}$ में कम कर दिया जाता है।
  3. परिणामी समीकरण: यह व्युत्पत्ति बिना किसी अतिरिक्त भौतिक सन्निकटन के मानक KM डिफरेंशियल समीकरणों को $K = 2\sigma_a$ और $S = \sigma_s \bar{p}_{-+}$ के साथ प्रदान करती है।
• त्रुटि विश्लेषण (Error Analysis): पेपर एक प्रोजेक्शन एरर मेट्रिक $\varepsilon = \|I^{(\perp)}\| / \|I\|$ पेश करता है, जहाँ $I^{(\perp)}$ तीव्रता का वह घटक है जो $P$ के अनंत-आयामी कर्नेल (kernel) में स्थित है (अर्थात वे कोणीय विवरण जिन्हें प्रोजेक्शन द्वारा त्याग दिया गया है)।

3. मुख्य योगदान (Key Contributions)

A. संरचनात्मक वर्गीकरण (Structural Classification)

यह पेपर KM सिद्धांत को परिवहन सन्निकटन के पदानुक्रम में एक विशिष्ट बिंदु के रूप में स्थापित करता है:

• KM सिद्धांत: $\text{span}\{\chi_+, \chi_-\}$ पर प्रोजेक्ट करता है। यह ज्यामितीय दिशात्मकता (फॉरवर्ड बनाम बैकवर्ड) को सुरक्षित रखता है लेकिन हेमिस्फीयर के भीतर सभी आंतरिक कोणीय रिज़ॉल्यूशन को त्याग देता है।
• डिफ्यूजन ($P_1$): $\text{span}\{1, \mu\}$ पर प्रोजेक्ट करता है। यह पहले कोणीय मोमेंट (ग्रेडिएंट) को सुरक्षित रखता है लेकिन सख्त फॉरवर्ड/बैकवर्ड अंतर को खो देता है।
• डिस्क्रीट ऑर्डिनेट्स ($S_N$): विशिष्ट कोणों को सैंपल करने वाली एक कोलोकेशन विधि, जो KM में उपयोग किए गए गैलरकिन एवरेजिंग से भिन्न है।

B. रैंक-संरक्षण प्रमेय (The Rank-Preservation Theorem)

एक महत्वपूर्ण सैद्धांतिक निष्कर्ष यह है कि प्रोजेक्शन का रैंक बहु-परत संरचना (multilayer composition) के तहत संरक्षित रहता है।

• चूंकि प्रोजेक्शन $P$ इडेम्पोटेंट ($P^2=P$) है और इसका रैंक 2 है, इसलिए $N$ परतों को स्टैक करना (ट्रांसफर मैट्रिसेस को गुणा करना) एक संयुक्त ऑपरेटर की ओर ले जाता है जो अभी भी 2D स्पेस को 2D स्पेस में मैप करता है।
• निहितार्थ: परतों की स्टैकिंग की मात्रा प्रारंभिक प्रोजेक्शन द्वारा त्याग दी गई कोणीय जानकारी को वापस नहीं ला सकती। 100-परत वाला KM कैलकुलेशन में भी वही कोणीय रिज़ॉल्यूशन (शून्य) होता है जो एक एकल-परत वाले कैलकुलेशन में होता है।

C. प्रोजेक्शन-एरर मानदंड (The Projection-Error Criterion)

पेपर उन स्थितियों को परिभाषित करता है जिनमें KM सिद्धांत सटीक है:

• सटीकता की स्थिति: KM तब अच्छा काम करता है जब वास्तविक समाधान $I$ का अधिकांश ऊर्जा $P$ के रेंज में होती है (अर्थात, हेमिस्फीयर के भीतर कोणीय वितरण लगभग सपाट होता है)।
• विफलता की स्थिति: KM तब विफल होता है जब महत्वपूर्ण ऊर्जा $\ker(P)$ में रहती है, जैसे कि स्ट्रॉन्गली फॉरवर्ड-पीक्ड प्रकीर्णन या ऑप्टिकली थिन (बैलिस्टिक) शासन में।
• शासी पैरामीटर (Governing Parameter): सटीकता रिड्यूस्ड ऑप्टिकल थिकनेस $\tau^* = \tau(1-g)$ द्वारा शासित होती है, जहाँ $g$ विषमता पैरामीटर (asymmetry parameter) है। जब $\tau^* \gg 1$ होता है, तो मल्टीपल स्कैटरिंग प्रकाश को आइसोट्रोपाइज कर देती है, जिससे समाधान $P$ के रेंज में आ जाता है और त्रुटि न्यूनतम हो जाती है।

4. परिणाम और संख्यात्मक चित्रण (Results and Numerical Illustrations)

लेखक एक उच्च-क्रम $S_{32}$ डिस्क्रीट-ऑर्डिनेट्स सॉल्वर के विरुद्ध संख्यात्मक तुलना का उपयोग करके सिद्धांत को मान्य करते हैं:

• आइसोट्रोपिक प्रकीर्णन ($g=0$): प्रोजेक्शन एरर $\varepsilon \approx 0.20$; रिफ्लेक्टेंस एरर नगण्य है ($\sim 0.001$)।
• मध्यम एनिसोट्रॉपी ($g=0.50$): प्रोजेक्शन एरर $\varepsilon \approx 0.25$।
• स्ट्रॉन्ग फॉरवर्ड पीकिंग ($g=0.85$): प्रोजेक्शन एरर $\varepsilon \approx 0.34$।
• अवलोकन: जबकि आंतरिक प्रोजेक्शन एरर $\varepsilon$, $g$ के साथ बढ़ता है, मल्टीपल स्कैटरिंग मीडिया में दृश्यमान रिफ्लेक्टेंस ($R$) में त्रुटि छोटी बनी रहती है। यह समझाता है कि क्यों कंपोजिशनल मॉडल (जैसे हेबर्ट और हर्ष द्वारा) मुद्रित कागज में उच्च रंग सटीकता प्राप्त करते हैं: सबस्ट्रेट (कागज, स्याही) की भौतिकी ने KM प्रोजेक्शन कार्य करने से पहले ही प्रकाश को आइसोट्रोपाइज कर दिया है, जिससे वह जानकारी प्रभावी रूप से नष्ट हो गई है जिसे KM नहीं देख सकता।

5. महत्व और निहितार्थ (Significance and Implications)

• सैद्धांतिक वैधता: यह पेपर KM सिद्धांत को "प्रपंचशास्त्र" (phenomenology) से "कठोर लो-रैंक सन्निकटन" (rigorous low-rank approximation) की ओर ले जाता है। यह सिद्ध करता है कि KM एक सटीक प्रोजेक्शन है, न कि एक त्रुटिपूर्ण सन्निकटन, बशर्ते उपयोगकर्ता इसके डोमेन की वैधता को समझता हो।
• सफलता की व्याख्या: यह उस विरोधाभास को हल करता है कि जटिल बहु-परत प्रणालियों (जैसे प्रिंटिंग) में एक रैंक-2 मॉडल क्यों सफल होता है। सफलता इसलिए नहीं है कि मॉडल परिष्कृत है, बल्कि इसलिए है क्योंकि भौतिक माध्यम (कागज, स्याही) स्वाभाविक रूप से प्रकाश को आइसोट्रोपाइज कर देता है, जिससे त्याग दी गई कोणीय जानकारी अप्रासंगिक हो जाती है।
• सीमाओं को परिभाषित करना: यह मात्रात्मक मानदंड ($\tau^*$ और $\varepsilon$) प्रदान करता है कि कब KM को छोड़ देना चाहिए। उन शासनों में जहाँ $\tau^*$ छोटा है (पतली फिल्में) या $g$ बहुत अधिक है (जैविक ऊतक, वायुमंडलीय एरोसोल), प्रोजेक्शन एरर बड़ा होता है, और वहां उच्च-क्रम विधियों ($P_N$, $S_N$, मोंटे कार्लो) की संरचनात्मक आवश्यकता होती है।
• इनवर्स प्रॉब्लम: यह पेपर स्पष्ट करता है कि बहु-परत KM के लिए इनवर्स प्रॉब्लम अनिर्धारित (underdetermined) है क्योंकि अवलोकन योग्य (कुल $R$ और $T$) $P$ के 2D रेंज में रहते हैं। प्रति परत दो से अधिक पैरामीटर को रिकवर करने के लिए कोणीय-रिज़ॉल्व्ड माप या स्पेक्ट्रल मल्टीप्लेक्सिंग की आवश्यकता होती है, लेकिन फिर भी, अनंत-आयामी कर्नेल अदृश्य बना रहता है।

निष्कर्ष:
यह पेपर सफलतापूर्वक कुबेल्का-मुंक सिद्धांत को "कोणीय रिज़ॉल्यूशन के बिना ज्यामितीय यथार्थवाद" (geometric realism without angular resolution) के सन्निकटन के रूप में वर्गीकृत करता है। यह उन उपकरणों को प्रदान करता है जिनसे यह अनुमान लगाया जा सके कि यह सन्निकटन कब पर्याप्त है और कब विफल होता है, और इसकी अनुभवजन्य सफलता को मल्टीपल स्कैटरिंग के भौतिकी और ऑर्थोगोनल प्रोजेक्शन के गणितीय गुणों में स्थापित करता है।