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

یک مدل شبکه بولتزمن دوفازی قدرتمند با قابلیت مدلسازی نسبت چگالی بالا برای شبیه‌سازی برخورد حباب به مانع متخلخل به‌کار گرفته شده است که می‌تواند جریان دوفازی مخلوط نشدنی با نسبت چگالی 1000 را با بقای جرم بسیار مطلوب شبیه‌سازی کند. در این مقاله با مدلسازی مانع در دامنه حل جریان دوفازی، نحوه تغییرات توابع توزیع در دیواره‌ها و گوشه‌های مانع بیان شده، نحوه اعمال شرط مرزی ترشوندگی سطح بر روی مانع تشریح گردیده و نشان داده شده است که با وجود اعمال شرط مرزی مربوط به مانع نیز بقای جرم مدل به‌شکل مطلوبی حفظ می‌شود. پس از ارزیابی دقت و قابلیت مدل و صحت پیاده‌سازی آن ابتدا با چند مساله پایه، برخورد یک حباب بالارونده با نسبت چگالی 1000 به یک مانع متخلخل مسطح شبیه‌سازی شده و تاثیر زاویه تماس، عدد اتوس و نسبت تخلخل و در تغییر شکل و نحوه عبور یا عدم عبور حباب از مانع به‌صورت سیستماتیک مورد بررسی قرار گرفته است. در برخورد حباب به مانع متخلخل با نسبت تخلخل پایین، در صورتی که زاویه تماس کم باشد حباب در زیر مانع باقی مانده و به سکون می‌رسد. در زوایای تماس بالا، خاصیت آب‌گریزی سطح مانع، حباب را به داخل تخلخل‌ها می‌کِشد و حباب به سطح بالای مانع منتقل می‌شود و بر روی سطح مانع به سکون می‌رسد. در سایر موارد حباب به طور کامل از مانع عبور می‌کند و از آن جدا می‌شود. خطای بقای جرم در عبور حباب از مانع متخلخل از مرتبه 1110 است که بسیار مطلوب ارزیابی می‌شود.

Numerical simulation of rising bubble striking a porous obstacle using mass-conserving lattice Boltzmann method

A powerful twophase lattice Boltzmann model with the ability of modeling high density ratio is applied to simulate a rising bubble striking a porous obstacle. This model is able to simulate immiscible twophase flow with density ratio of 1000 and result in desirable mass conservation. In present research, a porous obstacle is posed in twophase flow domain, bounce back and wetting boundary conditions at walls and corners is discussed and showed that after implementation of obstacle boundary conditions, mass conservation of the model is preserved. Accuracy and ability of the model firstly examined by some basic problems. Next, striking of a rising bubble with 1000 density ratio to a porous obstacle is simulated and the effect of contact angle, Eotvos number and porosity ratio in deformation and passing of the bubble from the obstacle is investigated systematically. Different porosity ratios and contact angles, result in different bubble behavior striking the porous obstacle; In low porosity ratios and low contact angles, the bubble remains below the obstacle. At high contact angles, the hydrophobicity of the obstacle surface draws the bubble into the porosities, and the bubble moves to the top of the obstacle and stays on the top surface of the obstacle. In other cases, the bubble completely passes through the obstacle and separates it. Mass conservation error of bubble passing the porous obstacle is of order of 1011 which is completely desirable.