# A Numerical Approach for an Unsteady Tangent Hyperbolic Nanofluid Flow past a Wedge in the Presence of Suction/Injection

Authors: U. Shahzad; M. Mushtaq; S. Farid; K. Jabeen; R.M.A. Muntazir
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

The analysis of unsteady tangent hyperbolic nanofluid flow past a wedge with injection-suction, because of its beneficial uses, has gained a lot of attention. The present study is mainly concerned with tangent hyperbolic nanofluid (non-Newtonian nanofluid). First, we have converted the system of partial differential equations (PDEs) to a system of ordinary differential equations (ODEs) with the help of appropriate similarity transformations. Boundary conditions are also transformed by utilizing suitable similarity transformation. Now, for the obtained ODEs, we have used the numerical techniquebvp4c and investigated the velocity, temperature, and concentration profiles. The accuracy of the flow model is validated by applying MAPLE d-solve command having good agreement while comparing the numerical results obtained by bvp4c for both suction and injection cases. The effects of distinct dimensionless parameters on the various profiles are being analyzed. The novel features such as thermophoresis and Brownian motion are also discussed to investigate the characteristics of heat and mass transfer. Graphical representation of the impact of varying parameters and the solution method for the abovementioned model is thoroughly discussed. It was observed that suction or injection can play a key role in controlling boundary layer flow and brings stability in the flow. It was also noticed that by increasing the Darcy number, velocity profile increases in both injection-suction cases.


---

## Body

## 1. Introduction

The heat transfer plays an important role in manufacturing industries. It has effect on cost and quantity of production and is also a challenge for effective heating and cooling effects in high-tech manufacturing industries and many other big industries. A conventional heat transfer fluid has poor thermal conductivity as compared with the solids like water and engine oils [1]. In all of these, thermal conductivity of heat transfer fluids plays a vital role in heat development equipment’s. Choi [2] firstly gave the concept of nanofluids; nanofluids are the fluids whose particles size is less than 100 nm. In many previous years, researchers gave more attention to the nanofluids because of various reasons, for example, heat transfer performance and enhancement of thermal conductivity. Buogiorno [3] explained the extraordinary thermal conductivity of nanofluids with boundary layer (BL) mass and heat transfer phenomena. Sarviya and Fuskele [4] reviewed the brief description of thermal conductivity of nanofluids. Prabhakar et al. [5] studied the effects of inclined Lorentz forces for zero normal fluxes of nanoparticles at the stretching sheet. A brief study on BL tangent hyperbolic flow of nanofluid is explained in these articles [6–17]. Fluids in which viscosity remain constant, no matter the shear applied on the fluid with constant temperature, e.g, water, air, gasoline, and alcohol On the other hand, non-Newtonian fluids are the fluids in which when stresses applied, viscosity will be changed, for example, quicksand, cornflour-water, silly putty, ketchup, cream, paint, glue, detergents, and molten plastics. Non-newtonian fluids are more complex to deal than the Newtonian fluids [18]. They are so difficult to deal with resulting equations and classical Newtonian models. The tangent hyperbolic model was the famous non-Newtonian model that was first introduced by Pop and Ingham [19] with the extraordinary uses in modern industries. Reddy et al. [20], Khan et al. [21], Mondal et al. [22], and Jabeen et al. [23] presented their numerical studies on MHD flows of tangent hyperbolic nanofluid over stretched sheet. Similarly, many applications of flow over a wedge-shaped surface are in scientific and manufacturing industries, for example, hydrodynamics, magnetohydrodynamics, and crude oil exploration and extraction.In 1937, Hartree [24] studied solutions on wedge. A brief study was devoted to analysis the differential parameters on the wedge flow. For example, the authors in [25–28] explained the heat transfer phenomena in wedge surface. Only few studies on unsteady MHD flow of tangent nanofluid along a wedge are in the literature, i.e., Kebede et al. [29] used a method to give the approximation of two-dimensional tangent hyperbolic Nanofluids over a stretched wedge in a porosity medium. We use the numerical method to approximate the solution of unsteady tangent hyperbolic flow over a wedge with suction-injection, and then we compared the results of numerical methods that we used to calculate the solution of nonlinear equations [30–32]. In this study, we have discussed numerical approximation of heat and mass transfer of a tangent hyperbolic nanofluid flow over a wedge shaped surface along with injection-suction. Taking into account rising demand of modern technology, every field of engineering and science uses nanoparticles to enhance their performance and minimize the cost of production. Fluids that have been traditionally used for heat transfer application in manufacturing industries have rather a low thermal conductivity. We used nanofluids which enhance the thermal conductivity behavior and rate of heat transfer as compared with the conventional fluids. Fluids, in general, may reveal both Newtonian and non-Newtonian behavior. Newtonian fluids behavior is observed when there is linear relationship in between the shear stress and the rate of deformation; otherwise, for any deviation for this behavior, fluids are treated as non-Newtonian fluids. It is so difficult to comprehend all the particular behaviors of non-Newtonian fluids as compared with the Newtonian fluids. Many of the fluids used in industries for the manufacturing of paint, honey, toothpaste, butter, shampoo, custard, soap solution, and so on do not obey Newton’s viscosity law. It is very difficult to analyze all the characteristics of such fluids. So, in many non-Newtonian formulations in the prior survey, tangent hyperbolic nanofluids are the famous non-Newtonian model that we considered here. A system of partially differential equations is converted to a system of ODEs under appropriate similarity variables. Numerical-based approaches are applied for numerical computations for the solution of nonlinear equations. It is expected that the convergence of these methods will give better results than the other methods.The numerical method always works in iterations and is easy to compute. When the analytical solution does not exist to solve a system, then numerical solution is preferred. These methods take short time to execute and cast less computational cost. MAPLE d solves basically a method to obtain the numerical solution for both initial and boundary value problems. This method is also an effective solver, with an approach towards the increasing Convergence criteria. Bvp4c is a fourth-order numerical method to solve boundary value ordinary differential equation in MATLAB software. The computed numerical results are demonstrated for a variety of nondimensional constraints in the equations. It is used for different parametric values, involving slip parameter, Prandtl number, Brownian movement variableNb and Nt, the thermophoresis variable, the Schmidt number Sc, and the ratio stretching parameter Sc. Also, M is the material fluids parameter, and n is the microrotation parametric influence. The MATLAB BVP solver bvp4c basically a typical finite difference code that deals with the three-stage LobattoIIIa formulational approach, which is used to obtain numerical results using bvp4c methodology.


