तकनीकी सारांश: गैर-अवशेष (Non-Residual) अनुमानकों द्वारा संचालित सममित द्वितीय-क्रम रैखिक दीर्घायु (Elliptic) PDE के लिए अनुकूलित FEM की इष्टतम जटिलता

1. समस्या विवरण

यह शोध पत्र सममित द्वितीय-क्रम रैखिक दीर्घायु आंशिक अवकल समीकरणों (PDEs) के लिए अनुकूलित परिमित तत्व विधि (AFEM) को संबोधित करता है, जो इस रूप में हैं:
$$-\text{div}(A \nabla u^\star) + c u^\star = f - \text{div}(\mathbf{f}) \quad \text{in } \Omega, \quad u^\star = 0 \text{ on } \partial\Omega$$
जहाँ $\Omega \subset \mathbb{R}^d$ एक परिबद्ध बहुफलकीय लिप्सचिट्ज़ डोमेन (polyhedral Lipschitz domain) है। मुख्य चुनौती दो स्रोतों से उत्पन्न त्रुटियों को एक साथ नियंत्रित करने की है:

1. विविक्तकरण त्रुटि (Discretization Error): जो एक मेष (mesh) पर परिमित तत्व सन्निकटन से उत्पन्न होती है।
2. बीजगणितीय त्रुटि (Algebraic Error): जो पुनरावृत्ति सॉल्वर (iterative solvers) के माध्यम से परिणामी रैखिक प्रणालियों के अनिश्चित समाधान से उत्पन्न होती है।

पिछले कार्यों के विपरीत, जो अक्सर सटीक विविक्त समाधानों को मानते हैं या विशेष रूप से अवशेष-आधारित (residual-based) त्रुटि अनुमानकों पर निर्भर करते हैं, यह कार्य गैर-अवशेष त्रुटि अनुमानकों (जैसे कि औसत-आधारित ZZ-प्रकार के अनुमानक और संतुलित फ्लक्स (equilibrated flux) अनुमानक) पर ध्यान केंद्रित करता है, जबकि एक अनिश्चित पुनरावृत्ति बीजगणितीय सॉल्वर को अनुकूलित लूप में एकीकृत करता है। लक्ष्य यह सिद्ध करना है कि अनुकूलन एल्गोरिदम अनवर्तनीय पूर्ण R-रैखिक अभिसरण (unconditional full R-linear convergence) और कुल कम्प्यूटेशनल लागत के संबंध में इष्टतम जटिलता (optimal complexity) प्राप्त करता है।

2. कार्यप्रणाली और ढांचा

2.1. अनुकूलन एल्गोरिदम

लेखक एक अनुकूलन एल्गोरिदम (एल्गोरिदम A) प्रस्तावित करते हैं जो चार मॉड्यूल के माध्यम से लूप करता है: SOLVE, ESTIMATE, MARK, और REFINE।

• SOLVE और ESTIMATE: ये मॉड्यूल आपस में जुड़े हुए हैं। एक पुनरावृत्ति बीजगणितीय सॉल्वर (संकुचन मानचित्र $\Psi_\ell$) तब तक लागू किया जाता है जब तक कि गणना योग्य बीजगणितीय त्रुटि अनुमानक $\zeta_\ell$ के आधार पर एक स्टॉपिंग मानदंड $\mu_\ell$ के सापेक्ष पूरा न हो जाए। विशेष रूप से, सॉल्वर तब रुकता है जब $\zeta_\ell(u^k_\ell) \leq \lambda \mu_\ell(u^k_\ell)$ हो जाता है।
• MARK: डॉर्लर मार्किंग मानदंड (Dörfler marking criterion) का उपयोग करके गैर-अवशेष अनुमानक $\mu_\ell$ के आधार पर तत्वों के एक सेट को चिह्नित किया जाता है।
• REFINE: चिह्नित तत्वों को न्यूएस्ट-वर्टेक्स बाइसेक्शन (Newest-Vertex Bisection - NVB) का उपयोग करके परिष्कृत (refine) किया जाता है।

