#### Abstract

This study focuses on the industrial and engineering interest quantities, such as drag force and rate of transmission of heat, for pseudo-plastic nanofluid flow. The attributes of natural convection of the pseudo-plastic nanofluid flow model over a vertical slender cylinder are explored. The pseudo-plastic flow is studied under the influence of concentration of nanoparticles, rate of heat transmission, and drag force. For the first time, the pseudo-plastic nanofluid flow model has been implemented over a vertical slender cylinder which is not yet investigated. The acquired model is based on thermophoresis and Brownian motion mechanisms. The governing equations of pseudo-plastic nanofluid in cylindrical coordinates are modelled. The developed system of nonlinear equations is tackled by boundary layer assumptions and similarity transformations. Moreover, the solution of the acquired system exhibited by using a new powerful numerical technique. A comprehensive debate on drag force and transmission of heat under the influence of various emerging parameters is illustrated in the table. Furthermore, the effects of dimensionless parameters over the velocity profile, temperature profile, and concentration of nanoparticle profile have been exhibited graphically.

#### 1. Introduction

Investigation of the non-Newtonian fluids gains prodigious attention of researchers over a half-century because naturally, most of the fluids used in the industrial applications are non-Newtonian fluids. This is the main cause of increased applications of non-Newtonian fluids in the industrial field and engineering such as petroleum production, molten plastics, food engineering, automobiles, polymer solution industry, chemical engineering, and power engineering. A solitary established equations cannot pronounce the attributes of such non-Newtonian fluids because these kinds of fluids have a nonlinear relationship among the rate of stress and strain. Therefore, several scientists and engineers have pronounced models for non-Newtonian fluids [1–10]. Ellahi et al. [11] explored thermally charged MHD biphase flow coating with non-Newtonian nanofluid over the slippery walls. The sustainable features of MHD Jaffrey fluid for bio-bi-phase flow are carried out by Zeeshan et al. [12]. Rehman et al. [13] examined MHD flow of the nanoparticle influenced near a stagnation point over an exponential stretched surface. Parand et al. [14] described the boundary layer flow of Powell–Eyring fluid for a stretching sheet. Generation/absorption of heat through the flow of axisymmetric Casson fluid over a stretched cylinder is addressed by Javed et al. [15]. Additional appropriate studies in this way are in [16, 17].

Over the last few years, we have adopted a pioneering procedure for refining the transmission of heat by utilizing ultrafine solid particles in the fluids, and these particles have been used widely called as nanofluid. The label “nanofluid” was used for the first time in 1995 by Choi and Eastman [18]. He cited that it is conceivable to boost convection of transmission of heat and thermal conductivity efficiency by using nanofluids. To handle nanofluids, Buongiorno [19] introduced a mathematical model and explored numerous transport mechanisms and applications about the nanofluids. Miscellaneous purposes of nanofluids have been uncovered by Das et al. [20]. Zhang et al. [21] scrutinized the nanofluid for MHD radiative flow over a variable heat flux and chemical reaction surface. Transport of heat for the ferromagnetic fluid with thermal stratification over a stretching sheet was examined by Muhammad et al. [22]. Saini and Sharma [23] extended the application of nanofluids through the investigation of double-diffusive bioconvection. Further extensive research and applications of nanofluids across numerous fields are heat transmission, energy, microequipment, and boiling applications [24–28].

Williamson’s fluid model has non-Newtonian behaviour in nature like pseudo-plastic fluid which defines the flow of shear thinning. The industrial, engineering, and biological fluids which observe as Williamson’s fluid are blood, glue, paint, ketchup, polymer solutions, and nail polish. Williamson [29] who discovered the model to communicate pseudo-plastic attributes along with features of extreme points of viscosity. Due to this invention, innovated researchers are motivated to discover more upfront classifications of non-Newtonian fluid. Ramzan et al. [30] analyzed the numerical solution of MHD flow over a stretched surface with convective boundary conditions using the shooting method. Two-dimensional flow of Williamson’s fluid film with heat diffusion under the influence of thermal radiation was inspected via an optimal approach by Shah et al. [31].