## 2. Constitutive Relation for Tangent Hyperbolic Nanofluid

The non-Newtonian fluids have number of formulations in the literature that exhibit different rheological fluid behaviors. Amongst them, the tangent hyperbolic fluid is one of an interesting non-Newtonian four constant fluid formulation which basically narrates a behavior of shear thinning. It means that the viscosity declines as the rate of shear is enhanced while the apparent viscosity varies cautiously in between zero and infinite rate of shear.So, for tangent hyperbolic fluid, the constitutive relation is given as follows:(1)τ=μ∞+μo+μ∞tanhΓω¯nω¯.In the above expression, ‘n’ indicates the power law constant, γ represents the time constant, τ denotes extra stress tensor whereas zero and infinite shear rate of viscosity is referred here as μo and μ∞, respectively. Also, ω¯ is defined as follows:(2)ω¯=12∑i∑jω¯ijω¯ji=12∏.In this expression,∏ is defined as follows:(3)∏=12trgradV+gradVT2.Since we are considering the behavior of shear thinning in case of tangent hyperbolic fluid, therefore, for the case ofΓω¯<1, then the extra stress tensor which is denoted by τ is typically reduced to the following:(4)τ=μoΓω¯nω¯=μo1+Γω¯−1nω¯=μo1+nΓω¯−1ω¯.


## 3. Mathematical Formulation

