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
লেখক: Shaokang Li, Livio Gibelli, Yonghao Zhang
শ্রেণিবিভাগ: cond-mat.soft physics.flu-dyn
প্রকাশকাল: ১৪ অক্টোবর ২০২৫ (arXiv v2)
গবেষণাপত্র লিংক: https://arxiv.org/abs/2508.15785

সারসংক্ষেপ

এনস্কগ-ভ্লাসভ সমীকরণ বর্তমানে অ-সাম্যাবস্থা বাষ্পীভবন প্রবাহ বর্ণনার সবচেয়ে নিখুঁত গতিবিদ্যা মডেল। যদিও এটি আণবিক গতিবিদ্যা (MD) পদ্ধতির চেয়ে অধিক দক্ষ, তবে প্রকৃত তরল প্রবাহের বৈশিষ্ট্য পুনরুৎপাদনে এর নির্ভুলতা ভ্লাসভ বলপ্রয়োগ পদের অনুমান এবং স্বল্প-পাল্লার আণবিক মিথস্ক্রিয়ার এনস্কগ সংঘর্ষ পদের আনুমানিকতার দ্বারা সীমাবদ্ধ। এই সীমাবদ্ধতা সমাধানের জন্য, এই গবেষণা প্রকৃত তরলের জন্য বিশেষভাবে ডিজাইন করা আণবিক গতিবিদ্যা মডেল প্রস্তাব করে, লেনার্ড-জোন্স তরলকে উদাহরণ হিসেবে নিয়ে। এই মডেলটি প্রথমে তরল-বায়বীয় সিস্টেমের সাম্যাবস্থা বৈশিষ্ট্য মূল্যায়নে ব্যবহৃত হয়, যার মধ্যে রয়েছে তরল-বায়বীয় সহাবস্থান বক্ররেখা, পরিবহন সহগ, বাষ্পচাপ এবং পৃষ্ঠটান সহগ। ফলাফল MD সিমুলেশন এবং পরীক্ষামূলক তথ্যের সাথে চমৎকার সামঞ্জস্য প্রদর্শন করে। উপরন্তু, এই মডেলটি অ-সাম্যাবস্থা বাষ্পীভবন অধ্যয়নে ব্যবহৃত হয়, বিশেষত তরল-বায়বীয় ইন্টারফেসের নিকটবর্তী বেগ বিতরণ ফাংশনের উপর মনোনিবেশ করে। ফলাফল নিশ্চিত করে যে বায়বীয় অঞ্চলে ম্যাক্সওয়েল বিতরণ থেকে বিচ্যুতি ক্রমাগত বিদ্যমান, যা উল্লেখযোগ্য অ-সাম্যাবস্থা পরিস্থিতিতে ক্লাসিক্যাল হার্টজ-নুডসেন সম্পর্কের সীমাবদ্ধতা নির্দেশ করে।

গবেষণা পটভূমি ও প্রেরণা

১. গবেষণা সমস্যা

এই গবেষণায় সমাধানযোগ্য মূল সমস্যা হলো ন্যানো-স্কেল বাষ্পীভবন প্রক্রিয়ায়, বিদ্যমান গতিবিদ্যা মডেলগুলি প্রকৃত তরলের (যেমন লেনার্ড-জোন্স তরল) অ-সাম্যাবস্থা বাষ্পীভবন আচরণ সঠিকভাবে বর্ণনা করতে অক্ষম।

২. সমস্যার গুরুত্ব

ন্যানো-স্কেল বাষ্পীভবন একাধিক অত্যাধুনিক প্রযুক্তি ক্ষেত্রে গুরুত্বপূর্ণ ভূমিকা পালন করে:

ন্যানো-ছিদ্রযুক্ত ঝিল্লিতে পদার্থ পৃথকীকরণ দক্ষতা বৃদ্ধি
উচ্চ-কর্মক্ষমতা ইলেকট্রনিক ডিভাইসের দক্ষ তাপ অপসারণ
প্রবাহ শীতলীকরণ ইত্যাদি ব্যবহারিক প্রয়োগ

৩. বিদ্যমান পদ্ধতির সীমাবদ্ধতা

