انتقال حرارت جابجایی یکی از مهمترین پدیدهها در بسیاری از مسائل مهندسی مانند خنککاری تجهیزات الکترونیکی، طراحی مبدلهای حرارتی، سیستمهای تهویه، توربینها و تجهیزات فرآیندی است. در بسیاری از این کاربردها، دانستن ضریب انتقال حرارت جابجایی یا همان hنقش کلیدی در تحلیل و طراحی سیستم دارد. این ضریب مشخص میکند که چه مقدار حرارت بین یک سطح و سیال اطراف آن مبادله میشود و معمولاً از رابطه معروف نیوتن برای سرمایش بهدست میآید.
در روشهای تحلیلی یا تجربی، ضریب انتقال حرارت اغلب از طریق روابط تجربی و عدد نوسلت محاسبه میشود؛ اما در شبیهسازیهای CFD مانند نرمافزار ANSYS Fluent میتوان این ضریب را به صورت دقیقتر و بر اساس میدان دما و شار حرارتی بهدست آورد. فلوئنت با حل معادلات بقای جرم، مومنتوم و انرژی، امکان بررسی توزیع دما، سرعت و شار حرارتی را فراهم میکند.
در این مقاله به بررسی روشهای بهدست آوردن ضریب انتقال حرارت جابجایی (h) در فلوئنت میپردازیم. همچنین نحوه استخراج شار حرارتی، محاسبه اختلاف دما و استفاده از دادههای خروجی نرمافزار برای تعیین مقدار h توضیح داده میشود تا بتوان از نتایج شبیهسازی برای تحلیلهای مهندسی و طراحی دقیق استفاده کرد.
مطالعه بیشتر: ترکیب انسیس با سالیدورک برای تحلیل مدلهای پیچیده؛ یک راهنمای عملی
ضریب انتقال حرارت جابجایی (h) چیست؟
ضریب انتقال حرارت جابجایی که با نماد h نمایش داده میشود، پارامتری است که شدت انتقال حرارت بین یک سطح جامد و سیال مجاور آن را توصیف میکند. این ضریب نشان میدهد که با وجود یک اختلاف دما بین سطح و سیال، چه مقدار شار حرارتی از سطح منتقل میشود. رابطه اصلی این پدیده از قانون سرمایش نیوتن بهدست میآید:
q = h (T_s – T_f)
در این رابطه، q شار حرارتی، T_s دمای سطح و T_f دمای سیال دور از دیواره است. مقدار h به عوامل مختلفی مانند سرعت سیال، خواص ترموفیزیکی سیال، نوع جریان (آرام یا آشفته)، هندسه سطح و شرایط مرزی وابسته است.
در مسائل واقعی، تعیین دقیق ضریب انتقال حرارت کار سادهای نیست. در بسیاری از منابع مهندسی، این ضریب از روابط تجربی که بر اساس آزمایشها استخراج شدهاند محاسبه میشود. این روابط معمولاً شامل اعداد بدون بعدی مانند عدد رینولدز، پرانتل و نوسلت هستند.
اما در شبیهسازیهای CFD، به جای استفاده از روابط تجربی، میدان جریان و دما بهطور کامل حل میشود. در نتیجه میتوان مقدار واقعی شار حرارتی و اختلاف دمای سطح و سیال را استخراج کرد و از آن برای محاسبه ضریب انتقال حرارت استفاده نمود. این موضوع باعث میشود که شبیهسازی عددی ابزار بسیار قدرتمندی برای تحلیل انتقال حرارت باشد.
اهمیت ضریب انتقال حرارت در تحلیلهای CFD
در تحلیلهای دینامیک سیالات محاسباتی، ضریب انتقال حرارت یکی از مهمترین پارامترهایی است که عملکرد سیستمهای حرارتی را مشخص میکند. بسیاری از تجهیزات صنعتی مانند مبدلهای حرارتی، رادیاتورها، سیستمهای خنککننده الکترونیک و حتی تجهیزات نیروگاهی به شدت به مقدار این ضریب وابسته هستند.
در شبیهسازی CFD، حل همزمان میدان سرعت و دما این امکان را فراهم میکند که مکانیزم انتقال حرارت با دقت بالاتری بررسی شود. بهعنوان مثال میتوان مشاهده کرد که در چه نواحی از سطح، انتقال حرارت بیشتر است یا کدام بخشها دارای مقاومت حرارتی بالاتری هستند.
همچنین با استفاده از نتایج شبیهسازی میتوان توزیع مکانی ضریب انتقال حرارت را روی سطح بررسی کرد. این موضوع در طراحی سیستمهای خنککاری اهمیت زیادی دارد، زیرا نشان میدهد کدام نقاط نیاز به بهبود طراحی دارند.
از سوی دیگر، مقایسه مقدار h حاصل از شبیهسازی با روابط تجربی میتواند به اعتبارسنجی مدل عددی نیز کمک کند. بنابراین استخراج ضریب انتقال حرارت از فلوئنت نهتنها برای تحلیل حرارتی بلکه برای ارزیابی صحت مدل نیز بسیار مفید است.
مراحل مدلسازی انتقال حرارت جابجایی در فلوئنت
برای محاسبه ضریب انتقال حرارت در فلوئنت ابتدا باید مسئله انتقال حرارت بهدرستی مدلسازی شود. این فرآیند شامل چند مرحله مهم مانند ایجاد هندسه، تولید شبکه، تعریف مدلهای فیزیکی و اعمال شرایط مرزی است.
در گام نخست هندسه مسئله در نرمافزارهای طراحی مانند DesignModeler یا SpaceClaim ایجاد میشود. سپس دامنه محاسباتی به شبکهای از سلولها تقسیم میشود تا معادلات حاکم روی آن حل شوند. کیفیت مش تأثیر زیادی بر دقت نتایج انتقال حرارت دارد، بنابراین انتخاب اندازه و نوع المانها بسیار مهم است.
پس از ایجاد مش، باید مدلهای فیزیکی مناسب فعال شوند. در مسائل انتقال حرارت جابجایی، معمولاً معادله انرژی فعال میشود و در صورت وجود جریان آشفته، یکی از مدلهای توربولانسی مانند k-ε یا k-ω نیز انتخاب میشود.
در مرحله بعد خواص سیال و شرایط مرزی تعریف میشوند. برای مثال دمای ورودی سیال، سرعت جریان و شرایط دمایی دیواره مشخص میشود. پس از تنظیم پارامترها، شبیهسازی اجرا شده و میدان دما و سرعت در کل دامنه محاسبه میشود.
فعالسازی معادله انرژی در فلوئنت
یکی از مهمترین مراحل در تحلیل انتقال حرارت در فلوئنت، فعالسازی معادله انرژی است. این معادله مسئول محاسبه توزیع دما در سیال و جامد بوده و بدون آن امکان بررسی انتقال حرارت وجود ندارد.
در بخش Models نرمافزار، گزینه Energy باید فعال شود. با فعال شدن این گزینه، فلوئنت معادله بقای انرژی را به مجموعه معادلات حل اضافه میکند. این معادله شامل اثرات رسانش حرارتی، جابجایی و در برخی موارد منابع حرارتی داخلی است.
پس از فعالسازی انرژی، لازم است خواص حرارتی مواد نیز تعریف شوند. پارامترهایی مانند ضریب هدایت حرارتی، ظرفیت گرمایی ویژه و چگالی در این مرحله وارد میشوند. دقت این دادهها تأثیر مستقیمی بر نتایج انتقال حرارت دارد.
همچنین در صورت وجود جریان آشفته، باید مدل توربولانسی مناسب انتخاب شود، زیرا توربولانس تأثیر قابل توجهی بر انتقال حرارت جابجایی دارد. مدلهای مختلف توربولانسی میتوانند نتایج متفاوتی برای ضریب انتقال حرارت ایجاد کنند.
روش محاسبه ضریب انتقال حرارت (h) در فلوئنت
پس از پایان شبیهسازی، میتوان از نتایج فلوئنت برای محاسبه ضریب انتقال حرارت استفاده کرد. سادهترین روش استفاده از رابطه قانون سرمایش نیوتن است:
h = q / (T_s – T_f)
در این رابطه q شار حرارتی سطح، T_s دمای دیواره و T_f دمای سیال مرجع است. هر سه این پارامترها را میتوان از نتایج شبیهسازی استخراج کرد.
در فلوئنت، شار حرارتی سطح از بخش Surface Integrals یا Reports قابل استخراج است. این مقدار نشاندهنده نرخ انتقال حرارت از سطح به سیال است. دمای سطح نیز بهراحتی از نتایج دمای دیواره قابل مشاهده است.
برای دمای سیال معمولاً از دمای میانگین جریان یا دمای ورودی سیال استفاده میشود. انتخاب این دما به نوع مسئله بستگی دارد. پس از استخراج این مقادیر، با قرار دادن آنها در رابطه بالا میتوان مقدار ضریب انتقال حرارت را محاسبه کرد.
استخراج شار حرارتی از نتایج شبیهسازی
برای محاسبه دقیق ضریب انتقال حرارت، باید مقدار شار حرارتی از سطح مورد نظر استخراج شود. این کار در فلوئنت از طریق ابزارهای گزارشگیری انجام میشود.
در بخش Reports → Surface Integrals میتوان گزینه Wall Heat Flux را انتخاب کرد. سپس سطح مورد نظر انتخاب شده و مقدار شار حرارتی متوسط یا کل محاسبه میشود. این مقدار نشان میدهد که چه میزان انرژی حرارتی از سطح عبور میکند.
در برخی مسائل، توزیع مکانی شار حرارتی نیز اهمیت دارد. در این حالت میتوان از کانتورهای حرارتی استفاده کرد تا مشخص شود کدام بخش از سطح دارای بیشترین انتقال حرارت است.
استفاده از این دادهها نهتنها برای محاسبه ضریب انتقال حرارت مفید است، بلکه به تحلیل عملکرد حرارتی سیستم نیز کمک میکند. به کمک این اطلاعات میتوان طراحی سطوح انتقال حرارت را بهینهسازی کرد.
نکات مهم برای افزایش دقت محاسبه h در فلوئنت
برای بهدست آوردن ضریب انتقال حرارت دقیق در شبیهسازی CFD باید چند نکته مهم رعایت شود. یکی از مهمترین عوامل، کیفیت شبکه نزدیک دیواره است. در مسائل انتقال حرارت جابجایی، لایه مرزی حرارتی نقش مهمی در تعیین مقدار شار حرارتی دارد.
بنابراین لازم است شبکه در نزدیکی دیواره بسیار ریز باشد تا گرادیان دما بهدرستی حل شود. استفاده از Inflation Layers در مشزنی یکی از روشهای رایج برای افزایش دقت در این ناحیه است.
انتخاب مدل توربولانسی مناسب نیز اهمیت زیادی دارد. برخی مدلها برای مسائل انتقال حرارت عملکرد بهتری دارند. برای مثال مدل k-ω SSTدر بسیاری از مسائل دیوارهای نتایج دقیقتری ارائه میدهد.
همچنین انجام Grid Independence Study یا مطالعه استقلال شبکه ضروری است. در این روش با افزایش تعداد سلولها بررسی میشود که آیا نتایج شبیهسازی تغییر میکند یا خیر. اگر نتایج پایدار شوند، میتوان اطمینان داشت که مقدار h بهدرستی محاسبه شده است.
اعتبارسنجی نتایج شبیهسازی
پس از محاسبه ضریب انتقال حرارت، بهتر است نتایج شبیهسازی با دادههای تجربی یا روابط تئوری مقایسه شوند. این کار به اعتبارسنجی مدل عددی کمک میکند و نشان میدهد که شبیهسازی تا چه حد قابل اعتماد است.
یکی از روشهای رایج، محاسبه عدد نوسلت از نتایج شبیهسازی و مقایسه آن با روابط تجربی موجود در منابع مهندسی است. اگر اختلاف نتایج در محدوده قابل قبول باشد، میتوان از مدل برای تحلیلهای بیشتر استفاده کرد.
همچنین بررسی توزیع دما و شار حرارتی در دامنه شبیهسازی کمک میکند تا از صحت فیزیکی نتایج اطمینان حاصل شود. گاهی ممکن است تنظیمات نادرست شرایط مرزی یا مدلهای فیزیکی باعث ایجاد نتایج غیرواقعی شوند.
در نهایت، ترکیب شبیهسازی CFD با دادههای تجربی بهترین روش برای تحلیل دقیق انتقال حرارت جابجایی و تعیین ضریب انتقال حرارت در سیستمهای مهندسی است.