Modelling and Simulation in Engineering

Modelling and Simulation in Engineering / 2021 / Article

Research Article | Open Access

Volume 2021 |Article ID 5583110 |

Asif Mehmood, Rashid Mahmood, Afraz Hussain Majeed, Farah Jabeen Awan, "Flow of the Bingham-Papanastasiou Regularized Material in a Channel in the Presence of Obstacles: Correlation between Hydrodynamic Forces and Spacing of Obstacles", Modelling and Simulation in Engineering, vol. 2021, Article ID 5583110, 14 pages, 2021.

Flow of the Bingham-Papanastasiou Regularized Material in a Channel in the Presence of Obstacles: Correlation between Hydrodynamic Forces and Spacing of Obstacles

Academic Editor: Houari Ameur
Received22 Feb 2021
Revised01 Apr 2021
Accepted13 Apr 2021
Published28 Apr 2021


The numerical modeling and simulation for the stationary Bingham fluid flow around two confined circular cylinders with various gap ratios are studied. The singularity in the model’s apparent viscosity is dealt by Papanastasiou’s regularization. The model equations are discretized by adopting the methodology based on finite element method (FEM) by choosing a mixed higher order LBB-stable finite element pair. The direct solver PARADISO has been utilized to solve the linearized system of equations. Hydrodynamic forces represented by drag and lift coefficients are computed, and a correlation coefficient is calculated for the gap ratios and for several values of the Bingham number . Line graphs for horizontal and vertical velocities are drawn. Moreover, velocity and pressure profiles are plotted for pertinent values of the parameters. Plug and shear zones are revealed via velocity snapshots in the domain. Pressure is nonlinear in the vicinity of the obstacles and becomes linear downstream in the cylinders as expected in channel flows.

1. Introduction

The theoretical hydrodynamics prospered in nineteenth century but was a bit distant from hydraulic science. The introduction of boundary layer concept in the 20th century reduced this gap to zero. The fluid mechanics suffering due to the complexity of fluid models in analytical solutions have crossed this barrier due to numerical procedures and invention of computer and technological advancement in the recent past. An important question is what should be the next step? The intellect is the one who looks beyond the senses. Non-Newtonian fluids, high Reynolds numbers, and turbulent flows [13] captured the attention of mathematicians and engineers in the past few decades. The mysteries of turbulence are still being revealed in order to understand the flow phenomenon. The notable contributions of Euler, Bernoulli, and d’ Alembert have been explained in detail by Darrigol and Frisch [3].

Over the past two decades, the flows of both Newtonian and non-Newtonian fluids have been investigated thoroughly [410]. For various values of Reynolds numbers, a variety of flow patterns can be observed even for Newtonian fluids [4, 5]. The impact of fluid interaction with bluff bodies is very complex. So far, less information is available in literature about such complexities; a few studies can be found in [47].

The hydrodynamic forces such as drag and lift are the key factors to be investigated in studying solid-liquid interaction. While investigating the flow characteristics when more than one obstacle is there, the arrangement and placement of obstacles are of particular interest, as it has meaningful impact on hydrodynamic forces. Two cylinders with different arrangements is the simplest of such cases [8]. The shear and dead zones are dependent not only on the placement of obstacles but on their sizes as well. Such zones are important to be investigated as they have practical significance in many fields such as placement of means for bridges, fluid flow in channels with submerged obstacles, and of course in food industry while dealing with mixing of different products. Tandem and side-by-side are the two basic configurations of obstacle arrangements, of which tandem arrangement captured more attention of researchers. The drag and lift have different impacts on up- and downstream obstacles in tandem arrangements. The gap spacing of obstacles is also of keen interest to control the impact of fluid forces especially drag and lift.

