Droplet Spreading on Nanowall-Textured Hydrophobic Surfaces: A Nonlinear Reduced-Order Model
تقدم هذه الورقة نموذجاً مختزلاً غير خطي بدرجة حرية واحدة ينجح في رصد التشوه الشعاعي وعوامل الانتشار القصوى لقطرات الماء المرتطمة بكل من الأسطح السيليكونية الكارهة للماء الملساء وتلك ذات الأسطح النانوية الجدارية، مما يثبت أن السطح ذو الملمس يقلل من الانتشار بشكل كبير مقارنة بالأسطح الملساء مع تحديد الصلابة غير الخطية ومقدار التثبيت كعوامل تحكم رئيسية.
عندما تسقط قطرة ماء على سطح ما، فإنها لا تكتفي بمجرد الارتطام والتوقف؛ بل تؤدي رقصة قصيرة ومعقدة من الطاقة. يصطدم اندفاع القطرة للأمام بالسطح، مما يؤدي إلى تسطح السائل وتحويله إلى قرص عريض ورقيق قبل أن تحاول قوة التوتر السطحي —وهي نفس القوة التي تجعل الماء يتكور على شكل حبات— سحبها مرة أخرى لتصبح كروية. ويعتمد ما إذا كانت القطرة ستنتشر على اتساع أم ستبقى صغيرة على ذلك التجاذب بين السرعة الأولية ومقاومة السطح الذي تصطدم به. ففي حالة وجود سطح أملس تماماً وطارد للماء، يمكن للقطرة أن تنتشر بشكل كبير. ولكن عندما يكون ذلك السطح مغطى بهياكل مجهرية دقيقة، مثل غابة من الجدران النانوية، تتغير القصة. فهذه الأنسجة يمكن أن تحبس الهواء، أو تثبت حافة القطرة في مكانها، أو تخلق احتكاكاً إضافياً، وكل هذه العوامل تغير مدى انتشار الماء. إن فهم هذا السلوك أمر بالغ الأهمية لكل شيء، بدءاً من تصميم النوافذ ذاتية التنظيف وصولاً إلى تحسين الطباعة بنفث الحبر، ومع ذلك، يظل التنبؤ بدقة بكيفية تفاعل القطرة مع مثل هذه الأسطح المعقدة تحدياً صعباً للعلماء.
في دراسة حديثة، سعى الباحثون إلى تبسيط هذه المشكلة الفيزيائية المعقدة إلى نموذج يمكن إدارته. فقد ركزوا على قطرات الماء، وتحديداً قطرات بحجم 20 ميكرولتر، تسقط من ارتفاع 20 مليمتراً على نوعين من أسطح السيليكون: أحدهما أملس وكاره للماء (طارد له)، والآخر مغطى بجدران نانوية بتباعد قدره 130 نانومتر فقط. لم يكن الهدف هو محاكاة كل جزيء ماء على حدة أو الشكل الدقيق للهواء المحبوس بين الجدران، لأن ذلك يتطلب قدرة حوسبية هائلة. بدلاً من ذلك، بنى الفريق نموذجاً "مختزلاً". فكر في الأمر كإنشاء خريطة مبسطة لمدينة تظهر فقط الطرق الرئيسية وإشارات المرور، متجاهلةً كل شارع جانبي ومبنى، وذلك فقط للتنبؤ بالوقت الذي ستستغرقه الرحلة. لقد تعامل نموذجهم مع القطرة بأكملها كنقطة واحدة متحركة يمكنها التمدد والتقلص، محكومة بعدة قواعد أساسية حول كيفية مقاومة السطح لهذا التمدد.
اختبر الباحثون نموذجهم أولاً على السطح الأملس. ووجدوا أنه إذا عاملوا قطرة الماء كزنبرك بسيط يقاوم التمدد في خط مستقيم، فإن النموذج سيتنبأ بأن القطرة ستنتشر بشكل مفرط، بما يتجاوز بكثير ما تظهره التجارب فعلياً. ولإصلاح ذلك، أضافوا قاعدة مفادها أن السطح يصبح أكثر صلابة كلما تمددت القطرة أكثر، تماماً مثل الشريط المطاطي الذي يصبح أصعب في السحب كلما مددته أكثر. ومع هذا التعديل، طابقت تنبؤات النموذج للسطح الأملس البيانات التجريبية بشكل مثالي تقريباً، حيث أظهرت معامل انتشار أقصى يبلغ حوالي 2.25. وهذا يعني أن قطر القطرة نما إلى حوالي 2.25 ضعف حجمها الأصلي قبل أن يتوقف عن التمدد.
بعد ذلك، طبق الفريق نفس النموذج على السطح ذي الملمس المكون من الجدران النانوية. وهنا، أدخلوا عنصراً جديداً: قوة "التثبيت". ويمثل هذا الجدران المجهرية التي تمسك بحافة قطرة الماء، مما يجعل انتشار القطرة أكثر صعوبة. وحسب النموذج أن هذه المقاومة الإضافية ستحد بشكل كبير من الانتشار. وكانت النتيجة هي معامل انتشار أقصى قدره 1.810، مما يعني أن القطرة على السطح ذي الملمس بقيت أصغر بكثير مقارنة بالسطح الأملس. وتوافقت هذه النتيجة مع الاتجاهات العامة الملحوظة في التجارب السابقة، مما أكد أن الملمس نجح في كبح القطرة. كما تحقق الباحثون من كيفية سلوك النموذج عند تغيير حجم القطرة أو الارتفاع الذي سقطت منه. وفي كل حالة، تنبأ النموذج بشكل صحيح بأن القطرة ستنتشر بشكل أكبر إذا اصطدمت بالسطح بقوة أكبر، لكنها ستنتشر دائماً بشكل أقل على السطح ذي الملمس مقارنة بالسطح الأملس.
ومع ذلك، كان المؤلفون واضحين جداً بشأن ما لا يستطيع نموذجهم القيام به. فبسبب تبسيطهم للقطرة إلى نقطة واحدة متحركة، لا يمكن لنموذجهم التنبؤ بما يحدث بعد توقف القطرة عن الانتشار. فهو لا يستطيع إخبارنا ما إذا كانت القطرة سترتد عن السطح، أو ما مدى ارتفاع قفزتها، أو متى ستنفصل عنه بالضبط. كما أنه لا يأخذ في الاعتبار الفيزياء المعقدة للهواء المحبوس تحت القطرة أو تسرب الماء في الفجوات الصغيرة بين الجدران النانوية. النموذج هو أداة حصرياً لفهم أقصى عرض تصل إليه القطرة خلال الاصطدام الأولي. كما أشار الباحثون إلى أنه بينما عمل نموذجهم بشكل جيد للظروف المحددة التي اختبروها، فإن الطريقة التي حسبوا بها قوة "التثبيت" كانت تقديراً بناءً على مرجع ثابت، وقد لا تلتقط بدقة كل تفاصيل الأنسجة السطحية المختلفة.
في نهاية المطاف، يوفر هذا العمل طريقة واضحة ومنخفضة التكلفة حوسبياً لفهم حدود مدى انتشار قطرة الماء على سطح خشن وطارد للماء. فمن خلال تجريد التفاصيل المعقدة لديناميكا السوائل والتركيز على القوى الرئيسية المؤثرة، أنشأ الباحثون أداة موثوقة تلتقط السلوك الأساسي للنظام. لقد أثبتوا أنه بينما تحدد سرعة القطرة وصلابة السطح مدى انتشارها، فإن الملمس الدقيق للسطح يعمل كمكبح، مما يبقي القطرة أصغر حجماً. يساعد هذا الفهم العلماء والمهندسين على تصميم أسطح أفضل لاحتياجات محددة، مع معرفة أن إضافة الملمس المجهري هي وسيلة فعالة للتحكم في كيفية سلوك السوائل، حتى لو ظل المشهد الكامل لحركة القطرة معقداً للغاية بحيث لا يمكن لهذا النموذج البسيط حله.
ملخص تقني: انتشار القطرات على الأسطح الكارهة للماء ذات التضاريس الجدارية النانوية
بيان المشكلة يتضمن اصطدام القطرات بالأسطح غير المبللة تحولاً عابراً معقداً للطاقة الحركية الانتقالية إلى طاقة بينية، وحركة داخلية، وتبدد. وبينما تلخص قوانين القياس الكلاسيكية (مثل βmax∼We1/4) اتجاهات أقصى تشوه، إلا أنها لا تحدد بشكل فريد كيف تعدل التضاريس السطحية توازنات القوة والطاقة الأساسية. وتحديداً، تُدخل التضاريس السطحية تأثيرات متنافسة: تغيرات في هندسة البلل المتوازن (حالات وينزل مقابل كاسي-باكستر)، وحواجز طاقة تعيق خط التماس (التثبيت/التعليق)، والقدرة المحتملة لتغلغل السائل داخل التضاريس مما يغير التبدد والطاقة المستردة.
الإطار المحدد الذي يتم تناوله هو اصطدام قطرات الماء بأسطح السيليكون الملساء والمعدلة بجدران نانوية بناءً على البيانات التجريبية من Yilbas وآخرون [11]. ويتمثل التحدي في تحديد أبسط نظام ديناميكي شعاعي قادر على إعادة إنتاج اتجاه واعتماد أقصى انتشار على التضاريس دون اللجوء إلى محاكاة متعددة الأطوار كاملة الدقة ومكلفة حاسوبياً (والتي تتطلب حل مقاييس طول تمتد عبر أربعة رتب مقدارية).
المنهجية صاغ المؤلفون نموذجاً مختزلاً غير خطي بـ درجة حرية واحدة (1-DOF) يعمل كبديل شعاعي مجمع. هذا النموذج ليس محاكاة حجم السائل (VOF) أو محاكاة حقل الطور (phase-field) دقيقة.
الديناميكيات الحاكمة: يتتبع النظام امتداداً شعاعياً فعالاً q(t)، مع قيمة مرجعية q0=D0/2. معادلة الحركة هي: mq¨=−kcap(q−q0)[1+κs(q0q−q0)2]−cviscq˙−Fpintanh(vsq˙)
الصلابة غير الخطية: يوفر حد تكعيبي (κs) تصلباً يعتمد على السعة لمنع الامتثال غير الواقعي عند التشوهات الكبيرة.
التخميد اللزوجي: يُنمذج بمعامل يعتمد على رقم رينولدز.
المقاومة الحركية: حد منتظم يعارض الحركة، ويقترب من مقدار Fpin عند السرعات العالية.
قاعدة التثبيت الهجينة: نظراً لأن قوة المقاومة المنتظمة تتلاشى عند السرعة صفر (مما يمنع التثبيت الاستاتيكي)، فقد تم فرض حدث انتقال من الحركة إلى السكون صراحةً. تنتهي المسارات عندما تنخفض السرعة تحت نطاق التقاط (vstick) وتكون قوة دفع الشعيرية غير كافية للتغلب على عتبة سعة استاتيكية (ηsFpin).
إغلاق التثبيت: يتم تقدير مقدار التثبيت Fpin باستخدام نموذج قوة على مستوى المسام مُدرج بمقياس عدد مواجهات مرجعي ثابت. ومن الجدير بالذكر أنه تحت افتراض أن القطر الهيدروليكي يساوي تباعد التضاريس (DH=2a)، فإن مقدار التثبيت غير البعدي يصبح مستقلاً عن كل من التباعد وحجم القطرة.
المعايرة والتحقق من الصحة:
تم ضبط المعامل التكعيبي κs ليتوافق مع نقطة بيانات رسومية واحدة (سطح أملس، V=20μL، h=20mm).
يتم التنبؤ باستجابة السطح ذي التضاريس باستخدام معاملات ثابتة، دون مزيد من الضبط.
تحليل الحساسية: تم إجراء حساسية حتمية محلية على معاملات الإغلاق (κs,Fpin,ζ0,vs) والإعدادات العددية للتمييز بين الضوابط الفيزيائية والآثار العددية.
المساهمات الرئيسية
تسلسل النماذج: يضع البحث تسلسلاً للنماذج (خطي غير مثبت، غير خطي غير مثبت، غير خطي مثبت) لعزل الأدوار المحددة للصلابة غير الخطية ومقاومة التثبيت.
الديناميكيات الهجينة الصريحة: يقدم قاعدة رياضية للانتقال من الحركة إلى السكون للتعامل مع التثبيت الاستاتيكي، مع الإقرار بأن التنعيم المنتظم وحده لا يمكنه الحفاظ على حالة مثبتة.
القابلية للتعريف والشفافية: يصنف العمل المعلمات صراحةً كقيم أدبية، أو مدخلات مُعاد بناؤها، أو معاملات مضبوطة، أو إغلاقات مفترضة. ويوضح أن الاتفاق عند حالة واحدة لا يمكنه التحقق من جميع الإغلاقات في آن واحد.
تحديد النطاق: النموذج محدد بصرامة للتنبؤ بـ عامل الانتشار الشعاعي الأقصى (βmax). وهو يستبعد صراحةً التنبؤ بالارتداد الرأسي، أو تحولات حالة البلل، أو الانفصال الفيزيائي.
النتائج
الحالة الأساسية (V=20μL, h=20mm):
النموذج الخطي غير المثبت يبالغ في تقدير الانتشار (βmax=3.955).
النموذج غير الخطي غير المثبت (المعاير) يطابق الهدف للسطح الأملس (βmax=2.253).
النموذج غير الخطي ذو التضاريس يتنبأ بانتشار منخفض (βmax=1.810)، وهو ما يتوافق مع الاتجاه النوعي لتقليل الانتشار على أسطح الجدران النانوية الملاحظ من قبل Yilbas وآخرون.
نقل الحالة: مع تثبيت المعلمات، يتنبأ النموذج بأن أقصى انتشار يزداد رتيباً مع رقم ويبر المُعاد بناؤه (We) لكلا السطحين. ويُظهر السطح ذو التضاريس باستمرار βmax أقل من السطح الأملس عبر نطاقات الحجم وارتفاع الإطلاق المستقصاة.
الحساسية:
تم تحديد الصلابة غير الخطية (κs) ومقدار التثبيت (Fpin) كضوابط رئيسية على أقصى انتشار.
مقدار التثبيت هو الضابط المهيمن الخاص بالسطح ذي التضاريس.
المعلمات العددية (تسامح الحل، تحديد الحدث، سرعة التنعيم الحركي) لها تأثيرات ضئيلة على أقصى انتشار متوقع عند الدقة المذكورة.
زمن النهاية: يتسم الزمن الذي تنتهي عنده المسارات بحساسية عالية تجاه عامل السعة الاستاتيكية (ηs) ومقدار التثبيت بالقرب من عتبات التبديل، مما يشير إلى أن زمن النهاية هو تشخيص لتصنيف الحالة وليس تنبؤاً فيزيائياً قوياً.
الأهمية والادعاءات يدعي البحث بتواضع تقديم بديل شعاعي قابل للتكرار ومنخفض التكلفة حاسوبياً يحدد ما يمكن وما لا يمكن استنتاجه من معلومات أقصى انتشار المتاحة.
نجح في إعادة إنتاج الاتجاه النوعي بأن التضاريس الجدارية النانوية تقلل من أقصى انتشار مقارنة بالأسطح الملساء.
أوضح أن تقليل الانتشار مدفوع بالتأثير المشترك لعمل التثبيت الخارجي وقاعدة السعة الاستاتيكية الهجينة.
يؤكد المؤلفون أن النموذج لا يتنبأ بالارتداد الرأسي، أو معاملات الاستعادة الطبيعية، أو أوقات التلامس الفيزيائية. إن "الصفر" في بروكسي الاسترداد الشعاعي للحالة ذات التضاريس هو عرف تصنيفي لعدم العودة الشعاعية، وليس تنبؤاً بارتفاع ارتداد تجريبي صفري.
يعمل العمل كخط أساس صريح للنماذج المستقبلية ذات الإحداثيات الثنائية أو النماذج متعددة الأطوار دقيقة الدقة، مسلطاً الضوء على قيود النماذج أحادية الإحداثي في التقاط تحولات حالة البلل المعقدة وآليات استرداد الطاقة.