محاكاة القاذف بطريقة CFD تعتمد كثيراً على نموذج الاضطراب. ما تذكره دراسة تحقق منشورة عن k-ω SST، وأين يضعف النموذج، وكيف تقارن بين التصاميم.

رسم توضيحي من DEEB-FLOW لحقل الضغط في قاذف
رسم تصوّري من مكتبة حالات DEEB-FLOW. وهو ليس نتيجة من هذا المقال.

القاذف يبدو بسيطاً ومحاكاته صعبة

ليس في القاذف (الإيجكتور) أجزاء متحركة، وأبعاده المؤثرة قليلة. يُفترض أن تكون محاكاته سهلة. ثم تشغّل الحالة نفسها بنموذجَي اضطراب، فتتغير نسبة السحب بأكثر من التحسين الذي كنت تأمل تحقيقه بالتصميم.

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

ما الذي يُراد التنبؤ به

الرقم الرئيسي هو نسبة السحب، أي التدفق الكتلي الثانوي مقسوماً على التدفق الكتلي الأولي. وللتوضيح، خذ فوهة أولية يمر فيها 0.050 كغ/ث وقاذفاً يسحب 0.021 كغ/ث من التدفق الثانوي.

ω = m_s / m_p = 0.021 / 0.050 = 0.42

الرقم الثاني هو الضغط الخلفي الحرج، أي أعلى ضغط مخرج يظل القاذف عنده يسحب بكامل طاقته.

ثلاثة أنظمة تشغيل ولماذا تهم المحاكاة

يتصرف القاذف تصرفاً مختلفاً بحسب ضغط مخرجه.

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

أي استنتاج من محاكاة قاذف ينبغي أن يذكر النظام الذي ينطبق عليه. النموذج الجيد في الأول قد يكون عديم الفائدة في الثاني.

ما تذكره دراسة تحقق

قارن بيزانيي وزملاؤه في جامعة بوليتكنيكو دي ميلانو عدة مقاربات لنمذجة الاضطراب في القواذف فوق الصوتية بالقياسات. أعطى نموذج k-ω SST أفضل توافق. بلغ خطؤه في نسبة السحب 14% في المتوسط و20% حداً أقصى.

هذه الأرقام تخص قواذف تلك الدراسة وظروفها، فاقرأها على أنها مرتبة مقدار. يمكن أن يُتوقع من نموذج مبني بعناية مع k-ω SST انحراف يقارب 15% عن نسبة السحب المقيسة في نظام التصميم. أما في النظام دون الحرج، حيث النتيجة أشد حساسية للمزج المتوقع، فتوقّع أقل من ذلك.

لماذا يهم النموذج هنا إلى هذا الحد

ثلاثة مواضع في الجريان تحسم الجواب، وكلها صعبة على نموذج ذي معادلتين.

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

طُوّر نموذج k-ω SST لتدرجات الضغط المعاكسة والانفصال، ولذلك يحسن الأداء هنا في الغالب. وتبقى الصعوبتان الأخريان.

كيف أبني نموذجاً كهذا

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

استخلاص استنتاج من رقم غير مؤكد

افترض أن الحالة الأساسية تعطي نسبة سحب 0.42 وأن إرجاع مخرج الفوهة 4 مم عكس اتجاه الجريان يعطي 0.45. مع ارتياب 15% قد تقع الحالة الأساسية وحدها في أي مكان بين 0.36 و0.48. وإذا أُخذت القيمتان مطلقتين مستقلتين كان التحسين داخل الضجيج.

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

يصح ذلك ما دامت الحالتان في النظام نفسه. إذا اقتربت الهندسة المعدلة من النقطة الحرجة دخلت المقارنة منطقة يضعف فيها النموذج، واحتاج الاستنتاج إلى اختبار.

عملياً، تحقق من نقطة واحدة بالقياس إن استطعت، ثم استخدم النموذج للاتجاهات حولها.

أخطاء شائعة

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

ما الذي تصلح له المحاكاة هنا

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

أما نسبة سحب مضمونة حتى منزلتين عشريتين فخارج قدرتها. القاذف الذي يُباع بضمان أداء ينبغي تأكيد ضمانه على منصة اختبار.

المصادر والمطالعة

خدمة ذات صلة

ديناميكا الموائع الحسابية والجريان

مقارنة الهندسة وظروف التشغيل قبل اعتماد التصميم.

استعرض الخدمات
كل المقالات