Consider a tangent hyperbolic fluid flow occurring in a two-dimensional laminar way of mode across a wedge with injection-suction. The origin of Cartesian coordinate systems is established at vertex of the wedge, by using the Cartesian coordinate system which is (x, y). The x-axis is lying along the surface of the wedge whereas the y-axis is lying perpendicular to the surface of wedge. Assume that the external flow moves with a free stream velocity of Uex,t=bxm/1−ct and the fluid enters the wedge boundary layer with a wall velocity of Uwx,t=axm/1−ct subject to the magnetic domain of B¯=0,Bo acting perpendicularly to the surface of wedge, where t representing the time variable while a, b, m, and c depict constants assuming a>0 and a<0 that demonstrates the rates of stretching and shirking of wedge, respectively. Here, m is Falkner–Skan constants defined as m=β/2−β and also described as a total angle of wedge Ω by β=Ω/π. Suppose the temperature of the surface Tw and concentration CW of the wedge vary in the form of power-law as Tw=T∞+bxm/1−ct and Cw=C∞+bxm/1−ct as shown in Figure 1.Figure 1
Physical demonstration of the flow model.The equations, provided as follows, express the conservation rules regulating flow phenomena when all of the following assumptions are taken into account [29]:(5)∂u∂x+∂v∂y=0,(6)∂u∂t+u∂u∂x+v∂u∂y=∂Ue∂t+Ue∂Ue∂x+ν1−n+2nΓ∂u∂y∂2u∂y2−σBo2ρf+μK1u−Ue+gβTT−T∞−βcC−C∞sinΩ2,(7)∂T∂t+u∂T∂x+v∂T∂y=α∂2T∂y2+τDB∂C∂y∂T∂y+DTT∞∂T∂y2+16σ∗T∞33k∗ρCpf∂2T∂y2+σBo2ρCpu−Ue2+νCp1−n∂u∂y2+nΓ2∂u∂y3+QoρCpfT−T∞,(8)∂C∂t+u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2.The components of velocity are (u, v) along x- and y-axis. C and T are the nanofluid concentration and temperature in the BL region. μ and ρ are the dynamic viscosity and density of the nanofluids. ν=μ/ρ is the kinematic viscosity, n is the power law index. The numbers Γ,g,K1,βT,andβC signify material constant with time dependency, the gravitational acceleration magnitude, the porous medium permeability, volumetric thermal, and the expansion coefficient of concentration, respectively. Mean absorption and Stefan–Boltzmann constants are k∗ and σ∗, respectively; and the coefficient Qo represents the generation of heat when Qo>0 and the absorption of heat when Qo<0.We apply some boundary conditions on the above four PDEs. The boundary conditions near and far from the wedge’s surface are used as follows [29]:(9)u=Uwx,t,v=Vwt,−kf∂T∂y=hfTw−T,−DB∂C∂y=hcCw−Caty=0,u⟶Uex,t=bxm1−ct,T⟶T∞,C⟶C∞aty⟶∞,where Vw=Vo/1−ct represents the transpiration velocity also representing the transmission of mass at the stretching wedge surface with Vo as the constant velocity value, and convective heat and mass transfers coefficients are represented by hf and hc, respectively.Te transformation functions are as follows [29]:(10)η=y1+mUe2ν1−ct,ψx,y,t=2νxUe1+mfη,θη=T−T∞Tw−T∞,ϕη=C−C∞Cw−C∞.Here,η denote the nondimensional quantity for transformation purpose and ψ is the stream function with following factors u=∂ψ/∂y and v=−∂ψ/∂x. The dimensionless velocity and thermal and concentration profiles are represented by fηandθηandϕη, respectively.By the similarity transformation into the governing equations where the continuity equation as (5) and the equations as (6) to (8) came up in the nondimensional system of ordinary differential equations as follows:(11)1−n+nWe1+mf″f′″+ff″+2mm+11−f′2+Am+12−2f′−nf″+2mm+1M+Da1−f′+Grθ+GcϕsinΩ2=0,1Pr1+43Rdθ″+Nbθ′ϕ′+Ntθ′2+fθ′−2mm+1f′θ−Am+1ηθ′+2θ+2mm+1MEcf′2−2f′+1+Ec1−nf″2+nWef″3+2mm+1Q=0,ϕ″+PrLefθ′−2mm+1f′ϕ−Am+1ηϕ′+2ϕ+NtNbθ″=0,where prime′ indicates the differentiation with respect to η. And we have some dimensionless parameters which are as follows:(12)Ec=uw2cpfTw−T∞,Da=νxK1Ue,Gr=gβTTw−T∞xU∞2,Gc=gβCCw−C∞xU∞2,where Gr and Gc are temperature and mass buoyancy parametric effects, respectively.(13)Nb=τDBCw−C∞ν,Nt=τDTTw−T∞νT∞,where Nb and Nt exhibit the Brownian movement and the thermophoresis parametric influences, respectively.(14)Q=xQoρCpfUw,where Q represents the heat generation effect when Q>0 and heat absorption effect when Q<0.(15)We=ΓUe3/2νx,A=cbxm−1,M=σBo2xρUe,Pr=να,Rd=4σ∗T∞3kfk∗,Le=αDB.With the help of abovementioned similarity transformation, boundary conditions are as follows:(16)f0=S,f′0=δ,θ′0+Bi11−θ0=0,ϕ′0+Bi21−ϕ0=0atη=0,f′∞⟶1,θ∞⟶0,ϕ∞⟶0atη⟶∞,where S=Vo/aν represents the suction when S<0 and injection when S>0 and δ=Ue/Uw represents the velocity ratio parameter.Biot numbers are as follows:(17)Bi1=hfkf2νUem+1,Bi2=hsDa2νUem+1.The wall frictionCf, the local Nusselt number Nux, and Sherwood number Shx are given as follows:(18)Cfx=τwρfUe2,Nux=xqwαTw−T∞,Shx=xqmDBCw−C∞,where(19)τw=μ1−n∂u∂y+nΓ2∂u∂y2y=0,qw=−α+16σ∗T∞33ρCpfk∗∂T∂yy=0,qm=−DB∂C∂yy=0.We get the following relationships by substituting and applying the similarity transformation:(20)ReCf=1+m1−n+1+m2nWef″0f″0,NuxRe=−1+m1+43Rdθ′0,ShxRe=−1+mϕ′0,where Re=Uex/ν is the local Reynold’s number.


