Towards $C^0$ finite element methods for fourth-order elliptic equation. Part I: general boundary conditions

तकनीकी सारांश: चतुर्थ-क्रम के एलिप्टिक समीकरणों के लिए $C^0$ परिमित तत्व विधियों की ओर, भाग I: सामान्य सीमा स्थितियाँ

समस्या विवरण
यह शोध पत्र चतुर्थ-क्रम के एलिप्टिक समीकरणों, विशेष रूप से बाइहारमोनिक समीकरण $\Delta^2 u = f$ के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है, जो कि बहुभुजीय डोमेन $\Omega \subset \mathbb{R}^2$ पर आधारित है, जिसमें सामान्य सीमा स्थितियाँ (boundary conditions) दी गई हैं। सीमा $\partial \Omega$ को विलगित डिरिचलेट ($\Gamma_D$) और न्यूमैन ($\Gamma_N$) खंडों में विभाजित किया गया है। प्राथमिक चुनौती बाइहारमोनिक समीकरण के लिए प्रभावी $C^0$ परिमित तत्व विधियों (finite element methods) को विकसित करने में निहित है। जबकि नेवियर सीमा स्थितियों ($\Gamma_D = \partial \Omega$) के साथ बाइहारमोनिक समीकरण के लिए मानक दृष्टिकोण अक्सर पॉइसन समीकरणों के एक निकाय में समस्या को विघटित (decouple) करने पर निर्भर करते हैं, यह "नाइव" (naive) मिश्रित सूत्रीकरण गैर-उत्तल (non-convex) डोमेन या विशिष्ट मिश्रित सीमा स्थितियों के लिए वास्तविक समाधान की ओर अभिसरण (converge) करने में विफल रहता है। यह विफलता सापोग्यन विरोधाभास (Sapongyan paradox) के कारण होती है, जहाँ विघटित मिश्रित सूत्रीकरण का समाधान स्थान मूल बाइहारमोनिक समस्या के समाधान स्थान के साथ मेल नहीं खाता है, जिसका कारण पुन: प्रवेश करने वाले कोनों (re-entrant corners) से उत्पन्न होने वाले निम्न-नियमितता वाले पद (low-regularity terms) हैं।

कार्यप्रणाली
लेखक एक संशोधित मिश्रित सूत्रीकरण (modified mixed formulation) प्रस्तावित करते हैं जो बाइहारमोनिक समस्या को पॉइसन समीकरणों के एक निकाय में विघटित करता है, जिसे अतिरिक्त मध्यवर्ती पॉइसन समस्याओं द्वारा संवर्धित किया जाता है ताकि यह सुनिश्चित किया जा सके कि समाधान सही सोबोलेव स्थान (Sobolev space, $H^2(\Omega)$) के भीतर रहे।

1. समाधान स्थान का विश्लेषण: शोध पत्र लैपलेस ऑपरेटर $\Delta$ के विश्लेषण को करता है जो $H^2(\Omega)$ (उपयुक्त सीमा स्थितियों के साथ) को $L^2(\Omega)$ में मैप करता है। यह पहचानता है कि बड़े आंतरिक कोणों वाले बहुभुजीय डोमेन के लिए, रेंज $S = \Delta(V^2)$ $L^2(\Omega)$ का एक उचित बंद उपसमूह (closed subspace) $S^\perp$ है।
2. आधार निर्माण (Basis Construction): लेखक $S^\perp$ के लिए एक आधार का निर्माण करते हैं जो सबसे बड़े आंतरिक कोण $\omega$ के पास परिभाषित फलनों $\xi_m$ से बना है। ये फलन सिंगुलर पदों $s_m(r, \theta)$ (जो लैपलेस के कोणीय भाग के आइजन मानों से प्राप्त होते हैं) और नियमित सुधार पदों $\zeta_m$ (जो सहायक पॉइसन समस्याओं को हल करके प्राप्त किए जाते हैं) से मिलकर बने हैं। $S^\perp$ का आयाम ($d_\perp$) सबसे बड़े आंतरिक कोण के परिमाण और आसन्न किनारों पर विशिष्ट सीमा स्थितियों (डिरिचलेट/न्यूमैन) के संयोजन पर निर्भर करता है।
3. संशोधित एल्गोरिदम: प्रस्तावित $C^0$ परिमित तत्व एल्गोरिदम (एल्गोरिदम 3.1) इस प्रकार आगे बढ़ता है:
   • एक सहायक फलन $w$ के लिए प्राथमिक पॉइसन समस्या को हल करना।
   • ऑर्थोगोनल पूरक के लिए आधार फलनों $\xi_{m,h}$ की गणना करने हेतु सहायक पॉइसन समस्याओं को हल करना।
   • $w$ को $S$ के उपसमूह पर प्रोजेक्ट करना (ग्राम मैट्रिक्स सिस्टम के माध्यम से गुणांकों की गणना करके $S^\perp$ में इसके घटक को हटाना)।
   • प्राप्त प्रोजेक्टेड स्रोत पद $w_S$ का उपयोग करके अंतिम सन्निकटन $u_h$ प्राप्त करने के लिए दूसरी पॉइसन समस्या को हल करना।
