Self-similar rupture of thin films of power-law fluid

ملخص تقني: التمزق ذاتي التشابه للأغشية الرقيقة من سوائل قانون القوة

صياغة المسألة
تتقصى هذه الدراسة التمزق في زمن محدد للأغشية السائلة الرقيقة المكونة من سوائل "قانون القوة" (power-law fluids)، والمدفوع بالمنافسة بين التوتر السطحي وقوى التجاذب بين الجزيئات (قوى فان دير فالس). وبينما يعد التمزق ذاتي التشابه للأغشية النيوتونية أمراً راسخاً، فإن سلوك السوائل غير النيوتونية، التي تتميز بأس قانون القوة $n$ (حيث يشير $n < 1$ إلى ترقق القص و $n > 1$ إلى تثخن القص)، يطرح تحديات تحليلية ورقمية كبيرة. المعادلة الحاكمة هي معادلة تفاضلية جزئية (PDE) من الدرجة الرابعة غير الخطية مشتقة من نظرية التزييت (lubrication theory)، وتتضمن حداً لضغط التجزئة لنمذجة تجاذب فان دير فالس. الهدف الأساسي هو حساب الحلول ذاتية التشابه التي تصف ملف الغشاء مع اقترابه من زمن التمزق $t_0$، ورسم خريطة لوجود واستقرار هذه الحلول عبر فضاء المعاملات $n$.

المنهجية
يستخدم المؤلفون مزيجاً من المحاكاة الرقمية المباشرة، والاستمرار العددي، والتحليل التقاربي:

1. المحاكاة الرقمية للمعادلة التفاضلية الجزئية: يتم حل معادلة الغشاء الرقيق المعتمدة على الزمن باستخدام طريقة الفروق المحدودة مع مخطط زمني ضمني (ode15s في MATLAB). ولتجنب المشكلات المتعلقة بالمشتقات غير الصحيحة المتأصلة في نماذج قانون القوة، يتم حساب التدفق مباشرة من ملف السماكة باستخدام نموذج خماسي النقاط للمشتقات، بدلاً من حساب المشتق الرابع المكاني للسماكة مباشرة.
2. بناء الحل ذاتي التشابه: يتم اختزال المعادلة التفاضلية الجزئية إلى نظام من المعادلات التفاضلية العادية (ODEs) عبر القياس ذاتي التشابه. ويتم حل مسألة القيم الحدية (BVP) الناتجة باستخدام حزمة الاستمرار العددي Auto07p. ويتتبع المؤلفون فروع الحلول عن طريق تغيير أس قانون القوة $n$، بدءاً من الحلول النيوتونية المعروفة ($n=1$).
3. التعامل مع حدود $n$ الصغيرة: تفشل عملية الاستمرار العددي القياسية عند قيم $n$ الصغيرة ($n \lesssim 0.2$) بسبب ظهور منطقة داخلية صغيرة أسياً حيث يتلاشى التدفق. وللتغلب على ذلك، قدم المؤلفون إعادة قياس للإحداثيات بناءً على متغير التدفق، مما أدى إلى تحويل النطاق لاستيعاب التغيرات السريعة في الطبقة الداخلية.
4. التحليل التقاربي: بالنسبة للحد $n \to 0$، يقوم المؤلفون بإجراء توسع تقاربي متوافق. يتضمن هذا منطقة خارجية، وطبقة داخلية يتلاشى فيها التدفق أسياً، ومنطقة داخلية للغاية. يكشف هذا التحليل عن بنية الحل ويبرر النهج العددي المعاد قياسه.

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

• بنية التفرع: يُظهر مخطط التفرع للحلول ذاتية التشابه، المفهرسة بـ $n$، بنية "تلوّي" (snaking) غير بديهية حول $n=1$.
  • بالنسبة لـ $n > 1$، يمكن فقط للفرع الأولي (الذي يمتلك أكبر سماكة مركزية $H(0)$ عند $n=1$) أن يستمر إلى قيم $n$ كبيرة جداً. بينما تتلاشى جميع الفروع الأخرى عبر تفرعات الطي (fold bifurcations) مع زيادة $n$.
  • بالنسبة لـ $n < 1$، يمكن فقط للفروع الأربعة الأولى (التي تمتلك أكبر $H(0)$ عند $n=1$) الاستمرار نحو $n \to 0$. أما الفروع المتبقية فتتلاشى عبر تفرعات الطي مع تناقص $n$، مما يخلق نمط اندماج معقداً بالقرب من $n=1$.
• اختيار الجاذب: تتلاقى المحاكاة الرقمية للمعادلة التفاضلية الجزئية المعتمدة على الزمن باستمرار نحو الفرع الأساسي الوحيد للحلول ذاتية التشابه. هذا الفرع هو الوحيد الذي يستمر في كل من حدود تثخن القص القصوى ($n \to \infty$) وترقق القص القصوى ($n \to 0$).
• تقارب $n$ الصغير: في الحد $n \to 0$، تمتلك الحلول ذاتية التشابه منطقة داخلية صغيرة أسياً حيث يصبح التدفق ضئيلاً. حدد التحليل التقاربي عدداً لا نهائياً من الحلول ذاتية التشابه في هذا الحد. المسألة الخارجية من الدرجة الأولى تتطلب حلاً عددياً، وتكشف شروط المطابقة عن بنية محددة تتضمن منطقة داخلية للغاية محكومة بمعادلة تفاضلية عادية من الدرجة الثالثة غير خطية ولا تملك حلاً مغلق الصيغة.
• خصائص الملف: تُظهر فروع التشابه من الرتب العليا سلوكاً تذبذبياً في المتغير ذاتي التشابه $\eta$. وتكون هذه التذبذبات أكثر وضوحاً لقيم $n$ الأصغر، وهي مرتبطة بالاختيار المنفصل للحلول المطلوب لتلبية الشروط الحدية، وهي ظاهرة تتعلق بظاهرة ستوكس (Stokes' phenomenon) في الحالة النيوتونية.

الأهمية والادعاءات
يزعم البحث توسيع نطاق حساب الحلول ذاتية التشابه من السوائل النيوتونية إلى سوائل قانون القوة غير النيوتونية، كاشفاً عن بنية تفرع معقدة لم تُستكشف من قبل. وتتمثل المساهمة الرئيسية في تحديد سلوك التفرع "المتلوّي" بالقرب من $n=1$ وإثبات أن الفرع الأولي هو الوحيد الذي يعمل كجاذب فيزيائي لديناميكيات التمزق.

يسلط المؤلفون الضوء على تحديين تقنيين هامين تمت معالجتهما:

1. عدم وجود المشتقات: تمت إدارة عدم وجود المشتقات المكانية العليا في تدفقات التزييت لقانون القوة عند نقاط تدرج الضغط الصفري من خلال صياغة المسألة بدلالة التدفق بدلاً من مشتقات السماكة العليا.
2. المناطق الصغيرة أسياً: كشفت الدراسة عن صعوبات رقمية شديدة في حد ترقق القص الشديد ($n \to 0$) بسبب المناطق الداخلية الصغيرة أسياً. ومن خلال اشتقاق البنية التقاربية، قدم المؤلفون مخططاً عددياً معاد قياسه قادراً على حساب الحلول في هذا النظام، وهي قدرة يزعمون أنها ضرورية لمحاكاة سوائل الترقق الشديد في المشكلات السطحية العامة.

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