Conforming/Non-conforming Virtual Elements and application to elasticity problems in curved three-dimensional domains

तकनीकी सारांश: अनुरूप/अनुरूप-रहित आभासी तत्व और वक्रित त्रि-आयामी डोमेन में प्रत्यास्थता समस्याओं का अनुप्रयोग

समस्या विवरण
वर्चुअल एलीमेंट मेथड (VEM) ने सामान्य बहुफलकीय मेष (polytopal meshes) पर आंशिक अंतर समीकरणों (PDEs) के विवेकीकरण (discretization) के लिए एक लचीला ढांचा स्थापित किया है। हालाँकि, वक्रित डोमेन वाले त्रि-आयामी (3D) समस्याओं में VEM को लागू करने में महत्वपूर्ण चुनौतियाँ बनी हुई हैं। 2D VEM, जो वक्रित सीमाओं के साथ विकसित हो चुका है, की तुलना में 3D कार्यान्वयन वक्रित सतहों पर आभासी स्थानों (virtual spaces) को परिभाषित करने और आवश्यक प्रक्षेपों (projections) की गणना करने में कठिनाइयों का सामना करते हैं। मौजूदा दृष्टिकोण अक्सर जटिल मास्टर-स्लेव निर्माणों, सीमा स्थितियों के कमजोर प्रवर्तन (weak enforcement), या विच्छिन्न चरों (discontinuous variables) वाले मिश्रित सूत्रीकरणों पर निर्भर करते हैं। ठोस यांत्रिकी (solid mechanics) के संदर्भ में, ऐसे तरीके जो विच्छिन्न विस्थापन क्षेत्रों (displacement fields) को जन्म देते हैं, भौतिक रूप से अवास्तविक होते हैं। इसके अतिरिक्त, सीधी भुजाओं वाले मानक पॉलीहेड्रल मेष, ज्यामितीय सन्निकटन त्रुटियों (geometric approximation errors) को पेश करते हैं जो उच्च-क्रम योजनाओं के समाधान त्रुटि पर हावी हो सकती हैं।

कार्यप्रणाली
लेखक एक नवीन अनुरूप/अनरूप-रहित (Conforming/Non-conforming - C-NC) वर्चुअल एलीमेंट मेथड प्रस्तावित करते हैं जो उपयोगकर्ता को मेष फलकों (mesh faces) को अनुरूप या अनरूप-रहित के रूप में चुनिweise रूप से नामित करने की अनुमति देता है।

• हाइब्रिड स्पेस निर्माण: विधि एक पॉलीहेड्रॉन $E$ पर एक स्थानीय आभासी स्थान $W^k_h(E)$ को इसके सीमा फलकों को एक अनुरूप सेट ($\mathcal{F}^c_E$) और एक अनरूप-रहित सेट ($\mathcal{F}^n_E$) में विभाजित करके परिभाषित करती है।
  • अनरूप फलकों पर, समाधान निरंतर ($C^0$) होता है, और यह स्थान मानक नोडल VEM परिभाषा का अनुसरण करता है (जो बहुपद क्षणों/polynomial moments के साथ संवर्धित है)।
  • अनरूप फलकों पर, समाधान को विच्छिन्न होने की अनुमति है, और अभिलंब अवकलज (normal derivative) को $k-1$ घात के बहुपद के रूप में प्रतिबंधित किया जाता है।
  • वैश्विक स्थान को केवल अनुरूप फलकों के माध्यम से निरंतरता लागू करके निर्मित किया जाता है। यह "हाइब्रिड" तत्व स्थानीय पॉइसन समस्याओं (Poisson problems) के माध्यम से बनाया जाता है जिनमें मिश्रित डिरिचलेट (अनरूप फलकों पर) और न्यूमैन (अनरूप फलकों पर) सीमा स्थितियाँ होती हैं।