2.2. प्रमुख धारणाएं

विश्लेषण $\mu_\ell$ के अमूर्त गुणों और बीजगणितीय सॉल्वर पर आधारित है:

• संकुचनकारी सॉल्वर (Contractive Solver): पुनरावृत्ति सॉल्वर एक संकुचन गुण $|||u^\star_\ell - \Psi_\ell(v_\ell)||| \leq q_{ctr} |||u^\star_\ell - v_\ell|||$ को संतुष्ट करता है जहाँ $0 < q_{ctr} < 1$, जो मेष आकार से स्वतंत्र है।
• गैर-अवशेष अनुमानक गुण: अनुमानक $\mu_\ell$ को निम्नलिखित का पालन करना चाहिए:
  1. स्थानीय तुल्यता (Local Equivalence): $\mu_\ell$ सटीक विविक्त समाधान $u^\star_\ell$ के लिए मानक अवशेष-आधारित अनुमानक $\eta_\ell$ के स्थानीय रूप से तुल्य है। विशेष रूप से, $\eta_\ell(U_\ell; u^\star_\ell) \lesssim \mu_\ell(T^m_\ell[U_\ell]; u^\star_\ell)$ और इसके विपरीत, जहाँ $T^m_\ell$ स्तर $m$ का एक पैच है।
  2. दुर्बल स्थिरता (Weak Stability): $\mu_\ell$ एक पैच स्तर $r$ वाले दो अलग-अलग विविक्त फलनों के बीच एक स्थिरता स्थिति को संतुष्ट करता है।
• मेष परिष्करण (Mesh Refinement): NVB का उपयोग आकार नियमितता (shape regularity), ओवरले अनुमान (overlay estimates), और मेष-समापन (mesh-closure) जैसे मानक गुणों को सुनिश्चित करता है।

2.3. क्वासी-त्रुटि (The Quasi-Error)

विश्लेषण की जाने वाली केंद्रीय मात्रा क्वासी-त्रुटि $M^k_\ell$ है, जिसे बीजगणितीय त्रुटि और विविक्तकरण त्रुटि अनुमानक के योग के रूप में परिभाषित किया गया है:
$$M^k_\ell := |||u^\star_\ell - u^k_\ell||| + \mu_\ell(u^\star_\ell)$$
ध्यान दें कि $u^\star_\ell$ (सटीक FE समाधान) को कभी भी नहीं निकाला जाता है; $\mu_\ell(u^\star_\ell)$ एक सैद्धांतिक निर्माण है जिसका उपयोग विश्लेषण के लिए किया जाता है, जबकि इसके गणना योग्य समकक्ष $\zeta_\ell(u^k_\ell) + \mu_\ell(u^k_\ell)$ को इसके तुल्य दिखाया गया है।

3. प्रमुख योगदान और परिणाम

3.1. अनवर्तनीय पूर्ण R-रैखिक अभिसरण

प्राथमिक सैद्धांतिक परिणाम (थ्योरम 8) स्थापित करता है कि क्वासी-त्रुटि $M^k_\ell$ अनवर्तनीय रूप से और R-रैखिक रूप से अभिसरित होती है। अर्थात, अनुकूलन इतिहास में किन्हीं भी दो सूचकांकों $(\ell, k)$ और $(\ell', k')$ के लिए, $C_{lin} > 0$ और $0 < q_{lin} < 1$ ऐसे स्थिरांक मौजूद हैं कि:
$$M^k_\ell \leq C_{lin} q_{lin}^{|(\ell, k)| - |(\ell', k')|} M^{k'}_{\ell'}$$
महत्वपूर्ण रूप से, यह अभिसरण अनुकूलन मापदंडों ($\theta, C_{mark}, \lambda$) के किसी भी चयन के लिए बना रहता है। यह अभिसरण सुनिश्चित करने के लिए "पर्याप्त छोटे" मापदंडों की आवश्यकता को समाप्त करता है, जो अनिश्चित सॉल्वर से संबंधित पूर्व साहित्य में एक सामान्य प्रतिबंध है।

