راهنمای شبیهسازی ذوب و انجماد (Solidification & Melting): آموزش کامل مدل Enthalpy-Porosity
چگونه یک انتخاب اشتباه برای “ثابت ناحیه خمیری” میتواند کل شبیهسازی انجماد شما را بیاعتبار سازد؟
شاید برایتان پیش آمده باشد که ساعتها و حتی روزها منتظر اتمام یک شبیهسازی تغییر فاز ماندهاید، اما در نهایت نتایجی گرفتید که هیچ شباهتی به واقعیت فیزیکی مسئله نداشتند. باور کنید این اتفاق بسیار رایج است. بعد از حدود ۷ سال کار تخصصی روی پروژههای CFD، یکی از تلخترین تجربیاتم مربوط به یک پروژه ریختهگری آلومینیوم بود که به خاطر انتخاب یک عدد اشتباه برای پارامتر Mushy Zone Constant، کل نتایج مربوط به سرعت پر شدن قالب و الگوی انجماد بیاعتبار شد و ما مجبور شدیم محاسبات سنگین رو از اول تکرار کنیم. این پارامتر به ظاهر ساده، میتواند تفاوت بین یک شبیهسازی دقیق و یک خروجی بیارزش باشد. تیم سیمومک در تمام مراحل انجام پروژه فلوئنت کنار شماست؛ چه برای انجام پروژه دانشجویی فلوئنت و مشاوره تخصصی انجام پایان نامه فلوئنت نیاز به راهنمایی داشته باشید ما راهکار دقیق را به شما ارائه میدهیم.
جدول راهنمای عیبیابی سریع(Troubleshooting Guide)
| مشکل رایج (Symptom) | علت احتمالی (Potential Cause) | راه حل پیشنهادی (Solution) |
| واگرایی حل در همان ایتریشنهای اول | گام زمانی (Time Step) خیلی بزرگ یا مش بیکیفیت | گام زمانی را یک مرتبه بزرگی (Order of Magnitude) کاهش دهید. کیفیت مش را (مخصوصاً Skewness) چک کنید. |
| دمای غیرفیزیکی (خیلی بالا یا پایین) | تعریف اشتباه خواص ماده یا شرایط مرزی | مقادیر گرمای نهان، ظرفیت گرمایی و شرایط مرزی حرارتی را دوباره بررسی کنید. |
| نوسانات شدید در نمودار باقیماندهها | مقدار نامناسب Under-Relaxation Factors یا Mushy Zone Constant | مقادیر Under-Relaxation را کمی کاهش دهید. مقدار Mushy Zone Constant را در محدوده 10⁴ تا 10⁷ تست کنید. |
| جبهه انجماد به صورت پلهای و غیرصاف | مش درشت در ناحیه تغییر فاز | مش را در ناحیهای که انتظار جبهه انجماد را دارید ریزتر کنید (Refine Mesh). |
در این راهنمای شبیهسازی ذوب و انجماد، ما فقط به شما نمیگوییم کدام دکمه را در فلوئنت بزنید؛ بلکه عمیقاً به چرایی هر انتخاب میپردازیم تا شما با اطمینان کامل پروژه خود را پیش ببرید. این مقاله بخشی از راهنمای جامع ما در زمینه آموزش کامل انسیس فلوئنت است که به شما کمک میکند بر تمام جنبههای این نرمافزار قدرتمند مسلط شوید.

مدل Enthalpy-Porosity واقعاً چگونه فرآیند پیچیده تغییر فاز را به یک مسئله قابل حل تبدیل میکند؟
بیایید خیلی ساده به قضیه نگاه کنیم. بزرگترین چالش در شبیهسازی ذوب یا انجماد، ردیابی مرز متحرک بین فاز جامد و مایع است. این مرز ثابت نیست و مدام جابجا میشود. مدل Enthalpy-Porosity یک ترفند هوشمندانه به کار میبرد: به جای ردیابی صریح این مرز، کل ناحیهای که در آن تغییر فاز رخ میدهد (که به آن ناحیه خمیری یا Mushy Zone میگویند) را مانند یک محیط متخلخل (Porous Media) در نظر میگیرد.
در این ناحیه، تخلخل (Porosity) از ۱ (مایع کامل) به ۰ (جامد کامل) تغییر میکند. فلوئنت با اضافه کردن یک جمله “میرایی” (Damping) به معادلات مومنتوم، سرعت سیال را در این ناحیه به تدریج به صفر نزدیک میکند. به این ترتیب، بدون نیاز به الگوریتمهای پیچیده ردیابی مرز، اثر فیزیکی انجماد به درستی شبیهسازی میشود. این روش، با وجود سادگی نسبی، برای طیف وسیعی از کاربردهای صنعتی مثل ریختهگری و جوشکاری دقت فوقالعادهای دارد و درک درست آن به فهم بهتر طرحهای گسستهسازی (Discretization Schemes) در فلوئنت نیز کمک میکند.

