Variational Leading-Edge Profiles for Centrifugal-Compressor Impellers under Prescribed Inlet Swirl

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

لفهم جوهر هذه المشكلة، يجب على المرء أولاً أن يستوعب أن الهواء الداخل إلى الضاغط لا يتحرك في خط مستقيم؛ بل غالباً ما يصل بدوران، وهو ما يُعرف بـ "الدوامة" (swirl)، وتوزيع محدد للسرعة عبر نصف قطر العجلة. والهدف من الشفرة هو توجيه هذا الهواء بسلاسة فوق سطحها. وفي العالم الحقيقي، يؤدي تغيير شكل الشفرة إلى تغيير التدفق، وتغيير التدفق يغير القوى المؤثرة على الشفرة، مما يخلق شبكة متشابكة من الأسباب والنتائج. ولتفكيك هذا التشابك، قرر باحث يدعى ف. ف. غوتسولينكو من معهد الفيزياء الحرارية الهندسية في كييف، أوكرانيا، النظر في نسخة مبسطة من المشكلة. فبدلاً من محاولة حل التفاعل الكامل والمعقد بين الهواء والمعدن دفعة واحدة، قام بعزل خط واحد أحادي البعد: وهو المنحنى الذي يحدد حافة الشفرة أثناء حركتها من المركز الداخلي إلى الحافة الخارجية. وقد تعامل مع ظروف دخول الهواء كحالة ثابتة، مثل لقطة لحالة تشغيل محددة، وطرح سؤالاً بسيطاً: بالنظر إلى هذه اللقطة المحددة للهواء الدوار، ما هو المنحنى المثالي رياضياً لحافة الشفرة لتقليل مقياس معين من "الجهد"؟

تركز الدراسة على نموذج رياضي يتم فيه تعريف "الجهد" بكيفية تفاعل الشفرة مع الزخم الزاوي للهواء الداخل. في هذا النموذج، يتم تجميد سرعة الهواء ودورانه في مكانها، حيث يتم أخذهما من حالة مرجعية، ويقوم الباحث بتغيير منحنى الشفرة فقط. الهدف هو إيجاد المنحنى الذي يقلل قيمة معينة مشتقة من طول الشفرة وموقعها بالنسبة لدوران الهواء. لا يدعي هذا النهج إيجاد أفضل شفرة لكل حالة ممكنة للمحرك، كما أنه لا يأخذ في الاعسبار الفيزياء المعقدة للاحتكاك أو الحرارة أو سمك المعدن. بدلاً من ذلك، يوفر معياراً تحليلياً صافياً — حلاً مثالياً لعالم مبسط — يمكن أن يكون بمثابة نقطة انطلاق أو حالة اختبار لمحاكاة حاسوبية أكثر تعقيداً.

اكتشف الباحث أن هذه المشكلة الصعبة يمكن اختزالها إلى مهمة أبسط بكثير. فباستخدام مبدأ كلاسيكي من مبادئ التفاضل والتكامل، تم تحويل المعادلة المعقدة التي تصف المنحنى إلى عملية حسابية واحدة سهلة تتضمن معاملاً يتحكم في مدى حدة دوران المنحنى. أثبتت الدراسة أنه بالنسبة لمجموعة واسعة من ظروف تدفق الهواء، يوجد بالضبط منحنى فريد واحد يلبي المتطلبات. علاوة على ذلك، أرست الدراسة قاعدة واضحة لمتى يمكن لوجود مثل هذا المنحنى أصلاً: وهي أن المسافة التي يجب أن تقطعها الشفرة من المركز الداخلي إلى الحافة الخارجية لا يمكن أن تكون طويلة جداً مقارنة بالفرق بين نصف قطريهما. فإذا كانت الشفرة طويلة جداً بالنسبة للمساحة المتاحة، فلن يوجد منحنى سلس ومثالي يحقق الشروط. ويوفر هذا الاكتشاف للمهندسين اختباراً نهائياً بـ "نعم أو لا" قبل أن يبدؤوا حتى في التصميم.

كما استكشفت الدراسة سيناريوهين شائعين لتوضيح الشكل الذي تبدو عليه هذه المنحنيات المثالية. في الحالة الأولى، حيث يدور الهواء بطريقة شائعة في العديد من المحركات، يتخذ حافة الشفرة المثالية شكل "المنحنى السلس" (catenary)؛ وهو نفس المنحنى الذي يشكله سلسلة ثقيلة تتدلى بحرية بين نقطتين. وفي الحالة الثانية، حيث يدور الهواء بشكل منتظم، يكون المنحنى أكثر تعقيداً، ويُوصف بدالة رياضية خاصة تُعرف باسم "التكامل الناقص" (elliptic integral)، والتي يمكن تصورها كقوس صاعد سلس يختلف قليلاً عن السلسلة المتدلية. وقد حسب الباحث هذه الأشكال لمثال محدد، باستخدام شفرة تبدأ من نصف قطر 50 ملم وتنتهي عند 90 ملم بطول 30 ملم. وأظهرت النتائج أن المنحنى الأمثل ليس خطاً مستقيماً، بل ينحني للداخل قليلاً، حيث يبقى قريباً من المركز لفترة أطول قبل أن يرتفع بحدة ليلتقي بالحافة الخارجية. وهذا السلوك هو استجابة رياضية مباشرة لكيفية توزيع زخم الهواء.

إن ما يجعل هذا العمل قيماً للغاية ليس فقط شكل المنحنى، بل اليقين الذي استُخلص به. فقد أثبت الباحث أن الحل ليس مجرد حد أمثل محلي، بل هو حد أدنى صارم وقوي، مما يعني أن أي انحراف طفيف عن هذا المنحنى سيؤدي إلى نتيجة أسوأ. وقد تم إثبات ذلك من خلال إظهار أن عائلة المنحنيات الممكنة لا تتقاطع مع نفسها، وأن الشروط الرياضية للحل المثالي قد استوفيت. إن الطريقة المستخدمة لإيجاد هذه المنحنيات قوية وفعالة، حيث تعتمد على بحث عددي مباشر يضمن العثور على الإجابة إذا كانت موجودة. وهذا يوفر أداة موثوقة للمهندسين لإنشاء تصميمات أولية بسرعة، والتي يمكن بعد ذلك صقلها باستخدام عمليات محاكاة أكثر تفصيلاً واستهلاكاً للموارد الحاسوبية.

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