A Staggered Discontinuous Galerkin Method for linear elasticity problem on Polytopal Meshes

तकनीकी सारांश: पॉलीटोपल मेश पर रैखिक लोच (Linear Elasticity) के लिए एक स्टैगर्ड डिस्कंटीन्यूअस गैलेरकिन विधि

समस्या विवरण
यह शोध पत्र हेलिंगर-रिएस्नर मिश्रित परिवर्तनशील सिद्धांत (Hellinger-Reissner mixed variational principle) का उपयोग करके रैखिक लोच की संख्यात्मक सन्निकटन (numerical approximation) समस्या को संबोधित करता है। इसका लक्ष्य सममित प्रतिबल टेंसर (symmetric stress tensor) $\sigma \in H(\text{div}, \Omega; S)$ और विस्थापन क्षेत्र (displacement field) $u \in L^2(\Omega; \mathbb{R}^d)$ का एक साथ सन्निकटन करना है। इस मिश्रित सूत्रीकरण (mixed formulation) के विवेकीकरण (discretization) में मुख्य चुनौतियाँ तीन प्रकार की हैं: (i) स्थिरता सुनिश्चित करने के लिए इन्फ-सुप (inf-sup) स्थिति को संतुष्ट करना; (ii) प्रतिबल टेंसर की समरूपता (symmetry) और उसके तत्वों की सीमाओं पर अभिलंब निरंतरता (normal continuity) दोनों को लागू करना; और (iii) इष्टतम अभिसरण दरें (optimal convergence rates) प्राप्त करना। मौजूदा विधियाँ अक्सर वर्टेक्स डिग्री ऑफ फ्रीडम (DoFs) पर निर्भर करती हैं, जो हाइब्रिडाइजेशन में बाधा डालती हैं, या कमजोर समरूपता (weak symmetry) का उपयोग करती हैं, जो स्थानीय रूप से कोणीय संवेग के संतुलन को सख्ती से संतुष्ट नहीं कर सकती हैं। इसके अलावा, कई विधियाँ असंपीड्य सीमा ($\lambda \to \infty$) में "लॉकिंग" (locking) की समस्या से जूझती हैं।

कार्यप्रणाली (Methodology)
लेखक एक नवीन स्टैगर्ड डिस्कंटीन्यूअस गैलेरकिन (SDG) विधि प्रस्तावित करते हैं जो सामान्य पॉलीटोपल मेश पर परिभाषित है। यह एक प्राइमल-डुअल मेश संरचना का उपयोग करता है:

1. मेश निर्माण: एक प्राइमल मेश $K_h$ सामान्य पॉलीटोपल तत्वों से बना है। एक डुअल मेश $K^*_h$ का निर्माण प्रत्येक प्राइमल तत्व के बैरीसेंटर (barycenter) को उसके शीर्षों (vertices) से जोड़कर एक परिष्कृत त्रिकोणीकरण (refined triangulation) $T_h$ के माध्यम से किया जाता है।
2. स्टैगर्ड निरंतरता (Staggered Continuity): यह विधि पिछले SDG लोच कार्यों में पाए जाने वाले निरंतरता असाइनमेंट को उलट देती है।
   • प्रतिबल ($\sigma_h$): परिष्कृत त्रिकोणीकरण $T_h$ पर खंडित बहुपदों (piecewise polynomials) के रूप में परिभाषित है। महत्वपूर्ण रूप से, प्रतिबल दृढ़ रूप से सममित (strongly symmetric) है और प्राइमल तत्वों ($K \in K_h$) की सीमाओं पर अभिलंब निरंतरता (normal continuity) रखता है।
   • विस्थापन ($u_h$): प्राइमल तत्वों $K_h$ पर खंडित बहुपदों द्वारा अनुमानित किया जाता है। विस्थापन डुअल तत्वों ($\omega_F \in K^*_h$) की सीमाओं पर निरंतर है।
3. परिमित तत्व स्थान (Finite Element Spaces):
   • प्रतिबल स्थान $\Sigma_h$ को दृढ़ समरूपता के साथ $H(\text{div})$-अनुरूप बनाया गया है। स्थानीय DoFs को प्राइमल फलक (faces) पर अभिलंब फ्लक्स और आंतरिक मोमेंट्स के माध्यम से परिभाषित किया जाता है, जिससे वर्टेक्स-आधारित अज्ञात चर (unknowns) से बचा जा सके।
   • विस्थापन स्थान $U_h$ प्राइमल सेल पर वेक्टर-मान वाले बहुपदों का समूह है।
   • निम्नतम क्रम ($k=0$) के लिए, सुव्यवस्थितता (well-posedness) सुनिश्चित करने के लिए प्रतिबल स्थान को रिजिड मोशन स्पेस के साथ समृद्ध किया गया है।
4. परिवर्तनशील सूत्रीकरण (Variational Formulation): योजना को रैखिक लोच के प्रथम-क्रम सिस्टम से व्युत्पन्न किया गया है। इसमें प्रतिबल के डाइवर्जेंस और परिष्कृत त्रिकोणीकरण के आंतरिक फलकों पर अभिलंब फ्लक्स के जंप (jumps) से संबंधित स्थिरीकरण पद (stabilization terms) शामिल हैं। ये पद संख्यात्मक लॉकिंग को रोकने के लिए डिज़ाइन किए गए हैं, न कि निरंतरता लागू करने के लिए।
5. हाइब्रिडाइजेशन (Hybridization): यह विधि एक हाइब्रिडाइज़ेबल संरचना को स्वीकार करती है। तत्व सीमाओं पर लैग्रेंज मल्टीप्लायर्स (Lagrange multipliers) को पेश करके, प्रतिबल चर को स्टैटिकली कंडेंस (statically condensed) किया जा सकता है, जिसके परिणामस्वरूप केवल विस्थापन के सीमा DoFs से संबंधित एक वैश्विक प्रणाली प्राप्त होती है।

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