3.2. इष्टतम जटिलता (Optimal Complexity)

शोध पत्र यह सिद्ध करता है कि कुल कम्प्यूटेशनल लागत (संचयी डिग्री ऑफ फ्रीडम और सॉल्वर स्टेप्स द्वारा मापी गई) के संबंध में क्वासी-त्रुटि की क्षय दर (decay rate) इष्टतम है।

• थ्योरम 15: यदि अनुकूलन मापदंड $\theta$ और $\lambda$ पर्याप्त रूप से छोटे चुने जाते हैं, तो एल्गोरिदम इष्टतम अभिसरण दर प्राप्त करता है। विशेष रूप से, क्वासी-त्रुटि की क्षय दर नॉन-लीनियर एप्रोक्सिमेशन क्लास $\mathcal{A}^s$ में सर्वोत्तम संभव एप्रोक्सिमेशन दर से मेल खाती है।
• परिणाम यह दर्शाता है कि एल्गोरिदम अनावश्यक सॉल्वर पुनरावृत्तियों या मेष परिष्करण पर अनावश्यक कम्प्यूटेशनल संसाधनों को नष्ट नहीं करता है, बशर्ते कि मापदंडों को सही ढंग से ट्यून किया गया हो।

3.3. विशिष्ट अनुमानकों पर अनुप्रयोग

इस अमूर्त ढांचे को दो विशिष्ट गैर-अवशेष अनुमानकों के वर्गों पर लागू किया गया है, यह सिद्ध करते हुए कि वे आवश्यक स्थानीय तुल्यता और स्थिरता मान्यताओं को पूरा करते हैं:

1. ZZ-प्रकार के औसत अनुमानक (ZZ-Type Averaging Estimators): ज़िएनिएव्स्की और ज़ू (Zienkiewicz and Zhu) के मौलिक कार्य पर आधारित। यह शोध पत्र किसी भी बहुपद डिग्री $p \geq 1$ के लिए अवशेष अनुमानक के साथ स्थानीय तुल्यता को सिद्ध करता है (थ्योरम 16)।
2. संतुलित फ्लक्स अनुमानक (Equilibrated Flux Estimators): स्थानीय फ्लक्स पुनर्निर्माण (जैसे रैवियट-थॉमस तत्वों) पर आधारित। यह शोध पत्र इन अनुमानकों के लिए स्थानीय तुल्यता और दुर्बल स्थिरता स्थापित करता है (थ्योरम 22), और उनकी $p$-रोबस्टनेस (p-robustness) को नोट करता है।

3.4. संख्यात्मक प्रयोग

अनुभाग 6 एक L-आकार के डोमेन (एक सिंगुलैरिटी वाली समस्या) पर 2D संख्यात्मक प्रयोग प्रस्तुत करता है। प्रयोग तुलना करते हैं:

• मानक अवशेष-आधारित अनुमानक।
• ZZ-प्रकार के अनुमानक।
• संतुलित फ्लक्स अनुमानक।

परिणाम पुष्टि करते हैं:

• तीनों अनुमानक तुलनीय मेष उत्पन्न करते हैं जहाँ परिष्करण सिंगुलैरिटी पर केंद्रित है।
• दोनों गैर-अवशेष अनुमानक डिग्री ऑफ फ्रीडम और संचयी रनटाइम के संबंध में इष्टतम अभिसरण दर प्राप्त करते हैं।
• संतुलित फ्लक्स अनुमानक बेहतर दक्षता सूचकांक (1 के करीब) और $p$-रोबस्टनेस प्रदर्शित करता है, हालांकि इसे सख्त स्टॉपिंग मानदंडों के कारण अधिक सॉल्वर पुनरावृत्तियों की आवश्यकता हो सकती है।

