Numerical analysis of the thermal relaxation of the dense gas between two parallel plates: the free energy monotonicity for the Enskog equation

यहाँ शोध पत्र "Numerical analysis of the thermal relaxation of the dense gas between two parallel plates: the free energy monotonicity for the Enskog equation" का विस्तृत तकनीकी सारांश दिया गया है।

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

यह शोध पत्र दो समानांतर, स्थिर प्लेटों के बीच सीमित एक सघन गैस (dense gas) के तापीय विश्रांति (thermal relaxation) की प्रक्रिया की जांच करता है, जिन्हें एक समान, स्थिर तापमान ($T_0$) पर बनाए रखा गया है। प्रणाली एक गैर-संतुलन अवस्था (non-equilibrium state) से शुरू होती है, जिसकी विशेषता एक समान तापमान लेकिन एक गैर-समान घनत्व प्रोफाइल (विशेष रूप से, एक साइनुसोइडल भिन्नता) है।

मुख्य उद्देश्य इस गैर-संतुलन सघन गैस प्रणाली में हेल्महोल्ट्ज़ मुक्त ऊर्जा (Helmholtz free energy - $F$) के समय विकास का विश्लेषण करना है। यह अध्ययन एन्स्कोग समीकरण (Enskog equation) के दो रूपों की तुलना करता है, जो आणविक आयतन और परिमित-आकार के टकराव प्रभावों को ध्यान में रखते हुए घने गैसों के लिए बोल्ट्ज़मैन समीकरण का विस्तार करता है:

1. मूल एन्स्कोग समीकरण (Original Enskog Equation - OEE): यह शास्त्रीय एन्स्कोग टकराव कारक का उपयोग करता है।
2. संशोधित एन्स्कोग समीकरण (Enskog Equation with a Slight Modification - EESM): यह ताकाता और ताकाहाशी (2025) द्वारा प्रस्तावित एक संशोधित एन्स्कोग कारक का उपयोग करता है ताकि ऊष्मप्रवैगिकी स्थिरता (thermodynamic consistency) सुनिश्चित की जा सके।

केंद्रीय प्रश्न यह है कि क्या इस गैर-संतquilibrium प्रणाली की मुक्त ऊर्जा समय के साथ एकदिष्ट (monotonically) रूप से घटती है, जैसा कि ऊष्मा सिंक (heat bath) के संपर्क में रहने वाली प्रणालियों के लिए द्वितीय नियम (विशेष रूप से, H-प्रमेय) द्वारा आवश्यक है।

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

A. गणितीय सूत्रीकरण (Mathematical Formulation)

• शासी समीकरण (Governing Equation): गैस के व्यवहार को हार्ड-स्फीयर अणुओं के लिए एन्स्कोग समीकरण का उपयोग करके मॉडल किया गया है। कोलिजन इंटीग्रल (collision integral) में एक एन्स्कोग कारक ($g$) शामिल है, जो स्थानीय घनत्व और आणविक आकार के आधार पर टकराव की आवृत्ति को ठीक करता है।
• $g$ के लिए दो मॉडल:
  • OEE: $g(X_1, Y_1) = 2S(\frac{4\pi\sigma^3}{3m}\rho(\frac{X_1+Y_1}{2}))$, जहाँ $S$ कार्नाहन-स्टार्लिंग समीकरण (Carnahan-Starling equation) से व्युत्पन्न है।
  • EESM: $g(X_1, Y_1) = S(R(X_1)) + S(R(Y_1))$, जहाँ $R(X_1)$ घनत्व का एक गैर-स्थानीय समाकल (non-local integral) है। यह संशोधन सैद्धांतिक रूप से सिद्ध किया गया है कि H-प्रमेय का पालन करता है।
• सीमा शर्तें (Boundary Conditions): प्लेटों पर विसरित परावर्तन (diffuse reflection) (अणु $T_0$ पर मैक्सवेलियन वितरण के साथ पुनः उभरते हैं)।
• मुक्त ऊर्जा परिभाषा: अध्ययन गैर-संतुलन अवस्थाओं के लिए मुक्त ऊर्जा ($F$) की एक विस्तारित परिभाषा का उपयोग करता है:
  $$F = RT_0 H + E$$
  जहाँ $H = H^{(k)} + H^{(c)}$ एक सामान्यीकृत H-फलन है (गतिक भाग + कोलिजनल/सहसंबंध भाग) और $E$ कुल ऊर्जा है। EESM के लिए, यह सैद्धांतिक रूप से सिद्ध है कि $F$ एकदिष्ट रूप से घटता है ($dF/dt \leq 0$)।

