A Thermodynamically Consistent Free Boundary Model for Two-Phase Flows in an Evolving Domain with Bulk-Surface Interaction

यहाँ पैट्रिक क्नोफ़ (Patrik Knopf) और याडोंग लिउ (Yadong Liu) के शोध पत्र "A Thermodynamically Consistent Free Boundary Model for Two-Phase Flows in an Evolving Domain with Bulk-Surface Interaction" का विस्तृत तकनीकी सारांश दिया गया है।

1. समस्या विवरण (Problem Statement)

यह शोध पत्र एक विकसित होते डोमेन (evolving domain) $\Omega(t)$ के भीतर दो-चरणीय प्रवाह (two-phase flows - दो श्यान, असंपीड्य तरल पदार्थों का मिश्रण) के गणितीय मॉडलिंग को संबोधित करता है। एक स्थिर डोमेन वाले पारंपरिक मॉडलों के विपरीत, यह कार्य उन परिदृश्यों पर विचार करता है जहाँ डोमेन की सीमा $\Gamma(t) = \partial\Omega(t)$ तरल वेग क्षेत्र (fluid velocity field) के जवाब में चलती और विकृत होती है।

संबोधित किए गए प्रमुख चुनौतियाँ शामिल हैं:

• बल्क-सरफेस इंटरेक्शन (Bulk-Surface Interaction): बल्क तरल और सीमा सतह के बीच सामग्री के आदान-प्रदान (जैसे, सतह पर अवशोषण, अधिशोषण, या रासायनिक अभिक्रियाएं) का मॉडलिंग करना।
• परिवर्तनीय संपर्क कोण (Variable Contact Angles): 90-डिग्री के निश्चित संपर्क कोण की प्रतिबंधात्मक धारणा से आगे बढ़कर, दो चरणों के बीच के इंटरफेस को गतिशील कोण पर सीमा से मिलने की अनुमति देना।
• ऊष्मागतिक स्थिरता (Thermodynamic Consistency): यह सुनिश्चित करना कि मॉडल ऊष्मागतिकी के नियमों, विशेष रूप से द्रव्यमान संरक्षण (mass conservation) और ऊर्जा अपव्यय (energy dissipation) का एक चलते हुए डोमेन सेटिंग में कड़ाई से पालन करता है।
• चलती संपर्क रेखाएँ (Moving Contact Lines): उस संपर्क रेखा की गति का सटीक वर्णन करना जहाँ चरण इंटरफेस (phase interface) चलती हुई सीमा से मिलता है, जिसके लिए मानक 'नो-स्लिप' (no-slip) सीमा स्थिति को शिथिल करने की आवश्यकता होती है।

2. कार्यप्रणाली (Methodology)

लेखक दो अलग-अलग गणितीय दृष्टिकोणों का उपयोग करके नियंत्रण समीकरणों को व्युत्पन्न करते हैं, जो दोनों स्थानीय द्रव्यमान संतुलन नियमों (local mass balance laws) से शुरू होते हैं:

A. लैग्रेंज मल्टीप्लायर दृष्टिकोण (सामान्य मामला)

यह दृष्टिकोण घनत्व या द्रव्यमान प्रवाह पर प्रतिबंधात्मक धारणाओं के बिना पूर्ण मॉडल को व्युत्पन्न करने के लिए उपयोग किया जाता है।

• स्थानीय ऊर्जा अपव्यय (Local Energy Dissipation): लेखक बल्क और सतह दोनों के लिए स्थानीय ऊर्जा अपव्यय असमानताओं को सूत्रबद्ध करते हैं।
• रचनात्मक धारणाएँ (Constitutive Assumptions): वे बाधाओं (असंपीड्यता, द्रव्यमान संरक्षण) को लागू करने के लिए लैग्रेंज मल्टीप्लायर्स (रासायनिक क्षमता $\mu, \theta$ और सतह दबाव $q$) पेश करते हैं।
• फ्लक्स की पहचान (Flux Identification): यह सुनिश्चित करके कि अपव्यय असमानता बनी रहे, वे स्ट्रेस टेंसर, द्रव्यमान प्रवाह और सीमा स्थितियों के लिए रचनात्मक संबंधों की पहचान करते हैं। यह एक सामान्य व्युत्पत्ति की अनुमति देता है जिसमें बेमेल घनत्व और सतह के बीच द्रव्यमान हस्तांतरण शामिल है।