4. महत्व और साहित्य से संबंध

लेखक अपने कार्य को मौजूदा साहित्य के एक एकीकरण और विस्तार के रूप में प्रस्तुत करते हैं:

• बनाम [KS11, CN12]: उन कार्यों के विपरीत जो सटीक सॉल्वर के साथ इष्टतम दरों पर ध्यान केंद्रित करते हैं, यह शोध पत्र अनिश्चित सॉल्वर को शामिल करता है और इष्टतम जटिलता (लागत बनाम त्रुटि) पर ध्यान केंद्रित करता है। इसके अलावा, यह कार्य [CN12] (पर्याप्त सूक्ष्म प्रारंभिक मेष, आंतरिक-नोड गुण) और [KS11] (निम्नतम-क्रम FEM, पड़ोसी परिष्करण) के प्रतिबंधात्मक अनुमानों से बचता है।
• बनाम [CKNS08, BM09, CFPP14]: ये कार्य आमतौर पर सटीक FE समाधानों की आवश्यकता रखते हैं या विचलन तर्कों (perturbation arguments) पर निर्भर करते हैं जो केवल छोटे मापदंडों के लिए अभिसरण की गारंटी देते हैं। यह शोध पत्र किसी भी पैरामीटर चयन के लिए अनवर्तनीय अभिसरण प्रदान करता है।
• बनाम [BFM+25]: जबकि [BFM+25] अनिश्चित सॉल्वर के लिए इष्टतम जटिलता का विस्तार करता है, यह अवशेष-आधारित अनुमानकों तक सीमित है। यह शोध पत्र इन परिणामों को गैर-अवशेष अनुमानकों (ZZ और संतुलित फ्लक्स) तक विस्तारित करने वाला पहला कार्य है, जो व्यवहार में व्यापक रूप से उपयोग किए जाते हैं लेकिन विश्लेषणात्मक रूप से अधिक चुनौतीपूर्ण हैं क्योंकि उनमें प्रत्यक्ष अवशेष संरचना का अभाव होता है।

मुख्य नवाचार: शोध पत्र इस कठिनाई को दूर करता है कि गैर-अवशेष और अवशेष अनुमानकों के बीच स्थानीय तुल्यता आमतौर पर केवल सटीक विविक्त समाधान के लिए होती है (जिसे कभी भी नहीं निकाला जाता है)। [BFM+25] में विश्लेषण में एक सूक्ष्म संशोधन का उपयोग करके और गैर-अवशेष अनुमानक की दुर्बल स्थिरता का लाभ उठाकर, लेखक गणना किए गए अनिश्चित समाधान और सैद्धांतिक सटीक विविक्त समाधान के बीच के अंतर को पाटते हैं, जिससे अनवर्तनीय अभिसरण और इष्टतम जटिलता सिद्ध होती है।

5. निष्कर्ष

यह कार्य अनिश्चित सॉल्रों के साथ अनुकूलित परिमित तत्व विधियों में गैर-अवशेष त्रुटि अनुमानकों का उपयोग करने के लिए एक कठोर गणितीय आधार प्रदान करता है। यह प्रदर्शित करता है कि, सामान्य धारणाओं के तहत, ये विधियाँ न केवल अनवर्तनीय रूप से अभिसरित होती हैं, बल्कि इष्टतम कम्प्यूटेशनल जटिलता भी प्राप्त करती हैं। यह उन लोकप्रिय अनुमानकों (जैसे ZZ-प्रकार और संतुलित फ्लक्स) के व्यावहारिक उपयोग को मान्य करता है जहाँ रैखिक प्रणालियों को सटीक रूप से हल करना कम्प्यूटेशनल रूप से अत्यधिक खर्चीला होता है।