• पॉलीटोप्स पर दृढ़ समरूपता (Strong Symmetry on Polytopes): यह शोध पत्र वर्टेक्स DoFs पर निर्भर किए बिना, सामान्य पॉलीटोपल मेश पर एक विशिष्ट ज्यामितीय अपघटन और DoFs के सावधानीपूर्ण चयन के माध्यम से सममित प्रतिबल स्थानों का निर्माण करता है।
• उलटी निरंतरता (Reversed Continuity): पिछली SDG लोच विधियों [31, 41] के विपरीत, जो कमजोर समरूपता लागू करती हैं, यह कार्य प्राइमल और डुअल मेश के सापेक्ष प्रतिबल और विस्थापन के निरंतरता पैटर्न को बदलकर दृढ़ समरूपता लागू करता है।
• लॉकिंग-मुक्त गुण (Locking-Free Property): इस विधि को लॉकिंग-मुक्त सिद्ध किया गया है। स्थिरता और अभिसरण असंपीड्य सीमा ($\lambda \to \infty$) में भी मजबूत रहता है क्योंकि स्थिरीकरण पद लैमे मापदंडों (Lamé parameters) से स्वतंत्र रूप से प्रतिबल के ट्रेस (trace) को नियंत्रित करते हैं।
• स्थानीय संरक्षण (Local Conservation): यह सूत्र रैखिक संवेग के स्थानीय संतुलन (डाइवर्जेंस समीकरण के माध्यम से) और कोणीय संवेग (प्रतिबल टेंसर की दृढ़ समरूपता के माध्यम से) को संरक्षित करता है।
• हाइब्रिडाइज़ेबल संरचना: वर्टेक्स DoFs की अनुपस्थिति एक हाइब्रिडाइजेशन प्रक्रिया की अनुमति देती है जो वैश्विक रैखिक प्रणाली के आकार को महत्वपूर्ण रूप से कम करती है, जिससे कार्यान्वयन सरल होता है और गणना लागत कम होती है।

परिणाम

• सैद्धांतिक विश्लेषण: लेखक एक व्यापक सैद्धांतिक विश्लेषण प्रदान करते हैं, जो डिस्क्रीट इन्फ-सुप स्थितियों और कोएर्सिविटी (coercivity) के माध्यम से योजना की सुव्यवस्थितता (well-posedness) को सिद्ध करते हैं। इष्टतम a priori त्रुटि अनुमान प्राप्त किए गए हैं, जो ऊर्जा नॉर्म में प्रतिबल के लिए $O(h^{k+1})$ और डिस्क्रीट $H^1$-समान नॉर्म में विस्थापन के लिए अभिसरण दर दिखाते हैं, और $L^2$ नॉर्म में विस्थापन के लिए $O(h^{k+2})$ दिखाते हैं।
• संख्यात्मक प्रयोग:
  • लॉकिंग-मुक्त सत्यापन: विभिन्न लैमे मापदंडों ($\lambda = 1, 10^2, 10^4, 10^6$) के साथ त्रिकोणीय मेश पर संख्यात्मक परीक्षणों से पता चलता है कि त्रुटि $\lambda$ से स्वतंत्र रहती है, जो विधि की लॉकिंग-मुक्त प्रकृति की पुष्टि करती है।
  • अभिसरण दरें: त्रिकोणीय और सामान्य बहुभुज (polygonal) मेश दोनों पर प्रयोग $k=0$ और $k=1$ के लिए सैद्धांतिक अभिसरण क्रमों की पुष्टि करते हैं।
  • पॉलीटोपल लचीलापन: विधि ने सामान्य बहुभुज मेश को सफलतापूर्वक संभाला है, जो मानक सिम्पलीशियल या चतुर्भुज ग्रिड से परे इसकी प्रयोज्यता को प्रदर्शित करता है।

महत्व
यह शोध पत्र दावा करता है कि यह कार्य पॉइसन समीकरण [16] पर पहले से लागू SDG ढांचे को रैखिक लोच के अधिक जटिल क्षेत्र तक विस्तारित करता है। इसका महत्व इस तथ्य में निहित है कि यह एक ऐसी विधि प्रदान करता है जो दृढ़ समरूपता, अभिलंब निरंतरता, और हाइब्रिडाइज़ेबिलिटी को सामान्य पॉलीटोपल मेश पर वर्टेक्स डिग्री ऑफ फ्रीडम के बिना जोड़ती है। यह लोच के लिए मिश्रित परिमित तत्व विधियों (mixed finite element methods) की एक दीर्घकालिक कठिनाई को संबोधित करता है, जहाँ समरूपता और निरंतरता को एक साथ संतुष्ट करना अक्सर जटिल तत्वों की ओर ले जाता है जिनमें वर्टेक्स DoFs होते हैं या जिनके लिए कमजोर समरूपता की आवश्यकता होती है। प्रस्तावित विधि एक मजबूत, लॉकिंग-मुक्त विकल्प प्रदान करती है जो सैद्धांतिक रूप से सुदृढ़ है और अपने हाइब्रिडाइज़ेबल ढांचे के कारण गणनात्मक रूप से कुशल है।