Abstract

In this work, a comprehensive study of fluid forces and thermal analysis of two-dimensional, laminar, and incompressible complex (power law, Bingham, and Herschel–Bulkley) fluid flow over a topological cross-sectional cylinder (square, hexagon, and circle) in channel have been computationally done by using finite element technique. The characteristics of nonlinear flow for varying ranges of power law index Bingham number , Prandtl number , Reynolds number , and Grashof number have been examined. Considerable evaluation for thermal flow field in the form of dimensionless velocity profile, isotherms, drag and lift coefficients, and average Nusselt number is done. Also, for a range of , the drag forces reduction is observed for circular and hexagonal obstacles in comparison with the square cylinder. At corresponding to Newtonian fluid, maximum reduction in drag force is reported.

1. Introduction

Since the last several decades, heat transfer from a single cylinder (circular, square) has been widely considered, given its significance in various fields of engineering mathematical problems. The key research factors for the design, maintenance of electric cooling system, evaporators, heat exchangers, thermal plants and automobile radiators, etc., can be reported. Moreover, thermal flow by natural, forced, and mixed convection from solid body to fluid stream has been simulated by many researchers in a variety of configurations. Ditchfield et al. [1] experimentally examined the effects of thermal flow on viscoplastic fluid by implementing a proposed model. Shyam et al. [2] considered the influence of nonisothermal and nonlinear viscous fluid flow over the pair of circular heaters at large Grashof number. Computations were reported by varying gap between cylinders to elucidate its impact on thermal flow field. Laidoudi et al. [3, 4] analyzed the characteristics of mixed convection thermal flow of dual circular cylinders at low Reynolds number. They also computed fluid forces and average Nusselt number and observed high Nusselt number on second obstacle compared to the first one due to thermal buoyancy. Kefayati et al. [5] have numerically investigated the effects of double diffusion and entropy generation of viscous fluid flow over a cold cylinder in a duct by Lattice Boltzmann technique. In addition, an enhancement in heat transfer was noticed due to increase in the length of diameter of the cylinder. The buoyancy effects on thermal flow characteristics of confined circular cylinder submerged laminar Poiseuille fluid are studied by Laidoudi et al. [6]. Masoumi et al. [7] have examined the influence of nonisothermal viscoplastic fluid flow inside the circular duck and also imposed the Galerkins Weighted Residual (GWR) scheme of FEM computation with hybrid mesh. The influence of Richardson number on flow and heat transfer topologies of an obstacle at high Reynolds number are examined by Zafar and Alam [8]. Laidoudi et al. [9] considered the effects of the free convection fluid flow between two cylinders at high Reynolds number in a circular duct and also provided the correlation impacts of the average Nusselt number.

In order to design structures for non-Newtonian fluid exposed in polymer processing activities, reliable knowledge about fluid forces is often necessary, and in this regard, the flow and heat transfer in the free convection regime have been considerably studied in extensive research resources in recent years. Bharti et al. [10, 11] investigated the characteristics of steady, incompressible, non-Newtonian fluid flow over a confined cylinder at low Reynolds number . They solved the model by using QUICK algorithm and finite volume scheme. Shyam et al. [12] numerically studied the influence of thermal flow of free convection over a constant heated block in a square duck and also reported the effects of average Nusselt number on Grashof and Prandtl number against power law index. Baranwal and Chhabra [13] considered an incompressible free convective thermal flow of power law fluid over a confined cylinder in square cavity. In addition, they examined the geometric effects of the cylinder on thermal flow fields. A. Pantokratoras [14] numerically reported the flow of power law fluid over a rotating cylinder solved by using finite volume scheme through commercial code FLUENT. Gangawane and Manikandan [15] have explored the thermal flow features over a hexagonal block with constant temperature, solved by using finite volume technique with QUICK discretization and SIMPLE algorithm. Asnaashari and Tohidi [16] investigated the effects of flow and heat transfer of power law fluid in an unconfined computational domain. Dhiman and Shyam [17] have examined the influence of Reynolds number on unsteady heat transfer over an equilateral triangular block and created grid by GAMBIT and solved it by using FLUENT.

Nirmalkar and Chhabra [18] studied the characteristics of Bingham fluid flow over a heated submerged cylinder using finite element computation. Patel and Chhabra [19] have extensively examined heat transfer, fluid forces, and yield stress results in steady flow regime. Mahmood et al. [20] have simulated a plastic fluid flow of square cavity by using the open-source software package FEATFLOW submerged in Papanastasiou regularization. Liu et al. [21] investigated the characteristics of Bingham fluid flow in deeply buried rock. Thakur et al. [22] have been considered a forced convective heat transport of a Bingham fluid over a rotating circular cylinder at low Reynolds number. Rem et al. [23] investigated the influence of unsteady Bingham fluid flow in a fractured channel.

