Entropy Scaling for Diffusion Coefficients in Fluid Mixtures

1. بيان المشكلة

تعد معاملات الانتشار بالغة الأهمية لتصميم عمليات فصل السوائل، والمفاعلات، وأنظمة الاحتراق. ومع ذلك، فإن البيانات التجريبية للانتشار في المخاليط نادرة، لا سيما في الأنظمة غير المثالية والحالات القصوى (فوق الحرجة، وشبه المستقرة).

• الفجوة: بينما يعد تدرج الإنتروبيا (الذي يربط الخصائص الانتقالية بالإنتروبيا التشكيلية) أمراً راسخاً للمكونات النقية ولزوجة المخاليط والناقلية الحرارية، إلا أنه لم يُطبق بنجاً في معاملات انتشار المخاليط.
• التحدي: النماذج التجريبية الحالية (مثل نماذج فيجينز، وداركن) غالباً ما تفشل في التعامل مع المخاليط شديدة عدم المثالية، ولا يمكنها وصف العلاقة بين الانتشار الذاتي ($D_i$)، وانتشار ماكسويل-ستيفان ($\mathcal{D}_{ij}$)، والانتشار الفيكي ($D_{ij}$) بشكل متسق عبر كامل نطاق المائع (الغاز، السائل، فوق الحرج، وشبه المستقر). علاوة على ذلك، لا يوجد علاقة عامة تربط بين معاملات الانتشار الذاتي ومعاملات الانتشار المتبادل.

2. المنهجية

يقترح المؤلفون إطار عمل لتدرج الإنتروبيا متسق ديناميكياً حرارياً يتنبأ بجميع معاملات الانتشار في المخاليط الثنائية دون الحاجة إلى بارامترات (معاملات) ضبط للمخلوط. يعتمد هذا النهج على ثلاثة مفاهيم أساسية:

أ. معالجة المكونات شبه النقية

يعامل الإطار معاملات انتشار التخفيف اللانهائي ($D^\infty_i$) كـ "مكونات شبه نقية".

• تبيّن أن معامل الانتشار عند التخفيف اللانهائي ($D^\infty_i$) يظهر سلوك تدرج أحادي المتغير فيما يتعلق بالإنتروبيا التشكيلية المختزلة ($\tilde{s}_{conf}$)، بشكل مشابه للانتشار الذاتي للمكونات النقية.
• يسمح هذا لـ $D^\infty_i$ بأن يتم نمذجته باستخدام نفس ارتباطات تدرج الإنتروبيا المطورة للموائع النقية، مما يتطلب فقط بضع نقاط بيانات مرجعية (أو حتى لا شيء إذا استُخدمت المعاملات العالمية).

ب. صياغة تدرج الإنتروبيا

يستخدم النموذج نهج روزنفيلد المعدل مدمجاً مع نظرية تشابمان-إنسكوغ (CE) لحد الكثافة المنخفضة.

• متغير التدرج: الإنتروبيا التشكيلية المختزلة $\tilde{s}_{conf} = -s_{conf}/(R m)$، حيث يتم اشتقاق $s_{conf}$ من معادلة حالة (EOS) قائمة على الجزيئات.
• الحالات الحدية: يحدد النموذج دوال التدرج لأربع حالات حدية:
  1. الانتشار الذاتي للمكون النقي ($D^{pure}_i$).
  2. الانتشار عند التخفيف اللانهائي للمكون النقي ($D^\infty_i$).
  3. الانتشار الذاتي للمكون النقي للنوع الثاني ($D^{pure}_j$).
  4. الانتشار عند التخفيف اللانهائي للمكون النقي للنوع الثاني ($D^\infty_j$).
• الشكل الدالي: يتم وصف معامل الانتشار المختزل $\hat{D}$ كدالة مستمرة لـ $\tilde{s}_{conf}$ تتضمن معاملات خاصة بالنظام ($\alpha_2, \alpha_3$) يتم ضبطها بناءً على الحالات الحدية.

ج. قواعد الخلط والتركيب

للتنبؤ بالانتشار عند أي تركيب كيميائي ($x_i$)، يستخدم الإطار قواعد خلط خطية للتدخل بين الحالات الحدية:

• الانتشار الذاتي ($D_i, D_j$): يُحسب عن طريق خلط معاملات التدرج للمكون النقي والتخفيف اللانهائي بناءً على الكسور المولية.
• انتشار ماكسويل-ستفان ($\mathcal{D}_{ij}$): يُحسب عن طريق خلط معاملات التخفيف اللانهائي لكلا المكونين.
• الانتشار الفيكي ($D_{ij}$): يُشتق من معامل ماكسويل-ستفان والعامل الديناميكي الحراري ($\Gamma_{ij}$)، والذي يُحسب مباشرة من معادلة الحالة ($\Gamma_{ij} = \frac{x_i x_j}{RT} \frac{\partial^2 G}{\partial x_i^2}$).

