🔬 condensed matter

Continuum granular flow model with restitution-derived viscoelastic damping

تقدم هذه الورقة إطاراً متصلاً موحداً يجمع بين اللزوجة المرنة واللزوجة اللدنة للتدفقات الحبيبية، والذي يؤسس رابطاً مباشراً بين معامل الارتداد واللزوجة المتصلة لنمذجة السلوك المعتمد على المعدل مع الحفاظ على قاعدة التدفق اللدن الكلاسيكية $\mu(I)$ ، والتي تم تنفيذها عبر طريقة نقطة المادة لمحاكاة العمليات العابرة المعقدة.

المؤلفون الأصليون: Bodhinanda Chandra, Sachith Dunatunga, Ken Kamrin

نُشر 2026-04-24

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Bodhinanda Chandra, Sachith Dunatunga, Ken Kamrin

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

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

هذه الورقة البحثية تدور حول بناء "كتاب قواعد" عالمي للمحاكاة الحاسوبية يمكنه التنبؤ بدقة بكيفية سلوك هذه المواد الحبيبية في جميع تلك الحالات المختلفة. لقد ابتكر المؤلفون، بوديناندا شاندرا، وساتشيث دوناتونجا، وكين كامرين، نموذجًا رياضيًا جديدًا يجسّر الفجوة بين العالم الصغير للحبيبات الفردية والعالم الكبير للانهيارات الجليدية والترابية.

إليك تفصيل هذا الإنجاز، مشروحًا ببعض التشبيهات من الحياة اليومية.

١. المشكلة: معضلة "الارتداد" مقابل "السيولة"

في الماضي، واجهت النماذج الحاسوبية صعوبة في التعامل مع شيئين محددين في وقت واحد:

"الارتداد" (الاسترداد - Restitution): عندما تصطدم حبتان ببعضهما البعض، فإنهما لا ترتدان بشكل مثالي مثل الكرات المطاطية الخارقة؛ بل تفقدان بعض الطاقة. يُسمى هذا "معامل الاسترداد". إذا فقدتا الكثير من الطاقة، فإنهما تلتصقان؛ وإذا فقدتا القليل، فإنهما ترتدان.
"السيولة" (اللدونة - Plasticity): عندما تضغط على كومة من الرمل بقوة كافية، تبدأ في التدفق مثل العسل. هذا التدفق يعتمد على سرعة ضغطك.

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

٢. الحل: "ممتص صدمات ذكي"

ابتكر المؤلفون نموذجًا جديدًا يعمل مثل ممتص الصدمات الذكي في السيارة.

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

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

هذا أمر ضخم لأن هذا يمنع الكمبيوتر من التخمين. فبدلاً من القول: "لنضف بعض الاحتكاك العشوائي لجعله يبدو صحيحًا"، يقول النموذج: "بما أن هذه الحبيبات ترتد بنسبة ٨٠٪، فإن الكومة بأكملها يجب أن تمتلك هذا القدر المحدد من الاحتكاك الداخلي".

٣. قاعدة "عدم التوتر": الخيط غير المرئي

للمواد الحبيبية ميزة غريبة: يمكن ضغطها، لكن لا يمكن شدها.

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

٤. نظام "إشارات المرور" للإجهاد

أحد أكبر التحديات في هذه المحاكاة هو التأكد من أن "الارتداد" (اللزوجة المرنة) لا يفسد "السيولة" (اللدونة).

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

٥. ماذا أثبتوا؟ (التجارب)

اختبروا نموذج "ممتص الصدمات الذكي" الخاص بهم عبر خمسة سيناريوهات مختلفة، مثل سائق يختبر سيارة جديدة على تضاريس مختلفة:

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

الصورة الكبيرة

تعد هذه الورقة البحثية خطوة كبيرة للأمام لأنها توحد فيزياء الاصطدامات (ارتطام الحبيبات ببعضها) مع فيزياء التدفق (انهيار الرمال).

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

إنه جسر بين العالم المجهري لارتداد حبة واحدة والعالم الماكروسكوبي للانهيارات الأرضية، وكل ذلك متصل بمعادلة رياضية بسيطة وأنيقة.

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

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

