Droplets sitting on thin elastic sheets: A study with the boundary element method

ملخص تقني: قطرات مستقرة على صفائح مرنة رقيقة: دراسة باستخدام طريقة العناصر الحدية

بيان المشكلة
تبحث الورقة في التفاعل المرن-السطحي (elasto-capillary) بين قطرة سائلة وصفيحة مرنة رقيقة. وخلافًا للركائز الصلبة، فإن المواد المرنة تتشوه تحت تأثير القوى السطحية، مما يؤدي إلى أشكال توازن معقدة حيث تغوص القطرة في الصفيحة، مكونةً بنية تشبه العدسة. وبينما استكشفت الدراسات السابقة أنماط التجعد، والتحرك نحو المناطق الأكثر ليونة (durotaxis)، والأشكال متباينة الخواص تحت الإجهاد، كانت هناك حاجة إلى إطار عددي موحد لمحاكاة أشكال التوازن للقطرات على الصفائح المرنة تحت ظروف حدودية وبروتوكولات شد مختلفة. وتحديدًا، يهدف المؤلفون إلى قياس مدى تأثير سمك الصفيحة، والشد متماثل المناحي (isotropic tension)، والشد أحادي المحور على مورفولوجيا (شكل) القطرة، وتوزيع الإجهاد، والخصائص البصرية.

المنهجية
قام المؤلفون بتوسيع "طريقة العناصر الحدية" (BEM) التي طوروها سابقًا لمحاكاة الديناميكيات المقترنة لقطرة وصفيحة مرنة. ويتضمن جوهر المنهجية ما يلي:

1. المعادلات الديناميكية: يتم نمذجة النظام باستخدام الديناميكيات المخمدة (overdamped dynamics) حيث ترتبط حقول السرعة لسطح القطرة ($\dot{R}_0$) والصفيحة المرنة ($\dot{R}_{sh}$) ببعضهما البعض. المعادلة الحاكمة هي نظام خطي $G\dot{R} = -\nabla_R F$، حيث $G$ هي مصفوفة الاحتكاك و $F$ هي إجمالي الطاقة الحرة.
2. صياغة الطاقة الحرة:
   • القطرة ($F_0$): تشمل التوتر السطحي للسائل-الغاز، والعمل المبذول ضد الواجهة بين الركيزة والسائل، والمحكوم بزاوية التلامس "يونغ" (Young's contact angle). تم تضمين حد التوتر الخطي ولكن تم ضبطه على الصفر لكونه مهملًا بالنسبة للصفائح المرنة المدروسة.
   • الصفيحة ($F_{sh}$): تم نمذجتها باستخدام نموذج "سكالاك" (Skalak model) للقص داخل المستوي وتمدد المساحة (المرونة الفائقة)، ونموذج "هيلفريتش" (Helfrich model) للانحناء خارج المستوي. تُشتق المعاملات المرنة (القص $\kappa_S$، والمساحة $\kappa_A$، والانحناء $\kappa_b$) من خصائص المادة الكتلية (معامل يونغ $E$ ونسبة بواسون $\nu$) وسمك الصفيحة $d$.
   • القيود ($F_{cons}$): يتم فرضها عبر مضاعفات لاجرانج للحفاظ على حجم القطرة الثابت وضمان الاتصال الهندسي بين شبكات (meshes) القطرة والصفيحة عند الواجهة.
3. مصفوفة الاحتكاك: تتضمن المصفوفة $G$ إجهادات القص اللزجة (المشتقة من موتر أوسين - Oseen tensor)، واحتكاك خط التلامس (نتيجة موفات - Moffatt's result)، واحتكاك السائل-الركيزة المحكوم بشرط انزلاق نافييه (Navier slip condition).
4. التنفيذ العددي: يتم تقسيم النظام باستخدام شبكات مثلثية تم إنشاؤها بواسطة برنامج gmsh. يستخدم المؤلفون مخطط استيفاء مرجح (باستخدام Surrogates.jl) لربط القوى والسرعات بين الشبكات المتميزة للقطرة والصفيحة. تم جعل المحاكاة غير أبعاد (non-dimensionalized) باستخدام نصف قطر خط التلامس $R_0$ والمقاييس الزمنية المميزة المستمدة من قانون "كوكس-فينوف" (Cox-Voinov law).

النتائج الرئيسية
تطبق الدراسة طريقة BEM الموسعة على ثلاثة بروتوكولات محددة:

1. الصفيحة المرنة المثبتة (Clamped):

   • شكل التوازن للقطرة هو عدسة تصبح غير متماثلة بشكل متزايد بالنسبة لمستوى خط التلامس مع زيادة سمك الصفيحة.
   • الاعتماد على السمك: مع زيادة السمك ($d$)، ينخفض الجزء السفلي من القطرة ($h_l$) بينما يزداد الجزء العلوي ($h_u$)، رغم أن الارتفاع الإجمالي يظل ثابتًا نسبيًا.
   • زوايا التلامس: تتفق زاوية التلامس المجهرية عند خط التلامس مع زاوية "يونغ"، مما يؤكد صحة النموذج. ومع ذلك، فإن زاوية التلامس العيانية (المجهرية/الكلية) (المحددة من خلال شكل العدسة العام) تختلف بشكل كبير عن زاوية "يونغ" وتعتمد على سمك الصفيحة.
   • فك الارتباط: يتحدد شكل النصف السفلي للقطرة بشكل أساسي من خلال مرونة الصفيحة وهو مستقل تقريبًا عن زاوية التلامس التوازنية ($\theta_{eq}$)، بينما تؤثر $\theta_{eq}$ بشكل أساسي على النصف العلوي.

2. الصفيحة المشدودة بشكل متماثل المناحي (Isotropically Stretched):

   • يؤدي تطبيق الشد متماثل المناحي إلى رفع القطرة، مما يقلل من تشوه الصفيحة حتى تصبح مستوية مرة أخرى عند التمدد العالي ($\epsilon \approx 0.1$).
   • ملفات الإجهاد: يحدد المؤلفون الملفات الشعاعية للإجهادات المرنة المحيطية ($T_{\phi\phi}$) والشعاعية ($T_{rr}$). تزداد الإجهادات مع التمدد، ويقع خط التلامس بالقرب من نقطة انقلاب في منحنيات الإجهاد.
   • الضبط البصري: تعمل القطرة ذات شكل العدسة كعدسة سائلة. يمكن ضبط البعد البؤري، المحسوب عبر معادلة صانع العدسات، بشكل عكسي بنسبة تقارب 15% عن طريق تغيير الشد المطبق. هذا النطاق يضاهي نطاق الضبط للعدسات في الهواتف الذكية الحديثة.

3. الصفيحة المشدودة أحادية المحور (Uniaxially Stretched):

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

الأهمية والادعاءات
تدعي الورقة أن طريقة BEM الموسعة توفر أداة عددية قوية لمحاكاة الأنظمة المرنة-السطحية، حيث تنجح في التقاط التفاعل بين ديناميكيات السوائل والتشوهات المرنة. وتتمثل المساهمات الرئيسية في:

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