Rate of transmission of heat, the transmission of mass rate, and skin friction coefficient play a dynamic character during coating of wires or polymer fibre coating. As wires have a thin cylindrical shape, miscellaneous researchers have launched several mathematical models. In axial incompressible flow, Seban and Bond [32] premeditated the attributes of drag force and rate of transmission of heat for a cylinder in 1951. Under the uniform surface heat flux, Mucoglu and Chen [33], for the first time, scrutinized the slender cylinder with the help of mixed convection regime. Nadeem et al. [34] analyzed the boundary layer flow and transmission of heat of a nanofluid in a vertical slender cylinder. Patil et al. [35] examined the mixed convection nanofluid boundary layer flow under the effect of surface roughness with a moving slender cylinder. Reddy et al. [36] explored the natural convection for the Casson fluid flow past over a hollow slender cylinder. With the help of Bejan’s heat function concept, Reddy et al. [37] investigated the unsteady MHD micropolar fluid flow in a homogeneously thermal radiative hollow slender cylinder with the radiative transmission of heat effect. Further latest stimulating work in this area can be found in [38–42].

The literature review replicates that, generally, the researchers engaged themselves to study non-Newtonian fluids’ behaviour by assuming different effects. The attributes of natural convection of Williamson’s nanofluid model over a vertical slender cylinder are not explored until now. Therefore, the developed model is simplified by the boundary layer approximation and similarity transformations. The governing coupled nonlinear system of equations is then solved by a new powerful numerical technique. Furthermore, physical behaviour for the industrial interests of the fluid will be examined through the table.

#### 2. Analysis of Flow and Mathematical Formulation

Suppose an incompressible Williamson’s nanofluid flow along with a permeable vertical slender cylinder with radius *a* having uniform ambient temperature . The coordinates are used, whereas *x* is acting along the surface of the cylinder and *r* along the radial direction. The velocity profile, temperature profile, and concentration profile are as follows:where are velocity components along the surface and radial direction. The boundary layer equations of motion, energy, and the nanoparticle concentration are

The accompanied boundary conditions are

Here, is the mainstream velocity, is called kinematic viscosity and is defined as , articulated the material constant for Williamson’s fluid, and denotes the density of the fluid. Now, the nondimensional variables and similarity transformations are defined as follows:in which and operated for the characteristic temperature and nanoparticle concentration . By using the above transformations, equation (2) is identically satisfied, and equations (3)–(5) will be articulated as

Dimensionless attached boundary conditions arewhere illustrates the curvature parameter, symbolizes Williamson’s dimensionless parameter, pronounces as the parameter of natural convection, means the buoyancy ratio, is identified as the Prandtl number, stands for the Brownian motion parameter, corresponds to the thermophoresis parameter, and is known as the Lewis number.

#### 3. Physical Quantities of the Industrial Interest

For the industrial interest, the physical quantities, i.e., drag force and transmission of heat, are stated aswhere stands for the shear stress tensor over the surface of the slender cylinder, while is called wall heat flux. These physical quantities can be articulated as

According to similarity transformation, equation (15) is transformed as

#### 4. Numerical Solution

The solution of the nonlinear system of ODEs (10)–(12) over the accompanied condition equations (13) and (14) is tackled through the numerical algorithm of MATLAB inherent scheme bvp4c. In order to apply this scheme, first of all, higher-order differential equations are transformed into first-order ODEs. The procedure is as follows.

Let us suppose that

Therefore, equations (10)–(14) are transformed as follows:with conditions

#### 5. Discussion

Natural convection of Williamson’s nanofluid over a slender cylinder has been analyzed. Numerical solution is obtained by using MATLAB scheme bvp4c. Physical properties of the nanofluid are portrayed by using the Buongiorno model. Likewise, transmission of heat and mass is differentiated by means of assorted parameters.

##### 5.1. Velocity Profiles

###### 5.1.1. Attributes of and for

Figures (1) and (2) show the attributes of curvature parameter and Williamson’s dimensionless parameter for . It is well defined from Figures (1) and (2) that intensifying the values of tends to devalue ; however, equivalent behaviour is anticipated for Williamson’s parameter . This is because after intensifying the curvature parameter , the radius of the cylinder with the fluid declines. Similarly, after intensifying Williamson’s parameter which will cause a decline in velocity because fluid opposes more resistance.