## 4. Results and Discussion

Under the boundary layer approximations, the dimensional form of a system of PDEs is yielded against the numerous distinct profiles which are transformed into ODEs by utilizing suitable similarity transformation. To explore the numerical results of velocity, temperature, and concentration profile, a numerical solution with the help of the bvp4c algorithm using MATLAB and MAPLE d-solve command has been constructed. Table1 shows that present results are in good agreement with previous studies reported in [11, 25, 29]. In Tables 2–4, the effects of dimensionless parameters upon local Sherwood #Shx, Nusselt #Nux, and skin friction Cfx are summarized.Table 1
Calculation ofRex1/2Cfx for different values of m.
Pop2013Ullah2017Tesfaye2020Present resultsf″0f″0f″0θ′0ϕ′0f″0θ′0ϕ′0m=0.00.46960.46960.46880.42280.23260.46970.42050.2323m=1/23−−0.56900.56930.45980.23480.56900.45420.2391m=1/110.65500.65500.65540.48180.24530.65500.48470.2457m=1/70.73200.73220.50900.25160.73190.51340.2523m=1/50.80210.80210.80230.54020.26690.80210.54030.2588m=1/30.92770.92770.92770.59040.26990.92770.59040.2716m=1/21.0389—1.03890.63690.28221.03890.63660.2842m=1.001.23261.23261.23260.72020.30631.23260.72040.3086m=5.001.55041.55030.86480.53171.55030.86480.3545m=100.01.67941.67940.92520.37261.67940.92520.3748Table 2
Calculation ofRex1/2Cfx for different auxiliary parameters.
Rex1/2CfxMATLABMAPLES=0.2S=−0.2S=0.2S=−0.2M0.01.12130.91801.12130.91801.51.36811.17261.36811.17263.51.63661.44581.63661.44585.51.86411.67581.86981.6757We0.21.14760.93721.14750.93580.61.14870.94531.14850.94501.21.16130.95721.16110.95711.81.17460.96841.17450.9682A0.01.10540.90061.10540.90060.21.17290.97291.17290.97290.41.23931.04381.23931.04380.61.30421.11281.30421.1128Da0.01.12130.91801.12130.91800.51.20951.00961.20951.00961.01.29141.09401.29141.09401.51.36801.17251.36801.1725δ0.11.24611.01741.24611.01740.21.14120.93811.14120.93810.31.02860.85131.02860.85130.40.90870.75710.90870.7571Gr0.11.10130.89551.10130.89550.41.12980.92661.12980.92660.81.16730.96721.16730.96721.21.20431.00711.20431.0071Gc0.11.12790.92201.12790.92200.41.13960.93711.13960.93710.81.15510.95701.15510.95701.21.17050.97681.17050.9768Table 3
Calculation ofRex−1/2Nux for different auxiliary parameters.
Rex−1/2ShxMATLABMAPLES=0.2S=−0.2S=0.2S=−0.2M0.00.40140.35060.40140.35061.50.40440.35600.40440.35603.50.40670.36000.40670.36005.50.40800.36230.40800.3623A0.00.38170.32370.38170.32370.20.41950.37490.41950.37490.40.45030.41450.45030.41450.60.47550.44580.47550.4458Q0.10.46060.41180.46060.41180.20.43110.38140.43110.38140.30.40160.35100.40160.35100.40.37220.32070.37220.3207Nb0.10.40160.35100.40160.35100.20.39960.34860.39960.34860.30.39760.34610.39760.34610.40.39570.34370.39570.3437Nt0.10.40160.35100.40160.35100.20.39960.34890.39960.34890.30.39760.34860.39760.34860.40.39560.34470.39560.3447Ec0.010.40490.35400.40490.35400.20.37390.32590.37390.32590.50.32460.28110.32460.28110.80.27460.23570.27460.2357Pr1.40.41700.35740.41700.35741.60.42280.35930.42280.35931.80.42770.36060.42770.36062.00.43200.36150.43200.3615Rd1.01.05750.94871.05750.94871.41.02830.93161.02830.93161.81.00430.91711.00430.91712.20.98430.90480.98430.9048Table 4
Calculation ofRex−1/2Shx for different auxiliary parameters.
MRex−1/2ShxMATLABMAPLES=0.2S=−0.2S=0.2S=−0.20.00.40520.36530.40520.36531.50.40790.37010.40790.37013.50.41040.37450.41040.37455.50.41240.37770.41240.3777A0.00.40200.35860.40200.35860.20.40880.37250.40880.37250.40.41570.38520.41570.38520.60.42230.39650.42230.3965Nb0.10.40540.36570.40540.36570.20.41430.36870.41430.36870.30.41720.36970.41720.36970.40.41870.37020.41870.3702Nt0.10.40540.36570.40540.36570.20.38930.36160.38930.36160.30.37420.35860.37420.35860.40.36010.35670.36010.3567Sc0.50.29160.28280.29160.28281.00.35140.32780.35140.32781.50.38390.35100.38390.35102.00.40540.36570.40540.3657A graphical demonstration for the behavior of emerging pertinent parameters associated different profiles is disclosed. Figure2 demonstrates the effects of M on velocity distribution. It can be seen that velocity profile increases by enhancing the values of M in the presence of suction which is physically fit with the argument that the decelerating force shows lesser domination in the entire flow region while an opposite trend followed up by temperature and concentration profile can be seen in Figures 3 and 4. We is the division of relaxation time and viscous forces of fluid. The effects of various profiles, i.e., f′, θ, and ϕ are exhibited in Figures 5–7. It is observed that the hiking (We) increases and the velocity of fluid starts to decline. On contrary, the temperature and concentration got increased which is due to the fact that the augmented Weissenberg number causes an increase in the relaxation time of the fluid flow which results in thicker nanofluid and enhanced resistance in the flow.Figure 2
Velocity profile versusM.Figure 3
Temperature profile versusM.Figure 4
Concentration profile versusM.Figure 5
Velocity profile versus We.Figure 6
Temperature profile versus We.Figure 7
Concentration profile versus We.Figures8–10 demonstrate that an increment in the (A) outcome is the hike in the velocity of fluid; however, it ultimately decreases the temperature and concentration profiles.Figure 8
Velocity profile versusA.Figure 9
Temperature profile versusA.Figure 10
Concentration profile versusA.The permeability of the medium is analyzed in the form ofDa whose effects are crucially important on the fluid’s velocity distribution as depicted in Figure 11. It can be measured that the increasing values of the Darcy number cause a hike in the velocity of nanoparticles. It is evident in the practical phenomenon as the increased values in this parameter yield hike in porous medium that turns out as the minimization in the hurdles based along the entire flow path. This permits the free flow of the fluid with increased fluid velocity in the BL region.Figure 11
Velocity profile versus Da.The effects of the non-Newtonian fluid are significantly important in the region of BL. These impacts are intimated in the form of (n), and results on the velocity profiles are exhibited in Figure 12. The velocity profile depicts an increasing trend with the gradual increment in distinct values of parameter. The reason behind this, as the parametric value start increasing, the fluid’s nature sharply changes, i.e., transmission between shear thinning and shear thickening state. The results given in Figures 13 and 14 show the influences of thermal and concentration distributions with the amendments in Nb parameter. The collisions between the nanoparticles increased due to which they start to move in an irregular manner that trivially rises the temperature and concentration profile.Figure 12
Velocity profile versusn.Figure 13
Temperature profile versusNb.Figure 14
Concentration profile versusNb.Figures15 and 16 show the directly proportional behavior of temperature and concentration profile in response to Nt (thermophoresis parameter). It can be observed that with the increasing behavior of Nt, the coming nanoparticles transfer from the hot surface towards cold ambient fluid state due to which the temperature rises in the boundary layer of fluid. Hence, the thickness of thermal boundary layer is improved. Also, the nanoparticles concentration is increased in the fluid, so the concentration boundary layer is enhanced significantly.Figure 15
Temperature profile versusNt.Figure 16
Concentration profile versusNt.Figure17 highlights the behavior of the Prandtl number (Pr) on thermal profile. It can be seen that with higher values of Pr, the temperature profile ultimately declines. Thus, it has an inverse relation with thermal diffusivity. So, increase in Pr delineates the thermal graph.Figure 17
Temperature profile versus Pr.Effects of behavior of different profiles against radiation parameter are demonstrated in Figure18. It is indicated that the increase in R inspires significant hike in the thermal BL region. Physical agreement can be seen is in the existence of thermal radiation relates to high kinetic energy of the molecules in the entire flow field.Figure 18
Temperature profile versusR.In Figure19, effects of the Ec can be observed on temperature profile. It is depicted that the increment in dissipation parameter produces frictional heating force which is to be stored in the fluid layer. Such frictional heating enhanced the thickness of temperature distribution at the surface of the wedge.Figure 19
Temperature profile versus Ec.Figure20 is plotted to emanate the effects on concentration profile for distinct parametric values of the Schmidt number (Sc). The different values for Sc are chosen to be 0.5, 1.0, 1.5, and 2.0 to indicate the diffusing chemical species in gasses like H2,H2O,andNH3. It is to be noted that the graph of incremented (Sc) tends to decrease fluid concentration.Figure 20
Concentration profile versus Sc.Figure21 illustrates that temperature profile is increasing functions of the heat source parameter.Figure 21
Temperature profile versusQ.The behavior of parameterδ on velocity profile is demonstrated in Figure 22. It is of highly consideration that velocity profile got enhanced by increasing the parameter δ.Figure 22
Velocity profile versusδ.The results portrayed in Figures23 and 24 represent the effects of thermal and concentration buoyancy parameters on velocity profile. It can be observed that the buoyancy parameters make substantial variation in velocity profile.Figure 23
Velocity profile versus Gr.Figure 24
Velocity profile versus Gc.The influences on the numerous distinct profiles are determined for of wedge angle in Figures25–27. Figures 25–27 depict that the velocity of fluid is increased while opposite trend is carried out by the other two profiles. Physically this corresponds to the fact that the increment in wedge angle parameter causes increase in the pressure of fluid. Moreover, changes in the hiking changing of wedge angle parameter influence the temperature near the surface of the wedge.Figure 25
Velocity profile versusm.Figure 26
Temperature profile versusm.Figure 27
Concentration profile versusm.


