Polymer additives are commonly utilized to manipulate bubbly flows in various applications. Here we investigate the effects of clean and contaminated bubbles driven upward (upflow) in Newtonian and viscoelastic turbulent channel flows. Interface-resolved direct numerical simulations are performed to examine sole and combined effects of soluble surfactant and viscoelasticity using an efficient three-dimensional finite-difference-front-tracking method. The incompressible flow equations are solved fully coupled with the FENE-P viscoelastic model and the equations governing interfacial and bulk surfactant concentrations. The latter coupling is accomplished by a nonlinear equation of state that relates the surface tension to the surfactant concentration. For Newtonian turbulent bubbly flows, the effects of Triton X-100 and 1-pentanol surfactant are examined. It is observed that the sorption kinetics highly affect the dynamics of bubbly flow. A minute amount of Triton X-100 is found to be sufficient to prevent the formation of bubble clusters restoring the single-phase behavior while even two orders of magnitude more 1-pentanol surfactant is not adequate to prevent the formation of layers. For viscoelastic turbulent flows, it is found that the viscoelasticity promotes formation of the bubble wall-layers and thus the polymer drag reduction is completely lost for the surfactant-free bubbly flows, while the addition of small amount of surfactant (Triton X-100) in this system restores the polymer drag reduction resulting in 25% drag reduction for the Wi=4 case.

Interface-resolved direct numerical simulations are performed to examine the combined effects of soluble surfactant and viscoelasticity on the structure of a bubbly turbulent channel flow. The incompressible flow equations are solved fully coupled with the FENE-P viscoelastic model and the equations governing interfacial and bulk surfactant concentrations. The latter coupling is achieved through a non-linear equation of state which relates the surface tension to the surfactant concentration at the interface. The two-fluid Navier-Stokes equations are solved using a front-tracking method, augmented with a very efficient FFT-based pressure projection method that allows for massively parallel simulations of turbulent flows. It is found that, for the surfactant-free case, bubbles move toward the wall due to inertial lift force, resulting in formation of wall layers and a significant decrease in the flow rate. Conversely, a high-enough concentration of surfactant changes the direction of lateral migration of bubbles, i.e., the contaminated bubbles move toward the core region and spread out across the channel. When viscoelasticity is considered, viscoelastic stresses counteract the Marangoni stresses, promoting formation of bubbly wall-layers and consequently strong decrease in the flow rate. The formation of bubble wall-layers for combined case depends on the interplay of the inertial and elastic, and Marangoni forces.