1. مقدمه: آشنایی با تحلیل آیرودینامیکی ایرفویل

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

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

• 2: حداکثر انحنای (camber) ایرفویل برابر با ۲٪ طول وتر (chord) است.
• 4: این انحنای حداکثری در فاصله ۴۰٪ از لبه حمله (leading edge) قرار دارد.
• 12: حداکثر ضخامت ایرفویل ۱۲٪ طول وتر آن است.

هدف اصلی این مقاله، ارائه یک راهنمای گام‌به‌گام و جامع برای شبیه‌سازی ایرفویل با فلوئنت است. این راهنما همه مراحل کلیدی، از آماده‌سازی و ساخت هندسه در محیط Ansys Design Modeler تا تولید یک شبکه محاسباتی باکیفیت، تنظیمات دقیق حلگر در Ansys Fluent، و در نهایت، تحلیل و استخراج نتایج معنادار را پوشش می‌دهد. در ادامه، هر یک از این مراحل مفصل بررسی خواهند شد تا خوانندگان بتوانند با اطمینان، تحلیل‌های مشابهی را برای پروژه‌های خود به انجام رسانند.

2. پیش‌نیازها و آماده‌سازی اولیه

هر پروژه موفق دینامیک سیالات محاسباتی با یک مرحله آماده‌سازی دقیق و منظم آغاز می‌شود. این مرحله شامل گردآوری داده‌های ورودی صحیح و اطمینان از دسترسی به ابزارهای نرم‌افزاری مناسب است. برای شبیه‌سازی یک ایرفویل، داشتن مختصات دقیق پروفایل آن و درک کلی از ابعاد و شکل دامنه حل (Fluid Domain) از اهمیت بالایی برخوردار است. آماده‌سازی صحیح در این گام، پایه‌ای محکم برای مراحل بعدی یعنی ساخت هندسه و تولید شبکه فراهم می‌کند.

برای اجرای این پروژه، پیش‌نیازهای زیر ضروری است:

• نرم‌افزار:

  دسترسی به مجموعه Ansys Workbench که شامل ابزارهای زیر است:

  • Design Modeler برای ساخت هندسه
  • Meshing برای تولید شبکه محاسباتی
  • Fluent برای تنظیمات فیزیکی، حل و تحلیل نتایج

• فایل مختصات ایرفویل:

  مختصات دقیق پروفایل ایرفویل NACA 2412 که از منبع معتبر ناسا (19950002355.pdf، شکل ۱) استخراج شده است. این مختصات باید در یک فایل متنی ساده (با پسوند .txt یا .dat) و با فرمت‌بندی مناسب برای وارد کردن به نرم‌افزار Design Modeler آماده شود. این فرمت معمولاً شامل ستون‌هایی برای شماره گروه، شماره نقطه، و مختصات X, Y, و Z است.

با آماده‌سازی این موارد، می‌توانیم اولین گام عملی پروژه، یعنی ساخت هندسه دقیق ایرفویل و دامنه سیال اطراف آن را آغاز کنیم.

3. گام اول: ساخت هندسه در Ansys Design Modeler

توپولوژی C-Grid استاندارد صنعتی برای شبیه‌سازی‌های آیرودینامیکی خارجی است، زیرا امکان ایجاد یک شبکه ساختاریافته و بسیار کارآمد را فراهم می‌کند. این ساختار به طور کامل با انحنای لبه حمله ایرفویل تطبیق می‌یابد و همزمان اجازه می‌دهد یک بلوک شبکه متراکم با نسبت ابعادی بالا به سمت پایین‌دست جریان گسترش یابد تا فیزیک دنباله (wake) با دقت بالا ثبت شود؛ همه اینها بدون اتلاف سلولی که یک دامنه مستطیلی ساده به همراه دارد، انجام می‌شود. دامنه محاسباتی باید به اندازه‌ای بزرگ باشد که مرزهای آن تأثیری بر میدان جریان نزدیک به ایرفویل نداشته باشند.

مراحل ساخت هندسه در محیط Design Modeler به شرح زیر است:

1. ورود به Design Modeler:

   ابتدا از طریق Ansys Workbench، یک پروژه Fluid Flow (Fluent) ایجاد کرده و با دو بار کلیک بر روی Geometry وارد محیط Design Modeler شوید. در ابتدای ورود، واحد طول را بر روی متر (Meter) تنظیم کنید.