4. न्यूमैन स्थितियों का विस्तार: कार्यप्रणाली को शुद्ध न्यूमैन मामले ($\Gamma_N = \partial \Omega$) के लिए विस्तारित किया गया है, जहाँ संगतता स्थिति (compatibility condition) $\int_\Omega f = 0$ आवश्यक है, और दिखाया गया है कि पुन: प्रवेश करने वाले कोनों के लिए ऑर्थोगोनल पूरक का आयाम 1 है।

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

• सापोग्यन विरोधाभास का समाधान: यह शोध पत्र एक कठोर सैद्धांतिक ढांचा प्रदान करता है जो यह दर्शाता है कि सामान्य बहुभुजीय डोमेन और मिश्रित सीमा स्थितियों के लिए बाइहारमोनिक समीकरण का नाइव विघटन विफल हो जाता है। यह एक सुधारात्मक सूत्रीकरण प्रस्तावित करता है जो वास्तविक समाधान की ओर अभिसरण सुनिश्चित करता है।
• सामान्यीकृत विघटन: नेवियर स्थितियों तक सीमित पिछले कार्यों के विपरीत, यह विधि डिरिचलेट और न्यूमैन स्थितियों के मनमाने संयोजनों को संभालती है। आवश्यक अतिरिक्त पॉइसन समीकरणों की संख्या (0, 1, या 2) सबसे बड़े आंतरिक कोण और स्थानीय सीमा स्थिति के प्रकारों द्वारा स्पष्ट रूप से निर्धारित होती है।
• त्रुटि विश्लेषण: लेखक प्रस्तावित $C^0$ रैखिक परिमित तत्व विधि के लिए कठोर त्रुटि अनुमान स्थापित करते हैं। वे सिद्ध करते हैं कि $H^1$ नॉर्म में समाधान $u_h$ की अभिसरण दर $O(h^{\min\{1, 2\alpha\}})$ है, जहाँ $\alpha$ आंतरिक कोण और सीमा स्थितियों पर निर्भर करता है। विशेष रूप से, $\omega > \pi$ वाले मिश्रित सीमा स्थितियों के लिए, दर सिंगुलैरिटी (singularity) द्वारा सीमित होती है, जबकि नेवियर/न्यूमैन मामलों के लिए, इष्टतम दरें अक्सर पुनः प्राप्त की जाती हैं।
• संख्यात्मक सत्यापन: विभिन्न बहुभुजीय डोमेन (L-आकार और पुन: प्रवेश कोने वाली ज्यामिति सहित) पर व्यापक संख्यात्मक प्रयोग प्रस्तावित सैद्धांतिक भविष्यवाणियों की पुष्टि करते हैं। परिणाम पुष्टि करते हैं कि जहाँ ऑर्थोगोनल पूरक गैर-तुच्छ (non-trivial) है, वहाँ नाइव विधि गलत समाधानों की ओर अभिसरण करती है, जबकि प्रस्तावित संशोधित विधि संदर्भ समाधान (जो $C^0$-IPDG द्वारा गणना किया गया है) की ओर अभिसरित होती है।

परिणाम और महत्व
शोध पत्र प्रदर्शित करता है कि प्रस्तावित संशोधित मिश्रित सूत्रीकरण सामान्य बहुभुजीय डोमेन पर बाइहारमोनिक समीकरण को हल करने के लिए सुव्यवस्थित (well-posed) और प्रभावी है। संख्यात्मक प्रयोग दिखाते हैं कि प्रस्तावित विधि अनुमानित अभिसरण दर प्राप्त करती है और उन विसंगतियों को सफलतापूर्वक हल करती है जो नाइव मिश्रित सूत्रीकरणों में देखी जाती हैं।

इस कार्य का महत्व चतुर्थ-क्रम की समस्याओं के लिए $C^0$ परिमित तत्व विधियों (जो मानक निरंतर टुकड़ों वाले बहुपदों का उपयोग करते हैं) की गणनात्मक दक्षता को सामान्य सीमा स्थितियों के साथ विस्तारित करने की क्षमता में निहित है, जो पहले $C^1$ तत्वों की जटिलता या सरल विघटन रणनीतियों की अभिसरण विफलताओं के कारण बाधित थी। सीमा स्थितियों और डोमेन ज्यामिति के विशिष्ट प्रभाव को अभिलक्षित करके, लेखक एक मजबूत एल्गोरिथम ढांचा प्रदान करते हैं जो कार्यान्वयन की सरलता से समझौता किए बिना गणितीय शुद्धता सुनिश्चित करता है। यह शोध पत्र एक श्रृंखला के भाग I के रूप में कार्य करता है, जो आगामी कार्यों में निम्न-क्रम के पदों वाली अधिक जटिल चतुर्थ-क्रम समीकरणों को संभालने के लिए सैद्धांतिक आधार तैयार करता है।