Fully discrete analysis of pressure-robust enriched Galerkin methods for the time-dependent Stokes equations

तकनीकी सारांश: समय-निर्भर स्टोक्स समीकरणों के लिए पूर्णतः विविक्त (Fully Discrete) दबाव-मजबूत समृद्ध गैलरकिन (Enriched Galerkin) विधियों का विश्लेषण

समस्या विवरण
यह शोध पत्र समय-निर्भर स्टोक्स समीकरणों के संख्यात्मक समाधान को संबोधित करता है, जो श्यान (viscous), असंपीड्य तरल प्रवाह को नियंत्रित करते हैं। शास्त्रीय मिश्रित परिमित तत्व विधियों (जैसे टेलर-हुड या क्रौज़िक-रिवियर्ट) में एक प्राथमिक चुनौती "दबाव-मजबूती" (pressure-robustness) का अभाव है। मानक स्वरूपों में, वेग त्रुटि सीधे निरंतर दबाव और बल (forcing term) के अघूर्णी घटक (irrotational component) पर निर्भर करती है। फलस्वरूप, जब दबाव अधिक होता है या उसमें तीव्र प्रवणता (gradients) होती है, तो वेग सन्निकटन की सटीकता काफी घट सकती है, भले ही सटीक वेग सुचारू (smooth) हो। यह समस्या कम-श्यानता वाले शासन (low-viscosity regimes) में विशेष रूप से गंभीर हो जाती है, जहाँ त्रुटि सीमाओं में श्यानता पैरामीटर ($\mu^{-1}$) के व्युत्क्रम घात शामिल होते हैं, जिससे अवास्तविक परिणाम प्राप्त होते हैं। जबकि स्थिर-अवस्था समस्याओं और विभिन्न अन्य स्वरूपों (जैसे $H(\text{div})$-सम्मत, वर्चुअल एलिमेंट्स) के लिए दबाव-मजबूत विधियाँ मौजूद हैं, समृद्ध गैलरकिन (EG) विधियों का उपयोग करके समय-निर्भर स्टोक्स प्रणालियों के लिए एक पूर्णतः विविक्त, पैरामीटर-स्पष्ट विश्लेषण—विशेष रूप से जिसमें समय-व्युत्पन्न पद में वेग पुनर्निर्माण (velocity reconstruction) को शामिल किया गया हो—स्थापित नहीं किया गया है।

कार्यप्रणाली
लेखक एक पूर्णतः विविक्त समृद्ध गैलरकिन (EG) ढांचा प्रस्तावित करते हैं जो निरंतर गैलरकिन (CG) विधियों की गणनात्मक दक्षता को विच्छिन्न गैलरकिन (DG) विधियों के स्थानीय द्रव्यमान संरक्षण गुणों के साथ जोड़ता है।

1. स्थानिक विवक्तीकरण (Spatial Discretization): वेग स्थान $V_h$ का निर्माण एक निरंतर खंडवार रैखिक लैग्रेंज स्पेस ($C_h$) को एक विच्छिन्न, खंडवार रैखिक, शून्य-माध्य मान वाले वेक्टर स्पेस ($D_h$) के साथ समृद्ध करके किया जाता है। दबाव स्थान $Q_h$ खंडवार स्थिरांकों (piecewise constants) से बना है। सख्त स्थानीय द्रव्यमान संरक्षण को लागू करने और दबाव-मजबूती प्राप्त करने के लिए, एक वेग पुनर्निर्माण ऑपरेटर $R: V_h \to BDM_1(\mathcal{T}_h)$ का उपयोग किया जाता है। यह ऑपरेटर EG आधार के निरंतर घटक को संरक्षित करता है और विच्छिन्न संवर्धन को $H(\text{div})$ के निम्नतम-क्रम वाले राविए-थॉमस ($RT_0$) उप-स्थान में मैप करता है।
2. कालिक विवक्तीकरण (Temporal Discretization): लेखक बैकवर्ड यूलर (प्रथम-क्रम) और क्रैंक-निकोल्सन (द्वितीय-क्रम) दोनों विधियों का उपयोग करके पूर्णतः विविक्त स्कीमें विकसित करते हैं। एक महत्वपूर्ण पद्धतिगत नवाचार यह है कि पुनर्निर्माण ऑपरेटर $R$ को न केवल बल पद में, बल्कि विविक्त समय-व्युत्पन्न पद में भी शामिल किया गया है। विशेष रूप से, मानक अंतर भाग $(u_h^{m+1} - u_h^m)/\tau$ को $(R u_h^{m+1} - R u_h^m)/\tau$ द्वारा प्रतिस्थापित किया जाता है।
3. विश्लेषण: लेखक बिना शर्त स्थिरता (unconditional stability) स्थापित करते हैं और 'ए प्रियोरी' (a priori) त्रुटि अनुमान प्राप्त करते हैं। विश्लेषण स्पष्ट रूप से स्थानिक जाल आकार ($h$), समय-चरण आकार ($\tau$), और श्यानता ($\mu$) पर निर्भरता को ट्रैक करता है।

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