• डिग्री ऑफ फ्रीडम (DoFs): स्थान को शीर्षों (vertices) और किनारों के बिंदुओं (अनरूप सीमाओं पर) पर बिंदुवार मूल्यांकन और फलन तथा उसके अभिलंब अवकलज के क्षणों द्वारा परिभाषित किया जाता है। महत्वपूर्ण रूप से, यह दिखाया गया है कि स्टिफनेस मैट्रिक्स (विशेष रूप से $L^2$ और ऊर्जा प्रक्षेपों) के लिए आवश्यक प्रक्षेप, इन DoFs से ही गणनीय हैं, भले ही सीमा स्थितियाँ मिश्रित हों।
• वक्रित डोमेन के लिए अनुप्रयोग: विधि को वक्रित डोमेन पर 3D रैखिक प्रत्यास्थता (linear elasticity) के लिए लागू किया गया है। रणनीति में शामिल है:
  • डोमेन के आंतरिक भाग के लिए समतल, अनुरूप फलकों का उपयोग करना ताकि वैश्विक विस्थापन निरंतरता सुनिश्चित हो सके।
  • सटीक ज्यामिति का प्रतिनिधित्व करने के लिए विशेष रूप से सीमा $\partial\Omega$ पर वक्रित, अनरूप फलकों का उपयोग करना।
  • एक सुचारू पैरामीट्राइजेशन $\gamma_F$ के माध्यम से मैप्ड बहुपदों का उपयोग करके वक्रित फलकों पर आभासी स्थान को परिभाषित करना।
  • प्रक्षेप ऑपरेटरों पर सीमा समाकल (boundary integrals) के सन्निकटन का उपयोग करना, यह स्वीकार करते हुए कि यह एक आवश्यक सन्निकटन है जब प्रक्षेप का अभिलंब अवकलज मैप्ड बहुपद स्थान में नहीं होता है।

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

1. नवीन हाइब्रिड तत्व: यह शोध पत्र एक ऐसा VEM सूत्रीकरण प्रस्तुत करता है जहाँ फलकों की अनुरूप/अनरूप प्रकृति एक उपयोगकर्ता-निर्धारित विकल्प है, जो एक "हाइब्रिड" तत्व को सक्षम बनाता है जो शास्त्रीय अनुरूप VEM, अनरूप VEM, या दोनों के मिश्रण को पुनः प्राप्त कर सकता है।
2. सैद्धांतिक विश्लेषण: लेखक C-NC स्थान के लिए एक कठोर स्थिरता और अभिसरण विश्लेषण प्रदान करते हैं। वे $H^1$-सेमीनॉर्म और DoF-आधारित स्थिरीकरण के बीच समानता स्थापित करते हैं, इंटरपोलेशन अनुमान प्राप्त करते हैं, और रैखिक प्रत्यास्थता समस्या के लिए इष्टतम त्रुटि सीमाएँ (error bounds) निकालते हैं।
3. वक्रित 3D प्रत्यास्थता: यह शोध पत्र सटीक ज्यामिति प्रतिनिधित्व के साथ वक्रित डोमेन पर 3D प्रत्यास्थता समस्याओं को हल करने के लिए एक व्यावहारिक कार्यप्रणाली प्रदर्शित करता है। पिछले 3D वक्रित VEM दृष्टिकोणों के विपरीत, जो विच्छिन्न विस्थापन उत्पन्न करते हैं या जटिल निर्माणों की आवश्यकता रखते हैं, यह विधि बल्क में एक निरंतर विस्थापन क्षेत्र की गारंटी देती है और साथ ही सटीक वक्रित सीमा को भी कैप्चर करती है।
4. संख्यात्मक सत्यापन: विधि का परीक्षण तीन परिदृश्यों पर किया गया है:
   • मानक सीधी ज्यामिति (अभिसरण दरों को मान्य करना)।
   • सटीक वक्रित सीमा प्रतिनिधित्व (सिलेंडर) के साथ वक्रित ज्यामिति।
   • रैखिक और द्विघाती ज्यामितीय सन्नकनों के साथ मनमाने, गैर-समतलीय बहुभुजों द्वारा टेसेलेटेड गोले (sphere) के उदाहरण।