2. ایجاد پروفایل ایرفویل:

   ابتدا، با استفاده از Concept > 3D Curve فایل مختصات را وارد کنید. پس از تولید Line Body منحنی، آن را انتخاب کرده و با استفاده از Concept > Surfaces from Edges یک سطح صفحه‌ای (planar surface) محصور در پروفایل ایرفویل ایجاد کنید. در نهایت، این سطح جدید را انتخاب کرده و با ابزار Tools > Fill یک بدنه سطحی صلب (solid surface body) بسازید.

3. طراحی دامنه سیال (C-Grid):

   یک دامنه استاندارد C-Grid باید تقریباً ۱۲ برابر طول وتر (12c) در بالادست و ۲۴ برابر طول وتر (24c) در پایین‌دست امتداد یابد تا اثرات مرزها به حداقل برسد. این دامنه را می‌توان با ترسیم یک نیم‌دایره با شعاع 12c به مرکزیت لبه حمله (با ابزار Arc by Center) و گسترش آن با یک بخش مستطیلی به طول 24c (با ابزار Rectangle) ایجاد کرد.

4. ایجاد حجم سیال:

   پس از ساخت سطح دامنه سیال و سطح ایرفویل، از عملیات بولین Create > Boolean استفاده کنید. عملیات را روی حالت Subtract تنظیم کنید. نکته حیاتی: اطمینان حاصل کنید که دامنه سیال به عنوان Target Body و ایرفویل به عنوان Tool Body انتخاب شده است. معکوس کردن این انتخاب به اشتباه دامنه را از ایرفویل کم خواهد کرد. با این کار، سطح ایرفویل از دامنه کلی کم شده و فضای نهایی جریان سیال ایجاد می‌شود.

5. تقسیم دامنه:

   این تقسیم‌بندی به چهار ربع (quadrants)، یک پیش‌نیاز حیاتی برای اعمال Mapped Face Meshing است. این کار به الگوریتم شبکه‌بندی، نواحی منطقی و چهارضلعی می‌دهد که می‌توانند با یک شبکه چهارضلعی کاملاً ساختاریافته پر شوند و المان‌های کج (skewed elements) که ممکن است در یک دامنه پیچیده توسط الگوریتم‌های بدون ساختار ایجاد شوند را حذف می‌کند. این کار معمولاً با ترسیم خطوطی از لبه‌های حمله و فرار ایرفویل به سمت مرزهای دامنه و سپس استفاده از ابزار Slice انجام می‌شود.

با تکمیل این مراحل، هندسه مسئله نهایی شده و آماده ورود به مرحله حیاتی بعدی، یعنی تولید شبکه محاسباتی است.

4. گام دوم: تولید شبکه محاسباتی در Ansys Meshing

کیفیت شبکه محاسباتی (Mesh) نقشی حیاتی در دقت، پایداری و هزینه محاسباتی نتایج CFD ایفا می‌کند. هدف، ایجاد یک شبکه ساختاریافته (Structured Mesh) با تراکم بالا در نزدیکی سطح ایرفویل برای ثبت دقیق لایه مرزی و در ناحیه دنباله (wake) است. دستیابی به مقدار Y+ کمتر از ۱ صرفاً یک توصیه نیست؛ بلکه یک الزام قطعی برای عملکرد صحیح مدل توربولانسی k-ω SST است که در ادامه انتخاب خواهد شد. این مدل برای حل مستقیم زیرلایه لزج (viscous sublayer) در لایه مرزی فرموله شده است. اگر Y+ بیش از حد بالا باشد، اولین گره شبکه خارج از این ناحیه قرار می‌گیرد و منجر به محاسبات نادرست تنش برشی دیواره و در نتیجه، پیش‌بینی غیردقیق نیروی پسا می‌شود.

فرآیند گام‌به‌گام تولید شبکه در محیط Ansys Meshing به شرح زیر است:

1. ورود به Ansys Meshing:

   از محیط Workbench، بر روی Mesh دو بار کلیک کنید تا محیط شبکه‌بندی باز شود.

2. ایجاد نام برای مرزها (Named Selections):

   با استفاده از ابزار انتخاب لبه (Edge Selection)، مرزهای اصلی دامنه را انتخاب کرده و با کلیک راست، برای آن‌ها نام‌های مشخصی تعریف کنید. این نام‌ها در Fluent برای اعمال شرایط مرزی استفاده می‌شوند و باید با دقت تعریف شوند:

   • inlet: کل مرز خارجی C-شکل که از بالای خروجی شروع شده، دور لبه حمله می‌چرخد و به پایین خروجی می‌رسد.
   • outlet: مرز عمودی در انتهای بخش مستطیلی دامنه.
   • airfoil_wall: دو لبه‌ای که پروفایل سطح ایرفویل را تشکیل می‌دهند.
   • symmetry: لبه افقی که روی محور X از لبه فرار ایرفویل تا مرز خروجی امتداد دارد. این شرط تنها در صورتی معتبر است که نیمی از دامنه مدل‌سازی شود و جریان‌متقارن فرض شود، که برای زاویه حمله صفر قابل قبول است.