• पूर्णतः दबाव-मजबूत EG स्वरूप: शोध पत्र एक दबाव-मजबूत EG स्थानिक विवक्तीकरण प्रस्तुत करता है जो एक $H(\text{div})$-सम्मत स्थान में समृद्ध विच्छिन्न घटक को मैप करके सख्त स्थानीय द्रव्यमान संरक्षण प्राप्त करता है।
• समय-व्युत्पन्न में पुनर्निर्माण: लेखक प्रदर्शित करते हैं कि समय-निर्भर समस्याओं के लिए पूर्णतः विविक्त स्तर पर दबाव-मजबूती बनाए रखने के लिए, वेग पुनर्निर्माण को केवल बल पद में ही नहीं, बल्कि विविक्त समय-व्युत्पन्न पद में भी लागू किया जाना चाहिए। यह क्षणिक (transient) समस्याओं के लिए स्थिर दबाव-मजबूत तकनीकों को विस्तारित करने की एक सीमा को सुधारता है।
• इष्टतम-क्रम त्रुटि अनुमान: लेखक बिना शर्त स्थिरता सिद्ध करते हैं और इष्टतम-क्रम, पैरामीटर-स्पष्ट 'ए प्रियोरी' त्रुटि अनुमान प्राप्त करते हैं। महत्वपूर्ण रूप से, वेग त्रुटि सीमाएं निरंतर दबाव और बल के अघूर्णी भाग से स्वतंत्र हैं, और इनमें कोई व्युत्क्रम-श्यानता कारक ($\mu^{-1}$) नहीं है।
• द्रव्यमान संरक्षण: पुनर्निर्मित वेग प्रत्येक तत्व पर बिंदुवार (pointwise) अपसरण-मुक्त (divergence-free) दिखाया गया है, जो प्रत्येक समय चरण पर सख्त स्थानीय द्रव्यमान संरक्षण सुनिश्चित करता है।

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

• अभिसरण दरें (Convergence Rates): स्कीमा अनुमानित अभिसरण दर प्राप्त करती हैं। $\tau \sim h$ के साथ बैकवर्ड युलर योजना के लिए, वेग के लिए प्रथम-क्रम अभिसरण देखा जाता है। $\tau \sim h^{1/2}$ के साथ क्रैंक-निकोल्सन योजना के लिए, वेग त्रुटि समय में द्वितीय-क्रम के करीब पहुँचते हुए अभिसरण करती है (प्रभावी रूप से युग्मन के आधार पर $O(h)$ या बेहतर), जबकि मजबूती बनाए रखती है।
• श्यानता मजबूती (Viscosity Robustness): कम-श्यानता शासन ($\mu = 10^{-6}$ से $10^{-10}$) में, प्रस्तावित स्कीमा सटीकता बनाए रखती है, जबकि मानक स्कीमा या वे जिनमें पुनर्निर्माण केवल बल पद पर लागू किया गया है (SR-BE), इष्टतम अभिसरण दर प्राप्त करने में विफल रहते हैं या त्रुटि क्षरण प्रदर्शित करते हैं।
• द्रव्यमान संरक्षण: 3D परीक्षणों में, पुनर्निर्मित वेग मशीन परिशुद्धता ($\approx 10^{-14}$) तक अपसरण-मुक्त है, जबकि कच्चा विविक्त वेग एक गैर-शून्य अपसरण दोष (divergence defect) रखता है जो समय के साथ संचित नहीं होता है।
• व्यावहारिक अनुप्रयोग: विधि को एक आयताकार बाधा के चारों ओर 3D प्रवाह पर सफलतापूर्वक लागू किया गया है, जो मिश्रित सीमा स्थितियों और कम श्यानता के साथ स्थिरता और जटिल प्रवाह विशेषताओं (पृथक्करण, पुनर्संयोजन, दबाव गिरावट) को पकड़ने की क्षमता प्रदर्शित करता है।

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