A Finite Analytic Method for Two-Dimensional Power-Law Fluid Flow in Heterogeneous Porous Media Using Unstructured Grids

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

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

शोधकर्ताओं ने एक विशिष्ट प्रकार के तरल पर ध्यान केंद्रित किया जिसे 'पावर-लॉ फ्लूइड' (power-law fluid) कहा जाता है। ये वे तरल पदार्थ हैं जहाँ उन्हें धकेलने वाले बल और उनकी गति के बीच का संबंध एक विशिष्ट, गैर-रैखिक नियम का पालन करता है। वास्तविक दुनिया में, यह पदार्थों की एक विस्तृत श्रृंखला को कवर करता है, जैसे कि गाढ़े पॉलीमर समाधान या कच्चे तेल के कुछ प्रकार। चुनौती तब उत्पन्न होती है जब ये तरल पदार्थ सरंध्र माध्यमों (porous media)—यानी छोटे छिद्रों वाली चट्टानों—से होकर बहते हैं, जहाँ पारगम्यता (permeability), या प्रवाह की सुगमता, एक स्थान से दूसरे स्थान पर अचानक बदल जाती है। ऐसे वातावरण में, तरल का दबाव सुचारू रूप से नहीं बदलता है; इसके बजाय, यह आसान-प्रवाह और कठिन-प्रवाह वाली चट्टान के बीच की सीमा पर एक तीक्ष्ण उछाल या "सिंगुलैरिटी" (singularity) बनाता है। पारंपरिक कंप्यूटर मॉडल, जो चट्टान को छोटे त्रिभुजों के ग्रिड में विभाजित करते हैं, इन उछालों के साथ संघर्ष करते हैं। वे मान लेते हैं कि दबाव बिंदुओं के बीच एक सीधी रेखा में बदलता है, जो इन सीमाओं पर भौतिक रूप से गलत है। परिणामस्वरूप, मॉडल प्रवाह दर का बहुत कम आकलन करते हैं, और उन्हें सत्य के एक मोटे अनुमान तक पहुँचने के लिए भी बहुत अधिक सूक्ष्म ग्रिड बिंदुओं की आवश्यकता होती है।

इसे ठीक करने के लिए, टीम ने एक 'फाइनाइट एनालिटिक मेथड' (Finite Analytic Method) विकसित किया। मुख्य विचार यह है कि सीमा पर तरल के व्यवहार का अनुमान लगाना बंद करें और इसके बजाय उस व्यवहार के ज्ञात, सटीक गणितीय विवरण का उपयोग करें। शोधकर्ताओं ने महसूस किया कि जहाँ चट्टान के गुण बदलते हैं, उस बिंदु के पास तरल का दबाव एक विशिष्ट वक्रीय पैटर्न (curved pattern) का अनुसरण करता है, जो तरल की पावर-लॉ प्रकृति पर निर्भर करता है। उन्होंने इस प्रवाह के सबसे महत्वपूर्ण हिस्से के लिए इस वक्र के सटीक आकार, जिसे 'इगेनफंक्शन' (eigenfunction) कहा जाता है, की गणना की। इस गणना में उन जटिल समीकरणों को हल करना शामिल था जो एक वेज-आकार (wedge-shaped) के क्षेत्र में तरल की गति का वर्णन करते हैं जहाँ चट्टान के गुण भिन्न होते हैं। इन समीकरणों के सबसे छोटे और सबसे महत्वपूर्ण समाधान को खोजकर, उन्होंने उस प्रमुख तरीके को पकड़ लिया जिससे तरल ठीक उस समस्या वाले स्थान पर व्यवहार करता है।

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

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

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

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