راهنمای کامل شرایط مرزی برای شبیه‌سازی توربوماشین‌ها (Turbomachinery)

مبانی و چالش‌های کلیدی شرایط مرزی

شبیه‌سازی توربوماشین‌ها، به جای تحلیل کل ماشین، اغلب بر روی یک یا چند ردیف از پره‌ها متمرکز می‌شود. این رویکرد نیازمند تعریف دقیق مرزهای دامنه محاسباتی است. چالش اصلی در این فرآیند، کوتاه کردن دامنه محاسباتی (truncation of the computational domain) است که منجر به ایجاد مرزهای باز مصنوعی (artificial open boundaries) می‌شود؛ مرزهایی که در دنیای واقعی وجود ندارند اما برای حل عددی مسئله ضروری هستند.

این مرزهای مصنوعی می‌توانند باعث ایجاد «بازتاب‌های عددی کاذب» شوند. این پدیده زمانی رخ می‌دهد که امواج فشار یا سایر اختلالات جریان به مرز دامنه می‌رسند و به جای عبور طبیعی از آن، به صورت غیرفیزیکی به داخل دامنه بازتابیده می‌شوند. این امواج ناخواسته می‌توانند میدان جریان داخلی را مختل کرده، پیش‌بینی بار آیرودینامیکی روی پره‌ها را مخدوش کرده و نتایج شبیه‌سازی را بی‌اعتبار سازند.

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

• مرز ورودی (Inlet): محل ورود جریان به دامنه محاسباتی که در آن مشخصات جریان ورودی تعریف می‌شود.
• مرز خروجی (Outlet): محل خروج جریان از دامنه که شرایط فشار پایین‌دست را تعیین می‌کند.
• دیواره‌ها (Walls): سطوح جامد مانند پره‌ها، هاب (hub) و شرود (shroud) که جریان از آن‌ها عبور نمی‌کند.
• مرزهای واسط (Interface): مرزهای بین بخش‌های ثابت (استاتور) و دوار (روتور) که اندرکنش بین آن‌ها را مدل می‌کنند (مانند Mixing Plane).
• مرزهای تناوبی (Periodic): برای شبیه‌سازی یک کانال جریان (passage) از بین چندین کانال مشابه به منظور کاهش هزینه‌های محاسباتی (مانند cyclicAMI).

در ادامه، به بررسی تخصصی شرایط مرزی ورودی و انواع آن می‌پردازیم.

انواع شرایط مرزی ورودی

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

 شرایط مرزی مبتنی بر فشار (Pressure-Based)

شرط مرزی «فشار کل» (Total Pressure) (که در نرم‌افزارهایی مانند ANSYS CFX به عنوان Total Pressure at Inlet و در Fluent به عنوان Pressure Inlet شناخته می‌شود) یکی از رایج‌ترین و پایدارترین روش‌ها برای تعریف ورودی است. در این حالت، فشار کل در ورودی مشخص می‌شود و حلگر بر اساس شرایط پایین‌دست (مانند فشار استاتیک خروجی)، سرعت و دبی جرمی را محاسبه می‌کند. این رویکرد به ویژه زمانی مفید است که شرایط عملکردی ماشین بر اساس نسبت فشار تعریف شود.

شرایط مرزی مبتنی بر جریان (Flow-Based)

شرط مرزی «دبی جرمی» (Mass Flow Inlet) زمانی به کار می‌رود که دبی جرمی عبوری از ماشین، یک پارامتر طراحی مشخص و ثابت باشد. مزیت این روش آن است که دبی جرمی را مستقیماً کنترل می‌کند و برای شبیه‌سازی نقاط عملکردی مشخص در منحنی مشخصه کمپرسور یا پمپ بسیار مناسب است.

 شرایط مرزی مبتنی بر سرعت (Velocity-Based)

شرط مرزی «سرعت ورودی» (Velocity Inlet) سرعت جریان را مستقیماً در مرز ورودی تعریف می‌کند. این روش کمتر در شبیه‌سازی توربوماشین‌های تراکم‌پذیر استفاده می‌شود، اما برای کاربردهای خاص مانند شبیه‌سازی توربین‌های آبی داخل لوله که دبی و در نتیجه سرعت متوسط ورودی مشخص است، کاربرد دارد.

علاوه بر متغیرهای اصلی جریان، تعریف مشخصات تلاطم در ورودی نیز برای دستیابی به نتایج دقیق ضروری است.

ملاحظات ویژه برای ورودی: تعریف تلاطم

تعریف شرایط تلاطم ورودی واقعی (realistic turbulence inlet conditions) یکی از بخش‌های مهم و معمولاً دشوار در تنظیم شبیه‌سازی است. این شرایط می‌توانند به طور قابل توجهی بر توسعه لایه مرزی روی پره‌ها، پدیده جدایش جریان و میزان تلفات آیرودینامیکی تأثیرگذار باشند.

روش‌های رایج برای تعریف تلاطم در ورودی بر اساس مدل‌های توربولانسی دو معادله‌ای (مانند k-ε یا k-ω) به شرح زیر است:

سطح تلاطم و مقیاس طول (Turbulence Level and Length-Scale): این دو پارامتر رایج‌ترین روش برای تعریف تلاطم ورودی هستند. سطح تلاطم (Turbulence Intensity) به صورت درصد بیان می‌شود و مقیاس طول (Length-Scale) اندازه مشخصه گردابه‌های ورودی را نشان می‌دهد. در صورت عدم وجود داده‌های تجربی، مقیاس طول را می‌توان بر اساس مشخصات هندسی قطعات بالادست (مانند ضخامت پره‌ها یا درصدی از ارتفاع کانال) تخمین زد.

تخمین‌های تجربی برای قطعات مختلف: در غیاب داده‌های اندازه‌گیری شده، می‌توان از مقادیر تخمینی زیر برای «سطح تلاطم» استفاده کرد:

• • توربین فشار بالا (بعد از محفظه احتراق): تا ۲۰٪ (به دلیل تلاطم بسیار بالای ناشی از احتراق)
  • کمپرسورهای فشار بالا یا توربین‌های فشار پایین: حدود ۵٪
  • فن‌ها و کمپرسورهای فشار پایین: تا ۱٪ (به دلیل وجود قطعات کمتر در بالادست)

نسبت ویسکوزیته گردابه‌ای (Eddy Viscosity Ratio): این روش جایگزینی برای مقیاس طول است و کنترل مستقیم‌تری بر تأثیر تلاطم ورودی اعمال می‌کند. این نسبت، ویسکوزیته توربولانسی را به ویسکوزیته مولکولی مرتبط می‌سازد. مقادیر نمونه برای این پارامتر از ۱۰ برای فن‌ها و کاربردهای با تلاطم پایین، تا ۱۰۰۰ یا در موارد نادر حتی ۱۰۰۰۰ برای توربین‌های فشار بالا متغیر است.

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

انواع شرایط مرزی خروجی

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

 شرایط مرزی مبتنی بر فشار استاتیک (Static Pressure-Based)

شرط مرزی «فشار استاتیک» (fixedValue Static Pressure) یکی از متداول‌ترین و پایدارترین روش‌ها برای مرز خروجی است. در این روش، یک مقدار ثابت برای فشار استاتیک در کل مرز خروجی تعریف می‌شود. این شرط به همراه ورودی فشار کل، یک ترکیب بسیار رایج و قوی برای شبیه‌سازی‌های توربوماشین است. یک جایگزین قوی‌تر، شرط مرزی fixedMeanValue است که در آن فشار به گونه‌ای تنظیم می‌شود تا میانگین مورد نظر در خروجی حفظ شود.

 شرایط مرزی گرادیان صفر (Zero Gradient)

شرط مرزی «گرادیان صفر» (zeroGradient) فرض می‌کند که گرادیان تمام متغیرهای جریان در جهت عمود بر مرز خروجی صفر است. این شرط به جریان اجازه می‌دهد تا به طور طبیعی و با توجه به فیزیک بالادست، از دامنه خارج شود. این رویکرد زمانی مناسب است که مرز خروجی به اندازه کافی از ناحیه مورد نظر (مانند ردیف پره‌ها) دور باشد و جریان در آنجا کاملاً توسعه‌یافته باشد. استفاده از شرط zeroGradient در فاصله‌ای بسیار نزدیک به لبه‌های فرار پره‌ها می‌تواند به صورت مصنوعی از توسعه کامل دنباله (wake) و گردابه‌های پایین‌دست جلوگیری کرده و منجر به تخمین کمتر از واقعیت تلفات ناشی از اختلاط شود.

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

شرایط مرزی پیشرفته: جلوگیری از بازتاب‌های عددی

شرایط مرزی استانداردی که پیش‌تر بررسی شدند، مانند فشار استاتیک ثابت، اغلب در جلوگیری از همان «بازتاب‌های عددی کاذب» که در مقدمه به عنوان منبع اصلی خطا معرفی شد، ناتوان هستند. برای تحلیل‌هایی که این بازتاب‌ها در آن‌ها غیرقابل تحمل است—به ویژه در شبیه‌سازی‌های گذرا (unsteady)، آیرواکوستیک یا آئروالاستیک—نیاز به یک فرمولاسیون پیشرفته‌تر است. در این موارد، امواج فشاری که به مرزها می‌رسند، به صورت غیرفیزیکی بازتاب شده و میدان حل را آلوده می‌سازند.

شرایط مرزی غیر بازتابنده (Non-Reflecting Boundary Conditions – NRBC)

این دسته از شرایط مرزی برای به حداقل رساندن بازتاب‌های عددی طراحی شده‌اند. رویکرد توسعه داده شده توسط Giles یکی از محبوب‌ترین روش‌ها در این زمینه است. این روش مبتنی بر تجزیه ریاضی میدان جریان به امواج مشخصه (characteristic waves) است. این امواج به صورت فیزیکی به موج‌های صوتی، گردابه‌ای (vorticity) و آنتروپی (entropy) مربوط می‌شوند. با شناسایی جهت انتشار هر یک از این امواج، NRBC اجازه خروج امواج از دامنه را می‌دهد، اما از ورود امواج ناخواسته و غیرفیزیکی از مرز به داخل دامنه جلوگیری می‌کند. اهمیت حیاتی این رویکرد صرفاً نظری نیست؛ همان‌طور که در مطالعه موردی توربین VKI LS89 در منابع فنی نشان داده شده است، استفاده از شرایط مرزی ریمان استاندارد یا تقریبی می‌تواند امواج ضربه‌ای (shock waves) را محو کرده و به پیش‌بینی کاملاً اشتباه بارگذاری روی پره منجر شود، در حالی که یک NRBC مناسب، فیزیک مسئله را به درستی ثبت می‌کند.

