راهنمای کامل توابع دیواره (Wall Functions) در فلوئنت: تفاوت‌ها و کاربردها

درک لایه مرزی آشفته: اساس کار توابع دیواره

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

ساختار لایه مرزی آشفته به طور کلاسیک به سه زیرلایه اصلی تقسیم می‌شود:

• زیرلایه لزج (Viscous Sublayer): این ناحیه که بی‌واسطه به دیواره چسبیده است، تحت سلطه کامل اثرات لزجت مولکولی قرار دارد. در این لایه، پروفیل سرعت تقریباً خطی است و با رابطه u+ = y+ توصیف می‌شود. این ناحیه معمولاً در محدوده y+ < 5 قرار دارد.
• لایه بافر (Buffer Layer): این لایه یک منطقه گذار است که در آن، اثرات لزجت مولکولی و اغتشاشات توربولانسی هر دو به یک اندازه اهمیت دارند. هیچ‌یک از روابط ساده خطی یا لگاریتمی به تنهایی قادر به توصیف دقیق پروفیل سرعت در این ناحیه نیستند. محدوده این لایه تقریباً 5 < +y+ < 30 & y است.
• لایه لگاریتمی (Log-law Layer): در این ناحیه که از دیواره فاصله بیشتری دارد، اغتشاشات توربولانسی کاملاً غالب هستند. پروفیل سرعت در اینجا از یک رابطه لگاریتمی معروف به «قانون دیواره» (Law of the Wall) پیروی می‌کند. محدوده رایج این ناحیه 30 < +y+ < 300 & y در نظر گرفته می‌شود، اما باید توجه داشت که حد بالای آن ثابت نیست و به شدت به عدد رینولدز بستگی دارد. در جریان‌های با عدد رینولدز بالا، این لایه می‌تواند تا y+ > 2000 نیز گسترش یابد.

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

y+ = y*Uτ/ν

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

دو رویکرد اصلی برای مدل‌سازی نزدیک دیواره (Near Wall Treatment)

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

 رویکرد حل کامل لایه مرزی (Wall-Resolving Approach)

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

الزامات: برای پیاده‌سازی صحیح این روش، لازم است که مقدار میانگین +y برای سلول‌های مجاور دیواره نزدیک به ۱ باشد. همچنین، باید حداقل ۵ تا ۱۰ لایه سلول در داخل زیرلایه لزج (یعنی در محدوده y+ < 5) قرار گیرند تا گرادیان‌های شدید سرعت به درستی ثبت شوند. البته باید توجه داشت که ریز کردن بیش از حد مش و رسیدن به مقادیر y+ < 0.1 معمولاً بهبود چشمگیری در نتایج ایجاد نمی‌کند و چنین مقادیری اغلب تنها برای مدل‌های گذار به آشفتگی (Transition Models) ضروری هستند.

معایب: عیب اصلی این رویکرد، هزینه محاسباتی بسیار بالا (computationally expensive) است. ایجاد یک شبکه بسیار ریز در نزدیکی دیواره منجر به افزایش چشمگیر تعداد کل سلول‌های محاسباتی می‌شود. برای مثال، در شبیه‌سازی‌های نمونه بر روی اشکال مهندسی، تغییر رویکرد از مدل‌سازی دیواره به حل کامل لایه مرزی، تعداد سلول‌ها را تقریباً دو برابر کرده است:

• • بدنه هواپیما: از حدود ۶۵۰ هزار به ۱.۱۶ میلیون سلول.
  • زیردریایی: از حدود ۶.۶ میلیون به ۱۱.۱ میلیون سلول.
  • خودروی فرمول یک: از حدود ۵۷.۸ میلیون به ۱۱۱.۱ میلیون سلول.

رویکرد مدل‌سازی با توابع دیواره (Wall-Modeling/Wall Function Approach)

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

• مزایا: مزیت اصلی این روش، کاهش چشمگیر هزینه محاسباتی است. با استفاده از شبکه درشت‌تر و قرار دادن اولین سلول در ناحیه لگاریتمی (y+ > 30)، دیگر نیازی به حل مستقیم زیرلایه لزج و لایه بافر نیست و تعداد کل سلول‌ها به شدت کاهش می‌یابد.
• محدودیت‌ها: دقت این توابع بر پایه فرضیاتی استوار است که در همه جریان‌ها صادق نیست. اعتبار فیزیکی و دقت توابع دیواره در شرایطی مانند جریان‌های دارای جدایش (separated flows)، گرادیان فشار نامطلوب شدید (strong adverse pressure gradient)، و جریان‌های با عدد رینولدز پایین، به طور جدی زیر سوال است.

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

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

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

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

