Numerical study of Taylor bubble breakup in counter-current flow using large eddy simulation

This paper investigates dynamics of Taylor bubble in counter-current flows, leveraging large eddy simulations combined with the volume of fluid method. Utilizing the OpenFOAM framework, we have implemented a high-order Runge-Kutta time-integration scheme, along with a piecewise linear interface calculation method for precise geometric reconstruction of the bubble interface. We examine the performance of algebraic vs geometric capturing techniques in the context of Taylor bubble breakup, focusing specifically on the transitional flow regime with a liquid Reynolds number of 1400. Our results reveal that the geometric capturing technique offers superior accuracy, improving our understanding of the breakup process and providing valuable insight for multiphase flow simulations in various engineering fields. Our study also reveals the emergence of a secondary vortex in the turbulent wake region behind the Taylor bubble, a phenomenon most prominent at finer mesh resolutions. This vortex represents a novel discovery in counter-current Taylor bubble flows.