A Divergence-Free Scott-Vogelius Finite Element Method for the Surface Stokes Problem

तकनीकी सारांश: सतह स्टोक्स समस्या के लिए एक डाइवर्जेंस-फ्री स्कॉट-वोगेलियसीय परिमित तत्व विधि (Finite Element Method)

समस्या विवरण
यह शोध पत्र एक चिकनी, बंद, उन्मुखी योग्य (orientable) द्वि-आयामी सतह $\gamma \subset \mathbb{R}^3$ पर स्थिर सतह स्टोक्स समस्या (stationary surface Stokes problem) को संबोधित करता है। इसके शासी समीकरणों में एक स्पर्शरेखीय वेग क्षेत्र (tangential velocity field) $u$ और एक दबाव $p$ शामिल हैं, जो इस प्रतिबंध के अधीन हैं कि $u$ सतह के स्पर्शरेखीय है ($u \cdot \nu = 0$) और प्रवाह असंपीड्य (incompressible) है ($\text{div}_\gamma u = 0$)।

इस समस्या के लिए परिमित तत्व विधियों का निर्माण दो प्राथमिक चुनौतियाँ प्रस्तुत करता है:

1. स्पर्शरेखीयता (Tangentiality): एक विविक्त ज्यामिति (discrete geometry) पर स्पर्शरेखीय प्रतिबंध को कड़ाई से लागू करने से मानक $H^1$-कन्फॉर्मिंग वेक्टर स्पेस का उपयोग करते समय "लॉकिंग" (locking) की घटना हो सकती है, जिससे सन्निकटन गुण (approximation properties) खराब हो जाते हैं।
2. असंपीड्यता (Incompressibility): मानक वेग-दबाव युग्म (जैसे टेलर-हुड) केवल विवर्तन-मुक्त (divergence-free) प्रतिबंध को दुर्बल रूप से (weakly) लागू करते हैं। जबकि ऐसी विधियाँ मौजूद हैं जो स्पर्शरेखीयता को कड़ाई से लागू करती हैं (जैसे $H(\text{div})$-कन्फॉर्मिंग BDM स्पेस), वे अक्सर बड़े बीजगणितीय सिस्टम का परिणाम देती हैं या सतहों पर पूर्ण स्थिरता विश्लेषण का अभाव होता है। इसके विपरीत, वे विधियाँ जो स्पर्शरेखीयता को कड़ाई से लागू करते हुए $H^1$-कन्फॉर्मिटी बनाए रखती हैं (जैसे सतहों पर MINI या टेलर-हुड), आमतौर पर असंपीड्यता को केवल दुर्बल रूप से ही लागू करती हैं।

लेखक एक ऐसी विधि का निर्माण करने का लक्ष्य रखते हैं जो स्पर्शरेखीयता और असंपीड्यता दोनों प्रतिबंधों को सटीक रूप से लागू करती है, मानक 2D यूक्लिडियन विविक्तीकरणों की आयामीता (dimensionality) को बनाए रखती है, और घुमावदार सतहों पर मैक्रो-एलिमेंट निर्माण की कार्यान्वयन जटिलताओं से बचती है।

कार्यप्रणाली (Methodology)
प्रस्तावित दृष्टिकोण Kone et al. [10, 11] के सरफेस फाइनाइट एलीमेंट फ्रेमवर्क को स्कॉट-वोगेलियस युग्म (Scott-Vogelius pair) तक विस्तारित करता है। कार्यप्रणाली के मुख्य घटक निम्नलिखित हैं:

• त्रिकोणीयकरण (Triangulation): विधि एक पॉलीहेड्रल सन्निकटन वाली सतह के क्लफ-टोचर (Clough-Tocher - CT) शोधन का उपयोग करती है। प्रारंभिक मेश के प्रत्येक त्रिभुज को उसके शीर्षों से उसके केंद्रक (barycenter) को जोड़कर विभाजित किया जाता है। यह शोधन स्थिरता के लिए बहुपद डिग्री की आवश्यकताओं को शिथिल करता है (यह अनुमति देता है कि $r \ge 2$) और "लगभग विलक्षण" (nearly singular) शीर्षों को समाप्त करता है।
• परिमित तत्व स्थान (Finite Element Spaces):
  • वेग ($V_h$): CT-परिष्कृत सतह त्रिकोणीयकरण पर डिग्री $r$ के निरंतर, खंडित बहुपदों (piecewise polynomials) के रूप में परिभाषित है। सतह पर $H(\text{div})$-कन्फॉर्मिटी के लिए खंडों के बीच दुर्बल निरंतरता (weak continuity) को लागू करने हेतु, लेखक लैग्रेंज डिग्री ऑफ फ्रीडम (DOFs) पर एक सहायक पियोला ट्रांसफॉर्म (Piola transform) का उपयोग करते हैं। यह इंटर-एलिमेंट इंटीग्रेशन या स्थिरीकरण पदों की आवश्यकता के बिना स्थानीय स्थानों को आपस में जोड़ता है।
  • दबाव ($Q_h$): उसी CT-परिष्कृत मेश पर खंडित बहुपदों के डिग्री $r-1$ के रूप में परिभाषित है।
• कार्यान्वयन रणनीति: मैक्रो-एलिमेंट परिभाषाओं (जहाँ आधार फलन एक विभाजित संदर्भ तत्व पर परिभाषित किए जाते हैं और वैश्विक रूप से मैप किए जाते हैं) पर आधारित पिछले आइसोपैरामेट्रिक स्कॉट-वोगेलियस विधियों के विपरीत, यह शोध पत्र परिमित तत्व स्थानों को सीधे परिष्कृत सतह त्रिकोणीयकरण पर परिभाषित करता है। यह दृष्टिकोण मैक्रो-एलिमेंट कार्यान्वयन के लिए आवश्यक जटिल डेटा संरचनाओं और विशिष्ट क्वाड्रचर नियमों से बचता है, जिससे एक अधिक व्यावहारिक योजना प्राप्त होती है जो मानक परिमित तत्व सॉफ्टवेयर के अनुकूल है।
• स्थिरता विश्लेषण (Stability Analysis): लेखक विविक्त इन्फ-सुप स्थिति (inf-sup condition या Ladyzhenskaya-Babuška-Brezzi स्थिति) को सिद्ध करते हैं। प्रमाण एक विक्षोभ तर्क (perturbation argument) का उपयोग करता है:
  1. स्थिरता को पहले एक थोड़े अलग उच्च-क्रम त्रिकोणीयकरण पर परिभाषित एक सहायक "मैक्रो-एलिमेंट" युग्म के लिए स्थापित किया जाता है।
  2. प्रस्तावित प्रत्यक्ष-शोधन स्थान और सहायक मैक्रो-एलिमेंट स्थान के बीच के अंतर को $O(h)$ विक्षोभ के रूप में मापा जाता है।
  3. इस विक्षोभ विश्लेषण के माध्यम से स्थिरता को सहायक युग्म से प्रस्तावित युग्म में स्थानांतरित किया जाता है।

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

1. सटीक डाइवर्जेंस-फ्री और स्पर्शरेखीय विविक्तीकरण: निर्मित क्लफ-टोचर घुमावदार त्रिकोणीयकरणों पर स्कॉट-वोगेलियस युग्म, स्पर्शरेखीयता ($u \cdot \nu = 0$) और असंपीड्यता ($\text{div}_\gamma u = 0$) दोनों प्रतिबंधों को सटीक रूप से लागू करता है।
2. सरलीकृत कार्यान्वयन: स्थानों को मैक्रो-एlement के बजाय सीधे परिष्कृत त्रिकोणीयकरण पर परिभाषित करके, लेखक एक ऐसी योजना प्रदान करते हैं जो स्कॉट-वोगेलियस युग्म के सैद्धांतिक लाभों को बनाए रखते हुए काफी सरल कार्यान्वयन प्रदान करती है।
3. स्थिरता और अभिसरण (Stability and Convergence):
   • विधि को मेश आकार से स्वतंत्र एक स्थिरांक के साथ इन्फ-सुप स्थिर सिद्ध किया गया है।
   • आइसोपैरामेट्रिक शासन ($k=r$) में इष्टतम-क्रम अभिसरण प्राप्त किया गया है। विशेष रूप से, वेग $H^1$-नॉर्म में $O(h^r)$ के क्रम से और $L^2$-नॉर्म में $O(h^{r+1})$ के क्रम से अभिसरित होता है, जबकि दबाव $L^2$-नॉर्म में $O(h^r)$ के क्रम से अभिसरित होता है।
   • विश्लेषण में सतह $P_r - P_0$ युग्म और संगत मैक्रो स्कॉट-वोगेलियस युग्म के लिए स्थिरता का प्रमाण भी शामिल है।