Giles دو فرمولاسیون متفاوت ارائه کرده است:

1. یک فرمولاسیون دقیق برای حالت پایا (steady) که غیرمحلی (non-local) است (نیازمند تبدیل فوریه در فضا است، که باعث می‌شود شرایط در یک نقطه به کل مرز وابسته باشد).
2. یک فرمولاسیون تقریبی برای حالت گذرا (unsteady) که برای محلی (local) بودن طراحی شده است (که آن را از نظر محاسباتی ارزان‌تر می‌کند اما برای انواع خاصی از امواج دقت کمتری دارد).

 شرایط مرزی مبتنی بر ناورداهای ریمان (Riemann Invariant-Based)

این روش، که در مرز ورودی استفاده می‌شود، نوع دیگری از شرایط مرزی یک‌بعدی غیر بازتابنده است. این شرایط مرزی بر اساس ناورداهای ریمان (Riemann Invariants) بنا شده‌اند تا از بازتاب امواج در مرز ورودی جلوگیری کنند و پایداری حل را بهبود بخشند.

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

شرایط مرزی واسط در تحلیل‌های چندمرحله‌ای

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

• صفحات اختلاط (Mixing Planes): این روش به طور گسترده برای شبیه‌سازی‌های پایا (steady) استفاده می‌شود. در این رویکرد، متغیرهای جریان در سطح واسط به صورت محیطی میانگین‌گیری می‌شوند تا یک پروفایل جریان یکنواخت شعاعی ایجاد شود. این کار عملاً تمام اثرات گذرا و غیریکنواختی‌های محیطی (مانند دنباله پره‌ها) را حذف می‌کند و یک شرایط مرزی پایدار برای ردیف پره بعدی فراهم می‌آورد. با وجود کارایی محاسباتی بالا برای تحلیل‌های پایا، این فرآیند میانگین‌گیری به این معناست که صفحات اختلاط اساساً قادر به ثبت پدیده‌های گذرا مانند انتشار دنباله (wake propagation)، پدیده‌های فرکانس عبور پره (blade passing frequency) یا اندرکنش‌های آکوستیکی بین مراحل نیستند؛ پدیده‌هایی که اغلب برای پیش‌بینی فلاتر پره یا نویز، حیاتی هستند.
• شرایط مرزی تناوبی/چرخشی (Periodic/Cyclic): این شرایط مرزی (مانند cyclicAMI در OpenFOAM) برای مدل‌سازی تنها یک قطاع (sector) از ماشین به جای کل آن (Full Annulus) استفاده می‌شوند. این کار هزینه‌های محاسباتی را به شدت کاهش می‌دهد. این شرایط، سطوح تناوبی یک کانال جریان را به یکدیگر متصل می‌کنند و فرض می‌کنند که جریان در تمام کانال‌های دیگر یکسان است.

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

نکات کاربردی و بهترین شیوه‌ها

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

 ترکیب‌های رایج شرایط مرزی

بر اساس تجربیات عملی و راهنماها، دو ترکیب زیر به طور گسترده برای شبیه‌سازی‌های توربوماشین استفاده می‌شوند و نتایج پایداری ارائه می‌دهند:

1. ورودی فشار کل، خروجی فشار استاتیک (Total pressure in, static pressure out): این ترکیب یک استاندارد صنعتی محسوب می‌شود و برای طیف وسیعی از شبیه‌سازی‌ها، به ویژه در مراحل اولیه طراحی، مناسب و پایدار است.
2. ورودی دبی جرمی، خروجی فشار استاتیک (Mass flow inlet, static pressure outlet): این ترکیب زمانی ایده‌آل است که دبی جرمی مشخصه اصلی عملکرد ماشین باشد و هدف، محاسبه افزایش فشار تولید شده توسط ماشین در آن دبی خاص باشد.

مدیریت خطاهای رایج: مثال پمپ

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

یک راه‌حل عملی برای این مشکل، افزایش تدریجی فشار خروجی (ramping up the outlet pressure) است. به جای اعمال مقدار نهایی فشار در تکرار صفر، می‌توان فشار خروجی را در یک بازه مشخص از تکرارها (مثلاً از تکرار ۱۰۰ تا ۲۰۰) به صورت خطی از یک مقدار اولیه (مانند صفر) به مقدار نهایی افزایش داد. این کار به میدان جریان اجازه می‌دهد تا به تدریج توسعه یافته، به فشار بالاتر عادت کرده و به یک راه حل پایدار و فیزیکی همگرا شود.

 جمع‌بندی 

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

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