Finite element and box-method discretizations for fractional elliptic problems with quadrature and mass lumping

तकनीकी सारांश: भिन्नात्मक दीर्घिक समस्याओं के लिए परिमित तत्व और बॉक्स-विधि विविक्तीकरण (Finite Element and Box-Method Discretizations)

समस्या विवरण
यह शोध पत्र एक सीमित उत्तल डोमेन $S \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$) पर भिन्नात्मक दीर्घिक समस्या $L^\beta u = f$ के संख्यात्मक सन्निकटन को संबोधित करता है, जो सजातीय डिरिचलेट (Dirichlet) या न्यूमैन (Neumann) सीमा शर्तों के अधीन है। यहाँ, $L = -\text{div}(A\nabla) + \kappa^2$ एक द्वितीय-क्रम का स्व-अdej (self-adjoint) दीर्घिक ऑपरेटर है, और $\beta > 0$ है। यद्यपि भिन्नात्मक घात वाले दीर्घिक ऑपरेटर गैर-स्थानीय प्रसार (nonlocal diffusion), परवलयिक नियमितता सिद्धांत (parabolic regularity theory), और गॉसियन क्षेत्र मॉडलिंग (जैसे व्हिटल-मैटर्न वर्ग) में मानक हैं, केंद्रीय संख्यात्मक चुनौती $L^{-\beta}$ के व्युत्क्रम भिन्नात्मक घात का सटीक और कुशल सन्निकटन करना है।

एक महत्वपूर्ण मुद्दा जो पहचाना गया है वह यह है कि व्यावहारिक कार्यान्वयन अक्सर आदर्श गैलरकिन (Galerkin) सूत्रीकरण से विचलित होते हैं। गणना में, स्टिफनेस मैट्रिसेस (stiffness matrices), मास मैट्रिसेस (mass matrices), और लोड वेक्टर्स (load vectors) को आमतौर पर संख्यात्मक क्वाड्रचर (numerical quadrature) का उपयोग करके संयोजित किया जाता है। विशेष रूप से, $L^2$ आंतरिक गुणनफल (inner product) को अक्सर क्वाड्रचर-आधारित सन्निकटन द्वारा प्रतिस्थापित किया जाता है, जैसे कि मास लम्पिंग (mass lumping - विकर्ण मास मैट्रिसेस)। गैर-भिन्नात्मक समस्याओं के लिए, इन्हें "वेरिएशनल अपराध" (variational crimes) के रूप में माना जाता है (स्ट्रैंग और फिक्स)। हालांकि, भिन्नात्मक घातों के लिए, जहाँ विविक्त ऑपरेटर (discrete operator) एक परिमित-आयामी स्थान पर स्पेक्ट्रल कैलकुलस के माध्यम से परिभाषित होता है, वहाँ आंतरिक गुणनफल का चयन और द्विपद रूप (bilinear form) संरचनात्मक रूप से ऑपरेटर को ही परिभाषित करते हैं। फलस्वरूप, यह स्पष्ट नहीं है कि कार्यान्वयन के इन विवरणों को बदलने पर गैर-भिצאत्मक त्रुटि सिद्धांत भिन्नात्मक सेटिंग में कैसे स्थानांतरित होते हैं।

कार्यप्रणाली (Methodology)
लेखक एक एकीकृत अनुरूप खंडीय रैखिक (piecewise linear) ढांचे को विकसित करते हैं जो परिमित तत्व विधि (FEM) और बॉक्स-विधि (एक शीर्ष-केंद्रित फाइनाइट वॉल्यूम विधि) विविक्तीकरणों को एक एकल ऑपरेटर-सैद्धांतिक संरचना के भीतर उपचारित करता है। इस पद्धति का मूल भाग तीन विशिष्ट घटकों को अलग करना है:

1. दीर्घिक ऑपरेटर का सन्निकटन करने वाला विविक्त द्विपद रूप $a_h$।
2. लोड पद का विविक्तीकरण (या तो सटीक $L^2$ युग्मन या बॉक्स-विधि नियंत्रण-आयतन औसत)।
3. परीक्षण स्थान $V_h$ पर विविक्त ऑपरेटर को साकार करने के लिए उपयोग किया जाने वाला आंतरिक गुणनफल $\langle \cdot, \cdot \rangle_h$।

यह ढांचा $V_h$ पर सक्षम आंतरिक गुणनफल (admissible inner products) की एक श्रेणी पेश करता है। ये ऐसे आंतरिक गुणनफल हैं जो $L^2$ नॉर्म के समान (uniformly equivalent) हैं और सटीक $L^2$ आंतरिक गुणनफल के सापेक्ष $O(h^2)$ के क्रम में एक निरंतरता अनुमान (consistency estimate) को संतुष्ट करते हैं। इस श्रेणी में शामिल हैं:

• सटीक $L^2$ आंतरिक गुणनफल।
• मास-लम्पड आंतरिक गुणनफल (विकर्ण मास मैट्रिक्स)।
• इन दोनों के बीच अंतरपोलेट करने वाले स्पष्ट बैंडेड (banded) परिवार।

विविक्त भिन्नात्मक ऑपरेटर $L_h^{-\beta}$ को हिल्बर्ट स्पेस $(V_h, \langle \cdot, \cdot \rangle_h)$ पर स्पेक्ट्रल कैलकुलस के माध्यम से परिभाषित किया जाता है। लेख दो विशिष्ट योजनाओं का विश्लेषण करता है:

• भिन्नात्मक FEM: लोड के लिए सटीक $L^2$ प्रोजेक्टर का उपयोग करता है।
• भिन्नात्मक बॉक्स-विधि: लोड के लिए बॉक्स लोड ऑपरेटर $P_h$ (जिसमें नियंत्रण-आयतन औसत $Q$ शामिल है) का उपयोग करता है।

एक प्रमुख सैद्धांतिक उपकरण डुनफोर्ड-टेलर (Dunford–Taylor) फंक्शनल कैलकुलस का उपयोग करना है ताकि कंटूर इंटीग्रल्स के माध्यम से भिन्नात्मक घातों का प्रतिनिधित्व किया जा सके। यह लेखकों को द्विपद रूप के व्यवधानों, आंतरिक गुणनफल के व्यवधानों और लोड पद के व्यवधानों में त्रुटि अनुमानों को विभाजित करने की अनुमति देता है।

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

1. एकीकृत ऑपरेटर ढांचा: यह शोध पत्र एक ऐसे ढांचे को तैयार करता है जहाँ द्विपद रूप, लोड विविक्तीकरण और आंतरिक गुणनफल को अलग रखा गया है। यह स्पष्ट करता है कि $\beta \neq 1$ के लिए, आंतरिक गुणनफल का चयन केवल एक प्रतिनिधित्व संबंधी विकल्प नहीं है, बल्कि मौलिक रूप से विविक्त भिन्नात्मक ऑपरेटर को बदल देता है।
2. सक्षम आंतरिक गुणनफल: लेखक सक्षम आंतरिक गुणनफलों की एक विस्तृत श्रेणी को परिभाषित करते हैं जो व्यावहारिक स्पार्स (sparse) कार्यान्वयन (जैसे मास लम्पिंग) और सैद्धांतिक सटीक योजनाओं को समाहित करती है। वे सिद्ध करते हैं कि इस श्रेणी में स्पार्स मैस मैट्रिसेस भी शामिल हैं जो SPDE मॉडल में कुशल सांख्यिकीय अनुमान के लिए प्रासंगिक हैं।
3. मास लम्पिंग का अंतर्निहित लक्षण वर्णन: एक महत्वपूर्ण संरचनात्मक परिणाम 'लम्प्ड-मास' आंतरिक गुणनफल की पहचान है, जो सक्षम आंतरिक गुणनफलों में एकमात्र ऐसा विकल्प है जो भिन्नात्मक बॉक्स विविक्तीकरण को "अंतर्निहित" (intrinsic) बनाता है। विशेष रूप से, जब लम्प्ड-मास आंतरिक गुणनफल का उपयोग किया जाता है, तो भिन्नात्मक बॉक्स समाधान ऑपरेटर को गैर-भिन्नात्मक बॉक्स समाधान ऑपरेटर $T_h^{\text{box}}$ से सीधे संबंध $T_{h,\beta}^Q = Q^{-1}(Q T_h^{\text{box}})^\beta$ के माध्यम से पुनर्गठित किया जा सकता है। इसका अर्थ है कि भिन्नात्मक बॉक्स विधि केवल इस विशिष्ट आंतरिक गुणनफल के तहत ही गैर-भिन्नात्मक बॉक्स विधि के साथ सुसंगत है।
4. त्रुटि अनुमान (Error Estimates): लेख प्राकृतिक द्वितीय-क्रम निरंतरता धारणाओं के तहत दोनों FEM और बॉक्स-विधि विविक्तीकरणों के लिए कठोर त्रुटि अनुमान स्थापित करता है।
   • FEM के लिए, अभिसरण दर शास्त्रीय क्रम $O(h^{\min\{2\beta-\delta, 2\}})$ (महत्वपूर्ण मामलों में लॉगरिदमिक कारकों के साथ) को बनाए रखती है, चाहे सक्षम आंतरिक गुणनफल का उपयोग किया जाए।
   • बॉक्स-विधि के लिए, अभिसरण दर लोड विविक्तीकरण द्वारा सीमित होती है। त्रुटि $O(h^{\min\{2\beta+\sigma-\delta, 1+\sigma\}})$ के रूप में व्यवहार करती है, जहाँ $\sigma$ डेटा की नियमितता को दर्शाता है। यह पुष्टि करता है कि बॉक्स विधि नियंत्रण-आयतन औसत के कारण अपेक्षित सटीकता की हानि का सामना करती है, लेकिन यह हानि सक्षम आंतरिक गुणनफल के चयन से स्वतंत्र है।