این تابع که بر اساس کار پیشگامانه Launder و Spalding بنا شده، پایه‌ای‌ترین مدل در میان توابع دیواره است.

• مکانیسم عملکرد: این تابع از یک شرط دوگانه استفاده می‌کند. نرم‌افزار ابتدا مقدار *y را محاسبه می‌کند. اگر y* < 11.225 باشد، پروفیل سرعت بر اساس رابطه خطی زیرلایه لزج (u* = y*) محاسبه می‌شود. اگر y* > 11.225 باشد، از رابطه لگاریتمی استفاده می‌شود.
• تفاوت کلیدی +y و *y: اگرچه در مباحث تئوری از پارامتر +y که بر اساس تنش برشی دیواره (τ) تعریف می‌شود استفاده می‌کنیم، نرم‌افزارهای CFD در عمل از پارامتر *y که بر اساس انرژی جنبشی آشفتگی (k) تعریف می‌شود، بهره می‌برند. دلیل این انتخاب، افزایش پایداری عددی است. در نقاط جدایش جریان، تنش برشی دیواره به صفر میل می‌کند که این امر منجر به ایجاد یک تکینگی (singularity) در تعریف +y می‌شود. استفاده از *y که به k وابسته است، این مشکل را برطرف می‌کند. این دو پارامتر تنها در شرایط ایده‌آل و تعادلی با یکدیگر برابر هستند.
• الزامات y+ و محدودیت‌ها: محدودیت اصلی این تابع در ناحیه بافر (5 < +y+ < 30 & y) خود را نشان می‌دهد. اگر اولین گره محاسباتی در این ناحیه قرار گیرد، هیچ‌یک از دو رابطه خطی یا لگاریتمی دقت کافی را ندارند و نتایج به شدت غیرقابل اعتماد خواهند بود. بنابراین، هنگام استفاده از این تابع، شبکه باید به گونه‌ای طراحی شود که y+ سلول‌های مجاور دیواره یا کاملاً در زیرلایه لزج (y+ < 5) یا کاملاً در ناحیه لگاریتمی (y+ > 30) قرار گیرند.
• کاربرد مناسب: این تابع برای جریان‌های برشی ساده و کاملاً توسعه‌یافته که به دیواره چسبیده‌اند (attached flows) و گرادیان فشار ندارند، مناسب است. به دلیل حساسیت بالا به y+، استفاده از آن در مسائل پیچیده امروزی کمتر توصیه می‌شود.

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

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

• مکانیسم عملکرد: این تابع با اعمال یک محدودیت پایین برای مقدار *y، حلگر را مجبور به استفاده از قانون لگاریتمی می‌کند. این کار با رابطه y* = max(y*, 11.225) انجام می‌شود. این مکانیزم هوشمندانه از بروز خطای فاحش جلوگیری می‌کند؛ زیرا اگر شبکه در ناحیه‌ای به صورت ناخواسته بیش از حد ریز شود و y+ به زیر ۱۱.۲۲۵ سقوط کند، تابع استاندارد به اشتباه از فرمول زیرلایه لزج استفاده می‌کند که منجر به نتایج نادرست می‌شود (مشکلی که در نمودار حساسیت به شبکه توابع استاندارد به وضوح دیده می‌شود).
• مزایا و کاربرد: این تابع نسبت به ریز شدن ناخواسته مش حساسیت کمتری دارد (y+ insensitive). این ویژگی آن را به یک انتخاب عمل‌گرایانه و قوی در شبیه‌سازی‌های صنعتی با هندسه‌های پیچیده تبدیل می‌کند، جایی که حفظ y+ > 30 در تمام سطوح عملاً غیرممکن است و این تابع از شکست شبیه‌سازی یا بروز خطاهای بزرگ جلوگیری می‌کند.

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

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

• مکانیسم عملکرد: این تابع اثرات گرادیان فشار و تغییرات شدید در جریان را در محاسبات خود لحاظ می‌کند و پروفیل سرعت را متناسب با این شرایط اصلاح می‌نماید.
• الزامات y+: همانند تابع دیواره استاندارد، این تابع نیز برای عملکرد صحیح نیازمند قرارگیری اولین گره محاسباتی در ناحیه لگاریتمی است (30 < +y+ < 300 & y).
• کاربرد مناسب: این تابع به طور خاص برای جریان‌های دارای گرادیان فشار شدید، جدایش و اتصال مجدد جریان (separation and reattachment) و جریان‌های برخوردی (impinging flows) توصیه می‌شود و نسبت به تابع استاندارد، پیش‌بینی دقیق‌تری در این شرایط ارائه می‌دهد.

