A Thermodynamically Consistent Free Boundary Model for Two-Phase Flows in an Evolving Domain with Bulk-Surface Interaction

فيما يلي ملخص تقني مفصل للورقة البحثية بعنوان "نموذج حدود حرة متسق ديناميكيًا حراريًا لتدفقات الطورين في نطاق متطور مع التفاعل بين الحجم والسطح" من إعداد باتريك كنوف ويادونج ليو.

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

تتناول الورقة البحثية النمذجة الرياضية لتدفقات الطورين (خليط من سائلين لزوجيين غير قابلين للانضغاط) التي تحدث داخل نطاق متطور زمنياً $\Omega(t)$. وخلافاً للنماذج التقليدية حيث يكون النطاق ثابتاً، ينظر هذا العمل في السيناريوهات التي يتحرك فيها حد النطاق $\Gamma(t) = \partial\Omega(t)$ ويتشوه استجابةً لمجال سرعة المائع.

تشمل التحديات الرئيسية التي تم تناولها ما يلي:

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

2. المنهجية

استخدم المؤلفون طريقتين رياضيتين مختلفتين لاستخلاص المعادلات الحاكمة، وكلاهما يبدأ من قوانين توازن الكتلة المحلية:

أ. نهج مضاعف لاجرانج (الحالة العامة)

يُستخدم هذا النهج لاستختاج النموذج الكامل دون افتراضات مقيدة على الكثافة أو تدفق الكتلة.

• تبدد الطاقة المحلي: صاغ المؤلفون متباينات تبدد الطاقة المحلية لكل من الحجم والسطح.
• الافتراضات التكوينية: قدموا مضاعفات لاجرانج (الجهد الكيميائي $\mu, \theta$ وضغط السطح $q$) لفرض القيود (عدم الانضغاط، وحفظ الكتلة).
• تحديد التدفق: من خلال ضمان تحقق متباينة التبدد، حددوا العلاقات التكوينية لموترات الإجهاد، وتدفقات الكتلة، والشروط الحدية. يتيح ذلك اشتقاقاً عاماً يتضمن كثافات غير متطابقة وانتقال كتلة بين الحجم والسطح.

ب. النهج المتغير الطاقي (EnVarA) (الحالة الخاصة)

يُطبق هذا النهج تحت افتراضات تبسيطية محددة: كثافات متطابقة ($\tilde{\rho}_1 = \tilde{\rho}_2$) وعدم وجود تدفق كتلة بين الحجم والسطح.

• مبدأ الفعل الأدنى: استُخدم لاستخلاص القوى العطالية والحافظة من الدوال الوظيفية للحركة والطاقة الحرة.
• مبدأ التبدد الأقصى (مبدأ أونسجر): استُخدم لاستخلاص القوى المبددة (اللزوجة، الاحتكاك) من دالة تشتت رايلي.
• توازن القوة: يتم الحصول على المعادلات النهائية من خلال موازنة مجموع القوى العطالية، والحافظة، والمبددة لتساوي الصفر.

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

تقدم الورقة نظام نافييه-ستوكس-كاهن-هيليارد (NS-CH) مبتكر مقترن بشروط حدية ديناميكية.

النظام الحاكم

يتكون النموذج من:

1. ديناميكيات الحجم:
   • معادلة نافييه-ستوكس: تصف توازن الزخم للخليط بكثافة تعتمد على الطور $\rho(\phi)$ ولزوجة $\nu(\phi)$. وهي تتضمن حداً إضافياً للإجهاد ناتجاً عن التوتر السطحي (إجهاد كورتويغ).
   • معادلة كاهن-هيليارد الحملية: تصف تطور حقل الطور في الحجم $\phi$، المقترن بسرعة المائع.
2. ديناميكيات السطح:
   • معادلة كاهن-هيليارد السطحية: شرط حدي ديناميكي يصف تطور حقل الطور السطحي $\psi$.
   • الاقتران: يتم اقتران معادلات الحجم والسطح عبر شروط حدية تسمح بانتقال المواد (يتم التحكم فيه بواسطة المعاملات $K$ و $L$).
3. الشروط الحدية:
   • انزلاق نافييه المعمم: يستبدل شرط عدم الانزلاق للسماح بخطوط التلامس المتحركة. السرعة المماسية تتناسب مع الإجهاد المماسي بالإضافة إلى حد قوة متعلق بتدرجات التوتر السطحي.
   • الشروط الحدية الديناميكية: يتم تحديد السرعة العمودية للحد من خلال المكون العمودي لسرعة المائع ($u \cdot n = V$)، مما يجعلها مسألة حدود حرة.
   • ديناميكيات زاوية التلامس: تسمح الشروط الحدية لـ $\phi$ و $\mu$ (المتحكم بها بواسطة المعاملات $K, L$) لزاوية التلامس بالتطور ديناميكياً بدلاً من كونها ثابتة عند 90 درجة.

الخصائص الفيزيائية الرئيسية

يستوفي النموذج المستمد ما يلي:

• حفظ الكتلة: يتم حفظ إجمالي الكتلة؛ إذا كان $L=\infty$، تُحفظ كتل الحجم والسطح بشكل منفصل.
• تبدد الطاقة: يحقق النظام قانون تبدد الطاقة العالمي، مما يثبت الاتساق الديناميكي الحراري. إجمالي الطاقة (الحركية + طاقة الحجم الحرة + طاقة السطح الحرة) يتناقص بمرور الوقت بسبب اللزوجة، والاحتكاك، والانتشار.
• حفظ الحجم والمساحة: نظراً لقيود عدم الانضغاط ($\text{div } u = 0$ و $\text{div}_\Gamma u = 0$)، يتم الحفاظ على حجم النطاق ومساحة الحد.

4. النتائج

• التحقق من الاشتقاق: نجح المؤلفون في اشتقاق نفس النظام باستخدام طريقتين مختلفتين، مما يؤكد قوة النموذج.
• التعميم: يُظهر النموذج أنه يعمم الأدبيات الموجودة:
  • عند ضبط $K=L=\infty$، يستعيد النماذج ذات زوايا التلامس الثابتة عند 90 درجة وعدم وجود انتقال كتلة سطحية.
  • عند تثبيت حد النطاق ($V=0$)، يستعيد نماذج NS-CH القياسية مع الشروط الحدية الديناميكية الموجودة في أعمال سابقة (مثل [50]).
  • يتم استعادة نموذج AGG (Abels-Garcke-Grün) كحالة خاصة في نطاق ثابت.
• تحليل المعاملات: تفصل الورقة كيف تتحكم المعاملات $K$ و $L$ في السلوك الفيزيائي:
  • $K=0$: يفرض علاقة خطية بين حقول الطور في الحجم والسطح (توازن سريع).
  • $L \in (0, \infty)$: يسمح بمعدل محدود لانتقال الكتلة (الامتصاص/الامتزاز).
  • $L=\infty$: يفرض صفراً لتدفق الكتلة عبر الحد.

5. الأهمية

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

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