A FENE-P k– ε turbulence model for low and intermediate regimes of polymer-induced drag reduction

► Two-equation turbulence model for FENE-P fluids for up to intermediate drag reduction. ► Modeled transport equation of rate of dissipation includes viscoelastic destruction. ► Closure for the nonlinear term in conformation equation based on exact equation. ► Improved closures for viscoelastic stress work. ► New contribution to viscoelastic turbulent diffusion. A new low-Reynolds-number k– ε turbulence model is developed for flows of viscoelastic fluids described by the finitely extensible nonlinear elastic rheological constitutive equation with Peterlin approximation (FENE-P model). The model is validated against direct numerical simulations in the low and intermediate drag reduction (DR) regimes (DR up to 50%). The results obtained represent an improvement over the low DR model of Pinho et al. (2008) [A low Reynolds number k– ε turbulence model for FENE-P viscoelastic fluids, Journal of Non-Newtonian Fluid Mechanics, 154, 89–108]. In extending the range of application to higher values of drag reduction, three main improvements were incorporated: a modified eddy viscosity closure, the inclusion of direct viscoelastic contributions into the transport equations for turbulent kinetic energy ( k) and its dissipation rate, and a new closure for the cross-correlations between the fluctuating components of the polymer conformation and rate of strain tensors ( NLT ij ). The NLT ij appears in the Reynolds-averaged evolution equation for the conformation tensor (RACE), which is required to calculate the average polymer stress, and in the viscoelastic stress work in the transport equation of k. It is shown that the predictions of mean velocity, turbulent kinetic energy, its rate of dissipation by the Newtonian solvent, conformation tensor and polymer and Reynolds shear stresses are improved compared to those obtained from the earlier model.