The Herschel–Bulkley fluid model narrates the rheological influence of viscoplastic fluids which occurred in innumerable applications in industry. Theses viscous fluids are characterized by yield stress which requires a finite stress for flow. The relation between the shear stress and the shear rate is nonlinear. Since thermal flow is important in industrial applications with viscoplastic fluid flow like paint, cosmetic, and petroleum, many researchers have examined the characteristics of thermal flow in such viscoplastic fluid flow regimes [2432].

Zdanski and Vaz [33] have implemented a second-order finite difference scheme to simulate a flow problem of incompressible 3D viscous fluid flow in a channel. Abbasi et al. [34] have numerically analyzed the impacts of hydrodynamic forces over a cylinder-based Lattice Boltzmann Method (LBM) at low Reynolds number. Also, they compared their results with numerical and experimental data for single square cylinder case. Mahmood et al. [35] studied the non-Newtonian flow in a channel driven cavity by using element pair for finite element computation. Abbasi et al. [36] examined the characteristics of fluid forces over square cylinders at varying range of gap spacing, solved by LBM computation. Khan et al. [37] imposed least square FEM computation of viscous fluid flow over a semicylinder block via COMSOL solver. Tomio et al. [38] implemented a numerical methodology, second-order finite difference scheme to solve a flow problem. Characterization of fluid forces of nonlinear material in a channel is based on FEM computation at low Reynolds number [3941].

This manuscript is organized as follows. in Section 1 and Section 2, the introduction of the problem and the flow configuration and governing equations with constitutive relations are explained. Numerical scheme and grid convergence are described in Section 3. Also, the results are compared with literature for code validation. In Section 4, we explain in detail all the results and discussion of the article. Finally, conclusion is mentioned in Section 5.

2. Flow Configuration and Governing Equations

A cylinder (square, hexagon, and circle) of diameter in an incompressible viscous polymer solution in laminar, steady, 2D flow regime is revealed in Figure 1. The surface of cylinders is maintained at a uniform temperature and flow with average velocity . The total length of given benchmark problem is where is the upstream length and downstream length is from the center of the cylinder and height of the channel is 0.41, respectively. Moreover, the thermophysical properties of fluids are considered to be free from temperature. Also, we have neglected the effects of viscous dissipation from heat equation.

The geometrical configuration of the problem is shown in Figure 1.

Continuity, momentum, and energy are the governing equations of steady, nonisothermal, and viscous incompressible fluid flow across an infinitely long channel in presence of heated cylinder. Their nondimensional forms of these equations are written as follows:where represents the category of the chosen fluid. Here we have used three models, namely, power law, Bingham, and Herschel–Bulkley fluid defined as follows.(i)Power law model:(ii)Bingham model:(iii)Herschel–Bulkley model:

The nondimensionalized equations (1)–(3) suggest that velocity, pressure, and temperature fields are dependent on dimensionless parameters, Reynolds number , Grashof number , Prandtl number , and Bingham number . These parameters are defined aswhere symbols have been defined in nomenclature. To compute the solution of equations (1)–(3), the system must be subject to some boundary conditions. At the inlet of this domain, fluid enters with parabolic profile in x-direction. At the outlet, the Neumann conditions for all variables have been imposed. On the surface of channel and obstacles there is no-slip boundary condition. Furthermore, fins are insulated and obstacles are heated uniformly while other walls of channel remain cold.

Ultimately, computational results from the preceding governing equations in the form of primitive variables , together with the above boundary conditions of the domain are postprocessed to evaluate the derived quantities like drag coefficient, lift coefficient, and , respectively.

The nondimensional net drag and lift forces act on the submerged cylinder in the flow and normal direction, while the drag and lift coefficient are

Here the reference velocity is where is the maximum velocity of given parabolic profile at the inlet. The local Nusselt number on the surface of duct and fins is estimated by . Such values have been more averaged over both surfaces (cylinder and fins) to obtain the average Nusselt number given aswhere are the surfaces of thermal region and normal direction of the surface; it is reasonable to postulate that the drag coefficient is function of whereas the depends on the . This work endeavors to exhibit and develop more efficient functional relationship for circular cylinder in future.

3. Numerical Scheme and Grid Convergence