B. संख्यात्मक दृष्टिकोण (Numerical Approach)

• विविक्तकरण (Discretization): लेखक एक हाइब्रिड संख्यात्मक विधि का उपयोग करते हैं:
  • स्थानिक अवकलज (Spatial Derivatives): फाइनाइट-डिफरेंस विधि (अपविंड स्कीम्स)।
  • कोलिजन इंटीग्रल: फास्ट फूरियर स्पेक्ट्रल मेथड (FFSM)। यह एन्स्कोग समीकरण के जटिल कोलिजन ऑपरेटर को संभालने में उच्च सटीकता प्रदान करता है।
• समय एकीकरण (Time Integration): कोलिजन टर्म के लिए दूसरे क्रम की सटीकता प्राप्त करने हेतु एक दूसरे क्रम के एक्सट्रपलेशन स्कीम (एडम्स-बाशफोर्थ) का उपयोग किया जाता है।
• विमाहीन पैरामीटर (Dimensionless Parameters): समस्या को चार विमाहीन मापदंडों द्वारा परिभाषित किया गया है:
  • $\hat{\sigma} = \sigma/L$ (अंतराल की चौड़ाई के सापेक्ष आणविक आकार)।
  • $\eta_0$ (संदर्भ आयतन अंश)।
  • $\lambda$ (प्रारंभिक घनत्व विक्षोभ की तरंग दैर्ध्य)।
  • $w$ (प्रारंभिक घनत्व विक्षोभ का आयाम)।
• द्रव्यमान संरक्षण (Mass Conservation): चूंकि संख्यात्मक योजना पूर्णतः संरक्षित (conservative) नहीं है, इसलिए कुल द्रव्यमान भिन्नता को नियंत्रित करने के लिए प्रत्येक समय चरण (time step) पर एक सामान्यीकरण चरण लागू किया जाता है।

3. मुख्य योगदान (Key Contributions)

1. EESM के लिए H-प्रमेय का संख्यात्मक सत्यापन: यह शोध पत्र पहला कठोर संख्यात्मक प्रमाण प्रदान करता है कि संशोधित एन्स्कोग समीकरण (EESM) विसरित परावर्तन सीमाओं वाले एक सीमित डोमेन में मुक्त ऊर्जा की एकदिष्ट कमी को संतुष्ट करता है।
2. OEE की विफलता का प्रदर्शन: अध्ययन संख्यात्मक रूप से प्रदर्शित करता है कि मूल एन्स्कोग समीकरण (OEE) मुक्त ऊर्जा की एकदिष्ट कमी की गारंटी देने में विफल रहता है। OEE सिमुलेशन में मुक्त ऊर्जा का व्यवहार गैर-एकदिष्ट (non-monotonic) होता है (ओवरशूट प्रदर्शित करता है), जो ऊष्मप्रवैगिकी अपेक्षा का उल्लंघन करता है।
3. गैर-गतिक योगदान की पहचान: विश्लेषण से पता चलता है कि EESM में एकदिष्टता काफी हद तक मुक्त ऊर्जा के गैर-गतिक (सहसंबंध) भाग ($H^{(c)}$) पर निर्भर करती है। "आदर्श गैस" वाला मुक्त ऊर्जा भाग ($F - RT_0 H^{(c)}$) दोनों मॉडलों में एकदिष्ट रूप से कम नहीं होता है, जो ऊष्मप्रवैगिकी स्थिरता के लिए गहन-गैस सहसंबंधों की महत्वपूर्ण भूमिका को रेखांकित करता है।
4. क्षणिक बनाम स्थिर अवस्था अंतर: अध्ययन दिखाता है कि जबकि OEE और EESM लगभग समान अंतिम स्थिर अवस्था घनत्व प्रोफाइल की ओर अभिसरित होते हैं, उनके क्षणिक समय विकास (घनत्व और तापमान प्रोफाइल) काफी भिन्न होते हैं, विशेष रूप से जब प्रारंभिक घनत्व विक्षोभ बड़ा होता है।

4. मुख्य परिणाम (Key Results)

