2025-11-20T16:49:15.210607

Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model

Li, Gibelli, Zhang
Enskog-Vlasov equation is currently the most sophisticated kinetic model for describing non-equilibrium evaporative flows. While it enables more efficient simulations than the molecular dynamics (MD) methods, its accuracy in reproducing the flow properties of real fluids is limited by both the assumptions underlying the Vlasov forcing term and the approximation introduced by the Enskog collision term for short-range molecular interactions. To address this limitation, this work proposes a molecular kinetic model specifically designed for real fluids, with the Lennard-Jones fluids as an example. The model is first applied to evaluate the equilibrium characteristics of a liquid-vapour system, including the liquid-vapour coexistence curve, transport coefficients, vapour pressure, and surface tension coefficient. The results show excellent agreement with the MD simulation and experimental data. Furthermore, the model is used to investigate non-equilibrium evaporation, with a particular focus on the velocity distribution function adjacent to the liquid-vapour interface. The results confirm that deviations from the Maxwellian distribution persist in the vapour region, indicating limitations of the classical Hertz-Knudsen relation under pronounced non-equilibrium conditions. This work represents a critical step towards the development of an accurate and efficient computational framework for modelling non-equilibrium liquid-vapour flows for real fluids, with direct relevance to practical applications such as flow cooling.
academic

लेनार्ड-जोन्स तरल पदार्थों का गैर-संतुलन वाष्पीकरण: एन्स्कोग-व्लासोव सिद्धांत और हर्ट्ज-क्नुडसन मॉडल

मूल जानकारी

  • पेपर ID: 2508.15785
  • शीर्षक: Non-equilibrium evaporation of Lennard-Jones fluids: Enskog-Vlasov theory and Hertz-Knudsen model
  • लेखक: शाओकांग ली, लिविओ गिबेली, योंगहाओ झांग
  • वर्गीकरण: cond-mat.soft physics.flu-dyn
  • प्रकाशन तिथि: 14 अक्टूबर 2025 (arXiv v2)
  • पेपर लिंक: https://arxiv.org/abs/2508.15785

सारांश

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

अनुसंधान पृष्ठभूमि और प्रेरणा

1. अनुसंधान समस्या

इस अनुसंधान द्वारा समाधान की जाने वाली मूल समस्या यह है कि नैनो-स्केल वाष्पीकरण प्रक्रिया में, मौजूदा गतिशील मॉडल वास्तविक तरल पदार्थ (जैसे लेनार्ड-जोन्स तरल पदार्थ) के गैर-संतुलन वाष्पीकरण व्यवहार का सटीक वर्णन नहीं कर सकते।

2. समस्या की महत्ता

नैनो-स्केल वाष्पीकरण कई अग्रणी तकनीकी क्षेत्रों में महत्वपूर्ण भूमिका निभाता है:

  • नैनो-छिद्रपूर्ण झिल्लियों में पदार्थ पृथक्करण दक्षता में वृद्धि
  • उच्च-प्रदर्शन इलेक्ट्रॉनिक उपकरणों का कुशल ताप निष्कासन
  • प्रवाह शीतलन जैसे व्यावहारिक अनुप्रयोग

3. मौजूदा विधियों की सीमाएं

  • हर्ट्ज-क्नुडसन संबंध: गैस चरण क्षेत्र में स्थानीय संतुलन मानता है, क्नुडसन परत में गतिशील संरचना को नजरअंदाज करता है
  • एन्स्कोग-व्लासोव समीकरण: सदरलैंड संभावित ऊर्जा मॉडल का उपयोग करता है, जो वास्तविक लेनार्ड-जोन्स संभावित ऊर्जा से महत्वपूर्ण रूप से भिन्न है
  • आणविक गतिविज्ञान विधि: अत्यधिक कम्प्यूटेशनल लागत, इंजीनियरिंग अनुप्रयोगों में लागू करना कठिन

4. अनुसंधान प्रेरणा

एक एकीकृत ढांचा विकसित करना जो:

  • तरल चरण, गैस चरण और उनके इंटरफेस का सटीक वर्णन कर सके
  • मनमाने वाष्पीकरण गुणांक पर निर्भरता को समाप्त करे
  • कम्प्यूटेशनल दक्षता बनाए रखते हुए पूर्वानुमान सटीकता में सुधार करे

