شبیه‌سازی ایرفویل در ANSYS Fluent: تحلیل جریان و نیروی آیرودینامیکی

ایرفویل‌ها (Airfoils) شکل‌های آیرودینامیکی هستند که اساس طراحی بال هواپیماها، پره‌های توربین‌های بادی و حتی پره‌های پمپ‌ها و کمپرسورها را تشکیل می‌دهند. هدف اصلی طراحی یک ایرفویل، بهینه‌سازی نسبت نیروی برآ (Lift) به نیروی پسا (Drag) است. شبیه‌سازی دینامیک سیالات محاسباتی (CFD) با استفاده از نرم‌افزاری مانند ANSYS Fluent، ابزاری قدرتمند برای تحلیل دقیق عملکرد ایرفویل در شرایط مختلف پروازی یا کاری است. این شبیه‌سازی به مهندسان اجازه می‌دهد تا توزیع فشار، میدان سرعت و نیروهای آیرودینامیکی را بدون نیاز به آزمایش‌های پرهزینه تونل باد، بررسی کنند.

۱. ایجاد هندسه و دامنه جریان (Domain Generation)

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

A. هندسه ایرفویل و مدل‌سازی دوبعدی/سه‌بعدی

• وارد کردن مختصات (Airfoil Coordinates): ایرفویل‌ها معمولاً با استفاده از مختصات استاندارد (مانند سری NACA) تعریف می‌شوند. این مختصات به صورت یک جدول از نقاط $(x,y)$(x,y) وارد نرم‌افزار طراحی هندسه (مانند ANSYS SpaceClaim یا DesignModeler) می‌شوند تا مقطع ایرفویل به صورت یک منحنی (Curve) ایجاد شود.
• ایجاد دامنه سیال (Fluid Domain): از آنجایی که جریان باید به طور کامل مدل شود، ایرفویل باید در داخل یک دامنه جریان بزرگتر قرار گیرد. این دامنه باید به اندازه‌ای بزرگ باشد که مرزهای خارجی آن بر روی جریان اطراف ایرفویل تأثیر نگذارند (معمولاً حداقل ۲۰ برابر طول وتر ایرفویل در جلو و پشت). در شبیه‌سازی‌های دوبعدی، این دامنه یک مستطیل یا دایره بزرگ است که حفره ایرفویل از آن خارج شده است.

B. تعریف نام‌گذاری‌ها (Named Selections)

نام‌گذاری صحیح مرزهای دامنه برای اعمال شرایط مرزی در Fluent حیاتی است:

• Inlet (ورودی): مرزهایی که جریان سیال به داخل دامنه وارد می‌شود (معمولاً سرعت ثابت یا فشار کل).
• Outlet (خروجی): مرزی که جریان سیال از دامنه خارج می‌شود (معمولاً فشار استاتیک صفر یا فشارسنج).
• Airfoil Wall (دیواره ایرفویل): سطح ایرفویل که در آن شرایط بدون لغزش (No-slip condition) اعمال می‌شود.
• Symmetry/Farfield (مرزهای دور): مرزهای بیرونی دامنه که باید به گونه‌ای تعریف شوند تا اثرات انسداد تونل باد را به حداقل برسانند.

۲. مش‌بندی دامنه جریان (Mesh Generation)

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

A. اهمیت لایه‌های مرزی (Inflation Layers)

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

• کنترل $y^{+}$y+: ارتفاع اولین لایه سلول (First Layer Height) باید بر اساس مقدار $y^{+}$y+ مورد نظر تنظیم شود. انتخاب $y^{+}$y+ به مدل توربولانس بستگی دارد:
• مدل‌های مبتنی بر دیوار (Wall-Function Models): $y^{+}$y+ باید بین ۳۰ تا ۳۰۰ باشد.
• مدل‌های با وضوح بالا (Low Reynolds Number Models یا $k-\omega$k−ω SST): $y^{+}$y+ باید نزدیک به ۱ باشد تا لایه مرزی به طور کامل حل شود.
• تراکم مش: مش باید در ناحیه لبه حمله (Leading Edge) و لبه فرار (Trailing Edge) که تغییرات جریان بسیار زیاد است، متراکم‌تر باشد.

B. استفاده از مش غیرساختاریافته/هیبرید

اغلب، مش‌بندی دامنه جریان ایرفویل به صورت هیبرید انجام می‌شود:

• ساختاریافته (Structured): استفاده از سلول‌های چهارضلعی (Quad) در نزدیکی دیواره ایرفویل و در ناحیه بیدار (Wake Region) برای حفظ کیفیت بالا.
• غیرساختاریافته (Unstructured): استفاده از سلول‌های مثلثی (Tri) در بقیه دامنه که ساختار جریان ساده‌تر است.

۳. تنظیمات حل‌گر در ANSYS Fluent

