Direct numerical simulation of turbulent non-Newtonian flow using OpenFOAM

• Published in: Applied Mathematical Modelling 2019-08, Vol.72, p.50-67
• Main Authors: Zheng, E.Z., Rudman, M., Singh, J., Kuang, S.B.
• Format: Article
• Language: English
• Subjects: Computational fluid dynamics; Computer simulation; Direct numerical simulation; Direct numerical simulation (DNS); Fluid flow; Generalized Newtonian fluid; Newtonian fluids; Non Newtonian flow; Non Newtonian fluids; OpenFOAM; Pipe flow; Shear flow; Shear thinning (liquids); Shear-thinning fluid; Turbulence; Turbulent flow; Turbulent pipe flow; Viscosity
• ISSN: 0307-904X; 1088-8691; 0307-904X
• DOI: 10.1016/j.apm.2019.03.003

•OpenFOAM predictions for turbulent Newtonian fluids match very well with DNS references for mean turbulence statistics.•OpenFOAM predicts turbulent shear-thinning flow to be more transitional compared to a high-order spectral element DNS code.•Turbulence statistics predicted by OpenFOAM differ by at most 16% for shear-thinning fluids, and usually be less than 10%.•OpenFOAM scales well from 8 to 512 CPUs but has poor intranode scalability for less than 8CPUs. Understanding transition and turbulence in the flow of shear-thinning non-Newtonian fluids remains substantially unresolved and additional research is required to develop better computational methods for wall-bounded turbulent flows of these fluids. Previous DNS studies of shear-thinning fluids mainly use purpose-built codes and simple geometries such as pipes and channels. However in practical application, the geometry of mixing vessels, pumps and other process equipment is far more complex, and more flexible computational methods are required. In this paper a general-purpose DNS approach for shear-thinning fluids is undertaken using the OpenFOAM CFD library. DNS of turbulent Newtonian and non-Newtonian flow in a pipe flow are conducted and the accuracy and efficiency of OpenFOAM are assessed against a validated high-order spectral element-Fourier DNS code – Semtex. The results show that OpenFOAM predicts the flow of shear-thinning fluids to be a little more transitional than the predictions from Semtex, with lower radial and azimuthal turbulence intensities and higher axial intensity. Despite this, the first and second order turbulence statistics differ by at most 16%, and usually much less. An assessment of the parallel scaling of OpenFOAM indicates that OpenFOAM scales very well for the CPUs from 8 to 512, but the intranode scalability is poor for less than 8CPUs. The present work shows that OpenFOAM can be used for DNS of shear-thinning fluids in the simple case of pipe flow, and suggests that more complex flows, where flow separation is often important, are likely to be simulated with accuracies that are acceptably good for engineering application.