Linear Stability Analysis of Two-phase, Two-Component Flow in Porous Media

तकनीकी सारांश: सरंध्र माध्यमों में द्वि-चरणीय, द्वि-घटक प्रवाह का रैखिक स्थिरता विश्लेषण

समस्या विवरण
सरंध्र माध्यमों में तरल विस्थापन के दौरान होने वाली श्यानता उंगली (viscous fingering) अस्थिरता, उन्नत तेल रिकवरी (enhanced oil recovery), CO2 अनुक्रमण (sequestration), और भूजल उपचार जैसे अनुप्रयोगों की दक्षता को महत्वपूर्ण रूप से बाधित करती है। जबकि पूर्णतः अविমিশ্রणीय (immiscible) और पूर्णतः मिश्रणीय (miscible) विस्थापन के लिए रैखिक स्थिरता विश्लेषण पर व्यापक साहित्य मौजूद है, आंशिक रूप से मिश्रणीय प्रवाह (partially miscible flow)—जहाँ चरणों के बीच सीमित द्रव्यमान स्थानांतरण होता है—का मध्यवर्ती मामला काफी हद तक अनछुआ रहा है। विशेष रूप से, एक द्वि-चरणीय, द्वि-घटक प्रणाली में गुरुत्वाकर्षण, केशिका बलों (capillary forces), यांत्रिक विक्षेपण (mechanical dispersion), और अंतरा-चरण द्रव्यमान स्थानांतरण के जटिल परस्पर प्रभाव को ध्यान में रखने वाले विश्लेषणात्मक ढांचे की कमी है। यह अध्ययन ऐसे आंशिक रूप से मिश्रणीय प्रणालियों में अस्थिर विक्षोभ मोड (perturbation modes) के पूर्वानुमान पर ध्यान केंद्रित करते हुए, एक गैसीय तरल द्वारा तरल के विस्थापन के परिदृश्य पर आधारित है।

कार्यप्रणाली
लेखक एक समदैशिक (isotropic) सरंध्र माध्यम में द्वि-आयामी प्रवाह के लिए एक गणितीय मॉडल विकसित करते हैं, जो एक द्वि-घटक (घटक "a" और "b") के लिए ऊष्मप्रवैगिकी साम्य (thermodynamic equilibrium) को मानता है। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