3. اعمال کنترل‌های شبکه:

   • Mapped Face Meshing: بر روی هر چهار سطح (quadrant) دامنه، کنترل Mapped Face Meshing را اعمال کنید. این کار به نرم‌افزار دستور می‌دهد تا یک شبکه کاملاً چهارضلعی و ساختاریافته ایجاد کند که برای این نوع مسائل ایده‌آل است.
   • Edge Sizing: با استفاده از این کنترل، تعداد تقسیمات المان‌ها را بر روی لبه‌های مختلف تعریف کنید. بر روی لبه‌های مربوط به پروفایل ایرفویل، تعداد تقسیمات بالایی (مثلاً ۲۰۰ تقسیم) و بر روی لبه‌های دورتر مانند مرز خروجی، تعداد کمتری (مثلاً ۵۰ تقسیم) در نظر بگیرید تا از افزایش بی‌رویه تعداد المان‌ها جلوگیری شود.
   • Bias Factor: Bias Factor یک ابزار حیاتی برای بهره‌وری محاسباتی است. با ایجاد یک تصاعد هندسی در اندازه سلول‌ها، می‌توانیم یک شبکه با وضوح بالا را در جایی که ضروری است – یعنی در لایه مرزی و دنباله – حفظ کنیم و همزمان آن را در میدان دور که گرادیان‌ها ناچیز هستند، به سرعت درشت‌تر کنیم. این استراتژی به شدت تعداد کل سلول‌ها و زمان حل را بدون به خطر انداختن دقت محاسبات نیرو کاهش می‌دهد.

4. تولید و بازبینی شبکه:

   پس از اعمال تمام کنترل‌ها، بر روی Generate Mesh کلیک کنید. شبکه تولید شده را به صورت چشمی بازبینی کرده و از کیفیت آن در نواحی کلیدی اطمینان حاصل کنید. به عنوان یک بررسی کیفی، حداقل Orthogonal Quality باید بالاتر از 0.15 و حداکثر Skewness باید زیر 0.85 باقی بماند تا از پایداری عددی و دقت حل اطمینان حاصل شود. شبکه نهایی باید دارای المان‌های کشیده و منظم در اطراف ایرفویل و در ناحیه دنباله باشد.

با آماده شدن یک شبکه محاسباتی باکیفیت، اکنون همه چیز برای تعریف فیزیک مسئله و تنظیمات حلگر در Ansys Fluent مهیا است.

5. گام سوم: تنظیمات فیزیکی مسئله در Ansys Fluent

این مرحله را می‌توان قلب شبیه‌سازی نامید، زیرا در اینجا معادلات حاکم بر جریان، خواص سیال و شرایط مرزی که رفتار فیزیکی واقعی مسئله را تعریف می‌کنند، تنظیم می‌شوند. دقت و صحت نتایج نهایی به طور مستقیم به درستی تنظیمات انجام شده در این بخش بستگی دارد. انتخاب مدل توربولانسی مناسب و تعریف دقیق شرایط مرزی، کلید دستیابی به یک شبیه‌سازی معتبر و قابل اتکا است.

راهنمای گام‌به‌گام تنظیمات در محیط Ansys Fluent به شرح زیر است:

1. راه‌اندازی Fluent:

   از محیط Workbench، با دو بار کلیک بر روی Setup وارد Fluent Launcher شوید. گزینه Double Precision را فعال کنید تا از اعداد ممیز شناور ۶۴-بیتی استفاده شود؛ این کار خطاهای گردکردن عددی که می‌توانند بر همگرایی محاسبات حساس نیرو تأثیر بگذارند را به حداقل می‌رساند. سپس حلگر را راه‌اندازی کنید.

2. تنظیمات عمومی حلگر (General Setup):

   در پنل General، حلگر Pressure-Based انتخاب می‌شود، زیرا برای جریان‌های با سرعت پایین و تراکم‌ناپذیر (عدد ماخ بسیار کمتر از 0.3) که در آن چگالی ثابت است، بهینه شده است. الگوریتم تصحیح فشار آن برای این رژیم جریان، پایدارتر و از نظر محاسباتی قوی‌تر از حلگر Density-Based است که برای جریان‌های پرسرعت و تراکم‌پذیر طراحی شده. حالت زمانی را نیز Steady (پایا) در نظر بگیرید.

