Transport properties in a model of confined granular mixtures at moderate densities

إليك ملخص تقني مفصل للورقة البحثية بعنوان "خصائص النقل في نموذج لخليط حبيبي محصور عند كثافات متوسطة" من تأليف ديفيد غونزاليس مينديز وفيسنتي غارزو.

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

تتناول الورقة البحثية تحدي وصف السلوك الهيدروديناميكي للخلائط الحبيبية المحصورة، ثنائية الأبعاد شبه المستوية عند كثافات متوسطة.

• السياق: غالبًا ما تكون الأنظمة الحبيبية محصورة في صناديق مهتزة حيث تقوم الاهتزازات الرأسية بضخ الطاقة، والتي تُنقل بعد ذلك إلى درجات الحرية الأفقية عبر التصادمات. هذا الإعداد يخلق حالة مستقرة غير متوازنة.
• الفجوة: كانت دراسات النظرية الحركية السابقة لهذه الأنظمة تقتصر إلى حد كبير على:
  • الأنظمة أحادية المكونات.
  • النظم الخفيفة (المتخلخلة) (معادلة بولتزمان).
  • حدود المتتبع (Tracer limits) (حيث يكون تركيز أحد المكونات ضئيلاً للغاية).
• الهدف: اشتقاق معادلات نافيير-ستوكس الهيدروديناميكية ومعاملات النقل لخليط مكون من $s$ من المكونات بـ تركيزات عشوائية عند كثافات متوسطة (حيث تكون الارتباطات بين الجسيمات مهمة)، باستخدام نموذج $\Delta$ لمراعاة تأثيرات الحصر.

2. المنهجية

يستخدم المؤلفون النظرية الحركية للغازات الحبيبية، وتحديداً توسيع نظرية إنسج (Revised Enslag Theory - RET) لتشمل التصادمات غير المرنة ونموذج $\Delta$.

• نموذج $\Delta$: هو نموذج خشن (coarse-grained) حيث يتم محاكاة حقن الطاقة الرأسية من الجدار المهتز عن طريق إضافة حد سرعة إضافي ($\Delta_{ij}$) إلى الحركة النسبية للجسيمات المتصادمة. يسمح هذا للنموذج بالوصول إلى حالة مستقرة رغم تبدد الطاقة.
• المعادلات الحاكمة:
  • تُستخدم معادلة إنسج غير المرنة لوصف تطور دوال توزيع السرعة أحادية الجسيم $f_i(\mathbf{r}, \mathbf{v}; t)$.
  • توسيع تشابمان-إنسج (Chapman-Enskog Expansion): يتم توسيع دوال التوزيع في قوى التدرجات المكانية للمجالات الهيدروديناميكية (الكثافة، سرعة التدفق، درجة الحرارة) حتى الدرجة الأولى.
  • الدرجة صفر: يُفترض أن تكون دوال التوزيع هي توزيع ماكسويل المحلي (المُقاس بدرجات الحرارة الجزئية) لتقدير معدلات التبريد والضغط.
  • الدرجة الأولى: يتم التعبير عن دوال التوزيع من خلال دوال مجهولة ($A_i, B_{ij}, C_{i,\lambda\beta}, D_i$) والتي تحقق مجموعة من المعادلات التكاملية الخطية المقترنة.
• التقريبات:
  • توسيع كثير حدود سونين (Sonine Polynomial Expansion): يتم تقريب الدوال المجهولة في المعادلات التكاملية بأول حدودها (أول كثير حدود سونين) للحصول على تعبيرات تحليلية لمعاملات النقل.
    يدرس التحليل حالة درجة حرارة مستقرة ($\zeta^{(0)} = 0$)، مما يبسط المعادلات التكاملية.
  • اللزوجة الحجمية: يتم إهمال مساهمة تصحيحات درجة الحرارة الجزئية من الدرجة الأولى ($\varpi_i$) في اللزوجة الحجمية، بناءً على نتائج سابقة تشير إلى أن هذا التأثير ضئيل.

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

1. التعميم على الخلائط متعددة المكونات: اشتقاق المجموعة الكاملة لمعاملات نقل نافيير-ستوكس لخليط مكون من $s$ من المكونات، متجاوزين بذلك حدود المتتبع والقيود الثنائية للأعمال السابقة.
2. نظام الكثافة المتوسطة: الحساب الصريح لكل من المساهمات الحركية والتصادمية في معاملات النقل، مع مراعاة دالة الارتباط الزوجي $\chi_{ij}$ (تأثيرات استبعاد الحجم).
3. التعبيرات التحليلية: اشتقاق تعبيرات تحليلية مغلقة (وإن كانت معقدة) لكل من:
   • معاملات الانتشار: الانتشار المتبادل ($D_{ij}$) والانتشار الحراري ($D_i^T$).
   • اللزوجة: لزوجة القص ($\eta$) واللزوجة الحجمية ($\eta_b$).
4. معيار الفصل (Segregation Criterion): تطوير معيار نظري يعتمد على عامل الانتشار الحراري ($\Lambda$) للتنبؤ بالانتقال بين تأثير البندق البرازيلي (BNE) وتأثير البندق البرازيلي العكسي (RBNE) في الخلائط الكثيفة المحصورة تحت تأثير الجاذبية والتدرجات الحرارية.