B. एनर्गेटिक वेरिएशनल दृष्टिकोण (EnVarA) (विशेष मामला)

यह दृष्टिकोण विशिष्ट सरलीकृत धारणाओं के तहत लागू किया जाता है: समान घनत्व ($\tilde{\rho}_1 = \tilde{\rho}_2$) और बल्क और सतह के बीच कोई द्रव्यमान प्रवाह नहीं।

• न्यूनतम क्रिया सिद्धांत (Least Action Principle): गतिज (kinetic) और मुक्त ऊर्जा फलन (free energy functionals) से जड़त्वीय (inertial) और संरक्षित बलों को व्युत्पन्न करने के लिए उपयोग किया जाता है।
• अधिकतम अपव्यय सिद्धांत (Onsager's Principle): रेले अपव्यय फलन (Rayleigh dissipation functional) से विपाची (dissipative) बलों (श्यानता, घर्षण) को व्युत्पन्न करने के लिए उपयोग किया जाता है।
• बल संतुलन (Force Balance): अंतिम समीकरण जड़त्वीय, संरक्षित और विपाची बलों के योग को शून्य के बराबर संतुलित करके प्राप्त किए जाते हैं।

3. मुख्य योगदान और मॉडल सूत्रीकरण (Key Contributions and Model Formulation)

यह शोध पत्र गतिशील सीमा स्थितियों के साथ युग्मित एक नवीन नेवियर–स्टोक्स–कैन–हिलियर्ड (NS-CH) प्रणाली प्रस्तुत करता है।

नियंत्रित प्रणाली (The Governing System)

मॉडल में निम्नलिखित शामिल हैं:

1. बल्क डायनेमिक्स (Bulk Dynamics):
   • नेवियर–स्टोक्स समीकरण: चरण-निर्भर घनत्व $\rho(\phi)$ और श्यानता $\nu(\phi)$ के साथ मिश्रण के संवेग संतुलन का वर्णन करता है। इसमें इंटरफेशियल तनाव (कोटरवेट स्ट्रेस) के कारण एक अतिरिक्त स्ट्रेस टर्म शामिल है।
   • संवहनी कैन–हिलियर्ड समीकरण (Convective Cahn–Hilliard Equation): तरल वेग के साथ युग्मित, बल्क फेज-फील्ड $\phi$ के विकास का वर्णन करता है।
2. सतह डायनेमिक्स (Surface Dynamics):
   • सतह कैन–हिलियर्ड समीकरण: एक गतिशील सीमा स्थिति जो सतह फेज-फील्ड $\psi$ के विकास का वर्णन करती है।
   • युग्मन (Coupling): बल्क और सतह के समीकरणों को सीमा स्थितियों के माध्यम से जोड़ा जाता है जो सामग्री हस्तांतरण (पैरामीटर $K$ और $L$ द्वारा नियंत्रित) की अनुमति देते हैं।
3. सीमा स्थितियाँ (Boundary Conditions):
   • सामान्यीकृत नेवियर स्लिप (Generalized Navier Slip): 'नो-स्लिप' स्थिति को प्रतिस्थापित करता है ताकि चलती संपर्क रेखाओं की अनुमति मिल सके। स्पर्शरेखीय वेग (tangential velocity), सतह तनाव प्रवणता (surface tension gradients) से संबंधित एक बल पद के साथ स्पर्शरेखीय स्ट्रेस के समानुपाती होता है।
   • गतिशील सीमा स्थितियाँ (Dynamic Boundary Conditions): डोमेन की सामान्य वेग (normal velocity) तरल वेग के सामान्य घटक ($u \cdot n = V$) द्वारा निर्धारित होती है, जिससे यह एक फ्री बाउंड्री समस्या बन जाती है।
   • संपर्क कोण गतिशीलता (Contact Angle Dynamics): $\phi$ और $\mu$ के लिए सीमा स्थितियाँ (जो $K, L$ द्वारा नियंत्रित हैं) यह अनुमति देती हैं कि संपर्क कोण 90 डिग्री पर स्थिर रहने के बजाय गतिशील रूप से विकसित हो।

प्रमुख भौतिक गुण (Key Physical Properties)

व्युत्पन्न मॉडल निम्नलिखित का पालन करता है:

• द्रव्यमान संरक्षण (Mass Conservation): कुल द्रव्यमान संरक्षित रहता है; यदि $L=\infty$ है, तो बल्क और सतह के द्रव्यमान अलग-अलग संरक्षित रहते हैं।
• ऊर्जा अपव्यय (Energy Dissipation): प्रणाली एक वैश्विक ऊर्जा अपव्यय नियम का पालन करती है, जो ऊष्मागतिक स्थिरता को सिद्ध करती है। कुल ऊर्जा (गतिज + बल्क मुक्त ऊर्जा + सतह मुक्त ऊर्जा) श्यानता, घर्षण और विसरण के कारण समय के साथ घटती है।
• आयतन और क्षेत्रफल संरक्षण (Volume and Area Preservation): असंपीड्यता बाधाओं ($\text{div } u = 0$ और $\text{div}_\Gamma u = 0$) के कारण, डोमेन का आयतन और सीमा का क्षेत्रफल संरक्षित रहता है।

4. परिणाम (Results)

• व्युत्पत्ति सत्यापन (Derivation Verification): लेखकों ने सफलतापूर्वक दो अलग-अलग विधियों का उपयोग करके एक ही प्रणाली को व्युत्पन्न किया, जो मॉडल की मजबूती को प्रमाणित करता है।
• सामान्यीकरण (Generalization): मॉडल मौजूदा साहित्य को सामान्यीकृत करता है:
  • $K=L=\infty$ सेट करने पर यह निश्चित 90-डिग्री संपर्क कोण और बिना सतह द्रव्यमान हस्तांतरण वाले मॉडलों को पुनः प्राप्त करता है।
  • डोमेन सीमा को स्थिर करने पर ($V=0$), यह गतिशील सीमा स्थितियों वाले मानक NS-CH मॉडलों को पुनः प्राप्त करता है।
  • AGG मॉडल (Abels-Garcke-Grün) को एक स्थिर डोमेन में एक विशेष मामले के रूप में पुनः प्राप्त किया जाता है।
• पैरामीटर विश्लेषण (Parameter Analysis): शोध पत्र विस्तार से बताता है कि कैसे पैरामीटर $K$ और $L$ भौतिक व्यवहार को नियंत्रित करते हैं:
  • $K=0$: बल्क और सतह के बीच फेज-फील्ड के बीच एक रैखिक संबंध लागू करता है (तीव्र साम्यावस्था)।
  • $L \in (0, \infty)$: द्रव्यमान हस्तांतरण (अवशोषण/अधिशोषण) की एक सीमित दर की अनुमति देता है।
  • $L=\infty$: सीमा के पार शून्य द्रव्यमान प्रवाह को लागू करता है।

5. महत्व (Significance)

• जैविक और भौतिक अनुप्रयोग: यह मॉडल जैविक प्रणालियों (जैसे, चलती बैक्टीरिया, कोशिका झिल्ली) और सॉफ्ट मैटर भौतिकी (जैसे, विकृत सबस्ट्रेट्स पर बूंदें) के अनुकरण के लिए विशेष रूप से प्रासंगिक है, जहाँ तरल गति के कारण डोमेन ज्यामिति गतिशील रूप से बदलती है।
• सैद्धांतिक प्रगति: यह विकसित होने वाले मैनिफोल्ड्स (manifolds) पर दो-चरणीय प्रवाह के लिए एक कठोर ऊष्मागतिक ढांचा प्रदान करता है, जो एक ऐसा विषय है जिसमें बल्क-सरफेस द्रव्यमान हस्तांतरण और परिवर्तनीय संपर्क कोण दोनों को एक साथ संभालने वाला सुसंगत स्वरूप की कमी थी।
• लचीलापन: $K$ और $L$ के समावेश से मॉडल को विभिन्न भौतिक परिदृश्यों, अभेद्य सीमाओं से लेकर अत्यधिक अवशोषक सतहों तक, के लिए ट्यून करने की अनुमति मिलती है।
• भावी दिशाएँ: यह शोध पत्र भविष्य के विश्लेषणात्मक कार्य (वेल-पोज्डनेस सिद्ध करना) और इन जटिल युग्मित गतिकी को विकसित होती ज्यामिति पर अनुमानित करने के लिए संख्यात्मक सिमुलेशन (numerical simulations) के लिए मार्ग प्रशस्त करता है।

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