###### 5.1.2. Attributes of and for

The attribution of and is shown in Figures (3) and (4). From Figures (3) and (4), it defines that, by intensifying natural convection and buoyancy ratio , the buoyancy force will cause the higher velocity attained by the fluid.

##### 5.2. Temperature Profiles

###### 5.2.1. Attributes of and for

The behaviour of temperature profile for distinct values of curvature parameter and Prandtl number is reported in Figures (5) and (6). Evidently, from Figure 5, it can be described that enhancement in the curvature parameter yields intensifying in . Moreover, in Figure 6, the temperature profile and the boundary layer thickness decrease due to enhancement in Prandtl number (Pr). Which uncovers the truth that, enlarging in Pr cause the reduction in the thermal diffusivity of the fluids accordingly.

###### 5.2.2. Attributes of and for

The contribution of thermophoresis parameter and Brownian motion is described in Figures (7) and (8). Evidently, larger and produce higher . Practically, movement of a small number of particles from higher temperature to the lower one is defined as the thermophoresis phenomenon. Hence, a greater number of nanoparticles are shifted from the hot region which raises the fluid temperature. However, due to the increase in Brownian motion parameter in the result, random motion of the nanoparticles will be raised which causes intensifying in the temperature of the fluid.

##### 5.3. Concentration Profiles

###### 5.3.1. Attributes of on

Figure 9 exhibits the consequences of curvature parameter on nanoparticle concentration . It is clear from the figure that intensifying curvature parameter diminishes the nanoparticle concentration of the fluid. According to the industrial view, after increase in the curvature parameter, it caused the reduction in the radius of the slender cylinder; hence, concentration of nanoparticles also reduced.

###### 5.3.2. Attributes of Lewis Number Le and Prandtl Number on

The behaviour of Lewis number Le is depicted for through Figure 10. It is scrutinized that is diminished by increasing Lewis number Le. Physically, due to the levitation of diffusivity of heat and mass, nanoparticle concentration diminished. Furthermore, scrutinization of the attribution of Prandtl number Pr for is portrayed through Figure 11. Physically, intensifying in Prandtl number will cause the more heat convection occur because thickness of the thermal boundary layer is greater than the velocity boundary layer. Hence, nanoparticle concentration diminished.

###### 5.3.3. Attributes of and on

The impact of thermophoresis parameter and Brownian motion parameter on is depicted in Figures (12) and (13). Due to the increase of , boundary layer thickness of nanoparticle concentration intensified. As a matter of fact, thermophoresis force increases for higher estimation of , which yields nanoparticle relocations from higher to lower temperature; hence, intensifies (see Figure 12). However, and its concentration boundary layer diminish when is amplified (see Figure 13).

##### 5.4. Attributes of Drag Force and Rate of Transmission of Heat

The attributes of physical parameters such as curvature parameter , Williamson’s parameter , natural convection , and buoyancy ratio on drag force are exhibited in Table 1. It is evidently clear from the table that drag force is intensified by enhancing , , , and . Table 2 expresses the behaviour of physical parameters on the rate of transmission of heat. It is inspected that the rate of transmission of heat enhances by increasing , , and , while it declines for inclination of , , and .

#### 6. Final Remarks

This exploration scrutinized the natural convection of Williamson’s nanofluid for boundary-layer stagnation point flow over a vertical slender cylinder. Hence, thermophoresis, Brownian motion, natural convection, Williamson’s parameter, and curvature parameter are utilized for modelling and analysis. The governing coupled nonlinear ODEs are then solved numerically by using MATLAB technique bvp4c. This research reveals the following outcomes:(i)The velocity profile is declining due to upsurge in the curvature parameter and Williamson’s parameter(ii)Inclination of the Prandtl number and curvature parameter causes decline in the temperature profile, whereas it is intensified by increments in , , and (iii)The concentration of nanoparticles is dropped via amplifying , , and (iv)The greater value of Williamson’s parameter and curvature parameter causes the drag force to increase(v)Rate of transmission of heat is reduced for larger Williamson’s parameter

#### Data Availability

The data used to support the findings of this study are included within the article.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.