شبیه سازی دو بعدی احتراق در انسیس فلوئنت

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

شبیه‌سازی دو بعدی احتراق، یک ابزار قدرتمند برای درک پدیده‌های اساسی احتراق با هزینه محاسباتی نسبتاً پایین‌تر نسبت به شبیه‌سازی‌های سه بعدی است. اگرچه شبیه‌سازی دو بعدی نمی‌تواند تمامی جزئیات پیچیده سه بعدی (مانند گردابه‌های سه‌بعدی یا اثرات چرخش) را به طور کامل مدل کند، اما برای مطالعه بسیاری از جنبه‌های کلیدی مانند انتشار شعله، توزیع دما و گونه‌ها، و اثرات اختلاط در هندسه‌های خاص (مانند نازل‌ها یا کانال‌های ساده) بسیار مفید است. نرم‌افزار ANSYS Fluent به دلیل قابلیت‌های گسترده در مدل‌سازی فیزیک‌های مختلف، یکی از ابزارهای پیشرو برای این منظور است.

مراحل شبیه‌سازی احتراق در ANSYS Fluent

شبیه‌سازی احتراق در Fluent یک فرآیند چند مرحله‌ای است که نیازمند دقت در انتخاب مدل‌ها، تعریف خواص مواد و شرایط مرزی است.

۱. تعریف هندسه و مش‌بندی (Geometry & Meshing)

• هندسه (Geometry): هندسه دو بعدی سیستم احتراق (مانند یک محفظه احتراق، مشعل، یا کانال) باید با استفاده از ابزارهای CAD مانند ANSYS DesignModeler یا SpaceClaim ایجاد شود. این هندسه شامل مرزهای ورودی سوخت و اکسیدکننده، دیواره‌های محفظه و خروجی محصولات احتراق است.
• مش‌بندی (Meshing): کیفیت مش در شبیه‌سازی احتراق از اهمیت فوق‌العاده‌ای برخوردار است. نواحی واکنش (Flame Zone) که در آن گرادیان‌های دما و غلظت گونه‌ها بسیار شدید هستند، نیازمند مش‌بندی بسیار ریز (Refined Mesh) هستند. استفاده از مش‌بندی تطبیقی (Adaptive Meshing) یا Refinement در نواحی مهم می‌تواند به ثبت دقیق شعله و کاهش خطای عددی کمک کند.
• مش‌بندی سازمان‌یافته (Structured Mesh): در صورت امکان، به دلیل کیفیت بالاتر، توصیه می‌شود.
• مش‌بندی بی‌سازمان (Unstructured Mesh): برای هندسه‌های پیچیده‌تر.
• لایه مرزی (Boundary Layer Meshing): در نزدیکی دیواره‌ها برای ثبت دقیق لایه مرزی حرارتی و سیالاتی ضروری است.

۲. تنظیمات اولیه حلگر (Solver Settings)

• نوع حلگر (Solver Type):
• Pressure-Based: معمولاً برای جریان‌های احتراق با عدد ماخ پایین و تراکم‌ناپذیر یا با تراکم‌پذیری کم استفاده می‌شود.
• Density-Based: برای جریان‌های احتراق با سرعت بالا و تراکم‌پذیر، به ویژه در جایی که نوسانات فشار قابل توجه است (مانند موتورهای راکت)، مناسب‌تر است.
• وابستگی زمانی (Time):
• Steady: برای شبیه‌سازی احتراق در حالت پایدار، جایی که شعله به موقعیت ثابتی می‌رسد و خواص جریان با زمان تغییر نمی‌کند.
• Transient: برای مطالعه فرآیندهای ناپایدار مانند آغاز احتراق، انتشار شعله، پایداری شعله یا نوسانات احتراق. بسیاری از پدیده‌های احتراقی ذاتاً ناپایدار هستند.
• ابعاد (2D Space): اطمینان از انتخاب گزینه “2D”.

۳. تعریف مواد و واکنش‌ها (Materials & Reactions)

