A Nonmatching-Grid Discontinuous Galerkin Method for Nonlinear Reactive Transport in Hybrid-Dimensional Fractured Media

हमारे पैरों के नीचे छिपी दुनिया में, पानी और घुले हुए रसायन चट्टान के एक जटिल भूलभुलैया के माध्यम से चलते हैं। यह चट्टान शायद ही कभी एकसमान होती है; यह अक्सर एक ठोस ब्लॉक होती है, जिसे 'मैट्रिक्स' कहा जाता है, जो महीन दरारों द्वारा क्रॉस की जाती है जिन्हें 'फ्रैक्चर' (fractures) कहा जाता है। ये फ्रैक्चर भूमिगत राजमार्गों की तरह होते हैं। वे चट्टान के माध्यम से तरल पदार्थों को तेजी से प्रवाहित करने वाले सुपर-फास्ट चैनलों के रूप में कार्य कर सकते हैं, या पूरी गति को रोकने वाले अवरोधों के रूप में भी, यह उनकी प्रकृति पर निर्भर करता है। क्योंकि ये दरारें अक्सर चौड़ाई में सूक्ष्म होती हैं लेकिन लंबाई में कई किलोमीटर तक फैली हो सकती हैं, इसलिए इनके माध्यम से पदार्थों के संचलन का मॉडल बनाना वैज्ञानिकों के लिए एक बड़ी चुनौती है। यदि वे एक ऐसा कंप्यूटर मॉडल बनाने की कोशिश करते हैं जो दरार के खुलेपन के हर छोटे विवरण को पकड़ सके, तो गणना इतनी भारी हो जाती है कि सबसे शक्तिशाली सुपरकंप्यूटर भी उसे हल नहीं कर पाते। इस समस्या से निपटने के लिए, शोधकर्ता 'हाइब्रिड-डायमेंशनल मॉडल' नामक एक चतुर शॉर्टकट का उपयोग करते हैं। दरार को तीन-आयामी पतली स्लाइस के रूप में मानने के बजाय, वे इसे एक दो-आयामी रेखा या एक एक-आयामी रेखा के रूप में समतल कर देते हैं जो चट्टान के भीतर स्थित होती है। यह उन्हें दरार की सूक्ष्म मोटाई को हल किए बिना प्रवाह का अनुकरण करने की अनुमति देता है।

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

अल्वारो एच. सालास और कार्लोस डी. अकोस्टा (यूनिवर्सिडाड नैशनल डी कोलंबिया के) द्वारा किया गया एक नया अध्ययन इस कठोरता को एक सरल लेकिन गहन प्रश्न पूछकर संबोधित करता है: क्या होगा यदि चट्टान और दरार को एक सीध में होने की आवश्यकता ही न हो? उन्होंने एक ऐसी विधि विकसित की जहाँ चट्टान और फ्रैक्चर को पूरी तरह से अलग आकारों और आकृतियों के टुकड़ों में विभाजित किया जा सकता है, जो एक दूसरे से स्वतंत्र रूप से उत्पन्न होते हैं, और फिर भी वे जानकारी का पूर्णतः आदान-प्रदान कर सकते हैं। कल्पना कीजिए कि आप एक बड़े, मोटे बाल्टी से एक बारीक, विस्तृत छलनी में पानी डालने की कोशिश कर रहे हैं। आमतौर पर, आपको पानी को बिना गिराए सुचारू रूप से स्थानांतरित करने के लिए एक विशेष फनल (कीप) की आवश्यकता होगी। शोधकर्ताओं ने सीधे पानी डालने का तरीका खोज लिया, यह गणना करते हुए कि बाल्टी के किनारे का छलनी के प्रत्येक बिंदु से जहाँ भी संपर्क होता है, वहां से ठीक कितना पानी गुजरता है, चाहे छेद कितने भी बेमेल क्यों न हों।

टीम ने एक गणितीय ढांचा तैयार किया जो चट्टान और फ्रैक्चर के बीच रसायनों के आदान-प्रदान को एक प्रत्यक्ष भौतिक अंतःक्रिया के रूप में मानता है। दोनों ग्रिडों को मिलाने या डेटा को अनुवादित करने के लिए किसी मध्यस्थ का उपयोग करने के बजाय, उन्होंने एक ऐसी प्रणाली बनाई जो उन सूक्ष्म, ओवरलैपिंग खंडों पर स्थानांतरण की गणना करती है जहाँ चट्टान का किनारा और फ्रैक्चर की सतह वास्तव में मिलते हैं। इस दृष्टिकोण को 'नॉनमैचिंग-ग्रिड विधि' (nonmatching-grid method) कहा जाता है। अपने सिमुलेशन में, उन्होंने इस प्रणाली का परीक्षण जटिल रासायनिक प्रतिक्रियाओं के साथ किया, जहाँ पदार्थ चलते समय बदलते रहते हैं, जो कि मॉडलिंग के लिए एक अत्यंत कठिन परिदृश्य है। उन्होंने गणितीय रूप से सिद्ध किया कि यह सीधा संबंध स्थिर और सटीक है, भले ही चट्टान को बहुत बड़े ब्लॉकों में विभाजित किया गया हो और फ्रैक्चर को बहुत महीन खंडों में, या इसके विपरीत। उनकी विधि की स्थिरता इन ब्लॉकों के आकार के अनुपात पर निर्भर नहीं करती है, जिसका अर्थ है कि यदि एक ग्रिड दूसरे की तुलना में बहुत महीन है, तो कंप्यूटर क्रैश नहीं होता या गलत परिणाम नहीं देता।

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

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

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

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