Rheology of dense vibrated granular flows: non-monotonic response controlled by granular temperature

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

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

تُظهر المواد الحبيبية سلوكيات ريولوجية معقدة، حيث تنتقل بين حالات تشبه الحالة الصلبة والحالات الشبيهة بالسوائل اعتمادًا على القوة الخارجية المؤثرة. وبينما تم توصيف ريولوجيا التدفقات الحبيبية غير المهتزة بشكل جيد (على سبيل المثال، عبر ريولوجيا $\mu(I)$)، فإن سلوك المواد الحبيبية الكثيفة تحت الاهتزاز الرأسي لا يزال أقل فهمًا، لا سيما فيما يتعلق بالتداخل بين معاملات الاهتزاز (السعة $A$، والتردد $f$) وضغط الحجز ($p$).

لقد حددت الدراسات السابقة اتجاهات متضاربة أو غير مكتملة:

• تشير بعض الأعمال إلى أن إضعاف الاحتكاك يتم التحكم فيه بواسطة النسبة $A^2/p$، مما يعني وجود عتبة لتعطيل سلاسل القوة.
• لاحظت أعمال أخرى اعتمادًا غير رتيب للزوجة والاحتكاك على تردد الاهتزاز $f$، حيث تزداد القدرة على التدفق لتصل إلى حد أقصى ثم تنخفض عند الترددات الأعلى.
• هناك نقص في وجود إطار عمل موحد يفسر كيف يحدد حقن الطاقة، والتبديد، والحجز، اللزوجة الفعالة ($\eta$) عبر النطاق الكامل للترددات، خاصة في النطاق الذي تزداد فيه اللزوجة مرة أخرى عند الترددات العالية.

٢. المنهجية

استخدم المؤلفون محاكاة طريقة العناصر المنفصلة (DEM) لمحاكاة إعداد تجريبي واقعي يُعرف باسم مقياس ريومتر الريشة المستحث بالإجهاد.

• إعداد المحاكاة:
  • النظام: نظام حبيبي أحادي التشتت ($N=2600$ جسيم، نصف القطر $r=2$ ملم) محصور في أسطوانة ذات قاعدة مخروطية وغطاء دائري كتلته $M$.
  • القوة المؤثرة: تخضع الحاوية لإزاحة جيبية رأسية $\Delta z(t) = A \sin(2\pi f t)$.
  • قياس القص: يتم تعريض ريشة مستطيلة مغمورة في الحبيبات لعزم دوران ثابت $T$. ويتم قياس السرعة الزاوية $\\Omega$ الناتجة لحساب اللزوجة الفعالة $\eta = T / |\langle \Omega \rangle|$.
  • نموذج الفيزياء: تم استخدام نموذج "هيرتز-ميندلين" (Hertz-Mindlin) لوصف تفاعلات الحبيبات مع بعضها، ومع الريشة، ومع الجدران، مع مراعاة القوى المرنة والمبددة.
• المعاملات المتغيرة:
  • سعة الاهتزاز ($A$): من $6$إلى$58$ ميكرومتر.
  • تردد الاهتزاز ($f$): تم تغييره لتغيير التسارع المعاد قياسه $\Gamma = (2\pi f)^2 A / g$.
  • ضغط الحجز ($p$): تم تغييره عن طريق تغيير كتلة الغطاء.
• المقاييس الرئيسية:
  • اللزوجة الفعالة ($\eta$).
  • درجة الحرارة الحبيبية ($K$): تُعرف بأنها متوسط الطاقة الحركية للحبيبات ($K \propto \langle v^2 \rangle$).
  • نسبة الطاقة ($\bar{K}$): نسبة طاقة اضطراب الحبيبات إلى مقياس الطاقة الناتج عن الحجز.

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

تقدم الورقة إطارًا فيزيائيًا موحدًا يوفق بين الملاحظات السابقة المتعلقة بإضعاف وتقوية اللزوجة. وتتمثل المساهمات الرئيسية فيما يلي:

1. تحديد المتغير الحاكم: أثبت المؤلفون أن اللزوجة الفعالة لا يتم التحكم فيها مباشرة بواسطة $A$ أو $f$ أو $p$ بشكل فردي، بل بواسطة النسبة بين درجة الحرارة الحبيبية ($K$) ومقياس طاقة الحجز ($K_p$).
   • تم تعريفها كـ $\bar{K} = K / K_p$، حيث $K_p \propto p$.
