🔬 condensed matter

Rheology of dense suspensions of granular spherocylinders by particle-based simulation

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

المؤلفون الأصليون: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

نُشر 2026-07-31

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: Alex Dixon, John Hone, Gavin Melaugh, Christopher Ness

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

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

في هذه الورقة البحثية، قام فريق من الباحثين ببناء محاكاة حاسوبية فائقة الذكاء لتعمل كمختبر افتراضي لهذه المعلقات من العصي المزدحمة. وبدلاً من خلط مواد كيميائية حقيقية في المختبر، أنشأوا عالماً رقمياً مليئاً بآلاف "الأسطوانات الكروية" (spherocylinders) — وهي جسيمات تشبه الكبسولات أو الحبوب (أسطوانة مع نصف كرة عند كل طرف). جعلوا هذه الكبسولات تطفو في سائل افتراضي، ثم بدأوا في عملية "القص" (shearing)، وهي طريقة فنية تعني أنهم سحبوا الطبقة العليا من السائل جانبياً بينما ظلوا يمسكون الطبقة السفلى ثابتة، تماماً مثل فرد الزبدة على قطعة توست. كان الهدف هو معرفة كيف تتغير لزوجة الخليط (viscosity) مع زيادة كثافة العصي في المساحة ومع زيادة طول العصي.

كشفت المحاكاة عن سلوكيات مثيرة للاهتمام. أولاً، عندما بدأوا عملية القص، أصبح الخليط سميكاً للغاية للحظة وجيزة — "طفرة في اللزوجة" (viscosity spike) — قبل أن يستقر. حدث هذا لأن العصي كانت في البداية مبعثرة في اتجاهات عشوائية، وتصطدم ببعضها البعض مثل حشد في حفلة "موش بيت" (mosh pit). وبمجرد بدء التدفق، بدأت العصي في الاصطفاف، مثل الجنود الذين يسيرون في تشكيل عسكري، وأصبح الخليط أسهل في التحريك. وجد الباحثون أنه كلما زاد عدد العصي التي يحشرونها (الكسر الحجمي، $\phi$ )، أصبح الخليط أكثر سمكاً، وصولاً إلى نقطة حيث يتكدس ويتوقف عن التدفق تماماً. ومن المثير للاهتمام أن العصي كلما كانت أطول (نسبة الجوانب، $A$ ، والتي تتراوح من 1 للكرة حتى 20 للعصي الطويلة جداً)، حدث هذا التكدس في وقت أبكر. على سبيل المثال، مع عصي بطول 5، كان نقطة التكدس عند كسر حجمي قدره 0.57، ولكن بالنسبة لعصي بطول 20، تكدس الخليط عند كثافة أقل بكثير تبلغ حوالي 0.41.

اكتشف الفريق أيضاً أن العصي لم تكتفِ بالاصطفاف فحسب، بل اصطفت بشكل محدد مع اتجاه التدفق. لقد قاسوا "معامل الترتيب" ( $S$ ) لمعرفة مدى جودة اصطفافها. في الكثافات المنخفضة، كانت العصي تتأرجح وتدور في حلقات (ما يعرف بمدارات جيفري - Jeffery orbits)، ولكن مع زيادة ازدحام الحشد، انغلقوا في خط مستقيم على طول اتجاه التدفق. ومع ذلك، إذا أصبح الحشد كثيفاً للغاية، بدأت العصي تشعر بالإحباط ولم تعد قادرة على البقاء مصطفة بشكل مثالي، مما أدى إلى انخفاض الترتيب قليلاً. كما قام الباحثون بتفكيك القوى المؤثرة، موضحين أنه عندما كان الخليط رقيقاً، كانت مقاومة السائل نفسه هي الأكثر أهمية، ولكن مع زيادة الازدحام، أصبحت الاصطدامات الفيزيائية والاحتكاك بين العصي هي القوة المهيمنة، حيث تولت مهمة جعل الخليط سميكاً.

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

