Error estimates for finite element discretizations of the instationary Navier-Stokes equations

समस्या विवरण
यह कार्य सजातीय डिरिचलेट (नो-स्लिप) सीमा स्थितियों के साथ द्वि-आयामी इंस्टेशनरी नेवियर-स्टोक्स समीकरणों के संख्यात्मक विश्लेषण को संबोधित करता है। समीकरणों को समय में एक डिस्कंटीन्यूअस गैलर्किन (dG) विधि और स्थान में इन्फ-सप (inf-sup) स्थिर परिमित तत्वों (finite elements) को संयोजित करते हुए एक पूर्णतः विविक्त (fully discrete) योजना का उपयोग करके विविक्त किया गया है। प्राथमिक उद्देश्य विविक्त समाधान के लिए त्रुटि अनुमान प्राप्त करना है, विशेष रूप से $L^\infty(I; L^2(\Omega))$ नॉर्म पर ध्यान केंद्रित करते हुए।

जबकि हाल ही में इंस्टेशनरी स्टोक्स समस्या के लिए $L^\infty(I; L^2(\Omega))$, $L^2(I; H^1(\Omega))$, और $L^2(I; L^2(\Omega))$ नॉर्म्स में 'बेस्ट एप्रोक्सिमेशन' प्रकार के त्रुटि अनुमान स्थापित किए गए हैं, इन परिणामों को गैररेखीय नेवियर-स्टोक्स समीकरणों तक विस्तारित करना महत्वपूर्ण चुनौतियाँ प्रस्तुत करता है। मौजूदा साहित्य अक्सर $L^\infty(I; L^2(\Omega))$ और $L^2(I; H^1(\Omega))$ में त्रुटियों का संयुक्त रूप से अनुमान लगाता है, जो स्थानिक $L^2$ त्रुटि अनुमान के लिए ऑर्डर रिडक्शन (क्रम में कमी) का कारण बनता है क्योंकि $H^1$ स्थानिक त्रुटि को ध्यान में रखना आवश्यक होता है। इसके अतिरिक्त, नेवियर-स्टोक्स के लिए पिछले परिणामों में अक्सर स्रोत पद $f$ (जैसे, $f \in W^{1,\infty}$ या $W^{2,\infty}$) पर मजबूत रेगुलैरिटी (नियमितता) धारणाओं की आवश्यकता होती है या वे उप-इष्टतम अभिसरण क्रम (जैसे, $O(\sqrt{k} + h)$) प्रदान करते हैं।

कार्यप्रणाली
लेखक त्रुट विभाजन (error splitting), द्वैत तर्क (duality arguments), और dG समय विवक्तीकरण के लिए अनुकूलित स्थिरता विश्लेषण के संयोजन का उपयोग करते हैं।

1. रेगुलैरिटी विश्लेषण: शोध पत्र पहले निरंतर नेवियर-स्टोक्स समाधान के लिए आवश्यक रेगुलैरिटी परिणाम स्थापित करता है। यह सिद्ध करता है कि उत्तल डोमेन (convex domain) की धारणाओं और उपयुक्त डेटा रेगुलैरिटी ($f \in L^2(I; L^2(\Omega)^2)$, $u_0 \in V$) के तहत, समाधान में $H^2$ स्थानिक रेगुलैरिटी होती है और गैररेखीय पद $(u \cdot \nabla)u$ पर्याप्त रेगुलैरिटी विरासत में प्राप्त करता है ताकि इसे स्टोक्स ऑपरेटर ढांचे में एक स्रोत पद के रूप में माना जा सके।
2. विवक्तीकरण (Discretization): स्थानिक विवक्तीकरण इन्फ-सप स्थिर युग्मों $(U_h, M_h)$ का उपयोग करता है जो विविक्त LBB स्थिति को संतुष्ट करते हैं। समय विवक्तीकरण $q=0$ या $q=1$ के साथ dG($q$) विधि का उपयोग करता है, जो समय नोड्स पर विच्छिन्नता (discontinuities) की अनुमति देता है।
3. स्थिरता और अस्तित्व: एक प्रमुख तकनीकी घटक एक विशेष रूप से तैयार किया गया विविक्त ग्रोनवाल लेम्मा (Lemma 4.11) का विकास है। यह लेम्मा $L^1$ में समय के आधार पर केवल राइट-हैंड साइड को संभालने के लिए अनुकूलित है और डेटा के सापेक्ष घातांकीय वृद्धि (exponential growth) से बचते हुए, अंतिम समय $T$ के संबंध में रैखिक डेटा वाले स्थिरता बाउंड प्राप्त करने की अनुमति देता है। इसका उपयोग विविक्त समाधानों के अस्तित्व और विशिष्टता को सिद्ध करने और मजबूत नॉर्म्स ($L^\infty(I; H^1(\Omega))$ और $L^2(I; H^2(\Omega))$) में उनकी स्थिरता को सिद्ध करने के लिए किया जाता है।
4. त्रुटि विभाजन और द्वैत (Duality): त्रुटि $e = u - u_{kh}$ को दो भागों में विभाजित किया गया है: $\xi = u - \tilde{u}_{kh}$ (एक इंस्टेशनरी स्टोक्स प्रोजेक्शन की त्रुटि) और $\eta_{kh} = \tilde{u}_{kh} - u_{kh}$ (शेष भाग)।
   • स्टोक्स भाग $\xi$ का अनुमान मौजूदा स्टोक्स समस्या के बेस्ट एप्रोक्सिमेशन परिणामों का उपयोग करके लगाया जाता है।
   • शेष भाग $\eta_{kh}$ का अनुमान एक द्वैत तर्क का उपयोग करके लगाया जाता है। एक विविक्त द्वैत समस्या को निरंतर और विविक्त समाधानों के औसत ($u_{ukh} = \frac{1}{2}(u + u_{kh})$) के आसपास रैखिकीकृत (linearized) किया गया है। इस द्वैत समस्या की स्थिरता को विविक्त ग्रोनवाल लेम्मा और पूर्व-सिद्ध प्राइमल समाधानों की रेगुलैरिटी का उपयोग करके स्थापित किया गया है।