• مدل گونه‌ها (Species Transport): این مدل باید فعال شود تا بتوان انتقال و واکنش شیمیایی گونه‌های مختلف (سوخت، اکسیدکننده، محصولات احتراق و گونه‌های میانی) را مدل‌سازی کرد.
• ایجاد مخلوط (Mixture Material): یک “Mixture Material” شامل سوخت، اکسیدکننده (مانند N2, O2) و محصولات احتراق (مانند CO2, H2O) باید در Fluent تعریف شود. خواص ترموفیزیکی هر گونه (مانند ویسکوزیته، رسانایی حرارتی، گرمای ویژه) و همچنین خواص مخلوط (مانند چگالی بر اساس قانون گاز ایده‌آل) باید مشخص شود.
• تعریف واکنش‌ها (Reactions):
• Single-Step Reaction: برای ساده‌سازی فرآیند احتراق به یک واکنش کلی (مانند واکنش متان با اکسیژن).
• Multi-Step Reaction: برای دقت بیشتر، می‌توان از مکانیسم‌های واکنش چند مرحله‌ای استفاده کرد که شامل گونه‌های میانی و واکنش‌های اولیه هستند. مکانیسم‌های اسکلتی (Skeletal Mechanisms) نسخه‌های ساده شده‌ای از مکانیسم‌های جامع (Detailed Mechanisms) هستند که پیچیدگی محاسباتی را کاهش می‌دهند.
• Finite-Rate Chemistry: این گزینه محاسبه سرعت واکنش‌ها را بر اساس قوانین سینتیک شیمیایی (مانند قانون آرنیوس) فعال می‌کند.

۴. مدل‌های احتراق (Combustion Models)

انتخاب مدل احتراق مناسب به نوع شعله (آمیخته یا پیش‌آمیخته) و میزان آشفتگی بستگی دارد.

• احتراق غیر پیش‌آمیخته (Non-Premixed Combustion): در این نوع احتراق، سوخت و اکسیدکننده به طور جداگانه وارد محفظه می‌شوند و احتراق در ناحیه اختلاط آنها رخ می‌دهد.
• Eddy-Dissipation Model (EDM): یک مدل ساده و مقاوم که فرض می‌کند سرعت واکنش توسط اختلاط آشفته کنترل می‌شود. برای شعله‌های پخششی آشفته که اختلاط عامل محدودکننده است، مناسب است.
• Eddy-Dissipation Concept (EDC): مدل پیچیده‌تر که امکان مدل‌سازی سینتیک شیمیایی را در داخل عناصر ریز آشفته فراهم می‌کند. دقت بالاتری نسبت به EDM دارد.
• Flamelet Model: برای احتراق با سینتیک شیمیایی پیچیده، جایی که شعله به عنوان مجموعه‌ای از شعله‌های لایه‌ای (Flamelets) توصیف می‌شود.
• PDF Transport Model: دقیق‌ترین مدل برای احتراق‌های آشفته غیر پیش‌آمیخته، که اثرات نوسانات آشفتگی بر واکنش‌ها را به طور مستقیم مدل می‌کند. از نظر محاسباتی پرهزینه است.
• احتراق پیش‌آمیخته (Premixed Combustion): در این حالت، سوخت و اکسیدکننده قبل از احتراق به طور کامل مخلوط شده‌اند و شعله به صورت یک جبهه واکنش (Reaction Front) منتشر می‌شود.
• G-Equation (Level Set): برای مدل‌سازی انتشار جبهه شعله نازک در جریان‌های آشفته و آرام.
• Flamelet Generated Manifold (FGM): برای کاهش هزینه محاسباتی مکانیسم‌های شیمیایی پیچیده در شعله‌های پیش‌آمیخته.
• احتراق نیمه‌پیش‌آمیخته (Partially Premixed Combustion): ترکیبی از دو حالت فوق.
• مدل تابش (Radiation Model): در احتراق‌های با دمای بالا، انتقال حرارت تابشی نقش مهمی ایفا می‌کند.
• P1 Model: ساده و نسبتاً ارزان، اما برای سیستم‌های اپتیکال نازک یا ضخیم مناسب نیست.
• Discrete Ordinates (DO) Model: دقیق‌تر و انعطاف‌پذیرتر، برای طیف وسیعی از کاربردها.

