شبیه‌سازی و تحلیل جابجایی اجباری و طبیعی در عبور از یک سیلندر نیم‌دایره‌ای در ANSYS Fluent

این شبیه‌سازی به تحلیل همزمان یا جداگانه دو مکانیزم انتقال حرارت جابجایی می‌پردازد:

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

در بسیاری از موارد، این دو پدیده با یکدیگر برهم‌کنش دارند (Mixed Convection).

1. مقدمه و فیزیک مسئله

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

2. ایجاد هندسه و مش‌بندی

2.1. ایجاد هندسه (Geometry Creation):

  • نرم‌افزار: ANSYS DesignModeler یا SpaceClaim.
  • مراحل:
  1. طراحی سیلندر نیم‌دایره‌ای: یک نیم‌دایره (معمولاً در فضای 2D برای شبیه‌سازی ساده‌تر یا 3D برای مدل واقعی) با ابعاد مشخص ایجاد کنید.
  2. ایجاد دامنه سیال (Fluid Domain): یک مستطیل (در 2D) یا مکعب/استوانه (در 3D) بزرگتر در اطراف سیلندر نیم‌دایره‌ای ایجاد کنید. ابعاد دامنه باید به اندازه‌ای بزرگ باشد که اثرات مرزی بر نتایج در اطراف سیلندر تأثیرگذار نباشد (معمولاً حداقل 5 تا 10 برابر قطر سیلندر در جهت بالا دست و 10 تا 20 برابر در جهت پایین دست).
  3. Boolean Operation: از ابزار Subtract یا Cut استفاده کنید تا سیلندر نیم‌دایره‌ای از دامنه سیال کم شود.
  4. نام‌گذاری سطوح (Named Selections):
  • inlet (ورودی جریان برای جابجایی اجباری)
  • outlet (خروجی جریان)
  • cylinder_wall (سطح سیلندر نیم‌دایره‌ای)
  • top_wall, bottom_wall (دیواره‌های بالای و پایین دامنه سیال، بسته به نوع شبیه‌سازی می‌توانند wall، symmetry یا outlet باشند)
  • symmetry (اگر از تقارن استفاده می‌کنید).

2.2. مش‌بندی (Meshing):

  • نرم‌افزار: ANSYS Meshing.
  • اصول مش‌بندی:
  1. تمرکز مش در اطراف سیلندر: ناحیه اطراف سیلندر نیم‌دایره‌ای، جایی که گرادیان‌های سرعت و دما بسیار شدید هستند، نیاز به مش بسیار ریز دارد.
  2. لایه‌های مرزی (Inflation Layers): این لایه‌ها در نزدیکی سطح سیلندر ضروری هستند تا لایه‌های مرزی هیدرودینامیکی و حرارتی به دقت حل شوند. تعداد لایه‌ها، ارتفاع اولین لایه و نرخ رشد لایه‌ها باید با توجه به عدد رینولدز و مدل آشفتگی (y+) تنظیم شوند.
  3. نوع مش:
  • 2D: Quad (چهارضلعی) ترجیح داده می‌شود، به خصوص با استفاده از ابزارهای Mapped Face Meshing برای ایجاد مش‌های ساختاریافته در نواحی اطراف سیلندر.
  • 3D: Hexahedral (هگزا) در نواحی قابل کنترل و Tetrahedral/Polyhedral در نواحی پیچیده‌تر. Polyhedral اغلب تعادلی بین دقت و تعداد سلول‌ها را ارائه می‌دهد.
  1. تراکم مش (Mesh Density): به تدریج با فاصله گرفتن از سیلندر، می‌توان اندازه سلول‌ها را افزایش داد.
  2. کیفیت مش (Mesh Quality): Skewness و Orthogonal Quality را بررسی کنید. Skewness کمتر از 0.8 و Orthogonal Quality بالای 0.2 برای پایداری و دقت حل توصیه می‌شود.

3. مدل‌سازی در ANSYS Fluent

3.1. فعال‌سازی مدل‌های فیزیکی:

  • General:
  • Solver Type: Pressure-Based (برای جریان‌های با سرعت کمتر از ماخ 0.3)
  • Time: Steady (برای حالت پایدار) یا Transient (برای حالت ناپایدار، مثلاً برای پدیده‌های ریزش گرداب).
  • Models:
  • Energy Equation: فعال شود.
  • Viscous Model (Turbulence):
  • Laminar: برای اعداد رینولدز پایین (Re < 200).
  • k-epsilon (Realizable/RNG) یا k-omega (SST): برای جریان‌های آشفته. SST k-omega برای حل لایه‌های مرزی و جداشدگی جریان در اطراف اجسام بسیار مناسب است.
  • Gravity: فعال شود و جهت آن (مثلاً -9.81 m/s² در جهت Y) مشخص شود تا نیروی بویانسی عمل کند.
  • Radiation Model (اختیاری): اگر دماهای بالا یا اهمیت تشعشع بین سطوح/محیط وجود دارد، می‌توان DO یا S2S را فعال کرد.

