A generalised Spalart-Allmaras model with blended near-wall and free-shear coefficients: multi-objective calibration in an open-source framework
تقدم هذه الورقة نموذج سبالارت-ألماراس للاضطراب المعمم الذي يستخدم دالة حجب لدمج معاملات القرب من الجدار والقص الحر، مما يحقق دقة محسنة في التدفقات المنفصلة وحرة القص مع الحفاظ على المتانة في التدفقات الملتصقة من خلال معايرة متعددة الأهداف ضمن إطار عمل مفتوح المصدر.
المؤلفون الأصليون:Tony Di Fabbio, Riccardo Vallino, Eric Segalerba, Viola Rossano, Joel Guerrero, Andrea Ferrero
تحيط بنا السوائل من كل جانب، من الهواء المتدفق فوق جناح طائرة إلى الماء الذي يدور داخل أنبوب. إن التنبؤ بكيفية حركة هذه السوائل أمر ضروري لتصميم مركبات أكثر كفاءة، ومبانٍ أكثر أماناً، وأنظمة طاقة أنظف. وللقيام بذلك، يعتمد المهندسون على محاكاة حاسوبية تحل معادلات معقدة تصف حركة كل جسيم. ومع ذلك، فإن هذه المعادلات صعبة للغاية بحيث لا يمكن حلها بدقة للمسائل الواقعية، لذا يستخدم العلماء نماذج مبسطة. ومن أكثر الأدوات موثوقية لعقود من الزمن نموذج سمي باسم مبتكريه، "سبالارت وألماراس" (Spalart and Allmaras). يعمل هذا النموذج بشكل جيد للغاية في التدفقات السلسة والمتصلة، مثل انزلاق الهواء على طول جناح، لأنه يستخدم مجموعة واحدة وثابتة من القواعد. ولكن عندما يصبح التدفق فوضوياً — كالانفصال عن سطح، أو الدوران في زاوية، أو الاختلاط بحرية في نفاث — فإن تلك القواعد الثابتة نفسها غالباً ما تفشل. وكان التحدي يكمن في جعل النموذج أكثر ذكاءً دون كسر البساطة التي تجعله مفيداً للغاية.
لقد اتخذ فريق من الباحثين من إيطاليا خطوة كبيرة نحو حل هذه المعضلة؛ حيث ابتكروا نسخة عامة من النموذج الكلاسيكي يمكنها تكييف سلوكها اعتماداً على مكان وجود السائل وما يفعله. فبدلاً من إجبار المحاكاة بأكملها على اتباع مجموعة واحدة جامدة من التعليمات، يسمح نموذجهم الجديد للقواعد بالتغير محلياً. ففي المناطق التي يلتصق فيها السائل بجدار، يستخدم مجموعة واحدة من المعايير؛ وفي المناطق التي ينفصل فيها السائل أو يختلط بحرية، ينتقل إلى مجموعة أخرى. ويتم دمج هاتين المجموعتين من القواعد بسلاسة بواسطة مستشعر رقمي يكتشف ظروف التدفق المحلية. ولم يخمن الباحثون هذه القواعد الجديدة، بل استخدموا عملية تحسين حاسوبية قوية لإيجاد أفضل القيم الممكنة لها، واختبروا النموذج مقابل مجموعة واسعة من التدفقات المرجعية، من الصفائح المسطحة إلى النفاثات المعقدة والمناطق المنفصلة.
تظهر النتائج أن هذا النهج نجح في تحسين دقة النموذج في المناطق الصعبة دون إفساد أدائه في المناطق السهلة. فعلى سبيل المثال، في الاختبارات التي تضمنت نفاثاً هوائياً حراً، قلل النموذج الجديد من خطأ التنبؤ بمقدار عشرين ضعفاً مقارنة بالنموذج الأصلي. وفي الاختبارات التي تضمنت نتوءاً على جدار يتسبب في انفصال التدفق، انخفض الخطأ بمقدار الخُمس. والأهم من ذلك، حافظ النموذج على دقته العالية في التدفقات السلسة والمتصلة حيث كان الإصدار الأصلي ممتازاً بالفعل. وهذا يشير إلى أن قصور النموذج القديم لم يكن ناتجاً عن خلل جوهري في بنيته، بل عن حقيقة أنه حاول استخدام نفس الإعدادات الثابتة لكل نوع من أنواع التدفق. ومن خلال السماح لهذه الإعدادات بالتغير، أطلق الباحثون العنان لمستوى أعلى بكوثير من الدقة.
ومع ذلك، كشفت الدراسة أيضاً عن حدود ما يمكن إصلاحه بمجرد تغيير الأرقام. ففي اختبار محدد يتعلق بقناة مربعة، أحدث التدفق حركة ثانوية في الزوايا لم يستطع النموذج القياسي إعادة إنتاجها، مهما تم تعديل المعاملات. ويعود ذلك إلى أن الفيزياء التي تقود تلك الحركة تتطلب علاقة أكثر تعقيداً بين إجهاد السائل وشكله، وهو ما لا يمكن لنموذج خطي بسيط توفيره. ولحل هذه المشكلة، أضاف الباحثون تصحيحاً غير خطي محدداً للنموذج. وعندما فعلوا ذلك، نجح النموذج في توليد الدوامات الزاوية الدوارة التي تُلاحظ في البيانات عالية الدقة، مما قلل الخطأ في تلك الحالة المحددة بنسبة 85 بالمائة. وتوضح هذه النتيجة الفرق بين المشكلات التي يمكن حلها عبر ضبط المعاملات وتلك التي تتطلب تغييراً أعمق في بنية النموذج.
كما اختبر الفريق مدى جودة أداء نماذجهم الجديدة في تدفقات لم يروها من قبل، مثل حالة فقدان جناح الطائرة لتوازنه (الانهيار) أو تفاعل موجة صدمة مع سطح ما. وقد أدى النموذج الذي تمت معايرته على التدفقات الخارجية أداءً جيداً جداً في هذه الحالات غير المرئية، محافظاً على دقة النموذج الأصلي مع تقديم تحسينات في مجالات محددة. ويشير هذا إلى أن النهج الجديد ليس مجرد طريقة محدودة بحالات الاختبار المستخدمة أثناء التدريب، بل هو تحسين قوي يمكن تطبيقه على التحديات الهندسية الواقعية. وقد أتاح الباحثون كود البرمجة والبيانات الخاصة بهم كبرمجيات مفتوحة المصدر، مما يسمح لعلماء آخرين بالبناء على هذا العمل. إن نجاحهم يثبت أن نمذجة الاضطراب يمكن أن تتطور من خلال أن تصبح أكثر مرونة وتكيفاً، لتردم الفجوة بين الأدوات البسيطة والمتينة وبين الواقع المعقد لحركة السوائل.
ملخص تقني: نموذج سبالارت-ألماراس المعمم بمعاملات مدمجة للاقتراب من الجدار والمجالات ذات القص الحر
بيان المشكلة يُعد نموذج سبالارت-ألماراس (SA) معياراً صناعياً لمحاكاة معادلات نافيير-ستوكس بمتوسط رينولدز (RANS) نظراً لمتانته وانخفاض تكلفته الحسابية. ومع ذلك، تعتمد دقته على مجموعة واحدة من المعاملات الثابتة عالمياً والتي تمت معايرتها بناءً على تدفقات مقيدة بالجدار في حالة التوازن. هذا التوحيد يحد من أدائه في التدفقات المنفصلة، ومجالات القص الحر، والهندسات المعقدة حيث تختلف فيزياء الاضطراب بشكل كبير عن خط الأساس الذي تمت المعايرة عليه. وبينما بذلت الجهود السابقة لإضافة تصحيحات مستهدفة (مثل تصحيحات الدوران/الانحناء) أو توسيع نطاق المعايرة لآليات محددة، إلا أنها غالباً ما تحتفظ بمعاملات عالمية ثابتة أو تتطلب تبديلاً مناطقياً. السؤال الجوهي الذي يعالجه البحث هو: إلى أي مدى يمكن تعميم نموذج SA ليتكيف مع الأنظمة التدفقية المحلية دون فقدان متانته المتأصلة أو الحاجة إلى معايير تبديل معقدة يحددها المستخدم؟
المنهجية يقترح المؤلفون نموذج سبالارت-ألماراس المعمم (GSA) حيث تمتلك معاملات معايرة مختارة (σν~، Cb1، Cb2) حالتين متمايزتين: حالة داخلية (Inner) للاضطراب القريب من الجدار والمتصل، وحالة خارجية (Outer) للمناطق المنفصلة وذات القص الحر. يتم دمج هاتين الحالتين باستمرار في الفراغ باستخدام دالة حجب (FSA) مستمدة من مفاهيم محاكاة الدوامات المتأخرة المنفصلة (DDES) والـ IDDES المحسنة. يضمن هذا الدمج تكيف النموذج مع نظام التدفق المحلي دون الحاجة إلى تقسيم مناطقِي.
تشمل المكونات المنهجية الرئيسية ما يلي:
دمج المعاملات: يتم حساب كل معامل ψ عبر المعادلة ψ(x)=FSA(ψI−ψO)+ψO. كما تمت إعادة بناء معامل تدمير الجدار Cw1 محلياً للحفاظ على علاقة اتساق الطبقة اللوغاريتمية في كلتا الحالتين الحدّيتين.
دالة الحجب: دالة هجينة تجمع بين فرع المسافة عن الجدار (لحماية الطبقات الحدودية المتصلة) وفرع رقم رينولدز الدوامي (لاكتشاف الانفصال والقص الحر). تم إدخال امتداد حساس للطوبولوجيا لرفع "أرضية" مستشعر الدوامات في المناطق التي يهيمن عليها الدوران، مما يعالج التعارض بين التنبؤ بالانفصال وتوليد التدفق الثانوي.
الامتداد غير الخطي: لمعالجة عدم قدرة نماذج اللزوجة الدوامية الخطية على إعادة إنتاج التدفقات الثانوية الناتمة عن تباين الخواص (مثل الممرات المربعة)، يدمج النموذج علاقة التكوين التربيعية (QCR2024) مع معاملات معايرة.
المعايرة متعددة الأهداف: تمت معايرة المعاملات باستخدام إطار عمل مفتوح المصدر يربط بين OpenFOAM-13 والمحسنات التطورية (PSO لتوليد البذور لحالة واحدة، وNSGA-II لمجموعات الحالات المتعددة). قللت عملية التحسين من متجه التكلفة حيث تم تطبيع كل حد خطأ بواسطة خطأ نموذج SA الأساسي، مما يضمن أن التحسينات في فئة تدفق معينة لا تخفي تدهور الأداء في فئة أخرى.
مجموعات الاختبار المرجعية: تمت المعايرة على مجموعتين متميزتين: مجموعة NASA لنمذجة الاضطراب (TMR) (لوح مسطح، قناة، نفث، وتلة) ومجموعة ERCOFTAC (تلال شديدة/لطيفة، ممر مربع، ودرجة خلفية منحنية).
المساهمات والنتائج الرئيسية تفصل الدراسة بشكل منهجي بين القيود الناجمة عن اختيار المعاملات وتلك المتأصلة في بنية الإغلاق.
التعميم دون تنازلات: يحافظ نموذج GSA على دقة نموذج SA الأساسي للتدفقات المقيدة بالجدار المتصلة (اللوح المسطح، القناة) مع تحسين الأداء بشكل كبير في الأنظمة الصعبة.
في مجموعة NASA، خفض متغير GSA-N هدف "النفث" بمقدار عشرين ضعفاً وهدف "التلة المثبتة على الجدار" بمقدار خمس واحد مقارنة بأساس SA.
في مجموعة ERCOFTAC، حسن متغير GSA-E التنبؤات للتلال اللطيفة والدرجات المنحنية.
صراع المعاملات والقابلية للنقل: تسلط الدراسة الضوء على أن فئات التدفق المختلفة تتطلب حالات خارجية (Outer) متضاربة. فعلى سبيل المثال، يتطلب النفث الحر إنتاجاً مخفضاً (Cb1O)، بينما يتطلب الانفصال الناتج عن الضغط على تلة زيادة في الإنتاج. وبناءً على ذلك، فإن النموذج الذي تمت معايرته على مجموعة NASA (GSA-N) ينتقل جيداً إلى حالات ERCOFTAC، ولكنه يفشل في حالة نفث NASA إذا تمت معايرته بناءً على ERCOFCT (GSA-E)، حيث ينمو الإنتاج الخارجي بشكل كبير مما يؤدي إلى انهيار النواة المحتملة للنفث.
القيود الهيكلية وQCR: تثبت الورقة أنه لا توجد مجموعة من المعاملات يمكنها توليد التدفق الثانوي في ممر مربع، لأن هذا ناتج عن تأثير تباين إجهاد رينولدز الذي يتجاوز الإغلاقات الخطية.
ساعدت أرضية الحجب الحساسة للطوبولوجيا جزئياً في تخفيف المقايضة بين التنبؤ بالانفصال والتدفق الثانوي.
نجح امتداد QCR2024 المعاير (GSA-E+QCR2024) في توليد الدوامات الزاوية في الممر المربع، مما قلل هدف الممر بنسبة 85% مقارنة بنموذج الدمج فقط، مع تحسين التنبؤات بالتلال أيضاً.
التحقق من حالات غير مرئية: تم اختبار النماذج على ثلاث تكوينات لم تُستخدم في المعايرة: قطب الجناح NACA 0012، وتلة عابرة للصوت متماثلة المحور (ATB)، ووصلة جناح-جسم (WBJ).
حافظ GSA-N على خصائص الانهر (Stall) لنموذج SA الأساسي على الجناح، بينما أخر GSA-E الانهر بسبب ارتفاع الإنتاج الخارري (وهو عجز تم تصحيحه بواسطة QCR2024).
في حالة ATB، لم تؤدِ نماذج الدمج إلى تدهور في فيزياء التفاعل بين الصدمة والطبقة الحدودية.
في حالة WBJ، قلل GSA-N من أخطاء السرعة بنسبة 15-18% مقارنة بأساس SA، رغم أن طوبولوجيا الدوامة ظلت دون تغيير، مما يشير إلى أن الخطأ المتبقي يعود إلى فرضية RANS المستقرة وليس إلى معايرة المعاملات.
الواقعية: أكد تحليل ثبات الخواص (Anisotropy-invariant) أن إجهاد رينولدز النموذجي يظل واقعياً (شبه محدد موجباً) في الممر المربع وفي مناطق القص/الدوران في الوصلة، حتى مع وجود حدود QCR غير الخطية.
الأهمية والادعاءات تدعي الورقة أن التعميم المقترح يوسع نطاق تطبيق نموذج SA دون التضحية بمتانته أو قابليته للنقل. تكمن الأهمية الأساسية في الفصل المنهجي بين أوجه القصور في النموذج:
المعامل مقابل البنية: يحدد العمل أي الأخطاء يمكن علاجها بمعايرة المعاملات (مثل انتشار النفث، طول الانفصال) وأيها هيكلي (مثل التدفقات الثانوية في الممرات)، مما يستلزم علاقات تكوينية غير خطية.
المعايرة المتكيفة مع النظام: يثبت العمل أن صياغة نموذج واحدة يمكنها التكيف مع أنظمة تدفق مختلفة عبر دمج المعاملات محلياً، مما يتجنب الحاجة إلى التبديل المناطقي.
القابلية للتفسير القائمة على البيانات: من خلال استخدام التحسين متعدد الأهداف على مجموعات الاختبار، تحدد الدراسة حالات المعاملات المحددة التي تقود المكاسب والمقايضات، مما يقدم بديلاً قابلاً للتفسير للتحسينات "الصندوق الأسود" القائمة على البيانات.
يخلص المؤلفون إلى أنه بينما يحسن نموذج GSA الدقة عبر نطاق واسع من التدفقات، فإن "أفضل" نموذج يعتمد على التطبيق المستهدف. يُفضل GSA-N للديناميكا الهوائية الخارجية حيث يكون التنبؤ بالانهر أمراً حاسماً، بينما يعد GSA-E+QCR2024 متفوقاً للتدفقات الداخلية ذات الحركة الثانوية القوية. يشير العمل إلى أن نمذجة الاضطراب المستقبلية يجب أن تركز على نماذج ذات أسس فيزيائية وقابلة للتفسير، حيث تتكيف المعاملات مع طوبولوجيا التدفق، بدلاً من مجرد التصحيحات القائمة على البيانات. تم إصدار جميع البرمجيات والبيانات وأدوات المعايرة كمصادر مفتوحة لتسهيل التطوير المستقبلي.