• मुक्त ऊर्जा विकास (Free Energy Evolution):

  • EESM: कुल मुक्त ऊर्जा $F$ समय के साथ सख्ती से एकदिष्ट रूप से घटती है, जो सैद्धांतिक H-प्रमेय की पुष्टि करती है।
  • OEE: $F$ शुरू में घटता है लेकिन अंततः स्थिर होने से पहले गैर-एकदिष्ट व्यवहार (ओवरशूट) प्रदर्शित करता है। यह इंगित करता है कि OEE इस गैर-संतुलन संदर्भ में दूसरे नियम का सख्ती से पालन नहीं करता है।
  • विक्षोभ पर निर्भरता: OEE और EESM के बीच विचलन प्रारंभिक घनत्व भिन्नता के आयाम ($w$) के बढ़ने के साथ अधिक स्पष्ट हो जाता है।

• घनत्व और तापमान प्रोफाइल:

  • स्थिर अवस्था (Steady State): दोनों मॉडल संतुलन पर समान गैर-समान घनत्व प्रोफाइल की भविष्यवाणी करते हैं (अपवर्जित आयतन प्रभावों के कारण दीवारों के पास उच्च घनत्व), जो साहित्य के अनुरूप है।
  • क्षणिक चरण (Transient Phase): घनत्व और तापमान के समय विकास में महत्वपूर्ण अंतर देखे गए हैं। उदाहरण के लिए, विशिष्ट समय चरणों पर (जैसे, $\hat{t} \approx 0.08$), OEE और EESM के बीच तापमान प्रोफाइल स्पष्ट विचलन दिखाते हैं, भले ही अंतिम अवस्थाएँ समान हों।

• आणविक आकार के प्रति संवेदनशीलता: OEE में "ओवरशूट" की घटना विशिष्ट समय पैमानों पर होती है, जो आणविक टकराव समय ($t_\sigma^0$) के सापेक्ष सुसंगत है, चाहे अंतराल की चौड़ाई के सापेक्ष आणविक आकार का अनुपात ($\sigma/L$) कुछ भी हो।

5. महत्व (Significance)

• गतिक सिद्धांत में ऊष्मप्रवैगिकी स्थिरता: यह कार्य मूल एन्स्कोग समीकरण की ऊष्मप्रवैगिकी स्थिरता के संबंध में लंबे समय से चले आ रहे मुद्दे को हल करता है। यह पुष्टि करता है कि हालांकि OEE व्यापक रूप से उपयोग किया जाता है, इसमें सीमित, गैर-संतुलन प्रणालियों में कठोर H-प्रमेय की गारंटी नहीं है।
• संशोधित मॉडलों का सत्यापन: ये परिणाम सघन गैसों के संख्यात्मक सिमुलेशन के लिए EESM (संशोधित एन्स्कोग समीकरण) को अपनाने का पुरजोर समर्थन करते हैं। EESM यह सुनिश्चित करता है कि संख्यात्मक समाधान मौलिक ऊष्मप्रवैगिकी नियमों (एकदिष्ट एन्ट्रॉपी/मुक्त ऊर्जा उत्पादन) का सम्मान करते हैं।
• माइक्रो/नैनो-स्केल प्रवाह के निहितार्थ: जैसे-जैसे माइक्रोफैब्रिकेशन तकनीक आगे बढ़ रही है, उप-माइक्रो स्केल में सघन गैस व्यवहार को समझना महत्वपूर्ण है। यह शोध पत्र एक सत्यापित ढांचा (EESM) प्रदान करता है, जो हीट ट्रांसफर और रिलैक्सेशन को सिम्युलेट करने के लिए, जिससे यह सुनिश्चित होता है कि सिमुलेशन अन भौतिक परिणाम (जैसे मुक्त ऊर्जा में वृद्धि) उत्पन्न न करें जो OEE के उपयोग से उत्पन्न हो सकते हैं।
• पद्धतिगत प्रगति: संशोधित एन्स्कोग समीकरण के लिए फास्ट फूरियर स्पेक्ट्रल मेथड का सफल अनुप्रयोग जटिल कोलिजन इंटीग्रल को संभालने के लिए एक मजबूत कम्प्यूटेशनल दृष्टिकोण का प्रदर्शन करता है।

निष्कर्षतः, यह शोध पत्र स्थापित करता है कि सघन गैस विश्रांति के अनुकरण के लिए संशोधित कारक वाला एन्स्कोग समीकरण (EESM) ऊष्मप्रवैगिकी रूप से सुसंगत विकल्प है, जबकि मूल एन्स्कोग समीकरण (OEE) सही स्थिर अवस्था की ओर अभिसरित होने के बावजूद, भौतिक रूप से असंगत क्षणिक व्यवहार उत्पन्न कर सकता है।