A fully diagonalized spectral method on the unit ball

तकनीकी सारांश: इकाई गोले (Unit Ball) पर एक पूर्णतः विकेंद्रीकृत स्पेक्ट्रल विधि (A Fully Diagonalized Spectral Method)

समस्या का विवरण
यह शोध पत्र इकाई गोले $B_d \subset \mathbb{R}^d$ पर एक स्थिर श्रोडिंगर-प्रकार (Schrödinger-type) समीकरण के लिए डिरिचलेट बाउंडरी वैल्यू प्रॉब्लम (BVP) के संख्यात्मक समाधान को संबोधित करता है। विशिष्ट समीकरण है:
$$-\Delta u + \lambda(1 - \|x\|^2)^\kappa u = f(x), \quad x \in B_d$$
जो सीमा शर्त $u(x) = 0$ (गोले $\partial B_d = S^{d-1}$ पर) को पूरा करता है, जहाँ $\lambda > 0$ और $\kappa$ एक गैर-ऋणात्मक पूर्णांक है। मानक रूपांतर вариेशनल (variational) फॉर्मूलेशन $H^1_0(B_d)$ में $u$ खोजने के लिए है ताकि एक विशिष्ट द्विपद रूप (bilinear form) एक रैखिक फलन के बराबर हो।

कार्यप्रणाली (Methodology)
लेखक एक स्पेक्ट्रल विधि प्रस्तावित करता है जो गोले पर आंशिक अंतर समीकरणों (PDEs) के विवेकीकरण (discretization) से जुड़े बड़े, सघन रैखिक प्रणालियों (dense linear systems) के समाधान से बचती है। कार्यप्रणाली इस प्रकार है:

1. वेरिएशनल फ्रेमवर्क और सोबोलेव आंतरिक उत्पाद (Sobolev Inner Product):
   बध्य रूप (bilinear form) से जुड़े रूपांतर вариेशनल समस्या को $H^1_0(B_d)$ में एक आंतरिक उत्पाद के रूप में पहचाना जाता है। सीमा शर्त $u=0$ को लागू करने के लिए, परीक्षण फलनों (test functions) को $R^d_n = \{ (1-\|x\|^2)p(x) : p(x) \in \Pi^d_n \}$ से चुना गया है। एक विशिष्ट सोबोलेव आंतरिक उत्पाद $\langle \cdot, \cdot \rangle_S$ को बहुपद स्थान $\Pi^d$ पर इस प्रकार परिभाषित किया गया है:
   $$\langle u, v \rangle_S = \lambda \int_{B_d} (1-\|x\|^2)^{\kappa+2} u v \, dx + \int_{B_d} \nabla((1-\|x\|^2)u) \cdot \nabla((1-\|x\|^2)v) \, dx$$
   यह आंतरिक उत्पाद रूपांतर вариेशनल फॉर्मूलेशन से स्वाभाविक रूप से उत्पन्न होता है।

2. सोबोलेव ऑर्थोगोनल बहुपदों का निर्माण:
   विधि का मूल आधार $R^d_n$ के लिए एक आधार (basis) का निर्माण करना है जो $\langle \cdot, \cdot \rangle_S$ के सापेक्ष परस्पर ऑर्थोगोनल (mutually orthogonal) बहुपदों से बना हो।

• संरचना: आधार का निर्माण गोलाकार हार्मोनिक्स (spherical harmonics) $Y^{n-2j}_\nu$ और रेडियल चर $r = \|x\|$ पर निर्भर एकवर्थी बहुपदों (univariate polynomials) $q_j$ का उपयोग करके किया जाता है। बहुपद का रूप $R^n_{j,\nu}(x) = q_j(2\|x\|^2 - 1)Y^{n-2j}_\nu(x)$ है।
• एकवर्थी मामले में न्यूनीकरण (Reduction to Univariate Case): यह सिद्ध किया गया है कि बहुचर (multivariate) सोबोलेव ऑर्थोगोनैलिटी स्थिति रेडियल बहुपदों $q_j$ के लिए एक एकवर्थी सोबोलेव आंतरिक उत्पाद में बदल जाती है। इस एकवर्थी आंतरिक उत्पाद में एक भारित $L^2$ पद और $(1-t)q(t)$ के अवकलज (derivative) से संबंधित एक पद शामिल है।

3. कनेक्शन सूत्र और पुनरावर्ती गुणांक (Recursive Coefficients):
   यह शोध पत्र इन नए सोबोलेव ऑर्थोगोनल बहुपदों और शास्त्रीय जैकोबी बहुपदों (विशेष रूप से $P^{(1, \beta)}_k$) के बीच एक संबंध स्थापित करता है।