ملخص تقني: ريولوجيا المعلقات الكثيفة من الحبيبات متطاولة الشكل (Spherocylinders) عبر المحاكاة القائمة على الجسيمات

بيان المشكلة
تنتشر المعلقات الكثيفة من الجسيمات الحبيبية ذات الشكل القضيب في سياقات طبيعية وصناعية، بدءاً من سدات جذوع الأشجار في الأنهار وصولاً إلى المعلقات البلورية في الصهارة وشبكات الألياف في تصنيع الورق. وبينما تم توصيف ريولوجيا المعلقات الكروية بشكل جيد، لا سيما من خلال توحيد المعلقات غير البراونية والتدفقات الحبيبية الجافة بواسطة "بوير وآخرون" (Boyer et al.)، إلا أن سلوك الجسيمات المستطيلة لا يزال غير مفهوم بشكل كامل. إن درجة الحرية الإضافية المرتبطة باتجاه الجسيمات تقدم ظواهر بنيوية مجهرية معقدة، مثل مدارات جيفري (Jeffery orbits) والمحاذاة الناجمة عن القص، والتي تغير اللزوجة الفعالة واستجابة الإجهاد بشكل كبير. وتعد البيانات التجريبية المتاحة للأنظمة الكثيفة (الكسور الحجمية الصلبة العالية، $\phi$ ) محدودة بسبب صعوبات التحضير والقياس. علاوة على ذلك، تواجه مناهج المحاكاة السابقة مقايضات: حيث تعد طريقة "ستوكسيان دايناميكس" (Stokesian Dynamics) مكلفة حاسوبياً للأنظمة الكثيفة بسبب عمليات قلب المصفوفات المطلوبة للهيدروديناميكا متعددة الأجسام، بينما غالباً ما تبسط نماذج "طريقة العناصر المنفصلة" (DEM) هندسة الجسيمات (مثل سلاسل الكرات) أو تهمل القوى الهيدروديناميكية الضرورية. لذا، هناك حاجة لنموذج قائم على الجسيمات يمكن معالجته حاسوبياً بفعالية، ويجسد بدقة الفيزياء المجهرية المهيمنة للمعلقات الكثيفة من الحبيبات غير البراونية ومتطاولة الشكل تحت تأثير القص.

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

معادلات الحركة: يتم تحديث مسارات الجسيمات باستخدام خوارزمية "فيليت-فيرلي" (Velocity-Verlet)، التي تدمج معادلات الحركة الانتقالية والدورانية. تأخذ هذه المعادلات في الاعتبار قوى التلامس بين الزوج، وقوى التزييت (lubrication) بين الزوج، والقوى الهيدروديناميكية أحادية الجسم (السحب، الرفع، وعزم الدوران).
الهيدروديناميكا: يتم نمذجة قوى الجسم الواحد باستخدام حلول تحليلية للأشكال البيضاوية المتطاولة (prolate spheroids) كتقريب للأسطوانات متطاولة الشكل. يتضمن ذلك السحب المعتمد على الاتجاه والرفع الناجم عن القص. يعيد عزم الدوران الهيدروديناميكي إنتاج مدارات جيفري.
التلامس والتزييت:
- التلامس: يتم نمذجته كنموذج "هوك" (Hookean) مع احتكاك "كولوم" (Coulombic)، بما في ذلك مكونات النابض والمخمد (spring-damper) العمودية والمماسية.
- التزييت: تم تطوير مخطط مبتكر لحساب قوى التزييت بين الأسطوانات متطاولة الشكل. ولمعالجة الانقطاع في الانحناء عند نقطة التقاء الجزء الأسطواني بالنهاية، يقوم النموذج بحساب ثلاثة مكونات تفاعل (جذع-جذع، نهاية-جذع، ونهاية-نهاية) في آن واحد. يتم وزن هذه المكونات بناءً على الموقع النسبي للجسيمات لضمان انتقالات سلسة بين أنواع التفاعل، مما يحل مشكلة الانقطاعات الموجودة في النماذج السابقة (مثل ياماني وآخرون، وباتلر وشافيه).
الخطوة الزمنية الديناميكية: يعد حساب الخطوة الزمنية الديناميكية ابتكاراً منهجياً رئيسياً. يتم تحديد الخطوة الزمنية بناءً على أصغر مقياس زمني مميز بين تفاعلات التلامس، وقوى التزييت، وعزوم دوران التزييت. يضمن هذا الاستقرار العددي مع تجنب الخطوات الزمنية الصغيرة غير الضرورية، مما يحسن الكفاءة الحاسوبية.
إعداد المحاكاة: يستخدم النظام حدوداً دورية في اتجاه التدفق واتجاه التدرج، وحدود "ليز-إدواردز" (Lees-Edwards) للقص. يُستخدم خليط ثنائي التشتت من أنصاف الأقطار (1:1.4) لمنع التبلور. تغطي المحاسبات نسب الاستطالة ( $A$ ) من 1 إلى 20 والكسور الحجمية ( $\phi$ ) من المخفف إلى مرحلة ما قبل التكدس (near-jamming).

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

