Rheotaxis of microswimmers in colloid-laden channel flow

ملخص تقني: الانتحاء التياري للسباحات الدقيقة في تدفق قنوات محملة بالمواد الغروانية

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

المنهجية
يستخدم المؤلفون إطار محاكاة متوسط المقياس يجمع بين ديناميكا التصادم متعدد الجسيمات (MPCD) للسائل الخلفي وديناميكا الجزيئات (MD) للأجسام الصلبة.

• نموذج السائل: يتم محاكاة السائل الخلفي باستخدام MPCD، الذي يلتقط التفاعلات الهيدروديناميكية والتقلبات الحرارية. يتم تطبيق قوة جسم ثابتة لتوليد تدفق "بوازيه" مستوٍ بين جدارين متوازيين.
• نموذج السباح الدقيق: تُمثل السباحات الدقيقة بنموذج "سكويرمر" (squirmer) كروي، يتميز بسرعة انزلاق سطحية. ويتم تصنيفها حسب معامل السكويمر $\beta$: الدافعات ($\beta < 0$)، المحايدات ($\beta = 0$)، والجاذبات ($\beta > 0$).
• المواد الغروانية: تم إدخال مواد غروانية كروية صلبة وأحادية التشتت في القناة مع اختلاف الكسور الحجمية ($\phi$).
• التفاعلات: تتم معالجة التفاعلات بين جسيمات السائل، والسكويرمر، والمواد الغروانية، والجدران عبر قاعدة "نصف وقت الارتداد" (half-time bounce-back) للهيدروديناميكا، وجهد "ويكس-تشاندلر-أندرسن" (WCA) المقطوع للتنافر الفراغي.
• التحقق من الصحة: تم التحقق من صحة تنفيذ برنامج CUDA-C الخاص بالدراسة مقابل معاملات الانتشار النظرية للمعلقات الغروانية، ومعاملات العزم/القوة للكرات في تدفق بوازيه، والأدبيات الموجودة حول الانتحاء التياري للسكويرمر في القنوات النظيفة.

المساهمات والنتائج الرئيسية
تبحث الدراسة في توزيع الاحتمالية المشترك لموقع السكويمر ($z$) وتوجهه ($\phi_a$) عبر مختلف شدات التدفق ($V_r$) والكسور الحجمية الغروانية ($\phi$).

1. تأثير المواد الغروانية على الانتحاء التياري عند شدة تدفق منخفضة ($V_r = 0.5$):

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

2. تأثير المواد الغروانية عند شدة تدفق عالية ($V_r = 4$):

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

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

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