হার্টজ-নুডসেন সম্পর্ক: বায়বীয় অঞ্চলে স্থানীয় সাম্যাবস্থা অনুমান করে, নুডসেন স্তরের গতিবিদ্যা কাঠামো উপেক্ষা করে
এনস্কগ-ভ্লাসভ সমীকরণ: সাদারল্যান্ড বিভব মডেল ব্যবহার করে, যা প্রকৃত লেনার্ড-জোন্স বিভবের সাথে উল্লেখযোগ্য পার্থক্য রয়েছে
আণবিক গতিবিদ্যা পদ্ধতি: গণনা ব্যয় অত্যধিক, প্রকৌশল বাস্তবায়নে প্রয়োগ করা কঠিন

৪. গবেষণা প্রেরণা

একটি একীভূত কাঠামো উন্নয়ন, যা সক্ষম:

তরল, বায়বীয় এবং তাদের ইন্টারফেস সঠিকভাবে বর্ণনা করতে
স্বেচ্ছাচারী বাষ্পীভবন সহগের উপর নির্ভরতা দূর করতে
গণনা দক্ষতা বজায় রেখে পূর্বাভাস নির্ভুলতা বৃদ্ধি করতে

মূল অবদান

১. উন্নত এনস্কগ-ভ্লাসভ গতিবিদ্যা মডেল প্রস্তাব: লেনার্ড-জোন্স তরলের জন্য বিশেষভাবে ডিজাইন করা, যুগল পারস্পরিক সম্পর্ক ফাংশন সংশোধনের মাধ্যমে প্রকৃত তরলের অবস্থা সমীকরণ সঠিকভাবে পুনরুৎপাদন করে

২. সুনির্দিষ্ট তাপগতিবিদ্যা প্যারামিটার ক্রমাঙ্কন পদ্ধতি প্রতিষ্ঠা: সমালোচনামূলক বিন্দু এবং ত্রি-দশা বিন্দু প্যারামিটার একযোগে মিলিয়ে, মডেল এবং প্রকৃত তরল বৈশিষ্ট্যের সামঞ্জস্য নিশ্চিত করে

৩. মডেলের সাম্যাবস্থা বৈশিষ্ট্য পূর্বাভাস ক্ষমতা যাচাই: তরল-বায়বীয় সহাবস্থান বক্ররেখা, পরিবহন সহগ, বাষ্পচাপ এবং পৃষ্ঠটান ইত্যাদিতে পরীক্ষামূলক তথ্যের সাথে উচ্চ সামঞ্জস্য

৪. অ-সাম্যাবস্থা বাষ্পীভবনে বেগ বিতরণ বৈশিষ্ট্য উন্মোচন: নিশ্চিত করে যে বায়বীয় অঞ্চলের বেগ বিতরণ ম্যাক্সওয়েল বিতরণ থেকে বিচ্যুত, হার্টজ-নুডসেন সম্পর্কের সীমাবদ্ধতার জন্য তাত্ত্বিক ভিত্তি প্রদান করে

পদ্ধতি বিস্তারিত

কার্য সংজ্ঞা

অ-সাম্যাবস্থা পরিস্থিতিতে লেনার্ড-জোন্স তরলের বাষ্পীভবন প্রক্রিয়া অধ্যয়ন, তরল-বায়বীয় ইন্টারফেসের নিকটবর্তী আণবিক বেগ বিতরণ বৈশিষ্ট্য বিশ্লেষণে মনোনিবেশ, এবং ক্লাসিক্যাল হার্টজ-নুডসেন সম্পর্কের প্রযোজ্যতা মূল্যায়ন।

মডেল কাঠামো

১. মৌলিক গতিবিদ্যা সমীকরণ

মডেল সংশোধিত এনস্কগ-ভ্লাসভ সমীকরণের উপর ভিত্তি করে:

$\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(\boldsymbol{x}, \boldsymbol{\xi}, t)$ : আণবিক বেগ বিতরণ ফাংশন
$\boldsymbol{F}$ : সাদারল্যান্ড বিভব দ্বারা উৎপন্ন স্ব-সামঞ্জস্যপূর্ণ বল ক্ষেত্র
$\Omega$ : সরলীকৃত সংঘর্ষ পদ

২. সরলীকৃত সংঘর্ষ পদ

সংঘর্ষ পদ টেইলর সম্প্রসারণ আনুমানিকতা গ্রহণ করে: $\Omega = J_S + J_e$

যেখানে $J_S$ শাখভ মডেল পদ, $J_e$ অতিরিক্ত পদ, সঠিক প্র্যান্ডটল সংখ্যা এবং আয়তন সান্দ্রতা নিশ্চিত করে।

৩. অবস্থা সমীকরণ সংশোধন