1. शासी समीकरण (Governing Equations): मॉडल एक द्वि-चरणीय, द्वि-घटक प्रणाली के लिए द्रव्यमान संतुलन समीकरणों को व्युत्पन्न करता है, जिसमें डार्सी का नियम (Darcy's law), भिन्नात्मक प्रवाह (fractional flow), केशिका दबाव और यांत्रिक विक्षेपण शामिल हैं। एक प्रमुख सरलीकृत धारणा यह है कि एकल-चरण क्षेत्र में तरल चरण का घनत्व स्थिर रहता है, जबकि श्यानता (viscosity) संरचना के साथ बदलती है।
2. आधार अवस्था निर्माण (Base State Formulation): आधार अवस्था को एक आयामी, यात्रा करने वाली तरंग समाधान (travạching-wave solution/shock front) के रूप में परिभाषित किया गया है, जिसे द्रव्यमान स्थानांतरण के लिए अनुकूलित बकल-लेवरेट (Buckley-Leverett) दृष्टिकोण से प्राप्त किया गया है। शॉक विन्यास की गणना विधि-चरों (method of characteristics) का उपयोग करके की जाती है, जो दो-चरणीय क्षेत्र से शुद्ध-तरल क्षेत्र में संक्रमण के दौरान प्रवाह कार्यों में विच्छिन्नता (discontinuity) को ध्यान में रखती है।
3. रैखिक स्थिरता विश्लेषण: आधार अवस्था पर छोटे-आयाम वाले, तरंग-नुमा विक्षोभ (perturbations) आरोपित किए जाते हैं। eigenvalue समस्या बनाने के लिए शासी समीकरणों को रैखिक बनाया जाता है जहाँ eigenvalue विक्षोभ वृद्धि दर ($\sigma$) का प्रतिनिधित्व करता है और eigenfunction विक्षोभ प्रोफाइल का प्रतिनिधित्व करते हैं।
4. संख्यात्मक समाधान: एक महत्वपूर्ण चुनौती जो पहचानी गई है वह दो-चरणीय और एकल-चरण प्रवाह के बीच संक्रमण बिंदु पर eigenfunction के डेरिवेटिव में विच्छिन्नता है। इसे हल करने के लिए, लेखक संक्रमण बिंदु पर डेरिवेटिव के लिए विशिष्ट जंप स्थितियों (jump conditions) को व्युत्पन्न करते हैं। eigenvalue समस्या को 'मैच्ड इनिशियल वैल्यू प्रॉब्लम' (MIVP) पद्धति का उपयोग करके संख्यात्मक रूप से हल किया जाता है। इसमें दूर-क्षेत्र सीमाओं (far-field boundaries) से संक्रमण बिंदु की ओर रैखिक साधारण अंतर समीकरणों (ODEs) को एकीकृत करना और व्युत्पन्न जंप स्थितियों का उपयोग करके समाधानों को मिलाना शामिल है।

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

• नवीन ढांचा: यह कार्य, लेखकों की जानकारी के अनुसार, पहला रैखिक स्थिरता विश्लेषण प्रस्तुत करता है जो अंतरा-चरण द्रव्यमान स्थानांतरण को स्पष्ट रूप से ध्यान में रखते हुए एक वास्तविक द्वि-चरणीय, द्वि-घटक प्रणाली का विश्लेषण करता है।
• जंप स्थितियाँ (Jump Conditions): चरण-परिवर्तन विच्छिन्नता पर eigenfunction के डेरिवेटिव के लिए जंप स्थितियों का व्युत्पन्न, उन प्रणालियों में eigenvalue समस्या के सटीक समाधान की अनुमति देता है जहाँ प्रवाह कार्य विच्छिन्न होते हैं।
• व्यापक पैरामीटर अध्ययन: अध्ययन विस्थापन की स्थिरता पर प्रारंभिक संरचना, विक्षेपण (अनुदैर्ध्य और अनुप्रस्थ), मिश्रणीयता स्तर (द्वि-घटक अंतःक्रिया गुणांकों के माध्यम से), श्यानता अनुपात और गुरुत्वाकर्षण के प्रभावों की व्यवस्थित जांच करता है।

परिणाम

• द्रव्यमान स्थानांतरण का स्थिर प्रभाव: प्राथमिक निष्कर्ष यह है कि द्रव्यमान स्थानांतरण का प्रभाव मुख्य रूप से विस्थापन पर स्थिर करने वाला होता है। यह साम्य श्यानता अनुपात ($M_e$) को कम करके और शॉक गुणों को बदलकर (शॉक वेग $v_s$ को कम करके और गैस संतृप्ति $S_{gs}$ को बढ़ाकर) अधिकतम वृद्धि दर ($\sigma_{max}$) को कम करता है। यह स्थिरीकरण उच्च श्यानता अंतरों के लिए विशेष रूप से स्पष्ट है।
• गुरुत्वाकर्षण अंतःक्रियाएं: ऊपर की ओर विस्थापन (upward displacements) में, द्रव्यमान स्थानांतरण गुरुत्वाकर्षण-प्रेरित अस्थिरता को कम करता है। घनत्व का अंतर जो गुरुत्वाकर्षण अस्थिरता को संचालित करता है, द्रव्यमान स्थानांतरण को भी बढ़ाता है, जो अस्थिर करने वाले प्रभाव का प्रतिकार करता है। नीचे की ओर विस्थापन (downward displacements) में, गुरुत्वाकर्षण और द्रव्यमान स्थानांतरण दोनों ही स्थिर करने वाले तंत्र के रूप में कार्य करते हैं।
• विक्षेपण गैर-रैखिकता (Dispersion Non-linearity): अध्ययन आयामहीन अनुदैर्ध्य विक्षेपण गुणांक ($D^*_{\ell,xx}$) और अस्थिरता के बीच एक गैर-रैखिक संबंध की पहचान करता है। एक महत्वपूर्ण "सबसे खतरनाक" मान $D^*_{\ell,xx}$ मौजूद है जहाँ वृद्धि दर और कटऑफ वेवनंबर ($n_{cut}$) दोनों अधिकतम होते हैं। यह एक जटिल अंतःक्रिया का सुझाव देता जहाँ यांत्रिक विक्षेपण का मध्यम स्तर, केशिका बलों के सापेक्ष, बहुत कम और बहुत उच्च विक्षेपण व्यवस्थाओं की तुलना में वास्तव में अस्थिरता को बढ़ा सकता है।
• कटऑफ वेवनंबर: जबकि द्रव्यमान स्थानांतरण वृद्धि दर को महत्वपूर्ण रूप से कम करता है, इसका कटऑफ वेवनंबर पर प्रभाव कम स्पष्ट और अधिक जटिल है, जो विशिष्ट परिस्थितियों में पूर्णतः अविमिट्य मामलों की तुलना में उच्च $n_{cut}$ मान भी दे सकता है।

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