آشنایی با معادله‌ی ناویر استوکس

زمان مطالعه: 4 دقیقه

در دینامیک سیالات، معادلات ناویر-استوکس Navier-Stokes معادلاتی هستند که حرکت سه بعدی سیال ویسکوز یا لزج را توصیف می­کنند. این معادلات به نام کلود لوئیس ناویر (۱۷۸۵-۱۸۳۶) و جورج گابریل استوکس (۱۸۱۹-۱۹۰۳) نام گذاری شده‌­اند. در شرایطی که گرادیان­‌های دمایی بزرگ در سیال وجود نداشته باشد، این معادلات تقریب بسیار خوبی از واقعیت ارائه می­دهند.

معادلات ناویر-استوکس شامل یک معادله پیوستگی وابسته به زمان برای بقای جرم، سه معادله‌ی بقای مومنتوم وابسته به زمان و یک معادله بقای انرژی وابسته به زمان است. در این دسته معادلات چهار متغیر مستقل وجود دارد، x، y و z مختصات فضایی دامنه و t متغیر زمان می­‌باشد.

آشنایی با معادله‌ی ناویر استوکس- مرجع تخصصی آموزش نرم افزار انسیس و فلوئنت
شکل عمومی معادلات ناویر استوکس

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

حل معادلات ناویر-استوکس

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

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

ویژگی­‌های جریان آشفته

  • جریان آشفته تمایل دارد در سرعت های بالاتر، ویسکوزیته کم و در ابعاد خطی مشخصه بالاتر رخ دهد.
  • اگر عدد رینولدز بیشتر از Re > 3500 باشد، جریان آشفته است.
  • بی نظمی: این جریان به واسطه حرکات نامنظم ذرات سیال شناسایی می­شود. حرکت ذرات سیال به شدت نامنظم است. چون در حال حاضر جواب قطعی برای جریان متلاطم در دسترس نیست، این پدیده به صورت آماری بررسی می‌شود.
  • پخش یا انتشار: در جریان آشفته، یک توزیع سرعت نسبتاً یکسان در سراسر بخش لوله وجود دارد، در نتیجه کل سیال در یک مقدار مشخص جریان می‌یابد و در نزدیکی دیواره‌ها این سرعت شدیداً کاهش پیدا می‌کند. مشخصه­‌ای که باعث افزایش اختلاط، نرخ انتقال جرم، مومنتوم و انرژی در یک جریان است، “پخش (Diffusivity)” نامیده می­شود.
  • چرخش: جریان آشفته با مکانیزم تولید گردابه سه بعدی قوی مشخص می­‌شود. این مکانیزم به کشش گردابی معروف است.
  • اتلاف یا اضمحلال: فرآیند اتلاف فرآیندی است که در آن انرژی جنبشی جریان آشفته توسط تنش برشی ویسکوز به انرژی داخلی تبدیل می‌شود.

مقیاسِ میکرو کولموگروف (Kolmogorov Microscales)

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

به عنوان مثال، یک ابر کومولوس را در نظر بگیرید. در مقیاس ماکرو ابر می­تواند در ابعاد کیلومتر باشد و ممکن است در مدت زمان طولانی رشد کند یا ثابت باقی بماند. در داخل ابر، گردابه‌ها ممکن است در مقیاس‌هایی از مرتبه میلی‌متر وجود داشته باشند. برای جریان­‌های کوچکتر مانند جریان داخل لوله­‌ها، مقیاس‌­های میکرو ممکن است بسیار کوچکتر باشند. بخش اعظم انرژی جنبشی جریان آشفته در مقیاس­‌های ماکرو موجود است. انرژی با استفاده از مکانیسم اینرسی از مقیاس ماکرو به مقیاس میکرو منتقل می‌­شود. این فرآیند به عنوان آبشار انرژی آشفته (turbulent energy cascade) شناخته می­شود.

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

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

آشنایی با معادله‌ی ناویر استوکس- مرجع تخصصی آموزش نرم افزار انسیس و فلوئنت

که در آن ε میانگین نرخ اتلاف انرژی جنبشی آشفته در واحد جرم است و واحد آن (m2/s3) است. ν ویسکوزیته سینماتیکی سیال است و واحد آن (m2/s ) است.

اندازه کوچکترین گردابه در جریان توسط ویسکوزیته تعیین می­شود. مقیاس طولی کلموگروف با کاهش ویسکوزیته کاهش می­‌یابد. برای جریان­‌های با عدد رینولدز بسیار بالا، نیروهای ویسکوز نسبت به نیروهای اینرسی کوچکتر هستند. الزاماً تا زمانی که اثرات ویسکوزیته چشم­گیر باشند و اتلاف انرژی صورت بگیرد، شاهد حرکات در مقیاس­های کوچک‌تر خواهیم بود. نسبت بزرگ‌ترین مقیاس طولی به کوچک‌ترین مقیاس طولی در جریان آشفته با عدد رینولدز (که با توان سه چهارم افزایش می‌یابد) متناسب است:

آشنایی با معادله‌ی ناویر استوکس- مرجع تخصصی آموزش نرم افزار انسیس و فلوئنت

این امر باعث می­‌شود که شبیه سازی عددی مستقیم جریان آشفته عملاً غیرممکن باشد. برای مثال، جریانی با عدد رینولدز ۱e+6 را در نظر بگیرید. در این حالت نسبت L/l متناسب با ۱۰ به توان ۱۸/۴ است. از آنجایی که مسئله باید به صورت سه بعدی تحلیل گردد، باید از شبکه­‌ای که حداقل از ۱e+14  نقطه تشکیل شده است استفاده کنیم. ایجاد شبکه و حل مسئله در این حالت بسیار فراتر از ظرفیت و امکانات رایانه‌­های موجود است.

با دوره اموزش انسیس فلوئنت به یادگیری کار با نرم افزار فلوئنت بپردازید.

دیدگاه‌ خود را بنویسید

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *

سبد خرید
error: محتوا محافظت شده است
کد تخفیف Flu_Student تا 16 آذر، روز دانشجو، فعال شد.
+ +
پیمایش به بالا