The application of interface tracking methods to bubbly flow modeling has grown in recent years due to improvements in computing performance and development of more efficient solvers. However, the standard formulation of most interface tracking methods is not designed to physically handle the interface interactions at reasonable grid sizes. Regardless of the method used, a high grid resolution is required in the liquid film region in order to properly model drainage process during bubble interaction, which in certain conditions prevents the coalescence. This makes large scale (many bubbles) simulations unaffordable. One of the popular interface tracking approached is the level-set (LS) method. To simulate realistic bubble coalescence behavior in the LS method an algorithm with the capability of delaying or preventing the process of multiple simultaneous coalescence events has been developed.
Bubble interaction plays a significant role in high void fraction flow behavior and affects the transition to other flow regimes (e.g. churn-turbulent or slug flows). The described algorithm allows to improve the accuracy of predicting coalescence events in these relevant cases and has been tested in a variety of conditions and computational meshes.
This novel algorithm uses the LS method field to detect when bubbles are in close proximity, indicating a potential coalescence event, and applies a subgrid scale force to simulate the unresolved liquid drainage force. The subgrid-model is introduced by locally modifying the surface tension force near the liquid film drainage area. The algorithm can also simulate the liquid drainage time of the thin film by controlling the length of time the increased surface tension has been applied. Thus a new method of modeling bubble coalescence has been developed.
Several test cases were designed to demonstrate the capabilities of the algorithm. The simulations, including a mesh study, confirmed the abilities to identify and prevent coalescence as well as implement the time tracking portion, with an additional 10–25% computational cost. Ongoing tests aim to verify the algorithm’s functionality for simulations with different flow conditions, a ranging number of bubbles, and both structured and unstructured computational mesh types. Specifically, a bubble rising towards a free surface provides a test of performance and demonstrates the ability to consistently prevent coalescence. In addition, a two bubble case and a seven bubble case provide a more complex demonstration of how the algorithm performs for larger simulations. These cases are compared to much more expensive simulations capable of resolving the liquid film drainage (through very high local mesh resolution), to investigate how the algorithm replicates the liquid film drainage process.