4. दबाव सुदृढ़ता विश्लेषण (Pressure Robustness Analysis): शोध पत्र दबाव सुदृढ़ता (जहाँ वेग त्रुटि दबाव सन्निकटन से स्वतंत्र होती है) की जांच करता है। लेखक तर्क देते हैं कि सतह स्टोक्स समस्याओं के लिए, ज्यामितीय सन्निकटन त्रुटियों और स्रोत पद के विस्तार के कारण, सटीक रूप से डाइवर्जेंस-फ्री विधियों के लिए भी दबाव सुदृढ़ता स्वचालित रूप से गारंटीकृत नहीं है। वे दबाव सुढ़ता को संरक्षित करने वाले स्रोत विस्तारों के निर्माण के लिए एक ढांचा प्रदान करते हैं और संख्यात्मक प्रयोगों के माध्यम से प्रदर्शित करते हैं कि मानक विस्तार इस गुण को नष्ट कर सकते हैं।

संख्यात्मक प्रयोग
लेखक यूनिट स्फीयर और टोरस पर तीन संख्यात्मक प्रयोगों के माध्यम से सिद्धांत को मान्य करते हैं:

• परीक्षण 1 (डाइवर्जेंस-फ्री प्रवाह): इष्टतम अभिसरण दर ($L^2$ वेग के लिए $O(h^{r+1})$, वेग ग्रेडिएंट और दबाव के लिए $O(h^r)$) की पुष्टि करता है और दिखाता है कि समाधान मशीन परिशुद्धता (machine precision) तक डाइवर्जेंस-फ्री है।
• परीक्षण 2 (नो-फ्लो/शून्य वेग): प्रदर्शित करता है कि जबकि दबाव इष्टतम रूप से अभिसरित होता है, ज्यामितीय त्रुटियों के कारण वेग सन्निकटन गैर-शून्य होता है, जो स्वचालित दबाव सुदृढ़ता की कमी को उजागर करता है। प्रयोग यह भी देखता है कि वेग $L^2$ त्रुटि में सुपरकन्वर्जेंस ($O(h^{r+2})$) होता है, जो यह सुझाव देता है कि सैद्धांतिक $L^2$ अनुमान सटीक (sharp) नहीं हो सकता है।
• परीक्षण 3 (टोरस): अधिक जटिल ज्यामिति पर इष्टतम अभिसरण दरों की पुष्टि करता है, हालांकि सिस्टम कंडीशनिंग और राउंड-ऑफ त्रुटियों के कारण कुछ गिरावट देखी गई है।

महत्व
यह शोध पत्र सतह द्रव गतिशीलता (surface fluid dynamics) के साहित्य में एक रिक्तता को भरने का दावा करता है, जो एक ऐसी विधि प्रदान करता है जो सटीक रूप से डाइवर्जेंस-फ्री विविक्तीकरणों (संरक्षण नियम, दबाव सुदृढ़ता क्षमता) के लाभों को मानक परिमित तत्व कार्यान्वयन की व्यावहारिकता के साथ जोड़ती है। मैक्रो-एलिमेंट निर्माण से दूर जाकर, लेखक एक ऐसा तंत्र प्रदान करते हैं जो स्कॉट-वोगेलियस युग्म की सैद्धांतिक कठोरता को बनाए रखते हुए "काफी सरल और अधिक व्यावहारिक कार्यान्वयन" प्रदान करता है। यह कार्य डाइवर्जेंस-फ्री विधियों की प्रयोज्यता को सतह समस्याओं तक विस्तारित करता है, जो लिपि बाइलेयर्स (lipid bilayers) और पतली फिल्मों (thin films) जैसी घटनाओं के मॉडलिंग के लिए महत्वपूर्ण हैं।