3. انتخاب مدل فیزیکی (Models):

   • مدل توربولانسی: به پنل Viscous بروید و مدل k-ω SST را انتخاب کنید. این مدل ترکیبی هوشمندانه از دقت بالای مدل k-ω در نواحی نزدیک به دیواره (برای ثبت دقیق لایه مرزی) و پایداری و کارایی مدل k-ε در جریان آزاد (دور از دیواره‌ها) است. به همین دلیل، مدل k-ω SST برای شبیه‌سازی جریان‌های خارجی همراه با گرادیان فشار معکوس و جدایش احتمالی جریان، مانند جریان حول ایرفویل، گزینه‌ای بسیار مناسب و پرکاربرد است.

4. تعریف مواد (Materials): در بخش Materials، اطمینان حاصل کنید که سیال air با خواص پیش‌فرض (چگالی 1.225 kg/m³ و ویسکوزیته 1.7894e-05 kg/m-s) تعریف شده است.

5. شرایط مرزی (Boundary Conditions):

   مقادیر دقیق را برای هر یک از مرزهایی که در مرحله شبکه‌بندی نام‌گذاری کردید، به شرح زیر تنظیم کنید:

   • inlet: نوع آن را velocity-inlet قرار دهید. سرعت ورودی (Velocity Magnitude) برابر با 40 m/s تنظیم شده است که با پارامترهای پروژه مرجع سازگار است. لازم به ذکر است که داده‌های تجربی مورد استفاده برای اعتبارسنجی ممکن است در اعداد رینولدز کمی متفاوت باشند. برای این شبیه‌سازی اولیه، زاویه حمله صفر درجه در نظر گرفته می‌شود.
   • outlet: نوع آن را pressure-outlet تنظیم کرده و مقدار فشار نسبی (Gauge Pressure) را برابر با صفر پاسکال قرار دهید (یعنی فشار اتمسفر).
   • airfoil_wall: نوع آن باید wall باشد و شرط مرزی عدم لغزش (No-Slip) به طور پیش‌فرض برای آن اعمال می‌شود.
   • symmetry: نوع این مرز را symmetry انتخاب کنید.

6. تنظیمات حل (Solution):

   • Solution Methods: برای الگوریتم فشار-سرعت، روش جفت‌شده (Coupled) را انتخاب کنید. این روش معادلات مومنتوم و فشار را به طور همزمان حل می‌کند و منجر به همگرایی سریع‌تر برای مسائل پایا مانند این مسئله می‌شود، هرچند هزینه حافظه در هر تکرار بالاتر است.
   • Monitors: به بخش Residuals بروید و معیار همگرایی (Convergence Criterion) را برای تمام معادلات بر روی 1e-6 تنظیم کنید. این کار تضمین می‌کند که حل تا رسیدن به دقت عددی بالا ادامه یابد.

پس از تکمیل دقیق تمامی این تنظیمات، شبیه‌سازی برای مرحله حل و محاسبه میدان جریان آماده است.

6. گام چهارم: اجرای شبیه‌سازی و تحلیل نتایج

پس از تنظیم کامل فیزیک مسئله، فرآیند حل آغاز می‌شود. حلگر Fluent از یک رویکرد تکراری (iterative) برای حل معادلات دیفرانسیل جزئی حاکم بر جریان استفاده می‌کند تا زمانی که باقیمانده (Residual) معادلات به زیر معیار همگرایی تعریف‌شده کاهش یابد. رسیدن به همگرایی به معنای آن است که حل به یک حالت پایدار و از نظر عددی ثابت رسیده است. پس از آن، نتایج خام باید به داده‌های بصری و عددی معنادار تبدیل شوند.

6.1. اجرای حل

• Initialization: پیش از شروع تکرارها، باید یک حدس اولیه برای میدان جریان تعریف شود. از روش Hybrid Initialization استفاده کنید. Hybrid Initialization به طور کلی قوی است، اما برای مواردی با گرادیان‌های سرعت بالا یا فیزیک پیچیده، توصیه می‌شود که میدان‌های اولیه (مانند فشار و سرعت) را بازرسی کنید تا از معقول بودن فیزیکی آنها قبل از شروع محاسبه اطمینان حاصل شود.
• Run Calculation: در پنل Run Calculation، تعداد تکرارها را روی یک مقدار مناسب (مثلاً ۵۰۰ تکرار) تنظیم کرده و بر روی Calculate کلیک کنید. در حین اجرای حل، نمودار باقیمانده‌ها را به دقت زیر نظر بگیرید تا از روند نزولی و رسیدن آن‌ها به زیر معیار 1e-6 اطمینان حاصل کنید.