قبل از شروع شبیهسازی، چه خواص ترموفیزیکی کلیدی برای ماده باید به درستی تعریف شوند؟
دقت شبیهسازی شما مستقیماً به دقت دادههای ورودی بستگی دارد. قبل از اینکه حتی نرمافزار را باز کنید، مطمئن شوید این خواص را برای ماده مورد نظرتان دارید. وارد کردن یک عدد اشتباه در این بخش، یعنی زباله ورودی و زباله خروجی! 🗑️
- چگالی (Density): هم برای فاز مایع و هم جامد.
- ظرفیت گرمایی ویژه (Specific Heat): این خاصیت هم میتواند با دما تغییر کند.
- هدایت حرارتی (Thermal Conductivity): تفاوت این پارامتر بین فاز جامد و مایع معمولاً قابل توجه است.
- گرمای نهان ذوب (Latent Heat of Fusion): این قلب مسئله تغییر فاز است. عددی که انرژی لازم برای شکستن پیوندهای کریستالی را مشخص میکند.
- دمای انجماد و ذوب (Solidus & Liquidus Temperature): برای آلیاژها این دو دما متفاوت است و ناحیه خمیری را تعریف میکند.
- ویسکوزیته (Viscosity): برای فاز مایع.
نکته حرفهای: در بسیاری از مسائل واقعی، این خواص ثابت نیستند و تابعی از دما هستند. در صورت امکان، حتماً از پروفایلهای دمایی برای تعریف این خواص استفاده کنید تا نتایج به واقعیت نزدیکتر شوند، مخصوصا در کاربردهایی که با دمای بالا سروکار دارند، مثل مدلسازی تشعشع (Radiation) در دماهای بالا.
| خاصیت (برای آلومینیوم خالص) | مقدار تقریبی | واحد SI |
| دمای ذوب | 933 K | K |
| گرمای نهان ذوب | 398,000 | J/kg |
| هدایت حرارتی (مایع) | ~95 | W/m·K |
| هدایت حرارتی (جامد) | ~237 | W/m·K |
برای پیشبینی دقیق جبهه ذوب و انجماد، چگونه باید مش محاسباتی را طراحی و بهینه کنیم؟
مشبندی در مسائل تغییر فاز شوخیبردار نیست. چون گرادیانهای دما و کسر مایع در ناحیه مرز فاز بسیار شدید است، شما به مش ریز و باکیفیتی در آن نواحی نیاز دارید. اگر مش شما به اندازه کافی ریز نباشد، ممکن است جبهه انجماد را به صورت پلهای و غیرفیزیکی ببینید یا حتی کلا فرآیند به درستی ثبت نشود.
یک استراتژی خوب این است که اگر حدوداً میدانید تغییر فاز در کدام نواحی اتفاق میافتد (مثلاً نزدیک دیوارههای سرد قالب)، از تکنیکهای Body Sizing برای ریز کردن مش در آن مناطق استفاده کنید. کیفیت مش هم اهمیت زیادی دارد. معیارهایی مثل Skewness باید در محدوده قابل قبول باشند. در این زمینه، انتخاب بین مش ششوجهی (Hex) و چهاروجهی (Tet) میتواند تأثیر زیادی بر دقت و هزینه محاسباتی داشته باشد، که باید با توجه به هندسه مسئله تصمیمگیری شود.
در نرمافزار انسیس فلوئنت، چگونه میتوانیم مدل Solidification & Melting را قدم به قدم فعال و تنظیم کنیم؟
فعالسازی این مدل در فلوئنت بسیار ساده است، اما تنظیمات آن نیاز به دقت دارد. مراحل کلی به این صورت است:
- از پنل Models در فلوئنت، گزینه Solidification & Melting را پیدا کرده و آن را فعال کنید.
- به بخش Materials بروید. ماده مورد نظر خود را بسازید یا ویرایش کنید. مهمترین بخش اینجا وارد کردن مقدار Latent Heat (گرمای نهان) و دماهای Solidus و Liquidus است.
- بقیه خواص مثل چگالی، هدایت حرارتی و ظرفیت گرمایی را به دقت وارد کنید (ترجیحاً به صورت تابعی از دما).
- تمام! حالا فلوئنت میداند که باید معادلات مربوط به تغییر فاز را حل کند. این مدل در واقع یک حالت خاص از مدلسازی جریان در محیطهای متخلخل (Porous Media) است که خود فلوئنت آن را مدیریت میکند.
مهمترین پارامتر کنترلی یعنی Mushy Zone Constant را بر چه اساسی و با چه دیدی باید انتخاب کرد؟
این همان پارامتری است که در ابتدای مقاله به آن اشاره کردم. این ثابت که با نماد C نمایش داده میشود، میزان میرایی سرعت در ناحیه خمیری را کنترل میکند. به زبان ساده، مشخص میکند که سیال در ناحیه نیمهمنجمد چقدر به سختی حرکت کند.
- مقدار خیلی کم: اگر این عدد کوچک باشد (مثلاً 10³)، میرایی کافی نخواهد بود و ممکن است مقداری سرعت غیرفیزیکی در ناحیه جامد دیده شود و حل شما واگرا شود.
- مقدار خیلی زیاد: اگر این عدد بیش از حد بزرگ باشد (مثلاً 10¹⁰)، میتواند باعث ناپایداری عددی در حل معادلات شود.
یک مقدار پیشنهادی خوب برای شروع، عددی بین 10⁴ تا 10⁷ است. این عدد به نوع ماده و فیزیک مسئله بستگی دارد و گاهی نیاز است کمی آن را تغییر دهید تا به بهترین همگراعی و نتایج فیزیکی برسید. درک این پارامترها وقتی با پدیدههای پیچیدهتری مثل شبیهسازی پدیده کاویتاسیون در فلوئنت که آن هم یک نوع تغییر فاز است، روبرو میشوید اهمیت دوچندان پیدا میکند.
جدول خلاصه پارامترهای کلیدی در فلوئنت
| پارامتر در فلوئنت | توضیح مختصر | نکته کلیدی |
| Solidus Temperature | دمایی که در آن انجماد آغاز میشود. | برای مواد خالص، برابر با Liquidus Temperature است. |
| Liquidus Temperature | دمایی که در آن ذوب کامل میشود. | فاصله بین این دو، ناحیه خمیری را مشخص میکند. |
| Pure Solvent Melting Heat | گرمای نهان ذوب ماده. | این عدد، قلب مسئله تغییر فاز است. با دقت وارد شود. |
| Mushy Zone Constant | ثابت میرایی در ناحیه خمیری. | با مقادیر بین 10⁴ تا 10⁷ شروع کنید. |
چرا تنظیم صحیح معیارهای همگرایی (Residuals) در مسائل تغییر فاز تا این حد حیاتی است؟ ⚠️
در شبیهسازیهای عادی، شاید به صاف شدن نمودار باقیماندهها (Residuals) در حد 1e-3 یا 1e-4 اکتفا کنید. اما در مسائل شبیهسازی ذوب و انجماد، این کافی نیست. به دلیل ماهیت شدیداً غیرخطی مسئله (به خصوص به خاطر گرمای نهان)، نمودار باقیماندهها ممکن است نوسانات زیادی داشته باشد یا روی یک مقدار بالا ثابت بماند، در حالی که حل شما هنوز به پایداری فیزیکی نرسیده است.
به جای اینکه فقط به نمودار باقیماندهها چشم بدوزید، حتماً متغیرهای فیزیکی مهم را مانیتور کنید:
- دمای نقاط کلیدی در دامنه
- مقدار کل حرارت منتقل شده از مرزها
- کسر مایع کلی در سیستم
زمانی میتوانید بگویید حل شما همگرا شده که این مقادیر فیزیکی به یک عدد ثابت و پایدار برسند، حتی اگر نمودار باقیماندهها کمی نوسان داشته باشد. اگر با مشکلات جدی در این زمینه مواجه هستید، مطالعه مقاله ما در مورد ۷ دلیل اصلی عدم همگرایی (Divergence) در فلوئنت میتواند بسیار راهگشا باشد. از پروژههای کلاسی و انجام پروژه دانشجویی فلوئنت گرفته تا سطوح پیشرفته مثل انجام پایان نامه فلوئنت و انجام پروژه انسیس فلوئنت با هندسههای پیچیده، تیم ما آماده انجام پروژه فلوئنت با تضمین کیفیت و آموزش کامل است.