الميزة الرئيسية: يربط الإطار نموذج تدرج الإنتروبيا مع معادلة حالة قائمة على الجزيئات (مثل PC-SAFT للموائع الحقيقية، أو Kolafa-Nezbeda لنموذج Lennard-Jones). يضمن هذا الربط الاتساق الديناميكي الحراري، مما يسمح للنموذج بالتعامل مع التوازن الطوري (VLE, LLE) والحالات شبه المستقرة بشكل طبيعي.

3. المساهمات الرئيسية

1. إطار عمل موحد: أول تطبيق ناجح لتدرج الإنتروبيا للتنبؤ بجميع معاملات الانتشار (الذاتي، وماكسويل-ستفان، والفيكي) في المخاليط بشكل متسق.
2. مفهوم المكون شبه النقي: إثبات أن معاملات الانتشار عند التخفيف اللانهائي تتبع قانون تدرج أحادي المتغير عالمي مع الإنتروبيا التشكيلية، من خلال معاملتها كمكونات شبه نقية.
3. التنبؤ بدون بارامترات ضبط: بمجرد توصيف الحالات الحدية (باستخدام الحد الأدنى من البيانات)، يتنبأ النموذج بمعاملات الانتشار عبر كامل نطاق التركيب وحيز الحالة دون أي بارامترات ضبط للمخلوط.
4. قابلية تطبيق واسعة: يعمل النموذج للمواد الغازية، والسوائل، والموائع فوق الحرجة، وحتى الحالات شبه المستقرة/غير المستقرة (مثل مناطق السبينودال)، وهي حالات لم تستطع النماذج السابقة التعامل معها.

4. النتائج

قام المؤلفون بالتحقق من صحة النموذج باستخدام محاكاة الديناميكا الجزيئية (MD) لمخاليط لينارد-جونز (LJ) وبيانات تجريبية لأنظمة المواد الحقيقية.

• مخاليط لينارد-جونز:

  • تنبأ النموذج بدقة بمعاملات الانتشار الذاتي وماكسويل-ستفان عبر كامل نطاق التركيب، بما في ذلك الأنظمة ذات عدم المثالية القوية، والنقاط القصوى، وفجوات التوازن سائل-سائل (LLE).
  • رصد بشكل صحيح السلوك عند النقطة الحرجة حيث يقترب معامل الانتشار الفيكي من الصفر.
  • تفوق على النماذج التجريبية (Vignes وDarken)، التي فشلت في رصد الاتجاهات في الأنظمة غير المثالية وLLE.

• أنظمة المواد الحقيقية:

  • تم اختبار النموذج على مخاليط مثل (n-hexane + n-dodecane، و acetone + chloroform، و toluene + n-hexane).
  • الدقة: كانت الانحرافات النسبية المتوسطة منخفضة بشكل عام (مثلاً ~1-5% للانتشار الذاتي في المخاليط غير القطبية).
  • عدم المثالية: نجح النموذج في رصد السلوكيات غير الخطية الناتجة عن الاختلافات الكبيرة في الكتلة المولية أو التفاعلات القوية.
  • القيود: لوحظت انحرافات طفيفة في الأنظمة المرتبطة (مثل 2-propanol + n-heptane) وبالقرب من القيم القصوى للانتشار الفيكي، ويرجع ذلك على الأرجح إلى الهيكلية المحلية المعقدة التي لا تستطيع معادلة الحالة الحالية استيعابها بالكامل.

• الحالات شبه المستقرة:

  • نجح الإطار في التنبؤ بمعاملات الانتشار في المناطق شبه المستقرة وغير المستقرة (مثل داخل منطقة السبينودال)، مما أظهر قدرات تتجاوز بكثير حدود النمذجة التجريبية أو التجريبية الحالية.

5. الأهمية

• القدرة التنبؤية: تتيح هذه الطريقة التنبؤ بمعاملات الانتشار في المناطق التي تفتقر للبيانات التجريبية (مثل الضغط العالي، الحالة فوق الحرجة، والحالة شبه المستقرة)، وهو أمر حيوي لتصميم العمليات في الظروف القاسية.
• الاتساق الديناميكي الحراري: من خلال الربط مع معادلة حالة (EOS)، يضمن النموذج أن تنبؤات الانتشار متسقة مع التوازن الطوري والاستقرار الديناميكي الحراري، وهي ميزة تفتقر إليها الارتباطات التجريبية.
• الكفاءة: يقلل من الحاجة إلى بيانات تجريبية مكثفة أو محاكاة جزيئية معقدة لكل مخلوط جديد؛ حيث يتطلب الأمر فقط بيانات الحالات الحدية (النقية والتخفيف اللانهائي) لتوليد ملف الانتشار الكامل.
• أساس للأعمال المستقبلية: توفر المنهجية أساساً فيزيائياً قوياً لتوسيع تدرج الإنتروبيا ليشمل الأنظمة متعددة المكونات والموائع المرتبطة الأكثر تعقيداً.