سلوك اللزوجة: يتنبأ النموذج بزيادة منهجية في اللزوجة المستقرة ( $\eta_s$ $η_{s}$ ) مع كل من الكسر الحجمي ( $\phi$ $ϕ$ ) ونسبة الاستطالة ( $A$ $A$ ).
- الاستجابة العابرة: عند بدء عملية القص، لوحظت طفرة في اللزوجة، تُعزى إلى التوزيع العشوائي الأولي لاتجاهات الجسيمات. تخبو هذه الطفرة مع محاذاة الجسيمات، مما يؤدي إلى حالة مستقرة. تزدلاد قيمة الطفرة بالنسبة للحالة المستقرة مع زيادة $\phi$ .
- الحالة المستقرة: تتبع اللزوجة في الحالة المستقرة تباعداً وفق قانون القوة $\eta_s \sim (1 - \phi/\phi_m)^{-\varsigma}$ . تنخفض نسبة التكدس ( $\phi_m$ ) مع زيادة $A$ (من $\approx 0.57$ لـ $A=5$ إلى قيم أقل لـ $A$ الأعلى). يؤكد النموذج استقلالية المعدل عن اللزوجة في الحالة المستقرة، وهو ما يتوافق مع إطار العمل للحبوب المعلقة.
- مساهمات القوة: عند $\phi$ منخفضة، تهيمن إجهادات المائع وقوى التزييت على اللزوجة. ومع زيادة $\phi$ ، تصبح قوى التلامس (العمودية والمماسية) هي المهيمنة، لتتجاوز مساهمات التزييت في النهاية حول $\phi \approx 0.5$ .
البنية المجهرية والمحاذاة:
- معامل الترتيب: يزداد معامل الترتيب الاتجاهي ( $S$ ) مع زيادة $\phi$ حتى نقطة حرجة، وبعدها ينخفض قليلاً بسبب الإحباط الهندسي الذي يمنع المحاذاة المثالية.
- اتجاه المحاذاة: تصطف الجسيمات تفضيلياً على طول اتجاه التدفق ( $x$ -axis)، مع زاوية ميل صغيرة ( $\phi_v$ ) بالنسبة لمستوى التدفق-التدرج. تزداد زاوية الميل هذه مع زيادة $\phi$ حتى حد معين ثم تنخفض.
- تأثير نسبة الاستطالة: تؤدي نسب الاستطالة العالية إلى محاذاة أقوى عند كسور حجمية منخفضة، ولكنها تؤدي أيضاً إلى نقطة تكدس أقل. بالنسبة لـ $A$ العالية جداً، تصبح زاوية الميل مهملة.
- العابر مقابل المستقر: ترتبط طفرة اللزوجة الأولية بالوقت اللازم لانتقال النظام من الحالة العشوائية إلى حالة المحاذاة الناتجة عن القص.

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

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

أبحاث مشابهة