TY - JOUR

T1 - One equation model for turbulent channel flow with second order viscoelastic corrections

AU - Pinho, F. T.

AU - Sadanandan, B.

AU - Sureshkumar, R.

N1 - Funding Information:
Acknowledgements FT Pinho acknowledges the funding from Fundação Calouste Gulbenkian (Scholarship 72259) and Fundação para a Ciência e a Tecnologia (Scholarship SFRH/BSAB/507/2005). Funding from FEDER and FCT via project POCI/56342/EQU/2004 is also gratefully acknowledged. Partial support from the NSF funding through grant CTS 9874813 (R.S.) is acknowledged.

PY - 2008/10

Y1 - 2008/10

N2 - A modified second order viscoelastic constitutive equation is used to derive a k-l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P w could be either positive or negative. When P w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P w ∈<∈0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k-l formulation rather than from the closure for P w . An alternative closure for P w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P w ∈>∈0 but has also a small region of negative P w in the buffer layer, was also successfully tested. This second model for P w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.

AB - A modified second order viscoelastic constitutive equation is used to derive a k-l type turbulence closure to qualitatively assess the effects of elastic stresses on fully-developed channel flow. Specifically, the second order correction to the Newtonian constitutive equation gives rise to a new term in the momentum equation involving the time-averaged elastic shear stress and in the turbulent kinetic energy transport equation quantifying the interaction between the fluctuating elastic stress and rate of strain tensors, denoted by P w , for which a closure is developed and tested. This closure is based on arguments of isotropic turbulence and equilibrium in boundary layer flows and a priori P w could be either positive or negative. When P w is positive, it acts to reduce the production of turbulent kinetic energy and the turbulence model predictions qualitatively agree with direct numerical simulation (DNS) results obtained for more realistic viscoelastic fluid models with memory which exhibit drag reduction. In contrast, P w ∈<∈0 leads to a drag increase and numerical breakdown of the model occurs at very low values of the Deborah number, which signifies the ratio of elastic to viscous stresses. Limitations of the turbulence model primarily stem from the inadequacy of the k-l formulation rather than from the closure for P w . An alternative closure for P w , mimicking the viscoelastic stress work predicted by DNS using the Finitely Extensible Nonlinear Elastic-Peterlin fluid model, which is mostly characterized by P w ∈>∈0 but has also a small region of negative P w in the buffer layer, was also successfully tested. This second model for P w leads to predictions of drag reduction, in spite of the enhancement of turbulence production very close to the wall, but the equilibrium conditions in the inertial sub-layer were not strictly maintained.

KW - Drag reduction

KW - K-l model

KW - Second order constitutive equation

KW - Viscoelastic

UR - http://www.scopus.com/inward/record.url?scp=51849116940&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=51849116940&partnerID=8YFLogxK

U2 - 10.1007/s10494-008-9134-6

DO - 10.1007/s10494-008-9134-6

M3 - Article

AN - SCOPUS:51849116940

VL - 81

SP - 337

EP - 367

JO - Flow, Turbulence and Combustion

JF - Flow, Turbulence and Combustion

SN - 1386-6184

IS - 3

ER -