خلاصه جامع پایان‌نامه

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

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

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

دوم، مدل دیواره با روش شبیه‌سازی گردابه‌های جداشده (DES) در چارچوب ترکیبی RANS/LES به کار گرفته شده است. در این حالت، مدل اجازه می‌دهد تا حل جریان بر روی شبکه‌هایی که در جهت عمود بر دیواره حدود دو مرتبه بزرگی درشت‌تر از شبکه‌های معمول DES هستند، به‌صورت قابل‌اعتماد انجام شود. این بدان معناست که هزینه محاسباتی DES به‌طور قابل‌توجهی کاهش می‌یابد، بدون آنکه دقت حل در نواحی دور از دیواره فدا شود. نتایج نشان می‌دهد که پروفیل‌های سرعت و تنش‌های برشی مدل‌شده در این روش، حساسیت بسیار کمتری به اندازه شبکه در نزدیکی دیواره دارند. علاوه بر این، بررسی‌ها نشان می‌دهد که حل‌های DES و WMDES به شرایط اولیه حساس هستند و انتخاب شرایط اولیه نامناسب می‌تواند منجر به همگرایی به حل‌های یک‌بعدی و غیرفیزیکی شود. این یافته اهمیت توجه به شرایط اولیه در شبیه‌سازی‌های hybrid RANS/LES را برجسته می‌کند.

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

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

با این حال، محدودیت‌هایی نیز وجود دارد؛ از جمله اینکه مدل بر فرض جریان تعادلی و اعداد رینولدز بالا بنا شده است و در اعداد رینولدز پایین‌تر یا در نواحی با گرادیان فشار شدید، ممکن است خطای مدل‌سازی قابل‌توجهی ایجاد شود. همچنین، در نزدیکی ناپیوستگی‌های شیب هندسی مانند لبه‌های تیز، رفتار غیرفیزیکی در فشار و چگالی مشاهده شده است که نیازمند اصلاحات آینده است. پژوهش‌های آینده می‌توانند بر توسعه اصلاحات برای اعداد رینولدز پایین، وارد کردن اثرات گرادیان فشار در مدل، بهبود پیش‌بینی ضخامت جابجایی، توسعه سپر لایه مرزی برای DES و تعیین شبکه بهینه برای WMLES متمرکز شوند. این جهت‌گیری‌ها می‌توانند دامنه کاربرد مدل را گسترده‌تر کرده و دقت آن را در شرایط پیچیده‌تر جریان بهبود بخشند.

وقتی دیواره دیگر مانع نیست: انقلابی در شبیه‌سازی جریان‌های آشفته

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

چرا نزدیک دیواره این‌قدر دردسر دارد؟

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

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

ایده‌ای که همه چیز را ساده می‌کند

نوآوری اصلی این پایان‌نامه در یک ایده ساده اما قدرتمند نهفته است: به جای اینکه سعی کنیم جریان پیچیده نزدیک دیواره را با تمام جزئیاتش حل کنیم، بیایید آن را عمداً ساده کنیم. نویسنده پایان‌نامه، خانم کارمن-یوانا اورساچی، با اصلاح مدل آشفتگی معروف اسپالارت-آلماراس، کاری کرد که در نزدیکی دیواره، سرعت جریان به‌جای تغییرات سریع و پیچیده، به‌صورت خطی و ساده تغییر کند. همزمان، ویسکوزیته گردابه‌ای – که معیاری از میزان اختلاط در جریان است – به یک مقدار ثابت می‌رسد. اما نکته کلیدی اینجاست که تنش برشی کل، همان مقدار واقعی خود را حفظ می‌کند.

این یعنی چه؟ یعنی می‌توانیم جریان را در نزدیکی دیواره با شبکه‌ای بسیار درشت‌تر حل کنیم، بدون اینکه دقت پیش‌بینی نیروهای آیرودینامیکی مثل برآ و پسا از دست برود. این دقیقاً همان چیزی است که مهندسان سال‌ها به دنبالش بودند.

پروفسور دیوید دارموفال، استاد دانشگاه MIT و از چهره‌های شاخص این حوزه، در مقدمه این پژوهش تأکید می‌کند که «روش‌های دقیق شبیه‌سازی جریان، اگرچه از نظر فیزیکی واقع‌گرایانه‌تر هستند، اما هزینه آن‌ها همچنان یک مانع جدی برای استفاده مؤثر از آن‌هاست.» این جمله به‌خوبی نشان می‌دهد که چرا چنین نوآوری‌هایی تا این حد حیاتی هستند.

سه میدان نبرد: RANS، DES و LES

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

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

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

چرا این کار مهم است؟

اهمیت این پژوهش را می‌توان در یک جمله خلاصه کرد: کاهش هزینه محاسباتی بدون فدا کردن دقت. در صنعت هوافضا، هر ساعت شبیه‌سازی روی ابررایانه‌ها هزینه‌های هنگفتی دارد. اگر بتوان با شبکه‌ای درشت‌تر به همان دقت رسید، یعنی می‌توان شبیه‌سازی‌های بیشتری انجام داد، طرح‌های بیشتری را بررسی کرد و در نهایت، هواپیماهای بهینه‌تری طراحی کرد. این دقیقاً همان فلسفه‌ای است که بنیاد ملی علوم آمریکا در گزارش «چشم‌انداز CFD 2030» بر آن تأکید کرده است.

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

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

آینده‌ای که پیش روست

نویسنده پایان‌نامه در بخش کارهای آینده، چند مسیر را پیشنهاد می‌کند: توسعه اصلاحات برای اعداد رینولدز پایین، وارد کردن اثرات گرادیان فشار در مدل، بهبود پیش‌بینی ضخامت جابجایی، توسعه سپر لایه مرزی برای DES و تعیین شبکه بهینه برای WMLES. هر یک از این مسیرها می‌تواند دامنه کاربرد این مدل را گسترده‌تر کند.

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