شبیه سازی جریان همرفت آزاد لیوان آب سرد با کامسول

همرفت آزاد (Natural Convection)، که گاهی اوقات همرفت طبیعی نیز نامیده می‌شود، یکی از مکانیزم‌های اساسی انتقال حرارت در سیالات (مایعات و گازها) است. این پدیده زمانی رخ می‌دهد که اختلاف دما در یک سیال منجر به اختلاف چگالی شود. سیال گرم‌تر، به دلیل چگالی کمتر، به سمت بالا حرکت می‌کند و سیال سردتر، با چگالی بیشتر، به سمت پایین می‌نشیند. این حرکت طبیعی سیال، بدون نیاز به هیچ نیروی خارجی (مانند پمپ یا فن)، باعث انتقال انرژی حرارتی می‌شود.

شبیه‌سازی همرفت آزاد در نرم‌افزاری مانند COMSOL Multiphysics ابزاری قدرتمند برای مهندسان و محققان است. این شبیه‌سازی نه تنها به ما امکان می‌دهد تا الگوهای پیچیده جریان سیال و توزیع دما را پیش‌بینی کنیم، بلکه در بهینه‌سازی طراحی‌ها در طیف وسیعی از کاربردها، از خنک‌سازی قطعات الکترونیکی و سیستم‌های تهویه مطبوع گرفته تا سیستم‌های انرژی خورشیدی و فرآیندهای زمین‌گرمایی، نقش حیاتی دارد.

COMSOL با استفاده از یکپارچه‌سازی فیزیک‌های مختلف (مانند انتقال حرارت و دینامیک سیالات) در یک محیط واحد، ابزاری ایده‌آل برای مدل‌سازی این پدیده‌های کوپل‌شده است.

---

اصول فیزیکی همرفت آزاد در COMSOL

مدل‌سازی همرفت آزاد در COMSOL بر پایه حل همزمان معادلات حاکم بر سیالات و انتقال حرارت استوار است:

معادله بقای جرم (معادله پیوستگی)
معادلات بقای تکانه (معادلات ناویر-استوکس)
معادله بقای انرژی (معادله انتقال حرارت)

1. معادله حالت:

برای کوپل کردن معادلات فوق، چگالی سیال باید تابعی از دما باشد: $\rho =\rho \left(T\right)$ρ=ρ(T).

• تقریب بوسینسک (Boussinesq Approximation): در بسیاری از موارد همرفت آزاد که تغییرات دما نسبتاً کوچک است، از این تقریب استفاده می‌شود. در این تقریب، چگالی در تمامی جملات معادلات (به جز جمله نیروی جاذبه) ثابت فرض می‌شود، اما در جمله جاذبه، تغییرات چگالی خطی با دما در نظر گرفته می‌شود.

• چگالی کامل وابسته به دما: برای تغییرات دمای بزرگ یا سیالاتی که تقریب بوسینسک برای آن‌ها دقیق نیست، از توابع چگالی واقعی (مانند قانون گازهای ایده‌آل برای گازها) استفاده می‌شود.

---

گام‌های شبیه‌سازی همرفت آزاد در COMSOL Multiphysics

برای مدل‌سازی موفقیت‌آمیز همرفت آزاد در COMSOL، مراحل زیر را دنبال کنید:

1. تنظیمات اولیه (Model Wizard)

• Space Dimension: بسته به هندسه، 2D، 2D Axisymmetric (برای هندسه‌های استوانه‌ای) یا 3D را انتخاب کنید. 2D Axisymmetric اغلب برای شبیه‌سازی‌های همرفت آزاد با هندسه متقارن کارایی بالایی دارد.
• Add Physics: فیزیک‌های اصلی مورد نیاز:
• Heat Transfer in Fluids (ht): برای انتقال حرارت در سیال.
• Laminar Flow (spf) یا Turbulent Flow (spf): بسته به عدد ریلی (Rayleigh number) و ماهیت جریان، نوع جریان را انتخاب کنید. برای همرفت‌های ضعیف، Laminar Flow مناسب است.
• Add Multiphysics: کوپلینگ Non-Isothermal Flow را اضافه کنید. این کوپلینگ حیاتی است؛ زیرا معادلات انتقال حرارت و جریان سیال را به هم متصل می‌کند و چگالی سیال را به دما وابسته می‌سازد.

2. تعریف هندسه (Geometry)

• هندسه را با دقت رسم کنید. دامنه سیال (جایی که همرفت اتفاق می‌افتد) باید به‌طور واضح مشخص باشد. اگر جامدی نیز در مدل دخیل است (مانند دیواره‌های محفظه یا یک جسم گرم/سرد)، دامنه جامد را نیز اضافه کنید.
• اطمینان حاصل کنید که مرزهای مشترک بین دامنه‌های سیال و جامد به درستی تعریف شده‌اند.

3. تعریف مواد (Materials)

• برای سیال: خواص وابسته به دما مانند چگالی ($\rho$ρ), ویسکوزیته دینامیکی ($\mathrm{<br>}$μ), ظرفیت حرارتی ویژه ($C_p$​) و رسانایی گرمایی ($k$) را تعیین کنید. ضریب انبساط حرارتی ($\beta$) برای تقریب بوسینسک بسیار مهم است. COMSOL کتابخانه‌های مواد غنی دارد که می‌توانید از آن‌ها استفاده کنید.
• برای جامد (در صورت وجود): خواص حرارتی مانند رسانایی گرمایی و ظرفیت حرارتی ویژه را تعریف کنید.

