مبانی شبیه‌سازی محیط‌های متخلخل و فیزیک مسئله کاتالیست این کارگاه آموزشی تخصصی با تمرکز بر یکی از پرکاربردترین مباحث دینامیک سیالات محاسباتی (CFD)، یعنی مدل‌سازی جریان در محیط‌های متخلخل (Porous Media) در نرم‌افزار ANSYS CFX طراحی شده است. پروژه محوری این آموزش، شبیه‌سازی جریان گاز نیتروژن (به عنوان گاز ایده‌آل) درون یک مبدل کاتالیستی (Catalytic Converter) است. گاز نیتروژن با سرعت یکنواخت ۱۰ متر بر ثانیه و دمای ۴۵۰ درجه سانتی‌گراد وارد دستگاه شده و از درون یک بستر سرامیکی یکپارچه (Ceramic Monolith Substrate) عبور می‌کند. این بستر دارای کانال‌هایی با مقطع مربع است که جریان را در راستای محور طولی (Z) عبور می‌دهد، اما در راستای محورهای عرضی (X و Y) غیرقابل نفوذ است. مخاطبان در این بخش یاد می‌گیرند که چگونه این ناهمسانگردی فیزیکی را با تعریف ضریب مقاومت جانبی که دو مرتبه بزرگ‌تر از مقاومت طولی است، به صورت ریاضی در نرم‌افزار شبیه‌سازی کنند. این آموزش برای دانشجویان و مهندسانی که می‌خواهند فرآیندهای پیچیده افت فشار و رفتار جریان در فیلترها و کاتالیزورها را مدل‌سازی کنند، بسیار ارزشمند است.

راهنمای گام‌به‌گام پیاده‌سازی پروژه در محیط CFX-Pre و حلگر مخاطبان پس از مطالعه و اجرای گام‌به‌گام این کارگاه، فرآیند کامل راه‌اندازی یک پروژه متخلخل را فرا می‌گیرند. این مراحل با ذخیره پروژه و وارد کردن شبکه محاسباتی (Mesh) با فرمت ICEM CFD آغاز می‌شود. سپس فراخوانی گاز نیتروژن از کتابخانه مواد و تنظیم دامنه‌ سیال پیش‌فرض با مدل آشفتگی k-epsilon و انتقال حرارت هم‌دما صورت می‌گیرد. گام کلیدی این آموزش، نحوه تعریف «دامنه متخلخل» (Porous Domain) برای بستر سرامیکی است. کاربران می‌آموزند که چگونه تخلخل حجمی را روی 0.5 تنظیم کنند، مدل افت جهت‌دار (Directional Loss) را فعال نمایند و ضرایب مقاومت خطی و درجه دوم را برابر با ۴۴۰ کیلوگرم بر متر به توان چهار (\(kg\cdot m^{-4}\)) اعمال کنند. علاوه بر این، نحوه فرمول‌نویسی برای مانیتور کردن افت فشار و دبی جرمی خروجی در حین حل مسئله آموزش داده می‌شود تا کاربر بتواند همگرایی فیزیکی مسئله را در Solver Manager به صورت لحظه‌ای ردیابی کند.

تحلیل مهندسی نتایج و محاسبات افت فشار در CFD-Post بخش پایانی این آموزش به تحلیل کیفی و کمی نتایج در محیط CFD-Post اختصاص دارد. کاربران یاد می‌گیرند که چگونه با ترسیم بردارهای سرعت و کانتورهای فشار روی صفحات برشی، رفتار بصری گاز را ارزیابی کنند. نتایج کمی حلگر نشان می‌دهد که افت فشار خالص در طول بستر متخلخل سرامیکی حدود ۲۸۵ پاسکال است. با این حال، با تعریف یک عبارت ریاضی (Expression) برای محاسبه اختلاف فشار کل بین ورودی و خروجی کل دستگاه، مقدار افت فشار کل حدود ۳۰۰ پاسکال به دست می‌آید. این تحلیل به مخاطب آموزش می‌دهد که چگونه با مقایسه این دو عدد نتیجه‌گیری کند که ۲۸۵ پاسکال از فشار به دلیل مقاومت بستر متخلخل افت کرده و ۱۵ پاسکال باقی‌مانده ناشی از تلفات جریان در سایر بخش‌های دستگاه (مانند نواحی انبساط و انقباض مجاری ورودی و خروجی) است. این تحلیل بودجه فشار، مهارتی بسیار حیاتی برای بهینه‌سازی طراحی و کاهش فشار برگشتی (Backpressure) موتورها به شمار می‌رود.

۱. مقدمه: چالش مدل‌سازی هزاران مجرای میکروسکوپی

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

اما یک سوال استراتژیک وجود دارد: چگونه می‌توان این حجم عظیم از جزئیات هندسی را شبیه‌سازی کرد بدون اینکه ابررایانه‌ها تحت فشار محاسباتی از کار بیفتند؟ پاسخ در جادوی مفهومی به نام «محیط متخلخل» (Porous Media) نهفته است که به ما اجازه می‌دهد بدون ترسیم تک‌تک سوراخ‌ها، فیزیک واقعی را مدل‌سازی کنیم.

