A reduced model for droplet dynamics with interfacial viscosity

ملخص تقني: نموذج مُخفّض لديناميكيات القطيرات مع اللزوجة البينية

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

المنهجية
يقترح المؤلفون نموذج "مافيتون-مينال الموسع" (EMM)، والذي يعمم معادلة (MM) الأصلية لتطور الموتر المورفولوجي (الشكلي) $S$ (الذي يصف شكل القطيرة كشكل إهليلجي).

• التوسيع النظري: يدمج نموذج (EMM) اللزوجات البينية من خلال أرقام بوزينيسك (Boussinesq numbers)، وهي $Bqs = \mu_s/\mu R$ و $Bqd = \mu_d/\mu R$. وبالاستناد إلى الحسابات الاضطرابية للواجهات اللزجة (تحديداً امتداد "نارسيمن" لعمل "راليشون")، اشتق المؤلفون تعبيرات عامة لمعاملات النموذج $f_1$ و $f_2$. وهذه المعاملات، التي كانت سابقاً دالة فقط في نسبة اللزوجة $\lambda$ للقطيرات النظيفة، أصبحت الآن تعتمد على $\lambda$ و $Bqs$و$Bqd$.
• التحقق العددي: لقياس دقة ونطاق صلاحية نموذج (EMM)، أجرى المؤلفون مقارنة منهجية مقابل المحاكاة كاملة الحل (FRS). تستخدم محاكاة (FRS) طريقة "الحدود المغمورة–لب التوازن" (immersed boundary–lattice Boltzmann method) لحل الديناميكيات المقترنة للسوائل الداخلية والخارجية والواجهة القابلة للتشوه دون افتراض تشوهات صغيرة.
• تحليل الخطأ: تقيم الدراسة متوسط الخطأ النسبي $\langle \epsilon \rangle_T$ بين التشوه الزمني $D(t)$ الذي يتنبأ به نموذج (EMM) ونتائج محاكاة (FRS) عبر فضاء بارامتري واسع محدد برقم كابيلاري ($Ca$) ورقمي بوزينيسك.

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

1. اشتقاق المعاملات العامة: يقدم البحث تعبيرات تحليلية صريحة لمعاملات (EMM) وهي $f_1^{(EMM)}$ و $f_2^{(EMM)}$ التي تأخذ في الاعتبار كلاً من لزوجة القص ولزوجة التمدد البينية. وتختزل هذه التعبيرات إلى معاملات (MM) الأصلية عندما تكون اللزوجات البينية صفراً.
2. توصيف الديناميكيات العابرة: تكشف الدراسة أن اللزوجة البينية تؤثر على زمن التشوه العابر بشكل غير خطي. فبينما تتصرف القطيرة ذات لزوجة القص السائدة أو لزوجة التمدد السائدة بشكل مشابه للقطيرة النظيفة فيما يتعلق بزمن الاسترخاء، فإن القطيرة التي تمتلك قيماً متقاربة لكلتا اللزوجتين تشهد تباطؤاً دراماتيكياً في ديناميكيات التشوه.
3. تحديد نطاق الصلاحية: من خلال الاختبار العددي المكثف، رسم المؤلفون حدود المناطق المحددة في الفضاء البارامتري $(Ca, Bqs, Bqd)$ التي يظل فيها نموذج (EMM) دقيقاً كمياً.

النتائج

• التشوه في الحالة المستقرة: يتنبأ نموذج (EMM) بدقة بالتشوه في الحالة المستقرة $D_s$ لمجموعة واسعة من البارامترات. ويلتقط النموذج الاتجاه غير المتوقع حيث أن زيادة لزوجة التمدد (مع وجود لزوجة قص صفرية) تزيد قليلاً من التشوه، بينما تقلل زيادة لزوجة القص (أو وجود لزوجة قص وتمدد متساويتين) من التشوه.
• السلوك العابر: ينجح النموذج في إعادة إنتاج الانتقال من الاسترخاء الرتيب إلى الاستجابات التذبذبية الملاحظة في (FRS)، خاصة عندما يكون $Bqs = Bqd$. ويزداد زمن التشوه المميز بشكل حاد فقط عندما تكون كلتا اللزوجتين البينيتين ملموستين.
• حدود الدقة:
  • بالنسبة للزوجات البينية الصغيرة ($Bqs, Bqd \ll 1$) وأرقام كابيلاري صغيرة ($Ca \le 0.2$)، يعيد نموذج (EMM) إنتاج بيانات (FRS) بأخطاء تقل عادة عن بضعة بالمائة.
  • تتدهور دقة النموذج بشكل أسرع في تكوينات اللزوجة غير المتماثلة بقوة (حيث يهيمن أحد رقمي بوزينيسك على الآخر) مقارنة بالحالة المتساوية ($Bqs = Bqd$).
  • مع زيادة $Ca$ (على سبيل المثال، $Ca = 0.5$)، تؤدي الاقترانات غير الخطية في الإجهاد البيني، والتي يقربها نموذج (EMM) القائم على فرضية التشوهات الصغيرة، إلى حدوث تباينات أكبر. وعند قيم $Ca$ المرتفعة و $Bqs$ المنخفضة، يقترب النموذج من عتبة التكسر، حيث تشير محاكاة (FRS) إلى عدم استقرار عددي بسبب التشوهات الكبيرة.

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