Navier-Stokes-Cahn-Hilliard system in a $3$D perforated domain with free slip and source term: Existence and homogenization

إليك ملخص تقني مفصل للورقة البحثية بعنوان "نظام نافييه-ستوكس-كان-هيلارد في مجال ثلاثي الأبعاد مثقوب مع انزلاق حر وحر عيني" من إعداد أمارتيا تشاكرابورتي، وهارادان دوتا، وهاري شانكار ماهاتو.

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

تتقصى الورقة البحثية التحليل الرياضي والسلوك التقاربي لنظام نافييه-ستوكس-كان-هيلارد (NSCH) الذي ينمذج خليط سوائل ثنائي غير قابل للانضغاط داخل وسط مسامي ثلاثي الأبعاد مثقوب دورياً.

السمات الفيزيائية والرياضية الرئيسية:

• المجال (Domain): مجال مثقوب $\Omega^\varepsilon_p \subset \mathbb{R}^3$ يمثل مسام السائل، محاط بكتل صلبة $\Omega^\varepsilon_s$. الهندسة دورية مع معامل مقياس صغير $\varepsilon > 0$.
• المعادلات الحاكمة:
  • الزخم (Momentum): معادلات نافييه-ستوكس المعتمدة على الزمن مع تنسور لزوجة متغير $A^\varepsilon(t,x)$ وحد كورتويغ (Korteweg-type) لقوة الشعيرية $\lambda_\varepsilon \phi^\varepsilon \nabla \mu^\varepsilon$.
  • حفظ الكتلة: قيد عدم الانضغاط $\nabla \cdot u^\varepsilon = 0$.
  • مجال الطور (Phase Field): معادلة كان-هيلارد مع حد مصدر غير محافظ $G(\phi^\varepsilon)$ (يمثل تفاعلات حجمية أو دفعاً خارجياً)، مقترنة بالجهد الكيميائي $\mu^\varepsilon = -\Delta \phi^\varepsilon + F'(\phi^\varepsilon)$.
• الشروط الحدية:
  • الحد الخارجي ($\partial \Omega$): شرط عدم الانزلاق ($u^\varepsilon = 0$).
  • الكتل الصلبة ($\Gamma^\varepsilon_s$): شروط الانزلاق الحر (انزلاق نافييه مع إجهاد مماس صفري): $u^\varepsilon \cdot n^\varepsilon = 0$ و $(A^\varepsilon D(u^\varepsilon)n^\varepsilon) \cdot \tau^\varepsilon = 0$.
• التحجيم (Scaling): تتبع قوة الشعيرية $\lambda_\varepsilon$ مقياساً بحيث $\lim_{\varepsilon \to 0} \lambda_\varepsilon = \lambda \in [0, +\infty)$. وتتحجم قوة الجسم الخارجية $g^\varepsilon$ كـ $\sqrt{\lambda_\varepsilon}$.

2. المنهجية

يستخدم المؤلفون نهجاً تحليلياً صارماً يتكون من خطوتين:

أ. الوجود والتقديرات الموحدة (على المستوى المجهري)

بالنسبة لـ $\varepsilon > 0$ ثابت، يثبت المؤلفون وجود حلول ضعيفة ويشتقون تقديرات موحدة في $\varepsilon$.

• تقريب غاليركين (Galerkin Approximation): يبنون الحلول التقريبية باستخدام مخطط غاليركين رباعي الخطوات يعتمد على الدوال الذاتية من نوع ستوكس المتكيفة مع الشروط الحدية المختلطة (عدم انزلاق على الحد الخارجي، وانزلاق حر على الكتل الصلبة).
• مؤثرات التمديد (Extension Operators): يمثل التحدي التقني الحرج هو المجال المعتمد على $\varepsilon$. قام المؤلفون ببناء:
  • مؤثرات بوغوفسكي (Bogovskiĭ Operators): معرفة على كل من مجال السائل $\Omega^\varepsilon_p$ والمجال الصلب المنفصل $\Omega^\varepsilon_s$ للتعامل مع قيد التباعد وإعادة بناء الضغط. ومن الأهمية بمكان أن معايير هذه المؤثرات مستقلة عن $\varepsilon$.
  • مؤثرات التقييد والتمديد: تُستخدم لتمديد بارامتر الترتيب $\phi^\varepsilon$ وسرعة المتجهة $u^\varepsilon$ من المجال المسامي إلى المجال الكامل $\Omega$ مع الحفاظ على حدود موحدة في فضاءات سوبوليف.