मुख्य योगदान

  1. सुधारे गए एन्स्कोग-व्लासोव गतिशील मॉडल का प्रस्ताव: लेनार्ड-जोन्स तरल पदार्थों के लिए विशेष रूप से डिज़ाइन किया गया, संबंधित कार्यों को संशोधित करके वास्तविक तरल पदार्थ के स्थिति समीकरण को सटीक रूप से पुनः प्रस्तुत करता है
  2. सटीक थर्मोडायनामिक पैरामीटर कैलिब्रेशन विधि की स्थापना: महत्वपूर्ण बिंदु और त्रिपल बिंदु मापदंडों को एक साथ मिलाकर, मॉडल और वास्तविक तरल पदार्थ गुणों के बीच संगति सुनिश्चित करता है
  3. मॉडल की संतुलन गुण पूर्वानुमान क्षमता का सत्यापन: तरल-गैस सह-अस्तित्व वक्र, परिवहन गुणांक, वाष्प दबाव और सतह तनाव के संदर्भ में प्रायोगिक डेटा के साथ उच्च समझौता
  4. गैर-संतुलन वाष्पीकरण में वेग वितरण विशेषताओं का प्रकटीकरण: गैस चरण क्षेत्र में वेग वितरण के मैक्सवेल वितरण से विचलन की पुष्टि करता है, हर्ट्ज-क्नुडसन संबंध की सीमा के लिए सैद्धांतिक आधार प्रदान करता है

विधि विवरण

कार्य परिभाषा

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

मॉडल आर्किटेक्चर

1. मूल गतिशील समीकरण

मॉडल संशोधित एन्स्कोग-व्लासोव समीकरण पर आधारित है:

∂f∂t+ξ⋅∂f∂x+Fm⋅∂f∂ξ=Ω\frac{\partial f}{\partial t} + \boldsymbol{\xi} \cdot \frac{\partial f}{\partial \boldsymbol{x}} + \frac{\boldsymbol{F}}{m} \cdot \frac{\partial f}{\partial \boldsymbol{\xi}} = \Omega

जहां:

  • f(x,ξ,t)f(\boldsymbol{x}, \boldsymbol{\xi}, t): आणविक वेग वितरण फ़ंक्शन
  • F\boldsymbol{F}: सदरलैंड संभावित ऊर्जा द्वारा उत्पन्न स्व-संगत बल क्षेत्र
  • Ω\Omega: सरलीकृत टकराव पद

2. सरलीकृत टकराव पद

टकराव पद टेलर विस्तार सन्निकटन को अपनाता है: Ω=JS+Je\Omega = J_S + J_e

जहां JSJ_S शाखोव मॉडल पद है, JeJ_e अतिरिक्त पद है, सही प्रांडटल संख्या और आयतन श्यानता सुनिश्चित करता है।

3. स्थिति समीकरण संशोधन

सामान्यीकृत वैन डेर वाल्स रूप अपनाता है: p=nkBT(1+ρbχ)−an2p = nk_BT(1 + \rho b\chi) - an^2

लेनार्ड-जोन्स तरल पदार्थ के थर्मोडायनामिक गुणों से मेल खाने के लिए संबंधित कार्य χ\chi को संशोधित करके।

तकनीकी नवाचार बिंदु

1. संबंधित कार्य का अनुभवजन्य संशोधन

अनुभवजन्य फ़ंक्शन Φ\Phi का परिचय: Φ=ρb+a1(ρb)2+a2(ρb)3+a3(ρb)4+a4(ρb)5\Phi = \rho b + a_1(\rho b)^2 + a_2(\rho b)^3 + a_3(\rho b)^4 + a_4(\rho b)^5

जहां χ=Φ′=dΦ/d(ρb)\chi = \Phi' = d\Phi/d(\rho b)

2. बहु-बाधा शर्त कैलिब्रेशन

एक साथ संतुष्ट करता है:

  • महत्वपूर्ण बिंदु शर्त: (∂p∂n)T=Tc,n=nc=0\left(\frac{\partial p}{\partial n}\right)_{T=T_c,n=n_c} = 0, (∂2p∂n2)T=Tc,n=nc=0\left(\frac{\partial^2 p}{\partial n^2}\right)_{T=T_c,n=n_c} = 0
  • त्रिपल बिंदु शर्त: दबाव और रासायनिक क्षमता समान

3. परिवहन गुणांक संशोधन

घनत्व प्रभाव पर विचार करते हुए श्यानता और ताप चालकता: μ=μ∗χ(1+25ρbχ)2+35ϖ\mu = \frac{\mu^*}{\chi}\left(1 + \frac{2}{5}\rho b\chi\right)^2 + \frac{3}{5}\varpiκ=κ∗χ(1+35ρbχ)2+cvϖ\kappa = \frac{\kappa^*}{\chi}\left(1 + \frac{3}{5}\rho b\chi\right)^2 + c_v\varpi

प्रायोगिक सेटअप

डेटासेट

  • संदर्भ तरल पदार्थ: आर्गन (Argon) लेनार्ड-जोन्स तरल पदार्थ के विशिष्ट प्रतिनिधि के रूप में
  • प्रायोगिक डेटा: तरल-गैस सह-अस्तित्व वक्र, परिवहन गुणांक, वाष्प दबाव, सतह तनाव के प्रायोगिक माप
  • MD सिमुलेशन डेटा: तुलनात्मक सत्यापन के लिए