4. النتائج الرئيسية

أ. معاملات النقل

• الانتشار: تعتمد معاملات الانتشار على معاملات الاستعادة ($\alpha$)، ونسب الكتل، ونسب الأقطار، والتركيزات، والكثافة (الكسر الحجمي الصلب $\phi$).
  • تزيد عدم المرونة عموماً من معاملات الانتشار مقارنة بالحالة المرنة، وإن كان هذا التأثير أقل وضوحاً في نموذج $\Delta$ مقارنة بنموذج الكرات الصلبة غير المرنة القياسي (IHS).
  • تأثيرات الكثافة كبيرة: تزداد معاملات الانتشار عموماً مع زيادة الكثافة عندما تكون نسب الكتلة/القطر $>1$، وتنخفض عندما تكون النسب $<1$.
• اللزوجة:
  • لزوجة القص ($\eta$): تتكون من أجزاء حركية وتصادمية. توضح الورقة أنه بالنسبة للأنظمة المحصورة، تنخفض لزوجة القص الموزونة مع زيادة عدم المرونة. الاعتماد على الكثافة أضعف في نموذج $\Delta$ مقارنة بنموذج IHS القياسي.
  • اللزوجة الحجمية ($\eta_b$): يهيمن عليها الانتقال التصادمي. تم تحليل نسبة اللزوجة الحجمية إلى قيمتها المرنة، مما أظهر اعتماداً ضعيفاً على عدم المرونة مقارنة بنموذج IHS.
• التحقق: تمت مقارنة النتائج النظرية للزوجة القص في حد المكون الواحد مع محاكاة الديناميكا الجزيئية (MD) (من المرجع 14). تعيد النظرية إنتاج الاتجاهات نوعياً وتظهر توافقاً كمياً جيداً (الانحرافات $<2\%$) للكثافات المتوسطة، حتى عند عدم مرونة عالية ($\alpha \approx 0.5$).

ب. فصل الانتشار الحراري

تطبق الورقة المعاملات المشتقة لتحليل الفصل الناتج عن التدرج الحراري والجاذبية.

• عامل الانتشار الحراري ($\Lambda$): يُعرف بحيث يقابل $\Lambda > 0$ تأثير البندق البرازيلي (BNE) (الجسيمات الكبيرة ترتفع إلى الأعلى البارد) ويقابل $\Lambda < 0$ تأثير البندق البرازيلي العكسي (RBNE) (الجسيمات الكبيرة تغوص إلى الأسفل الساخن).
• المخططات الطورية: يرسم المؤلفون خط الانتقال ($\Lambda = 0$) في فضاء المعاملات لنسبة الحجم ($\sigma_1/\sigma_2$) ونسبة الكتلة ($m_1/m_2$).
  • غياب الجاذبية: تؤدي زيادة الكثافة قليلاً إلى توسيع منطقة RBNE (تجمع الجسيمات الكبيرة بالقرب من اللوح الساخن).
  • هيمنة الجاذبية: تكون منطقة RBNE أكبر بشكل عام، ولكن زيادة الكثافة يمكن أن توسع منطقة BNE اعتماداً على معايير الخليط المحددة.
  • عدم المرونة: يؤثر معامل الاستعادة بشكل كبير على حدود الفصل.

5. الأهمية والآثار المترتبة

• التقدم النظري: يجسّر هذا العمل الفجوة بين نظرية الغاز الحبيبي الخفيف ونظرية المادة الحبيبية الكثيفة، مما يوفر إطاراً حركياً صارماً للخلائط المحصورة التي يمكن تطبيقها في التجارب (مثل الصناديق المهتزة).
• القدرة التنبؤية: تسمح معايير الفصل المشتقة للباحثين بالتنبؤ بما إذا كان خليط معين من الجسيمات سينفصل أو يختلط تحت ظروف معينة من الاهتزاز، وتدرج الحرارة، والجاذبية.
• الصلة التجريبية: النموذج قابل للتطبيق مباشرة على تجارب الحبيبات ثنائية الأبعاد شبه المستوية. يشير الاتفاق مع محاكاة الديناميكا الجزيئية (MD) إلى أن النظرية موثوقة للكثافات المتوسطة ($\phi \lesssim 0.25$) وعدم المرونة المتوسطة.
• الاتجاهات المستقبلية: يشير المؤلفون إلى أنه لم يتم حل معاملات تدفق الحرارة بالكامل في هذه الورقة، ولكن الإطار قائم للعمل المستقبلي. كما يسلطون الضوء على إمكانية إجراء تحليل الاستقرار الخطي للحالة المستقرة المتجانسة والمقارنة مع نظريات السوائل ميكرو-بولارية (micropolar fluid theories).

باختصار، تقدم هذه الورقة وصفاً شاملاً للنظرية الحركية للنقل والفصل في الخلائط الحبيبية الكثيفة والمحصورة، وتوفر أدوات تحليلية لفهم الظواهر غير المتوازنة المعقدة في التدفقات الحبيبية الصناعية والطبيعية.