• التعامل مع حد المصدر: على عكس نماذج كان-هيلارد الكلاسيكية، فإن حد المصدر $G(\phi^\varepsilon)$ يكسر حفظ الكتلة. استنتج المؤلفون تقديرات محددة للمتوسط المكاني لـ $\phi^\varepsilon$ وللجهد الكيميائي للتعامل مع غياب خصائص المتوسط الصفري، مستفيدين من البنية المحددة لـ $G$ (bi-Lipschitz).
• تقديرات الطاقة: يثبتون قانون تبديد طاقة موحد، حيث تكون الطاقة الإجمالية محدودة بـ $C\sqrt{\lambda_\varepsilon}$.

ب. التجانس الدوري (المستوى العياني)

يقوم المؤلفون بعملية الحد $\varepsilon \to 0$ باستخدام طريقة التفكيك الدوري (Periodic Unfolding Method) ($T^\varepsilon_*$).

• التقارب ثنائي المقياس (Two-Scale Convergence): يحددون الحقول العيانية ($u, \phi, \mu$) والمصححات المجهرية ($\hat{u}, \hat{\phi}, \hat{\mu}$) المعرفة على الخلية المرجعية $Y_p$.
• التمييز بين الحالات: يعتمد الحد المجانس بشكل حاسم على القيمة الحدية $\lambda = \lim_{\varepsilon \to 0} \lambda_\varepsilon$:
  1. تلاشي الشعيرية ($\lambda = 0$): تصبح قوة الشعيرية مهملة بالنسبة للحدود الأخرى في معادلة الزخم الموزونة.
  2. النظام المتوازن ($\lambda \in (0, +\infty)$): تظل قوة الشعيرية، والتبديد اللزوجي، والتحرك غير الخطي في حالة توازن.

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

1. وجود الحلول الضعيفة

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

2. اشتقاق نموذجين فعالين متميزين

بناءً على تحجيم قوة الشعيرية، يتم اشتقاق نموذجين مختلفين على المستوى العياني:

• الحالة الأولى: $\lambda = 0$ (حد ستوكس-كان-هيلارد)

  • تفقد معادلة الزخم العيانية حد الحمل غير الخطي $(u \cdot \nabla)u$.
  • تفقد معادلة مجال الطور حد النقل الإعلاتي $u \cdot \nabla \phi$.
  • النتيجة: نظام ستوكس-كان-هيلارد (SCH) حيث يتم توجيه التدفق فقط بواسطة قوة الشعيرية المجانسة وقوى الجسم الخارجية، دون وجود عطالة أو حمل عياني.
  • المعادلة (1.5): $\partial_t u - \nabla \cdot (A^{hom} D(u)) + \nabla \tilde{p} + \phi C^{hom} \nabla \mu = g$.

• الحالة الثانية: $\lambda \in (0, +\infty)$ (حد نافييه-ستوكس-كان-هيلارد)

  • يظل الحمل غير الخطي والحدود الإعلاتية باقية في الحد، موازنة بـ $\sqrt{\lambda}$.
  • النتيجة: نظام نافييه-ستوكس-كان-هيلارد (NSCH) كامل على المقياس العياني.
  • المعادلة (1.7): $\partial_t u + \sqrt{\lambda}(u \cdot \nabla)u - \nabla \cdot (A^{hom} D(u)) + \nabla \tilde{p} + \phi C^{hom} \nabla \mu = g$.
  • تتضمن معادلة مجال الطور أيضاً الحد الإعلاتي $\sqrt{\lambda} u \cdot \nabla \phi$.

3. التنسورات المجانسة ومشكلات الخلية

يحدد المؤلفون صراحةً التنسورات الفعالة (المجانسة) $A^{hom}$ و $B^{hom}$ و $C^{hom}$ عبر مشكلات الخلية على المسام المرجعية $Y_p$. تلتقط هذه التنسورات تأثير الهندسة المجهرية وشروط الانزلاق الحر على خصائص النقل العيانية.

4. تقارب الطاقة

تثبت الورقة تقارب الطاقة المجهرية الإجمالية إلى دالة طاقة مجانسة. تحقق الطاقة الحدية قانون تبديد مماثل، مما يؤكد الاتساق الديناميكي الحراري للنماذج العيانية المستمدة.

4. الأهمية

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

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