6.2. تحلیل نتایج (Post-Processing)

پس از اتمام حل، می‌توانید نتایج را به روش‌های مختلفی تحلیل کنید:

• ایجاد کانتورهای فشار و سرعت: از بخش Graphics > Contours برای نمایش کانتورهای فشار استاتیک (Static Pressure) و سرعت (Velocity Magnitude) در اطراف ایرفویل استفاده کنید. همانطور که انتظار می‌رود، باید یک ناحیه با فشار استاتیک بالا در نقطه سکون (لبه حمله) و یک ناحیه با فشار پایین و سرعت بالا بر روی سطح بالایی ایرفویل مشاهده شود. این اختلاف فشار، عامل اصلی تولید نیروی برآ (Lift) است.
• ترسیم خطوط جریان (Streamlines): برای درک بهتر مسیر حرکت ذرات سیال و مشاهده نحوه عبور جریان از روی ایرفویل، از ابزار Streamline استفاده کنید. این خطوط به وضوح نشان می‌دهند که جریان به سطح ایرفویل چسبیده باقی می‌ماند یا در نقاطی دچار جدایش می‌شود.
• محاسبه ضرایب آیرودینامیکی: برای استخراج مقادیر عددی ضرایب برآ (Lift Coefficient) و پسا (Drag Coefficient)، به بخش Reports > Forces > Print بروید. در این پنجره، جهت بردارهای لیفت (عمود بر جریان) و درگ (موازی با جریان) را مشخص کرده و airfoil_wall را به عنوان سطحی که نیرو بر آن محاسبه می‌شود، انتخاب کنید. با کلیک بر روی Print، مقادیر عددی این ضرایب در کنسول نرم‌افزار نمایش داده می‌شود.

در پایان، باید به این نکته توجه داشت که نتایج به دست آمده از شبیه‌سازی باید با داده‌های معتبر تجربی یا تحلیلی مقایسه شوند تا از صحت و اعتبار مدل CFD اطمینان حاصل شود.

7. جمع‌بندی و نتیجه‌گیری

در این مقاله، یک راهنمای جامع و گام‌به‌گام برای شبیه‌سازی جریان حول ایرفویل NACA 2412 با استفاده از مجموعه نرم‌افزاری Ansys ارائه شد. مراحل کلیدی پروژه، از جمله ساخت هندسه دقیق و دامنه C-Grid، تولید شبکه ساختاریافته با کنترل‌های پیشرفته، تنظیم مدل فیزیکی k-ω SST در Fluent و استخراج نتایج کلیدی مانند ضرایب آیرودینامیکی، به تفصیل پوشش داده شد. این فرآیند بار دیگر نشان داد که شبیه‌سازی ایرفویل با فلوئنت یک ابزار مهندسی قدرتمند است که به طراحان اجازه می‌دهد عملکرد آیرودینامیکی را با دقت و سرعت بالا پیش‌بینی و بهینه کنند.

مهم‌ترین یافته‌های این شبیه‌سازی برای زاویه حمله صفر، شامل مشاهده توزیع فشار نامتقارن حول ایرفویل (به دلیل وجود انحنا یا camber) بود که منجر به تولید نیروی برآی مثبت حتی در غیاب زاویه حمله می‌شود. همچنین کانتورهای سرعت، شتاب گرفتن جریان بر روی سطح بالایی و کاهش سرعت در لبه حمله را به وضوح به تصویر کشیدند.

به عنوان یک توصیه نهایی و گامی برای اعتبارسنجی مدل، پیشنهاد می‌شود که این شبیه‌سازی برای زوایای حمله مختلف (مانند ۸، ۱۲، و ۱۶ درجه) تکرار شود. سپس، نتایج ضرایب لیفت و درگ به دست آمده را با داده‌های تجربی موجود در منبع 19950002355.pdf (شکل 6) مقایسه کنید. این مقایسه نه تنها به شما درک عمیق‌تری از رفتار آیرودینامیکی ایرفویل در شرایط نزدیک به واماندگی می‌دهد، بلکه اعتبار و دقت مدل CFD شما را نیز به طور عملی تأیید می‌کند.