Y+ چیست و چگونه آن را در شبیه‌سازی CFD تنظیم کنیم؟

مفاهیم بنیادی: لایه مرزی آشفته و پارامتر Y+

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

Y+ دقیقاً چیست؟

به زبان ساده، Y+ یک فاصله بی‌بعد است که رزولوشن شبکه را نسبت به زیرلایه‌های لایه مرزی اندازه‌میگیرد. می‌توان آن را به عنوان یک “خط‌کش مخصوص” برای ارزیابی مش در نزدیکی دیواره در نظر گرفت. این پارامتر مشخص می‌کند که اولین لایه سلول‌های ما برای کدام رویکرد مدل‌سازی توربولانس مناسب است. مقدار Y+ که از رابطه زیر به دست می‌آید، فاصله بی‌بعد از دیواره را نشان می‌دهد:

y+ = y * Uτ / ν

در این رابطه، y فاصله عمودی از دیواره، Uτ سرعت برشی (friction velocity) و ν ویسکوزیته سینماتیکی سیال است.

زیرلایه‌های لایه مرزی آشفته

لایه مرزی آشفته به طور کلی به سه زیرلایه اصلی تقسیم می‌شود که رفتارجریان در هر یک از آن‌ها منحصر به فرد است. مقدار Y+ تعیین می‌کند که اولین سلول محاسباتی در کدام یک از این نواحی قرار دارد.

اهمیت این ساختار در این است که لایه بافر (5 < y+ < 30) اغلب به عنوان “منطقه ممنوعه” (forbidden zone) برای شبکه‌بندی در نظر گرفته می‌شود. دلیل این امر آن است که نه فرض پروفیل سرعت خطیِ زیرلایه لزج و نه فرض پروفیل لگاریتمی لایه خارجی، هیچ‌کدام نمی‌توانند فیزیک پیچیده گذار در این ناحیه را به درستی توصیف کنند. این عدم تطابق، منشأ اصلی خطا در محاسبات تنش برشی دیواره و انتقال حرارت است.

در نهایت، مقدار Y+ هدف تعیین می‌کند که اولین سلول شبکه در کدام یک از این نواحی فیزیکی قرار می‌گیرد، که این تصمیم مستقیماً به دو استراتژی اصلی مدل‌سازی که در ادامه بحث می‌شود، منجر خواهد شد.

استراتژی‌های اصلی مدل‌سازی: Wall-Resolving در مقابل Wall-Modeling

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

رویکرد تفکیک کامل دیواره (Wall-Resolving)

در این رویکرد، شبکه به صورت عمدی آنقدر ریز تولید می‌شود که مرکز اولین سلول در داخل زیرلایه لزج قرار گیرد (ایده‌آل y+ < 1). این استراتژی به حلگر اجازه می‌دهد تا معادلات حاکم را مستقیماً تا سطح دیواره انتگرال‌گیری کرده و گرادیان‌های شدید سرعت و دما را با دقت بالا ثبت کند. این روش دقیق‌ترین نتایج را ارائه می‌دهد، اما هزینه محاسباتی بسیار بالایی دارد.

این رویکرد با دقت بالا در سناریوهای زیر ضروری است:

• شبیه‌سازی‌های دقیق انتقال حرارت: برای محاسبه دقیق ضریب انتقال حرارت و شار حرارتی.
• تحلیل نیروی درگ اصطکاکی (Friction Drag): برای پیش‌بینی دقیق اصطکاک پوستی که مستقیماً به گرادیان سرعت روی دیواره وابسته است.
• پیش‌بینی جدایش جریان (Flow Separation): در مسائلی مانند واماندگی (stall) ایرفویل که نقطه دقیق جدایش جریان اهمیت دارد.
• مدل‌سازی پدیده گذار (Transition): برای ثبت صحیح نقطه شروع آشفتگی (گذار از جریان آرام به آشفته).

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

رویکرد استفاده از توابع دیواره (Wall-Modeling)