সাধারণীকৃত ভ্যান ডার ওয়ালস রূপ গ্রহণ: $p = nk_BT(1 + \rho b\chi) - an^2$

যুগল পারস্পরিক সম্পর্ক ফাংশন $\chi$ সংশোধনের মাধ্যমে লেনার্ড-জোন্স তরলের তাপগতিবিদ্যা বৈশিষ্ট্য মিলানো হয়।

প্রযুক্তিগত উদ্ভাবন বিন্দু

১. যুগল পারস্পরিক সম্পর্ক ফাংশনের অভিজ্ঞতামূলক সংশোধন

অভিজ্ঞতামূলক ফাংশন $\Phi$ প্রবর্তন: $\Phi = \rho b + a_1(\rho b)^2 + a_2(\rho b)^3 + a_3(\rho b)^4 + a_4(\rho b)^5$

যেখানে $\chi = \Phi' = d\Phi/d(\rho b)$

২. বহু-সীমাবদ্ধতা শর্ত ক্রমাঙ্কন

একযোগে সন্তুষ্ট করে:

সমালোচনামূলক বিন্দু শর্ত: $\left(\frac{\partial p}{\partial n}\right)_{T=T_c,n=n_c} = 0$ , $\left(\frac{\partial^2 p}{\partial n^2}\right)_{T=T_c,n=n_c} = 0$
ত্রি-দশা বিন্দু শর্ত: চাপ এবং রাসায়নিক বিভব সমান

৩. পরিবহন সহগ সংশোধন

ঘনত্ব প্রভাব বিবেচনা করে সান্দ্রতা এবং তাপ পরিবাহিতা: $\mu = \frac{\mu^*}{\chi}\left(1 + \frac{2}{5}\rho b\chi\right)^2 + \frac{3}{5}\varpi$ $\kappa = \frac{\kappa^*}{\chi}\left(1 + \frac{3}{5}\rho b\chi\right)^2 + c_v\varpi$

পরীক্ষামূলক সেটআপ

ডেটাসেট

রেফারেন্স তরল: আর্গন (Argon) লেনার্ড-জোন্স তরলের প্রতিনিধি উদাহরণ হিসেবে
পরীক্ষামূলক তথ্য: তরল-বায়বীয় সহাবস্থান বক্ররেখা, পরিবহন সহগ, বাষ্পচাপ, পৃষ্ঠটানের পরীক্ষামূলক পরিমাপ মান
MD সিমুলেশন তথ্য: তুলনামূলক যাচাইয়ের জন্য ব্যবহৃত

মূল্যায়ন সূচক

১. সাম্যাবস্থা বৈশিষ্ট্য:

তরল-বায়বীয় সহাবস্থান ঘনত্ব
শিয়ার সান্দ্রতা এবং তাপ পরিবাহিতা
সম্পৃক্ত বাষ্পচাপ
পৃষ্ঠটান সহগ

২. অ-সাম্যাবস্থা বৈশিষ্ট্য:

বেগ বিতরণ ফাংশন ম্যাক্সওয়েল বিতরণ থেকে বিচ্যুতির মাত্রা
নুডসেন স্তর কাঠামো

তুলনামূলক পদ্ধতি

মূল এনস্কগ-ভ্লাসভ সমীকরণ
পরীক্ষামূলক তথ্য
আণবিক গতিবিদ্যা সিমুলেশন ফলাফল

বাস্তবায়ন বিস্তারিত

গণনা ডোমেইন দৈর্ঘ্য: $L$ , পর্যায়ক্রমিক সীমানা শর্ত (সাম্যাবস্থা) বা শোষণকারী সীমানা শর্ত (অ-সাম্যাবস্থা) গ্রহণ
তাপমাত্রা পরিসীমা: ত্রি-দশা বিন্দু থেকে সমালোচনামূলক বিন্দু পর্যন্ত
আর্গন প্যারামিটার: $a = 162.5 \text{ J·l·mol}^{-2}$ , $n_c = 13.41 \text{ mol·l}^{-1}$ , $T_c = 150.69 \text{ K}$

পরীক্ষামূলক ফলাফল

প্রধান ফলাফল

১. সাম্যাবস্থা বৈশিষ্ট্য যাচাই