मूल्यांकन संकेतक

  1. संतुलन गुण:
    • तरल-गैस सह-अस्तित्व घनत्व
    • कतरनी श्यानता और ताप चालकता
    • संतृप्त वाष्प दबाव
    • सतह तनाव गुणांक
  2. गैर-संतुलन विशेषताएं:
    • वेग वितरण फ़ंक्शन के मैक्सवेल वितरण से विचलन की डिग्री
    • क्नुडसन परत संरचना

तुलनात्मक विधि

  • मूल एन्स्कोग-व्लासोव समीकरण
  • प्रायोगिक डेटा
  • आणविक गतिविज्ञान सिमुलेशन परिणाम

कार्यान्वयन विवरण

  • कम्प्यूटेशनल डोमेन लंबाई: LL, आवधिक सीमा शर्तों (संतुलन) या अवशोषण सीमा शर्तों (गैर-संतुलन) का उपयोग करते हुए
  • तापमान रेंज: त्रिपल बिंदु से महत्वपूर्ण बिंदु तक
  • आर्गन पैरामीटर: a=162.5 J\cdotpl\cdotpmol−2a = 162.5 \text{ J·l·mol}^{-2}, nc=13.41 mol\cdotpl−1n_c = 13.41 \text{ mol·l}^{-1}, Tc=150.69 KT_c = 150.69 \text{ K}

प्रायोगिक परिणाम

मुख्य परिणाम

1. संतुलन गुण सत्यापन

  • तरल-गैस सह-अस्तित्व वक्र: संशोधित मॉडल पूरे तापमान रेंज में प्रायोगिक डेटा के साथ उच्च समझौता दिखाता है, जबकि मूल EV समीकरण कम तापमान पर महत्वपूर्ण विचलन दिखाता है
  • परिवहन गुणांक: कतरनी श्यानता प्रायोगिक डेटा के साथ पूरी तरह से सहमत है, ताप चालकता औसत त्रुटि लगभग 10%
  • वाष्प दबाव: व्यापक तापमान रेंज में प्रायोगिक मानों के साथ अच्छा समझौता
  • सतह तनाव: औसत सापेक्ष त्रुटि लगभग 9%

2. असंतृप्त तरल गुण

संतृप्ति वक्र से दूर की स्थितियों में:

  • उच्च तापमान पर परिवहन गुणांक पूर्वानुमान सटीक है
  • कम तापमान उच्च घनत्व पर विचलन मौजूद है (अपेक्षित है, क्योंकि मॉडल मुख्य रूप से संतृप्ति गुणों के लिए कैलिब्रेट किया गया है)

गैर-संतुलन वाष्पीकरण विश्लेषण

1. मैक्रोस्कोपिक पैरामीटर वितरण

  • तापमान केंद्र से सीमा की ओर घटता है
  • संख्या घनत्व शीतलन प्रभाव के कारण थोड़ा बढ़ता है
  • शरीर वेग इंटरफेस और गैस चरण क्षेत्र में बढ़ता है

2. वेग वितरण फ़ंक्शन विशेषताएं

विभिन्न स्थानों (x=11,13,14.5,15,17x = 11, 13, 14.5, 15, 17) पर विश्लेषण दिखाता है:

  • तरल चरण और इंटरफेस क्षेत्र: वेग वितरण मैक्सवेल वितरण के करीब है
  • गैस चरण क्षेत्र: नकारात्मक वेग रेंज में स्पष्ट विचलन, आणविक टकराव की कमी के कारण
  • सीमा के पास क्षेत्र: शरीर गति प्रभाव महत्वपूर्ण है, वितरण में स्पष्ट विस्थापन

प्रायोगिक निष्कर्ष

  1. हर्ट्ज-क्नुडसन संबंध की सीमा: गैस चरण क्षेत्र में वेग वितरण की गैर-मैक्सवेल विशेषता मजबूत गैर-संतुलन स्थितियों के तहत HK संबंध की अप्रयोज्यता की पुष्टि करती है
  2. दो महत्वपूर्ण कारक:
    • तरल-गैस इंटरफेस के पास गैस चरण क्षेत्र में वेग वितरण का मैक्सवेल वितरण से विचलन
    • गैस चरण पक्ष पर शरीर गति प्रभाव क्रमशः महत्वपूर्ण हो जाता है

संबंधित कार्य

मुख्य अनुसंधान दिशाएं

  1. शास्त्रीय सिद्धांत: हर्ट्ज-क्नुडसन संबंध, श्रेज समीकरण
  2. गतिशील विधि: बोल्ट्जमैन समीकरण, DSMC विधि, एन्स्कोग-व्लासोव सिद्धांत
  3. आणविक सिमुलेशन: आणविक गतिविज्ञान, मोंटे कार्लो विधि