While studying the rheological aspects of nonlinear viscosity fluids, there are two types of fluid models to deal with: shear rate-dependent models [1117] and yield stress models [18, 19]. The materials like toothpaste and tomato ketchup are yield stress oriented. The rheology of such is investigated by many researchers [1824]. There are some recent discussion and hypothesis about the existence of yield stress. One statement is that the fluid behaves like solid below a certain threshold value of yield stress (of course fluid flows beyond that threshold value). The second statement summarizes that, all solids and fluids can flow over a long period of time [25, 26]. The first statement looks more practical as described in [27]. Beris et al. [28] were the first to investigate yield stress of the Bingham fluids using FEM. Adachi and Yoshioka [29] also established some results of viscous fluids past a cylinder. Mitsoulis [30] also used FEM to investigate the hydrodynamic forces of the Bingham fluid flowing around the cylinder. Some studies are available in literature which confined to Newtonian fluid past a couple of cylinders; one of these is proposed by Vakil and Green [31] and produced some results for Reynolds number between 1 and 20. Jossic and Magnin produced some results [32] of two cylinders by fixing the range of the Bingham number up to 40. Recently, Koblitz et al. [33] investigated the hydrodynamic forces for very high values of Bn. Calder and Yezzi [34] considered a rate of convergence and application for obstacle problems. In some studies mentioned therein [3538], authors have numerically investigated the flows of incompressible non-Newtonian fluids in various computational domains. Ameur and Menni [36] performed simulations at high Reynolds numbers giving rise to a turbulent flow in a backward facing step configuration, and a comparison with the experimental studies is performed. Houssem and Mohamed [3941] analyzed a mixed convective thermal flow over a single cylinder and two cylinders. They have measured the effects of wall proximity of obstacles and gap ratios for varying Reynolds and Richardson number. In addition, they have also considered the effects of the Reynolds and Richardson numbers on fluid forces and the Nusselt number.

We organize the rest of the manuscript as follows. “Mathematical Modeling” is concerned with the mathematical modeling of flow equations for a regularized Bingham fluid. The nondimensional reformulation of model equations is also presented in the same section. Problem setup and numerical approach including element type and solvers are described in “Problem Setup.” “Results and Discussions” is devoted to display the results of numerical simulations for the involved parameters and the analysis of gap ratio on hydrodynamic forces. The correlation coefficient between the forces versus gap between the cylinders is computed. The conclusion of the present study is revealed in “Conclusions.”

2. Mathematical Modeling

The set of coupled partial differential equations for stationary incompressible flows in dimensional form are defined generally as below where the all terms have their specific meanings. A simplified rheological relationship is proposed by Bingham [4] fluid for viscoplastic materials as where , , , and , represent the yield stress, stress tensor, plastic viscosity, and shear rate, respectively. The strain tensor is defined as

Here, velocity vector is denoted as . The shear rate and stress magnitude are defined as

Unyielded regions produced by (3) do not create any problem as long as , but apparent viscosity gives rise to a singularity in the limit . Computational schemes such as finite element method implemented in our study do not allow such singularities. To circumvent this issue, some regularization schemes are present in the literature. The purpose of regularization schemes is to replace discontinuous apparent viscosity for with such an expression that remains bounded for arbitrarily small and approximating the rheological behavior at the same time. In this direction, we removed the singularity with the help of the regularization introduced by Papanastasiou [20] as

Here, the parameter is showing the stress growth. Owing to Equation (4), the viscosity can be written as which is applicable in the whole flow domain. In order to introduce the nondimensional numbers for the simulation purposes, we choose as reference length and velocity, respectively, and nondimensional variables and obtain the following formulation: in which where and are the Reynolds and Bingham number. The stress growth parameter is now given by . The viscosity in dimensionless form is

is the nondimensional analogue of .

3. Problem Setup

3.1. Flow Configuration

Two equally sized cylinders with diameter are placed in a rectangular channel at several relative positions. The dimension of the computational domain is . The circular obstacle is fixed at , and is placed at different locations to alter spacing between the cylinders . The up and down walls of channel are set with no slip conditions. An inflow parabolic profile with is exposed to the inlet of the channel, and a do-nothing boundary condition at the outlet is chosen (see Figure 1).

3.2. Numerical Approach and Solvers

