Numerical analysis of the Biot equations coupled to frictional contact mechanics

तकनीकी सारांश: घर्षण संपर्क यांत्रिकी से जुड़े बायोट समीकरणों का संख्यात्मक विश्लेषण

समस्या का निरूपण
यह शोध पत्र एक तरल पदार्थ से संतृप्त पोरो-विस्को-इलास्टिक माध्यम के गणितीय मॉडलिंग और संख्यात्मक सन्निकटन को संबोधित करता है, जो एक कठोर बाधा के साथ घर्षण संपर्क (frictional contact) के अधीन है। भौतिक प्रणाली अर्ध-स्थैतिक (quasi-static) बायोट समीकरणों द्वारा शासित होती है, जो रैखिक लोच (विशेष रूप से केल्विन-वोइग्ट विस्कोइलास्टिक मॉडल) को डार्सी प्रवाह (Darcy flow) के साथ जोड़ती है। संपर्क यांत्रिकी को सीमा के एक विशिष्ट उपसमुच्चय ($\Gamma_c$) पर दो विशिष्ट बाधाओं के माध्यम से मॉडल किया गया है:

1. नॉर्मल कंप्लायंस (Normal Compliance): एक संरचनात्मक संबंध जहाँ सामान्य संपर्क तनाव (normal contact traction), प्रवेश गहराई (penetration depth) के समानुपाती होता है, जो गैर-सुचारू साइनोर्नी (Signorini) स्थिति से बचता है।
2. कूलम्ब घर्षण (Coulomb Friction): एक नियम जो स्पर्शरेखीय गति (tangential motion) को नियंत्रित करता है, जहाँ शरीर तब तक चिपका रहता है जब तक कि स्पर्शरेखीय तनाव, घर्षण सीमा (सामान्य तनाव के समानुपाती) से नीचे रहता है, और अन्यथा फिसलता है।

परिणामी गणितीय मॉडल एक रैखिक आंशिक अंतर समीकरण (द्रव्यमान संतुलन समीकरण) के साथ एक गैर-रैखिक परिवर्तनशील असमानता (मोमेंटम संतुलन समीकरण जिसमें संपर्क स्थितियाँ शामिल हैं) के युग्मन का एक तंत्र है। प्राथमिक चुनौती प्रवाह और यांत्रिकी के बीच का युग्मन, समस्या की समय-निर्भरता, और घर्षण संपर्क स्थितियों द्वारा उत्पन्न गैर-रैखिकता है।

कार्यप्रणाली
लेखक इस समस्या के लिए एक पूर्णतः विविक्त (fully discrete) संख्यात्मक योजना प्रस्तावित करते हैं और उसका विश्लेषण करते हैं:

• स्थानिक विविक्तीकरण (Spatial Discretization): अनुरूप परिमित तत्वों (conformal finite elements) का उपयोग किया जाता है। विस्थापन स्थान का अनुमान क्रम $k_1$ के लैग्रेंजियन तत्वों का उपयोग करके लगाया जाता है, और दबाव स्थान का अनुमान क्रम $k_2$ के तत्वों का उपयोग करके लगाया जाता है।
• कालिक विविक्तीकरण (Temporal Discretization): समय-चरण (time stepping) के लिए अंतर्निहित यूलर विधि (implicit Euler method) का उपयोग किया जाता है।
• अस्तित्व प्रमाण रणनीति (Existence Proof Strategy): विविक्त समाधान के अस्तित्व और अद्वितीयता को स्थापित करने के लिए, लेखक एक नवीन सैद्धांतिक दृष्टिकोण अपनाते हैं। वे पहले एक सरलीकृत समस्या पर विचार करते हैं जहाँ संपर्क पदों को "लैग्ड" (lagged) किया जाता है (प्रभावी रूप से घर्षण मॉडल को ट्रेस मॉडल में बदल देता है)। इस सरलीकृत समस्या के लिए, वे प्रवाह और यांत्रिकी को अलग करने के लिए आमतौर पर उपयोग किए जाने वाले 'फिक्स्ड-स्ट्रेस स्प्लिटिंग स्कीम' को एक सैद्धांतिक उपकरण के रूप में उपयोग करते हैं, ताकि कंट्रैक्शन तर्क के माध्यम से समाधान के अस्तित्व को सिद्ध किया जा सके। इसके बाद, पूर्ण कूलम्ब घर्षण मॉडल को पुनः प्राप्त करने के लिए संपर्क यांत्रिकी साहित्य से एक फिक्स्ड-पॉइंट तर्क लागू किया जाता है।
• त्रुटि विश्लेषण (Error Analysis): पूर्ववर्ती त्रुटि अनुमानों (a priori error estimates) का व्युत्पन्न तकनीक बायोट समीकरणों (गिरौट एट अल., फिलिप्स और व्हीलर) और समय-निर्भर परिवर्तनशील असमानताओं (सोफोना और हान) से तकनीकों को जोड़ता है। एक प्रमुख विश्लेषणात्मक चरण यह दिखाना है कि प्रवाह और यांत्रिकी के बीच का युग्मन पद स्थिरता और त्रुटि विश्लेषण में रद्द हो जाता है, जिससे यह उनके अनुमानों को खराब होने से रोकता है।