این رویکرد یک جایگزین مقرون‌به‌صرفه است که در آن از یک شبکه درشت‌تر استفاده می‌شود و مرکز اولین سلول در لایه لگاریتمی قرار می‌گیرد (y+ > 30). در این روش، به جای حل‌مستقیم زیرلایه لزج و بافر، از فرمول‌های نیمه‌تجربی (توابع دیواره یا Wall Functions) برای “پل زدن” بین اولین سلول محاسباتی و دیواره استفاده می‌شود. این توابع، تنش برشی و شار حرارتی دیواره را بر اساس مقادیر متغیرها در اولین سلول تخمین می‌زنند.

این روش در شرایط زیر مناسب بوده و نتایج قابل قبولی ارائه می‌دهد:

• آیرودینامیک خارجی در مقیاس بزرگ: مثل شبیه‌سازی جریان هوا در اطراف ساختمان‌ها یا یک خودرو کامل.
• جریان‌های داخلی کاملاً توسعه‌یافته: در لوله‌ها و کانال‌های طولانی که پروفیل جریان به یک حالت پایدار رسیده است.

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

تحلیل توابع دیواره (Wall Functions) در ANSYS Fluent

برای مهندسانی که رویکرد Wall-Modeling را انتخاب می‌کنند، انتخاب صحیح تابع دیواره در Ansys Fluent برای دستیابی به نتایج معتبر بسیار حیاتی است. هر تابع دارای روش‌شناسی، محدوده Y+ ایده‌آل و کاربردهای خاص خود است.

تابع دیواره استاندارد (Standard Wall Function)

این تابع، مدل پیش‌فرض برای خانواده مدل‌های k-ε در Fluent است. مکانیسم آن بر یک رویکرد دوبخشی استوار است: اگر فاصله بی‌بعد (*y) کمتر از 11.225 باشد، از پروفیل‌سرعت خطی استفاده می‌کند و در غیر این صورت، پروفیل لگاریتمی را به کار می‌گیرد.

در این نقطه، معرفی یک تمایز فنی ضروری است: در حالی که Y+ (مبتنی بر تنش برشی دیواره، τ_w) معیار شناخته‌شده برای طراحی مش توسط کاربر است، حلگر Fluent در محاسبات داخلی خود از پارامتر مشابهی به نام Y* (مبتنی بر انرژی جنبشی توربولانسی، k) برای پایداری عددی استفاده می‌کند. دلیل این امر آن است که در نقاط جدایش یا رکود جریان، τ_w می‌تواند به صفر نزدیک شود و باعث ایجاد یک تکینگی (singularity) در محاسبه Y+ گردد، در حالی که k همواره مقداری مثبت دارد. این موضوع، زمینه را برای درک توابع دیواره پیشرفته‌تر فراهم می‌کند.

بزرگترین نقطه ضعف تابع دیواره استاندارد، نداشتن دقت در لایه بافر است. بنابراین، هنگام استفاده از آن، باید اطمینان حاصل کرد که مرکز سلول‌ها به طور قطعی خارج از محدوده 5 < y+ < 30 قرار دارد.

تابع دیواره مقیاس‌پذیر (Scalable Wall Function)

این تابع برای کاهش حساسیت به مقدار Y+ طراحی شده است و به آن “عدم حساس به وای پلاس” (y+ insensitive) نیز گفته می‌شود. منطق اصلی آن، اعمال یک محدودکننده است که از استفاده از پروفیل خطی جلوگیری می‌کند (y* = max(y*, 11.225)). این تابع همیشه محاسبات را بر اساس قانون لگاریتمی انجام می‌دهد، که آن را برای هندسه‌های پیچیده که رسیدن به Y+ ثابت در تمام سطوح دشوار است، به گزینه‌ای قوی‌تر تبدیل می‌کند.

تابع دیواره غیرتعادلی (Non-Equilibrium Wall Function)

این تابع به طور خاص برای پدیده‌های جریانی پیچیده‌تر مانند جریان‌های با گرادیان فشار بالا، جدایش جریان (flow separation) و اتصال مجدد جریان (reattachment) طراحی شده است. این تابع با در نظر گرفتن اثرات غیرتعادلی در لایه مرزی، نتایج را در چنین شرایطی نسبت به مدل استاندارد بهبود می‌بخشد. برای استفاده موثر از این تابع، شبکه باید به‌گونه‌ای طراحی شود که مرکز اولین سلول به طور قطعی در ناحیه لگاریتمی (y+ > 30) قرار گیرد.

