Embedded Trefftz DG method for steady Navier-Stokes flow. Part I: Oseen linearization

तकनीकी सारांश: स्थिर नेवियर-स्टोक्स प्रवाह के लिए एम्बेडेड ट्रेफ़्ट्ज़ DG विधि (भाग I: ओसीन लीनियराइजेशन)

समस्या विवरण
यह शोध पत्र स्थिर ओसीन समस्या (steady Oseen problem) के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है, जो एक निर्धारित संवहन क्षेत्र (convection field) $w$ के चारों ओर स्थिर असंपीड्य नेवियर-स्टोक्स समीकरणों के लीनियराइजेशन (linearization) के रूप में कार्य करता है। नियंत्रित समीकरण एक वेग $u$ और दबाव $p$ की खोज करते हैं जो एक परिबद्ध डोमेन $\Omega \subset \mathbb{R}^d$ ($d=2,3$) में समरूप डिरिचलेट सीमा स्थितियों (homogeneous Dirichlet boundary conditions) को संतुष्ट करते हैं:
$$-\nu\Delta u + (w \cdot \nabla)u + \nabla p = f, \quad \text{div } u = 0$$
मुख्य चुनौती एक एम्बेडेड ट्रेफ़्ट्ज़ डिस्कंटीन्यूअस गैलेरकिन (DG) विधि विकसित करना है जो कम डिग्री ऑफ फ्रीडम (DoFs) के साथ उच्च सन्निकटन गुणवत्ता बनाए रखती है और संवहन पद (convection term) $(w \cdot \nabla)u$ द्वारा पेश किए गए परिवर्तनीय गुणांकों (variable coefficients) को संभाल सकती है। स्टोक्स समस्या (जहाँ $w=0$) के विपरीत, ओसीन ऑपरेटर में परिवर्तनीय गुणांक शामिल होते हैं, जिससे मानक स्थिर-गुणांक ट्रेफ़्ट्ज़ निर्माण अपर्याप्त हो जाते हैं।

कार्यप्रणाली (Methodology)
लेखक एक एम्बेडेड ट्रेफ़्ट्ज़-DG दृष्टिकोण प्रस्तावित करते हैं जो स्थानीय ट्रेफ़्ट्ज़ आधारों (local Trefftz bases) के स्पष्ट निर्माण से बचता है। इसके बजाय, विधि एक मानक बहुपद DG स्थान $X$ से शुरू होती है और स्थानीय प्रक्षेपण बाधाओं (projection constraints) के माध्यम से PDE को लागू करती है।

1. स्थानीय ऑपरेटर और ट्रेफ़्ट्ज़ स्थान:
   एक स्थानीय ऑपरेटर $K_{T,w}$ प्रत्येक तत्व $T$ पर परिभाषित किया जाता है, जो समाधान स्थान को एक द्वैत परीक्षण स्थान (dual test space) $Q(T)'$ पर मैप करता है। ट्रेफ़्ट्ज़ स्थान $T_w$ को $X$ के उस उप-स्थान के रूप में परिभाषित किया गया है जहाँ स्थानीय अवशेष (residual) $Q(T)$ के विरुद्ध एक दुर्बल अर्थ (weak sense) में शून्य होता है। विशेष रूप से, $T_w$ के फलन केवल तब PDE को संतुष्ट करते हैं जब उन्हें निम्न-क्रम के बहुपदों (वेग के लिए डिग्री $k-2$, दबाव के लिए $k-1$) के विरुद्ध परीक्षण किया जाता है।

2. स्थानीय-वैश्विक विभाजन और पूरक स्थान (Local-Global Splitting and Complement Spaces):
   एक महत्वपूर्ण सैद्धांतिक चरण एक स्थानीय पूरक स्थान $L(T)$ का निर्माण है ताकि $X(T) = T_w(T) \oplus L(T)$ हो सके।

