A scalable preconditioning framework for stabilized contact mechanics with hydraulically active fractures

तकनीकी सारांश: हाइड्रोलिक रूप से सक्रिय फ्रैक्चर के साथ स्थिर संपर्क यांत्रिकी के लिए एक स्केलेबल प्रीकंडीशनिंग फ्रेमवर्क

समस्या विवरण
यह शोध पत्र छिद्रयुक्त माध्यम (porous medium) में अंतर्निहित फ्रैक्चर नेटवर्क के भीतर घर्षण संपर्क यांत्रिकी (frictional contact mechanics) और द्रव प्रवाह (fluid flow) के युग्मित भौतिकी (coupled physics) के अनुकरण की कम्प्यूटेशनल चुनौती को संबोधित करता है। यह युग्मित समस्या भू-तापीय ऊर्जा उत्पादन, हाइड्रोलिक फ्रैक्चरिंग और CO2 अनुक्रमण (sequestration) जैसे उपसतह अनुप्रयोगों के लिए महत्वपूर्ण है। भौतिक प्रक्रियाएं आपस में गहराई से जुड़ी हुई हैं: फ्रैक्चर का एपर्चर (aperture) और फिसलन (slippage) द्रव प्रवाह को संचालित करते हैं, जबकि द्रव दबाव के परिवर्तन स्ट्रेस अवस्था को प्रभावित करते हैं और संपर्क यांत्रिकी को प्रभावित करते हैं।

लेखक एक डिस्क्रीट फ्रैक्चर मॉडल (DFM) का उपयोग करते हैं जहाँ फ्रैक्चर को स्पष्ट रूप से दर्शाया गया है। विविक्तीकरण (discretization) मैकेनिकल विस्थापन क्षेत्र (mechanical displacement field) के लिए लो-ऑर्डर निरंतर परिमित तत्वों (low-order continuous finite elements) और द्रव दबाव के लिए सेल-सेंटर्ड परिमित आयतन योजना (cell-centered finite volume scheme) का उपयोग करता है। संपर्क बाधाओं (contact constraints) को सेल-सेंटर्ड लैग्रेंज मल्टीप्लायर्स (जो ट्रैक्शन बलों का प्रतिनिधित्व करते हैं) का उपयोग करके लागू किया जाता है जो दबाव क्षेत्र के समान ग्रिड साझा करते हैं, जिससे क्षेत्रों के बीच इंटरपोलेशन (interpolation) की आवश्यकता समाप्त हो जाती है। हालांकि, यह विशिष्ट सूत्रीकरण समान रूप से इन्फ-सुप स्थिर (inf-sup stable/LBB stable) नहीं है और इसके लिए एक वैश्विक बीजगणितीय स्थिरीकरण (global algebraic stabilization) की आवश्यकता होती है।

परिणामी रैखिक प्रणाली, जो गैर-रैखिक समस्या के न्यूटन-क्रायलोव लीनियराइजेशन (Newton-Krylov linearization) से उत्पन्न होती है, एक बड़ी, विरल (sparse), गैर-सममित $3 \times 3$ ब्लॉक जैकोबियन मैट्रिक्स है। मैट्रिक्स की संरचना गतिशील है, जो फ्रैक्चर नेटवर्क के स्टिक (stick), स्लिप (slip) और ओपन (open) क्षेत्रों के माध्यम से विकसित होने के दौरान अपना स्वरूप बदलती रहती है (सिमेट्रिक सैडल-पॉइंट से नॉन-सिमेट्रिक इंडेफिनेट तक)। बड़े पैमाने की समस्याओं के लिए मानक डायरेक्ट सॉल्वर कम्प्यूटेशनल रूप से निषेधात्मक हैं, जिससे मजबूत, स्केलेबल इटरेटिव सॉल्वर की आवश्यकता होती है।

कार्यप्रणाली (Methodology)
मुख्य योगदान भौतिक रूप से सूचित ब्लॉक विभाजन (physically-informed block partitioning) और अत्याधुनिक बीजगणितीय मल्टीग्रिड (Algebraic Multigrid - AMG) तकनीकों पर आधारित एक स्केलेबल प्रीकंडीशनिंग फ्रेमवर्क का डिज़ाइन है। फ्रेमवर्क पूरी तरह से युग्मित मल्टी-फिजिक्स सिस्टम को एकल-भौतिकी उप-समस्याओं (single-physics sub-problems) तक सीमित करने, उन्हें अनुमानित रूप से हल करने और फिर समाधान को पूर्ण प्रणाली में प्रोलोंग (prolong) करने के माध्यम से कार्य करता है।

