آیا کاهش باقیمانده‌ها (Residuals) برای همگرایی کافی است؟ راهنمای نهایی سیمومک برای قضاوت همگرایی در CFD 💡

شاید براتون پیش اومده باشه که ساعت ۲ نصف شب، خیره شدید به مانیتور و می‌بینید که خطوط نمودار باقیمانده‌ها (Residuals) توی نرم‌افزار فلوئنت یا CFX قشنگ سرازیر شدن و به زیر ۱۰ به توان ۳- رسیدن. یه نفس راحت می‌کشید، سیستم رو خاموش می‌کنید و می‌رید که بخوابید. صبح که بیدار می‌شید و کانتورها رو چک می‌کنید، می‌بینید یه جای کار میلنگه! مثلاً توزیع دما اصلا با منطق جور در نمیاد یا نیروی درگ بدست اومده پرت و پلاست.

این دقیقاً همون نقطه‌ای هست که تجربه من و شما رو از هم جدا می‌کنه. توی این ۷ سالی که درگیر پروژه‌های صنعتی بودم، از طراحی مبدل‌های حرارتی گرفته تا آیرودینامیک پهپاد، بارها دیدم که دانشجوها یا حتی مهندسای تازه‌کار گول همین نمودار سبز رنگ رو می‌خورن. آیا کاهش باقیمانده‌ها (Residuals) برای همگرایی کافی است؟ راهنمای تفسیر نمودار حل دقیقاً قراره همین “تله‌ی سبز” رو براتون باز کنه. ما توی سیمومک قرار نیست فقط بهتون بگیم دکمه Calculate رو بزنید، قراره یاد بگیریم چطور مثل یک کارآگاه رفتار حلگر رو زیر نظر بگیریم.

چرا صرفاً سبز شدن وضعیت همگرایی در نرم‌افزارهایی مثل ANSYS Fluent تضمین‌کننده صحت نتایج فیزیکی نیست؟

بیاید روراست باشیم؛ نرم‌افزار “فیزیک” رو نمی‌فهمه، اون فقط داره ریاضی حل می‌کنه. وقتی چراغ همگرایی توی فلوئنت سبز میشه، صرفاً به این معنیه که معادلات جبریِ خطی شده (Linearized Algebraic Equations) به یه تلورانس مشخصی رسیدن. اما آیا این یعنی گردابه‌های پشت سیلندر درست تشکیل شدن؟ یا انتقال حرارت توی لایه مرزی درست مدل شده؟ ابدا! ❌

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

معادلات ریاضی پشت نمودار باقیمانده‌ها واقعاً چه چیزی را درباره خطای حل عددی به ما نشان می‌دهند؟

خیلی‌ها فکر می‌کنن این خطوط رنگی توی مانیتورینگ، خطای واقعی حل نسبت به جواب مطلق هستن. کاش اینطور بود! ولی واقعیت تلخ‌تره. باقیمانده‌ها در واقع “عدم تعادل” (Imbalance) در معادلات بقا برای هر سلول هستن.
فرض کنید معادله ما به صورت

        A⋅x=bA \cdot x = bA⋅x=b

باشه. حلگر یه حدسی برای

        xxx

میزنه. باقیمانده

        RRR

میشه:

        R=∣b−A⋅x∣R = |b – A \cdot x|R=∣b−A⋅x∣

.

این یعنی چی؟ یعنی نمودار داره داد میزنه که “چقدر معادله جبری

        Ax=bAx=bAx=b

رو خوب حل کردم”. حالا اگه خود ماتریس

        AAA

(که از گسسته‌سازی مش و فیزیک مسئله اومده) از بیخ و بن ایراد داشته باشه چی؟ رزیدیوال هیچی بهت نمیگه. یه نکته ریز که خیلی‌ها بهش دقت نمیکنن اینه که این مقادیر معمولاً “اسکیل” (Scaled) میشن. یعنی اگه توی ستاپ اولیه، مقیاس‌دهی رو اشتباه زده باشید، ممکنه رزیدیوال‌ها به صورت مصنوعی پایین یا بالا نشون داده بشن. پس نمودار رزیدیوال فقط داره بهت میگه تکرارها (Iterations) دارن به یه ثبات می‌رسن، نه لزوماً به حقیقت.

