۵ نکته کلیدی و غافلگیرکننده در شبیه‌سازی جریان تراکم‌پذیر بر روی ایرفویل NACA 0012

مقدمه شبیه‌سازی پرواز در رژیم جریان انتقالی یا ترانسونیک (Transonic)، همواره یکی از جذاب‌ترین و در عین حال چالش‌برانگیزترین حوزه‌ها برای مهندسان آیرودینامیک بوده است. عبور از مرزهای سرعت‌های پایین و ورود به دنیای اثرات تراکم‌پذیری، جایی است که معادلات ساده سیالات دیگر پاسخگو نیستند. اما چقدر می‌توان به نتایج عددی نرم‌افزارهایی مانند ANSYS CFX در مقایسه با واقعیت‌های فیزیکی ثبت شده در تونل‌های باد ناسا (NASA) اعتماد کرد؟ در این مقاله، با نگاهی تحلیلی به شبیه‌سازی ایرفویل NACA 0012، به بررسی ۵ درسی می‌پردازیم که فراتر از آموزش‌های معمول، دقت شبیه‌سازی‌های شما را تضمین می‌کنند.

درس اول: خطای مقیاس؛ وقتی ۱ میلی‌متر به ۱۰۰۰ متر تبدیل می‌شود! یکی از رایج‌ترین و در عین حال مخرب‌ترین اشتباهات در شبیه‌سازی، نادیده گرفتن واحدهای اندازه‌گیری هنگام وارد کردن مش (Mesh) است. در پروژه‌های آیرودینامیکی، وتر ایرفویل (Chord Length) مبنای تمام محاسبات و اعداد بدون بعد است. اگر مش در محیطی با واحد میلی‌متر طراحی شده باشد و بدون تغییر در CFX وارد شود، نرم‌افزار آن را بر حسب متر تفسیر می‌کند. نتیجه؟ ایرفویل یک متری شما ناگهان به یک غول ۱۰۰۰ متری تبدیل می‌شود که تمام معادلات رینولدز را به هم می‌ریزد.

بررسی اولیه در CFX-Pre تنها یک پیشنهاد نیست، بلکه یک ضرورت حیاتی است. برای اصلاح این خطا، استفاده از ابزار Transform Mesh و اعمال ضریب مقیاس 0.001 (Scale) تنها راه بازگرداندن ابعاد به واقعیت فیزیکی است.

“پس از باز کردن CFX-Pre، مش را می‌توان بررسی کرد و مشخص است که مقیاس اشتباه است، زیرا وتر ایرفویل به جای ۱ متر، ۱۰۰۰ متر است؛ این نشان می‌دهد که مش به جای متر، با واحد میلی‌متر ساخته شده است.”

درس دوم: آیرودینامیک سرعت بالا؛ فیزیک پنهان در لایه‌مرزی در عدد ماخ ۰.۷، با جریانی روبرو هستیم که به شدت تحت تاثیر اثرات تراکم‌پذیری است. در این رژیم، مدل‌سازی دقیق لایه‌مرزی و میدان‌های حرارتی تعیین‌کننده است. برای این شبیه‌سازی، انتخاب مدل توربولانسی Shear Stress Transport (SST) یک انتخاب استراتژیک است؛ چرا که این مدل در مدیریت گرادیان‌های فشار و پیش‌بینی جدایش جریان، عملکردی به مراتب برتر از مدل‌های استاندارد k-epsilon دارد.

در تنظیمات Domain، استفاده از گزینه Total Energy ضروری است، اما نکته کلیدی‌تر فعال‌سازی Incl. Viscous Work Term است. در سرعت‌های بالا، اصطکاک در لایه‌مرزی گرمای قابل توجهی تولید می‌کند (Viscous Heating) که مستقیماً بر چگالی محلی و توزیع فشار اثر می‌گذارد. نادیده گرفتن این ترم، یعنی نادیده گرفتن بخشی از فیزیک واقعی جریان ترانسونیک. همچنین برای تطابق با داده‌های تجربی (عدد رینولدز 9×10^6 در دمای ۲۸۸ کلوین)، فشار مرجع (Reference Pressure) باید دقیقاً روی ۵۶۸۶۷ پاسکال تنظیم شود تا چگالی سیال در شرایط ایده‌آل حفظ شود.