दो अलग-अलग रिडक्शन सीक्वेंस (दृ pendekatan) प्रस्तावित और विश्लेषित किए गए हैं:

1. $t-p-u$ दृष्टिकोण (ट्रैक्शन-प्रेशर-डिस्प्लेसमेंट):

   • रणनीति: विस्थापन ($u$) के लिए एक एकल साम्यावस्था समीकरण (equilibrium equation) तक सिस्टम को कम करने के लिए ट्रैक्शन ($t$) और दबाव ($p$) चरों को एक साथ समाप्त करता है।
   • शूर कॉम्प्लीमेंट (Schur Complement): इसके परिणामस्वरूप एक शूर कॉम्प्लीमेंट $S = A + C_1 \tilde{H}^{-1}C_2 - Q_1 T^{-1}Q_2$ प्राप्त होता है।
   • अनुमान: ट्रैक्शन ब्लॉक $H$ को एक ब्लॉक-डायगोनल मैट्रिक्स $\tilde{H}$ द्वारा और दबाव ब्लॉक $T$ को उसके डायगोनल $\tilde{T}$ द्वारा अनुमानित किया जाता है। परिणामी शूर कॉम्प्लीमेंट $S$ के व्युत्क्रम (inverse) को एक AMG ऑपरेटर (विशेष रूप से PETSc से GAMG) का उपयोग करके अनिश्चित रूप से (inexactly) लागू किया जाता है।
   • विशेषताएं: यह दृष्टिकोण प्रति इटरेशन कम्प्यूटेशनल रूप से सस्ता (एक AMG कॉल) है लेकिन $S$ की पॉजिटिव डेफिनिटनेस (positive definiteness) पर निर्भर करता है। लेख नोट करता है कि $S$ फ्रैक्चर की स्थिति के आधार पर (विशेष रूप से जब ओपन क्षेत्र मौजूद होते हैं) इंडेफिनेट और नॉन-सिमेट्रिक हो सकता है, जिससे AMG सॉल्वर की प्रभावशीलता से समझौता हो सकता है।

2. $t-u-p$ दृष्टिकोण (ट्रैक्शन-डिस्प्लेसमेंट-प्रेशर):

   • रणनीति: पहले ट्रैक्शन ($t$) को समाप्त करके एक प्रथम-स्तरीय शूर कॉम्प्लीमेंट $S_1 = A + C_1 \tilde{H}^{-1}C_2$ (एक काल्पनिक कठोरता वाला इलास्टिसिटी मैट्रिक्स) बनाया जाता है। फिर, विस्थापन ($u$) को समाप्त करके एक द्वितीय-स्तरीय शूर कॉम्प्लीमेंट $S_2 = T - Q_2 S_1^{-1} Q_1$ बनाया जाता है।
   • अनुमान: $S_1$ के व्युत्क्रम को AMG का उपयोग करके लागू किया जाता है। दूसरे स्तर के शूर कॉम्प्लीमेंट $S_2$ को एक डायगोनल फिक्स्ड-स्ट्रेस मैट्रिक्स $S_K$ द्वारा अनुमानित किया जाता है, जिसे बीजगणितीय रूप से या स्थानीय वॉल्यूम अनुमानों के माध्यम से गणना किया जा सकता है।
   • विशेषताएं: इस दृष्टिकोण में प्रति प्रीकंडीशनर एप्लिकेशन में दो AMG कॉल की आवश्यकता होती है (प्रति इटरेशन उच्च लागत)। हालांकि, यह अधिक मजबूती (robustness) प्रदान करता है क्योंकि पहला स्तर का शूर कॉम्प्लीमेंट $S_1$ फ्रैक्चर की स्थिति की परवाह किए बिना सैद्धांतिक रूप से पॉजिटिव डेफिनेट होने की गारंटी देता है। संभावित इंडेफिनेटनेस को दूसरे स्तर के शूर कॉम्प्लीमेंट में सीमित किया गया है, जो कि सिमेट्रिक है और इसे नेस्टेड डायरेक्ट सॉल्वर द्वारा प्रभावी ढंग से संभाला जा सकता है।

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