A wide range of flow problems can be described with the incompressible Navier-Stokes equations provided in (2). These equations describe real-life processes and help to understand nature. However, there are many difficulties associated with the numerical treatment of the system given in (1) and (2). Firstly, the nonlinearity of the convective term gives rise to numerical instabilities at larger values of Reynolds numbers. The second problem is caused by the incompressibility constraint that gives rise to checkerboard pressure modes resulting in pressure oscillations. The third problem is associated with the nonlinear algebraic system that results after the discretization of these equations. The size of the global matrix is very large and in general has bad condition number, so to tackle these issues, efficient iterative solvers with preconditioning techniques [42] are required.

Keeping in view all these difficulties, the methodology was adopted in the finite element method. The LBB-stable FEM pair is chosen to approximate the discrete velocity and pressure. A hybrid coarse mesh is first created and then is successively refined to meet the grid independent results. For the solution of nonlinear algebraic systems, Newton’s method is selected to linearize the system, and a direct linear solver PARADISO with special reordering of unknowns is applied for the linearized system. To stop the nonlinear iterative process, we adopt the following criteria:

Here, represents the iteration number, and denotes a component of solution. Figure 2 shows the computational grids at refinement level 1.

Following this strategy of refinement, Table 1 shows the number of elements and degrees of freedom at different refinement levels.

Refinement levels


Domain discretization of channel with pair of circular obstacles in tandem arrangement (several gap spacing) at different refinement levels is provided in Table 1. From the evaluated data about the degrees of freedom at , it is concluded that for the high refinement levels, the degree of freedom is 190551 at 190941 at , and 218349 at . Table 1 also provides the number of domain elements as well as boundary elements for all the generated meshes.

3.3. Quantities of Interest

The important quantities are values of drag and lift which represent the forces in the horizontal and vertical directions, respectively, defined by the surface integration over the cylinders:

These drag and lift forces combine pressure and viscous forces as evident from the above surface integrals. Once the hydrodynamic forces are available, then their nondimensional analogue namely the drag coefficient and lift coefficient are readily available as postprocessing by

Here, is the average velocity, and D is the diameter of the obstacle.

4. Results and Discussions

The incompressible fluid is considered in the present study to explore some new dynamics of hydrodynamic forces. For this purpose, a pair of obstacles of circular shapes is placed in a channel to observe the impact of the Bingham fluid flow. The upstream cylinder is fixed at the position (0.2, 0.2), and the position of downstream cylinder is changed at three different locations to establish the gap spacing . For the fixed value of Reynolds number , the flow of the Bingham fluid is observed. Starting from (the Newtonian case) up to , the drag and lifts are calculated at both the cylinders (see Tables 2 and 3). Fluid velocity, pressure, and viscosity are observed for gap spacing () and for different values of . For all values, it is observed that velocity decreases with the increase in (see Figures 3, 4, and 5(a)5(c)) and plug zone extends from the center of channel towards the solid walls, and shear zone is restricted to the neighborhood of obstacles.

for left cylinder for right cylinder


Left cylinder liftRight cylinder lift


Figures 6, 7, and 8(a)8(c) show the pressure profiles representing the stagnation zones. Stagnation point appears as the fluid hits upstream cylinder. There is lesser impact of pressure on downstream cylinder for (Figure 6(a)). However, it increases with the increase in (Figure 6(c)), and it is also observed that the pressure increases on downstream cylinder with the increase of gap spacing (Figure 8).

The impact of pertinent parameters on viscosity is plotted in Figures 9, 10, and 11. It is observed that viscosity increases with the increase in . Higher values of viscosity observed even between the two cylinders as increases. Moreover, small islands of high viscosity appear around the cylinders for all gap spacing and their size grows with an increase in .

The steady flow in two-dimensional space is considered for analysis. Cut lines are marked exactly at the middle of two cylinders () to observe the velocity impact in the center of the two cylinders, and and components are considered for minute analysis. Plots show the dependence of stream lines on , , and. For , the visual impact of velocity is observed for the component in Figure 12(a). For (non-Newtonian case), the fluidity is much higher than the rest of the . The gradual increase in nonlinearity due to increase in can also be observed. The increase in gap ratio also has an impact on flow pattern in the middle of cylinders that can be seen in Figures 12(b) and 12(c).

The impact of spacing on vertical velocity is negligible; however, for Newtonian case , it gets some peak as described in Figure 13.