• वेग घटक ($L_u$): इसे एक सब-बॉल पर शून्य होने वाले बबल फलनों (bubble functions) का उपयोग करके बनाया गया है, जिसे विशिष्ट बहुपद पदों के साथ समृद्ध किया गया है ताकि फैसेट फ्लक्स (facet fluxes) को नियंत्रित किया जा सके, और $d$ फैसेट्स पर शून्य माध्य फ्लक्स सुनिश्चित करने के लिए प्रतिबंधित किया गया है।
• दबाव घटक ($L_p$): इसे डिग्री $k-1$ के शून्य-माध्य बहुपदों के रूप में परिभाषित किया गया है।
• व्युत्क्रमणीयता (Invertibility): शोध पत्र यह सिद्ध करता है कि स्थानीय ऑपरेटर $K_{T,w}$, $L(T)$ पर बाउंडेडली व्युत्क्रमणीय (boundedly invertible) है, बशर्ते कि संवहन क्षेत्र $w$ एक समाधान स्थिति (resolution condition) को संतुष्ट करे (विशेष रूप से श्यानता (viscosity) के सापेक्ष $|w|_{h,d}$ पर एक लघुता की स्थिति)। यह व्युत्क्रमणीयता यह सुनिश्चित करती है कि समाधान के स्थानीय भाग को प्रत्येक तत्व पर स्वतंत्र रूप से संगणना किया जा सके।

3. न्यूनीकृत सूत्रीकरण (Reduced Formulation):
   उच्च-क्रम दबाव घटकों को स्थानीय स्तर पर हटाने के बाद (जो स्थानीय ऑपरेटर के माध्यम से वेग से जुड़े होते हैं), लेखक एक न्यूनीकृत सूत्रीकरण प्राप्त करते हैं जो केवल वेग अज्ञात (velocity unknown) के संदर्भ में प्रस्तुत किया जाता है। यह सूत्रीकरण:

• स्थानीय स्तर पर सैडल-पॉइंट संरचना (saddle-point structure) को हटा देता है।
• स्टोक सेटिंग से कोएर्सिविटी (coercivity) को रिड्यूस्ड ओसीन-ट्रेफ़्ट्ज़ स्थानों में स्थानांतरित करता है।
• केवल पीसवाइस-कांस्टेंट दबाव घटक के लिए वैश्विक युग्मन (global coupling) बनाए रखता है ताकि असंपीड्य डाइवर्जेंस बाधा को लागू किया जा सके।

4. स्थिरता विश्लेषण (Stability Analysis):
   विश्लेषण निम्नलिखित को स्थापित करता है:

• द्विक-रूपों (bilinear forms) की निरंतरता (Continuity)।
• रिड्यूस्ड फॉर्म का ट्रेफ़्ट्ज़ उप-स्थान पर कोएर्सिविटी (Coercivity)।
• मानक DG रूपों में पूर्ण युग्मित प्रणाली के लिए इन्फ-सुप स्थिरता (Inf-sup stability)।
  ये परिणाम अमूर्त स्थानीय-वैश्विक ढांचे और $L(T)$ के विशिष्ट निर्माण पर निर्भर करते हैं, जो सु-संरचित (well-posedness) और क्वासी-ऑप्टिमलिटी (quasi-optimality) सुनिश्चित करते हैं।

प्रमुख योगदान
यह शोध पत्र स्टोक्स समस्या से आगे बढ़कर स्थिर ओसीन समीकरणों तक एम्बेडेड ट्रेफ़्ट्ज़-DG कार्यप्रणाली का विस्तार करता है, जो बाद के गैर-रेखीय नेवियर-स्टोक्स विश्लेषण (भाग II) के लिए रैखिक आधार प्रदान करता है। विशिष्ट योगदान हैं:

• स्थानीय पूरक स्थानों का निर्माण: स्टोक्स मामले के विपरीत, जो स्थिर-गुणांक टेलर-बहुपद तर्कों पर निर्भर था, यह कार्य संवहन पद से उत्पन्न परिवर्तनीय-गुणांक स्थानीय ऑपरेटर को संभालने के लिए $L(T)$ का निर्माण करता है।
• विषम दाएँ-पक्ष (Inhomogeneous Right-Hand Sides) का उपचार: विधि सुधार तर्क (correction argument) के बजाय सीधे स्थानीय समस्या प्रारूप में विषम पदों को शामिल करती है, जो अमूर्त स्थानीय-वैश्विक विभाजन ढांचे में स्वाभाविक रूप से फिट बैठती है।
• न्यूनीकृत सूत्रीकरण और स्थिरता: लेखक एक न्यूनीकृत वेग-मात्र सूत्रीकरण प्राप्त करते हैं और पूर्ण ओसीन ऑपरेटर के लिए न्यूनीकृत इन्फ-सुप स्थिरता सिद्ध करते हैं। यह स्थानीय स्तर पर सैडल-पॉइंट संरचना को समाप्त करता है और संवहन क्षेत्र पर एक उपयुक्त लघुता की स्थिति के तहत कोएर्सिविटी के हस्तांतरण की अनुमति देता है।