4. تنظیمات فیزیک و شرایط مرزی (Physics Settings & Boundary Conditions)

• Heat Transfer (ht):
• Initial Values: دمای اولیه سیستم را (مثلاً دمای محیط) تنظیم کنید.
• Boundary Conditions:
• Temperature: دمای ثابت (مثلاً دمای یک دیواره گرم یا سرد).
• Heat Flux یا Heat Source: نرخ انتقال حرارت ثابت یا چشمه حرارتی.
• Heat Flux (همرفتی/تابشی): Convective Heat Flux یا Surface-to-Ambient Radiation در مرزهای خارجی برای تعامل با محیط اطراف.
• Laminar/Turbulent Flow (spf):
• Gravity: Gravity را فعال کرده و جهت آن را (مثلاً در جهت Y منفی برای همرفت رو به بالا) مشخص کنید. این مهمترین عامل محرک در همرفت آزاد است.
• Initial Values: سرعت اولیه سیال را صفر در نظر بگیرید.
• Boundary Conditions:
• No Slip: در تمام مرزهای جامد که سیال به آن چسبیده است.
• Pressure, No Viscous Stress یا Open Boundary: برای مرزهایی که سیال می‌تواند آزادانه به داخل و خارج دامنه جریان یابد (مثلاً در یک محفظه باز به محیط). این برای تنظیم فشار مرجع نیز استفاده می‌شود.
• Non-Isothermal Flow (coupling): در تنظیمات این کوپلینگ، می‌توانید انتخاب کنید که از Boussinesq approximation استفاده شود یا Full compressibility برای چگالی سیال.

5. مش‌بندی (Meshing)

• لایه‌های مرزی: برای دقت بیشتر در مناطق با گرادیان‌های دمایی و سرعتی بالا (نزدیک سطوح جامد)، از Boundary Layers استفاده کنید.
• ریز کردن مش: مش را در نواحی که انتظار تغییرات سریع دما یا سرعت را دارید، مانند مناطق نزدیک به چشمه‌های حرارتی یا سطوح، ریزتر کنید.

6. مطالعه و حل (Study & Compute)

• Study Type:
• Stationary: برای پیدا کردن حالت پایدار (steady-state) جریان و دما، در صورتی که سیستم به تعادل برسد.
• Time-Dependent: برای مشاهده تکامل پدیده در طول زمان (مثلاً خنک شدن یا گرم شدن یک سیستم). این برای بررسی پایداری و دینامیک‌های گذرا ضروری است.
• Compute: مدل را حل کنید.

7. پس‌پردازش و تحلیل نتایج (Results & Post-processing)

• میدان دما: از Contour Plot یا Surface Plot برای نمایش توزیع دما (متغیر ht.T) در دامنه‌های سیال و جامد استفاده کنید.
• میدان سرعت: Streamline Plot یا Arrow Plot برای نمایش الگوهای جریان سیال (متغیرهای spf.u, spf.v, spf.w). حلقه‌های همرفتی را در اینجا مشاهده خواهید کرد.
• فشار: Contour Plot از فشار (متغیر spf.p).
• شار حرارتی: با Surface Integration یا Line Integration می‌توانید نرخ انتقال حرارت همرفتی یا کلی را در مرزهای مشخصی محاسبه کنید.
• عدد ناسلت (Nusselt Number): برای کمی‌سازی کارایی انتقال حرارت همرفتی، می‌توان عدد ناسلت را در مرزهای تبادل حرارت محاسبه کرد.
• انیمیشن: برای مطالعات وابسته به زمان، ایجاد انیمیشن به درک بهتر دینامیک پدیده کمک می‌کند.

کلیدواژه ها : Free Convection-Natural Convection-COMSOL Multiphysics-CFD Simulation-Heat Transfer-Non-Isothermal Flow-Laminar Flow-Turbulent Flow-Boussinesq Approximation-Navier-Stokes Equations-Energy Equation-Continuity Equation-Temperature Dependent Density-Gravity-Temperature Field-Velocity Field-Heat Flux-Nusselt Number-Boundary Layers-Boundary Conditions-2D Axisymmetric Geometry-Time-Dependent Study-Stationary Study-Thermal جریان همرفت آزاد-همرفت طبیعی-کامسول Multiphysics-شبیه‌سازی CFD-انتقال حرارت-غیرایزوترمال-جریان آرام-جریان آشفته-تقریب بوسینسک-معادلات ناویر-استوکس-بقای انرژی-بقای جرم-چگالی وابسته به دما-گرانش-میدان دما-میدان سرعت-شار حرارتی-عدد ناسلت-لایه‌های مرزی-شرایط مرزی-هندسه دو بعدی متقارن-مطالعه وابسته به زمان-حالت پایدار-پایداری حرارتی-دینامیک سیالات-فیزیک حرارتی-Stability-Fluid Dynamics-Thermal Physics-Multiphysics Modeling.-