• फ्रेमवर्क डिज़ाइन: लेखक एक एकीकृत प्रीकंडीशनिंग फ्रेमवर्क प्रस्तुत करते हैं जो हाइड्रोलिक रूप से सक्रिय फ्रैक्चर के साथ स्टेबलाइज्ड कॉन्टैक्ट मैकेनिक्स से उत्पन्न होने वाले नॉन-सिमेट्रिक, इंडेफिनेट और समय-परिवर्तनीय $3 \times 3$ ब्लॉक जैकोबियन को संभालने में सक्षम है।
• तुलनात्मक विश्लेषण: यह शोध पत्र $t-p-u$ और $t-u-p$ रिडक्शन सीक्वेंस का एक कठोर सैद्धांतिक और संख्यात्मक तुलना प्रदान करता है। यह स्थापित करता है कि जबकि $t-p-u$ दृष्टिकोण सैद्धांतिक रूप से कुशल है, $t-u-p$ दृष्टिकोण जटिल, वास्तविक दुनिया के सिमुलेशन में, जहाँ फ्रैक्चर की स्थिति अप्रत्याशित होती है, अधिक मजबूती के लिए बेहतर है।
• स्केलेबिलिटी प्रदर्शन: कार्य यह प्रदर्शित करता है कि प्रस्तावित प्रीकंडिशनर एल्गोरिथमिक स्केलेबिलिटी प्राप्त करते हैं, जिससे समस्या के आकार के बढ़ने और फ्रैक्चर-आधारित अज्ञात (unknowns) और बल्क अज्ञात के अनुपात के बदलने पर भी क्रायलोव इटरेशन की संख्या लगभग स्थिर रहती है।

संख्यात्मक परिणाम
फ्रेमवर्क का परीक्षण तीन प्रयोगों के सेट पर किया गया था:

1. मजबूती (परीक्षण 1): एक सिंगल-फ्रैक्चर समस्या ने खुलासा किया कि जबकि $t-p-u$ दृष्टिकोण डायरेक्ट सॉल्वर के साथ अच्छा प्रदर्शन करता है, यह AMG का उपयोग करने पर विफल हो सकता है यदि शूर कॉम्प्लीमेंट इंडेफिनेट हो जाता है। $t-u-p$ दृष्टिकोण ने सभी परीक्षण की गई कॉन्फ़िगरेशन में अभिसरण (convergence) बनाए रखा।
2. वीक स्केलेबिलिटी (परीक्षण 2): समान रूप से रिफाइंड ग्रिड (~760,000 अज्ञात तक) और बदलते 2D-से-3D वेरिएबल अनुपात पर परीक्षणों ने दिखाया कि रैखिक सॉल्वर इटरेशन काउंट स्थिर रहा (औसतन ~17–24 GMRES इटरेशन), चाहे मेश रिफाइनमेंट या फ्रैक्चर नेटवर्क का सापेक्ष आकार कुछ भी हो।
3. वास्तविक दुनिया का अनुप्रयोग (परीक्षण 3): हाइड्रोलिक फ्रैक्चरिंग परिदृश्य में नौ फ्रैक्चरों के साथ एक टिल्टेड वेल (tilted well) के सिमुलेशन (1 मिलियन से अधिक अज्ञात) ने कम्प्यूटेशनल दक्षता की पुष्टि की। जटिल स्टिक/स्लिप/ओपन ट्रांज़िशन के कारण न्यूटन इटरेशन की कई एक्टिव-सेट स्टेप्स की आवश्यकता के बावजूद, रैखिक सॉल्वर स्थिर और कुशल रहा, जिसमें प्रति स्टेप औसतन 24.2 GMRES इटरेशन दर्ज किए गए।

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

लेखक विनम्रतापूर्वक निष्कर्ष निकालते हैं कि जबकि $t-p-u$ दृष्टिकोण सैद्धांतिक दक्षता और कम एप्लिकेशन लागत में मामूली लाभ प्रदान करता है, $t-u-p$ दृष्टिकोण जटिल फ्रैक्चर विकास के दौरान उत्पन्न होने वाले शूर कॉम्प्लीमेंट मैट्रिसेस की इंडेफिनेटनेस के खिलाफ अपनी बेहतर मजबूती के कारण व्यावहारिक अनुप्रयोगों के लिए पसंदीदा विकल्प है। यह कार्य बढ़ते बड़े और जटिल जियोमैकेनिकल समस्याओं को हल करने के लिए एक आधार स्थापित करता है, जिसका भविष्य का कार्य हाई-परफॉर्मेंस कंप्यूटिंग कार्यान्वयन और मल्टी-फेज फ्लो विस्तार की ओर निर्देशित है।