پس از اتمام حل، کانتورهای کسر مایع (Liquid Fraction) را چگونه باید تحلیل و تفسیر کرد؟
تبریک! شبیهسازی شما تمام شده. حالا وقت تحلیل نتایج است. یکی از مهمترین خروجیها در این نوع تحلیل، کانتور Liquid Fraction است. این کانتور به شما یک تصویر واضح از وضعیت فاز ماده در دامنه میدهد:
- مقدار 1 (معمولاً قرمز): نشاندهنده فاز مایع کامل است.
- مقدار 0 (معمولاً آبی): نشاندهنده فاز جامد کامل است.
- مقادیر بین 0 و 1: این ناحیه، همان Mushy Zone یا منطقه خمیری است که در آن جامد و مایع در کنار هم وجود دارند.
با بررسی این کانتور در زمانهای مختلف (برای حل گذرا)، میتوانید پیشروی جبهه انجماد یا ذوب را به وضوح ببینید، مناطقی که دچار انجماد زودهنگام میشوند را شناسایی کنید و عیوب احتمالی مثل حفرات گازی در ریختهگری را پیشبینی نمایید. این تحلیلها در حوزههایی مانند مدلسازی احتراق (Combustion) در فلوئنت که در آن هم با تغییرات فازی (شیمیایی) سروکار داریم، بسیار حیاتی هستند.
چگونه میتوانیم نتایج شبیهسازی انجماد را با دادههای آزمایشگاهی برای اطمینان از صحت اعتبارسنجی کنیم؟
هیچ شبیهسازیای بدون اعتبارسنجی (Validation) کامل نیست. حتی دقیقترین تنظیمات هم باید با دنیای واقعی مقایسه شوند. در مسائل ذوب و انجماد، یکی از بهترین روشهای اعتبارسنجی، مقایسه موقعیت جبهه انجماد در شبیهسازی با نتایج تجربی است که در مقالات معتبر علمی منتشر شدهاند. مثلاً برای مسائل کلاسیک مانند “انجماد آب در یک حفره مربعی”، دادههای آزمایشگاهی فراوانی وجود دارد.
شما میتوانید موقعیت مرز Liquid Fraction = 0.5 را در زمانهای مختلف از شبیهسازی خود استخراج کرده و آن را روی نموداری در کنار دادههای تجربی رسم کنید. اگر این دو نمودار تطابق خوبی با هم داشته باشند، میتوانید با اطمینان بالایی به نتایج خود اعتماد کنید. این فرآیند، که در مقاله راهنمای جامع اعتبارسنجی نتایج CFD به تفصیل توضیح داده شده، برای هر پروژه آکادمیک یا صنعتی یک امر ضروری است.
رایجترین دلایل واگرایی حل در شبیهسازی ذوب و انجماد و راه حلهای عملی آنها چیست؟
واگرایی (Divergence) در این نوع مسائل یک کابوس رایج است. اگر با خطای معروف “Floating Point Exception” مواجه شدید، قبل از اینکه کامپیوتر را از پنجره بیرون بیندازید، این موارد را چک کنید. اینها دلایل اصلی و راهحلهایشان هستند که اغلب در پروژههای پیچیده با آن مواجه میشویم:
- گام زمانی (Time Step) بیش از حد بزرگ: در حلهای گذرا، اگر گام زمانی بزرگ باشد، تغییرات شدید در ناحیه خمیری باعث ناپایداری حل میشود. گام زمانی را کوچکتر کنید.
- مش بیکیفیت: المانهای با Skewness بالا یا تغییر اندازه ناگهانی مش، قاتل همگرایی هستند.
- مقدار نامناسب Mushy Zone Constant: همانطور که گفتیم، این عدد نقش کلیدی دارد. کمی آن را تغییر دهید.
- تنظیمات Under-Relaxation Factors: گاهی اوقات کاهش مقادیر این فاکتورها برای متغیرهایی مثل مومنتوم و انرژی میتواند به پایداری حل کمک کند.
برای درک عمیقتر این خطاها و راهحلهای دیگر، میتوانید به راهنمای ما در مورد خطای Floating Point Exception در فلوئنت مراجعه کنید.
آیا سرعت پایین همگرایی در شبیهسازی شما را کلافه کرده است؟ راهکارهای افزایش سرعت حل کدامند؟
شبیهسازیهای تغییر فاز ذاتاً زمانبر هستند. اما چند تکنیک وجود دارد که میتواند فرآیند را سریعتر کند. یکی از موثرترین روشها، استفاده از محاسبات با عملکرد بالا (HPC) و حل موازی است. با تقسیم کردن محاسبات بین چندین هسته پردازنده (CPU)، میتوانید زمان حل را به شکل چشمگیری کاهش دهید.
یک ترفند دیگر، شروع حل با یک مش درشتتر برای رسیدن به یک حالت نیمهپایدار و سپس ادامه حل روی مش نهایی و ریزتر است. این کار کمک میکند تا حل از یک نقطه شروع بهتر آغاز شود. همچنین استفاده از طرحهای گسستهسازی مرتبه اول (First Order) در ابتدای حل و سپس سوییچ کردن به مرتبه دوم (Second Order) میتواند به پایداری و سرعت اولیه کمک کند.
چه زمانی مدل Enthalpy-Porosity دیگر پاسخگو نیست و باید به سراغ مدلهای پیشرفتهتر برویم؟
با تمام مزایایی که این مدل دارد، محدودیتهایی هم دارد. مدل Enthalpy-Porosity فرض میکند که سرعت فاز جامد و مایع در ناحیه خمیری یکسان است (که همان سرعت صفر است). این فرض برای اکثر کاربردهای رایج مثل ریختهگری در قالب ثابت کاملاً قابل قبول است.
اما در مسائلی که فاز جامد میتواند حرکت کند، مثل انجماد یک دوغاب در حال حرکت یا فرآیندهایی که در آن کریستالهای جامد توسط جریان حمل میشوند، این مدل دیگر دقت کافی را ندارد. در چنین شرایطی، باید به سراغ مدلهای پیچیدهتر چندفازی مثل مدل Eulerian-Eulerian برای جریانهای چندفازی در فلوئنت برویم که به ما اجازه میدهند برای هر فاز معادلات مومنتوم و انرژی جداگانهای را حل کنیم.
چگونه سیمومک از شبیهسازی انجماد برای کاهش عیوب در فرآیندهای ریختهگری صنعتی استفاده میکند?
در “سیمومک”، ما از این شبیهسازیها فقط برای تولید تصاویر زیبا استفاده نمیکنیم. یکی از پروژههای چالشبرانگیز ما، بهینهسازی سیستم خنککاری باتری خودروی الکتریکی با استفاده از مواد تغییر فاز دهنده (PCM) بود. با شبیهسازی دقیق فرآیند ذوب و انجماد PCM، توانستیم الگوی چیدمان کانالهای خنککننده را طوری بهینه کنیم که دمای باتری در شرایط پیک کاری، ۵ درجه سانتیگراد کاهش پیدا کند. این نتیجه مستقیماً به افزایش عمر و ایمنی باتری منجر شد. این نوع تحلیلهای عمیق، ستون فقرات خدمات انجام پروژه فلوئنت در تیم ما است.