For a better analysis, the data presented in Tables 2 and 3 have been plotted in Figure 14 to see the trend of gap ratio versus . A linear increasing profile for drag on both left and right cylinders is observed. Furthermore, the drag coefficient is higher for left cylinder as compared with the right one since the fluid forces are dominant on the upstream cylinder. The lift coefficient varies in a nonlinear fashion for all gap ratios; however, for left cylinder, it increases with an augmentation in whereas for right cylinder, the lift coefficient drops for after a certain threshold of .

A correlation analysis has also been done to check the strength of relationship between drag and lift coefficients at various gap ratios and for all the values of ranging from 0 to 50. The software package SPSS 23.0 is used for this purpose.

In Table 4, the correlation between drag and lift is evaluated for for left and right cylinders, respectively, with different spacing. For left cylinder, the drag and lift are positively correlated, and the relationship is very strong. It gets even stronger as the gap increases. For right cylinder, the drag and lift are inversely related for all gap spacing. For downstream cylinder at , the drag and lifts are almost independent. But as the gap increases, the correlation becomes strong.

Left cylinder (upstream)Right cylinder (downstream)


5. Conclusions

We have simulated the flow of the Bingham fluid around obstacles using the Papanastasiou regularization and the finite element method. The effects of the Bingham number and gap spacing on the drag and lift coefficient of the cylinders have been investigated in detail. It has been observed that the separation plays a key role in determining the drag and lift coefficients, and it influences the shape of unyielded zones near cylinders. An increment in the Bingham number results in the enhancement of plug zone towards the walls of the channel downstream the obstacles. Increasing the separation between the cylinders gives rise to more stagnation pressure on the downstream cylinder. The drag and lift coefficients are positively correlated for the upstream cylinder for all gap spacing while the correlation becomes negative for the downstream cylinder.


