Quadrature error estimates on non-matching grids in a fictitious domain framework for fluid-structure interaction problems

तकनीकी सारांश: एक काल्पनिक डोमेन फ्रेमवर्क में गैर-मिलान वाले ग्रिड्स पर क्वाड्रैचर एरर अनुमान (Fluid–Structure Interaction समस्याओं के लिए)

समस्या विवरण
यह शोध पत्र डिस्ट्रिब्यूटेड लैग्रेंज मल्टीप्लायर्स (distributed Lagrange multipliers) पर आधारित एक काल्पनिक डोमेन फॉर्मुलेशन का उपयोग करके फ्लूइड-स्ट्रक्चर इंटरेक्शन (FSI) समस्याओं के संख्यात्मक विविक्तकरण (numerical discretization) को संबोधित करता है। इस फ्रेमवर्क में, फ्लूइड और सॉलिड को दो पूरी तरह से स्वतंत्र, गैर-मिलान वाले मेषों (non-matching meshes) का उपयोग करके विविक्तकृत किया जाता है। सॉलिड को एक संदर्भ डोमेन $B$ पर लैग्रेंजियन फ्रेमवर्क में वर्णित किया गया है, जिसे एक स्थिर यूलेरियन फ्लूइड डोमेन $\Omega$ के भीतर उसकी वर्तमान स्थिति में मैप किया गया है। फ्लूइड वेलोसिटी और सॉलिड विरूपण (deformation) के बीच गतिज युग्मन (kinematic coupling) को सॉलिड संदर्भ डोमेन पर परिभाषित लैग्रेंज मल्टीप्लायर $\lambda$ के माध्यम से लागू किया जाता है।

एक महत्वपूर्ण कम्प्यूटेशनल चुनौती कपलिंग (इंटरफेस) मैट्रिक्स को असेंबल करने में आती है। इस मैट्रिक्स के लिए फ्लूइड बेसिस फंक्शन्स (जो सॉलिड मैपिंग के साथ संयोजित हैं) और सॉलिड बेसिस फंक्शन्स के गुणनफल को सॉलिड एलिमेंट्स पर एकीकृत करने की आवश्यकता होती है। चूंकि ये फलन गैर-मिलान वाले ग्रिड्स पर परिभाषित हैं, इसलिए एक कंपोजिट क्वाड्रैचर नियम बनाने के लिए फ्लूइड और मैप किए गए सॉलिड मेष के ज्यामितीय प्रतिच्छेदन (geometric intersection) की गणना करना सटीक एकीकरण (exact integration) के लिए आवश्यक है। यह प्रक्रिया कम्प्यूटेशनल रूप से महंगी है, विशेष रूप से क्योंकि कपलिंग टर्म को प्रत्येक टाइम स्टेप पर असेंबल किया जाना चाहिए। वैकल्पिक रूप से, एक अनुमानित दृष्टिकोण बिना मेष प्रतिच्छेदन के, मानक क्वाड्रैचर नियमों का उपयोग करके सीधे सॉलिड एलिमेंट्स पर एकीकरण करता है, जिससे एक क्वाड्रैचर एरर (quadrature error) उत्पन्न होता है। यह शोध पत्र इस अनुमानित एकीकरण रणनीति से जुड़े सैद्धांतिक वैधता और एरर अनुमानों की जांच करता है।

कार्यप्रणाली (Methodology)
लेखक एक स्टेडी-स्टेट स्टोक्स-आधारित FSI मॉडल का विश्लेषण करते हैं, हालांकि वे उल्लेख करते हैं कि उनके परिणाम अधिक सामान्य गैररेखीय मामलों के प्रतिनिधि हैं। विवरणात्मक फॉर्मुलेशन (variational formulation) में फ्लूइड वेलोसिटी $u$, प्रेशर $p$, सॉलिड विरूपण मैप $X$, और लैग्रेंज मल्टीप्लायर $\lambda$ खोजने की प्रक्रिया शामिल है जो एक सैडल-पॉइंट सिस्टम (saddle-point system) को संतुष्ट करते हैं।

अध्ययन कपलिंग को परिभाषित करने के लिए उपयोग किए जाने वाले द्विरेखीय रूप (bilinear form) $c(\cdot, \cdot)$ के दो विशिष्ट विकल्पों पर केंद्रित है:

1. $c_0$ (Dual pairing): $L^2(B)$ इनर प्रोडक्ट पर आधारित, जहाँ मल्टीप्लायर स्पेस $H^1(B)$ का ड्यूल है।
2. $c_1$ (Scalar product): $H^1(B)$ इनर प्रोडक्ट पर आधारित, जहाँ मल्टीप्लायर स्पेस स्वयं $H^1(B)$ है।

कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

• विविक्तिकरण (Discretization): सभी चरों के लिए परिमित तत्व स्थानों (finite element spaces) को पेश किया जाता है। फ्लूइड के लिए बेर्वियर-पिरोनो एलिमेंट ($P1$-iso-$P2/P1$) का उपयोग किया जाता है, जबकि सॉलिड मैपिंग और मल्टीप्लायर निरंतर पीसवाइज लीनियर एलिमेंट्स का उपयोग करते हैं।
• अनुमानित एकीकरण (Approximate Integration): लेखक उन डिस्क्रीट द्विरेखीय रूपों को परिभाषित करते हैं जहाँ कपलिंग टर्म्स को सॉलिड मेष पर मानक क्वाड्रैचर नियमों का उपयोग करके कंप्यूट किया जाता है, यह नजरअंदाज करते हुए कि संयोजित फ्लूइड बेसिस फंक्शन्स सॉलिड मेष के सापेक्ष विच्छिन्न (discontinuous) पीसवाइज पॉलीनोमियल्स हैं।
• अमूर्त त्रुटि विश्लेषण (Abstract Error Analysis): इन्फ-सफ (inf-sup) स्थितियों के ढांचे का पालन करते हुए, लेखक एक सामान्य त्रुटि अनुमान (थ्योरम 2) व्युत्पन्न करते हैं, जो सटीक डिस्क्रीट समाधान और अपूर्ण एकीकरण के साथ प्राप्त समाधान के बीच के अंतर को बांधता है। यह सीमा क्वाड्रैचर एरर फंक्शनल्स पर निर्भर करती है।
• क्वाड्रैचर एरर अनुमान (Quadrature Error Estimation): मुख्य तकनीकी योगदान कपलिंग टर्म्स के लिए विशिष्ट अनुमानों को सिद्ध करना है। लेखक इस तथ्य का उपयोग करते हैं कि संयोजित फलन $v_h(X)$ सॉलिड एलिमेंट्स पर निरंतर लेकिन पीसवाइज लीनियर है, जिससे भिन्नात्मक सोबोलेव स्पेस ($H^{1+s}$) में विशिष्ट रेगुलैरिटी गुण प्राप्त होते हैं। वे एरर को बाउंड करने के लिए इंटरपोलेशन थ्योरी और इनवर्स इनइक्वेलिटीज (inverse inequalities) का उपयोग करते हैं।
• स्थिरता प्रमाण (Stability Proof): लेखक यह सिद्ध करते हैं कि डिस्क्रीट समस्या तब भी सुव्यवस्थित (well-posed) रहती है (इन्फ-सफ स्थितियों को संतुष्ट करती है), जब तक कि विशिष्ट मेष साइज की शर्तें पूरी हों, भले ही कपलिंग टर्म का निर्माण अनुमानित रूप से किया गया हो।

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

1. अनुमानित कपलिंग के साथ सुव्यवस्थितता (Well-Posedness): पेपर यह सिद्ध करता है कि अनुमानित क्वाड्रैचर का उपयोग करके कपलिंग मैट्रिक्स को असेंबल करने पर भी डिस्क्रीट समस्या स्थिर और सुव्यवस्थित रहती है, जिसके लिए मेष प्रतिच्छेदन की आवश्यकता नहीं है।
2. क्वाड्रैचर एरर अनुमान:
   • $L^2$-आधारित कपलिंग ($c_0$) के लिए, क्वाड्रैचर एरर $O(h_B^{3/2} |\log h_B|)$ द्वारा सीमित है। लेखक दिखाते हैं कि अनुमानित विधि इष्टतम अभिसरण दर (optimal convergence rates) प्रदान करती है, बशर्ते $h_B \to 0$ हो।
   • $H^1$-आधारित कपलिंग ($c_1$) के लिए, ग्रेडिएंट टर्म के कारण एरर विश्लेषण अधिक जटिल है। एरर बाउंड में मेष साइज के अनुपात पर निर्भर टर्म शामिल हैं। विशेष रूप से, इष्टतम अभिसरण के लिए न केवल $h_B \to 0$ की आवश्यकता है, बल्कि यह भी कि $h_B/h_\Omega \to 0$ हो (जहाँ $h_\Omega$ फ्लूइड मेष साइज है)। यदि $h_B/h_\Omega$ को स्थिर रखा जाता है, तो विधि अपनी इष्टतमता खो सकती है।