परिणाम
संख्यात्मक प्रयोग सैद्धांतिक भविष्यवाणियों की पुष्टि करते हैं:

• अभिसरण दरें (Convergence Rates): समतल और वक्रित दोनों डोमेन के लिए, विधि इष्टतम अभिसरण दर प्राप्त करती है: बहुपद घात $k$ के लिए $L^2$-नॉर्म में $O(h^{k+1})$ और $H^1$-सेमीनॉर्म में $O(h^k)$।
• वक्रित सीमा प्रदर्शन: सटीक वक्रित सीमा प्रतिनिधित्व वाले सिलेंडर के मामले में, विधि इष्टतम दर बनाए रखती है, जो टेट्राहेड्रल मेष पर मानक द्वितीय-क्रम परिमित तत्व (finite element) योजनाओं से बेहतर प्रदर्शन करती है, जो ज्यामितीय सन्निकटन त्रुटियों ($O(h^2)$ और $O(h^{3/2})$) से ग्रस्त होती हैं।
• ज्यामितीय सन्निकटन सीमाएँ: गोले के उदाहरण के लिए, विधि मनमाने पॉलीगोनल टेसेलेशन को सफलतापूर्वक संभालती है। हालाँकि, लेखक देखते हैं कि $k=2$ के लिए, इष्टतम अभिसरण दर तभी प्राप्त होती है जब ज्यामितीय सन्निकटन पर्याप्त रूप से सटीक हो। द्विघाती ज्यामितीय सन्नकन (Quad) का उपयोग करने पर, $L^2$ त्रुटि प्रारंभ में अभिसरण का उच्च क्रम दिखाती है लेकिन अंततः स्थिर (plateau) हो जाती है। लेखक इसका श्रेय ज्यामितीय सन्निकटन त्रुटि (विशेष रूप से सतह क्षेत्र और आयतन त्रुटि) को देते हैं, जो प्रभावी हो जाती है, यह नोट करते हुए कि सेरेंडिपिटी आधार (serendipity basis) $k=2$ तत्वों के साथ इष्टतम $L^2$ अभिसरण के लिए आवश्यक $O(h^3)$ ज्यामितीय त्रुटि की गारंटी नहीं देता है।

महत्व और दावे
शोध पत्र का दावा है कि प्रस्तावित C-NC VEM वक्रित डोमेन पर 3D प्रत्यास्थता समस्याओं के लिए एक लचीला और मजबूत ढांचा प्रदान करता है। इसका प्राथमिक महत्व एक निरंतर विस्थापन समाधान प्रदान करने (जो ठोस यांत्रिकी के लिए आवश्यक है) की क्षमता में निहित है, जबकि अनरूप फलकों का उपयोग करके सटीक वक्रित सीमाओं का प्रतिनिधित्व किया जाता है—एक ऐसा संयोजन जिसे मानक VEM या FEM दृष्टिकोणों के माध्यम से प्राप्त करना कठिन था। लेखक इस बात पर जोर देते हैं कि हालांकि विधि वक्रित फलकों पर सीमा समाकल के प्रक्षेप में एक सन्निकटन पेश करती है, लेकिन यह विधि की स्थिरता या इष्टतम अभिसरण दरों से समझौता नहीं करती है, बशर्ते ज्यामितीय सन्निकटन पर्याप्त हो। यह कार्य जटिल 3D ज्यामिति को गैर-भौतिक विच्छिन्न समाधानों या अत्यधिक जटिल मास्टर-स्लेव निर्माणों के बिना संभालने के लिए एक प्रमाण-अवधारणा (proof of concept) के रूप में कार्य करता है।