इस पेपर के लाभ

  • शास्त्रीय सिद्धांत की तुलना में: क्नुडसन परत की गतिशील संरचना पर विचार करता है
  • मूल EV समीकरण की तुलना में: वास्तविक तरल पदार्थ गुणों का अधिक सटीक वर्णन
  • MD विधि की तुलना में: कम्प्यूटेशनल दक्षता में महत्वपूर्ण सुधार

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. लेनार्ड-जोन्स तरल पदार्थों के लिए सुधारे गए गतिशील मॉडल का सफल विकास
  2. मॉडल संतुलन गुण पूर्वानुमान में प्रायोगिक डेटा के साथ उच्च समझौता
  3. गैर-संतुलन वाष्पीकरण में वेग वितरण की गैर-मैक्सवेल विशेषता का प्रकटीकरण
  4. हर्ट्ज-क्नुडसन संबंध की सीमा के लिए सैद्धांतिक आधार प्रदान करता है

सीमाएं

  1. प्रयोज्यता रेंज: मुख्य रूप से संतृप्ति गुणों के लिए अनुकूलित, गैर-संतृप्ति स्थितियों में सटीकता सीमित
  2. तरल पदार्थ प्रकार: वर्तमान में केवल निष्क्रिय गैसों के लिए लागू, बहु-परमाणु अणुओं तक विस्तार के लिए आंतरिक ऊर्जा मोड पर विचार करने की आवश्यकता है
  3. ज्यामितीय जटिलता: मुख्य रूप से समतल इंटरफेस स्थिति का सत्यापन किया गया है

भविष्य की दिशाएं

  1. अन्य निष्क्रिय गैसों तक विस्तार
  2. जटिल बहु-परमाणु तरल पदार्थों के लिए उपयुक्त मॉडल विकास
  3. जटिल ज्यामितीय स्थितियों के तहत अनुप्रयोग अनुसंधान

गहन मूल्यांकन

लाभ

  1. विधि नवाचार: एन्स्कोग सिद्धांत और व्लासोव सन्निकटन को चतुराई से संयोजित करता है, अनुभवजन्य संशोधन के माध्यम से वास्तविक तरल पदार्थ का सटीक वर्णन प्राप्त करता है
  2. सत्यापन पूर्णता: संतुलन गुणों और गैर-संतुलन विशेषताओं सहित कई कोणों से मॉडल की सटीकता का सत्यापन करता है
  3. सैद्धांतिक योगदान: गैर-संतुलन वाष्पीकरण को समझने के लिए महत्वपूर्ण अंतर्दृष्टि प्रदान करता है, विशेष रूप से HK संबंध की सीमा के प्रकटीकरण में
  4. व्यावहारिक मूल्य: कम्प्यूटेशनल दक्षता बनाए रखते हुए पूर्वानुमान सटीकता में महत्वपूर्ण सुधार

कमियां

  1. पैरामीटर निर्भरता: मॉडल में कई अनुभवजन्य पैरामीटर हैं, जो इसकी सामान्यीकरण क्षमता को प्रभावित कर सकते हैं
  2. सैद्धांतिक आधार: कुछ संशोधनों में कठोर सैद्धांतिक व्युत्पत्ति का अभाव है
  3. सत्यापन रेंज: मुख्य रूप से आर्गन के आधार पर सत्यापित, अन्य तरल पदार्थों की प्रयोज्यता को आगे की पुष्टि की आवश्यकता है

प्रभाव

  1. शैक्षणिक योगदान: गैर-संतुलन तरल पदार्थ यांत्रिकी के लिए नई मॉडलिंग सोच प्रदान करता है
  2. अनुप्रयोग संभावनाएं: नैनो तकनीक, ताप प्रबंधन आदि क्षेत्रों में महत्वपूर्ण अनुप्रयोग मूल्य
  3. पुनरुत्पादनीयता: विधि विवरण स्पष्ट है, अच्छी पुनरुत्पादनीयता है

प्रयोज्य परिदृश्य

  • नैनो-स्केल वाष्पीकरण शीतलन प्रणाली डिजाइन
  • छिद्रपूर्ण माध्यम में चरण परिवर्तन ताप स्थानांतरण
  • सूक्ष्म-नैनो उपकरणों का ताप प्रबंधन
  • आसवन और पृथक्करण प्रक्रिया अनुकूलन

संदर्भ

पेपर 41 महत्वपूर्ण संदर्भों का हवाला देता है, जो गतिशील सिद्धांत, प्रायोगिक डेटा और संख्यात्मक विधियों सहित कई पहलुओं को कवर करते हैं, अनुसंधान के लिए एक मजबूत सैद्धांतिक आधार प्रदान करते हैं।