3. इनवर्स इनइक्वेलिटीज (Inverse Inequalities): लेखक ड्यूल नॉर्म के मामले में लैग्रेंज मल्टीप्लायर के नॉर्म को संभालने के लिए विशिष्ट इनवर्स इनइक्वेलिटीज व्युत्पन्न करते हैं और उनका उपयोग करते हैं।
4. संख्यात्मक सत्यापन (Numerical Validation): सैद्धांतिक अनुमानों को मान्य करने के लिए चार संख्यात्मक परीक्षण प्रस्तुत किए गए हैं:
   • टेस्ट 1 और 2: विभिन्न मेष रिफाइनमेंट रणनीतियों (कॉन्स्टेंट रेश्यो बनाम डिक्रीजिंग रेश्यो) के तहत $c_0$ और $c_1$ के लिए अभिसरण दरों की तुलना करते हैं। परिणाम पुष्टि करते हैं कि $c_0$ अनुमानित एकीकरण के साथ मजबूत (robust) है, जबकि $c_1$ को सटीकता बनाए रखने के लिए $h_B/h_\Omega \to 0$ की आवश्यकता होती है।
   • टेस्ट 3: इंटरफेस पर वेलोसिटी ग्रेडिएंट में जंप (जो FSI में सामान्य है) वाले समाधान का उपयोग करता है, जो अपेक्षित सब-ऑप्टिमल अभिसरण दर $O(h^{1/2})$ की पुष्टि करता है।
   • टेस्ट 4: उच्च-क्रम तत्वों ($P2/P1$) की जांच करता है, जो दिखाता है कि जबकि $c_0$ मजबूत बना रहता है, $c_1$ को इष्टतम दर प्राप्त करने के लिए अभी भी सटीक एकीकरण या विशिष्ट मेष रेश्यो की आवश्यकता होती है।

महत्व और दावे
पेपर का दावा है कि यह काल्पनिक डोमेन FSI विधियों में अनुमानित एकीकरण का उपयोग करने के लिए एक कठोर सैद्धांतिक औचित्य प्रदान करता है, जो अक्सर कम्प्यूटेशनल दक्षता के लिए अपनाया जाता है। लेखक प्रदर्शित करते हैं कि:

• "अपूर्ण" दृष्टिकोण संख्यात्मक योजना की स्थिरता से समझौता नहीं करता है।
• $L^2$-आधारित कपलिंग ($c_0$) अनुमानित एकीकरण के लिए विशेष रूप से उपयुक्त है, क्योंकि इसके लिए सटीकता बनाए रखने के लिए फ्लूइड और सॉलिड मेष को एक विशिष्ट अनुपात में रिफाइन करने की आवश्यकता नहीं होती है।
• $H^1$-आधारित कपलिंग ($c_1$), हालांकि अलग स्थिरता गुण प्रदान कर सकती है, लेकिन इसमें एरर को समाधान एरर पर हावी होने से रोकने के लिए सख्त मेष आवश्यकताएं ($h_B/h_\Omega \to 0$) शामिल हैं।

लेखक निष्कर्ष निकालते हैं कि उनके सैद्धांतिक परिणाम उनके पिछले कार्य [12, 13] के प्रयोगात्मक अवलोकनों की पुष्टि करते हैं, जो विशिष्ट कपलिंग फॉर्म और मेष कॉन्फ़िगरेशन के आधार पर सटीक (मेष-इंटरसेक्टिंग) बनाम अनुमानित (नॉन-इंटरसेक्टिंग) असेंबली तकनीकों के बीच चयन करने के लिए एक ठोस गणितीय आधार प्रदान करते हैं।