آیا می‌دانستید تفاوت “همگرایی ریاضی” با “همگرایی فیزیکی” می‌تواند کل اعتبار پروژه پایانی شما را زیر سوال ببرد؟

اینجا دقیقاً جاییه که تفاوت یه اپراتور نرم‌افزار با یه مهندس مکانیک مشخص میشه.

• همگرایی ریاضی: اعداد توی ماتریس‌ها دیگه تغییر نمی‌کنن (رزیدیوال‌ها صاف میشن).
• همگرایی فیزیکی: پارامترهای مهندسی مسئله (مثل دبی خروجی، افت فشار، نیروی لیفت) ثابت شدن.

ممکنه رزیدیوال شما روی

        10−410^{-4}10−4

ثابت شده باشه، ولی نیروی درگ روی ایرفویل هنوز داره نوسان می‌کنه یا با یه شیب ملایم داره تغییر می‌کنه. اگر اینجا حل رو متوقف کنید، فاجعه رخ میده. توی یکی از پروژه‌های توربوماشین که انجام می‌دادیم، رزیدیوال‌ها فلت شده بودن ولی راندمان کمپرسور هنوز داشت هر ۵۰۰ تکرار، یک درصد بالا می‌رفت. اگه همونجا قطع می‌کردیم، طراحی کلاً رد می‌شد. پس همیشه یادتون باشه، همگرایی ریاضی شرط لازمه، ولی کافی نیست.

چرا در پروژه‌های پیچیده مهندسی مکانیک نباید به تنظیمات پیش‌فرض همگرایی (۱۰ به توان ۳-) اعتماد کنیم؟

این

        10−310^{-3}10−3

که دیفالت اکثر نرم‌افزاراست، بیشتر شبیه یه شوخی میمونه برای کارهای حساس! 📉 برای یه جریان ساده لمینار توی لوله شاید بد نباشه، اما برای جریان‌های پیچیده، آشفته، یا احتراقی، این مقدار یعنی “هنوز اول راهیم”.

تجربه نشون داده برای مقاله‌های ISI و پروژه‌های صنعتی سنگین، ما معمولاً نیاز داریم باقیمانده‌ها رو تا

        10−510^{-5}10−5

یا حتی

        10−610^{-6}10−6

پایین بیاریم. مخصوصاً معادله پیوستگی (Continuity) و انرژی که خیلی حساسن. وقتی رزیدیوال روی

        10−310^{-3}10−3

هست، یعنی هنوز خطای موضعی قابل توجهی توی دامین داریم. یه بار توی شبیه‌سازی نازل سوخت، با

        10−310^{-3}10−3

همه چی آروم بود، ولی وقتی ادامش دادیم تا

        10−610^{-6}10−6

، تازه جدایش جریان (Flow Separation) توی دهانه نازل خودشو نشون داد. پس لطفاً تنبلی نکنید و اون معیار Convergence Criterion رو دستکاری کنید.

چگونه تعریف مانیتورهای سطحی مثل ضریب درگ یا دمای خروجی به قضاوت صحیح ما از پایداری حل کمک می‌کند؟

این بخش مهم‌ترین توصیه منه به شما: مانیتور بسازید!
بدون مانیتور کردن کمیت‌های فیزیکی، شما دارید توی تاریکی رانندگی می‌کنید. قبل از اینکه دکمه Run رو بزنید، باید از خودتون بپرسید: “هدف من از این شبیه‌سازی چیه؟”

• اگه دنبال نیروی درگی: مانیتور
•         CdC_dCd​

بساز.

• اگه مبدل حرارتیه: مانیتور دمای میانگین خروجی (Mass-weighted Average Temperature) بساز.