তরল-বায়বীয় সহাবস্থান বক্ররেখা: সংশোধিত মডেল সম্পূর্ণ তাপমাত্রা পরিসীমায় পরীক্ষামূলক তথ্যের সাথে উচ্চ সামঞ্জস্য প্রদর্শন করে, যেখানে মূল EV সমীকরণ নিম্ন তাপমাত্রায় উল্লেখযোগ্য বিচ্যুতি দেখায়
পরিবহন সহগ: শিয়ার সান্দ্রতা পরীক্ষামূলক তথ্যের সাথে সম্পূর্ণ সামঞ্জস্যপূর্ণ, তাপ পরিবাহিতা গড় ত্রুটি প্রায় ১০%
বাষ্পচাপ: বিস্তৃত তাপমাত্রা পরিসীমায় পরীক্ষামূলক মানের সাথে ভাল সামঞ্জস্য
পৃষ্ঠটান: গড় আপেক্ষিক ত্রুটি প্রায় ৯%

২. অ-সম্পৃক্ত তরল বৈশিষ্ট্য

সম্পৃক্ত বক্ররেখা থেকে দূরবর্তী পরিস্থিতিতে:

উচ্চ তাপমাত্রায় পরিবহন সহগ পূর্বাভাস সঠিক
নিম্ন তাপমাত্রা উচ্চ ঘনত্বে বিচ্যুতি বিদ্যমান (প্রত্যাশিত, কারণ মডেল প্রধানত সম্পৃক্ত বৈশিষ্ট্য ক্রমাঙ্কনের জন্য লক্ষ্যবদ্ধ)

অ-সাম্যাবস্থা বাষ্পীভবন বিশ্লেষণ

১. সামগ্রিক প্যারামিটার বিতরণ

তাপমাত্রা কেন্দ্র থেকে সীমানার দিকে হ্রাস পায়
সংখ্যা ঘনত্ব শীতলীকরণ প্রভাবের কারণে সামান্য বৃদ্ধি পায়
আয়তন বেগ ইন্টারফেস এবং বায়বীয় অঞ্চলে বৃদ্ধি পায়

২. বেগ বিতরণ ফাংশন বৈশিষ্ট্য

বিভিন্ন অবস্থানে ( $x = 11, 13, 14.5, 15, 17$ ) বিশ্লেষণ প্রদর্শন করে:

তরল এবং ইন্টারফেস অঞ্চল: বেগ বিতরণ ম্যাক্সওয়েল বিতরণের কাছাকাছি
বায়বীয় অঞ্চল: ঋণাত্মক বেগ পরিসীমায় স্পষ্ট বিচ্যুতি দেখা যায়, আণবিক সংঘর্ষের অপর্যাপ্ততার কারণে
সীমানা-নিকটবর্তী অঞ্চল: আয়তন গতি প্রভাব উল্লেখযোগ্য, বিতরণে স্পষ্ট স্থানান্তর ঘটে

পরীক্ষামূলক আবিষ্কার

১. হার্টজ-নুডসেন সম্পর্কের সীমাবদ্ধতা: বায়বীয় অঞ্চলের বেগ বিতরণের অ-ম্যাক্সওয়েলীয় বৈশিষ্ট্য নিশ্চিত করে যে HK সম্পর্ক শক্তিশালী অ-সাম্যাবস্থা পরিস্থিতিতে অপ্রযোজ্য

২. দুটি মূল কারণ:

তরল-বায়বীয় ইন্টারফেসের নিকটবর্তী বায়বীয় অঞ্চলের বেগ বিতরণ ম্যাক্সওয়েল বিতরণ থেকে বিচ্যুত
বায়বীয় পার্শ্বের আয়তন গতির প্রভাব ক্রমবর্ধমান উল্লেখযোগ্য

উপসংহার ও আলোচনা

প্রধান উপসংহার

১. লেনার্ড-জোন্স তরলের জন্য প্রযোজ্য উন্নত গতিবিদ্যা মডেল সফলভাবে উন্নয়ন ২. সাম্যাবস্থা বৈশিষ্ট্য পূর্বাভাসে মডেল পরীক্ষামূলক তথ্যের সাথে উচ্চ সামঞ্জস্যপূর্ণ ৩. অ-সাম্যাবস্থা বাষ্পীভবনে বেগ বিতরণের অ-ম্যাক্সওয়েলীয় বৈশিষ্ট্য উন্মোচন ৪. হার্টজ-নুডসেন সম্পর্কের সীমাবদ্ধতার জন্য তাত্ত্বিক ভিত্তি প্রদান

সীমাবদ্ধতা