پردازش دیواره پیشرفته (Enhanced Wall Treatment – EWT)

EWT پیشرفته‌ترین رویکرد در میان توابع دیواره فلوئنت است. این روش از یک تابع ترکیبی (blended function) استفاده می‌کند که یک پروفیل واحد و پیوسته را در تمام زیرلایه‌های لایه‌مرزی (لزج، بافر و لگاریتمی) ارائه می‌دهد. این ویژگی EWT را برای شبکه‌هایی که دارای طیف وسیعی از مقادیر Y+ هستند، بسیار مقاوم می‌کند. با این حال، توصیه رسمی Ansys این است که برای دستیابی به ماکزیمم دقت، همچنان شبکه برای یکی از دو حالت y+ ≈ 1 یا y+ > 30 طراحی شود.

برای کسب حداکثر دقت با EWT، راهنمای Ansys توصیه می‌کند که حداقل ۱۰ سلول در ناحیه متأثر از لزجت (تقریباً y+ < 40) قرار داده شود تا پروفیل‌های سرعت میانگین و کمیت‌های توربولانسی به درستی تفکیک شوند.

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

راهنمای انتخاب: تطبیق مدل توربولانس با استراتژی نزدیک دیواره

ارتباط بین مدل توربولانس RANS و استراتژی پردازش نزدیک دیواره، یک پیوند حیاتی است. عدم تطابق این دو، مانند “ریختن بنزین در موتور دیزل” است و قطعاً به نتایج نادرست منجر خواهد شد. هر مدل توربولانس بر اساس فرضیات خاصی توسعه یافته و برای عملکرد صحیح، به نوع خاصی از اطلاعات از ناحیه نزدیک دیواره نیاز دارد.

جدول زیر یک راهنمای جامع برای انتخاب ترکیب صحیح مدل توربولانس و محدوده Y+ در ANSYS Fluent ارائه می‌دهد:

با تثبیت این چارچوب نظری، اکنون می‌توانیم به مراحل عملی مورد نیاز برای پیاده‌سازی این بهترین شیوه‌ها بپردازیم.

گردش کار عملی: از محاسبه تا اعتبارسنجی

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

گام اول: تخمین ارتفاع اولین لایه مش (First Layer Thickness)

قبل از شروع فرآیند شبکه‌بندی، باید ارتفاع مورد نیاز برای اولین سلول مجاور دیواره (Δy) را برای رسیدن به Y+ هدف خود تخمین بزنید. این فرآیند معمولاً شامل مراحل زیر است:

• • محاسبه عدد رینولدز (Re)
  • تخمین ضریب اصطکاک پوستی (Cf)
  • محاسبه تنش برشی دیواره (τw)
  • محاسبه سرعت برشی (Uτ)
  • محاسبه نهایی ارتفاع لایه اول (Δy) ابزارهای آنلاین متعددی (Y+ Calculator) برای ساده‌سازی این محاسبات وجود دارند.

گام دوم: تنظیم لایه‌های مرزی (Inflation) در نرم‌افزار مش‌بندی

پس از محاسبه Δy، این مقدار را در ابزار Inflation در ANSYS Meshing وارد کنید. تنظیمات کلیدی عبارتند از:

• • انتخاب متد (Method): گزینه First Layer Thickness را انتخاب کنید.
  • ارتفاع لایه اول (First Layer Height): مقدار Δy محاسبه شده را در این قسمت وارد نمایید.
  • تعداد لایه‌ها (Number of Layers): برای اطمینان از یک گذار نرم به شبکه اصلی، مقداری بین 10 تا 20 لایه توصیه می‌شود.
  • نرخ رشد (Growth Rate): مقدار 1.2 یک انتخاب استاندارد و مناسب برای ایجاد یک گذار هموار است.

گام سوم: بازرسی و اعتبارسنجی مقادیر Y+ پس از حل

پس از همگرایی شبیه‌سازی، مرحله حیاتی اعتبارسنجی Y+ آغاز می‌شود. برای مصورسازی توزیع Y+ در ANSYS Fluent:

1. 1. به مسیر Graphics > Contours بروید.
   2. از منوی کشویی Contours of، گزینه Turbulence... و سپس Wall Yplus را انتخاب کنید.
   3. در بخش Surfaces، تمام سطوح دیواره مدل را انتخاب کرده و بر روی Display کلیک کنید.

نکته کلیدی:

برای وضوح بیشتر، تیک گزینه Node Values را بردارید تا فقط کانتور مبتنی بر سلول نمایش داده شود.

1. در کانتور نمایش داده شده، به دنبال لکه‌های بزرگ رنگی در “منطقه ممنوعه” (5 < Y+ < 30) یا یک الگوی “پوست گوری” (zebra stripe pattern) (که نشان‌دهنده تغییرات ناگهانی و شدید در اندازه سلول‌های لایه مرزی است) باشید.

اشتباهات رایج که باید از آنها اجتناب کرد (Common Mistakes to Avoid)

• خطای مرگبار (The Fatal Flaw): هرگز به مقدار میانگین Y+ که توسط نرم‌افزار گزارش می‌شود، به تنهایی اعتماد نکنید. یک مقدار میانگین 30 می‌تواند حاصل ترکیبی از مقادیر نامناسب (مانند 1 و 59) باشد. همیشه کانتور توزیع Y+ را به صورت بصری بررسی کنید.
• تله شماره دو (The Second Trap): محاسبات اولیه بر اساس فرضیات ساده است. نواحی با انحنای زیاد یا گوشه‌های تیز نیازمند شبکه محلی ریزتری برای رسیدن به Y+ هدف هستند. این نواحی را نادیده نگیرید.

۷. نکته پیشرفته: پیامدهای عددی تفاوت بین Y+ و Y*

همانطور که در بخش ۴ اشاره شد، حلگرهای مدرن CFD مانند Fluent به دلیل پایداری عددی، از Y* (مبتنی بر انرژی جنبشی k) به جای Y+ (مبتنی بر تنش برشی τ_w) در محاسبات داخلی خود استفاده می‌کنند. این بخش به پیامدهای این انتخاب می‌پردازد.

رویکرد مبتنی بر Y+ ذاتاً تکراری (iterative) است. برای محاسبه Y+، به τ_w نیاز داریم و برای محاسبه τ_w با استفاده از توابع دیواره، به Y+ نیازمندیم. این وابستگی متقابل باید برای هر سلول مرزی در هر تکرار حل شود که از نظر محاسباتی پرهزینه است.

در مقابل، رویکرد Y* مستقیم (direct) است. حلگر مقدار k را از تکرار قبلی در اختیار دارد، با استفاده از آن Y* را مستقیماً محاسبه می‌کند و سپس τ_w را به دست می‌آورد. این فرآیند بدون تکرار، سریع‌تر و در جریان‌های دارای جدایش یا رکود (که در آن‌ها τ_w به صفر میل می‌کند) بسیار پایدارتر است.

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

۸. جمع‌بندی و توصیه‌های نهایی

این راهنما مفاهیم کلیدی مرتبط با Y+، دو استراتژی اصلی مدل‌سازی نزدیک دیواره (Wall-Resolving در مقابل Wall-Modeling)، و گردش کار عملی برای تنظیم Y+ در فلوئنت را پوشش داد. با درک این اصول، مهندسان می‌توانند تصمیمات آگاهانه‌تری بگیرند و از اعتبار نتایج شبیه‌سازی خود اطمینان حاصل کنند.

توصیه نهایی و کلیدی یک اصل غیرقابل تخطی است: همیشه و تحت هر شرایطی، کانتور توزیع Y+ را بر روی تمام سطوح دیواره به صورت بصری بازرسی کنید. اکتفا به مقدار میانگین گزارش‌شده توسط حلگر، یک خطای مهندسی فاحش است که می‌تواند اعتبار کل تحلیل شما را زیر سوال ببرد. این بازرسی بصری، مهم‌ترین و موثرترین سپر دفاعی شما در برابر خطاهای بحرانی مدل‌سازی نزدیک دیواره است.