نرم‌افزار بهت اجازه میده اینا رو روی نمودار ببینی. قانون طلایی اینه: حل زمانی تمام شده که این مانیتورها کاملاً افقی و بی تغییر بشن. حتی اگه رزیدیوال‌ها نوسان داشته باشن، اگه مانیتورهای فیزیکی ثابت باشن، یعنی به احتمال ۹۹٪ به جواب رسیدید. برای اینکه بدونید دقیقاً چه پارامترهایی رو باید مانیتور کنید و چطور گزارش‌گیری کنید، مقاله استخراج داده‌های کمی دقیق (مانند ضرایب درگ و لیفت) راهنمای فوق‌العاده‌ایه که نباید از دستش بدید.

آیا بالانس جرم و انرژی در ورودی و خروجی دامین محاسباتی شما به طور کامل برقرار شده است؟

بعد از اینکه حل تموم شد، نرید سریع سراغ کانتور رنگی گرفتن! اول برید سراغ منوی Reports > Fluxes.
چک کنید ببینید آیا “جرم ورودی” با “جرم خروجی” برابره؟ (Net Mass Flow Rate). این عدد باید خیلی خیلی نزدیک به صفر باشه.
توی شبیه‌سازی‌های حرارتی هم همینطور؛ کل حرارت تولیدی یا ورودی باید با حرارت خروجی بالانس باشه.

اگه دیدید Net Imbalance شما مثلاً ۵٪ کل جریان ورودیه، یعنی ۵٪ از جرم سیال توی محاسبات غیب شده! 🤯 این یعنی حل شما هنوز همگرا نشده، حتی اگه نمودار رزیدیوال کف زمین باشه. معمولاً این مشکل وقتی پیش میاد که مش بی‌کیفیت باشه یا “تعداد تکرار” کافی نبوده. این چک‌لیستِ منه برای هر پروژه قبل از تحویل:

1. بررسی رزیدیوال‌ها (زیر حد مجاز؟)
2. بررسی مانیتورها (ثابت شدن؟)
3. بررسی بالانس جرم و انرژی (خطای زیر ۱٪؟)

چرا نوسانی شدن نمودار باقیمانده‌ها همیشه نشانه بدی نیست و چه ارتباطی با فیزیک جریان‌های ناپایا دارد؟

بعضی وقتا هر کاری می‌کنید نمودار صاف نمیشه و شروع میکنه به سینوسی شدن یا زیگزاگ رفتن. دانشجوها اینجا معمولاً می‌ترسن و فکر می‌کنن مدل غلطه. اما صبر کنید! شاید فیزیک مسئله شما ذاتاً ناپایا (Unsteady) هست.
مثلاً جریان پشت یک استوانه (Vortex Shedding) ذاتاً نوسانیه. هرچقدر هم سعی کنید به زورِ حلگر Steady (پایا) این رو حل کنید، بازم رزیدیوال‌ها نوسان می‌کنن چون گردابه‌ها دارن کنده میشن.

اینجا نوسان رزیدیوال یعنی حلگر داره سعی میکنه یه پدیده وابسته به زمان رو توی یه فریم ثابت گیر بندازه و نمیتونه. اگر نوسانات منظم و پریودیک دیدید، تبریک میگم! شما احتمالاً با یک جریان گذرا طرف هستید و باید حلگر رو به Transient تغییر بدید. البته گاهی هم نوسان بخاطر مش بد یا تایم‌استپ نامناسبه که باید جداگانه بررسی بشه، ولی همیشه نوسان مساوی با خطا نیست.

چگونه با قرار دادن پروب‌های نقطه‌ای در نواحی پرگردابه از عدم تغییرات محلی جریان مطمئن شویم؟