۵. مدل آشفتگی (Turbulence Model)

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

• k-epsilon (k-ε) Models: پرکاربردترین مدل‌ها (Standard, Realizable, RNG). مقاوم و مناسب برای بسیاری از جریان‌های صنعتی.
• k-omega (k-ω) SST Model: برای جریان‌های با جدایش و گرادیان فشار معکوس، به ویژه نزدیک دیواره‌ها، دقت بالایی دارد و اغلب برای احتراق توصیه می‌شود.
• Large Eddy Simulation (LES): برای شبیه‌سازی با وضوح بالا، که به منابع محاسباتی زیادی نیاز دارد. معمولاً برای شبیه‌سازی‌های سه بعدی و در صورت نیاز به جزئیات دقیق‌تر ساختار شعله استفاده می‌شود.

۶. شرایط مرزی (Boundary Conditions)

شرایط مرزی باید به دقت و با توجه به فیزیک مسئله تعریف شوند.

• ورودی (Inlet):
• Mass Flow Inlet / Velocity Inlet: برای جریان ورودی سوخت و اکسیدکننده. باید مقادیر غلظت گونه‌ها (Mass Fractions)، دما و پارامترهای آشفتگی (مانند شدت آشفتگی و مقیاس طول) را مشخص کنید.
• خروجی (Outlet):
• Pressure Outlet: فشار استاتیک در خروجی (معمولاً فشار اتمسفریک).
• دیواره‌ها (Walls):
• Thermal Conditions: دمای دیواره (Temperature)، شار حرارتی ثابت (Heat Flux)، یا آدیاباتیک (Adiabatic) برای دیواره‌های عایق.
• No-Slip: شرط سرعت صفر در دیواره.

۷. روش‌های حل و همگرایی (Solution Methods & Convergence)

• Schemes (طرح‌های تفاضل‌گیری): برای معادلات مومنتوم، انرژی، گونه‌ها و آشفتگی، استفاده از طرح‌های Second Order Upwind برای دقت بالاتر توصیه می‌شود. برای جفت‌شدگی فشار-سرعت (Pressure-Velocity Coupling)، طرح‌های SIMPLEC یا PISO (برای مسائل Transient) یا Coupled (برای سرعت و مقاومت بیشتر) مناسب هستند.
• Residuals (باقیمانده‌ها): معیار همگرایی حل. مقادیر باقیمانده برای معادلات مومنتوم و گونه‌ها معمولاً ${10}^{-5}$10−5 و برای انرژی ${10}^{-7}$10−7 در نظر گرفته می‌شود. همچنین، نظارت بر توازن جرم و انرژی (Mass and Energy Balance) نیز برای اطمینان از همگرایی مهم است.

۸. مقداردهی اولیه و اجرا (Initialization & Run Calculation)

• مقداردهی اولیه (Initialization): مقادیر اولیه برای میدان‌های جریان (سرعت، فشار، دما، غلظت گونه‌ها) را تعیین کنید. اغلب از “Standard Initialization” استفاده می‌شود.
• Patching: در برخی موارد، به ویژه برای شروع احتراق، ممکن است نیاز باشد یک ناحیه کوچک با دمای بالا یا غلظت گونه‌های محصولات احتراق (Ignition Source) را به صورت دستی در هندسه “Patch” کنید تا فرآیند احتراق آغاز شود.
• Run Calculation: تعداد تکرارها (Iterations) یا گام‌های زمانی (Time Steps) را برای رسیدن به همگرایی یا پایان شبیه‌سازی ناپایدار مشخص کنید.

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

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