3.2. Materials:

  • Fluid: معمولاً Air را با ویژگی‌های مناسب (چگالی، ویسکوزیته، گرمای ویژه، رسانایی حرارتی، ضریب انبساط حرارتی) تعریف کنید.
  • Density:
  • Boussinesq: برای جابجایی طبیعی با اختلاف دمای کم تا متوسط.
  • Incompressible Ideal Gas: برای اختلاف دمای زیاد یا گازهای غیر هوا.
  • User-Defined Function (UDF): برای مدل‌های چگالی پیچیده‌تر.

3.3. Boundary Conditions:

  • inlet (برای جابجایی اجباری):
  • Velocity Inlet یا Mass Flow Inlet: مقدار سرعت و دما.
  • Pressure Inlet: اگر جریان طبیعی در ورودی وجود دارد.
  • outlet:
  • Pressure Outlet: برای خروج جریان، معمولاً با فشار صفر نسبی.
  • cylinder_wall:
  • Wall: نوع No-Slip (دیواره ثابت).
  • Thermal Conditions: Temperature (دمای ثابت) یا Heat Flux (شار حرارتی ثابت).
  • Radiation: Emissivity را برای مدل تشعشع تنظیم کنید.
  • top_wall, bottom_wall:
  • می‌توانند Wall (عایق یا با دمای مشخص)، Symmetry (اگر دامنه سیال در آن جهت نامحدود فرض شود) یا Pressure Outlet (اگر این مرزها نیز برای خروج جریان باز باشند).

3.4. Solution:

  • Methods:
  • Pressure-Velocity Coupling: SIMPLE, SIMPLEC, یا Coupled. Coupled معمولاً سریع‌تر است.
  • Spatial Discretization: Second Order Upwind برای Momentom, Energy, Turbulence.
  • Controls: ضرایب آرامش (Under-Relaxation Factors) را برای پایداری حل تنظیم کنید.
  • Monitors:
  • Residuals: حداقل تا 1e-6 برای Energy و 1e-4 برای سایر معادلات.
  • Surface Monitors: نرخ انتقال حرارت روی سیلندر (Heat Transfer Rate), ضرایب درگ (Drag Coefficient) و لیفت (Lift Coefficient) روی سیلندر، دمای متوسط سیال در خروجی. این‌ها نشانه‌های بهتری از همگرایی فیزیکی هستند.

4. تحلیل نتایج (Post-Processing and Analysis)

  • میدان‌های سرعت و دما:
  • Contours: توزیع دما و سرعت در اطراف سیلندر.
  • Vectors/Streamlines: الگوهای جریان، تشکیل لایه‌های مرزی، جداشدگی جریان و تشکیل گردابه‌ها (برای جریان اجباری).
  • انتقال حرارت:
  • Local and Average Nusselt Number: عدد ناسلت محلی و میانگین روی سطح سیلندر برای کمی‌سازی انتقال حرارت جابجایی.
  • Heat Flux: شار حرارتی از سطح سیلندر.
  • نیروهای آیرودینامیکی:
  • Drag Coefficient (Cd) و Lift Coefficient (Cl): برای جریان اجباری (در صورت وجود).
  • برهم‌کنش جابجایی اجباری و طبیعی: با مقایسه نتایج حالاتی که فقط جابجایی اجباری، فقط جابجایی طبیعی، و هر دو با هم را شبیه‌سازی می‌کنید، می‌توانید سهم هر کدام را در انتقال حرارت کلی تعیین کنید.
  • اعتبارسنجی (Validation): مقایسه نتایج با داده‌های تجربی یا عددی موجود برای حالت‌های مشابه (مثلاً جریان روی سیلندر دایره‌ای) برای اطمینان از صحت مدل.

این شبیه‌سازی به دلیل وجود دو پدیده و هندسه خاص، نیاز به دقت بالایی در مراحل مش‌بندی و انتخاب مدل‌ها دارد.

کلیدواژه ها : جابجایی اجباری- Forced Convection- جابجایی طبیعی- Natural Convection- سیلندر نیم‌دایره‌ای- Semi-Circular Cylinder- ANSYS Fluent- CFD- هندسه- Geometry- دامنه سیال- Fluid Domain- Named Selections- مش‌بندی- Meshing- لایه‌های مرزی- Inflation Layers- کیفیت مش- Mesh Quality- Skewness- Orthogonal Quality- مدل آشفتگی- Turbulence Model- k-omega SST- Boussinesq Approximation- Incompressible Ideal Gas- شرایط مرزی- Boundary Conditions- Velocity Inlet- Pressure Outlet- Wall- Temperature- Heat Flux- Emissivity- ضرایب درگ- Drag Coefficient- ضرایب لیفت- Lift Coefficient- عدد ناسلت- Nusselt Number- همگرایی- Convergence- تحلیل نتایج- Post-Processing.