مانیتورهای سطحی (که بالا گفتم) میانگین می‌گیرن. گاهی میانگین ثابته، ولی یه گوشه‌ای از دامین جریان داره دیوونه‌بازی درمیاره! برای پروژه‌های خیلی حساس، ما توی سیمومک از “Point Monitors” استفاده می‌کنیم.
یعنی چی؟ یعنی مختصات چند تا نقطه حساس (مثلاً پشت مانع، داخل ناحیه ویک، یا نزدیک دیواره) رو میدیم و سرعت یا فشار اون نقطه خاص رو مانیتور می‌کنیم.

این روش مثل اینه که یه سنسور بذارید وسط جریان. اگه نمودار سرعت اون نقطه بعد از ۲۰۰۰ تکرار هنوز داشت بالا پایین می‌شد، یعنی جریان هنوز اونجا شکل نگرفته. این تکنیک برای وقتهایی که شک دارید “آیا گردابه‌های کوچک هم همگرا شدن یا نه” عالیه و دید خیلی عمیق‌تری نسبت به کل دامین بهتون میده. 📍

آیا ثابت ماندن (Stall) نمودار باقیمانده‌ها در مقادیر بالا لزوماً به معنای شکست شبیه‌سازی و نیاز به مش‌بندی مجدد است؟

این یکی از اون صحنه‌های اعصاب‌خردکنه: نمودار باقیمانده‌ها میاد پایین، میاد پایین، و یهو روی عدد

        10−210^{-2}10−2

یا

        10−310^{-3}10−3

گیر می‌کنه و تبدیل میشه به یه خط صاف. نه بالا میره، نه پایین. اینجور وقتا اکثر بچه‌ها سریع دکمه Stop رو می‌زنن و میرن سراغ تغییر مش. اما دست نگه دارید! ✋

تجربه من میگه همیشه مشکل از مش نیست. گاهی وقتا حلگر توی یه “مینیمم محلی” گیر کرده. یه ترفند که ما توی سیمومک استفاده می‌کنیم اینه که برای ۵۰ تا ۱۰۰ تکرار، نوع حلگر رو تغییر میدیم (مثلاً از coupled به SIMPLE) یا یه تغییر کوچیک توی ضرایب زیر-تغییر (Relaxation Factors) میدیم تا به اصطلاح یه “شوک” به حل وارد بشه. خیلی وقتا با همین کار نمودار دوباره شروع میکنه به ریزش. البته اگر بعد از هر ترفندی تکون نخورد، اون وقت باید شک کنید که شاید فیزیک مسئله رو اشتباه تعریف کردید یا کیفیت مش توی نواحی حساس افتضاحه.

چگونه کیفیت مش و پارامترهایی مثل Skewness می‌توانند باعث ایجاد همگرایی کاذب یا واگرایی در حلگر شوند؟

ببینید، حلگرهای CFD هرچقدر هم قدرتمند باشن، وقتی به یه المان (Cell) کج و کوله با Skewness بالای 0.95 برسن، قاطی می‌کنن. وقتی المان بیش از حد اریبه، معادلات گسسته‌سازی شده دقتشون رو از دست میدن. نتیجه؟ ممکنه نمودار باقیمانده‌ها واگرا بشه (Error معروف Floating Point Exception) یا بدتر از اون: همگرایی کاذب.

همگرایی کاذب یعنی نمودار میاد پایین و سبز میشه، ولی جوابی که میده کاملاً پرته. انگار که با یه خط‌کش کج بخواید دیوار صاف بسازید؛ دیوار ساخته میشه ولی کجه! 🧱 پس قبل از اینکه وقتتون رو روی تنظیمات حلگر بذارید، توی مرحله مش‌زنی (Meshing) مطمئن بشید که حداقل کیفیت متعامد (Orthogonal Quality) بالای 0.1 هست. یه بار روی یه پروژه منیفولد دود سه روز وقت گذاشتیم تا حل همگرا بشه، آخرش فهمیدیم سه تا سلولِ هرمیِ خراب گوشه مدل باعث تمام بدبختی‌ها بود!

تغییر طرحواره‌های گسسته‌سازی از مرتبه اول به دوم چه تاثیری بر رفتار نمودار باقیمانده‌ها و دقت نهایی دارد؟

