Second-gradient models for incompressible viscous fluids and associated cylindrical flows

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

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

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

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

يركز المؤلفان على إطار فريد-غورتين (Fried-Gurtin) للسوائل اللزجة الخطية غير القابلة للانضغاط ذات التدرج الثاني. وبينما يقدم هذا الإطار موتر إجهاد فائق (hyperstress tensor) من الدرجة الثالثة ($G$) لمحاكاة المشتقات المكانية العليا للسرعة، فإنه يعاني من فجوة نظرية حرجة: وهي أن الضغط الفائق ($\pi$)، وهو حقل متجهي يظهر في علاقة التكوين لـ $G$، يظل غير محدد. وبدون علاقة تكوين لـ $\pi$، يفتقر النموذج إلى القدرة التنبؤية، لا سيما في الشروط الحدية المتعلقة بالشد الفائق (hypertractions). علاوة على ذلك، فإن المحاولات السابقة لنمذجة اللزوجة المعتمدة على الضغط ضمن الميكانيكا الكلاسيكية غالباً ما تؤدي إلى مسائل سيئة التحديد حيث تفقد المعادلات الحاكمة خاصية الإهليلجية (ellipticity).

2. المنهجية

يطور المؤلفان إطاراً رياضياً صارماً لحل مشكلة عدم تحديد الضغط الفائق وتوسيع النموذج ليشمل اللزوجات المعتمدة على الضغط.

• الإغلاق التكويني للضغط الفائق: بناءً على الحجج المستمدة من الأدبيات السابقة (تحديداً [37])، يقترح المؤلفان علاقة تكوين محددة تربط الضغط الفائق بتدرج متوسط الضغط ($p$):
  $$\pi = \ell_1^2 \nabla p$$
  هنا، $\ell_1$ هو مقياس طول مميز مستمد من المعلمات الجوهرية للنموذج. هذا الاختيار يقلل من عدد الثوابت الحرة ويسمح بإغلاق معادلة الضغط.

• المعادلات الحاكمة: يتضمن النموذج توازن الزخم الخطي مع موازنة قدرة عطالية غير قياسية (بإدخال مقياس طول ثانٍ $\ell_0$) وعلاقة الإجهاد الفائق المقترحة. وتكون معادلة الزخم الناتجة هي:
  $$\rho \left( \dot{v} - \ell_0^2 \text{div}(\dot{\nabla}v) \right) = -\nabla(p - \ell_1^2 \Delta p) + \mu \Delta v - \mu \ell_1^2 \Delta^2 v + \rho b$$
  وبتطبيق التباعد (divergence) على هذه المعادلة، نحصل على معادلة إهليلجية من الدرجة الرابعة للضغط، على عكس معادلة بويسون من الدرجة الثانية في نافييه-ستوكس الكلاسيكية.

• توسيع اللزوجة المعتمدة على الضغط: يوسع المؤلفان النموذج ليشمل السوائل التي تكون فيها اللزوجة $\mu$ دالة في الضغط، $\hat{\mu}(p)$. ويثبتان أن إدراج حد التدرج الثاني ($\ell_1^2 \Delta^2 p$) يضمن بقاء معادلة الضغط إهليلجية بغض النظر عن سلوك $\hat{\mu}(p)$، مما يحل مشكلات التحديد الجيد (well-posedness) الموجودة في نماذج اللزوجة المعتمدة على الضغط الكلاسيكية.

• الشروط الحدية: تحلل الدراسة نوعين من شروط الالتصاق:

  • الالتصاق القوي: $v = 0$ و $\partial v / \partial n = 0$ عند الحدود.
  • الالتصاق الضعيف: $v = 0$ و الشد الفائق $m = 0$ عند الحدود.
    لتحديد الضغط بشكل فريد في النظام من الدرجة الرابعة، تم اقتراح شرط حدي إضافي وهو: $\partial p / \partial n = 0$ عند الحدود.

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

1. حل عدم تحديد الضغط الفائق: توفر الورقة علاقة تكوين مبررة فيزيائياً ($\pi = \ell_1^2 \nabla p$) تغلق نموذج فريد-غورتين، مما يجعله جيد التحديد رياضياً وقادراً على التنبؤ.
2. ضمان الإهليلجية للزوجة المعتمدة على الضغط: يثبت المؤلفان أنه بالنسبة للسوائل ذات اللزوجة المعتمدة على الضغط، فإن صياغة التدرج الثاني تضمن إهليلجية معادلة الضغط. وهذا يتناقض بشكل صارخ مع النماذج الكلاسيكية حيث يمكن لتدرجات الضغط العالية أن تغير نوع المعادلة (تصبح زائديّة/hyperbolic)، مما يؤدي إلى عدم استقرار رياضي.
3. الحلول التحليلية الصريحة: يستنتج المؤلفان حلولاً دقيقة للتدفقات الأسطوانية المستقرة (Poiseuille و Taylor-Couette) تحت كل من شروط الالتصاق القوي والضعيف. تتضمن هذه الحلول دالات بيسل المعدلة ($I_0, K_0, I_1, K_1$).
4. التقارب التقاربي: تم تقديم براهين صارمة تظهر أنه عندما تؤول مقاييس الطول المميزة ($\ell_1, \ell_0$) إلى الصفر، فإن ملفات السرعة والضغط المستمدة تتقارب نقطياً نحو حلول نافييه-ستوكس الكلاسيكية.

4. النتائج

تطبق الورقة النموذج على مسألتين كلاسيكيتين للتدفق الأسطواني:

• تدفق بوازيه (التدفق المستقر في أنبوب):

  • الالتصاق القوي: يتضمن ملف السرعة حداً مكافئاً كلاسيكياً بالإضافة إلى حد تصحيحي يتضمن $I_0(r/\ell_1)$. كما يتم تعديل معدل التصريف بواسطة مقياس الطول.
  • الالتصاق الضعيف: يحقق الحل الشرط $G[n \otimes n] = 0$. يختلف ملف السرعة عن حالة الالتصاق القوي ولكنه لا يزال يتقارب نحو الملف المكافئ مع $\ell_1 \to 0$.
  • التقارب: باستخدام التوسعات التقاربية لدالات بيسل، يثبت المؤلفون أن السرعة غير البعدية $u(\sigma)$ تتقارب إلى $1-\sigma^2$ (الحل الكلاسيكي) عندما $\lambda_1 \to 0$.

• تدفق تايلور-كويت (أسطوانة دوارة):

  • الالتصاق القوي: يتضمن ملف السرعة $I_1(r/\ell_1)$ ويتقارب نحو الملف الخطي $v = \Omega r$ مع $\ell_1 \to 0$.
  • الالتصاق الضعيف: ومن المثير للاهتمام، بالنسبة للالتصاق الضعيف، يحقق الحل الخطي الكلاسيكي $v = \Omega r$ الشروط الحدية تماماً، مما يعني أن تأثيرات التدرج الثاني تتلاشى لهذا التكوين المحدد للتدفق تحت هذه الشروط.
  • حقل الضغط: يتم تحديد الضغط بواسطة معادلة تفاضلية عادية من الدرجة الرابعة. وتؤكد الحلول العددية (باستخدام bvp4c في MATLAB) أن حقل الضغط غير البعدي يتقارب نحو الحل الكلاسيكي مع تلاشي مقاييس الطول.

5. الأهمية

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

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