A quasi-incompressible Cahn-Hilliard-Darcy model for two immiscible fluids in porous media

तकनीकी सारांश: दो असंगत तरल पदार्थों के लिए एक अर्ध-असंपीड्य (Quasi-Incompressible) काह्न-हिलार्ड-डार्सी मॉडल

समस्या विवरण
यह अध्ययन सरंध्र माध्यमों (porous media) के भीतर दो-चरणीय असंगत तरल पदार्थों के गणितीय मॉडलिंग को संबोधित करता है। जबकि शास्त्रीय मुस्काट (Muskat) समस्या तरल इंटरफेस को अनंत रूप से पतली विच्छेदन रेखाओं के रूप में प्रभावी ढंग से वर्णित करती है, यह परिमित-समय विलक्षणताओं (finite-time singularities) के प्रति संवेदनशील है, जैसे कि चिकनी प्रोफाइल का अनियंत्रित विकास और सुव्यवस्थित प्रतीत होने वाले शासन में भी अस्थिरता। इन सीमाओं को दूर करने के लिए, डिफ्यूज-इंटरफेस (फेज-फील्ड) सूत्रीकरण का उपयोग तीक्ष्ण विच्छेदन को संकीर्ण संक्रमण परतों से बदलने के लिए किया जाता है, जिससे गतिशीलता को नियमित किया जा सके।

मौजूदा साहित्य ने समान और असमान घनत्वों के लिए काह्न-हिलार्ड-नेवियर-स्टोक्स (CHNS) प्रणाली, और काह्न-हिलार्ड-हीले-शॉ/डार्सी (CHHS/CHD) मॉडलों को व्यापक रूप से कवर किया है। हालाँकि, एक अर्ध-असंपीड्य काह्न-हिलार्ड-डार्सी (qCHD) मॉडल के लिए एक पूर्ण सैद्धांतिक ढांचा जो एक विलक्षण लघुगणकीय फ्लॉरी-हगिंस्स क्षमता (singular logarithmic Flory–Huggins potential) और गैर-समान तरल गुणों (विशेष रूप से भिन्न घनत्व) को समाहित करता है, एक अनसुलझा मुद्दा बना हुआ है। अधिकांश मौजूदा डार्सी-प्रकार के मॉडल या तो सख्त असंपीड्यता ($\nabla \cdot u = 0$) मानते हैं या बहुपद मुक्त ऊर्जा क्षमता (polynomial free energy potentials) का उपयोग करते हैं, जिनमें चरण सीमाओं के पास लघुगणकीय क्षमता की भौतिक यथार्थता का अभाव होता है।

कार्यप्रणाली
लेखक qCHD मॉडल को व्युत्पन्न करते हैं और निम्नलिखित कार्यप्रणालीगत चरणों के माध्यम से कमजोर समाधानों (weak solutions) के वैश्विक अस्तित्व को स्थापित करते हैं:

1. वैरिएशनल व्युत्पत्ति (Variational Derivation):

   • मॉडल को ओन्सेगर के एक्सट्रीमम सिद्धांत (रेले के न्यूनतम ऊर्जा अपव्यय सिद्धांत के समकक्ष) का उपयोग करके व्युत्पन्न किया गया है।
   • प्रणाली आयतन अंश अंतर $\phi$ को ऑर्डर पैरामीटर के रूप में और द्रव्यमान-औसत मिश्रण वेग $u$ का उपयोग करती है, जो घनत्व के अंतर के कारण गैर-सोलेनोइडल ($\nabla \cdot u \neq 0$) है।
   • कुल ऊर्जा कार्यात्मक $E(\phi)$ में फ्लॉरी-हगिंस्स समरूप मुक्त ऊर्जा $F(\phi)$, एक ग्रेडिएंट पद (जो इंटरफेशियल ऊर्जा का प्रतिनिधित्व करता है), और एक गुरुत्वाकर्षण क्षमता शामिल है।
   • निरंतरता बाधा को शिथिल करने के लिए एक लैग्रेंज मल्टीप्लायर $p$ पेश करके और रेलेियन (Rayleighian) को न्यूनतम करके, लेखक गति संतुलन (डार्सी का नियम), द्रव्यमान संरक्षण, और ऑर्डर पैरामीटर के विकास को नियंत्रित करने वाली युग्मित समीकरणों की प्रणाली व्युत्पन्न करते हैं।