این نکته خیلی کلیدیه. وقتی تنظیمات رو روی First Order Upwind می‌ذارید، حلگر خیلی سریع و راحت همگرا میشه. نمودارها مثل سرسره میان پایین. چرا؟ چون این روش ذاتاً “پخش‌کننده” (Diffusive) است و نوسانات رو توی خودش حل میکنه و از بین میبره. اما دقتش چی؟ پایینه.

برای کارهای دقیق مهندسی، حتماً باید برید سراغ Second Order. وقتی این سوییچ رو انجام میدید، یهو می‌بینید نمودار باقیمانده‌ها می‌پره بالا! 📈 نترسید، این نشونه خوبیه. یعنی حلگر تازه داره متوجه گرادیان‌های واقعی و جزئیات جریان میشه. پیشنهاد من؟ همیشه ۱۰۰۰ تکرار اول رو با First Order برید تا جریان شکل بگیره، بعد سوییچ کنید روی Second Order برای دقت نهایی. اگر دیدید روی مرتبه دوم واگرا شد، یعنی کیفیت مشتون برای حل دقیق مناسب نیست.

در شبیه‌سازی‌های وابسته به زمان (Transient) چه معیاری برای تعیین تعداد تکرار در هر گام زمانی مناسب است؟

توی حل‌های زمانی، نمودار باقیمانده‌ها رفتار دندان‌اره‌ای داره. توی هر گام زمانی (Time Step) باید بیاد پایین و دوباره برای گام بعدی میره بالا. سوال اینه: توی هر گام چقدر باید صبر کنیم؟
قانون نانوشته اینه: باقیمانده‌ها باید توی هر گام زمانی حداقل ۲ تا ۳ مرتبه بزرگی (Order of Magnitude) کاهش پیدا کنن.

اگه می‌بینید توی ۵ تکرار همگرا میشه، احتمالاً گام زمانیتون خیلی کوچیکه و دارید وقت تلف می‌کنید. اگه ۵۰ تکرار میگذره و هنوز پایین نیومده، گام زمانیتون بزرگه و دقت حل داره قربانی میشه. معمولاً تنظیم “۲۰ تکرار در هر گام زمانی” یه نقطه شروع خوبه، ولی وحی منزل نیست. همه اینا به عدد کورانت(Courant Number) بستگی داره که باید حواستون باشه (معمولاً زیر ۱ برای حل‌های دقیق).

آیا بررسی کانتورهای سرعت و فشار در حین حل می‌تواند خطاهایی را نشان دهد که در نمودار باقیمانده‌ها پنهان مانده‌اند؟

گاهی وقتا باقیمانده‌ها عالی‌ان، ولی وقتی کانتور فشار رو می‌بینید، انگار صفحه شطرنجیه! (Checkerboarding). یعنی یه سلول فشارش بالاست، سلول بغلی پایین. این یه خطای عددی معروفه که نمودار رزیدیوال اصلاً نشونش نمیده.

برای همین همیشه میگم “چشمی” هم کار رو چک کنید. کانتورها باید نرم و پیوسته باشن (مگر در محل شاک یا شوک). اگر تغییرات ناگهانی و غیرفیزیکی دیدید، به حل شک کنید. اینجا بحث Validaiton خیلی مهم میشه. برای اینکه بدونید چطور نتایج رو با واقعیت تطبیق بدید، حتماً مقاله جامع ما در مورد چگونه نتایج شبیه‌سازی را در فلوئنت اعتبارسنجی(Validation) کنیم؟ رو بخونید. اونجا کامل توضیح دادیم که چطور مچِ نرم‌افزار رو بگیرید.

چگونه استفاده از ضرایب زیر-تغییر (Under-Relaxation Factors) می‌تواند به کنترل نوسانات شدید باقیمانده‌ها کمک کند؟