प्रमुख परिणाम

1. सुव्यवस्थितता (Well-posedness): शोध पत्र विविक्त समाधान के अस्तित्व और अद्वितीयता को स्थापित करता है, बशर्ते कि समय-चरण का आकार द्विपक्षीय रूपों (bilinear forms) के कोर्सिविटी स्थिरांकों और संपर्क फलन के लिप्सचिट्ज़ स्थिरांक से संबंधित एक विशिष्ट शर्त को पूरा करता हो।
2. स्थिरता (Stability): एक स्थिरता अनुमान निकाला गया है, जो यह दर्शाता है कि विविक्त समाधान मेष आकार ($h$) और समय चरण ($\Delta t$) से स्वतंत्र रूप से सीमित रहता है, यह मानते हुए कि समय चरण संपर्क पदों की गैर-रैखिकता को संभालने के लिए पर्याप्त छोटा है।
3. पूर्ववर्ती त्रुटि अनुमान (A Priori Error Estimates): लेखक विस्थापन के लिए $H^1$-नॉर्म और दबाव के लिए $L^2$-नॉर्म में त्रुटि अनुमान प्राप्त करते हैं। सैद्धांतिक अभिसरण दरें हैं:
   • समय: प्रथम-क्रम अभिसरण ($O(\Delta \text{t})$)।
   • स्थान: दबाव अनुकूल क्रम ($O(h^s)$) के साथ अभिसरण करता है। विस्थापन अभिसरण का क्रम मानक लोच समस्याओं की तुलना में $\sqrt{h}$ के कारक से कम हो जाता है, जिसके परिणामस्वरूप रैखिक तत्वों के लिए लगभग $O(h^{r-0.5})$ का क्रम प्राप्त होता है। यह कमी संपर्क स्थितियों के कारण समाधान की नियमितता (regularity) कम होने के कारण है।
4. संख्यात्मक प्रयोग: FEniCSx सॉफ़्टवेयर का उपयोग करके एक संख्यात्मक प्रयोग किया गया। परिणामों ने सैद्धांतिक भविष्यवाणियों की पुष्टि की:
   • दबाव और विस्थापन दोनों के लिए प्रथम-क्रम का समय अभिसरण देखा गया।
   • दबाव के लिए स्थानिक अभिसरण इष्टतम था।
   • विस्थापन के लिए स्थानिक अभिसरण थोड़ा कम हुआ (रैखिक तत्वों के लिए लगभग 0.8 की दर देखी गई), जो संपर्क के कारण कम नियमितता की सैद्धांतिक अपेक्षा के अनुरूप है।

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

लेखक जोर देते हैं कि उनका कार्य विस्कोइलास्टिक संपर्क समस्याओं (जो गैर-पोरस मीडिया के लिए पहले से स्थापित हैं) के सुलभ विश्लेषण को युग्मित पोरो-विस्को-इलास्टिक सेटिंग तक विस्तारित करता है। इस कार्य का महत्व है:

• युग्मित प्रवाह और घर्षण संपर्क वाले छिद्रपूर्ण माध्यमों के संख्यात्मक सिमुलेशन के लिए एक कठोर गणितीय आधार प्रदान करना।
• यह प्रदर्शित करना कि पोरोइलास्टिसिटी और संपर्क यांत्रिकी की तकनीकों को संयोजित करते समय युग्मन पद त्रुटि अनुमानों के व्युत्पन्न में बाधा नहीं डालता है।
• एक पूर्णतः विविक्त योजना प्रदान करना जिसकी स्थिरता और अभिसरण दरें प्रमाणित हैं और संख्यात्मक प्रयोगों द्वारा सत्यापित हैं।

यह शोध पत्र अपने अनुप्रयोगों के संबंध में विनम्र बना हुआ है, जो व्युत्पन्न अनुमानों के दायरे से परे नए भौतिक अनुप्रयोगों या भविष्य के विस्तार के बजाय पूरी तरह से प्रस्तावित संख्यात्मक योजना के गणितीय विश्लेषण और सत्यापन पर केंद्रित है।