5. कार्यान्वयन क्षमता: इस विश्लेषण को एक्सपोनेंशियल रूप से अभिसरण करने वाले रैशनल सन्निकटन (जैसे sinc quadrature) के साथ स्थानिक विविक्तीकरण को जोड़कर पूरी तरह से लागू करने योग्य विधियों तक विस्तारित किया गया है, जिससे यह सुनिश्चित होता है कि कुल त्रुटि स्थानिक विविक्तीकरण त्रुटि द्वारा नियंत्रित हो।

परिणाम
एक और दो आयामों में संख्यात्मक प्रयोग सैद्धांतिक निष्कर्षों का समर्थन करते हैं:

• एकसमान मेष (uniform meshes) और खंडीय स्थिरांक लोड के साथ 1D में, FEM और बॉक्स-विधि की त्रुटियां मेल खाती हैं, और सैद्धांतिक दरों को प्राप्त करती हैं।
• विलक्षण (singular) लोड के मामलों में, बॉक्स-विधि अभिसरण क्रम में अनुमानित संतृप्ति (saturation) प्रदर्शित करती है (लोड सन्निकटन द्वारा सीमित), जबकि FEM उच्च सटीकता बनाए रखती है।
• 2D में, आंतरिक गुणनफल (सटीक $L^2$, मास-लम्पड, और बैंडेड इंटरपोलेशन) को बदलने पर समान अभिसरण ढलान (convergence slopes) प्राप्त होते हैं, जो यह पुष्टि करता है कि सक्षम आंतरिक गुणनफल का चयन विविक्त ऑपरेटर को प्रभावित करता है लेकिन इसके आनुपातिक अभिसरण दर को नहीं।

महत्व
लेखक मुख्य रूप से दो क्षेत्रों में महत्व का दावा करते हैं:

1. सांख्यिकीय मॉडलिंग: यह कार्य सांख्यिकी में "SPDE दृष्टिकोण" के लिए एक विश्लेषणात्मक आधार प्रदान करता है, जहाँ गॉसियन रैंडम फील्ड्स को भिन्नात्मक SPDEs के समाधान के रूप में मॉडल किया जाता है। यह मास लम्पिंग और स्पार्स प्रिसिजन ऑपरेटरों के उपयोग को न्यायोचित ठहराता है, यह दिखाते हुए कि भले ही आंतरिक गुणनफल विविक्त ऑपरेटर को बदल देता है, फिर भी यह सन्निकटन के अभिसरण क्रम को कम नहीं करता है।
2. संख्यात्मक विश्लेषण: यह इस अस्पष्टता को सुलझाता है कि "वेरिएशनल अपराध" (क्वाड्रचर और मास लम्पिंग) भिन्नात्मक ऑपरेटरों को कैसे प्रभावित करते हैं। विविक्तीकरण के घटकों को अलग करके, लेखक प्रदर्शित करते हैं कि भिẪत्मक बॉक्स विधि की संरचनात्मक विशेषताएं केवल विशिष्ट आंतरिक गुणनफल विकल्पों के तहत संरक्षित रहती हैं, जो भिẪत्मक दीर्घिक समस्याओं के लिए कुशल और सुसंगत संख्यात्मक योजनाओं के निर्माण के लिए एक कठोर मार्गदर्शिका प्रदान करती हैं।

लेखक निष्कर्ष निकालते हैं कि विकसित ढांचा भविष्य के कार्यों के लिए दिशाएँ खोलता है, जैसे कि भिẪत्मक परवलयिक समस्याओं (जैसे भिẪत्मक ऊष्मा समीकरण) और उनसे जुड़े स्टोकेस्टिक पार्शियल डिफरेंशियल इक्वेशंस (SPDEs) का विश्लेषण।