آیا میدانستید این مدل در بهینهسازی سیستمهای ذخیره انرژی حرارتی (PCM) نیز کاربرد دارد؟
بله! یکی از جذابترین کاربردهای مدرن شبیهسازی ذوب و انجماد، در طراحی سیستمهای ذخیره انرژی حرارتی با استفاده از مواد تغییر فاز دهنده (Phase Change Materials – PCM) است. این مواد میتوانند مقدار زیادی انرژی را در دمای تقریباً ثابت ذخیره و آزاد کنند. از این تکنولوژی در سیستمهای خنککننده ساختمانها، مدیریت حرارتی باتریها و حتی لباسهای هوشمند استفاده میشود. شبیهسازی دقیق به ما کمک میکند تا بهترین ماده PCM و طراحی هندسی را برای حداکثر کردن نرخ شارژ و دشارژ انرژی پیدا کنیم.
چکلیست نهایی سیمومک قبل از اجرای یک شبیهسازی طولانیمدت ذوب و انجماد چیست؟
قبل از فشردن دکمه “Calculate”، همیشه این لیست را مرور کنید تا از هدر رفتن وقت و هزینه محاسباتی جلوگیری کنید:
- [✅] صحت خواص ماده: آیا تمام خواص ترموفیزیکی، به خصوص گرمای نهان و دماهای تغییر فاز، به درستی وارد شدهاند؟
- [✅] کیفیت مش: آیا کیفیت مش (Skewness, Orthogonal Quality) در محدوده استاندارد است؟ آیا مش در نواحی حساس به اندازه کافی ریز است؟
- [✅] شرایط مرزی: آیا شرایط مرزی حرارتی (دما، شار حرارتی) به درستی تنظیم شدهاند؟
- [✅] مقداردهی اولیه (Initialization): آیا دامنه با شرایط اولیه منطقی مقداردهی شده است؟
- [✅] تنظیمات حلگر: آیا گام زمانی (برای حل گذرا) و Under-Relaxation Factors مناسب انتخاب شدهاند؟
- [✅] مانیتورینگ: آیا مانیتور کردن متغیرهای فیزیکی کلیدی را برای بررسی همگرایی واقعی تنظیم کردهاید؟
رعایت این چکلیست ساده، تفاوت بین یک تحلیل موفق و هفتهها کار بیهوده است.
برای برونسپاری پروژههای شبیهسازی تغییر فاز و دریافت مشاوره تخصصی، چگونه با مهندسان سیمومک تماس بگیرید؟
انجام یک شبیهسازی ذوب و انجماد (Solidification & Melting) دقیق و معتبر، نیازمند ترکیبی از دانش تئوری عمیق، تسلط بر نرمافزار و تجربه عملی در مواجهه با چالشهای واقعی است. اگر با یک پروژه پیچیده صنعتی یا دانشگاهی روبرو هستید و به دنبال نتایجی قابل اعتماد در کوتاهترین زمان ممکن میگردید، تیم مهندسی سیمومک آماده است تا تخصص و تجربه خود را در اختیار شما قرار دهد.
ما خدمات زیر را به صورت تخصصی ارائه میدهیم:
- انجام کامل پروژههای شبیهسازی صنعتی و آکادمیک
- مشاوره در زمینه راهاندازی و عیبیابی شبیهسازیهای شما
- اعتبارسنجی نتایج CFD با دادههای تجربی
- آموزشهای خصوصی و گروهی متناسب با نیاز شما
برای بحث در مورد پروژه خود و دریافت راهنماییهای اولیه، میتوانید از طریق صفحه تماس با ما، با کارشناسان ما در ارتباط باشید. برای اطمینان از کیفیت و دقت نتایج، میتوانید از خدمات انجام پروژه انسیس فلوئنت ما استفاده کنید. همچنین برای پروژههای حساس، امکان عقد قرارداد و انجام پروژه فلوئنت در تهران به صورت حضوری و یا انجام پروژه فلوئنت به صورت آنلاین برای سراسر کشور فراهم است.
سوالات متداول (FAQ) به همراه پاسخ
- تفاوت اصلی بین دمای Solidus و Liquidus چیست و چرا در شبیهسازی مهم است؟
- دمای Solidus نقطهای است که در آن انجماد ماده شروع میشود (آخرین قطره مایع منجمد میشود) و دمای Liquidus نقطهای است که در آن ذوب ماده کامل میشود (آخرین ذره جامد ذوب میشود). فاصله بین این دو دما، ناحیه خمیری (Mushy Zone) را تعریف میکند. برای آلیاژها این دو دما متفاوت است اما برای مواد خالص یکسان است. وارد کردن دقیق این دو دما برای تعریف صحیح رفتار ماده در ناحیه تغییر فاز حیاتی است.
- آیا میتوان از این مدل برای شبیهسازی جوشکاری استفاده کرد؟
- بله، مدل Enthalpy-Porosity یکی از روشهای رایج برای شبیهسازی حوضچه مذاب در فرآیندهای جوشکاری است. این مدل میتواند فرآیند ذوب و انجماد سریع فلز را به خوبی شبیهسازی کند. البته برای تحلیلهای دقیقتر، ممکن است نیاز به کوپل کردن آن با مدلهای دیگر (مانند مدلهای منبع حرارتی متحرک با UDF) باشد.
- آیا برای شبیهسازی ذوب و انجماد باید از حلگر گذرا (Transient) استفاده کرد؟
- در ۹۹٪ موارد، بله. چون ذوب و انجماد فرآیندهایی وابسته به زمان هستند و ما میخواهیم پیشروی جبهه انجماد یا ذوب را در طول زمان ببینیم، استفاده از حلگر گذرا ضروری است. حل پایا (Steady-State) فقط در شرایط خاصی که به دنبال توزیع نهایی دما در یک حالت تعادلی هستیم، کاربرد دارد.
- چگونه میتوان اثر جابجایی طبیعی (Natural Convection) ناشی از تغییر چگالی را در فاز مایع مدل کرد؟
- برای این کار باید در تنظیمات ماده (Materials)، چگالی را به صورت تابعی از دما تعریف کنید (مثلاً با استفاده از تقریب بوسینسک – Boussinesq Approximation). سپس با فعال کردن گرانش (Gravity)، فلوئنت به طور خودکار نیروهای شناوری که باعث ایجاد جریان جابجایی طبیعی در حوضچه مذاب میشوند را محاسبه خواهد کرد.
- چرا نمودار باقیمانده انرژی (Energy Residual) در این شبیهسازیها اغلب بالاتر از بقیه باقیماندهها باقی میماند؟
- این یک پدیده کاملاً رایج است. به دلیل آزاد شدن یا جذب مقدار زیادی انرژی (گرمای نهان) در سلولهایی که دچار تغییر فاز میشوند، معادله انرژی به شدت غیرخطی میشود. این امر باعث میشود که باقیمانده انرژی به سختی به مقادیر خیلی پایین (مثل 1e-6) برسد. به همین دلیل است که باید همگرایی را با مانیتور کردن متغیرهای فیزیکی (مثل دما در نقاط خاص) قضاوت کنید، نه فقط با نمودار باقیماندهها.
- آیا مدل Enthalpy-Porosity انقباض یا انبساط ناشی از تغییر فاز را در نظر میگیرد؟
- خیر. این مدل به طور پیشفرض تغییر حجم ناشی از تفاوت چگالی بین فاز جامد و مایع را در نظر نمیگیرد. این یکی از محدودیتهای اصلی آن است. برای شبیهسازی دقیق پدیدههایی مثل ایجاد حفرات انقباضی (Shrinkage Cavity) در ریختهگری، نیاز به استفاده از مدلهای پیشرفتهتر چندفازی (مانند Volume of Fluid – VOF) است که تغییرات چگالی را به درستی مدل میکنند.
- مقدار Mushy Zone Constant را برای آب باید چند در نظر گرفت؟
- برای آب که یک ماده خالص است و ناحیه خمیری بسیار کوچکی دارد، معمولاً از مقادیر بزرگتری برای این ثابت استفاده میشود (مثلاً 10⁶ تا 10⁹). این کار به فلوئنت کمک میکند تا یک مرز مشخص و تیز بین آب و یخ ایجاد کند.
- آیا میتوان این مدل را با UDF ترکیب کرد؟
- بله، قطعاً. شما میتوانید با استفاده از UDF (User-Defined Function) خواص ماده را به صورت توابع پیچیدهتری از دما تعریف کنید، یا یک منبع حرارتی خاص (مثلاً برای شبیهسازی لیزر) به دامنه اضافه کنید. این کار انعطافپذیری شبیهسازی را بسیار بالا میبرد.
- حداقل کیفیت مش مورد نیاز برای این تحلیل چقدر است؟
- به عنوان یک قانون سرانگشتی، سعی کنید حداکثر Skewness مش شما زیر 0.85 و میانگین آن زیر 0.3 باشد. همچنین Orthogonal Quality باید تا حد امکان به ۱ نزدیک باشد. مهمتر از آن، داشتن حداقل ۵ تا ۱۰ لایه سلول در ناحیهای که انتظار دارید مرز فاز در آن قرار بگیرد، برای ثبت دقیق گرادیانها ضروری است.
- چرا نتایج شبیهسازی من نشان میدهد که مایع در داخل ناحیه جامد باقی مانده است (Liquid Entrapment)؟
- این پدیده میتواند فیزیکی باشد، مخصوصاً در فرآیندهای انجماد سریع آلیاژها که ممکن است مقداری مایع در بین دندریتهای جامد محبوس شود. اما اگر این اتفاق به صورت غیرفیزیکی رخ میدهد، ممکن است به دلیل مش درشت، گام زمانی بزرگ یا مقدار نامناسب Mushy Zone Constant باشد که اجازه نمیدهد حرارت به درستی از آن ناحیه خارج شود.