Energy non-equipartition in vibrofluidized particles

ملخص تقني: عدم تساوي توزيع الطاقة في الجسيمات المتميعة بالاهتزاز

بيان المشكلة
غالبًا ما تنحرف نظرية تساوي توزيع الطاقة، التي تفترض أن الطاقة الحركية تتوزع بالتساوي بين جميع درجات الحرية في المائع، في الأنظمة الحبيبية. وبينما وثقت الدراسات السابقة عدم تساوي توزيع الطاقة في الأسِرّة المتميعة بالاهتزاز التي تتكون من جسيمات ملساء (مما يؤدي إلى تباين في الطاقة الحركية المتذبذبة) وجسيمات خشنة، فإن التأثير المحدد لمعلمات تفاعل الجسيمات الواقعية يظل موضوعًا للبحث. وتحديدًا، يتناول هذا العمل كيف تؤثر نسبة معامل صلابة النابض المماسي إلى النابض العمودي ($\kappa = k_t/k_n$) على تقسيم الطاقة بين الأنماط الانتقالية والدورانية. تسعى الدراسة إلى تحديد ما إذا كانت النماذج الواقعية لتلامس الجسيمات، لا سيما تلك التي تتضمن احتكاكًا عاليًا ونسب صلابة محددة، تؤدي إلى تساوي توزيع الطاقة أم إلى استمرار عدم تساوي التوزيع.

المنهجية
استخدم المؤلفون طريقة العناصر المنفصلة (DEM) باستخدام برنامج LAMMPS مفتوح المصدر لمحاكاة نظام حبيبي ثلاثي الأبعاد مهتز رأسيًا.

• نموذج التفاعل: تم نمذجة قوى التلامس باستخدام نموذج النابض المخمد الخطي (LSD). تم اختيار ثابت صلابة النابض العمودي ($k_n = 10^8 mg/d$) ليكون كبيرًا بما يكفي لضمان أن التصادمات هي تصادمات ثنائية في الغالب ($t_c/t_f \ll 1$).
• المعلمات: تباينت الدراسة في معامل الاحتكاك ($\mu$) من الجسيمات شبه الملساء ($\mu \approx 0$) إلى الجسيمات الخشنة تمامًا ($\mu = 10$). وتم استقصاء قيمتين رئيسيتين لنسبة الصلابة $\kappa$ وهما: $2/7$ و $3/4$. كما أُجريت عمليات محاكاة إضافية لـ $0.67 \le \kappa < 1$.
• بروتوكول المحاكاة: جرت عمليات المحاكاة على مرحلتين: مرحلة استرخاء أولية ($7 \times 10^7 t_c$) للوصول إلى حالة الاستقرار، تلتها مرحلة جمع البيانات بخطوة زمنية أصغر ($\Delta t = t_c/100$).
• التحليل: تم تحليل السرعات الخطية والزاوية اللحظية لتحديد ملامح متوسط الطاقة الحركية الانتقالية ($T$) والدورانية ($T_R$) المتذبذبة. وتم حساب المتوسطات الجماعية لأكثر من $10^3$ من التكوينات. واستُخدمت النسبة $K = K_{ET}/K_{ER}$ (الطاقة الحركية الانتقالية إلى الدورانية) كمعيار أساسي لتساوي توزيع الطاقة.

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

1. الاعتماد على نسبة الصلابة ($\kappa$): يعتمد سلوك نسبة الطاقة $K$ بشكل حاسم على $\kappa$.

   • عند $\kappa = 2/7$: تنخفض النسبة $K$ بشكل رتيب مع زيادة معامل الاحتكاك $\mu$. وعند الاحتكاك العالي ($\mu = 10$)، تقترب $K$ من الواحد الصحيح، مما يشير إلى تحقيق تساوي توزيع الطاقة بين الأنماط الانتقالية والدورانية.
   • عند $\kappa = 3/4$: السلوك غير رتيب؛ فبعد انخفاض أولي مع زيادة $\mu$، تزداد $K$ وتستقر عند قيمة تقارب 5. يشير هذا إلى استمرار عدم تساوي توزيع الطاقة، حتى عند معاملات الاحتكاك العالية جدًا.
   • أكدت عمليات المحاكاة لـ $0.67 \le \kappa < 1$ أن عدم تساوي توزيع الطاقة يستمر للجسيمات ذات معاملات الاحتكاك العالية في هذا النطاق.

2. ديناميكيات التصادم وتبادل الطاقة:

   • كشف تحليل التصادمات الفردية أنه بالنسبة لـ $\kappa = 2/7$، تزداد نسبة التلامسات في نظام "الالتصاق" (حيث يكون معامل الارتداد الدوراني $\beta \approx 1$) مع زيادة $\mu$. وعند $\mu = 1$، تكون حوالي 75% من التلامسات في هذا النظام الحافظ للطاقة، مما يسهل تساوي توزيع الطاقة.
   • بالنسبة لـ $\kappa = 3/4$، تستقر نسبة التلامسات في نظام "الالتصاق-الانزلاق" عند حوالي 55% لـ $\mu > 1$، حيث يكون $\beta < 1$. وبناءً على ذلك، تظل نسبة تبادل الطاقة $\Theta = E_{exch}/\Delta E$ أقل من 1 وتنخفض للجسيمات شديدة الخشونة. إذ يتجاوز التبديد عملية تبادل الطاقة بين الأنماط، مما يمنع تساوي توزيع الطاقة.

3. التحقق من النظرية: أظهرت نتائج المحاكاة لدرجة الحرارة الحبيبية عند التوازن ($T_o$) توافقًا ممتازًا مع التنبؤات النظرية (المبنية على توسيعات العزم من الرتبة الأولى) لكل من الجسيمات الملساء تمامًا والجسيمات الخشنة في الأنظمة المتخلخلة، بشرط استخدام معلمات القياس المناسبة (التي تتضمن $\beta$).

الأهمية والادعاءات
يزعم البحث أن اختيار نسبة الصلابة $\kappa$ هو عامل حاسم في التنبؤ بسلوك التدفق العياني للجسيمات الحبيبية الواقعية. وتوضح الدراسة أن تساوي توزيع الطاقة ليس نتيجة مضمونة في الأنظمة الحبيبية الخشنة المتميعة بالاهتزاز؛ بل هو مرهون بالتفاعل المحدد بين الاحتكاك ونسبة الصلابة المماسية إلى العمودية.

وتحديدًا، خلص المؤلفون إلى ما يلي:

• لوحظ تساوي توزيع الطاقة عند $\kappa = 2/7$ في حالات الاحتكاك العالي بسبب هيمنة تصادمات الالتصاق الحافظة للطاقة.
• لوحظ عدم تساوي توزيع الطاقة عند $0.67 \le \kappa < 1$ (بما في ذلك $\kappa = 3/4$) لأن تصادمات "الالتصاق-الانزلاق" في هذه الأنظمة ليست حافظة للطاقة، مما يؤدي إلى تفاوت كبير بين الطاقات الحركية الانتقالية والدورانية.

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