Density evolution at fluid-fluid interfaces: A generalized Gibbs-Duhem theory

तकनीकी सारांश: तरल-तरल इंटरफेस पर घनत्व का विकास: एक सामान्यीकृत गिब्स-डुहेम सिद्धांत

समस्या विवरण
यह शोध पत्र शास्त्रीय गिब्स ऊष्मागतिकी (Gibbs thermodynamics) और न्यूटनियन यांत्रिकी के बीच एक मौलिक विच्छेद की पहचान करता है, विशेष रूप से इंटरफेशियल प्रणालियों के संबंध में। शास्त्रीय न्यूटनियन यांत्रिकी पदार्थ की गति को $\rho \frac{du}{dt} = -\nabla p$ के माध्यम से वर्णित करती है, जबकि शास्त्रीय गिब्स-डुहेम संबंध ($-dp + \phi d\mu + s_\tau d\tau = 0$) कड़ाई से अर्ध-स्थैतिक (quasi-static), साम्यावस्था प्रक्रियाओं पर लागू होता है, जो गतिज ऊर्जा योगदान की उपेक्षा करता है। इसके अतिरिक्त, मानक तरल गतिकी अक्सर असंपीड्यता (incompressibility) की धारणा पर निर्भर करती है ($\frac{d\rho}{dt} = 0$), जो स्थानीय घनत्व विविधताओं और अतिरिक्त-आयतन प्रभावों (जैसे, जल-गैस मिश्रण) को पकड़ने में विफल रहती है। द्रव्यमान संरक्षण से व्युत्पन्न मौजूदा निरंतरता समीकरण (continuity equations), एंट्रॉपी विकास से जुड़े ऊष्मागतिक बाधाओं को स्पष्ट रूप से शामिल नहीं करते हैं, जो बहु-घटक इंटरफेशियल प्रणालियों में घनत्व की गतिशीलता का वर्णन करने की उनकी क्षमता को सीमित करते हैं।

कार्यप्रणाली
लेखक गतिज ऊर्जा योगदान को शामिल करने के लिए गिब्स-डुहेम संबंध का विस्तार करके एक सामान्यीकृत ढांचा विकसित करते हैं। कार्यप्रणाली इस प्रकार है:

1. प्रणाली परिभाषा: दो अविभाज्य तरल चरणों वाले एक बंद निकाय का मॉडल बनाया गया है। इंटरफेस को एक परिमित मोटाई वाले क्षेत्र के रूप में माना गया है जहाँ अणु सह-अस्तित्व में होते हैं।
2. आयतन और घनत्व का पुनर्गठन: लेखक स्थानीय असंपीड्यता बाधा को शिथिल करते हैं। शून्य अतिरिक्त आयतन मानने के बजाय, वे एक अतिरिक्त आयतन पद ($\Delta v_e$) पेश करते हैं और इसे प्रभावी आंशिक आयतनों ($\Delta v'_w, \Delta v'_a$) में पुनर्वितरित करते हैं। यह आंशिक घनत्व को स्थानीय दबाव परिवर्तनों के साथ बदलने की अनुमति देता है, जिससे आयतन संरचना ($\phi$) और कुल घनत्व ($\rho$) में अतिरिक्त आयतन की जानकारी समाहित हो जाती है।
3. मुक्त ऊर्जा फलन (Free Energy Functional): एक मुक्त ऊर्जा फलन $E$ का निर्माण किया गया है, जिसमें हेल्महोल्ट्ज़-प्रकार का मुक्त ऊर्जा घनत्व, गतिज मुक्त ऊर्जा और वॉल फ्री एनर्जी शामिल है। महत्वपूर्ण रूप से, एंट्रॉपी $s$ को तीन घटकों में विभाजित किया गया है: बोल्ट्ज़मैन मिश्रण एंट्रॉपी ($s_\phi$), तापीय एंट्रॉपी ($s_\tau$), और वेग एंट्रॉपी ($s_u$)।
4. वेरिएशनल कैलकुलस: स्थिर आयतन बाधाओं के तहत मुक्त ऊर्जा फलन पर वेरिएशनल कैलकुलस लागू करके, लेखक तापमान, संरचना और वेग के संयुग्मी विभव (conjugate potentials) प्राप्त करते हैं।
5. समीकरणों का व्युत्पन्न:
   • एक संशोधित गिब्स-डुहेम समीकरण व्युत्पन्न किया गया है, जो रासायनिक विभव, वेग विभव और तापमान विभव के परिवर्तनों से दबाव परिवर्तनों को जोड़ते हुए मुक्त ऊर्जा के कुल अवकलज (total derivative) को लागू करता है।
   • ऊष्मा, द्रव्यमान और वेग प्रसार के लिए गतिज समीकरण विक्षोभ सिद्धांत (dissipation principle) का उपयोग करके प्राप्त किए गए हैं।
   • संवेग विकास समीकरण को संशोधित गिब्स-ड्यूहेम संबंध से व्युत्पन्न किया गया है।
   • एक नया घनत्व विकास समीकरण संवेग विकास समीकरण से वेग विक्षोभ समीकरण को घटाकर तैयार किया गया है।

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

