Deformation and orientation of a capsule with viscosity contrast in linear flows: a theoretical study

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

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

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

• تباين اللزوجة ($\lambda$): حيث يُسمح لنسبة لزوجة السائل الداخلي إلى لزوجة السائل الخارجي بأن تختلف عن الواحد الصحيح.
• التوتر السطحي ($\sigma$): وجود التوتر السطحي، وهو أمر ذو صلة بالأغشية المسامية أو أنظمة الزيت في الماء.
• صلابة الانحناء ($\kappa$): مقاومة الغشاء لتغيرات الانحناء، وهو أمر حاسم لتحليل الاستقرار (البعج/التجعد).
• القوانين التكوينية: يتم نمذجة الاستجابة المرنة للغشاء باستخدام ثلاثة قوانين متميزة: سكالاك (Skalak)، وهوكي (Hookean)، ونيو-هوكيان (Neo-Hookean).

الهدف الأساسي هو اشتقاق تعبيرات تحليلية لشكل الكبسولة (التشوه) وزاوية توجهها حتى الرتبة الثانية في رقم كابيلاري ($Ca$)، وهو نظام توقفت فيه الدراسات السابقة غالباً عند الرتبة الأولى.

2. المنهجية

يستخدم المؤلفون نهج نظرية الاضطراب (Perturbation theory) ضمن نظام تدفق ستوكس (رقم رينولدز منخفض).

• المعادلات الحاكمة: يوصف النظام بمعادلات ستوكس لكل من السوائل الداخلية والخارجية، مع ربطها بشروط حدودية عند واجهة الكبسولة. وتشمل هذه الشروط استمرارية السرعة، وتوازن الإجهاد (الذي يدمج القوى الهيدروديناميكية، والمرنة، والانحناء، والكبيلارية)، والشرط الكينماتيكي.
• القوانين التكوينية:
  • المرونة: يتم اشتقاق إجهاد الغشاء من دوال طاقة الانفعال ($w$) الخاصة بنماذج سكالاك، وهوكي، ونيو-هوكيان، مما يقدم معاملات مثل معامل القص ($G_s$)، ونسبة بواسون ($\nu$)، ومعامل تصلد الانفعال ($C$).
  • الانحناء والكبيلارية: تُشتق القوى من طاقة هيلفريک للانحناء وطاقة السطح، مما يقدم أرقاماً لا بُعدية لصلابة الانحناء ($B = \kappa / G_s R^2$) ونسبة المرونة الكبيلارية ($\Sigma = \sigma / G_s$).
• استراتيجية التوسع:
  • يتم توسيع نصف قطر الكبسولة كـ $r = 1 + F^{(1)} + F^{(2)}$، حيث $F^{(1)}$ هو التشوه من الرتبة الأولى، و$F^{(2)}$ هو التصحيح من الرتبة الثانية.
  • يتم توسيع جميع الكميات الفيزيائية (السرعة، الضغط، الإجهاد) كقوى من رقم كابيلاري $Ca$.
  • يستخدم الحل التوافقيات الصلبة (حل لامب) لحل معادلات ستوكس. ويتم توسيع تينسورات التشوه باستخدام التوافقيات الكروية من الرتبة الثانية (لتشوه الرتبة الأولى) والرتبة الرابعة (لتأثيرات الرتبة الثانية).
• الحل التحليلي: يقوم المؤلفون بحل نظام المعادلات الخطية الناتج عن طريق التكامل عبر سطح الكرة وتطبيق مؤثرات التماثل. وهذا يسمح بفصل المعادلات الخاصة بتينسور التشوه ($F$) وتينسور التوجه ($K$).
• التحقق من الصحة: تم التحقق من النتائج التحليلية مقابل المحاكاة العددية التي أجريت باستخدام كود مخصص يعتمد على طريقة التكامل الحدودي (BIM) وطريقة العناصر المحدودة (FEM).

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

1. نظرية الرتبة الثانية للتوجه: تشتق الورقة نظير تشافي-برينر (Chaffey-Brenner) للكبسولات. وبينما تتنبأ نظرية الرتبة الأولى بزاوية توجه ثابتة (غالباً $\pi/4$ في تدفق القص)، فإن نظرية الرتبة الثانية توفر التصحيح لهذه الزاوية ($\delta$)، وتوضح كيفية اعتمادها على تباين اللزوجة، وصلابة الانحناء، والتوتر السطحي.
2. تعميم القوانين التكوينية: تقدم الدراسة تعبيرات تحليلية موحدة لثلاثة نماذج مرنة مختلفة (سكالاك، هوكي، نيو-هوكيان) ضمن إطار عمل واحد يتضمن الانحناء والتوتر السطحي.
3. استقلال التشوه عن تباين اللزوجة في الرتبة الرائدة: يؤكد المؤلفون أنه، على غرار نتائج بارثس-بيسيل وراليسون (1981) الكلاسيكية، فإن التشوه في الرتبة الرائدة مستقل عن تباين اللزوجة $\lambda$ عند غياب التوتر السطحي والانحناء. ومع ذلك، يوضحون أنه في الرتب الأعلى، تؤثر $\lambda$ على التوجه.
4. دور حدود الرتب العليا: تثبت الدراسة أن مساهمة الرتبة الثانية في التشوه (تغير الشكل) هي صفر لجميع القوانين التكوينية. وهذا يعني أن السلوك الميكانيكي غير الخطي لشكل الكبسولة يظهر فقط في الرتبة الثالثة، مما يفسر سبب إظهار البيانات التجريبية غالباً لاستجابة خطية عبر نطاق واسع من $Ca$.
5. الترميز بالمعاملات اللابعدية: يتم التعبير عن النتائج بدلالة أربعة أرقام رئيسية لا بُعدية:
   • رقم كابيلاري ($Ca$)
   • تباين اللزوجة ($\lambda$)
   • رقم الكابيلارية المرنة ($\Sigma$)
   • صلابة الانحناء اللابعدية ($B$)

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

• التشوه (معامل تايلور):
  • في النظام الخطي (الرتبة الرائدة)، يتناسب التشوه $D$ طردياً مع $Ca$.
  • وجود التوتر السطحي ($\Sigma$) وصلابة الانحناء ($B$) يقلل من التشوه. وفي حدود التوتر السطحي العالي أو صلابة الانحناء العالية، يصبح التشوه مستقلاً عن معامل القص $G_s$.
  • تعبيرات التشوه لنماذج سكالاك، وهوكي، ونيو-هوكيان متميزة ولكنها تتقارب تحت خرائط معاملات محددة (على سبيل المثال $\nu = C/(1+C)$).
• زاوية التوجه ($\phi_{TT}$):
  • تنحرف زاوية الميل عن $\pi/4$ بمقدار يتناسب مع $Ca$.
  • على عكس التشوه، يعتمد التوجه على تباين اللزوجة $\lambda$ حتى في رتبة التصحيح الأولى.
  • تنجح المعادلات المستمدة لـ $\phi_{TT}$ في استعادة علاقة تشافي-برينر لقطرات السائل في حد انعدام مرونة الغشاء.
• المقارنة مع النتائج العددية:
  • تظهر التوقعات التحليلية لكل من التشوه والتوجه توافقاً ممتازاً مع محاكاة BIM العددية عبر مجموعة واسعة من المعاملات ($\lambda, \Sigma, B$).
  • يظل هذا التوافق قائماً حتى في الحدود التقاربية (التوتر السطحي العالي أو صلابة الانحناء العالية)، مما يؤكد قوة توسع الاضطراب.

5. الأهمية والتداعيات

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

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