تُظهر التدفقات الحبيبية الجافة سلوكيات معقدة تعتمد على المعدل ناتجة عن آليتين مجهريتين متمايزتين:

القصور الذاتي المجهري (Micro-inertia): يسود في التدفقات القصية البلاستيكية الكثيفة، مما يؤدي إلى احتكاك يعتمد على المعدل موصوف بـ ريولوجيا $\mu(I)$ .
الارتداد (التشتت التصادمي): يحكم فقدان الطاقة أثناء التصادمات المرنة بين الجسيمات. وهذا يؤثر على الاستجابات المرنة الكلية، مثل تخميد موجات الإجهاد والاهتزازات في الحالات الشبيهة بالصلبة، وفقدان الطاقة في الحالات الشبيهة بالغاز.

التحدي: بينما تم دمج القصور الذاتي المجهري بشكل جيد في النماذج المستمرة (على سبيل المثال، عبر نموذج $\mu(I)$ )، لا يوجد إطار عمل منهجي لدمج الارتداد المرن في الميكانيكا المستمرة دون تغيير سلوك التدفق البلاستيكي بشكل اصطناعي. غالبًا ما تربط النهج الحالية التشتت التصادمي مباشرة بالتشوه البلاستيكي أو تعتمد على معاملات تخميد عشوائية (مثل تخميد رايلي) تفتقر إلى رابط فيزيائي مباشر بمعامل الارتداد المجهري ( $e$ ). الهدف هو اشتقاق إطار عمل لزوجي مرن-لدن (viscoelastic-viscoplastic) موحد يربط معامل الارتداد المجهري باللزوجة المستمرة الكلية، مما يضمن أن التخميد اللزجي يؤثر على انتشار الموجات والتصادمات دون تعديل قاعدة التدفق البلاستيكي.

2. المنهجية

أ. الصياغة النظرية

يقترح المؤلفون إطار عمل مستمر لزوجي مرن-لدن موحد بناءً على المبادئ التالية:

التفكيك الكينماتيكي: يتم تفكيك تدرج السرعة ( $L$ ) جمعيًا إلى مكونات لزوجية مرنة ( $L_{ve}$ ) ولدن ( $L_p$ ).
تفكيك الإجهاد: باتباع ترتيب كيلفن-فويغت (Kelvin-Voigt)، يتم تقسيم الإجهاد الكلي إلى أجزاء مرنة ( $\sigma_e$ $σ_{e}$ ) ولزجة ( $\sigma_v$ $σ_{v}$ )، والتي يتم مساواتها بعد ذلك بالإجهاد البلاستيكي ( $\sigma_p$ $σ_{p}$ ) في الترتيب المتوالي.
- تمييز جوهري: يتم تطبيق التخميد اللزجي فقط على المكون المرن ( $\sigma_v$ تعتمد على $D_{ve}$ ). يضمن هذا أن التشتت اللزجي يخمد الموجات والتصادمات ولكن لا يضيف تخميدًا اصطناعيًا إلى التدفق البلاستيكي، الذي يظل محكومًا بدقة بريولوجيا $\mu(I)$ .
شرط عدم وجود شد (No-Tension Condition): يتضمن النموذج قاعدة فصل حيث تصبح المادة خالية من الإجهاد ( $\sigma=0$ ) إذا انخفضت الكثافة الحجمية عن قيمة حرجة ( $\rho_c$ )، مما يحاكي انفصال الحبيبات (الحالة الشبيهة بالغاز).

ب. اشتقاق الرابط بين الارتداد واللزوجة

المساهمة المركزية هي اشتقاق علاقة صريحة بين معامل الارتداد ( $e$ ) واللزوجة الحجمية المستمرة ( $\vartheta$ ).