درس سوم: ترفند زاویه حمله؛ بهره‌وری مهندسی با Expressions تغییر زاویه حمله (AOA) در مطالعات پارامتریک معمولاً به معنای چرخاندن هندسه و تولید مجدد مش است که فرآیندی زمان‌بر و مستعد خطاست. یک روش هوشمندانه برای اعمال زاویه حمله ۱.۴۹ درجه، استفاده از عبارات محاسباتی (Expressions) برای تجزیه سرعت ورودی است.

در این روش، به جای تغییر فیزیکی مش، جهت جریان ورودی را تغییر می‌دهیم. با تعریف سرعت جریان آزاد (Uinf = 238.12 m/s)، مولفه‌های سرعت را به صورت زیر تعریف می‌کنیم:

• Ux = Uinf * cos(AOA)
• Uy = Uinf * sin(AOA)

این “هک بهره‌وری” باعث می‌شود بدون نیاز به Re-meshing، تنها با تغییر یک عدد در بخش Expressions، زاویه حمله را تغییر دهید و در زمان ارزشمند پیش‌پردازش (Pre-processing) صرفه‌جویی کنید.

درس چهارم: فریب باقی‌مانده‌ها را نخورید؛ هنر همگرایی بسیاری از تحلیل‌گران تصور می‌کنند رسیدن به باقی‌مانده‌های (RMS Residuals) در سطح 1e-4 به معنای پایان کار است. اما در جریان‌های پیچیده، همگرایی ریاضی لزوماً به معنای پایداری فیزیکی نیست. در این پروژه، معیار 1e-4 کافی نبود؛ زیرا در این سطح، مقادیر ضریب برا (CL) هنوز نوسان داشتند.

نویسنده با سخت‌گیرانه‌تر کردن معیار همگرایی به 1e-6 و افزایش Timescale Factor از ۱ به ۱۰، فرآیند رسیدن به پاسخ پایدار را تسریع و تثبیت کرد. افزایش این فاکتور یک ترفند حرفه‌ای برای رسیدن سریع‌تر به حالت Steady-state در حل‌گرهای کوپل شده است. همگرایی واقعی تنها زمانی تایید می‌شود که Monitor Points مربوط به CL و CD در Solver Manager به یک خط کاملاً مستقیم تبدیل شوند.

درس پنجم: ظهور مناطق مافوق صوت در قلب جریان زیرصوت یکی از هیجان‌انگیزترین نتایج این شبیه‌سازی، مشاهده مناطق مافوق صوت (Supersonic) بر روی سطح ایرفویل است. با وجود اینکه جریان ورودی ۰.۷ ماخ است، کانتورهای عدد ماخ نشان می‌دهند که سرعت محلی در نزدیکی لبه حمله از سرعت صوت (۳۴۰ متر بر ثانیه) فراتر رفته و به حدود ۳۵۰ متر بر ثانیه می‌رسد.

دقت این شبیه‌سازی با دو معیار اصلی اثبات می‌شود: ۱. ضریب برا (CL): مقدار عددی ۰.۲۳۶ که در تطابق بسیار نزدیک با داده تجربی ناسا (۰.۲۴۱) قرار دارد. ۲. توزیع ضریب فشار (Cp): نمودار توزیع فشار (Pressure Loop) بر روی ایرفویل، انطباق تقریباً کاملی با داده‌های آزمایشگاهی نشان داد. این یعنی مکان دقیق تشکیل شوک و توزیع فشار روی سطوح بالا و پایین به درستی توسط نرم‌افزار تسخیر شده است.

جمع‌بندی و نگاه به آینده شبیه‌سازی NACA 0012 به ما یادآوری می‌کند که شبیه‌سازی CFD، فراتر از فشردن چند دکمه، درک عمیق فیزیک و هندسه است. از اصلاح مقیاس مش گرفته تا انتخاب دقیق فاکتورهای زمانی و مدل‌های حرارتی، هر پارامتری وزنی در نتیجه نهایی دارد.

در پایان، به عنوان یک متخصص این سوال را مطرح می‌کنم: “با توجه به اینکه افزودن جزئیات بیشتر (مانند لبه فرار با ضخامت واقعی یا مدل‌های گذار جریان) هزینه محاسباتی را بالا می‌برد، مرز میان دقت کافی و پیچیدگی بیهوده در یک پروژه مهندسی کجاست؟“ پاسخ به این سوال، تفاوت میان یک تکنسین نرم‌افزار و یک مهندس تحلیل‌گر را مشخص می‌کند.