مدل‌سازی جریان و انتقال حرارت دوفازی در داخل کانال و اطراف لوله با روش لتیس بولتزمن

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

Simulation of Two-Phase Flow and Heat Transfer in a Channel and around a Tube by Lattice-Boltzmann Method

Determination of multiphase flow dynamics and thermal behavior of twophase flow in a channel are of importance. The smallscale surface tension effect and related simulation efficiency, precision, and stability, have caused mesoscopic Lattice Boltzmann method broadening application. In the current study, the thermalhydraulic behavior of subcooled falling flow in a vertical channel and around a single horizontal tube is simulated by using the Lee method and phasefiled model, and thermal passive scalar model. The modified curved boundary conditions and two different boundary conditions for side boundaries are investigated. The density ratio is 20 and other property’s ratios of water are applied, and the outside diameter of the tube is 28.9mm. The flow, temperature, and pressure fields are presented and a detailed understanding of the movement of the threephase contact line, circulating flow and local and average Nusselt numbers are determined. The film thickness, thermal boundary layer variation by the film thickness, Reynold number effect on Nusselt number and mass conservation are investigated as verification. The results have shown good consistency and high effectiveness in the simulation of multiphase gasliquid flows in the presence of a circular obstacle, and for viscosity and thermal diffusivity ratios of water.