{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2025,7,30]],"date-time":"2025-07-30T11:22:18Z","timestamp":1753874538631,"version":"3.41.2"},"reference-count":61,"publisher":"AIP Publishing","issue":"9","license":[{"start":{"date-parts":[[2023,9,11]],"date-time":"2023-09-11T00:00:00Z","timestamp":1694390400000},"content-version":"vor","delay-in-days":41,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"},{"start":{"date-parts":[[2023,9,11]],"date-time":"2023-09-11T00:00:00Z","timestamp":1694390400000},"content-version":"tdm","delay-in-days":41,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100000266","name":"Engineering and Physical Sciences Research Council","doi-asserted-by":"publisher","award":["EP\/N5059681\/1"],"award-info":[{"award-number":["EP\/N5059681\/1"]}],"id":[{"id":"10.13039\/501100000266","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["pubs.aip.org"],"crossmark-restriction":true},"short-container-title":[],"published-print":{"date-parts":[[2023,8,1]]},"abstract":"<jats:p>A novel open-source anisotropic k\u2212\u03b5\u2212v2\u2212f model is presented for turbulent viscoelastic duct flow with dilute polymeric solutions described by the finitely extensible nonlinear elastic-Peterlin constitutive model. The turbulence model for channel and square duct flow of Newtonian fluids is adapted to incorporate the polymeric terms within the governing equations. All the required non-linear terms are validated with simple closure models and are assessed a priori against independent direct numerical simulation data in fully developed channel flow. The NLTij term, which accounts for the interaction between fluctuating components of the conformation tensor and the velocity gradient tensor, is modeled with the mean flow direction, ti, and wall-normal, ni, present in the Newtonian model, based on the streamwise alignment of mean polymer stretch. The implicit polymer effects on pressure\u2013strain are assessed with a simple ad hoc closure accounting for the reduced near-wall production of turbulent kinetic energy. The same closure is also adapted for the spanwise Reynolds stress predictions of polymer-enhanced secondary flow. The model performs well in channel flow and captures low, intermediate, and high drag reduction features for a wide range of rheological parameters. The capabilities are extended for square ducts (or any regular polygon) due to the symmetric modeling of the closure models, which can predict the mean streamwise and secondary flow features associated with second normal Reynolds stress differences. Accessible codes and models are crucial for the advancement and improvement of turbulent viscoelastic models, and an OpenFOAM C++ code package is developed and freely available on GitHub (https:\/\/github.com\/MikeMcDermott-Code\/v2f).<\/jats:p>","DOI":"10.1063\/5.0159668","type":"journal-article","created":{"date-parts":[[2023,9,11]],"date-time":"2023-09-11T10:42:19Z","timestamp":1694428939000},"update-policy":"https:\/\/doi.org\/10.1063\/aip-crossmark-policy-page","source":"Crossref","is-referenced-by-count":1,"title":["An open-source anisotropic k\u2212\u03b5\u2212v2\u2212f model for turbulent viscoelastic duct flows"],"prefix":"10.1063","volume":"35","author":[{"ORCID":"https:\/\/orcid.org\/0000-0002-5091-7170","authenticated-orcid":false,"given":"M.","family":"McDermott","sequence":"first","affiliation":[{"name":"School of Mechanical Engineering, University of Leeds 1 , Woodhouse Lane, Leeds LS2 9JT, United Kingdom"}]},{"ORCID":"https:\/\/orcid.org\/0009-0007-1474-8449","authenticated-orcid":false,"given":"T. A. E.","family":"Riou","sequence":"additional","affiliation":[{"name":"School of Mechanical Engineering, University of Leeds 1 , Woodhouse Lane, Leeds LS2 9JT, United Kingdom"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-3410-6482","authenticated-orcid":false,"given":"P. R.","family":"Resende","sequence":"additional","affiliation":[{"name":"Escola Superior de Tecnologia e Gest\u00e3o, Instituto Polit\u00e9cnico de Viana do Castelo 2 , 4900\u2013347 Viana do Castelo, Portugal"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-1058-2003","authenticated-orcid":false,"given":"M. C. T.","family":"Wilson","sequence":"additional","affiliation":[{"name":"School of Mechanical Engineering, University of Leeds 1 , Woodhouse Lane, Leeds LS2 9JT, United Kingdom"}]},{"ORCID":"https:\/\/orcid.org\/0000-0003-2825-0709","authenticated-orcid":false,"given":"A. M.","family":"Afonso","sequence":"additional","affiliation":[{"name":"Transport Phenomena Research Center, Faculty of Engineering, University of Porto 3 , Rua Dr. Roberto Frais s\/n, 4200\u2013465 Porto, Portugal"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-5647-1771","authenticated-orcid":false,"given":"G.","family":"de Boer","sequence":"additional","affiliation":[{"name":"School of Mechanical Engineering, University of Leeds 1 , Woodhouse Lane, Leeds LS2 9JT, United Kingdom"}]}],"member":"317","published-online":{"date-parts":[[2023,9,11]]},"reference":[{"first-page":"135","article-title":"Some observations on the flow of linear polymer solutions through straight tubes at large Reynolds numbers","year":"1948","key":"2023102617535174800_c1"},{"key":"2023102617535174800_c2","doi-asserted-by":"publisher","first-page":"367","DOI":"10.1146\/annurev.fl.01.010169.002055","article-title":"Drag reduction by additives","volume":"1","year":"1969","journal-title":"Annu. Rev. Fluid Mech."},{"key":"2023102617535174800_c3","unstructured":"J. W.\n            Hoyt\n          , \u201c\nTurbulent flow of drag-reducing suspensions,\u201d Technical Report No. NUC-TP-299 (\nNaval Undersea Center, \nSan Diego, CA, 1972)."},{"key":"2023102617535174800_c4","doi-asserted-by":"publisher","first-page":"305","DOI":"10.1017\/S0022112067001442","article-title":"The Toms phenomenon: Turbulent pipe flow of dilute polymer solutions","volume":"30","year":"1967","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c5","doi-asserted-by":"publisher","first-page":"625","DOI":"10.1002\/aic.690210402","article-title":"Drag reduction fundamentals","volume":"21","year":"1975","journal-title":"AIChE J."},{"key":"2023102617535174800_c6","doi-asserted-by":"publisher","first-page":"1016","DOI":"10.1063\/1.1345882","article-title":"Budgets of Reynolds stress, kinetic energy and streamwise enstrophy in viscoelastic turbulent channel flow","volume":"13","year":"2001","journal-title":"Phys. Fluids"},{"key":"2023102617535174800_c7","doi-asserted-by":"publisher","first-page":"251","DOI":"10.1017\/S0022112003005305","article-title":"Turbulent channel flow near maximum drag reduction: Simulations, experiments and mechanisms","volume":"490","year":"2003","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c8","doi-asserted-by":"publisher","first-page":"311","DOI":"10.1007\/s10494-005-9002-6","article-title":"New answers on the interaction between polymers and vortices in turbulent flows","volume":"74","year":"2005","journal-title":"Flow, Turbul. Combust."},{"key":"2023102617535174800_c9","doi-asserted-by":"publisher","first-page":"23","DOI":"10.1016\/j.jnnfm.2005.12.012","article-title":"Influence of rheological parameters on polymer induced turbulent drag reduction","volume":"140","year":"2006","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c10","doi-asserted-by":"publisher","first-page":"N19","DOI":"10.1080\/14685248.2012.685522","article-title":"Some dynamical features of the turbulent flow of a viscoelastic fluid for reduced drag","volume":"13","year":"2012","journal-title":"J. Turbul."},{"key":"2023102617535174800_c11","doi-asserted-by":"publisher","first-page":"235","DOI":"10.1146\/annurev.fluid.40.111406.102156","article-title":"Mechanics and prediction of turbulent drag reduction with polymer additives","volume":"40","year":"2008","journal-title":"Annu. Rev. Fluid Mech."},{"key":"2023102617535174800_c12","doi-asserted-by":"publisher","first-page":"121302","DOI":"10.1063\/1.5129619","article-title":"Turbulent drag reduction by polymer additives: Fundamentals and recent advances","volume":"31","year":"2019","journal-title":"Phys. Fluids"},{"key":"2023102617535174800_c13","doi-asserted-by":"publisher","first-page":"213","DOI":"10.1016\/0377-0257(80)85007-5","article-title":"Polymer solution rheology based on a finitely extensible bead\u2013spring chain model","volume":"7","year":"1980","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c14","doi-asserted-by":"publisher","first-page":"281","DOI":"10.1017\/S0022112007006611","article-title":"Effects of polymer stresses on eddy structures in drag-reduced turbulent channel flow","volume":"584","year":"2007","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c15","doi-asserted-by":"publisher","first-page":"99","DOI":"10.1016\/j.jnnfm.2013.09.007","article-title":"A viscoelastic k\u2212\u03b5\u2212v2\u2212f turbulent flow model valid up to the maximum drag reduction limit","volume":"202","year":"2013","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c16","doi-asserted-by":"publisher","first-page":"52","DOI":"10.1016\/j.ijheatfluidflow.2013.05.016","article-title":"Analysis of polymer drag reduction mechanisms from energy budgets","volume":"43","year":"2013","journal-title":"Int. J. Heat Fluid Flow"},{"key":"2023102617535174800_c17","doi-asserted-by":"publisher","first-page":"A11","DOI":"10.1017\/jfm.2020.402","article-title":"Direct numerical simulations of turbulent viscoelastic jets","volume":"899","year":"2020","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c18","doi-asserted-by":"publisher","first-page":"A26","DOI":"10.1017\/jfm.2022.559","article-title":"Turbulent planar wakes of viscoelastic fluids analysed by direct numerical simulations","volume":"946","year":"2022","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c19","doi-asserted-by":"publisher","first-page":"041301","DOI":"10.1103\/PhysRevFluids.5.041301","article-title":"Increase of turbulent drag by polymers in particle suspensions","volume":"5","year":"2020","journal-title":"Phys. Rev. Fluids"},{"key":"2023102617535174800_c20","doi-asserted-by":"publisher","first-page":"1057","DOI":"10.1017\/jfm.2018.858","article-title":"Turbulent duct flow with polymers","volume":"859","year":"2019","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c21","doi-asserted-by":"publisher","first-page":"513","DOI":"10.1134\/S0015462821040091","article-title":"Prandtl's secondary flows of the second kind. problems of description, prediction, and simulation","volume":"56","year":"2021","journal-title":"Fluid Dyn."},{"key":"2023102617535174800_c22","unstructured":"L. C.\n            Hoagland\n          , \u201c\nFully developed turbulent flow in straight rectangular ducts: Secondary flow, its cause and effect on the primary flow,\u201d Ph.D. thesis (\nMassachusetts Institute of Technology, 1962)."},{"key":"2023102617535174800_c23","doi-asserted-by":"publisher","first-page":"306","DOI":"10.1007\/BF02079937","article-title":"Untersuchungen \u00fcber turbulente str\u00f6mungen in nicht kreisf\u00f6rmigen rohren","volume":"1","year":"1930","journal-title":"Ing.-Arch."},{"key":"2023102617535174800_c24","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1061\/JYCEAJ.0000866","article-title":"Turbulent flow in rectangular ducts","volume":"89","year":"1963","journal-title":"J. Hydraul. Div."},{"key":"2023102617535174800_c25","doi-asserted-by":"publisher","first-page":"A36","DOI":"10.1017\/jfm.2021.1139","article-title":"Turbulent entrainment in viscoelastic fluids","volume":"934","year":"2022","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c26","doi-asserted-by":"publisher","first-page":"863","DOI":"10.1016\/0020-7225(82)90008-8","article-title":"On turbulent secondary flows in pipes of non-circular cross-section","volume":"20","year":"1982","journal-title":"Int. J. Eng. Sci."},{"key":"2023102617535174800_c27","doi-asserted-by":"publisher","first-page":"459","DOI":"10.1017\/S0022112087001319","article-title":"On nonlinear k \u2013 l and k\u2212\u03b5 models of turbulence","volume":"178","year":"1987","journal-title":"J. Fluid Mech."},{"key":"2023102617535174800_c28","doi-asserted-by":"publisher","first-page":"61","DOI":"10.1007\/BF00271429","article-title":"A numerical study of turbulent square-duct flow using an anisotropic k\u2212\u03b5 model","volume":"2","year":"1990","journal-title":"Theor. Comput. Fluid Dyn."},{"key":"2023102617535174800_c29","doi-asserted-by":"publisher","first-page":"331","DOI":"10.1017\/S0022112075003382","article-title":"A more general effective-viscosity hypothesis","volume":"72","year":"1975","journal-title":"J. Fluid Mech."},{"year":"2008","key":"2023102617535174800_c30","article-title":"Predictions of turbulent secondary flows using the v2\u2212f model"},{"key":"2023102617535174800_c31","first-page":"267","article-title":"RANS modelling for compressible and transitional flows","volume-title":"Summer Program Centre for Turbulence Research","year":"1998"},{"key":"2023102617535174800_c32","doi-asserted-by":"publisher","first-page":"713","DOI":"10.1007\/s00162-020-00545-9","article-title":"A priori tests of eddy viscosity models in square duct flow","volume":"34","year":"2020","journal-title":"Theor. Comput. Fluid Dyn."},{"key":"2023102617535174800_c33","doi-asserted-by":"publisher","first-page":"169","DOI":"10.1007\/s10494-005-1974-8","article-title":"A robust formulation of the v2\u2212f model","volume":"73","year":"2005","journal-title":"Flow, Turbul. Combust."},{"key":"2023102617535174800_c34","doi-asserted-by":"publisher","first-page":"332","DOI":"10.1016\/j.ijheatmasstransfer.2016.04.053","article-title":"A RANS model for heat transfer reduction in viscoelastic turbulent flow","volume":"100","year":"2016","journal-title":"Int. J. Heat Mass Transfer"},{"first-page":"735","article-title":"Reynolds stress modelling for drag-reducing viscoelastic flows","year":"2003","key":"2023102617535174800_c35"},{"key":"2023102617535174800_c36","doi-asserted-by":"publisher","first-page":"1581","DOI":"10.1063\/1.1689971","article-title":"Characteristic scales and drag reduction evaluation in turbulent channel flow of non-constant viscosity viscoelastic fluids","volume":"16","year":"2004","journal-title":"Phys. Fluids"},{"key":"2023102617535174800_c37","doi-asserted-by":"publisher","first-page":"89","DOI":"10.1016\/j.jnnfm.2008.02.008","article-title":"A low Reynolds number turbulence closure for viscoelastic fluids","volume":"154","year":"2008","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c38","doi-asserted-by":"publisher","first-page":"177","DOI":"10.1016\/j.jnnfm.2006.04.012","article-title":"Turbulent channel flow of dilute polymeric solutions: Drag reduction scaling and an eddy viscosity model","volume":"139","year":"2006","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c39","doi-asserted-by":"publisher","first-page":"639","DOI":"10.1016\/j.jnnfm.2011.02.012","article-title":"A FENE-P k\u2212\u03b5 turbulence model for low and intermediate regimes of polymer-induced drag reduction","volume":"166","year":"2011","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c40","doi-asserted-by":"publisher","first-page":"69","DOI":"10.1007\/s10494-012-9424-x","article-title":"Development of a low-Reynolds-number k\u2212\u03c9 model for FENE-P fluids","volume":"90","year":"2013","journal-title":"Flow, Turbul. Combust."},{"key":"2023102617535174800_c41","doi-asserted-by":"publisher","first-page":"30","DOI":"10.1016\/j.ijheatfluidflow.2018.07.004","article-title":"An improved k\u2212\u03b5 turbulence model for FENE-P fluids capable to reach high drag reduction regime","volume":"73","year":"2018","journal-title":"Int. J. Heat Fluid Flow"},{"key":"2023102617535174800_c42","first-page":"221","article-title":"Improvement of the energy distribution in isotropic turbulent viscoelastic fluid models","year":"2014","journal-title":"Proceedings of MEFTE"},{"key":"2023102617535174800_c43","doi-asserted-by":"publisher","first-page":"025508","DOI":"10.1088\/1873-7005\/aa9e37","article-title":"New developments in isotropic turbulent models for FENE-P fluids","volume":"50","year":"2018","journal-title":"Fluid Dyn. Res."},{"key":"2023102617535174800_c44","doi-asserted-by":"publisher","first-page":"8140","DOI":"10.3390\/app10228140","article-title":"A FENE-P k\u2212\u03b5 viscoelastic turbulence model valid up to high drag reduction without friction velocity dependence","volume":"10","year":"2020","journal-title":"Appl. Sci."},{"key":"2023102617535174800_c45","doi-asserted-by":"publisher","first-page":"108799","DOI":"10.1016\/j.ijheatfluidflow.2021.108799","article-title":"An improved k\u2212\u03c9 turbulence model for FENE-P fluids without friction velocity dependence","volume":"90","year":"2021","journal-title":"Int. J. Heat Fluid Flow"},{"key":"2023102617535174800_c46","doi-asserted-by":"publisher","first-page":"376","DOI":"10.1016\/j.jnnfm.2010.01.013","article-title":"Reynolds-averaged modelling of polymer drag reduction in turbulent flows","volume":"165","year":"2010","journal-title":"J. Non-Newtonian Fluid Mech."},{"key":"2023102617535174800_c47","doi-asserted-by":"publisher","first-page":"337","DOI":"10.1007\/s10494-008-9134-6","article-title":"One equation model for turbulent channel flow with second order viscoelastic corrections","volume":"81","year":"2008","journal-title":"Flow, Turbul. Combust."},{"key":"2023102617535174800_c48","doi-asserted-by":"publisher","first-page":"1338","DOI":"10.1016\/j.physd.2009.07.013","article-title":"A short review on drag reduction by polymers in wall bounded turbulence","volume":"239","year":"2010","journal-title":"Phys. D"},{"key":"2023102617535174800_c49","doi-asserted-by":"publisher","first-page":"220","DOI":"10.1016\/j.ijheatfluidflow.2015.05.017","article-title":"A Reynolds stress model for turbulent flow of homogeneous polymer solutions","volume":"54","year":"2015","journal-title":"Int. J. Heat Fluid Flow"},{"key":"2023102617535174800_c50","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1080\/14685248.2013.851385","article-title":"A Reynolds stress model for turbulent flows of viscoelastic fluids","volume":"14","year":"2013","journal-title":"J. Turbul."},{"first-page":"V01AT03A010","article-title":"Reynolds-averaged simulation on turbulent drag-reducing flows of viscoelastic fluid based on user-defined function in FLUENT package","year":"2014","key":"2023102617535174800_c51"},{"key":"2023102617535174800_c52","doi-asserted-by":"publisher","first-page":"65","DOI":"10.22059\/JCAMECH.2017.232031.138","article-title":"Parametric study of a viscoelastic RANS turbulence model in the fully developed channel flow","volume":"48","year":"2017","journal-title":"J. Comput. Appl. Mech."},{"key":"2023102617535174800_c53","doi-asserted-by":"publisher","first-page":"509","DOI":"10.1146\/annurev-fluid-010719-060107","article-title":"Numerical methods for viscoelastic fluid flows","volume":"53","year":"2021","journal-title":"Annu. Rev. Fluid Mech."},{"key":"2023102617535174800_c54","doi-asserted-by":"publisher","first-page":"035106","DOI":"10.1063\/1.1850920","article-title":"Viscoelastic effects on higher order statistics and on coherent structures in turbulent channel flow","volume":"17","year":"2005","journal-title":"Phys. Fluids"},{"key":"2023102617535174800_c55","doi-asserted-by":"publisher","first-page":"869","DOI":"10.1002\/1097-0363(20000730)33:6&lt;869::AID-FLD36&gt;3.0.CO;2-E","article-title":"Prediction of turbulent oscillatory flows in complex systems","volume":"33","year":"2000","journal-title":"Int. J. Numer. Methods Fluids"},{"key":"2023102617535174800_c56","doi-asserted-by":"publisher","first-page":"943","DOI":"10.1063\/1.869966","article-title":"Direct numerical simulation of turbulent channel flow up to","volume":"11","year":"1999","journal-title":"Phys. Fluids"},{"first-page":"303","article-title":"Assessment of non-local effect on pressure term in RANS modeling using a DNS database","year":"1998","key":"2023102617535174800_c57"},{"key":"2023102617535174800_c58","doi-asserted-by":"publisher","first-page":"2345","DOI":"10.1063\/1.1287517","article-title":"A new form of the elliptic relaxation equation to account for wall effects in RANS modeling","volume":"12","year":"2000","journal-title":"Phys. Fluids"},{"first-page":"63","article-title":"A turbulence model for polymer flows","year":"2004","key":"2023102617535174800_c59"},{"key":"2023102617535174800_c60","unstructured":"M.\n            McDermott\n          , \u201c\nRANS turbulence models for FENE-P viscoelastic fluids,\u201d Ph.D. thesis (\nUniversity of Leeds, 2022)."},{"key":"2023102617535174800_c61","doi-asserted-by":"publisher","first-page":"743","DOI":"10.1063\/1.869229","article-title":"Direct numerical simulation of the turbulent channel flow of a polymer solution","volume":"9","year":"1997","journal-title":"Phys. Fluids"}],"container-title":["Physics of Fluids"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/pubs.aip.org\/aip\/pof\/article-pdf\/doi\/10.1063\/5.0159668\/18188052\/095116_1_5.0159668.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"text-mining"},{"URL":"https:\/\/pubs.aip.org\/aip\/pof\/article-pdf\/doi\/10.1063\/5.0159668\/18188052\/095116_1_5.0159668.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/pubs.aip.org\/aip\/pof\/article-pdf\/doi\/10.1063\/5.0159668\/18188052\/095116_1_5.0159668.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,5,9]],"date-time":"2025-05-09T18:59:41Z","timestamp":1746817181000},"score":1,"resource":{"primary":{"URL":"https:\/\/pubs.aip.org\/pof\/article\/35\/9\/095116\/2910549\/An-open-source-anisotropic-k-v-2-f-model-for"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2023,8,1]]},"references-count":61,"journal-issue":{"issue":"9","published-print":{"date-parts":[[2023,8,1]]}},"URL":"https:\/\/doi.org\/10.1063\/5.0159668","relation":{},"ISSN":["1070-6631","1089-7666"],"issn-type":[{"type":"print","value":"1070-6631"},{"type":"electronic","value":"1089-7666"}],"subject":[],"published-other":{"date-parts":[[2023,8]]},"published":{"date-parts":[[2023,8,1]]},"article-number":"095116"}}