प्रमुख योगदान और परिणाम
शोध पत्र ने नेवियर-स्टोक्स समीकरणों के लिए "बेस्ट एप्रोक्सिमेशन टाइप" त्रुटि अनुमान प्राप्त किए हैं, जो $H^1$ त्रुटि के साथ जुड़ने के बिना $L^\infty(I; L^2(\Omega))$ नॉर्म में त्रुटि को अलग करता है।

मुख्य सैद्धांतिक परिणाम हैं:

• $L^2(I; H^1(\Omega))$ अनुमान: $f \in L^2(I; L^2(\Omega)^2)$ और $u_0 \in V$ के लिए, त्रुटि संतुष्ट करती है $\|u - u_{kh}\|_{L^2(I; H^1(\Omega))} \leq C(\sqrt{k} + h)$।
• $L^\infty(I; L^2(\Omega))$ अनुमान: $f \in L^\infty(I; L^2(\Omega)^2)$ और $u_0 \in V \cap H^2(\Omega)^2$ के लिए, त्रुटि संतुष्ट करती है:
  $$\|u - u_{kh}\|_{L^\infty(I; L^2(\Omega))} \leq C \ln(T/k)^2 (k + h^2)$$
  यह परिणाम एक बेस्ट एप्रोक्सिमेशन प्रकार के अनुमान के रूप में प्रस्तुत किया गया है जिसमें विविक्त स्थानों पर इन्फिमम (infimum) और स्थिर स्टोक्स प्रोजेक्शन त्रुटि शामिल है।
• $L^2(I; L^2(\Omega))$ अनुमान: $f \in L^2(I; L^2(\Omega)^2)$ और $u_0 \in V$ के लिए, त्रुटि संतुष्ट करती है $\|u - u_{kh}\|_{L^2(I; L^2(\Omega))} \leq C(k + h^2)$।

ये अनुमान अभिसरण क्रम (लॉग कारकों तक) और आवश्यक रेगुलैरिटी के संदर्भ में इष्टतम (optimal) हैं, जो पिछले परिणामों में सुधार करते हैं जो या तो संयोजित नॉर्म्स (जिससे ऑर्डर रिडक्शन होता है) का उपयोग करते थे या डेटा पर मजबूत रेगुलैरिटी धारणाओं की आवश्यकता रखते थे।

महत्व
लेखक बताते हैं कि इस कार्य का प्राथमिक उद्देश्य राज्य बाधाओं (state constraints) के अधीन अनुकूलन नियंत्रण समस्याओं (optimal control problems) का विश्लेषण करना है जो समय में बिंदुवार (pointwise in time) होती हैं। ऐसी समस्याओं में, संगत लैग्रेंज मल्टीप्लायर समय में एक माप (measure) होता है, जिसके लिए $L^\infty(I; L^2(\Omega))$ नॉर्म में त्रुटि अनुमानों की आवश्यकता होती है। प्राप्त अनुमान dG समय विवक्तीकरण का उपयोग करके इन नियंत्रण समस्याओं के विश्लेषण के लिए आवश्यक सैद्धांतिक आधार प्रदान करते हैं।

यह कार्य यह भी रेखांकित करता है कि विकसित की गई तकनीकें, विशेष रूप से अनुकूलित विविक्त ग्रोनवाल लेम्मा और गैररेखीय पद के लिए द्वैत तर्क, इष्टतम अभिसरण दर और न्यूनतम रेगुलैरिटी आवश्यकताओं को बनाए रखते हुए स्टोक्स त्रुटि अनुमानों को गैररेखीय नेवियर-स्टोक्स समीकरणों तक विस्तारित करने की अनुमति देती हैं। शोध पत्र इस बात पर जोर देता है कि ये परिणाम अक्सर पिछले साहित्य में आवश्यक स्रोत पद $f$ की उच्च रेगुलैरिटी को मान बिना प्राप्त किए गए हैं।