تابع دیواره بهبودیافته (Enhanced Wall Treatment – EWT)

این گزینه که در فلوئنت با نام Enhanced Wall Treatment شناخته می‌شود، پیشرفته‌ترین و انعطاف‌پذیرترین رویکرد در میان توابع دیواره است.

• مکانیسم عملکرد: EWT به جای استفاده از یک شرط دوگانه، از یک تابع ترکیبی (blending function) استفاده می‌کند تا یک پروفیل واحد، پیوسته و از نظر فیزیکی سازگار برای کل ناحیه نزدیک دیواره—شامل زیرلایه لزج، بافر و لگاریتمی—ایجاد کند. این رویکرد، گذار نرمی بین لایه‌ها فراهم کرده و به طور خاص، مدل‌سازی ناحیه چالش‌برانگیز بافر را که سایر توابع در آن شکست می‌خورند، به شکل قابل توجهی بهبود می‌بخشد.
• مزایا: مزیت کلیدی EWT حساس نبودن آن به مقدار y+ است. این تابع می‌تواند هم برای مش‌های بسیار ریز که زیرلایه لزج را حل می‌کنند (y+ ~ 1) و هم برای مش‌های درشت که از قانون لگاریتمی استفاده می‌کنند  (y+ > 30) نتایج دقیقی ارائه دهد.
• کاربرد مناسب: EWT تقریباً برای تمام انواع جریان‌ها گزینه برتری است، به ویژه زمانی که محدوده وسیعی از y+ در دامنه محاسباتی وجود دارد یا پدیده‌هایی مانند انتقال حرارت دقیق، جدایش جریان و جریان‌های با عدد رینولدز پایین اهمیت بالایی دارند.

جمع‌بندی و توصیه‌های کلیدی

انتخاب صحیح روش مدل‌سازی نزدیک دیواره و «توابع دیواره در فلوئنت» یکی از مهم‌ترین تصمیمات در فرآیند شبیه‌سازی CFD است که مستقیماً بر دقت نتایج و هزینه محاسباتی تأثیر می‌گذارد. در ادامه، مهم‌ترین نکات این راهنما به صورت توصیه‌های کاربردی مرور می‌شوند.

• اصل موازنه را درک کنید: انتخاب بین رویکرد حل کامل لایه مرزی (y+ ~ 1) و استفاده از توابع دیواره (y+ > 30) یک تصمیم راهبردی بین حداکثر دقت و هزینه محاسباتی قابل قبول است.
• از ‘منطقه ممنوعه’ دوری کنید: از قرار دادن اولین گره محاسباتی در لایه بافر (‘منطقه ممنوعه’) (5 < +y+ < 30 & y)، به خصوص هنگام استفاده از توابع دیواره استاندارد، جداً خودداری کنید. مدل‌سازی در این ناحیه با خطای قابل توجهی همراه است.
• تابع دیواره مناسب را برگزینید: برای جریان‌های پیچیده با جدایش، Non-Equilibrium (با y+ > 30) گزینه مناسبی است. برای دستیابی به حداکثر انعطاف‌پذیری و دقت در گستره‌های مختلف y+، تابع Enhanced Wall Treatment (EWT) برترین گزینه محسوب می‌شود.
• دقت را در اولویت قرار دهید: برای مسائلی که به پارامترهای دیواره بسیار حساس هستند، مانند محاسبه دقیق انتقال حرارت یا نیروی درگ، رویکرد حل کامل لایه مرزی (y+ ~ 1) همیشه ارجحیت دارد و نتایج قابل اعتمادتری ارائه می‌دهد.

توصیه نهایی: همواره مقدار y+ را به صورت یک کانتور بر روی سطوح یا یک نمودار توزیع بررسی کنید. هرگز به مقدار میانگین y+ اکتفا نکنید، زیرا این مقدار می‌تواند نواحی بحرانی با مقادیر نامناسب y+ (بسیار بالا، بسیار پایین یا در ناحیه بافر) را از دید شما پنهان کند و منجر به تحلیل نادرست نتایج شود.