المنهج: يحلل المؤلفون استجابة تجمع حبيبي كروي لنبضة ضغط حجمي. يعرّفون مقياس ارتداد مستمر ( $E$ ) كنسبة لمعدلات الانفعال الحجمي قبل وبعد التصادم.
تعريف زمن التلامس: اتباعًا لـ شواجر وبوشل، يتم تعريف زمن التلامس ( $t_c$ ) كالوقت الذي تتلاشى فيه قوة التلامس (الإجهاد)، وليس عندما تعود الإزاحة إلى الصفر. هذا يأخذ في الاعتبار الطبيعة أحادية الجانب لتلامس الحبيبات (لا يوجد شد).
النتيجة: من خلال حل معادلة الموجة المخمدة للتجمع الكروي، يشتقون علاقة مغلقة الشكل بين اللزوجة الحجمية غير البعدية ( $\tilde{\vartheta}$ $\tilde{ϑ}$ ) والارتداد المستمر ( $E$ $E$ ).
- علاقة تقريبية: $\tilde{\vartheta} \approx 0.237 \frac{|\ln E|}{\pi^2}$ .
- يسمح هذا بتحويل معامل الارتداد التجريبي ( $e$ ) مباشرة إلى اللزوجة الحجمية المستمرة ( $\vartheta$ ).
لزوجة القص: ترتبط لزوجة القص ( $\eta$ ) باللزوجة الحجمية عبر تخميد من نوع رايلي يتناسب مع الصلابة: $\eta = (G/K)\vartheta$ .

ج. التنفيذ العددي

المحلل (Solver): يتم تنفيذ النموذج ضمن طريقة نقطة المادة (MPM)، وهي إطار عمل هجين (أويلري-لاغرانجي) مثالي للتشوهات الكبيرة، وانفصال المواد، وإعادة التكتل.
خوارزمية تحديث الإجهاد: يتم تطوير مخطط تنبؤ-تصحيح (predictor-corrector).
1. حساب الإجهاد المرن التجريبي بما في ذلك الحدود اللزجة.
2. التحقق من شرط الخضوع ( $\tau \leq \mu(I)p$ ).
3. في حالة حدوث الخضوع، يتم إجراء خوارزمية "العودة إلى الخريطة" (return-mapping) لتحديث معدل الانفعال البلاستيكي والإجهاد، مما يضمن أن التدفق البلاستيكي مدفوع حصريًا بريولوجيا $\mu(I)$ بينما تعمل الحدود اللزجة على تخميد الاستجابة المرنة.

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

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

4. النتال والتحقق من الصحة

تم التحقق من صحة النموذج من خلال خمسة أمثلة عددية:

الضغط الكروي: تم التحقق من العلاقة التحليلية المشتقة بين اللزوجة الحجمية والارتداد المستمر. طابقت النتائج العددية التنبؤ النظري للاضمحلال الأسي في نسب السرعة.
تدفق المستوى المائل: أظهر أن تغيير معامل الارتداد (وبالتالي اللزوجة) لا يغير ملف سرعة باجنولد (Bagnold) في الحالة المستقرة. هذا يؤكد أن النموذج يحدد بشكل صحيح عملية فصل التخميد اللزجي عن قاعدة التدفق البلاستيكي، بما يتوافق مع نظرية $\mu(I)$ .
تدفق الصومعة ذو القاع المسطح: أظهر أنه بينما معدل التفريغ (المحكوم بالتدفق البلاستيكي) غير حساس لـ $e$ ، فإن سلوك إعادة التكتل حساس للغاية. أدى الارتداد المنخفض (التخميد العالي) إلى تقليل الارتداد غير الفيزيائي وأدى إلى مسافات جري أقصر وزوايا استقرار ديناميكية أكثر حدة، مما يطابق ملاحظات DEM.
اصطدام السرير الحبيبي: حقق المحاكي اصطدام جسم بكتلة حبيبية نجاحًا في تخميد موجات الإجهاد عالية التردد والترددات النابضة المرتبطة بالاهتزازات المرنة غير المخمدة، مع الحفاظ على التشوه البلاستيكي الصحيح وأنماط المقذوفات.
الأنماط الحبيبية المهتزة: أعاد النموذج إنتاج تكوين الأنماط المربعة والمعينية في الطبقات الحبيبية المهتزة رأسيًا. ومن المهم أن هذه الأنماط ظهرت فقط عند تضمين التخميد اللزج المرن (ارتداد محدود $e < 1$ )، مما يطابق الملاحظات التجريبية. الفشل في تكوين أنماط مستقرة في المحاكاة ذات المرونة التامة ( $e=1$ ).

5. الأهمية

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

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

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