DOI: https://doi.org/10.1017/jfm.2024.605
تاريخ النشر: 2024-09-10
المؤلف: Shaokang Li وآخرون
الموضوع الرئيسي: ديناميات الغاز ونظرية الحركة
نظرة عامة
تقدم البحث نموذجًا جديدًا لنظرية الحركة الجزيئية مصمم لمحاكاة تدفقات السائل-البخار، مع التركيز بشكل خاص على عمليات التبخر على النانو. يعالج هذا النموذج قيود أدوات المحاكاة الحالية من خلال التقاط التغيرات الطورية غير المتوازنة بدقة مع الحفاظ على الكفاءة الحسابية. على عكس النماذج السابقة التي افترضت ظروفًا متساوية الحرارة، يأخذ النهج المقترح في الاعتبار تقلبات درجة الحرارة أثناء التبخر دون المطالب الحسابية العالية لطرق الحركة المتوسطة الأكثر تعقيدًا. تم التحقق من صحة النموذج مقابل حلول مرجعية في ثلاثة سيناريوهات: توازن السائل-البخار، التبخر في بيئة شبه فراغية، والتبخر في البخار، مما يظهر تخفيضات كبيرة في وقت المحاكاة – تصل إلى مرتبتين من حيث الحجم – مع الحفاظ على الدقة.
يستخدم النموذج الحركي إمكانات ساذرلاند لتقريب تفاعلات الجسيمات، حيث يجمع بين إمكانات الطرد الصلبة والإمكانات الجذابة اللينة. يتم نمذجة التفاعلات الطاردة من خلال تكامل تصادم تقريبي، بينما يتم التعامل مع القوى الجذابة باستخدام تقريب المجال المتوسط. تظهر التحقق من صحة النموذج مقابل الحلول التحليلية والعددية توافقًا ممتازًا، خاصة عند درجات الحرارة العالية، مع بعض الانحرافات عند درجات الحرارة المنخفضة لا تزال تلتقط السلوكيات الفيزيائية الأساسية. ستركز الأبحاث المستقبلية على توسيع النموذج إلى أبعاد أعلى وتطوير شروط حدودية لمحاكاة أفضل للتفاعلات مع الجدران الصلبة، مما يعزز قابليته للتطبيق في دراسة ظواهر ديناميكا السوائل المعقدة مثل تكتل القطرات وخصائص ترطيب السوائل. بشكل عام، يضع هذا العمل النموذج الحركي كأداة واعدة لتقدم تطبيقات تكنولوجيا النانو المتعلقة بالتبخر على النانو.
مقدمة
تسلط المقدمة الضوء على أهمية التبخر في سياقات طبيعية وتقنية متنوعة، خاصة على النانو، حيث أدت التقدمات في تكنولوجيا النانو إلى تحسين التطبيقات مثل التبريد بالتبخر وفصل المواد باستخدام أغشية نانوية مسامية. يتم التأكيد على الحاجة إلى أدوات حسابية دقيقة لمحاكاة تدفق التبخر على النانو، حيث تتضمن هذه العمليات تفاعلات معقدة عبر مقاييس مكانية وزمنية متعددة. يتميز مجال تدفق السائل المتبخر بأربعة مناطق متميزة: الكتلة السائلة، الواجهة، طبقة كنوذسن، وكتلة البخار، حيث تقدم كل منها تحديات فريدة للنمذجة بسبب الحالات غير المتوازنة والتقلبات الكبيرة في المتغيرات الكلية.
تناقش هذه الفقرة السياق التاريخي لنمذجة التبخر، مشيرة إلى صيغة هيرتز-كنوذسن (H-K) وقيودها، خاصة عدم قدرتها على حساب سرعة البخار في الاتجاه السفلي. على الرغم من أن التعديلات اللاحقة، مثل معادلة شراج، قد حسنت التنبؤات، إلا أنها لا تزال تعاني من نقص تحت ظروف معينة. تلاحظ المقدمة أيضًا أنه بينما تم تطوير طرق نظرية وعددية، بما في ذلك طرق اللحظات ونهج محاكاة مونت كارلو المباشرة (DSMC) لدراسة ديناميات التبخر، فإنها غالبًا ما تعتمد على شروط حدودية ظاهرة تقدم عدم يقين كبير بسبب تباين معامل التبخر المبلغ عنه في الأدبيات. وهذا يبرز ضرورة وجود نماذج أكثر دقة يمكنها التقاط الديناميات بدقة عند واجهة السائل-البخار.
النتائج
في قسم النتائج، يقدم البحث شروط حدودية متنوعة تم تطبيقها على نظام تم نمذجته تحت سيناريوهات مختلفة. على وجه الخصوص، تستخدم الحالة (ب) شروط حدودية لجدار ماص لمحاكاة بيئة شبه فراغية، مما يؤدي إلى الحالة \( f \mp L_2, \xi_x = 0 \) لـ \( \xi_x \gtrless 0 \). في الحالة (ج)، يتم استخدام شروط حدودية بعيدة لتمثيل منطقة البخار، المميزة بتوزيع ماكسويل المنجرف مع سرعة \( u_\infty \)، المعبر عنها كالتالي:
\[
f \mp L_2, \xi_x = \frac{n_\infty m}{\sqrt{2\pi k_B T_\infty}} \exp\left(-\frac{m(\xi_x – u_\infty)^2}{2k_B T_\infty}\right), \quad \xi_x \gtrless 0.
\]
تبسط التحليل من خلال التركيز على اتجاه x، مما يؤدي إلى إدخال دوال توزيع مختزلة تم الحصول عليها من خلال التكامل على مكونات السرعة \( \xi_y \) و \( \xi_z \).
يتم جعل النظام المترابط الناتج من معادلات التوزيع المختزلة بعدم الأبعاد باستخدام القطر الجزيئي \( \ell_0 \) وأعلى سرعة جزيئية محتملة عند درجة الحرارة المرجعية \( T_0 \). يتم تعريف الوقت المرجعي \( t_0 \)، وحقل القوة المتوسطة، ودالة التوزيع، وكثافة العدد وفقًا لذلك. يتم حل المعادلات باستخدام طريقة السرعة المنفصلة، التي تستخدم شبكات موحدة في المجالات المكانية والسرعة، مع اختلافات منتهية من الدرجة الثانية للمشتقات ومخطط أويلر الصريح لتكامل الزمن، مما يحقق دقة من الدرجة الثانية في الفضاء ودقة من الدرجة الأولى في الزمن. حددت دراسات استقلال الشبكة حجم شبكة مكانية قدره \( \ell_0/10 \) وحجم شبكة سرعة قدره \( u_0/20 \) لمحاكاة عبر النطاق \([-4u_0, 4u_0]\). يتم الإشارة إلى الكميات غير البعدية بعلامة التلدة لزيادة الوضوح.
مناقشة
في هذا القسم، يناقش المؤلفون تطوير والتحقق من صحة نموذج حركي يهدف إلى محاكاة تدفقات السائل-البخار، مع التركيز بشكل خاص على معامل التبخر. يستخدم النموذج معادلة إنسكوغ-فلاسوف (EV)، التي تشمل كل من الطور السائل والبخار، بالإضافة إلى المنطقة الواجهة. تمتد معادلة EV معادلة بولتزمان من خلال تضمين مصطلح تصادم إنسكوغ لأخذ تأثيرات السوائل الكثيفة في الاعتبار ومصطلح قوة ذاتية متسقة لفلاسوف للتفاعلات الجزيئية بعيدة المدى. على الرغم من قوتها النظرية، تواجه معادلة EV قيودًا في التطبيقات العملية بسبب كثافتها الحسابية وفرضية تفاعلات الكرات الصلبة. يقترح المؤلفون نموذجًا حركيًا مبسطًا يقرب مصطلح تصادم إنسكوغ باستخدام توسيع سلسلة تايلور ومصطلح استرخاء شاكوف، مما يعزز الكفاءة الحسابية مع الحفاظ على الدقة.
يتم تقييم أداء النموذج من خلال محاكاة ثلاثة سيناريوهات: توازن السائل-البخار، التبخر في ظروف شبه فراغية، والتبخر في البخار. تشير النتائج إلى أن النموذج الحركي يلتقط بدقة الفيزياء الأساسية لتفاعلات السائل-البخار، مع توافق جيد مع كل من معادلة EV والحلول التحليلية، خاصة عند درجات الحرارة العالية. ومع ذلك، تظهر بعض التباينات عند درجات الحرارة المنخفضة، تُعزى إلى التبسيطات في تكامل التصادم. يسمح النموذج أيضًا بتقدير معامل التبخر، مما يكشف أن اختيار درجة الحرارة عند واجهة السائل-البخار يؤثر بشكل كبير على النتائج. بشكل عام، يظهر هذا العمل نهجًا واعدًا لنمذجة تدفقات السائل-البخار مع تقليل المطالب الحسابية، مما يمهد الطريق للأبحاث المستقبلية والتطبيقات في نظرية الحركة.
DOI: https://doi.org/10.1017/jfm.2024.605
Publication Date: 2024-09-10
Author(s): Shaokang Li et al.
Primary Topic: Gas Dynamics and Kinetic Theory
Overview
The research presents a novel molecular kinetic theory model designed to simulate liquid-vapour flows, particularly focusing on nanoscale evaporation processes. This model addresses the limitations of existing simulation tools by accurately capturing non-equilibrium phase changes while maintaining computational efficiency. Unlike previous models that assumed isothermal conditions, the proposed approach accounts for temperature variations during evaporation without the high computational demands of more complex mean-field kinetic methods. The model has been validated against benchmark solutions in three scenarios: liquid-vapour equilibrium, evaporation into near-vacuum, and evaporation into vapour, demonstrating significant reductions in simulation time—up to two orders of magnitude—while maintaining accuracy.
The kinetic model employs a Sutherland potential to approximate particle interactions, combining hard-sphere repulsive and soft-attractive potentials. The repulsive interactions are modeled through an approximated collision integral, while attractive forces are treated with a mean-field approximation. Validation against analytical and numerical solutions shows excellent agreement, particularly at higher temperatures, with some deviations at lower temperatures still capturing essential physical behaviors. Future research will focus on extending the model to higher dimensions and developing boundary conditions to better simulate interactions with solid walls, thereby enhancing its applicability in studying complex fluid dynamics phenomena such as droplet coalescence and fluid wetting properties. Overall, this work positions the kinetic model as a promising tool for advancing nanotechnology applications related to nanoscale evaporation.
Introduction
The introduction highlights the significance of evaporation in various natural and technological contexts, particularly at the nanoscale, where advancements in nanotechnology have enhanced applications such as evaporative cooling and material separation using nanoporous membranes. The need for accurate computational tools to simulate nanoscale evaporative flow is emphasized, as these processes involve complex interactions across multiple spatial and temporal scales. The flow field of an evaporating liquid is characterized by four distinct regions: the liquid bulk, the interface, the Knudsen layer, and the vapour bulk, each presenting unique challenges for modeling due to non-equilibrium states and substantial variations in macroscopic variables.
The section discusses the historical context of evaporation modeling, referencing the Hertz-Knudsen (H-K) formula and its limitations, particularly its inability to account for downstream vapour velocity. Although subsequent modifications, such as Schrage’s equation, have improved predictions, they still fall short under certain conditions. The introduction also notes that while theoretical and numerical methods, including moment methods and the Direct Simulation Monte Carlo (DSMC) approach, have been developed to study evaporation dynamics, they often rely on phenomenological boundary conditions that introduce significant uncertainty due to the variability of the evaporation coefficient reported in the literature. This underscores the necessity for more refined models that can accurately capture the dynamics at the liquid-vapour interface.
Results
In the results section, the study presents various boundary conditions applied to a system modeled under different scenarios. Specifically, case (b) employs absorbing wall boundary conditions to simulate a near-vacuum environment, leading to the condition \( f \mp L_2, \xi_x = 0 \) for \( \xi_x \gtrless 0 \). In case (c), far-field boundary conditions are utilized to represent the vapor region, characterized by a drifting Maxwellian distribution with velocity \( u_\infty \), expressed as:
\[
f \mp L_2, \xi_x = \frac{n_\infty m}{\sqrt{2\pi k_B T_\infty}} \exp\left(-\frac{m(\xi_x – u_\infty)^2}{2k_B T_\infty}\right), \quad \xi_x \gtrless 0.
\]
The analysis simplifies by focusing on the x-direction, leading to the introduction of reduced distribution functions obtained by integrating over the velocity components \( \xi_y \) and \( \xi_z \).
The resulting coupled system of reduced distribution equations is made dimensionless using the molecular diameter \( \ell_0 \) and the most probable molecular velocity at the reference temperature \( T_0 \). The reference time \( t_0 \), mean force field, distribution function, and number density are defined accordingly. The equations are solved using the discrete velocity method, which employs uniform meshes in spatial and velocity domains, with second-order finite differences for derivatives and an explicit Euler scheme for time integration, achieving second-order accuracy in space and first-order accuracy in time. Mesh independence studies determined a spatial grid size of \( \ell_0/10 \) and a velocity grid size of \( u_0/20 \) for simulations across the range \([-4u_0, 4u_0]\). Dimensionless quantities are denoted with a tilde for clarity.
Discussion
In this section, the authors discuss the development and validation of a kinetic model aimed at simulating liquid-vapour flows, particularly focusing on the evaporation coefficient. The model utilizes the Enskog-Vlasov (EV) equation, which incorporates both liquid and vapour phases, along with the interfacial region. The EV equation extends the Boltzmann equation by including an Enskog collision term to account for dense fluid effects and a Vlasov self-consistent force term for long-range molecular interactions. Despite its theoretical robustness, the EV equation faces limitations in practical applications due to its computational intensity and the assumption of hard-sphere interactions. The authors propose a simplified kinetic model that approximates the Enskog collision term using a Taylor series expansion and a Shakhov-like relaxation term, thus enhancing computational efficiency while maintaining accuracy.
The model’s performance is evaluated through simulations of three scenarios: liquid-vapour equilibrium, evaporation into near-vacuum conditions, and evaporation into vapour. Results indicate that the kinetic model accurately captures the essential physics of liquid-vapour interactions, with good agreement with both the EV equation and analytical solutions, particularly at higher temperatures. However, some discrepancies arise at lower temperatures, attributed to the simplifications in the collision integral. The model also allows for the estimation of the evaporation coefficient, revealing that the choice of temperature at the liquid-vapour interface significantly influences the results. Overall, this work demonstrates a promising approach to modeling liquid-vapour flows with reduced computational demands, paving the way for future research and applications in kinetic theory.