Universal mapping of drop impact spreading from wetting to Leidenfrost regimes
تضع هذه الورقة البحثية مخططاً عالمياً قائماً على تبدد الطاقة يتنبأ بدقة بالانتشار الأقصى للقطرات عبر كل من نظامي البلل وظاهرة ليدنفروست عبر نطاق واسع من أرقام أوهنسور وبربر ويبر، مما يتيح تحويل تنبؤات تأثير البلل إلى نظيراتها في حالة ليدنفروست.
المؤلفون الأصليون:Shushan Hu, Nan Hu, Zirui Li, Yifei Sun, Xiang Gao, Liwu Fan
عندما تسقط قطرة من سائل على سطح ما، فإنها لا تتوقف ببساطة؛ بل تتسطح، وتنتشر، ثم غالباً ما ترتد. هذا الحدث اليومي، من المطر الذي يصطدم بالنافذة إلى الحبر الذي يضرب الورق، محكوم بصراع بين الزخم الأمامي للقطرة، ولزوجة سطحها الخاص، والاحتكاك الداخلي للسائل. لقد أدرك العلماء منذ زمن طويل كيف يحدث هذا عندما يلمس السائل السطح مباشرة، وهي حالة تُعرف باسم "الترطيب" (wetting). ومع ذلك، يحدث سيناريو مختلف وأكثر دراماتيكية عندما يكون السطح شديد السخونة. في هذه الحالة، لا يلمس السائل الصلب على الإطلاق، بل يطفو بدلاً من ذلك على وسادة رقيقة غير مرئية من بخاره الخاص، وهي ظاهرة تُعرف باسم "تأثير ليدن فروست" (Leidenfrost effect). وتعد هذه الحالة حاسمة لتقنيات تتراوح من أنظمة التبريد إلى حقن الوقود، ومع ذلك ظل التنبؤ بمدى انتشار القطرة بدقة في حالة الطفو هذه لغزاً صعباً، ويرجع ذلك إلى حد كبير إلى أن فيزياء القطرة العائمة تختلف بشكل كبير عن تلك الملتصقة بالأرض.
لقد استطاع باحثون من جامعة تشجيانغ وجامعة برينستون الآن الربط بين هذين العالمين. فمن خلال الجمع بين التجارب عالية السرعة والمحاكاة الحاسوبية التفصيلية، وضعوا قاعدة عالمية تنقل سلوك الانتشار من قطرة "رطبة" إلى سلوك قطرة "ليدن فروست". ويكشف عملهم أن الفرق الجوهري بين الحالتين لا يكمن في التوتر السطحي أو سرعة القطرة، بل في كيفية فقدان السائل للطاقة. فعندما تصطدم قطرة بسطح مبلل، يعمل الاحتكاك ضد الجدار مثل مكبح، مما يبطئ الانتشار ويشتت الطاقة بالقرب من القاع. أما عندما تطفو القطرة على البخار، فإن احتكاك الجدار هذا يختفي، مما يسمح للقطرة بالانتشار لمسافة أبعد بكما. ووجد الفريق أن هذه المسافة الإضافية ليست عشوائية؛ بل هي مرتبطة مباشرة بمقدار الطاقة التي تم توفيرها بتجنب احتكاك الجدار ذاك.
ولكشف هذا الارتباط، راقب العلماء قطرات مليمترية من الماء ومزيج الجلسرين وهي تصطدم بقرص من الياقوت المصقول. في بعض التجارب، كان القرص في درجة حرارة الغرفة، مما سمح للقطرات بترطيب السطح. وفي تجارب أخرى، تم تسخين القرص إلى 550 درجة مئوية، مما خلق وسادة البخار. وقد سجلوا الاصطدامات بكاميرات تلتقط 4000 إطار في الثانية، متتبعين بدقة مدى اتساع القطرات قبل أن تتوقف عن التمدد. وقد قرنوا هذه الملاحظات بنماذج حاسوبية يمكنها حساب التدفق غير المرئي للطاقة داخل القطرة. وأظهرت عمليات المحاكاة أنه بينما تجنبت القطرات العائمة بالفعل الكبح الشديد للجدار، بدأ السائل بداخلها في الاضطراب بقوة أكبر، مما خلق نوعاً جديداً من الاحتكاك الداخلي الذي عوّض جزئياً ذلك المكسب. وبالرغم من هذا التعويض الداخلي، كان إجمالي الطاقة التي وفرتها القطرة العائمة كبيراً ويمكن التنبؤ به.
اكتشف الباحثون أن توفير هذه الطاقة يتبع نمطاً دقيقاً عبر نطاق واسع من الظروف. فقد اختبروا قطرات ذات سماكات وسرعات مختلفة، مغطين مجالاً تباينت فيه مقاومة السائل للتدفق بعامل قدره مئة، واختلفت فيه قوة الاصطدام بعشرة أضعاف. وفي كل حالة، كان الفرق في مدى انتشار القطرة بين حالتي الترطيب والطفو يمكن حسابه ببساطة من خلال النظر في مقدار الطاقة المفقودة بسبب الاحتكاك في الحالة الرطبة. ويسمح هذا الاكتشاف للعلماء بأخذ أي تنبؤ موجود لقطرة رطبة وتحويله فوراً إلى تنبؤ دقيق لقطرة عائمة، دون الحاجة إلى بناء نموذج جديد ومعقد من الصفر.
ويظل هذا الربط صحيحاً حتى عندما تكون القطرات لزجة جداً أو تتحرك بسرعة كبيرة، بشرط ألا تتفتت عند الاصطدام. وتؤكد الدراسة أن "قاعدة النصف"، وهي فكرة سابقة تقترح أن القطرات العائمة تحول دائماً نصف طاقتها إلى مساحة سطحية، ليست دقيقة إلا للقطرات الرقيقة والسريعة جداً. أما بالنسبة للقطرات الأكثر سمكاً أو الأبطأ، يصبح الاحتكاك الداخلي مؤثراً، ويأخذ هذا الربط الجديد في الاعتبار هذا الأمر من خلال إظهار مدى قلة انتشار القطرة مقارنة بالحد المثالي. يوفر هذا العمل طريقة موحدة لفهم اصطدام القطرات، جسر الفجوة بين الارتطام المألوف لسطح مبلل والانزلاق الغامض لقطرة ليدن فروست، مما يوفر مساراً أوضح للتحكم في هذه الاصطدامات في التطبيقات الصناعية والحرارية.
ملخص تقني: رسم خرائط شامل لانتشار اصطدام القطرات من أنظمة البلل إلى أنظمة ليدنفروست
بيان المشكلة يُعد اصطدام قطرة سائلة بسطح ما مسألة نموذجية في ديناميكا السوائل، تحكمها التفاعلات بين القصور الذاتي، والخاصية الشعرية، والتبدد اللزوجي. وبينما تم توصيف عامل الانتشار الأقصى (βm=Dm/D0) لعمليات الاصطدام في حالة البلل (حيث يحدث تلامس مباشر بين السائل والصلب) بشكل جيد عبر نطاقات واسعة من أرقام ويبر (We) وأوهنسجور (Oh)، فإن نظام ليدنفروست (LF) (حيث تطفو القطرة فوق طبقة بخار مستقرة) يفتقر إلى إطار تنبؤي عالمي مماثل. تعتمد نماذج ليدنفروست المستخدمة على نطاق واسع، مثل "قاعدة الـ 1/2" (وايلديمان وآخرون)، على افتراض تدفق غير لزج وتتجاهل التبدد اللزوجي في الكتلة، مما يؤدي إلى مبالغة منهجية في تقدير الانتشار عندما تكون Oh≳0.05. ولم يتم بعد وضع رسم خرائط شامل يمكنه تحويل تنبؤات البلل المستقرة إلى نظيراتها في نظام ليدنفروست.
المنهجية يضع المؤلفون رسم خرائط شاملاً يعتمد على منظور تبدد الطاقة، حيث يربط أقصى انتشار في نظام ليدنفروست بالانتشار في نظام البلل من خلال الفرق في إجمالي التبدد اللزوجي.
النهج التجريبي والعددي: تجمع الدراسة بين تجارب التصوير عالي السرعة والمحاكاة العددية المباشرة (DNS) باستخدام برنامج حل حجم السوائل مفتوح المصدر Basilisk.
التجارب: تم اصطدام قطرات مليمترية من الماء ومحاليل الجلسرين المائية بسطح من الياقوت الأزرق في ظروف غير مسخنة (بلل) وظروف مسخنة (550∘ مئوية، ليدنفروست).
المحاكاة: أُجريت محاكاة متماثلة محورياً لطورين؛ حيث استخدمت حالة البلل شرط حد عدم الانزلاق، بينما استخدمت حالة ليدنفروست شرط حد عدم وجود إجهاد قص لتمثيل كبح إجهاد القص بالقرب من الجدار بواسطة طبقة البخار.
تحليل ميزانية الطاقة: تقوم الدراسة بتفكيك إجمالي التبدد اللزوجي (Ed) إلى تبدد الطبقة الحدودية القريبة من الجدار (EBLd) وتبدد الكتلة (Ebulkd). ومن خلال مقارنة موازنات الطاقة لاصطدامات البلل وليدنفروست تحت ظروف أولية متطابقة، يعزل المؤلفون النقص في إجمالي التبدد اللزوجي (ΔEd=Ed,w−Ed,LF) المتاح للتوسع السطحي الإضافي في نظام ليدنفروست.
اشتقاق القياس: يشتق المؤلفون قانون قياس لفرق التبدد، حيث يوضح أن ΔEd يتناسب مع مساهمة الطبقة الحدودية للبلل (EBLd∼We−1/4Oh1/2βm,w2(βm,w−1)1/2). يؤدي هذا إلى معادلة رسم خرائط مباشرة تربط βm,LF بـ βm,w.
المساهمات والنتائج الرئيسية
معادلة رسم الخرائط الشاملة: يشتق البحث معادلة رسم خرائط مباشرة (المعادلة 4) تتنبأ بعامل انتشار ليدنفروست الأقصى من أي قانون انتشار بلل محدد مسبقاً: βm,LF=βm,w[1+607We3/4Oh1/2(βm,w−1)1/2]1/2 يتيح هذا الإطار تحويل نماذج البلل الراسخة (مثل نماذج سانجاي ولوه، أو ليو وآخرون، أو أبوت وبون) إلى تنبؤات ليدنفروست دون الحاجة لبناء نموذج مستقل ومنفصل لليدنفروست.
التحقق من الصحة: تم التحقق من صحة رسم الخرائط مقابل البيانات التجريبية والمحاكاة التي تغطي نطاقات 100≤We≤103 و 10−3≤Oh≤0.2. تغطي النت نیز نطاقين إلى ثلاثة أوامر قدرية في We وأمرين قدريين في Oh، وتتضمن بيانات للماء، ومزيج الجلسرين والماء، والإيثانول، والسوائل المفلورة.
الرؤية الفيزيائية: توضح الدراسة أنه بينما تعمل طبقة البخار في نظام ليدنفروست على كبح التبدد اللزوجي بالقرب من الجدار بقوة، فإن تبدد الكتلة يزداد ويعوض هذا الفقد جزئياً. ومع ذلك، يظل صافي النقص في إجمالي التبدد مقارباً لمساهمة الطبقة الحدودية للبلل، وهو ما يحدد المقياس المميز لرسم الخرائط.
الأداء: ينجح رسم الخرائط في التقاط الانتشار المعزز في نظام ليدنفروست والانحراف التدريجي عن "قاعدة الـ 1/2" غير اللزجة مع زيادة اللزوجة. كما يظل قوياً حتى عند قيم We المنخفضة، حيث يتلاشى حد التصحيح طبيعياً عندما βm,w→1.
الأهمية والادعاءات يدعي المؤلفون أنهم أرسوا اتصالاً طاقياً موحداً بين اصطدام قطرات البلل واصطدام قطرات ليدنفروست. ومن خلال إثبات أن الفرق في إجمالي التبدد اللزوجي بين النظامين ينعكس مباشرة في الفرق في عوامل الانتشار، يوفر البحث طريقاً بسيطاً لتوسيع الإطار النظري الواسع لاصطدامات البلل ليشمل نظام ليدنفروست.
يُقدم العمل كأداة عملية للتنبؤ بالاصطدام والتحكم فيه في تطبيقات مثل الإدارة الحرارية، وحقن الوقود، والدفع الذاتي، وحصاد الطاقة. ويشير المؤلفون إلى أن الإطار صالح ضمن نطاق المعلمات المدروسة؛ وقد تظهر قيود عند قيم We القصوى (بسبب الرذاذ أو عدم استقرار طبقة البخار) أو قيم We المنخفضة جداً (حيث تشكل التشوهات الضعيفة تحدياً للتقريبات الهندسية)، ولكن أظهرت هذه التأثيرات تأثيراً ضئيلاً على الانتشار الأقصى ضمن النطاق المدروس. كما يوفر رسم الخرائط طريقة لتقدير فرق إجمالي التبدد من خلال قياسات مزدوجة للبلل وليدنفروست، وهي كمية يصعب الوصول إليها مباشرة بطرق أخرى.