2. मॉडल सूत्रीकरण (Model Formulation):

   • परिणामी qCHD प्रणाली में शामिल है:
     • एक संशोधित डार्सी समीकरण जो वेग, दबाव, रासायनिक क्षमता और गुरुत्वाकर्षण को जोड़ता है।
     • मिश्रण घनत्व के लिए एक निरंतरता समीकरण।
     • ऑर्डर पैरामीटर $\phi$ के लिए एक काह्न-हिलार्ड प्रकार का समीकरण, जो रासायनिक क्षमता $\mu_p = \mu + \alpha p$ द्वारा संचालित होता है, जहाँ $\alpha$ एटवुड नंबर (Atwood number) है।
     • मुक्त ऊर्जा घनत्व $F(\phi)$ एक विलक्षण लघुगणकीय फलन है: $F(\phi) = \frac{\theta_0}{2}[(1+\phi)\ln(1+\phi) + (1-\phi)\ln(1-\phi)] - \frac{\theta_c}{2}\phi^2 + C$।

3. अस्तित्व के लिए विश्लेषणात्मक ढांचा (Analytical Framework for Existence):

   • समय-विविक्तकरण (Time-Discretization): मुक्त ऊर्जा क्षमता ($F = F_v + F_n$) के कॉनवेक्स-कॉनकेव विभाजन का उपयोग करके एक टाइम-मार्चिंग स्कीम बनाई गई है।
     कड़ाई से एकात्मकता (strict monotonicity) और कोएर्सिविटी (coercivity) पर निर्भर करते हुए, संबद्ध ऑपरेटर के समाधान के लिए एक फिक्स्ड पॉइंट तर्क का उपयोग किया जाता है।
   • अ प्रायोर अनुमान (A Priori Estimates): लेखक ऊर्जा अपव्यय नियम और डार्सी समीकरणों से प्राप्त दबाव के विशिष्ट अनुमानों का लाभ उठाते हुए सन्निकट समाधानों के लिए समान अनुमान प्राप्त करते हैं।
   • संकुचन (Compactness): ऑर्डर पैरामीटर की प्रबल अभिसरण (strong convergence) को स्थापित करने के लिए ऑबिन-लियंस-साइमन लेम्मा का उपयोग किया जाता है, जिससे गैर-रेखीय पदों, विशेष रूप से विलक्षण लघुगणकीय क्षमता में सीमा की ओर जाने की अनुमति मिलती है।

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

• qCHD मॉडल की व्युत्पत्ति: यह शोध पत्र एक व्यवस्थित qCHD मॉडल प्रस्तुत करता है जो स्थानीय द्रव्यमान संरक्षण और एक ऊर्जा-अपव्यय नियम को संतुष्ट करता है। यह मॉडल एटवुड नंबर के माध्यम से घनत्व अंतर को ध्यान में रखता है और एक द्रव्यमान-औसत वेग का उपयोग करता है।
• कमजोर समाधानों का वैश्विक अस्तित्व: प्राथमिक गणितीय परिणाम 2D और 3D डोमेन में qCHD प्रणाली के लिए कमजोर समाधानों के वैश्विक अस्तित्व का प्रमाण है। यह भौतिक रूप से यथार्थवादी फ्लॉरी-हगिंस्स लघुगणकीय क्षमता के लिए प्राप्त किया गया है, जो यह सुनिश्चित करता है कि ऑर्डर पैरामीटर भौतिक सीमाओं $(-1, 1)$ के भीतर रहे।
• तीक्ष्ण इंटरफेस सीमा (Sharp Interface Limit): औपचारिक मैच की गई एसिम्प्टोटिक एक्सपेंशन (matched asymptotic expansion) के माध्यम से, लेखक यह प्रदर्शित करते हैं कि qCHD मॉडल तीक्ष्ण इंटरफेस सीमा में एक संशोधित मुस्काट समस्या की ओर अभिसरित होता है, जो शास्त्रीय इंटरफेशियल गतिशीलता के साथ इसकी निरंतरता की पुष्टि करता है।
• ऊर्जा स्थिरता: व्युत्पन्न प्रणाली एक कठोर ऊर्जा-अपव्यय नियम का पालन करती है, जो यह सुनिश्चित करती है कि श्यान अपव्यय (viscous dissipation) और रासायनिक प्रसार के कारण कुल ऊर्जा समय के साथ घटती है।

महत्व
यह शोध पत्र दावा करता है कि यह पहला अध्ययन है जो एक अर्ध-असंपीड्य काह्न-हिलार्ड-डार्सी प्रणाली प्रस्तुत करता है जो साथ ही एक विलक्षण लघुगणकीय क्षमता और गैर-समान तरल गुणों (असमान घनत्व) को समाहित करता है।

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