A semi-implicit double-point Material Point Method for both free-surface flow and seepage in deformable porous media

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

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

قصة الماء "الذكي" والرمل "الزلق"

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

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

لحل ذلك، طبخ الفريق وصفة جديدة بخمسة مكونات سرية:

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

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

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

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

5. مزيج "المثبت" (TPIC 1% و Delta-Correction)
عندما يتحرك الماء بسرعة، يمكن للنقاط الصغيرة (نقاط المادة) التي يستخدمها الحاسوب لتتبعه أن تتداخل، فتتكتل معاً أو تتطاير بعيداً، مما يؤدي إلى انفجار المحاكاة. استخدم المؤلفون مزيجاً من تقنيتين: إحداهما تحافظ على الطاقة (FLIP) والأخرى تنعم الحركة (TPIC). وجدوا أن استخدام 99% من طريقة الحفاظ على الطاقة و1% فقط من طريقة التنعيم هو التوازن المثالي. الأمر يشبه قيادة السيارة: تريد محركاً قوياً (طاقة)، ولكنك تحتاج إلى القليل جداً من امتصاص الصدمات (التنعيم) للحفاظ على الرحلة من الاهتزاز. كما أضافوا "تصحيح دلتا"، وهو بمثابة دفعة لطيفة تنشر نقاط الماء المتكتلة حتى لا تتراكم وتكسر القواعد الرياضية.

الدليل في "النتيجة النهائية" (وفي التسونامي)

لم يكتفِ المؤلفون بكتابة هذه الأفكار، بل اختبروها مقابل تجارب واقعية. لقد أجروا عمليات محاكاة لـ:

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

الخلاصة

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

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