۲. نکته اول: جادوی ساده‌سازی؛ وقتی کل از اجزا فراتر می‌رود

در شبیه‌سازی CFD، همیشه هوشمندانه‌ترین راه، بازسازی دقیق هندسه نیست. برای مدل‌سازی بستر سرامیکی (Ceramic Monolith)، استفاده از یک «دامنه متخلخل» یا Porous Domain راهکاری است که پیچیدگی را به شدت کاهش می‌دهد.

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

“این کارگاه مدل‌سازی یک مبدل کاتالیستی را نمایش می‌دهد که در آن گاز نیتروژن از میان یک بستر منولیت سرامیکی با کانال‌های مربع‌شکل عبور می‌کند.”

نکته فنی برجسته در اینجا، تنظیم پارامتر «تخلخل حجمی ۰.۵» (Volume Porosity 0.5) است. از دیدگاه تحلیلی، این عدد به معنای آن است که نیمی از حجم این ناحیه توسط دیواره‌های سرامیکی مسدود شده است. این موضوع مستقیماً بر سرعت محلی (Local Velocity) تاثیر می‌گذارد؛ به طوری که سرعت سیال در هنگام عبور از مجاری بستر، به دلیل کاهش سطح مقطع موثر، عملاً دو برابر می‌شود. این ساده‌سازی به ما اجازه می‌دهد بدون درگیر شدن با مش‌بندی میلیون‌ها کانال، تاثیر واقعی بستر بر جریان را محاسبه کنیم.

۳. نکته دوم: مهندسی جهت‌دهی؛ چگونه جریان را در مسیر دلخواه مهار کنیم؟

در دنیای واقعی، کانال‌های بستر کاتالیست اجازه حرکت جانبی به گاز را نمی‌دهند. برای بازسازی این رفتار فیزیکی در نرم‌افزار، مهندسان از تکنیک «افت جهت‌دار» (Directional Loss) استفاده می‌کنند.

در تنظیمات این پروژه (مطابق با)، جریان در راستای محور Z با مولفه‌های (0,0,1) هدایت شده است. اما ترفند اصلی در بخش Transverse Loss نهفته است:

• ضریب مقاومت درجه دوم: 440 kg/m^4.
• ضریب ضرب‌کننده جهت جانبی (Multiplier): ۱۰۰.

با اعمال این ضریب ۱۰۰ در جهت‌های X و Y (یعنی ۲ مرتبه بزرگی بیشتر از جهت اصلی)، بستر در برابر جریان‌های عرضی عملاً غیرقابل نفوذ می‌شود. تایید بصری این مهندسی دقیق را می‌توان در [SOURCE_IMAGE_8] مشاهده کرد؛ جایی که وکتورهای سرعت کاملاً مستقیم و بدون کوچک‌ترین انحراف جانبی از بستر عبور می‌کنند، دقیقاً همان‌گونه که دیواره‌های صلب سرامیکی در واقعیت عمل می‌کنند.

۴. نکته سوم: توزیع نابرابر فشار؛ گلوگاه اصلی کجاست؟

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

طبق پایش لحظه‌ای حلگر ثبت شده است، مقدار دقیق افت فشار در بستر متخلخل بر روی عدد ۲۸۵.۸۹ پاسکال تثبیت می‌شود. با توجه به اینکه افت فشار کل دستگاه حدود ۳۰۰ پاسکال تخمین زده شده است، تحلیل داده‌ها نشان می‌دهد:

• ۹۵٪ از کل افت فشار (۲۸۵ پاسکال) منحصراً در بستر متخلخل رخ می‌دهد.
• تنها ۵٪ باقی‌مانده (۱۵ پاسکال) سهم سایر بخش‌های دستگاه شامل ورودی‌ها و خروجی‌هاست.

این یعنی بستر سرامیکی، Bottleneck یا گلوگاه اصلی سیستم است و هرگونه تلاش برای بهبود کارایی باید از همین نقطه آغاز شود.

۵. نتیجه‌گیری: فراتر از مبدل‌های کاتالیستی

ابزارهایی مانند CFX به ما می‌آموزند که در سیستم‌های پیچیده مهندسی، همیشه شهود اولیه ما درست نیست. وقتی شبیه‌سازی با چنین دقتی ثابت می‌کند که ۹۵٪ مقاومت در یک بخش کوچک متمرکز شده است، استراتژی مهندسی تغییر می‌کند.

اندیشه نهایی: اگر بدانید سهم بدنه فلزی در اتلاف انرژی ناچیز است، آیا باز هم زمان خود را صرف آیرودینامیک کردن بدنه می‌کنید یا تمام تمرکزتان را بر روی «بهینه‌سازی توزیع جریان در ورودی» (Inlet Flow Distribution) و مهندسی تخلخل می‌گذارید؟ شناخت گلوگاه‌ها، تفاوت میان یک طراحی معمولی و یک شاهکار مهندسی است.