🔬 physics

Network modelling of yield-stress fluid flow in randomly disordered porous media

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

المؤلفون الأصليون: Cláudio P. Fonte, Elliott Sutton, Kohei Ohie, Eleanor Doman, Yuji Tasaka, Anne Juel

نُشر 2026-03-11

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

CC BY 4.0

المؤلفون الأصليون: Cláudio P. Fonte, Elliott Sutton, Kohei Ohie, Eleanor Doman, Yuji Tasaka, Anne Juel

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

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

إليك قصة اكتشافهم، مقسمة إلى مفاهيم بسيطة وتشبيهات من الحياة اليومية.

1. المشكلة: معجون الأسنان "العالق"

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

عندما تحاول هذه السوائل الانضغاط عبر متاهة فوضوية وعشوائية من الصخور (الأوسط المسامية)، تصبح الأمور معقدة:

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

2. الطريقة القديمة مقابل الطريقة الجديدة

حاول العلماء محاكاة هذا بطريقتين:

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

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

3. الحل الجديد: خريطة قائمة على الفيزياء

بنى مؤلفو هذه الورقة البحثية خريطة أكثر ذكاءً:

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

4. الاكتشاف الكبير: الأمر يتعلق بـ "العنق"، وليس "الجسم"

اختبر الباحثون خريطتهم الجديدة مقابل طريقة "الكمبيوتر الخارق" ووجدوا أنها دقيقة بشكل مذهل، رغم بساطتها.

لكن المفاجأة الكبرى جاءت عندما حاولوا فهم ما الذي يتحكم في الضغط المطلوب لدفع السائل عبر الوسط.

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

التشبيه: تخيل حشداً من الناس يحاول مغادرة ملعب رياضي.

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

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

5. لماذا هذا مهم؟

هذا النموذج الجديد يعد تغييراً جذرياً في قواعد اللعبة لأنه:

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

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

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

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

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

لا يبدأ التدفق إلا فوق عتبة ضغط محددة.
l. العلاقة بين معدل التدفق وهبوط الضغط هي علاقة غير خطية.
بالقرب من نقطة الخضوع (Yield point)، يصبح النقل غير متجانس للغاية، ويتميز بظاهرة التشكيل القنوي (Channelization) حيث يتم تحريك مجموعة فرعية فقط من المسام.
كما أن الانزلاق عند الجدار (Wall slip) (الشائع في المواد الحبيبية والمواد الزجاجية اللينة) يزيد من تعقيد الفيزياء من خلال تقليل الاحتكاك وتغيير عتبة الخضوع.

القيود الحالية:

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

2. المنهجية

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

بناء الهندسة:
- يتم نمذجة الوسط المسامي كدومين ثنائي الأبعاد ( $L \times L$ ) يحتوي على عوائق دائرية عشوائية غير متداخلة بنصف قطر $R$ .
- يحدد تجزئة فورونوي (Voronoi tessellation) لمراكز العوائق الشبكة: حيث تمثل الرؤوس المسام، وتمثل الحواف الممرات (Throats).
- تتم نمذجة هندسة الممر كـ فجوة متقاربة-متباعدة (Converging-diverging gap) بين العوائق المتجاورة، وتُعرف بنصف عرض أدنى قدر له $h_{min}$ .
فيزياء مقياس الممر (الإغلاق):
- بدلاً من استخدام معاملات مُعدة، يشتق النموذج علاقة الضغط بالتدفق ( $P(Q)$ ) مباشرة من ميكانيكا السوائل لممر متقارب-متباعد.
- المعادلات الحاكمة:
  - الريولوجيا: قانون هيرشل-باكلي ( $\eta(\dot{\gamma}) = \tau_0/\dot{\gamma} + K\dot{\gamma}^{n-1}$ ).
  - شرط الانزلاق: يتم فرض سرعة انزلاق $u_s$ عند الجدران، تعتمد على تجاوز إجهاد خضوع الانزلاق $\tau_s$ لإجهاد القص عند الجدار.
  - توازن الزخم: يتم حله لتدفق أحادي الاتجاه مستقر لربط تدرج الضغط المحلي $G(x)$ بمعدل التدفق $Q$ .