پس از مش‌بندی، نوبت به تعریف مدل‌های فیزیکی و تنظیمات حل در Fluent می‌رسد.

A. مدل‌های فیزیکی (Physical Models)

• انتخاب حل‌گر (Solver): برای جریان‌های تراکم‌ناپذیر و سرعت‌های زیرصوت (مانند هواپیماهای تجاری)، حل‌گر مبتنی بر فشار (Pressure-Based) معمولاً انتخاب می‌شود.
• مدل توربولانس (Turbulence Model): انتخاب مدل توربولانس برای جریان‌های دارای جدایش و لایه‌های مرزی بسیار حیاتی است.
• $k-\omega$k−ω SST: این مدل به طور گسترده‌ای برای شبیه‌سازی ایرفویل‌ها استفاده می‌شود زیرا در حل لایه مرزی (ناحیه $k-\omega$k−ω) و ناحیه دوردست (ناحیه $k-ϵ$k−ϵ) عملکرد قوی دارد و برای پیش‌بینی جدایش جریان دقیق‌تر است.

B. شرایط عملیاتی و مرزی (Operating and Boundary Conditions)

• شرایط عملیاتی: تعریف خواص سیال (چگالی، ویسکوزیته) و فشار مرجع.
• زاویه حمله (Angle of Attack - AoA): در Fluent، زاویه حمله با تنظیم مؤلفه‌های سرعت ورودی (Velocity Components) در مرز ورودی (Inlet) تعیین می‌شود.
• اگر سرعت $V$V باشد و زاویه حمله $\alpha$α باشد، مؤلفه‌های سرعت به صورت $V_x=V\cos \left(\alpha \right)$Vx​=Vcos(α) و $V_y=V\sin \left(\alpha \right)$Vy​=Vsin(α) تعریف می‌شوند.

C. مانیتورینگ (Monitoring)

برای اطمینان از صحت و همگرایی نتایج، باید پارامترهای کلیدی زیر را مانیتور کرد:

• نیروهای آیرودینامیکی: مانیتور کردن ضرایب برآ ($C_L$CL​) و پسا ($C_D$CD​) در هر تکرار حل (Iteration). حل زمانی همگرا شده است که این ضرایب به مقداری ثابت برسند.
• باقیمانده‌ها (Residuals): باقیمانده‌های معادلات پیوستگی، مومنتوم و انرژی باید به مقادیر پایین (معمولاً ${10}^{-4}$10−4 تا ${10}^{-6}$10−6) کاهش یابند.

۴. تحلیل و تفسیر نتایج

پس از همگرایی حل، می‌توان نتایج CFD را استخراج و تحلیل کرد.

A. محاسبه ضرایب آیرودینامیکی

ANSYS Fluent به طور خودکار ضرایب برآ ($C_L$CL​) و پسا ($C_D$CD​) را بر اساس توزیع فشار و تنش برشی روی سطح ایرفویل محاسبه می‌کند. این مقادیر می‌توانند به صورت تابعی از زاویه حمله در نمودار $C_L$CL​ در مقابل $\alpha$α رسم شوند.

B. تحلیل توزیع فشار و میدان سرعت

• توزیع فشار (Pressure Distribution): بررسی کانتور فشار روی سطح ایرفویل. فشار کمتر در سطح بالایی (Suction side) نیروی برآ را تولید می‌کند.
• کانتور سرعت (Velocity Contours): بررسی سرعت جریان در اطراف ایرفویل و به ویژه در ناحیه بیدار (Wake) پشت ایرفویل.
• بررسی جدایش جریان (Flow Separation): جدایش جریان در زوایای حمله بالا (Stall) با مشاهده خطوط جریان (Streamlines) در سطح ایرفویل مشخص می‌شود، که به شدت نیروی برآ را کاهش و پسا را افزایش می‌دهد.

C. تحلیل لایه مرزی

بررسی پروفایل‌های سرعت و مقادیر $y^{+}$y+ در نزدیکی دیواره برای اطمینان از اینکه مش‌بندی به درستی لایه مرزی را حل کرده است.

کلیدواژه ها : ایرفویل-Airfoil-انسیس فلوئنت-ANSYS Fluent-شبیه سازی CFD-CFD Simulation-نیروی برآ-Lift-نیروی پسا-Drag-زاویه حمله-Angle of Attack-AoA-لایه مرزی-Boundary Layer-اینفلیشن لایر-Inflation Layer-کیفیت مش-Mesh Quality-yplus-مدل توربولانس-Turbulence Model-k-omega SST-جدایش جریان-Flow Separation-ضریب برآ-Lift Coefficient-ضریب پسا-Drag Coefficient-هندسه ایرفویل-Airfoil Geometry-دامنه جریان-Fluid Domain-DesignModeler-SpaceClaim-NACA