• संरचना: $q^{(\beta)}_k$ को शास्त्रीय जैकोबी बहुपदों के एक रैखिक संयोजन के रूप में व्यक्त किया जाता है।
• महत्वपूर्ण रूप से, यह संबंध "बैंडेड" (banded) है: $P^{(1, \beta)}_k$ को $q^{(\beta)}_k$ और केवल उससे पिछले $\kappa+1$ बहुपदों से जोड़ा गया है।
• एक एल्गोरिदम प्रदान किया गया है जो इनके कनेक्शन गुणांकों को पुनरावर्ती रूप से (recursively) गणना करने के लिए उपयोग किया जाता है। $\kappa=0$ के विशिष्ट मामले के लिए, गुणांकों का एक बंद-रूप (closed-form) प्रतिनिधित्व और एक गैर-रैखिक पुनरावृत्ति संबंध व्युत्पन्न किया गया है, जिससे उन्हें निरंतर भिन्नों (continued fractions) के माध्यम से व्यक्त किया जा सकता है।

4. पूर्णतः विकेंद्रीकृत समाधान (Fully Diagonalized Solution):
   समाधान $u(x)$ को सीमा कारक (boundary factor) $(1-\|x\|^2)$ के साथ सोबोलेव ऑर्थोगोनल बहुपदों के आधार में विस्तारित किया जाता है। चूंकि आधार ऑपरेटर के प्राकृतिक आंतरिक उत्पाद के सापेक्ष ऑर्थोगोनल है, इसलिए गैलेरकिन प्रोजेक्शन एक विकेंद्रित (diagonal) प्रणाली की ओर ले जाता है। समाधान के फूरियर-सोबोलेव गुणांकों को शास्त्रीय ऑर्थोगोनल आधार में स्रोत पद $f(x)$ के फूरियर गुणांकों से पुनरावर्ती रूप से प्राप्त किया जाता है, जिससे एक मैट्रिक्स को उलटने (invert) की आवश्यकता समाप्त हो जाती है।

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

• विकेंद्रीकरण (Diagonalization): प्राथमिक योगदान एक ऐसी स्पेक्ट्रल विधि है जो गोला पर डिरिचलेट समस्या के लिए एक पूर्णतः विकेंद्रित प्रणाली में परिणत होती है। यह मानक स्पेक्ट्रल विधियों के विपरीत है जो अक्सर पूर्ण, खराब-सशर्त (ill-conditioned) मैट्रिसेस प्रदान करती हैं।
• सोबोलेव ऑर्थोगोनल बहुपद: यह शोध पत्र एक बहुचर सोबोलेव ऑर्थोगोनल बहुपदों के रचनात्मक सिद्धांत को प्रदान करता है जो ग्रेडिएंट्स और एक विशिष्ट भार $(1-\|x\|^2)^\kappa$ वाले आंतरिक उत्पाद के संदर्भ में इकाई गोले पर निर्मित हैं।
• पुनरावर्ती एल्गोरिदम: यह इन सोबोलेव बहुपदों और शास्त्रीय जैकोबी बहुपदों के बीच स्पष्ट कनेक्शन सूत्र स्थापित करता है, जो शास्त्रीय ऑर्थोगोनल आधार में विस्तार गुणांकों की कुशल, पुनरावर्ती गणना को सक्षम बनाता है।
• संख्यात्मक सत्यापन: $d=2, \kappa=0$ के लिए एक संख्यात्मक प्रयोग इस विधि की सटीकता को प्रदर्शित करता है। त्रुटि विश्लेषण (error analysis) दिखाता है कि जैसे-जैसे सन्निकटन की डिग्री बढ़ती है, समाधान की एनालिटिसिटी को देखते हुए घातीय (सुपर-एक्सपोनेंशियल) अभिसरण (convergence) प्राप्त होता है।

महत्व और दावे
लेखक का दावा है कि यह दृष्टिकोण समस्या की रूपांतर вариेशनल संरचना का लाभ उठाकर संख्यात्मक कार्यान्वयन को महत्वपूर्ण रूप से सरल बनाता है। सोबोलेव ऑर्थोगोनल बहुपदों का उपयोग करके, यह विधि बड़े रैखिक सिस्टम को हल करने की कम्प्यूटेशनल बाधा को दूर करती है, जो गोले पर PDE के लिए सामान्य स्पेक्ट्रल विधियों की एक सीमा है। शोध पत्र का तर्क है कि व्युत्पन्न कनेक्शन सूत्र समाधान के गुणांकों के कुशल, पुनरावर्ती व्युत्पन्न को सक्षम करते हैं, जिससे यह विधि "पूर्णतः विकेंद्रित" बन जाती है। संख्यात्मक परिणाम पुष्टि करते हैं कि विधि विश्वसनीय और सटीक है, जो सुचारू समाधानों के लिए अपेक्षाकृत कम-डिग्री के सन्निकटन के साथ उच्च परिशुद्धता प्राप्त करती है।