১. প্রযোজ্য পরিসীমা: প্রধানত সম্পৃক্ত বৈশিষ্ট্যের জন্য অনুকূলিত, অ-সম্পৃক্ত পরিস্থিতিতে নির্ভুলতা সীমিত ২. তরল প্রকার: বর্তমানে শুধুমাত্র নিষ্ক্রিয় গ্যাসের জন্য প্রযোজ্য, বহু-পরমাণুক অণুতে সম্প্রসারণের জন্য অভ্যন্তরীণ শক্তি মোড বিবেচনা প্রয়োজন ৩. জ্যামিতিক জটিলতা: প্রধানত সমতল ইন্টারফেস পরিস্থিতি যাচাই করা হয়েছে

ভবিষ্যৎ দিকনির্দেশ

১. অন্যান্য নিষ্ক্রিয় গ্যাসে সম্প্রসারণ ২. জটিল বহু-পরমাণুক তরলের জন্য প্রযোজ্য মডেল উন্নয়ন ৩. জটিল জ্যামিতিক পরিস্থিতিতে প্রয়োগ অধ্যয়ন

গভীর মূল্যায়ন

সুবিধা

১. পদ্ধতিগত উদ্ভাবন: এনস্কগ তত্ত্ব এবং ভ্লাসভ আনুমানিকতা দক্ষতার সাথে সমন্বিত, অভিজ্ঞতামূলক সংশোধনের মাধ্যমে প্রকৃত তরলের সঠিক বর্ণনা অর্জন ২. যাচাইয়ের পর্যাপ্ততা: একাধিক দৃষ্টিকোণ থেকে মডেলের নির্ভুলতা যাচাই, সাম্যাবস্থা বৈশিষ্ট্য এবং অ-সাম্যাবস্থা বৈশিষ্ট্য সহ ৩. তাত্ত্বিক অবদান: অ-সাম্যাবস্থা বাষ্পীভবন বোঝার জন্য গুরুত্বপূর্ণ অন্তর্দৃষ্টি প্রদান, বিশেষত HK সম্পর্কের সীমাবদ্ধতা উন্মোচন ৪. ব্যবহারিক মূল্য: গণনা দক্ষতা বজায় রেখে পূর্বাভাস নির্ভুলতা উল্লেখযোগ্যভাবে বৃদ্ধি

অপূর্ণতা

১. প্যারামিটার নির্ভরতা: মডেলে একাধিক অভিজ্ঞতামূলক প্যারামিটার রয়েছে, যা এর সাধারণীকরণ ক্ষমতাকে প্রভাবিত করতে পারে ২. তাত্ত্বিক ভিত্তি: কিছু সংশোধনের কঠোর তাত্ত্বিক উদ্ভাবনের অভাব ৩. যাচাই পরিসীমা: প্রধানত আর্গনের উপর ভিত্তি করে যাচাই, অন্যান্য তরলের প্রযোজ্যতা আরও নিশ্চিতকরণ প্রয়োজন

প্রভাব

১. একাডেমিক অবদান: অ-সাম্যাবস্থা তরল গতিবিদ্যার জন্য নতুন মডেলিং চিন্তাধারা প্রদান ২. প্রয়োগ সম্ভাবনা: ন্যানো প্রযুক্তি, তাপ ব্যবস্থাপনা ইত্যাদি ক্ষেত্রে গুরুত্বপূর্ণ প্রয়োগ মূল্য রয়েছে ৩. পুনরুৎপাদনযোগ্যতা: পদ্ধতি বর্ণনা স্পষ্ট, ভাল পুনরুৎপাদনযোগ্যতা রয়েছে

প্রযোজ্য পরিস্থিতি

ন্যানো-স্কেল বাষ্পীভবন শীতলীকরণ সিস্টেম ডিজাইন
ছিদ্রযুক্ত মাধ্যমে দশা পরিবর্তন তাপ স্থানান্তর
মাইক্রো-ন্যানো ডিভাইসের তাপ ব্যবস্থাপনা
পাতন এবং পৃথকীকরণ প্রক্রিয়া অনুকূলকরণ

তথ্যসূত্র

গবেষণাপত্রটি ৪১টি গুরুত্বপূর্ণ তথ্যসূত্র উদ্ধৃত করেছে, যা গতিবিদ্যা তত্ত্ব, পরীক্ষামূলক তথ্য এবং সংখ্যাসূচক পদ্ধতি ইত্যাদি একাধিক দিক কভার করে, গবেষণার জন্য দৃঢ় তাত্ত্বিক ভিত্তি প্রদান করে।