حل عددی معادلات ناویر استوکس در حالت پایای تراکم‌ناپذیر آشفته با استفاده از روش تابع پایه شعاعی چند ربعی

در روش‌های عددی هزینه و انرژی قابل توجهی صرف ایجاد و در مراحل بعدی اعمال تغییرات لازم در شبکه می‌شود. به همین دلیل روش‌های بدون شبکه به سرعت در حال توسعه و به کارگیری در مسائل فیزیکی و مهندسی هستند. یکی از انواع این روش‌ها، روشهای تابع پایه شعاعی هستند که روش چندربعی یکی از توانمندترین آنهاست. در این پژوهش، روش تابع پایه شعاعی چندربعی برای حل معادلات تراکم‌ناپدیر جریان پایا شامل معادلات پیوستگی، ناویراستوکس و مدل آشفتگی استاندارد، در یک میدان دوبعدی مورد ارزیابی قرار گرفته‌است. این میدان شامل یک هندسهِ حفره با درپوش متحرک مربعی، به ابعاد m0.5×m0.5 می‌باشد که در پنج عدد رینولدز 10^5×2.5، 10^5×5، 10^6×1، 10^6×2 و 10^6×5.5 تحلیل گردیده است. دامنه مذکور دو بار با تعداد نقاط داخلی 36 و 121مورد حل قرار گرفته و متغیر‌های مختلف محاسبه شده‌اند. طی این فرآیند، دو کمیت مهم متغیر شکل c بهینه و مجموعه ضرایب λ بهینه برای هر میدان جریان، مورد بحث و بررسی قرار گرفته است. نتایج نشان می‌دهد، با اتخاذ فرض استقلال مقادیر پارامتر c، برای می‌توان به یک الگوی قابل پیش‌بینی برای مجموعه λ دست‌یافت. الگوهای مذکور به‌همراه توابع پیش‌بین میدان‌های جریان مورد نظر، با نتایج روش حجم محدود (نرم افرار انسیس فلوئنت)، مورد مقایسه قرار گرفتند. ضرایب نشساتکلیف 93 الی 99 درصدِ و بیشینه خطای جذر میانگین مربعات نسبی در حد یک درصد، بدست آمده از این مقایسه برای پنج متغیر مستقل نشان‌دهنده قابل اعتماد بودن ترکیب فرض مذکور، الگوهای تکرارپذیر مجموعه λ و نیز توابع پیش‌بینی‌کننده آنها، می‌باشند.

Numerical Solution of Steady Incompressible Turbulent Navier–Stokes Equations using Multiquadric Radial Basis Function (MQ-RBF) Method

The inconveniences of introducing and modifying the mesh grids in meshbased numerical methods lead the researchers to meshfree methods, among which the RBF methods are probably the most interesting and powerful ones. In this research, the numerical solution of the steadystate incompressible continuity and Navier–Stokes equations, and the standard kƐ turbulence model was investigated in a 2D domain. The computational domain consisting of a 0.5 m×0.5 m square liddriven cavity was analyzed for five Reynolds numbers of 2.5×10^5, 5×10^5, 10×10^5, 2×10^6, and 5.5×10^6. The Multiquadric Radial Basis Function (MQRBF), as the most successful RBF, was employed with 36 and 121 domain computational nodes to solve the PDEs. The velocity fields in two directions, the static pressure, the turbulent kinetic energy and the turbulent energy dissipation, were computed. A try–and–error algorithm was used for solving a set of nonlinear equations, and the optimal values of the shape parameter c and the λ set coefficients were evaluated and discussed for each flow field. According to the results, assuming the independence of the values of the shape parameter c for each flow field at different Reynolds numbers, a predictable pattern can be obtained for the λ set for different Reynolds’ numbers in the studied range. These patterns with the predictor functions of the flow fields were compared to existing benchmark results of the finite volume method (ANSYS Fluent). The NashSutcliffe coefficients of 9399% and RRSME of about %1 obtained from this comparison indicated the reasonable accuracy of the assumption concerning the independence of the shape parameter c of the Reynolds’ numbers, the repeatable patterns of the normalized λ set, and polynomial predictor functions in the MQRBF method for each flow field.