Due to nonlinearity of the governing equations as well as of non-Newtonian fluids, having plastic viscosity applied to channel flow, the energy and momentum equations cannot be supported with analytical solutions, so we implemented numerical scheme finite element method (FEM) for computing the optimized solution with different shapes of cylinder. The underlying discrete nonlinear system of equations has been solved by Newton method and the linearized inner system is solved with a direct solver PARDISO. The PARDISO solver utilizes LU matrix factorization and reduces the number of iterations required for the desired level of convergence. This reduces the number of the iterations required. The Lagrangian system was used for the velocity discretization (second order) of elements and pressures and temperatures (linear) as is shown in Figure 2.

Grid independence study has been performed for the , , and for the different cylinders and corresponding results are given in Table 1–Table 3 for all the three shapes of cylinders at , and . The benchmark results reveal the fact that results at levels 7 and 8 are in close agreement so to save computational cost, and we performed the rest of the simulation at level 7.

Finally, the current code has also been validated against the Schaefer and Turek [42] in Table 4 by setting the temperature on circular cylinder as zero.

4. Results and Discussions

This work endeavors in three different cross-sectional cylinders with significant boundary conditions to elucidate the impact of non-Newtonian fluid flow and heat transfer behavior. The results shown in this study cover the range of parameters as: , and respectively. Here is a comment on the choice of these ranges. The value of the power law index is encountered in many products such as suspension and polymer fluids found in industrial practices. Moreover, are considered adequate for determining their influence on flow and heat transfer characteristics. The lowest value of is Newtonian and the fluid in the very wide is plug-like and heat is transmitted only via conduction which is slightly increased by piston-like fluid motion. The selection of the range of parameters for present work is based on a combination of the actual material characteristics and the anticipated fluid limiting conditions. Extensive results are discussed in the ensuring sections on lengths entrance, velocity profile and temperature, yield surfaces, and . It is, however, edifying first to determine the accuracy of the existing numeric assumptions and this objective is accomplished by presenting a number of benchmarks for a number of limiting circumstances and comparing the results available in literature.

4.1. Power Law Model

The velocity distribution illuminated in Figure 3 with varying cylinder shapes at Reynolds with Newtonian fluid is restricted. Since the velocity of the parabola is induced at inlet, and other boundaries are kept in no-slip conditions, velocity variation close to channel obstacles is observed. A conclusion is reached that the parabolic nature of velocity profiles is similar in the entry. Figure 4 is plotted to evaluate pressure variation across the physical domain, particularly in vicinity of circular obstacle and constant . From the figure, it is seen that pressure exhibits optimized nonlinear (power law, Bingham, and Herschel–Bulkley) fluids behavior close to the obstacle and gets linear downstream in channel flow as expected. Here, the optimum pressure is observed to have maximum value in front of the obstruction that interacts with fluid. It also appears that its pressures become linear when the fluid passes over the obstacle.

Figure 5 shows the effects of the shape of cylinder on and against from 0.4 to 1.6 covering a wide range of Newtonian, shear-thinning (thixotropic), and shear-thickening (rheopectic) fluids. The reductions in the values of and are observed from square to circle cylinders. Moreover, reduction in power law index results in lowering the apparent fluid’s viscosity beside the cylinder, thereby reducing the frictional drag strength. The shape effects for Newtonian fluid against are illustrated in Figure 6. Furthermore, the power law index influences the total drag coefficient which is decreasing gradually as the Re increases from 10 to 50. In Figure 7, the effect of drag and lift coefficients on obstacle shapes was shown; with increasing Prandtl numbers, and decrease. Both coefficients produce lower values for the case of circular cylinder as compared with the other two shapes.

The changes in the flow field because of the parametric values of , Re, and significantly impact the over the surface of cylinders. Towards this end, we presented in this section the effects on the for the different shapes of the cylinder for , , and . Figure 8(a) shows that increasing the gives rise to an increment in the . The dependence of on and for three different cylinders is shown in Figures 8(b) and 8(c). increases as the increases because the thermal boundary layer is progressively thinning. In addition, an increase in shear-thinning (decreased value of n) also improves the rate of thermal flow, due to the reduction of the apparent fluid viscosity because of steep gradients close to the submerged cylinders.

4.2. Bingham Model

In this case, the flow field comprises the yielded and unyielded zones, which vary with the choice of Bingham fluid which is nondimensional analog of yield stress. Figure 9 illustrates the areas where fluid moves like a rigid body in unyielding (plug) zones. The viscosity increases as the yield stress in the Bingham fluid gradually increases. In the vicinity of the submerged cylinder there is always a yielded region. Bingham fluid acts in the center of the channel like a solid. Figure 10 shows the isotherms around the cylinders for the Newtonian case at fixed Re. The isothermal profile reflects the physical phenomena observed through the analysis of patterns of streamlines. More thermal contours are crowded near the curved surface than on the flat surface, which leads to a higher on the curved surface. The crowds of temperature contours depend on the and