## 5. Conclusion

Unsteady, two-dimensional tangent hyperbolic fluid of nanofluid with injection-suction is studied. Following results were investigated:(i)
With the increase inM, the velocity profile shows behavior of increasing velocity with injection suction. And temperature and concentration profiles show the behavior of decreasing temperature and concentration with injection-suction.(ii)
By increasing the Weissenberg number parameter (We), the velocity profile decreases and temperature and concentration profile will increase with injection-suction. The Weissenberg number is the significance in the study of flow of non-Newtonian fluids, depending upon the reheology and ratio of elastic to viscous forces.(iii)
With the increase in the unsteadiness parameter (A), the velocity profile shows that velocity increases and temperature concentration profiles show that both temperature and concentration decrease with injection-suction.(iv)
With the increase in the Darcy number (Da), velocity profile shows that velocity increases with injection-suction.(v)
With the increase in the power law index parameter (n), velocity profile increases with injection-suction.(vi)
By increasing the Brownian motion parameterNb, temperature profile shows that temperature increases and concentration profile shows that concentration decreases with injection-suction.(vii)
By increasing the thermophoresis parameterNt, temperature and concentration profiles show that both are increased with injection-suction.(viii)
By increasing radiation parameter (R), temperature profile shows that temperature decreases with injection-suction.(ix)
By increasing the Eckert number parameter (Ec), the temperature profile shows that temperature will increase with injection-suction.(x)
By increasing the Schmidt number parameter, the concentration profile shows that concentration decreases with injection-suction.(xi)
By increasing the heat source (Q), the temperature profile shows that temperature increases with injection-suction.(xii)
By increasing the velocity ratio parameter (d), velocity profile will show that velocity increases with injection-suction.(xiii)
With the increment inm, the velocity profile shows that velocity increases with injection-suction while temperature and concentration profile show that both temperature and concentration decrease with injection-suction.



---

*Source: 8653091-2021-10-15.xml*