On the assessment of a VOF based compressive interface capturing scheme for the analysis of bubble impact on and bounce from a flat horizontal surface, International Journal of Multiphase Flow (2014), doi: http://dx.doi.org/10. 1016/ j.ijmultiphaseflow.2014.05.017 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.On the assessment of a VOF based compressive interface capturing scheme for the analysis of bubble impact on and bounce from a flat horizontal surface
AbstractThe process of free rise, collision on and bounce from a solid horizontal surface for a single isolated bubble is investigated by numerical simulations based on the Volume of Fluid method (VOF). The volume fraction advection equation is solved algebraically using the compressive scheme implemented in the CFD open source library (OpenFOAM R ) using both axi-symmetrical and three dimensional domains. The solution sensitivity to the mesh refinement towards the solid boundary and the contact angle formulation (static and dynamic) are assessed with two different fluid mixtures for a range of Bond numbers [0.298 − 1.48] and two different surface hydrophilicities. Numerical results are assessed against published as well as new experiments to include both axi-symmetrical and three dimensional rise trajectories. The investigation addresses the liquid microfilm formation and drainage considering both flow and pressure fields and bubble dynamic characteristics over successive rebounds. Results highlight the importance of resolving the liquid micro layer at the interface between the gas and solid surface in particular in the case of superhydrophobic surfaces. A coarse mesh is shown to precipitate the liquid film drainage. This results in early formation of a triple phase contact line (TPCL) which can occur as soon as the first rebound whereas physical observations indicate that this typically happens much later at a stage when * Corresponding authors a significant part of the bubble kinetic energy has been dissipated following several rebounds. As a result numerical predictions are shown to be much more sensitive to the contact angle formulation than when a refined mesh allows a more accurate representation of the film drainage. In this case, static and dynamic contact angle models give broadly similar rebound characteristics. Following validation, the numerical simulations are used to provide some useful insight in the mechanisms driving the film drainage and the gas liquid interface as it interacts with the solid surface.