آموزش روش MRF (دستگاه مختصات متحرک) در فلوئنت برای شبیهسازی فن (راهنمای ۰ تا ۱۰۰)
اگر تا به حال سعی کرده باشید جریان هوا را اطراف یک فن یا ملخ هواپیما شبیهسازی کنید، حتماً با این دوراهی کلافهکننده روبرو شدهاید: یا باید هفتهها وقت بگذارید تا یک حل “مش لغزان” (Sliding Mesh) با گام زمانی کوچک جلو برود، یا اینکه دنبال راهی باشید که در چند ساعت، یک نتیجه مهندسی قابل قبول بگیرید. واقعیت این است که در صنعت، همیشه زمان لوکسترین دارایی ماست و اینجاست که آموزش روش MRF (دستگاه مختصات متحرک) در فلوئنت برای شبیهسازی فن مثل یک فرشته نجات ظاهر میشود.
مقایسه فنی روشMRF وSliding Mesh (راهنمای تصمیمگیری)
| ویژگی | روش MRF (دستگاه مختصات متحرک) | روش Sliding Mesh (مش لغزان) |
| نوع حلگر | پایا (Steady-State) | گذرا (Transient) |
| هزینه محاسباتی | پایین (چند ساعت) | بسیار بالا (چند روز تا چند هفته) |
| دقت نتایج | خوب برای نقاط طراحی و راندمان کلی | عالی برای تمام شرایط حتی استال و نویز |
| حرکت مش | مش ثابت است (فریم میچرخد) | مش واقعاً میچرخد (Physically Rotating) |
| کاربرد اصلی | استخراج منحنی عملکرد، طراحی اولیه | تحلیل آکوستیک، تداخل روتور-استاتور، راهاندازی |
| پیشنیاز سختافزاری | سیستم متوسط مهندسی | سیستمهای قوی یا سرورهای HPC |
قبل از اینکه وارد جزئیات فنی شویم، اگر تازه کار با این نرمافزار را شروع کردهاید یا میخواهید یک دید کلی نسبت به ماژولهای مختلف آن داشته باشید، پیشنهاد میکنم حتماً سری به راهنمای کامل انسیس فلوئنت (Ansys Fluent): راهنمای جامع از مقدماتی تا پیشرفته بزنید تا نقشهی راه دستتان بیاید. اما اگر آمادهاید تا یکبار برای همیشه پرونده شبیهسازی فن در حالت پایا (Steady) را ببندید، با من همراه باشید.
۱. چرا برای شبیهسازی فن در فلوئنت باید به جای مش لغزان از روش MRF استفاده کنیم؟
بیایید روراست باشیم؛ همه ما دوست داریم دقیقترین شبیهسازی ممکن را انجام دهیم. اما وقتی مدیر پروژه بالای سرتان ایستاده و نتایج را “برای دیروز” میخواهد، اجرای یک حل گذرا (Transient) که هر دور چرخش فن را به ۳۶۰ گام زمانی تقسیم میکند، عملاً خودکشی است! روش MRF (Multiple Reference Frame) دقیقاً برای همین روزهای مبادا ساخته شده.
در این روش، ما به جای اینکه مش را واقعاً بچرخانیم (کاری که در Sliding Mesh انجام میشود و پدر CPU را درمیآورد)، معادلات را در یک دستگاه مختصات چرخان حل میکنیم. این یعنی شما میتوانید با یک حلگر Steady-state، اثرات چرخش را با دقت بسیار خوبی (معمولاً بالای ۹۰٪ برای نقاط طراحی) بدست آورید. البته برای درک بهتر اینکه چه زمانی باید سراغ کدام روش بروید، مرور [مراحل گام به گام انجام یک پروژه در فلوئنت] دید خوبی به شما میدهد تا استراتژی درستی بچینید. در تجربه من، برای ۹۰ درصد پروژههای صنعتی مثل تهویه مطبوع یا فنهای رادیاتور، MRF کاملاً کفایت میکند و نیازی به هزینههای سنگین محاسباتی نیست.
۲. تفاوت فیزیکی روش دستگاه مختصات متحرک (MRF) با روشهای واقعی چرخش پره چیست؟
تصور کنید یک دوربین عکاسی دارید و از یک پنکه روشن با سرعت شاتر بسیار بالا عکس میگیرید. چه میبینید؟ پرهها در جای خودشان “منجمد” شدهاند، اما ما میدانیم که هوا در حال حرکت است. روش MRF دقیقاً همین کار را میکند؛ به آن روش “Frozen Rotor” هم میگویند.
در واقعیت، پره دارد محیط سیال را میشکافد. اما در MRF، ما پره را ثابت نگه میداریم و به جایش به کل فضای اطرافش (در معادلات ریاضی) سرعت زاویهای میدهیم. این یعنی نیروهای کوریولیس و گریز از مرکز به معادلات تکانه اضافه میشوند. درک این مفهوم خیلی مهم است چون اگر دنبال بررسی گردابههایی هستید که با زمان تغییر میکنند (مثلاً در بحث آکوستیک)، این روش شاید کمی گولت بزند. برای اینکه بهتر متوجه شوید سیال چطور رفتار میکند، خواندن مقاله [آیرودینامیک به زبان ساده برای مهندسان] میتواند پایههای تئوری شما را محکمتر کند.
۳. چگونه باید ناحیه دوار (Rotating Zone) را در هندسه برای اعمال روش MRF جدا کنیم؟
اینجا همانجایی است که خشت اول پروژه گذاشته میشود. اگر در مرحله هندسه (Geometry) خرابکاری کنید، هیچ تنظیماتی در فلوئنت نجاتتان نمیدهد. برای استفاده از آموزش روش MRF (دستگاه مختصات متحرک) در فلوئنت برای شبیهسازی فن، شما نمیتوانید کل دامنه حل را یکپارچه بگیرید.
باید یک حجم مجزا (معمولاً استوانهای یا مخروطی) دقیقاً اطراف پرههای فن ایجاد کنید. این حجم، همان ناحیهای است که ما بعداً به آن “Frame Motion” میدهیم. نکته طلایی که خیلیها رعایت نمیکنند این است: این ناحیه نباید آنقدر کوچک باشد که به نوک پرهها بچسبد (چون گرادیانها شدید میشود) و نباید آنقدر بزرگ باشد که به دیوارههای ساکن کانال برخورد کند. یک فاصله ایمن (مثلاً ۱۰ تا ۲۰ درصد قطر پره) در نظر بگیرید. 📐
۴. چرا کیفیت مش در مرز اینترفیس (Interface) بین ناحیه ساکن و متحرک حکم مرگ و زندگی را برای پروژه دارد؟
من در طول ۷ سال فعالیتم در حوزه CFD، بارها دیدم که دانشجوها ساعتها روی تنظیمات توربولانس وقت میگذارند اما از کیفیت مش در مرز اتصال ناحیه چرخان به ساکن غافل میشوند. وقتی شما دو ناحیه جداگانه دارید، نرمافزار باید اطلاعات جریان (سرعت، فشار و…) را از سلولهای ناحیه چرخان به سلولهای ناحیه ساکن پاس بدهد.
اگر مش شما در این مرز نامنظم باشد یا اختلاف سایز شدیدی داشته باشد، خطای درونیابی (Interpolation Error) بیداد میکند! نتیجهاش میشود پرشهای غیرعادی در نمودار فشار یا حتی واگرایی حل. سعی کنید تا حد امکان در این نواحی مش ریز و هماندازه داشته باشید. استفاده از [راهنمای کنترل کیفیت مش (Mesh Quality) در فلوئنت] میتواند چکلیست خوبی باشد تا قبل از ایمپورت کردن مش به فلوئنت، از سلامت شبکه مطمئن شوید. یک مش بد مثل فونداسیون کج برای یک ساختمان است؛ هرچقدر هم مصالح (تنظیمات) خوب باشد، باز هم میریزد.
۵. چگونه باید تنظیمات Cell Zone Conditions را برای فعالسازی Frame Motion در ناحیه پره انجام دهیم؟
خب، وارد محیط فلوئنت شدیم. این مرحله قلب تپنده تنظیمات MRF است. به تب Cell Zone Conditions بروید و ناحیهای که برای پره ساخته بودید (همان استوانه دور فن) را انتخاب کنید. اینجا یک اشتباه رایج وجود دارد: تیک Mesh Motion را نزنید! آن برای وقتی است که مش واقعاً حرکت میکند (مثل روش Sliding Mesh).
برای MRF، باید تیک Frame Motion را فعال کنید. با این کار به فلوئنت میگویید: “آقای نرمافزار، این ناحیه از نظر مشبندی ثابت است، اما دستگاه مختصاتش دارد میچرخد.” سپس باید سرعت دورانی (RPM) و محور چرخش (Rotation-Axis Origin/Direction) را تعیین کنید. 🛠️ اگر در تنظیم سایر شرایط مرزی هم شک دارید، مقاله [راهنمای کامل شرایط مرزی (Boundary Conditions) در فلوئنت] مرجع خوبی برای مرور است.
راهنمای تنظیمات پیشنهادی برای شبیهسازی فن محوری(Best Practices)
| پارامتر تنظیماتی | مقدار یا گزینه پیشنهادی | دلیل انتخاب |
| Solver Type | Pressure-Based | مناسب برای جریانهای تراکمناپذیر و سرعتهای زیر صوت |
| Turbulence Model | k-ω SST | پیشبینی دقیق جدایش جریان روی ایرفویل پره |
| Wall Treatment | Enhanced Wall Treatment | بهبود دقت در لایه مرزی (نیاز به y+ زیر 5) |
| Pressure-Velocity Coupling | Coupled | همگرایی پایدارتر در جریانهای چرخشی |
| Initialization | Hybrid Initialization | ایجاد میدان اولیه بهتر برای جلوگیری از واگرایی |
| Discretization | Second Order Upwind | کاهش خطاهای پخش عددی و افزایش دقت |
۶. سرعت دورانی (RPM) و جهت چرخش فن را بر اساس قانون دست راست چگونه در فلوئنت تعیین کنیم؟
شاید باورتان نشود ولی یکی از پرتکرارترین خطاهای مهندسی، حتی بین باتجربهها، اشتباه زدن جهت چرخش است! یادم میآید در یک پروژه تهویه تونل، همه چیز عالی بنظر میرسید اما فن به جای اینکه هوا را بیرون بکشد، داشت داخل میدمید، فقط بخاطر یک علامت منفی!
فلوئنت از قانون دست راست استفاده میکند. انگشت شست خود را در جهت محور چرخش (مثلاً محور Z مثبت) بگیرید؛ جهت بسته شدن انگشتان، جهت چرخش مثبت است. اگر فن شما باید ساعتگرد بچرخد اما محور Z شما رو به بیرون است، احتمالاً باید عدد RPM را منفی وارد کنید. حتماً قبل از ران گرفتن، یک بار بردار سرعت را چک کنید تا مطمئن شوید هوا به سمت درستی هل داده میشود.
۷. تنظیمات دیوارههای پره (Wall) باید روی حالت Moving Wall باشند یا Stationary Wall؟
این بخش کمی گیجکننده است و نیاز به دقت دارد. پرههای فن داخل ناحیه چرخان (Fluid-Rotating) قرار دارند. سوال اینجاست: دیوار پره باید چه تنظیمی داشته باشد؟
از آنجایی که خود پره همراه با هوای داخل آن ناحیه میچرخد، سرعت نسبی (Relative Velocity) پره نسبت به آن ناحیه “صفر” است. پس باید دیواره پره را روی Moving Wall بگذارید اما سرعتش را نسبت به “Adjacent Cell Zone” (ناحیه همسایه) برابر با صفر تنظیم کنید (Rotational Speed = 0). این یعنی پره همپا با دستگاه مختصات میچرخد. این ریزهکاریها در مدلسازی دقیق خیلی موثرند. برای درک بهتر تاثیر این تنظیمات روی جریان، مطالعه [مدل سازی جریان های آشفته: از تئوری تا عمل] به شما کمک میکند بفهمید لایه مرزی روی پره چطور تشکیل میشود.
۸. چگونه شرایط مرزی اینترفیس (Mesh Interface) را برای انتقال صحیح جریان بین نواحی تنظیم کنیم؟
وقتی مش شما Non-Conformal است (یعنی گرههای مش دو ناحیه دقیقا روبروی هم نیستند)، باید دستی اینترفیس تعریف کنید. در نسخههای جدید فلوئنت (۲۰۲۲ به بعد)، معمولاً خود نرمافزار به صورت هوشمند این سطوح تماس را تشخیص میدهد و رابطهای از نوع Mismatched یا Interface میسازد.
اما همیشه خودتان چک کنید. به بخش Domain > Interfaces بروید. باید مطمئن شوید که مرز بیرونی ناحیه دوار به مرز داخلی ناحیه ساکن کوپل شده است. اگر این اتصال برقرار نشود، انگار یک دیوار نامرئی جلوی جریان هوا را گرفته و هیچ بادی از فن خارج نمیشود! در جدول زیر تنظیمات پیشنهادی برای انواع اینترفیس در شبیهسازی فن را برایتان خلاصه کردهام:
| نوع اتصال (Interface Type) | کاربرد اصلی | نکته کلیدی تنظیمات |
| Periodic | وقتی فقط یک پره را شبیهسازی میکنیم | حتماً زاویه تناوب (مثلاً ۳۶۰ تقسیم بر تعداد پره) دقیق باشد. |
| Interior | وقتی مشها کاملاً Conformal (منطبق) هستند | نیازی به تنظیم خاصی نیست، جریان خودکار عبور میکند. |
| Mesh Interface | برای مشهای Non-Conformal (روش رایج) | گزینههای Matching و Mapping را بر اساس کیفیت مش تنظیم کنید. |
تا اینجا نیمی از مسیر آموزش روش MRF (دستگاه مختصات متحرک) در فلوئنت برای شبیهسازی فن را طی کردیم و پیکربندی اصلی مدل انجام شده است. اما هنوز انتخاب مدل توربولانس مناسب و تحلیل نتایج مانده که در ادامه به آن میپردازیم…
۹. کدام مدل آشفتگی (Turbulence Model) بهترین دقت را برای شبیهسازی جریانهای چرخشی شدید ارائه میدهد؟
اینجا جایی است که خیلی از مهندسان تازهکار به دیوار میخورند. اگر فکر میکنید مدل استاندارد k-epsilon برای همه چیز جواب میدهد، سخت در اشتباهید! در جریانهای چرخشی (Swirling Flows) که در فنها داریم، مدل استاندارد به شدت دسیپیشن (تلاف انرژی) را بالا تخمین میزند و نتایج را خراب میکند.
تجربه من و خیلی از همکاران نشان داده که برای روش MRF، مدل Realizable k-epsilon یا k-omega SST بهترین انتخابها هستند. مدل Realizable فرمولاسیون اصلاحشدهای برای لزجت توربولانسی دارد که چرخش را خیلی بهتر میبیند. اگر میخواهید دقیقاً بدانید تفاوت اینها در معادلات ریاضی چیست و چرا یکی بر دیگری برتری دارد، مقاله [بررسی کامل انواع مدل k-epsilon در فلوئنت (Standard, RNG, Realizable)] را حتما بخوانید. البته اگر جدایش جریان روی پره خیلی شدید باشد (مثلاً در زوایای حمله بالا)، مدل SST عملکرد بهتری در نزدیکی دیواره دارد.
۱۰. برای جلوگیری از واگرایی حل در ثانیههای اول شبیهسازی فن چه استراتژیهایی در مقداردهی اولیه (Initialization) وجود دارد؟
لحظهای که دکمه “Calculate” را میزنید و در عرض ۲ ثانیه ارور “Divergence detected” میگیرید، یکی از بدترین حسهای دنیاست! 😤 معمولاً دلیلش این است که سرعت اولیه را صفر گذاشتهاید و ناگهان میخواهید فن را با ۳۰۰۰ دور بر دقیقه بچرخانید؛ خب معلوم است حلگر شوکه میشود.
بهترین راهکار استفاده از Hybrid Initialization است که یک میدان جریان اولیه تقریبی میسازد. اما اگر باز هم جواب نداد، یک ترفند که من همیشه استفاده میکنم این است: ابتدا سرعت چرخش (RPM) را کم کنید (مثلاً ۱۰٪ سرعت نهایی)، حل را همگرا کنید و سپس پلهپله سرعت را زیاد کنید. همچنین دستکاری ضرایب زیر-آرامش (Under-Relaxation Factors) هم میتواند کمک کند. اگر مدام با خطاهای عجیب و غریب مواجه میشوید، پیشنهاد میکنم نگاهی به [۷ دلیل اصلی عدم همگرایی (Divergence) در فلوئنت و راهحل آنها] بیندازید تا ریشه مشکل را پیدا کنید.
۱۱. چگونه میتوانیم گشتاور (Torque) و افت فشار فن را پس از حل جریان محاسبه کنیم؟
شبیهسازی قشنگ است، اما کارفرما از شما “عدد” میخواهد. برای محاسبه گشتاور، باید قبل از حل (یا بعد از آن) به بخش Report Definitions > Moment بروید. نکتهای که اینجا باید حواستون باشه (و خیلیها فراموش میکنند) تنظیم درست “Moment Center” است. اگر مرکز گشتاور را دقیقاً روی محور چرخش فن نگذارید، اعداد فضایی و غلط تحویل میگیرید.
برای افت فشار یا افزایش فشار استاتیک هم باید یک گزارش Surface Integral > Area-Weighted Average روی ورودی و خروجی تعریف کنید. در جدول زیر یک راهنمای سریع برای گزارشگیریهای مهم در فن آوردهام:
| پارامتر مورد نظر | نوع گزارش در فلوئنت | ناحیه انتخاب شده | کاربرد |
| گشتاور (Torque) | Moment | تمام سطوح پره (Blades) | محاسبه توان مصرفی |
| تراست (Thrust) | Force (در جهت محور) | سطوح پره | نیروی جلوبرنده |
| افزایش فشار | Surface Integral (Pressure) | Inlet & Outlet | محاسبه راندمان |
| دبی جرمی | Mass Flow Rate | Outlet | نقطه کاری فن |
۱۲. چرا بررسی کانتورهای سرعت نسبی (Relative Velocity) در ناحیه پره مهمتر از سرعت مطلق است؟
وقتی کانتور سرعت مطلق (Velocity Magnitude) را نگاه میکنید، معمولاً ناحیه پره کاملاً قرمز است چون سرعت خطی بالایی دارد. این چیز زیادی به شما نشان نمیدهد. هنر یک متخصص CFD این است که “سرعت نسبی” (Relative Velocity) را بررسی کند.
در فریمورک [آموزش تکنیکهای حرفهای پسپردازش در CFD-Post]، همیشه تاکید میشود که سرعت نسبی نشان میدهد هوا از دید ناظر روی پره چطور حرکت میکند. با این دید، شما میتوانید نواحی “جدایش جریان” (Separation) روی سطح مکش پره یا گردابههای نوک پره را ببینید. اگر سرعت نسبی در پشت پره نزدیک صفر یا منفی شد، یعنی طراحی ایرودینامیکی شما مشکل دارد و فن در آن ناحیه کارایی ندارد.
۱۳. آیا نتایج روش MRF برای تمام نقاط کاری فن (Fan Curve) قابل اعتماد است یا محدودیتهایی دارد؟
باید با خودمان صادق باشیم؛ روش MRF یک تقریب “پایا” (Steady) از یک پدیده ذاّتاً گذراست. این روش در نقطه طراحی (Design Point) که جریان روی پره چسبیده و منظم است، فوقالعاده دقیق عمل میکند.
اما اگر فن را در شرایط Stall (واماندگی) یا Surge شبیهسازی کنید که جریان ناپایدار است و گردابههای بزرگ جدا میشوند، MRF نمیتواند نوسانات فرکانسی را ببیند و خطا بالا میرود. در این شرایط، نتایج MRF ممکن است کمی خوشبینانه باشد. پس همیشه این محدودیت را در گزارش نهاییتان ذکر کنید تا اعتبار فنیتان حفظ شود.
۱۴. چه زمانی باید از شبیهسازیهای پیچیدهتر مانند Sliding Mesh یا FSI به جای روش MRF استفاده کنیم؟
اگر پروژه شما نیاز به بررسی نویز (آکوستیک) دارد، یا میخواهید بدانید وقتی فن از جلوی یک مانع رد میشود چه ضربهای (Impulse) به آن وارد میکند، MRF دیگر جوابگو نیست. اینجا باید آستینها را بالا بزنید و سراغ روش Sliding Mesh بروید که حل گذرا است.
همچنین اگر پرههای فن شما از جنس پلاستیک نازک است و در اثر نیروی هوا خم میشود، دیگر با یک مسئله سیالاتی خالص طرف نیستید. اینجا باید وارد دنیای [اندرکنش سیال و سازه (FSI) در فلوئنت چیست و چگونه مدلسازی میشود؟] شوید تا تغییر شکل پره و تاثیر آن روی دبی هوا را همزمان حل کنید. این نوع پروژهها معمولاً زمانبر هستند و سیستمهای قدرتمند میطلبند.
۱۵. تیم مهندسی سیمومک چگونه در پروژههای طراحی توربوماشینها صحت نتایج شبیهسازی را تضمین میکند؟
ما در سیمومک (Simumech) فقط نرمافزار ران نمیکنیم؛ ما روی فیزیک مساله وسواس داریم. وقتی صحبت از [انجام پروژه فلوئنت] در سطح صنعتی میشود، تیم ما یک پروتکل سختگیرانه را دنبال میکند تا خیال کارفرما از بابت نتایج راحت باشد:
- بررسی استقلال از مش (Grid Independence): ما پروژه را با سه سایز مش مختلف حل میکنیم تا مطمئن شویم نتایج وابسته به تعداد سلول نیست.
- اعتبارسنجی با دادههای تجربی: مقایسه نتایج با منحنیهای استاندارد (Fan Curve) یا مقالات معتبر ISI.
- چک کردن Y+: اطمینان از اینکه مدل آشفتگی با شبکهبندی لایه مرزی سازگار است.
- تحلیل دقیق عدم قطعیت: بررسی اینکه تغییرات کوچک در ورودیها چقدر خروجی را تغییر میدهد.
این رویکرد باعث شده تا بتوانیم در پروژههای حساسی مثل خنککاری تجهیزات الکترونیکی یا طراحی پرههای کمپرسور، اعتماد صنایع را جلب کنیم.
۱۶. چگونه میتوانید برای مشاوره در زمینه شبیهسازی فن و تحلیل جریانهای چرخشی با متخصصان سیمومک ارتباط بگیرید؟
اگر در حال کار روی یک پایاننامه پیچیده هستید یا یک پروژه صنعتی دارید که نیاز به دقت و تخصص بالا دارد، تیم ما آماده است تا در کنار شما باشد. چه نیاز به [انجام پروژه دانشجویی فلوئنت] داشته باشید و چه بخواهید بخشی از فرآیند تحقیق و توسعه (R&D) شرکتتان را برونسپاری کنید، میتوانید روی دانش فنی ما حساب کنید.
امیدواریم این مقاله در خصوص آموزش روش MRF (دستگاه مختصات متحرک) در فلوئنت برای شبیهسازی فن توانسته باشد گرههای ذهنی شما را باز کند و مسیر را برای شبیهسازیهای دقیقتر هموار سازد. یادتان باشد، CFD ترکیبی از علم، هنر و تجربه است؛ پس از خطا کردن نترسید و همیشه نتایجتان را با دید انتقادی بررسی کنید.