2. تفسير عدم الرتابة: توضح الورقة أن السلوك غير الرتيب للزوجة بالنسبة للتردد مدفوع بـ كفاءة نقل الطاقة.
   • عند الترددات المنخفضة، يؤدي زيادة $f$ إلى زيادة حقن الطاقة، مما يرفع $K$ ويخفض $\eta$ (التمييع).
   • عند الترددات العالية، يزداد معدل التبديد (معدل التصادم) بشكل أسرع من امتصاص الطاقة، مما يؤدي إلى انخفاض $K$، وهو ما يؤدي بدوره إلى زيادة $\eta$ (التغلظ).
3. التوفيق بين قوانين القياس: تتحقق الدراسة من قانون القياس $A^2/p$ الملاحظ في دراسات الاحتكاك السابقة، لكنها تقيد صلاحيته بـ نطاق الحد الأدنى للزوجة (حالة التمييع). وتوضح أن هذا القياس يفشل في التنبؤ بالسلوك عند الترددات العالية حيث يهيمن الاتجاه غير الرتيب.
4. النموذج الأدنى: تم تطوير نموذج ذي كتلتين (الملحق ج) لإعادة إنتاج قياس $A^2/p$، مما يظهر أنه ينشأ من التفاعلات المرنة غير الخطية (تلامسات هيرتز) وليس فقط من تمزق سلاسل القوة.

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

• اللزوجة مقابل التردد ($\eta$ مقابل $f$):

  • تُظهر اللزوجة الفعالة اتجاهًا غير رتيب: فهي تنخفض مع زيادة التردد حتى تصل إلى حد أدنى ($\eta_{min}$)، ثم تزداد عند الترددات الأعلى.
  • يرتبط هذا الاتجاه تمامًا بـ درجة الحرارة الحبيبية ($K$)، والتي تصل أيضًا إلى ذروتها عند تردد متوسط $f^*$.
  • توصف العلاقة بقانون القوة: $\eta \propto K^\alpha$ حيث $\alpha \in [-2, -1.5]$.

• اللزوجة مقابل ضغط الحجز ($\eta$ مقابل $p$):

  • تؤدي زيادة ضغط الحجز $p$ إلى زيادة اللزوجة الفعالة (مما يجعل النظام أكثر استقرارًا).
  • ومع ذلك، ليس لـ $p$ تأثير يذكر على درجة الحرارة الحبيبية $K$ نفسها؛ إذ يكمن تأثير $p$ في كونه حاجز الطاقة المطلوب لتحريك الحبيبات.

• المنحنى الرئيسي ($\eta$ مقابل $\bar{K}$):

  • عند رسم $\eta$ مقابل نسبة الطاقة عديمة الأبعاد $\bar{K} = K / (p \pi r^3)$، تندمج جميع نقاط البيانات لمختلف السعات والترددات والضغوط في منحنى رئيسي واحد.
  • يؤكد هذا أن الحالة الريولوجية تتحدد من خلال التنافس بين طاقة الاضطراب (التي تحاول كسر أقفاص الحبيبات) وطاقة الحجز (التي تحافظ على سلامة الأقفاص).

• محدودية $A^2/p$:

  • بينما يتناسب $\eta_{min}$ مع $A^2/p$، إلا أن هذا المتغير لا يمكنه التنبؤ باللزوجة عند الترددات العالية (الفرع الصاعد من المنحنى).
  • يعد قياس $A^2/p$ نتيجة لوصول النظام إلى حالة تمييع محددة، لكنه لا يستطيع التقاط آلية تبديد الطاقة المعتمدة على التردد.

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

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

في الختام، نجح المؤلفون في إثبات أن ريولوجيا التدفقات الحبيبية الكثيفة المهتزة محكومة بالتوازن بين الطاقة الحركية لاضطراب الحبيبات وطاقة الحجز المثبتة، مما يوفر نموذجًا تنبؤيًا قويًا للزوجة الفعالة عبر ظروف الدفع المتنوعة.