:Yield stress
:Plastic viscosity
:Rate of strain tensor
:Shear rate
:Stress tensor
:Gap ratio
:Dimensional velocity vector
:Dimensionless velocity vector
:Reference velocity
:Dimensional pressure
:Dimensionless pressure
:Dimensional stress growth parameter
M:Dimensionless stress growth parameter
:Reynolds number
:Bingham number
#EL:Number of elements
DOF:Number of degrees of freedom
:Drag coefficient
:Lift coefficient

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors are very thankful for the anonymous referees whose suggestions and comments helped to improve the quality of the manuscript.


  1. M. Gad-el-Hak and P. R. Bandyopadhyay, “Reynolds number effects in wall-bounded turbulent flows,” Applied Mechanics Reviews, vol. 47, no. 8, pp. 307–365, 1994. View at: Publisher Site | Google Scholar
  2. P. R. Spalart, “Direct simulation of a turbulent boundary layer up toRθ= 1410,” Journal of Fluid Mechanics, vol. 187, pp. 61–98, 1988. View at: Publisher Site | Google Scholar
  3. O. Darrigol and U. Frisch, “From Newton's mechanics to Euler’s equations,” Physica, vol. 237, no. 14–17, pp. 1855–1869, 2008. View at: Publisher Site | Google Scholar
  4. M. M. Zdravkovich, Flow around Circular Cylinders. Volume 1: Fundamentals, Oxford University Press, New York, 1997.
  5. M. M. Zdravkovich, Flow around Circular Cylinders. Volume 2: Applications, Oxford University Press, New York, 2003.
  6. M. M. Zdravkovich, “Review of flow interference between two circular cylinders in various arrangements,” ASME J. Fluids Eng., vol. 199, pp. 618–633, 1977. View at: Google Scholar
  7. M. M. Zdravkovich, “The effects of interference between circular cylinders in cross flow,” Journal of Fluids and Structures, vol. 1, no. 2, pp. 239–261, 1987. View at: Publisher Site | Google Scholar
  8. R. P. Chhabra, “Hydrodynamics of Non-Spherical Particles in Non-Newtonian Fluids,” in Handbook of Applied Polymer Processing Technology, N. P. Cheremisinoff and P. N. Cheremisinoff, Eds., Marcel Dekker, New York, 1996, Chapter 1. View at: Google Scholar
  9. R. P. Chhabra, “Heat and Mass Transfer in Rheologically Complex Systems,” in AdVances in the Rheology and Flow of Non-Newtonian Fluids, D. Siginer, D. Kee, and R. P. Chhabra, Eds., Elsevier, Amsterdan, 1999, Chapter 39. View at: Google Scholar
  10. R. P. Chhabra, Bubbles, Drops, and Particles in Non-Newtonian Fluids, CRC Press, Boca Raton, FL, 2006. View at: Publisher Site
  11. R. E. Powell and H. Eyring, “Mechanism for relaxation theory of viscosity,” Nature, vol. 154, pp. 422–428, 1944. View at: Publisher Site | Google Scholar
  12. Y. I. Cho and K. R. Kensey, “Effects of the Non-Newtonian Viscosity of Blood on Hemodynamicsof Diseased Arterial Flows,” Advances in Bioengineering, vol. 15, pp. 147-148, 1989. View at: Google Scholar
  13. M. M. Cross, “Rheology of non-Newtonian fluids: A new flow equation for pseudoplastic systems,” Journal of Colloid Science, vol. 20, no. 5, pp. 417–437, 1965. View at: Publisher Site | Google Scholar
  14. H. Steffan, W. Brandstatier, G. Bachler, and R. Pucher, “Comparison of Newtonian and non-Newtonian blood flow in stenotic vessels using numerical simulation,” in Proc. 2nd Int. Symp. On Bioflud Mechanics & Biorheology, pp. 695–701, Munich, Germany, 1989. View at: Google Scholar
  15. R. B. Bird, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids, Fluid Mech., vol. 1, John Wiley and Sons, New York, 1987.
  16. M. Nakamura and T. Sawada, “Numerical study on the flow of a non-Newtonian fluid through an axisymmetric stenosis,” Journal of Biomechanical Engineering, vol. 110, no. 2, pp. 137–143, 1988. View at: Publisher Site | Google Scholar
  17. Y. I. Cho and K. R. Kensey, “Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: steady flows,” Biorheology, vol. 28, no. 3-4, pp. 241–262, 1991. View at: Publisher Site | Google Scholar
  18. E. C. Bingham, An Investigation of the Laws of Plastic Flow, US Government Printing Office, Pueblo, CO, USA, 1917.
  19. W. Herschel and R. Bulkley, “Konsistenzmessungen von Gummi-Benzollösungen,” Kolloid Zeitschrift, vol. 39, no. 4, pp. 291–300, 1926. View at: Publisher Site | Google Scholar
  20. T. C. Papanastasiou, “Flows of materials with yield,” Journal of Rheology, vol. 31, no. 5, pp. 385–404, 1987. View at: Publisher Site | Google Scholar
  21. D. D. KÉ and G. Turcotte, “Viscosity of biomaterials,” Chemical Engineering Communications, vol. 6, no. 4-5, pp. 273–282, 1980. View at: Publisher Site | Google Scholar
  22. R. Mahmood, N. Kousar, M. Yaqub, and K. Jabeen, “Numerical simulations of the square lid driven cavity flow of Bingham fluids using nonconforming finite elements coupled with a direct solver,” Adv. Math. Phys., vol. 2017, pp. 1–10, 2017. View at: Publisher Site | Google Scholar
  23. R. Mahmood, N. Kousar, K. Usman, and A. Mehmood, “Finite element simulations for stationary Bingham fluid flow past a circular cylinder,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, vol. 40, no. 9, 2018. View at: Publisher Site | Google Scholar
  24. A. Mehmood, W. A. Khan, R. Mahmood, and K. U. Rehman, “Finite element analysis on Bingham-Papanastasuou viscoplastic flow in a channel with circular/square obstacles: a comparative benchmarking,” PRO, vol. 7, p. 779, 2020. View at: Google Scholar
  25. H. Barnes and K. Walters, “The yield stress myth?” Rheologica Acta, vol. 24, no. 4, pp. 323–326, 1985. View at: Publisher Site | Google Scholar
  26. H. A. Barnes, “The yield stress- a review or “” everything flows?” Journal of Non-Newtonian Fluid Mechanics, vol. 81, pp. 133–178, 1999. View at: Google Scholar
  27. R. C. Gupta, “Developing Bingham fluid flow in a channel,” Mathematical and computer modeling, vol. 21, no. 8, pp. 21–28, 1995. View at: Publisher Site | Google Scholar
  28. A. Beris, J. Tsamopoulos, R. Armstrong, and R. Brown, “Creeping motion of a sphere through a Bingham plastic,” Journal of Fluid Mechanics, vol. 158, pp. 219–244, 1985. View at: Publisher Site | Google Scholar
  29. K. Adachi and N. Yoshioka, “On creeping flow of a visco-plastic fluid past a circular cylinder,” Chemical Engineering Science, vol. 28, no. 1, pp. 215–226, 1973. View at: Publisher Site | Google Scholar
  30. E. Mitsoulis, “On creeping drag flow of a viscoplastic fluid past a circular cylinder: wall effects,” Chemical Engineering Science, vol. 59, no. 4, pp. 789–800, 2004. View at: Publisher Site | Google Scholar
  31. A. Vakil and S. I. Green, “Two-dimensional side-by-side circular cylinders at moderate Reynolds numbers,” Computers and Fluids, vol. 51, no. 1, pp. 136–144, 2011. View at: Publisher Site | Google Scholar
  32. L. Jossic and A. Magnin, “Drag of an isolated cylinder and interactions between two cylinders in yield stress fluids,” Journal of Non-Newtonian Fluid Mechanics, vol. 164, no. 1-3, pp. 9–16, 2009. View at: Publisher Site | Google Scholar
  33. A. R. Koblitz, S. Lovett, and N. Nikiforakis, “Viscoplastic squeeze flow between two identical infinite circular cylinders,” Phys. Rev. Fluids, vol. 3, no. 2, 2018. View at: Publisher Site | Google Scholar
  34. J. Calder and A. Yezzi, “PDE acceleration: a convergence rate analysis and applications to obstacle problems,” Res. Math. Sci., vol. 6, no. 4, pp. 1–30, 2019. View at: Publisher Site | Google Scholar
  35. A. Mahammedi, H. Ameur, and A. Ariss, “Numerical investigation of the performance of kenics static mixers for the agitation of shear thinning fluids,” Journal of Applied Fluid Mechanics, vol. 10, no. 3, pp. 989–999, 2017. View at: Publisher Site | Google Scholar
  36. H. Ameur and Y. Menni, “Non-Newtonian fluid flows through backward facing steps,” S.N. Applied sciences, vol. 1, no. 12, 2019. View at: Publisher Site | Google Scholar
  37. H. Ameur, “Pressure drop and vortex size of power law fluids flow in branching channels with sudden expansion,” Journal of Applied Fluid Mechanics, vol. 11, no. 6, pp. 1739–1749, 2018. View at: Publisher Site | Google Scholar
  38. H. Ameur, “3D hydrodynamics involving multiple eccentric impellers in unbaffled cylindrical tank,” Chinese Journal of Chemical Engineering, vol. 24, no. 5, pp. 572–580, 2016. View at: Publisher Site | Google Scholar
  39. H. Laidoudi, “Upward flow and heat transfer around two heated circular cylinders in SQUARE duct under aiding thermal buoyancy,” Journal of the Serbian Society for Computational Mechanics, vol. 14, no. 1, pp. 113–123, 2020. View at: Publisher Site | Google Scholar
  40. H. Laidoudi, “Mixed convection heat transfer from confined tandem circular cylinders in cross flow at low Reynolds number,” Mechanica, vol. 23, 2017. View at: Google Scholar
  41. H. Laidoudi and M. Bouzit, “Mixed convection in Poiseuille fluid an asymmetrically confined heated circular cylinder,” Thermal Science, vol. 22, no. 2, pp. 821–834, 2016. View at: Publisher Site | Google Scholar
  42. M. Ferronato, “Preconditioning for sparse linear systems at the down of the 21st Century, History, Current Developments, and Future Perspectives,” Int. Sch. Research Net., vol. 49, 2012. View at: Google Scholar

Copyright © 2021 Asif Mehmood et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.