Equivalent Roughness Height for CFD Studies on NACA Airfoils
تقترح هذه الورقة إجراءً لتقدير ارتفاع خشونة حبيبات الرمل المكافئة من خلال تغطية جزيئات السطح للأجنح الهوائية من طراز NACA، وتقيم دقة نماذج الاضطراب المختلفة المزودة بتصحيحات الخشونة في التنبؤ بالأداء الديناميكي الهوائي تحت تدرجات ضغط معقدة، مما يعالج أوجه القصور في المعايرات الحالية القائمة على التدفقات النموذجية.
المؤلفون الأصليون:Alexis De Urraza, Juan M. Torres Zanardi, Ana Scarabino, Federico Bacchi
تخيل السماء كمحيط عملاق غير مرئي. تماماً كما ينزلق سمك عبر الماء، تشق الطائرات وتوربينات الرياح طريقها عبر هذا الهواء، معتمدة على طبقات انسيابية غير مرئية تلتصق بأجنحتها لتبقى فعالة. ولكن ماذا يحدث عندما لا تكون تلك الأجنحة ملساء تماماً؟ فكر في طبقة طلاء جديدة على سيارة؛ إذا رششت القليل من الرمل أو الغبار عليها بينما لا تزال رطبة، يصبح السطح خشناً. في عالم الديناميكا الهوائية، هذه "الخشونة" هي لص صامت؛ فهي تعطل التدفق السلس للهواء، مما يخلق سحباً إضافياً (مثل محاولة الجري عبر مياه يصل عمقها إلى الخصر بدلاً من الجري على مضمار سباق) وتسرق الطاقة من توربينات الرياح. لقد عرف المهندسون منذ زمن طويل أن الشفرات المتسخة أو المتآكلة تفقد كفاءتها، وأحياناً بنسبة تصل إلى النصف، لكن تحديد كيفية محاكاة هذا الاتساخ بدقة في جهاز الكمبيوتر يشبه محاولة وصف طعم قطعة بسكويت معينة باستخدام وصفة دقيق فقط. أنت بحاجة إلى طريقة لترجمة الحبيبات الفوضوية من العالم الحقيقي إلى أرقام رياضية نظيفة يمكن للكمبيوتر فهمها.
هنا يأتي دور فريق من الباحثين من جامعة لابلاتا الوطنية (Universidad Nacional de La Plata) بتقديم وصفة جديدة للمطبخ الرقمي. لقد عالجوا مشكلة معقدة: كيف يتم تحويل الوصف الغامض لـ "الخشونة القياسية" الموجودة في تجارب الأنفاق الهوائية القديمة إلى رقم دقيق يمكن لعمليات المحاكاة الحاسوبية الحديثة استخدامه. في الماضي، كان العلماء يختبرون الأشكال الانسيابية (المقطع العرضي للجناح) المغطاة بنوع معين من المواد التي تشبه ورق الصنفرة الخشن، لكن النماذج الحاسوبية لم تكن تعرف كيفية التعامل مع تلك الخشونة المحددة. كانوا بحاجة إلى "مترجم". طور المؤلفون طريقة ذكية لربط النسبة المئوية لتغطية سطح الجناح بهذه الجزيئات الصغيرة بـ "ارتفاع خشونة مكافئ" — وهو رقم واحد يخبر الكمبيوتر مدى وعورة السطح. ومن خلال إجراء آلاف الاختبارات الافتراضية في نفق هوائي على ثلاثة أشكال مختلفة من الأجنحة (NACA 0012، وNACA 2418، وNACA 633-418)، وجدوا أن كمية السطح المغطاة بالحصى المجهري مهمة للغاية. واكتشفوا أنه بالنسبة لبعض الأجنحة، كانت نسبة تغطية حوالي 10% من السطح هي الأفضل لمطابقة البيانات الواقعية، بينما كانت لغيرها نسبة 5% أو 7.5% هي النقطة المثالية. أظهرت عمليات المحاكاة الخاصة بهم أن طريقة الترجمة الجديدة هذه تعمل بشكل جيد بشكل مدهش، حيث تبقي تنبؤات الكمبيوتر ضمن نطاق 5% تقريباً من النتائج التجريبية الحقيقية لمعظم زوايا الطيران، رغم أن الكمبيوتر لا يزال يعاني قليلاً عندما يكون الجناح على وشك الانهيار (فقدان قوة الرفع). جوهرياً، لقد منحوا المهندسين مسطرة جديدة وأكثر دقة لقياس "الأوساخ الرقمية"، مما يساعدنا في بناء توربينات رياح وطائرات أفضل يمكنها التعامل مع واقع السماء الفوضوي.
ملخص تقني: ارتفاع الخشونة المكافئة لدراسات ديناميكا السوائل الحسابية (CFD) على الأجنحة الهوائية من طراز NACA
بيان المشكلة يعد التنبؤ الدقيق بتدهور الأداء الديناميكي الهوائي الناتج عن خشونة السطح أمراً بالغ الأهمية لريش توربينات الرياح والأنظمة الأخرى القائمة على الأجنحة الهوائية. تؤدي العوامل البيئية مثل الغبار، والجليد، وتأثير الحشرات، والتآكل إلى تغيير الحافة الأمامية، مما يؤدي إلى تعطيل الطبقة تحت اللزوجة (viscous sublayer) وزيادة معاملات احتكاك السطح. وبينما توفر البيانات التجريبية، لا سيما من كتاب "نظرية مقاطع الأجنحة" لأبوت ودون هوف (Abbott and Doenhoff)، معياراً لـ "الخشونة القياسية" (المعرفة بأحجام حبيبات الكاربوراندوم ونسب تغطية محددة)، فإن ترجمة هذه الظروف الفيزيائية إلى محاكاة ديناميكا السوائل الحسابية (CFD) يمثل تحدياً. تتطلب نماذج الاضطراب في الـ CFD عادةً ارتفاع خشونة مكافئ لحبيبات الرمل (ks)، بينما تُعرف الخشونة التجريبية غالباً بخصائص هندسية أو نسبة تغطية السطح. ولا يوجد تحويل مباشر بين حجم الحبيبات الفيزيائي المستخدم في اختبارات نفق الرياح ومعامل الخشونة الرقمي المنتظم المطلوب من قبل برامج الحل مثل ANSYS Fluent.
المنهجية تقترح الدراسة إجراءً لتقدير ارتفاع الخشونة المكافئ (ks) بناءً على النسبة المئوية لمساحة السطح المغطاة بالجسيمات. تسير المنهجية كما يلي:
توليف البيانات: قام المؤلفون بسد الفجوة بين التعريفات التجريبية لأبوت ودون هوف والعمل النظري لـ "شليختينج" (Schlichting) و"كولمان" (Coleman). قدم شليختينج جداول تربط بين هندسة الجسيمات وتوزيعها وبين الخشونة المكافئة، بينما قام كولمان وآخرون لاحقاً بتصحيح هذه القيم، موضحين أن التقديرات غير المصححة قد تصل إلى ضعف قيمة الخشونة المكافئة الفعلية.
تطوير الارتباط: من خلال تحليل توزيعات الجسيمات في بيانات شليختينج المصححة (تحديداً الجسيمات الكروية)، قام المؤلفون بحساب النسبة المئوية لمساحة السطح المغطاة لكل حالة. وتم إنشاء علاقة خطية بين النسبة المئوية لمساحة السطح المغطاة ونسبة الارتفاع المكافئ للخشونة إلى الخشونة الهندسية (ks,eq/K). والارتباط الناتج هو: ks,eq/K=19.19970⋅Covered Surf. [%]−0.12448
إعداد الـ CFD: أُجريت عمليات المحاكاة باستخدام نموذج الاضطراب k−ω SST مع تصحيحات الخشونة (Hellsten and Laine) لثلاثة أجنحة هوائية من طراز NACA: وهي NACA 0012، وNACA 2418، وNACA 633-418. حاكت عمليات المحاكاة ظروف اختبارات نفق الرياح لأبوت ودون هوف: رقم رينولدز يقارب 6×106 وكثافة اضطراب منخفضة (0.1–0.2%).
المعايرة: أُجريت عمليات محاكاة أولية لثلاث سيناريوهات تغطية (5%، 7.5%، و10%) لتحديد أي قيمة للخشونة المكافئة تحقق أفضل مطابقة لمنحنيات الرفع والسحب التجريبية لكل هندسة جناح محددة.
النتائج الرئيسية
الاعتماد على التغطية: وجدت الدراسة أن الخشونة المكافئة المطلوبة لمطابقة البيانات التجريبية تختلف باختلاف هندسة الجناح. فبالنسبة لـ NACA 0012، وفرت تغطية مساحة بنسبة 10% أفضل مطابقة؛ ولـ NACA 2418، كانت نسبة 5% هي المثلى؛ ولـ NACA 633-418، كانت نسبة 7.5% هي الأفضل مطابقة. يشير هذا إلى أن الخشونة الفعالة ليست مجرد دالة للتلوث الاسمي، بل تعتمد أيضاً على الخصائص الديناميكية الهوائية للجناح.
الدقة: ضمن نطاق زاوية الهجوم المنخفضة إلى المتوسطة (حوالي -6° إلى 5°)، أعادت عمليات محاكاة الـ CFD باستخدام قيم ks المستمدة إنتاج معاملات السحب التجريبية بفروقات تقل عموماً عن 5% لـ NACA 2418 وNACA 633-418. وأظهرت NACA 0012 انحرافات أعلى قليلاً، وصلت إلى 10% في بعض الظروف.
نظام التدفق: أكد التحليل الأولي أن جميع التكوينات المدروسة تعمل في نظام الخشونة الكاملة (ks+>140)، مما يؤكد صحة استخدام نماذج خشونة رمل الحبيبات المكافئة حيث يتم تعطيل تفاعلات الطبقة تحت اللزوجة.
محددات الانهيار (Stall): كما هو متوقع مع محاكاة RANS المستقرة، زادت الفوارق بين النتائج العددية والتجريبية بشكل كبير بالقرب من منطقة الانهيار. ويعزو المؤلفون ذلك إلى القيود المتأصلة في نموذج k−ω SST في التقاط الانفصال غير المستقر واسع النطاق، بدلاً من وجود قصور في إجراء تقدير الخشونة.
الأهمية والادعاءات يدعي البحث توفير إطار عمل عملي لترجمة التلوث السطحي المحدد تجريبياً (المتميز بتغطية الجسيمات) إلى معاملات خشونة متوافقة مع الـ CFD. ومن خلال إنشاء ارتباط خطي بين نسبة مئوية لتغطية السطح وارتفاع الخشونة المكافئ، تسمح المنهجية للباحثين بتجاوز الحاجة إلى الأوصاف الهندسية التفصيلية للتلوث عند نمذجة التدهور الديناميكي الهوائي.
يؤكد المؤلفون أنه بينما تم اشتقاق الارتباط من مجموعة محدودة من نقاط البيانات التجريبية ويفترض توزيعاً منتظماً للجسيمات الكروية، إلا أنه يوفر نهجاً قابلاً للتطبيق للتنبؤات على المستوى الهندسي. وتؤكد الدراسة أن الإجراء المقترح يمكنه إعادة إنتاج منحنيات السحب بدقة معقولة لظروف الخشونة القياسية، مما يسمح بمحاكاة عددية أكثر واقعية لفقدان الأداء الناتج عن تراكم الجسيمات أو تآكل السطح في تطبيقات طاقة الرياح والفضاء.
المحددات يشير المؤلفون صراحةً إلى أن الارتباط يعتمد على عدد محدود من نقاط البيانات ويفترض توزيعاً منتظماً للجسيمات الكروية، والتي قد لا تعكس الأشكال غير المنتظمة والتجمعات الموجودة في التلوث الواقعي. علاوة على ذلك، فإن الاعتماد على محاكاة RANS المستقرة يحد من دقة التنبؤ في منطقة الانه