A potential alternative CFD simulation for steady Carreau–Bird law-based shear thickening model: Part-I

• Published in: Journal of the Brazilian Society of Mechanical Sciences and Engineering Vol. 41; no. 4; pp. 1 - 13
• Main Authors: Rehman, Khalil Ur, Malik, M. Y., Mahmood, R., Kousar, N., Zehra, Iffat
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer Berlin Heidelberg 01.04.2019; Springer Nature B.V
• Subjects: Aerodynamic coefficients; Circular cylinders; Computational fluid dynamics; Computer simulation; Drag; Engineering; Finite element method; Fluid flow; Lift; Mathematical models; Mechanical Engineering; Newtonian fluids; Non Newtonian fluids; Nonlinear differential equations; Nonlinear equations; Power law; Reynolds number; Shear thickening (liquids); Technical Paper; Thickening; Velocity; Velocity distribution; Finite element method; Linear and parabolic profiles; Drag and lift coefficients; Power law fluid model
• ISSN: 1678-5878; 1806-3691
• DOI: 10.1007/s40430-019-1664-4

The current article provides the numerical investigation into an infinite-length circular cylinder placed as an obstacle in the flow of non-Newtonian fluid. To be more specific, a channel of length 2.2 m and height 0.41 m is considered. The Power law fluid model is carried out with Carreau–Bird law as a non-Newtonian fluid model, and both the Power law linear (constant) and parabolic velocity profiles are initiated simultaneously at an inlet of the channel. The right wall as an outlet is carried with Neumann condition. The relative velocity of Power law fluid particles with both the lower and upper walls is taken zero. A mathematical model is structured in terms of nonlinear differential equations. A well-trusted numerical technique named finite element method is adopted commercially. The LBB-stable element pair is utilized to approximate the velocity and pressure. The nonlinear iterations are stopped when the residual is dropped by 10 - 6 . The impact of Power law index and Reynolds number on the primitive variables is inspected. The obtained observations in this direction are provided by means of both the contour plots and line graphs. Due to a circular obstacle, both the drag and lift coefficients are evaluated around the outer surface of an obstacle towards the higher values of the Power law index. The numerical values of drag and lift coefficients up to various refinement mesh levels of domain are provided by way of tables. It is noticed that the parabolic velocity profile at an inlet of channel is compatible as compared to the linear velocity profile. Further, both the drag and lift coefficients are increasing function of Power law index.