परिणाम

• सैद्धांतिक: शोध पत्र मानक DG रूपों में पूर्ण स्थिरता और क्वासी-ऑप्टिमलिटी को सिद्ध करता है। यह स्थापित करता है कि विधि सु-संरचित है और त्रुटि (error) सर्वोत्तम सन्निकटन त्रुटि द्वारा सीमित है, जो कि डिस्क्रीट स्पेस में स्थिरांकों द्वारा स्केल की गई है।
• संख्यात्मक:
  • ओसीन-कोवासज़ना (Oseen–Kovasznay) परीक्षण: निर्धारित ओसीन-कोवासज़ना प्रवाह पर संख्यात्मक प्रयोगों से पता चलता है कि ट्रेफ़्ट्ज़-DG विधि अपेक्षित एसिम्प्टोटिक अभिसरण दर ($u$ के लिए $O(h^{k+1})$, $p$ के लिए $O(h^k)$) प्राप्त करती है।
  • दक्षता: जबकि मानक DG विधि निश्चित बहुपद डिग्री के लिए थोड़ा छोटा त्रुटि स्थिरांक देती है, ट्रेफ़्ट्ज़-DG विधि उच्च बहुपद क्रमों के लिए, विशेष रूप से, काफी कम डिग्री ऑफ फ्रीडम के साथ वांछित त्रुटि स्तर प्राप्त करती है।
  • संगणनात्मक लागत: टाइमिंग परिणाम संकेत देते हैं कि ट्रेफ़्ट्ज़-DG विधि लगातार मानक DG विधि से तेज़ है, और महीन मेष (finer meshes) तथा उच्च क्रमों के लिए इसका लाभ बढ़ता जाता है।
  • सीमाएँ: विधि संवहन-प्रधान शासन (convection-dominated regime - कम श्यानता) में संवेदनशीलता दिखाती है, जहाँ त्रुटियाँ मानक DG की तुलना में जल्दी खराब हो जाती हैं। यह स्थानीय एम्बेडेड ट्रेफ़्ट्ज़ निर्माण के रिज़ॉल्यूशन और कंडीशनिंग के कारण है।

महत्व
यह शोध पत्र दावा करता है कि यह भाग II में विकसित गैर-रेखीय नेवियर-स्टोक्स सिद्धांत के लिए आवश्यक रैखिक विश्लेषणात्मक आधार प्रदान करता है। परिवर्तनीय गुणांकों वाले ओसीन समस्या तक एम्बेडेड ट्रेफ़्ट्ज़ दृष्टिकोण को सफलतापूर्वक विस्तारित करके, लेखक प्रदर्शित करते हैं कि ट्रेफ़्ट्ज़ तकनीकों को स्पष्ट आधार निर्माण के बिना विषम, परिवर्तनीय-गुणांक ऑपरेटरों पर लागू किया जा सकता है। न्यूनीकृत सूत्रीकरण को गैर-रेखीय समस्या के विश्लेषण के लिए एक महत्वपूर्ण कदम के रूप में रेखांकित किया गया है, क्योंकि यह स्थानीय दबाव को वैश्विक प्रणाली से अलग करता है, जिससे इटरेटिव सॉल्वर के लिए स्थिरता विश्लेषण सरल हो जाता है। यह कार्य ओसीन समस्या को शास्त्रीय DG विविक्तकरणों (discretizations) और ट्रेफ़्ट्ज़-प्रकार के मॉडल रिडक्शन के बीच एक सेतु के रूप में स्थापित करता है, जो श्यान असंपीड्य प्रवाह के लिए उच्च-क्रम, निम्न-DoF सिमुलेशन का मार्ग प्रशस्त करता है।