Non-Monotonic Marangoni Suppression of Hydrodynamic Coarsening in Bicontinuous Liquid-Liquid Phase Separation

إليك ملخص تقني مفصل للورقة البحثية بعنوان "التثبيط غير الرتيب لظاهرة مارانغوني في التخشين الهيدروديناميكي في الفصل الطوري السائل-السائل ثنائي الاستمرارية" (Non-Monotonic Marangoni Suppression of Hydrodynamic Coarsening in Bicontinuous Liquid-Liquid Phase Separation) من إعداد ليو وآخرون.

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

غالباً ما يؤدي الفصل الطوري السائل-السائل (LLPS) في السوائل متعددة المكونات إلى تكوين نطاقات ثنائية الاستمرارية تخضع لعملية تخشين هيدروديناميكي (hydrodynamic coarsening)، حيث تزداد أحجام النطاقات بمرور الوقت. وبينما يُعرف أن إضافة المواد الخافضة للتوتر السطحي القابلة للذوبان تعمل على تأخير هذه العملية، إلا أن الدور الميكانيكي المحدد لـ إجهادات مارانغوني (Marangoni stresses) — وهي الإجهادات المماسية الناتجة عن تدرجات التوتر السطحي — مقابل مجرد الخفض البسيط لمتوسط التوتر السطحي، لا يزال غير مفهوم بشكل كافٍ.

هناك سؤال جوهري لم يُحل بعد، وهو كيف ينظم نقل المواد الخافضة للتوتر السطحي القابلة للذوبان — والذي يحكمه التوازن بين الانتشار والحمل (والذي يتميز برقم بيكليه للمادة الخافضة للتوتر السطحي، $Pe_\psi$) — إجهادات مارانغوني هذه. وتحديداً، من غير الواضح لماذا يعتمد تثبيط التخشين على علاقة غير رتيبة مع $Pe_\psi$، بدلاً من الزيادة أو النقصان ببساطة مع حركية المادة الخافضة للتوتر السطحي.

2. المنهجية

استخدم المؤلفون نموذج مجال طوري (phase-field model) موثق ذو معلمتين مرتبط بمعادلات نافيير-ستوكس غير القابلة للانضغاط لمحاكاة النظام.

• المعادلات الحاكمة:
  • مجال الطور ($\phi$) والمادة الخافضة للتوتر السطحي ($\psi$): يتم وصفهما بمعادلات "كان-هليارد" الحملية. وتُشتق الجهود الكيميائية من دالة طاقة حرة تتضمن حدود "جينزبرغ-لانداو"، وإنتروبيا "فلوري-هجينج"، وحدود الامتزاز/التفاعل من نوع "لانغميور".
  • الهيدروديناميكا: تتضمن معادلات نافيير-ستوكس قوى جسمية تمثل القوة الشعيرية (العمودية على الواجهة) وقوة مارانغوني (المماسية للواجهة، والمدفوعة بـ $\nabla \ln(1-\psi)$).
• المخطط العددي:
  • طريقة الحجم المحدود على شبكة منتظمة من نوع (Marker-and-Cell - MAC).
  • تم استخدام مخططات التشتت لتمثيل حدودات الحمل (advection) من الدرجة الخامسة (WENO)، بينما استخدمت طرق "أدامز-باشفوت" الصريحة (للحمل) و"أويلر/كرانك-نيكولسون" (للانتشار واللزوجة) للتقدم الزمني.
  • تستخدم طريقة الإسقاط لحل معادلات الزخم والاستمرارية المقترنة.
• أنماط المحاكاة:
  1. النموذج B: انتشار محض (بدون هيدروديناميكا).
  2. النموذج H: تخشين هيدروديناميكي (بدون مواد خافضة للتوتر السطحي).
  3. النموذج HS: تخشين هيدروديناميكي مع مواد خافضة للتوطة سطحية قابلة للذوبان واقتران كامل (بما في ذلك قوى مارانغوني).
• التحقق من الصحة: تم التحقق من صحة النموذج مقابل قوانين القياس الكلاسيكية ($R(t) \sim t^{1/3}$ للنموذج B؛ وأنظمة $t, t^{2/3}, t^{1/2}$ للنموذج H) وتمت مقارنته بالحلول التحليلية لتشوه القطرات المحملة بالمادة الخافضة للتوتر السطحي تحت تأثير القص.

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

1. التمييز الميكانيكي: تفصل هذه الدراسة بشكل قاطع بين تأثيرات خفض متوسط التوتر السطحي وبين إجهادات مارانغوني. وقد أثبتت أن تثبيط التخشين مدفوع أساساً بإجهادات مارانغوني، وليس مجرد انخفاض في متوسط التوتر السطحي.
2. الاعتماد غير الرتيب: حددت الورقة البحثية وفسرت وجود علاقة غير رتيبة بين رقم بيكليه للمادة الخافدة للتوتر السطحي ($Pe_\psi$) وتثبيط التخشين. ويحدث أقوى تثبيط عند قيمة متوسطة ($Pe_\psi = 10$)، وليس عند القيم المنخفضة ($Pe_\psi = 1$) أو العالية ($Pe_\psi = 100$).
3. الآلية المتحكم بها عبر النقل: يقترح المؤلفون آلية "يتحكم فيها النقل" حيث تعتمد فعالية تثبيط مارانغوني على المنافسة بين تعويض المادة الخافضة للتوتر السطحي (الانتشار من الكتلة) والاحتفاظ بالتدرج (الحمل الذي يحافظ على عدم التجانس).

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

أ. إجهادات مارانغوني مقابل خفض متوسط التوتر السطحي

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

ب. الاعتماد غير الرتيب على $Pe_\psi$

حللت الدراسة ثلاثة أنظمة لنقل المادة الخافضة للتوتر السطحي ($Pe_\psi = 1, 10, 100$):

• $Pe_\psi$ منخفض (الانتشار هو المهيمن، $Pe_\psi = 1$): يكون الانتشار فعالاً في تعويض المواد الخافضة للتوتر السطحي عند الواجهة، ولكنه أيضاً يعمل بسرعة على تسوية تدرجات التركيز عند الواجهة. ويؤدي غياب التدرجات المستمرة إلى ضعف إجهادات مارانغوني وضعف التثبيط.
• $Pe_\psi$ مرتفع (الحمل هو المهيمن، $Pe_\psi = 100$): يحافظ الحمل على عدم التجانس في الواجهة (التدرجات)، ولكن الانتشار يكون بطيئاً جداً لتعويض نقص المواد الخافضة للتوتر السطح من الكتلة إلى الواجهة. مما يؤدي إلى استنزاف الواجهة من المادة الخافضة للتوتر السطحي، وبالتالي ضعف إجهادات مارانغوني رغم وجود التدرجات.
• $Pe_\psi$ متوسط (التوازن الأمثل، $Pe_\psi = 10$): يحقق هذا النظام "النقطة المثالية". فالانتشار سريع بما يكفي للحفاظ على متوسط تحميل عالٍ للمادة الخافضة للتوتر السطحي عند الواجهة، بينما يكون الحمل قوياً بما يكفي للحفاظ على تدرجات تركيز كبيرة. هذا التعايش ينتج أقوى إجهادات مارانغوني وأقصى تثبيط للتخشين الهيدروديناميكي.

ج. التطور المورفولوجي

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

5. الأهمية والآثار المترتبة

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

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