Strongly non-linear continuous model for wave propagation in 2D granular media: derivation and hyperbolicity analysis

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

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

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

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

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

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

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

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