- التنفيذ العددي: يتم عكس العلاقة المحلية عددياً لإيجاد $G(x)$ لقيمة $Q$ معطاة، ويتم الحصول على إجمالي هبوط الضغط عبر الممر عن طريق تكامل $G(x)$ على طول الممر. يُستخدم معامل تسوية (Regularization parameter) لضمان المتانة العددية بالقرب من التدفق الصفري.
حل الشبكة:
- يقوم النظام بحل حفظ الكتلة عند كل مسام وعلاقات الضغط-التدفق غير الخطية عبر جميع الممرات.
- تستخدم طريقة من نوع نيوتن (Newton-type method) مع تعميم منطقة الثقة (Trust-region globalization) (منفذة بلغة Julia) لحل النظام المقترن.
- الاستمرار (Continuation): يتم حساب الحلول عن طريق مسح ضغط المدخل، باستخدام الحلول السابقة كتقديرات أولية للتعامل مع الصلابة (Stiffness) بالقرب من حد الخضوع.
التحقق من الصحة:
- تمت مقارنة النتائج مع المحاكاة العددية المباشرة (DNS) من تشاباريان وتامييسولا (2021) ومع محاكاة الحد النيوتوني الخاصة بالمؤلفين.
- تم اختبارها عبر هندسات متعددة (الحالات CT, G1, G2, G3) مع تغير المسامية ( $\phi$ ).

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

نموذج شبكة مسام خالٍ من المعاملات (Parameter-Free): هو أول نموذج شبكة مسام لتدفق اللزوجة البلاستيكية يدمج انزلاق الجدران وريولوجيا إجهاد الخضوع دون الحاجة لمعايرة تجريبية. يتم اشتقاق مقاومة الممر بالكامل من المبادئ الأولى (الهندسة والريولوجيا).
إدراج انزلاق الجدران: ينمذج صراحةً تأثير انزلاق الجدران على كل من هبوط الضغط الكلي وطوبولوجيا التدفق، وهي ميزة كانت غائبة سابقاً في نماذج الشبكة التنبؤية لهذا النظام.
تحديد مقياس الطول الحاكم: تتحدى الدراسة استخدام نصف قطر العائق ( $R$ ) أو مقاييس المسامية كمقياس طول أساسي للتدفق بالقرب من الخضوع. وقد حددت متوسط عرض الممر الأدنى ( $\bar{h}_{min}$ ) كالمقياس الهندسي الصحيح الذي يحكم التبديد في نظام الهيمنة البلاستيكية.

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

المقارنة والتحقق من الدقة:
- يعيد النموذج إنتاج الانحراف الماكروسكوبي عن سلوك دارسي بدقة، ملتقطاً الزيادة الرتيبة لتدرج الضغط اللابعدي ( $\mathcal{G}$ ) مع رقم بينجس ( $\mathcal{B}$ ).
- يتنبأ بشكل صحيح بالانتقال من التدفق الموزع مكانياً (رقم بينجس منخفض) إلى التدفق القنوي القوي (رقم بينجس مرتفع).
- التباين: يقلل النموذج من تقدير المقاومة قليلاً في نظام التدفق القنوي الشديد (بنسبة تصل إلى ~31% عند رقم بينجس مرتفع). ويُعزى ذلك إلى افتراض التدفق أحادي البعد في التضيقات الواسعة وإهمال التبديد داخل أجسام المسام (العقد). ومع ذلك، يظل أداة فعالة للغاية لتقليل الرتبة (Reduced-order tool).
تأثير انزلاق الجدران:
- يقلل الانزلاق بشكل كبير من تدرج الضغط المطلوب لاستمرار التدفق.
- تغيير الطوبولوجيا: في حالة عدم وجود انزلاق، يقتصر التدفق على مسارات قليلة جداً. ومع وجود الانزلاق، يتم تنشيط مسارات ثانوية، وتصبح شبكة التدفق أكثر اتصالاً.
- التحجيم (Scaling): يتغير التحجيم التقاربي لتدرج الضغط من $\mathcal{G} \sim \mathcal{B}$ (بدون انزلاق) إلى $\mathcal{G} \sim (\tau_s/\tau_0)\mathcal{B}$ (مع الانزلاق)، مما يشير إلى أن إجهاد خضوع الانزلاق يصبح عتبة مهيمنة.
قوانين التحجيم (الانهيار - The Collapse):
- يفشل التحجيم التقليدي باستخدام نصف قطر العائق $R$ في جعل البيانات تنهار (Collapse) عبر المسامية المختلفة.
- التحجيم باستخدام مقاييس المسامية ( $\phi/(1-\phi)$ ) يحسن الانهيار ولكنه يترك فروقات متبقية عند المسامية العالية.
- التحجيم الجديد: يؤدي إعادة تحجيم رقم بينجس وتدرج الضغط باستخدام متوسط عرض الممر الأدنى ( $\bar{{h}_{min}}$ ) (المعرف بـ $\mathcal{B}^*$ و $\mathcal{G}^*$ ) إلى انهيار عالمي (Universal collapse) للبيانات.
- في نظام الهيمنة البلاستيكية، تصبح العلاقة منحنى رئيسي: $\mathcal{G}^* = \mathcal{B}^*$ . وهذا يؤكد أن خسائر الضغط بالقرب من الخضوع محكومة بإحصاءات التضيقات بدلاً من مقياس العائق.

5. الأهمية

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