• सामान्यीकृत गिब्स-डुहेम संबंध: यह पत्र एक एकीकृत संबंध व्युत्पन्न करता है जो गतिज योगदान को शामिल करता है, जो शास्त्रीय ऊष्मागतिकी और न्यूटनियन यांत्रिकी के बीच गणितीय रूप से सेतु बनाता है।
• एकीकृत ढांचा: यह स्थापित करता है कि सामान्यीकृत गिब्स-डुहेम संबंध संवेग विकास समीकरण के गणितीय रूप से समकक्ष है, जो ऊष्मागतिक प्रेरक बलों और यांत्रिक गति के बीच एक सीधा संबंध प्रदान करता है।
• वैकल्पिक घनत्व विकास समीकरण: लेखक तरल-तरल इंटरफेस पर घनत्व के लिए एक नया विकास समीकरण प्रस्तावित करते हैं:
  $$u \frac{d\rho}{dt} = -\nabla p - \sum_i \phi_i \nabla \mu_i - \tau \nabla \mu_\tau - \rho \nabla \cdot [D_1 \nabla \mu_u + D_2 (\nabla \cdot \mu_u)I]$$
  यह समीकरण स्पष्ट रूप से दबाव, रासायनिक विभव प्रवणता (मिश्रण), तापीय प्रवणता और वेग एंट्रॉपी प्रसार द्वारा संचालित घनत्व विविधताओं को ध्यान में रखता है।
• श्यानता उपचार: नेवियर-स्टोक्स फॉर्मूलेशन के विपरीत जहाँ श्यान प्रभाव स्पष्ट यांत्रिक बल होते हैं, यह ढांचा प्रसार समीकरणों में वेग एंट्रॉपी के माध्यम से विक्षोभ (dissipation) को ऊष्मागतिक विवरण के भीतर समाहित करता है।

परिणाम
प्रस्तावित घनत्व विकास समीकरण की वैधता यह दिखाकर प्रदर्शित की गई है कि यह उपयुक्त सीमा मामलों में स्थापित भौतिक नियमों को स्वाभाविक रूप से पुनः प्राप्त करता है:

• ध्वनि की गति: आइसेंट्रोपिक स्थितियों के तहत, समीकरण ध्वनि की गति की क्लासिक परिभाषा ($u = \sqrt{dP/d\rho}$) को पुनः प्राप्त करता है।
• बरनौली का नियम: समीकरण एक सामान्यीकृत बरनौली समीकरण देता है। यह मानक बरनौली सिद्धांत ($p + \frac{1}{2}\rho u^2 = \text{constant}$) को पुनः प्राप्त करता है जब वेग एंट्रॉपी स्थिर होती है, लेकिन ऊर्जा परिवर्तन को भी ध्यान में रखता है जब वेग एंट्रॉपी परिवर्तित होती है।
• अवस्था का समीकरण (EOS): बिना किसी मैक्रोस्कोपिक संचलन के स्थिर अवस्था में एक आदर्श गैस के लिए, यह आदर्श गैस समीकरण (एक लघुगणकीय घनत्व कारक के साथ) को पुनः प्राप्त करता है।
• उच्च घनत्व अनुपात: यह दावा किया गया है कि यह ढांचा "क्लासिक हाई डेंसिटी रेश्यो समस्या" को संबोधित करता है, विशेष रूप से जल-गैस इंटरफेस पर जहाँ घनत्व अनुपात लगभग 1000 है।

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

1. जब घनत्व समान होता है तो शास्त्रीय निरंतरता समीकरण को एक सीमा मामले के रूप में पुनः प्राप्त करता है।
2. उन इंटरफेशियल प्रणालियों के लिए लागू रहता है जिनमें संरचना प्रवणता, एंट्रॉपी विकास और घनत्व विविधताएं हैं जहाँ शास्त्रीय निरंतरता समीकरण अपर्याप्त है।
3. गैर-साम्यावस्था ऊष्मागतिकी और यांत्रिक प्रक्रियाओं का एक एकीकृत विवरण प्रदान करता है, जो चरण इंटरफेस और बहु-घटक तरल गतिकी के ऊष्मागतिकी के लिए नई अंतर्दृष्टि प्रदान करने की क्षमता रखता है।

लेखक उल्लेख करते हैं कि इस कार्य को कोई बाहरी वित्त पोषण प्राप्त नहीं हुआ है और वे आशा करते हैं कि यह ढांचा गैर-साम्यावस्था ऊष्मागतिकी के लिए व्यापक रूप से प्रासंगिक होगा, हालांकि वे इन सैद्धांतिक विस्तारों के परे विशिष्ट प्रयोगात्मक सत्यापन या भविष्य के अनुप्रयोगों का प्रस्ताव नहीं करते हैं।