The trend of and based on the shape of cylinders versus Bn can be seen in Figure 11. There is a linear growing profile for drag coefficient for all shapes. However, for a square cylinder the is higher as compared with the other cases because the fluid forces dominate at the surface of square cylinder. The varies nonlinearly for all cylinders; for the square cylinder, however, the decreases after a certain threshold with an increase in the . However, quantitative analysis of and for different is presented in Table 5. Table 6 shows that, in the case of three shaped cylinders, maximum reduction is observed for all values. In addition, the maximum reduction can be seen at for all cylinder combinations. In the range , for all configurations that approach their corresponding minimum values, the values of a percent reduction decrease.

Figure 12 reveals the disparity of the based on the shapes of cylinder and with the Bn. The  has been shown to dominate in case of square cylinder. 

4.3. Herschel–Bulkley Model

The fluid model Herschel–Bulkley combines the characteristics of Power Law and Bingham models. The fluid flow of Herschel–Bulkley thus shares the features of power law and of Bingham fluid. The effect on the flow components is similar from the power index and Bingham number. For this reason, only those flow characteristics over the cylinder that were not described above will be discussed in this case. Such as the Bingham fluid, a fluid region is formed around the cylinder when there is a flow over the cylinder, while the medium is a rigid substance on the periphery. The fluid region is mainly based on the number of Bingham: the larger the amount of Bingham, the lesser the region of fluid.

The size depends simultaneously on the power index: the lower the, the smaller the fluid region. It can thus be concluded that a perturbation of the cylinders is reduced if the rheology of the fluids is non-Newtonian, if the number of the Bingham and/or the power index is lower. In strongly non-Newtonian fluids, dependence on the drag and kinematic features of the flow is less significant than in Newtonian fluids. The responsiveness of the flow properties to relative velocity changes is weakened as the non-Newtonian fluid properties increase. Figure 13 reveals the characteristics of drag and lift coefficient against Bingham number for different fluid properties over the different cylinders. The drag trend for Bn is the same across the fluid, but the drag values for the rheopectic fluid are higher as the boundaries are nonsmooth compared to the thixotropic fluid. The lift values demonstrate an increasing trend for from 0 to 10 and reveal a decreasing trend for a further increase in the Bingham number for each fluid case. Figure 14 reveals the effects of cylinder shapes on average Nusselt number against Bn for complex fluids. Figure 15 shows the variation of the with on both shape cylinders and the. The increases with as usual.

A line graph is shown in Figure 16 to depict the impact of rheology of fluid on velocity profiles. The velocity profile changes as the yield stress increases, from a parabola to a flat plateau in the channel central region, which represents the unyielded material.

5. Conclusions

In this manuscript, extensive computations have been carried out to elucidate the thermal flow features based on the topology of the obstacles. The results revealed herein enhance over the following varying ranges of nondimensional parameters: from 0.4 to 1.6, from 0 to 50, from 10 to 50, from 0.7 to 10, and from 1 to 10, respectively. In nonlinear fluid, decreases from square to circle by order to optimize cylinder. The results for the development and completion of the velocity profile and thermal flow and for three typical shapes of partially insulated cylinders are detailed. The main findings include the following:(i)Drag and lift coefficients increase against and while they decrease along over all the topologies of cylinders(ii)For the Newtonian case, reduction in is maximum while reduction is decreasing for increasing the (iii)The rate of thermal flow tends to increase as the Re values increase due to a decline in the (iv)The increases with increasing both and (v)The velocity profile shifts from a parabolic to a smooth plateau in the center along the channel when the yield stress increases indicating a plug zone in the center of the channel

Nomenclature

:Dimensionless velocity
:Dimensionless temperature
:Dimensionless pressure
m:Dimensional stress growth parameter
n:Flow behavior index
:Density
:Gravitational acceleration
:Reference length

Greek

:Yield stress
:Plastic viscosity
:Stress tensor
:Shear rate
:Reference velocity
:Reynolds number
:Bingham number
:Prandtl number
:Grashof number
:Drag force
:Drag coefficient
:Lift force
:Lift coefficient
:Average Nusselt number
:Local Nusselt number
:Gradient operator
:Thermal diffusivity
:Kinematic viscosity

Abbreviations

#EL:Number of elements
#DOF:Number of degrees of freedom.

Data Availability

The data underlying the results presented in the study are available within the article.

Conflicts of Interest

The authors state that there are no conflicts of interest.

Acknowledgments

This work was supported by HEC-Pakistan through NRPU Project No. 14038.