مطالعه عددي تبخير در يك كانال متخلخل قائم با روش حجم سيال در اپن‌فوم

Numerical Study of Evaporation in a Vertical Porous Channel By the volume of fluid method in OpenFOAM

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

In this study, using volume of fluid method in open source software OpenFOAM, the phenomenon of evaporation in the porous medium was analyzed. At the beginning of the solution, the system consists of a water phase and a porous copper environment. In the next steps of numerical simulation and as a result of partial evaporation of water, the vapor phase appears as the second fluid phase. Water and vapor are assumed to be incompressible and incompatible, and the phenomenon of evaporation occurs unevenly. The interface between phases is modeled by the VOF method and the Lee model has been used to mass transfer between two phases of water and vapor. For surface tension between phases, the continuous surface force (CSF) method was considered. The comparison of simulation results with experimental results showed that the combined solver of porous medium evaporation would well estimate the rate of evaporation at different sections of the channel. In addition, the results of the wall temperature indicate that the channel is divided into two zones of heating and evaporation. In the region of heating, the temperature increases linearly with the channel length to reach saturation temperature. After the point of saturation, the wall temperature first remains constant and eventually forms an oscillatory shape, in which locally there are temperature jumps. The evaporated flow rate also increases at high intensity first, but in the end regions of the porous channel, its growth rate is slow.