سلول‌های بنارد (Bénard Cells) الگوهای جریان همرفتی منظم و زیبایی هستند که در یک سیال محبوس بین دو صفحه افقی، هنگامی که صفحه پایینی داغ‌تر از صفحه بالایی باشد، شکل می‌گیرند. این پدیده کلاسیک همرفت طبیعی، نمونه‌ای بارز از بی‌ثباتی حرارتی (Thermal Instability) است و در سیالات مختلف از جمله هوا، آب و روغن قابل مشاهده است. شبیه‌سازی این پدیده‌ها در COMSOL Multiphysics به درک بهتر مکانیسم‌های انتقال حرارت و دینامیک سیالات کمک می‌کند.

1. تعریف هندسه و مواد

گام اول در شبیه‌سازی، ایجاد هندسه است. این هندسه شامل یک فضای مستطیلی یا مکعبی (بسته به شبیه‌سازی دو بعدی یا سه بعدی) بین دو صفحه موازی می‌شود. ضخامت سیال (فاصله بین دو صفحه) و ابعاد افقی دامنه نقش مهمی در شکل‌گیری و اندازه سلول‌ها دارند. سپس، ماده سیال را به این دامنه اختصاص می‌دهیم. خواص سیال مانند چگالی، ویسکوزیته، ظرفیت حرارتی و هدایت حرارتی باید به درستی تعریف شوند. برای مدل‌سازی همرفت طبیعی، لازم است تغییرات چگالی سیال با دما در نظر گرفته شود. معمولاً از تقریب بوسینسک (Boussinesq Approximation) استفاده می‌شود که در آن چگالی سیال به جز در ترم گرانش، ثابت فرض شده و در ترم گرانش تابعی خطی از دما در نظر گرفته می‌شود. ضریب انبساط حرارتی (Coefficient of Thermal Expansion) نیز یک پارامتر کلیدی در این مدل‌سازی است.

2. انتخاب فیزیک‌ها

برای شبیه‌سازی سلول‌های بنارد، دو فیزیک اصلی مورد نیاز است:

  • جریان سیال (Fluid Flow): معمولاً فیزیک “Laminar Flow” (جریان آرام) انتخاب می‌شود، زیرا سلول‌های بنارد در اعداد رایلی (Rayleigh Number) پایین تا متوسط که متناظر با جریان آرام است، شکل می‌گیرند. این فیزیک معادلات ناویه-استوکس (Navier-Stokes Equations) را حل می‌کند.
  • انتقال حرارت (Heat Transfer): فیزیک “Heat Transfer in Fluids” برای مدل‌سازی توزیع دما در سیال و انتقال حرارت ناشی از هدایت و همرفت استفاده می‌شود.

این دو فیزیک باید با یکدیگر کوپل (Couple) شوند. در COMSOL، معمولاً از کوپلینگ “Non-Isothermal Flow” (جریان غیرهم‌دما) استفاده می‌شود که به طور خودکار اثرات تغییرات دما بر چگالی سیال و در نتیجه نیروی شناوری (Buoyancy Force) را در معادلات جریان لحاظ می‌کند.

3. شرایط مرزی

تعریف دقیق شرایط مرزی برای تحریک و توسعه سلول‌های بنارد ضروری است:

  • مرزهای دما (Temperature Boundaries):
  • صفحه پایینی: دمای ثابت و بالاتر (Thot).
  • صفحه بالایی: دمای ثابت و پایین‌تر (Tcold).
  • مرزهای جریان سیال (Fluid Flow Boundaries):
  • دیواره‌ها (صفحات بالا و پایین و دیواره‌های جانبی): شرط “No-Slip” (عدم لغزش) اعمال می‌شود، به این معنی که سرعت سیال در مجاورت دیواره‌ها صفر است.
  • برای مدل‌سازی یک دامنه بی‌نهایت در جهت افقی، می‌توان از شرایط مرزی “Periodic” (تناوبی) در دیواره‌های جانبی استفاده کرد، که به سیال اجازه می‌دهد به طور طبیعی الگوی سلول‌های همرفتی را توسعه دهد. در غیر این صورت، دیواره‌های جامد کناری نیز No-Slip خواهند بود.

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

کیفیت مش‌بندی تأثیر زیادی بر دقت نتایج دارد. مش باید به اندازه کافی ریز باشد تا گرادیان‌های دما و سرعت را به خوبی در سراسر سیال، به ویژه در نزدیکی دیواره‌ها، ثبت کند. در همرفت طبیعی، لایه‌های مرزی حرارتی و سرعت در نزدیکی دیواره‌ها می‌توانند بسیار نازک باشند، بنابراین استفاده از مش‌های متراکم‌تر یا لایه‌های مرزی (Boundary Layers) در این مناطق توصیه می‌شود.

5. نوع مطالعه (Study Type)

برای مشاهده شکل‌گیری و تکامل سلول‌های بنارد، معمولاً از مطالعه “Time-Dependent” (وابسته به زمان) استفاده می‌شود. این به شما امکان می‌دهد تا دینامیک توسعه الگوهای همرفتی را از حالت اولیه تا رسیدن به یک حالت پایدار یا شبه‌پایدار مشاهده کنید.

  • شرایط اولیه (Initial Conditions): برای آغاز همرفت، معمولاً لازم است یک اختلاف دمایی اولیه بین صفحات وجود داشته باشد. برای شکستن تقارن و تحریک شکل‌گیری الگوها، می‌توان یک اغتشاش کوچک (Small Perturbation) در میدان دما یا سرعت در شرایط اولیه اعمال کرد. این اغتشاش می‌تواند یک نویز تصادفی یا یک گرادیان دمایی بسیار کوچک باشد.

6. پارامترهای کلیدی و نتایج مورد انتظار

پارامتر بی‌بعد اصلی که رفتار سلول‌های بنارد را کنترل می‌کند، عدد رایلی (Rayleigh Number) است. این عدد نسبت نیروهای شناوری به نیروهای اتلاف‌کننده (ویسکوزیته و نفوذ حرارتی) را نشان می‌دهد. با افزایش عدد رایلی، انتقال از هدایت خالص به همرفت و سپس به شکل‌گیری الگوهای سلولی منظم اتفاق می‌افتد.

نتایج مورد انتظار شامل:

  • نقشه‌های توزیع دما (Temperature Distribution): نشان‌دهنده نواحی گرم و سرد و نحوه انتقال حرارت.
  • نقشه‌های میدان سرعت (Velocity Field): نمایش جریان‌های صعودی (Upwelling) و نزولی (Downwelling) سیال که سلول‌های همرفتی را تشکیل می‌دهند.
  • بردارهای سرعت (Velocity Vectors): الگوهای گردابی شکل (Vortex Patterns) که مشخصه سلول‌های بنارد هستند.

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

کلیدواژه ها : COMSOL Multiphysics-سلول‌های بنارد-Bénard Cells-همرفت طبیعی-Natural Convection-بی‌ثباتی حرارتی-Thermal Instability-فضای بین دو صفحه-تقریب بوسینسک-Boussinesq Approximation-خواص سیال-Laminar Flow-Heat Transfer in Fluids-Non-Isothermal Flow-شرایط مرزی-No-Slip-Periodic-مش‌بندی-Boundary Layers-Time-Dependent-اغتشاش اولیه-Initial Perturbation-عدد رایلی-Rayleigh Number-توزیع دما-میدان سرعت-بردارهای سرعت