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.

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.