1. تنظیمات شبیهسازی در انسیس فلوئنت (Fluent Setup)
این مرحله شامل وارد کردن تمام اطلاعات فیزیکی و عددی لازم برای حل معادلات جریان است.
1.1. تنظیمات کلی حلگر (General Solver Settings)
- نوع حلگر (Solver Type): برای جریانهای دوفازی با سیالات تراکمناپذیر مانند آب و روغن، معمولاً از حلگر مبتنی بر فشار (Pressure-Based Solver) استفاده میشود.
- حالت حلگر (Time-Dependent/Steady):
- اگر هدف بررسی وضعیت نهایی پایدار (Steady-State) باشد که در آن متغیرها با زمان تغییر نمیکنند، میتوان از حلگر حالت پایدار استفاده کرد.
- اما برای ردیابی دقیق حرکت رابط بین آب و روغن، که اغلب یک پدیده وابسته به زمان است (مثلاً نحوه پر شدن مخزن، یا حرکت یک حباب)، استفاده از حلگر وابسته به زمان (Transient Solver) ضروری است. این حالت امکان مشاهده دینامیک جریان را فراهم میکند.
1.2. انتخاب و تنظیم مدلهای فیزیکی (Physics Model Selection)
- مدل چندفازی (Multiphase Model): همانطور که قبلاً ذکر شد، برای جریان دوفازی آب و روغن با رابط مشخص، مدل Volume of Fluid (VOF) بهترین گزینه است.
- در تنظیمات VOF، باید فازهای اولیه و ثانویه را تعریف کنید (مثلاً آب فاز اولیه و روغن فاز ثانویه).
- کشش سطحی (Surface Tension): این پارامتر حیاتی است. در این بخش، مقدار عددی کشش سطحی بین آب و روغن باید وارد شود. فلوئنت از یک مدل خاص برای محاسبه کشش سطحی (Continuum Surface Force - CSF) استفاده میکند که اثرات آن را بر حرکت رابط اعمال میکند.
- مدل فعالسازی رابط (Interface Tracking Schemes): برای حفظ وضوح رابط، فلوئنت روشهای دیسکرتیزاسیون خاصی مانند Geo-Reconstruct، CICSAM یا HRIC را برای کسری حجم ارائه میدهد. انتخاب این روشها برای جلوگیری از “پخش شدن” رابط و حفظ وضوح آن بسیار مهم است.
- خواص مواد (Material Properties): برای هر فاز (آب و روغن) باید خواص فیزیکی دقیق وارد شود:
- چگالی (Density): بسته به نوع شبیهسازی (ایزوترمال یا همراه با انتقال حرارت)، میتواند ثابت یا وابسته به دما/فشار باشد.
- ویسکوزیته (Viscosity): ثابت یا وابسته به دما.
- خواص حرارتی (در صورت نیاز): گرمای ویژه، رسانایی حرارتی، و ضریب انبساط حرارتی.
- مدل تلاطم (Turbulence Model): اگر جریان در دامنه آشفته باشد، انتخاب یک مدل تلاطم مناسب (مانند k-epsilon، k-omega، SST) و تنظیمات مربوط به آن ضروری است. مدل SST اغلب برای جریانهای شامل جدایش و گرادیانهای فشار نامطلوب عملکرد خوبی دارد.
1.3. شرایط مرزی (Boundary Conditions)
این مرحله شامل تعریف دقیق رفتار جریان در مرزهای دامنه محاسباتی است:
- ورودیها (Inlets):
- نوع ورودی: سرعت (Velocity Inlet) یا دبی جرمی (Mass Flow Inlet).
- مقدار فازها: برای هر ورودی، باید کسری حجم فازهای وارد شده (مثلاً در ورودی آب، کسری حجم آب 1 و روغن 0 است) مشخص شود.
- خروجیها (Outlets):
- نوع خروجی: فشار خروجی (Pressure Outlet) رایجترین گزینه است.
- دیوارهها (Walls):
- شرایط عدم لغزش (No-Slip): این شرط استاندارد برای سیالات چسبناک است که سرعت سیال در سطح دیواره برابر با سرعت دیواره (معمولاً صفر) است.
- زاویه تماس (Contact Angle): برای مدل VOF، تعریف زاویه تماس بین رابط آب و روغن و سطح دیوارهها (بسته به خاصیت ترشوندگی) حیاتی است. این زاویه بر شکلگیری و حرکت رابط در نزدیکی دیواره تأثیر میگذارد.
1.4. تنظیمات روشهای حل (Solution Methods)
- کوپلینگ فشار-سرعت (Pressure-Velocity Coupling): الگوریتمهایی مانند PISO، SIMPLE یا COUPLED برای اتصال معادلات مومنتوم و فشار استفاده میشوند. PISO اغلب برای حلهای Transient و COUPLED برای همگرایی سریعتر مناسب هستند.
- دیسکرتیزاسیون فضایی (Spatial Discretization): انتخاب روشهای دیسکرتیزاسیون برای هر متغیر:
- برای مومنتوم، انرژی و تلاطم، معمولاً از روش مرتبه دوم (Second Order Upwind) استفاده میشود.
- برای کسری حجم (Volume Fraction)، به دلیل اهمیت حفظ وضوح رابط، استفاده از روشهای مرتبه بالاتر و با قابلیت فشردگی (Compressive Schemes) مانند Geo-Reconstruct، HRIC یا CICSAM بسیار توصیه میشود.
2. اجرای شبیهسازی (Running the Simulation)
2.1. مقداردهی اولیه (Initialization)
- قبل از شروع محاسبات، باید شرایط اولیه برای همه متغیرها در کل دامنه محاسباتی تعریف شود.
- برای جریان دوفازی، میتوان یک فاز را به عنوان فاز پرکننده اولیه کل مخزن در نظر گرفت (مثلاً مخزن کاملاً پر از آب است) و سپس فاز دیگر (روغن) را از طریق ورودیها وارد کرد.
- مقداردهی اولیه میتواند به سرعت همگرایی و پایداری حل کمک کند.
2.2. اجرای محاسبات (Run Calculation)
- گام زمانی (Time Step Size): برای حلهای Transient، انتخاب اندازه گام زمانی بسیار مهم است. گام زمانی کوچکتر دقت را افزایش میدهد اما زمان محاسبات را به شدت زیاد میکند. گام زمانی بزرگتر میتواند منجر به ناپایداری و واگرایی حل شود. برای جریانهای VOF، گام زمانی اغلب توسط معیار Courant-Friedrichs-Lewy (CFL) کنترل میشود.
- تعداد گامهای زمانی (Number of Time Steps): تعیین تعداد گامها برای پوشش بازه زمانی مورد نظر از شبیهسازی.
- نظارت بر همگرایی (Monitoring Convergence):
- باقیماندهها (Residuals): نمودارهای باقیمانده نشاندهنده میزان برآورده شدن معادلات حاکم در هر تکرار هستند. باید تا حد معینی (مثلاً 1e-4 یا 1e-5) کاهش یابند تا حل همگرا محسوب شود.
- مانیتورینگ متغیرهای مورد علاقه: نظارت بر تغییرات متغیرهای مهم (مانند کسری حجم در یک نقطه خاص، یا دبی خروجی) در طول زمان میتواند نشاندهنده رسیدن به یک حالت پایدار یا دینامیک ثابت باشد.
3. تحلیل نتایج (Post-Processing and Analysis)
پس از اتمام محاسبات، نتایج باید استخراج و تفسیر شوند تا به درک پدیدههای فیزیکی و پاسخ به سؤالات مهندسی کمک کنند. این کار در ماژول Ansys CFD-Post یا مستقیماً در فلوئنت انجام میشود.
3.1. نمایش بصری (Visualization)
- نمایش کسری حجم (Volume Fraction Contours/Iso-surfaces): این مهمترین ابزار برای مشاهده موقعیت، شکل و حرکت رابط بین آب و روغن است. میتوان تصاویر یا انیمیشنهایی از توزیع فازها در زمانهای مختلف ایجاد کرد.
- نمایش میدانهای سرعت (Velocity Vectors/Contours): مشاهده الگوهای جریان برای هر فاز، شناسایی مناطق گردابهای، جریان معکوس یا نواحی با سرعت بالا.
- نمایش میدان فشار (Pressure Contours): بررسی توزیع فشار در مخزن و شناسایی مناطق پرفشار یا کمفشار.
- خطوط جریان (Streamlines): برای تجسم مسیرهای حرکت سیالات.
3.2. استخراج دادههای کمی (Quantitative Data Extraction)
- دبی جرمی یا حجمی (Mass/Volume Flow Rates): محاسبه دبی ورودی و خروجی هر فاز.
- نیروها (Forces): محاسبه نیروهای وارد بر دیوارهها (مثلاً نیروی درگ).
- تنش برشی (Shear Stress): محاسبه تنش برشی بر روی دیوارهها.
- نمودارهای تغییرات زمانی (Time-History Plots): ترسیم تغییرات متغیرهای خاص در یک نقطه یا در یک سطح مشخص در طول زمان (مثلاً ارتفاع رابط در یک مخزن).
- مقادیر متوسط یا انتگرالی (Averaged or Integrated Quantities): محاسبه مقادیر متوسط یا انتگرالی برای یک ناحیه یا کل دامنه (مثلاً میانگین سرعت فازها در یک مقطع).
3.3. تفسیر نتایج و اعتبار سنجی (Interpretation and Validation)
- درک فیزیک جریان: با بررسی تصاویر و دادهها، میتوان پدیدههایی مانند جدایش فازها، اختلاط، شکلگیری حبابها یا قطرات، و تأثیر کشش سطحی را درک کرد.
- اعتبار سنجی (Validation): مقایسه نتایج شبیهسازی با دادههای تجربی، نتایج سایر مدلهای عددی، یا حتی مشاهدات کیفی برای اطمینان از صحت مدل. این گام برای افزایش اعتماد به نتایج شبیهسازی حیاتی است.
- تحلیل حساسیت (Sensitivity Analysis): بررسی اینکه چگونه تغییر در پارامترهای ورودی (مثلاً کشش سطحی، دبی ورودی، یا خواص مواد) بر نتایج تأثیر میگذارد.
خلاصه چالشها در شبیهسازی و تحلیل
- پایداری عددی: به دلیل حضور رابط و کشش سطحی، حل میتواند ناپایدار باشد. نیاز به تنظیمات دقیق و گام زمانی کوچک.
- هزینه محاسباتی بالا: شبیهسازیهای Transient و VOF، به خصوص در سهبعدی و با مش ریز، منابع محاسباتی زیادی نیاز دارند.
- دشواری در اعتبار سنجی: جمعآوری دادههای تجربی دقیق برای جریانهای دوفازی میتواند چالشبرانگیز باشد.
- تفسیر پیچیده نتایج: حجم بالای دادههای تولید شده نیازمند مهارت بالا در تحلیل و تجسم است.
با درک عمیق این مراحل و ملاحظات، میتوان یک شبیهسازی موفق و معتبر از جریان دوفازی آب و روغن را در انسیس فلوئنت انجام داد و بینشهای ارزشمندی از پدیدههای فیزیکی مربوطه کسب کرد.
کلیدواژه ها : شبیهسازی و تحلیل جریان دوفازی آب و روغن- انسیس فلوئنت- حلگر مبتنی بر فشار- حلگر وابسته به زمان- مدل Volume of Fluid (VOF)- کشش سطحی بین فازها- فاز اولیه و ثانویه- مدل CSF- روشهای دیسکرتیزاسیون رابط مانند Geo-Reconstruct CICSAM HRIC- خواص مواد شامل چگالی ویسکوزیته خواص حرارتی- مدلهای تلاطم مانند k-epsilon k-omega SST- شرایط مرزی شامل ورودی سرعت ورودی جرم خروجی فشار دیواره شرط عدم لغزش- زاویه تماس- کوپلینگ فشار سرعت شامل PISO SIMPLE COUPLED- دیسکرتیزاسیون فضایی مرتبه دوم برای کسری حجم- مقداردهی اولیه- گام زمانی Courant Friedrichs Lewy- معیار CFL- نظارت بر همگرایی نمودار باقیماندهها- تحلیل نتایج- نمایش کسری حجم- نمایش میدان سرعت- نمایش میدان فشار- خطوط جریان- استخراج دادههای کمی- دبی جرمی- دبی حجمی- نیروهای وارد بر دیواره- تنش برشی- نمودارهای تغییرات زمانی- تفسیر نتایج- اعتبار سنجی مدل- تحلیل حساسیت- پایداری عددی- هزینه محاسباتی بالا- دشواری اعتبار سنجی- تفسیر پیچیده نتایج