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

The analysis of unsteady tangent hyperbolic nanoﬂuid ﬂow past a wedge with injection-suction, because of its beneﬁcial uses, has gained a lot of attention. The present study is mainly concerned with tangent hyperbolic nanoﬂuid (non-Newtonian nanoﬂuid). First, we have converted the system of partial diﬀerential equations (PDEs) to a system of ordinary diﬀerential 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 technique bvp4c and investigated the velocity, temperature, and concentration proﬁles. The accuracy of the ﬂow 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 eﬀects of distinct dimensionless parameters on the various proﬁles 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 ﬂow and brings stability in the ﬂow. It was also noticed that by increasing the Darcy number, velocity proﬁle increases in both injection-suction cases.


Introduction
e 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][7][8][9][10][11][12][13][14][15][16][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. Nonnewtonian fluids are more complex to deal than the Newtonian fluids [18]. ey are so difficult to deal with resulting equations and classical Newtonian models. e 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][26][27][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 twodimensional 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][31][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. Numericalbased 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. e 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. ese 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.
is 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. e 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 variable N b and N t , 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.
e MATLAB BVP solver bvp4c basically a typical finite difference code that deals with the threestage LobattoIIIa formulational approach, which is used to obtain numerical results using bvp4c methodology.

Constitutive Relation for Tangent
Hyperbolic Nanofluid e 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) In the above expression, 'n' indicates the power law constant, c 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: In this expression, is defined as follows: 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:

Mathematical Formulation
Consider a tangent hyperbolic fluid flow occurring in a twodimensional laminar way of mode across a wedge with injection-suction. e origin of Cartesian coordinate systems is established at vertex of the wedge, by using the 2 Mathematical Problems in Engineering Cartesian coordinate system which is (x, y). e 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 U e (x, t) � (bx m /1 − ct) and the fluid enters the wedge boundary layer with a wall velocity of U w (x, t) � (ax m /(1 − ct)) subject to the magnetic domain of B � (0, B o ) 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 T w and concentration C W of the wedge vary in the form of power-law as as shown in Figure 1. e equations, provided as follows, express the conservation rules regulating flow phenomena when all of the following assumptions are taken into account [29]: zC zt e 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. e numbers Γ, g, K 1 , β 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 Q o represents the generation of heat when Q o > 0 and the absorption of heat when Q o < 0.
We apply some boundary conditions on the above four PDEs. e boundary conditions near and far from the wedge's surface are used as follows [29]: where represents the transpiration velocity also representing the transmission of mass at the stretching wedge surface with V o as the constant velocity value, and convective heat and mass transfers coefficients are represented by h f and h c , respectively.
Te transformation functions are as follows [29]: Here, η denote the nondimensional quantity for transformation purpose and ψ is the stream function with following factors u � (zψ/zy) and v � − (zψ/zx). e 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: where prime ′ indicates the differentiation with respect to η. And we have some dimensionless parameters which are as follows:   Mathematical Problems in Engineering where Gr and Gc are temperature and mass buoyancy parametric effects, respectively.
where N b and N t exhibit the Brownian movement and the thermophoresis parametric influences, respectively.
where Q represents the heat generation effect when Q > 0 and heat absorption effect when Q < 0.
With the help of abovementioned similarity transformation, boundary conditions are as follows: represents the suction when S < 0 and injection when S > 0 and δ � U e /U w represents the velocity ratio parameter.
Biot numbers are as follows: , e wall friction C f , the local Nusselt number Nu x , and Sherwood number Sh x are given as follows: where We get the following relationships by substituting and applying the similarity transformation: where Re � (U e x/]) is the local Reynold's number.

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. Table 1 shows that present results are in good agreement with previous studies reported in [11,25,29]. In Tables 2-4, the effects of dimensionless Mathematical Problems in Engineering parameters upon local Sherwood # Sh x , Nusselt # N ux , and skin friction C fx are summarized. A graphical demonstration for the behavior of emerging pertinent parameters associated different profiles is disclosed. Figure 2 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  Table 2: Calculation of (Re x ) 1/2 C fx for different auxiliary parameters. 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. e 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. Figures 8-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. e permeability of the medium is analyzed in the form of (D a ) 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.
is permits the free flow of the fluid with increased fluid velocity in the BL region. e effects of the non-Newtonian fluid are significantly important in the region of BL. ese impacts are intimated in the form of (n), and results on the velocity profiles are exhibited in Figure 12.
e velocity profile depicts an Table 3: Calculation of (Re x ) − 1/2 N ux for different auxiliary parameters.
( increasing trend with the gradual increment in distinct values of parameter. e 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. e results given in Figures 13 and 14 show the influences of thermal and concentration distributions with the amendments in N b parameter. e collisions between the nanoparticles increased due to which they start to move in an irregular manner that trivially rises the temperature and concentration profile. Figures 15 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 17 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. us, it has an inverse relation with thermal diffusivity. So, increase in Pr delineates the thermal graph.
Effects of behavior of different profiles against radiation parameter are demonstrated in Figure 18. 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.
In Figure 19, 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 20 is plotted to emanate the effects on concentration profile for distinct parametric values of the Schmidt number (Sc). e 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 H 2 , H 2 O, and NH 3 . It is to be noted that the graph of incremented (Sc) tends to decrease fluid concentration. Figure 21 illustrates that temperature profile is increasing functions of the heat source parameter. e 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 δ.
e results portrayed in Figures 23 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.

Conclusion
Unsteady, two-dimensional tangent hyperbolic fluid of nanofluid with injection-suction is studied. Following results were investigated: (i) With the increase in M, 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. e 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 parameter (N b ), temperature profile shows that temperature increases and concentration profile shows that concentration decreases with injection-suction. (vii) By increasing the thermophoresis parameter (N t ), 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 in m, 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.

Data Availability
e data set used to support the results and conclusion is included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.