• کانتورها (Contours):
• Temperature Contours: برای مشاهده توزیع دما و محل جبهه شعله.
• Species Mass Fraction Contours: توزیع غلظت سوخت، اکسیدکننده، محصولات احتراق (CO2, H2O) و آلاینده‌ها (CO, NOx).
• Heat Release Rate Contours: نشان‌دهنده نواحی اصلی واکنش شیمیایی و تولید گرما.
• Velocity Contours / Vectors: برای بررسی الگوهای جریان، مناطق رکود یا جدایش.
• پلات‌های X-Y (X-Y Plots): رسم پروفایل دما و غلظت گونه‌ها در امتداد خطوط مشخص در محفظه احتراق.
• گزارش‌ها (Reports):
• Integral Quantities: محاسبه شار جرمی گونه‌ها در ورودی و خروجی، تولید گرمای کلی.
• Average Quantities: دمای میانگین، غلظت میانگین گونه‌ها در خروجی (برای تحلیل آلایندگی).
• تجزیه و تحلیل پایداری شعله و آلایندگی: بررسی شکل و پایداری شعله و همچنین میزان تولید آلاینده‌های حرارتی (NOx) و ناقص (CO).

چالش‌ها و نکات مهم در شبیه‌سازی احتراق دو بعدی

• تعامل آشفتگی-شیمی (Turbulence-Chemistry Interaction): این بزرگترین چالش در شبیه‌سازی احتراق است، زیرا آشفتگی بر سرعت اختلاط و در نتیجه سرعت واکنش تأثیر می‌گذارد، و از سوی دیگر گرما و محصولات واکنش بر خواص سیال و آشفتگی تأثیر می‌گذارند. انتخاب مدل احتراق و آشفتگی مناسب برای مدیریت این تعامل حیاتی است.
• پیچیدگی مکانیسم واکنش (Reaction Mechanism Complexity): مکانیسم‌های جامع می‌توانند بسیار بزرگ و از نظر محاسباتی سنگین باشند. انتخاب یک مکانیسم اسکلتی مناسب یا استفاده از رویکردهای Flamelet می‌تواند این چالش را کاهش دهد.
• اثرات تابش (Radiation Effects): در دماهای بالا، تابش گرمایی می‌تواند بخش قابل توجهی از انتقال حرارت را شامل شود و باید به درستی مدل‌سازی شود.
• پایداری عددی (Numerical Stability): معادلات واکنش‌های شیمیایی اغلب “Stiff” هستند و ممکن است نیاز به گام‌های زمانی بسیار کوچک در حل ناپایدار داشته باشند.
• محدودیت‌های دو بعدی (Limitations of 2D): شبیه‌سازی دو بعدی نمی‌تواند پدیده‌های ذاتی سه بعدی مانند چرخش (Swirl)، نوسانات سه بعدی شعله یا گردابه‌های پیچیده را به طور کامل نشان دهد. در صورت اهمیت این پدیده‌ها، شبیه‌سازی سه بعدی ضروری است.

شبیه‌سازی احتراق دو بعدی در ANSYS Fluent ابزاری قدرتمند برای تحلیل و بهینه‌سازی سیستم‌های احتراق است، اما نیازمند دانش دقیق از اصول احتراق، مکانیک سیالات محاسباتی و قابلیت‌های نرم‌افزار است.

کلیدواژه ها : شبیه-سازی-احتراق-Combustion-Simulation-انسیس-فلوئنت-ANSYS-Fluent-دو-بعدی-2D-مدل-احتراق-Combustion-Model-Eddy-Dissipation-Model-EDM-Eddy-Dissipation-Concept-EDC-Flamelet-Model-PDF-Transport-مدل-گونه-ها-Species-Transport-مدل-آشفتگی-Turbulence-Model-k-epsilon-k-omega-SST-مکانیسم-واکنش-Reaction-Mechanism-مقدمه-Introduction-هندسه-Geometry-مش-بندی-Meshing-حلگر-Solver-خواص-مواد-Materials-شرایط-مرزی-Boundary-Conditions-مقداردهی-اولیه-Initialization-پس-پردازش-Post-processing-تحلیل-نتایج-Result-Analysis-توزیع-دما-Temperature-Distribution-گونه-ها-Species-جبهه-شعله-Flame-Front-آلاینده-ها-Pollutants-NOx-CO-چالش-ها-Challenges-تعامل-آشفتگی-شیمی-Turbulence-Chemistry-Interaction-مدل-تابش-Radiation-Model-پایداری-شعله-Flame-Stability