ضرایب URF مثل پدال ترمز و گاز ماشین می‌مونن. 🚗 وقتی جاده لغزنده‌ست (جریان پیچیده یا مش بد)، نباید گاز رو تا ته فشار بدید. مقادیر پیش‌فرض (مثلاً 0.3 برای فشار و 0.7 برای مومنتوم) معمولاً خوبن، ولی اگه نمودارها نوسان شدید داشتن، این اعداد رو کم کنید (مثلاً فشار 0.2 و مومنتوم 0.5).

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

تیم فنی سیمومک در پروژه‌های صنعتی حساس از چه چک‌لیست نهایی برای تایید اعتبار همگرایی استفاده می‌کند؟

ما در سیمومک فقط کار رو انجام نمیدیم که تحویل بدیم؛ ما کیفیت رو “مهندسی” می‌کنیم. خدمات ما طیف وسیعی از نیازهای مهندسی مکانیک رو پوشش میده، اما برای هر کدوم، یه پروتکل سفت و سخت داریم.
خدمات تخصصی ما شامل موارد زیر هست:

• شبیه‌سازی سیالاتی (CFD): تحلیل‌های چندفازی، احتراق، و آیرودینامیک.
• تحلیل‌های جامداتی (FEA): تحلیل تنش، خستگی و دینامیک.
• کدنویسی اختصاصی (UDF/Scripting): برای مدل کردن فیزیک‌های خاص که توی نرم‌افزار نیست.
• مشاوره و بهینه‌سازی: کاهش هزینه ساخت با شبیه‌سازی پیش از تولید.

توی تمام این پروژه‌ها، ما تا وقتی “تست استقلال از مش” (Grid Independence Test) و “تست استقلال از دامنه” رو پاس نکنیم، خروجی رو تایید نمی‌کنیم. این یعنی مشتری مطمئن میشه نتایجش به خاطر شانس و اقبال نبوده!

آیا کاهش باقیمانده‌ها تا ۱۰ به توان ۶- برای همه پروژه‌های مهندسی ضروری است یا باعث اتلاف هزینه محاسباتی می‌شود؟

مدیرای پروژه همیشه عجله دارن. “مهندس چی شد؟” جمله‌ایه که زیاد می‌شنویم. واقعیت اینه که برای یه پروژه دانشجویی یا یه تخمین اولیه صنعتی، شاید

        10−310^{-3}10−3

هم کار راه انداز باشه و نیاز نباشه سه روز سیستم روشن بمونه تا به

        10−610^{-6}10−6

برسه.

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

چرا وجود جریان‌های برگشتی (Reversed Flow) در خروجی می‌تواند نمودارهای همگرایی را فریب دهد و راه حل آن چیست؟

دیدید توی کنسول فلوئنت یهو مینویسه Reversed flow in xxx faces on pressure-outlet؟ ⚠️ خیلیا اینو نادیده می‌گیرن، ولی این سمه!
این یعنی دامین محاسباتی شما کوتاه بوده و گردابه‌ها رسیدن به خروجی و دارن برمی‌گردن داخل. این باعث میشه همگرایی به شدت سخت بشه و حتی اگه رزیدیوال بیاد پایین، جواب غلطه چون شرط مرزی شما نقض شده. راه حل ساده‌ست: طول لوله یا دامین خروجی رو بیشتر کنید تا جریان قبل از خروج، آرام و توسعه‌یافته بشه.

در نهایت، یادتون باشه که نمودار، فقط یه ابزاره. این شما هستید که به عنوان مهندس باید قضاوت کنید. امیدوارم دفعه بعد که اون خطوط رنگی رو دیدید، با نگاهی انتقادی‌تر بررسیشون کنید. و سوال اصلی رو همیشه گوشه ذهنتون داشته باشید:  آیا کاهش باقیمانده‌ها (Residuals) برای همگرایی کافی است؟ راهنمای تفسیر نمودار حل به ما یاد داد که جواب قطعاً “خیر” است؛ باید عمیق‌تر نگاه کرد.