حل عددي مدل‌هاي دوسيالي تراكم‌پذير هم‌دما با استفاده از روش‌هاي بقايي در امتداد مسير مركزي

Numerical solution of compressible, isothermal two-fluid models using path-conservative central schemes

دقيق‌ترين و پيچيده‌ترين مدل‌هاي موجود براي تحليل جريآنهاي دوفازي، مدل‌هاي دوسيالي مي‌باشند. براي جريآنهاي دوفازي هم‌دماي تراكم‌پذير، دو مدل تك‌فشاري و دو فشاري موجود مي‌باشد. علي‌رغم قابليت بالاي اين مدل‌ها در تحليل جريان، به دليل حضور عبارات غيربقايي در معادلات مومنتم فازها نمي‌توان آنها را در قالب بقايي بيان كرد. بنابراين، شرط كلاسيك رنكين-هوگونيوت در محل ناپيوستگي‌هاي ميدان جريان، براي اين‌گونه معادلات قابل اعمال نيست و اعمال روشهاي عددي كلاسيك براي حل اين معادلات غالباً با مشكل مواجه بوده است. در مقاله حاضر براي غلبه بر اين مشكل، از شيوه جديد انتگرال‌گيري در امتداد مسير استفاده شده‌ است. در اين شيوه مي‌توان شرط عمومي رنكين-هوگونيوت را در امتداد مسير واصل ميان حالت‌هاي سمت چپ و راست ناپيوستگي اعمال كرد. پس از ارايه شكل بقايي در امتداد مسير روش‌هاي عددي مركزي لكس-فردريكس، لكس-وندروف و روزانوف، مساله شير آب و لوله ضربه با سرعت نسبي بزرگ با استفاده از اين روش‌ها حل شده‌اند و با بكارگيري شبكه‌هاي مختلف، استقلال حل از شبكه حاصل شد. مقايسه نتايج براي شير آب نشان مي‌دهد كه تطابق خوبي با حل تحليلي وجود دارد و نتايج براي لوله ضربه نيز نشان‌دهنده توانايي بالاي روش بكار رفته براي تسخير ناپيوستگي هاي ميدان جريان در مسايل دوفازي است.

Two-fluid models are the most accurate and complex models for analysis of two-phase flows. There are two different two-fluid models for analyzing compressible isothermal two-phase flows which are Single Pressure Model (SPM) and Two-Pressure Model (TPM). In spite of capabilities of these models in capturing two-phase flow behavior, it is not possible to express them in conservative form due to existence of non-conservative term in momentum equation of phases. Therefore, the classical Rankine-Hugoniot condition across discontinuities in the flow filed is not applicable for these equations and there would be difficulty in using classical numerical methods for solving these equations. In this paper a new path-conservative method is used to overcome this difficulty. In this method, one can apply general Rankine-Hugoniot condition along a path connecting left and right states of the discontinuity. After expressing path-conservative form of the employed central numerical methods which are Lax-Fridriches, Lax-Wendroff and Rusanove, water faucet and large relative velocity shock tube problems are solved by using these schemes. Grid independence was achieved using different grid sizes. For water faucet problem, comparison of numerical results with analytical solution show good agreement and for shock tube problem, the results indicate that this method is highly capable in capturing discontinuities in two-phase flow.