Scrutinization of Waste Discharge Concentrations in Eyring-Powell Nanoﬂuid Past a Deformable Horizontal Plane Surface

: Nanomaterials have been the focus of intense study and growth in the modern era across the globe because of their outstanding qualities, which are brought about by their nanoscale size; for instance, increased adsorption and catalysis capabilities plus signiﬁcant reactivity. Multiple investigations have veriﬁed the fact that nanoparticles may successfully remove a variety of pollutants from water, and, as a result, they have been utilized in the treatment of both water and wastewater. Therefore, the current research intent is to examine the nonlinear heat source/sink inﬂuence on the 3D ﬂow of water-based silver nanoparticles incorporated in an Eyring–Powell ﬂuid across a deformable sheet with concentration pollutants. Silver particles have been used intensively to ﬁlter water, due to their potent antibacterial properties. The leading equations involving partial differential equations are renewed into the form of ordinary ordinary differential equations through utilizing the appropriate similarity technique. Then, these converted equations are solved by utilizing an efﬁcient solver bvp4c. Visual displays and extensive exploration of the different impacts of the non-dimensional parameters on the concentration, temperature, and velocity proﬁles are provided. Also, the important engineering variables including skin friction, the rate of heat, and mass transfer are examined. The ﬁndings suggest that the mass transfer rate declines due to pollutant parameters. Also, the results suggest that the friction factor is uplifted by about 15% and that the heat transfer rate, as well as the mass transfer rate, declines by about 21%, due to the presence of the nanoparticle volume fraction. We believe that these results may improve the ﬂow rate of nanoﬂuid systems, improve heat transfer, and reduce pollutant dispersal.


Introduction
Environmental harms caused by water pollution also contribute to the dangerous consequences of pollution in the air on human health.Water contamination also has a negative effect on the social and economic advancement of the afflicted cultures and nations.In the 21st century, it is increasingly difficult to find clean, fresh water; dirty water is dangerous to human life, according to the latest UN study [1,2].When undesired materials penetrate into water reservoirs or bodies and when water is polluted or contaminated, it cannot be used for other purposes or for consumption by humans.With the objective of solving this rising issue, numerous chemical, mechanical, and physical techniques exist.Additionally, researchers are still exploring many new procedures to improve inexpensive processes of water filtration [3,4].The recently created field of nanotechnology is currently offering an intriguing way of purifying water with affordable prices, significant operating efficacy in removing contaminants, and recycling ability [5].Cintolesi et al. [6] focused on the interface of thermal and inertial forces in urban canyons and how these forces affect turbulent properties and pollution removal.Chinyoka and Makinde [7] investigated the kinetics of polymeric pollutant dispersal of Newtonian fluid flow through a rectangular channel.The numerical model offers assessments for identifying contamination levels and offers perspectives on pollution hazards arising from an inappropriate release of hydrocarbon products.Recently, Vinutha et al. [8] inspected the flow of H 2 O-based ternary nanofluids from divergent/convergent channels by incorporating concentrated contaminants and a nonlinear heat source/sink.
Materials with structural elements between 1 and 100 nm in size are referred to as nanomaterials [9].Nanomaterials are radically distinct from normal materials in terms of their magnetic, optical, electrical, and mechanical properties, because of their nanoscale dimension.Many different types of nanomaterials exhibit catalysis, high reactivity, and adsorption.Over the past few decades, nanomaterials have undergone significant research and study and have successfully been applied in a number of fields, such as medicine [10], catalysis [11], biology [12], and sensing [13].Silver nanomaterials are employed in a variety of industries, including those of medicines, medical devices, food and beverage, cosmetics, and water purification, due to their antibacterial qualities.More specifically, the application of nanoparticles to the purification of water and wastewater has attracted a lot of interest.Due to being tiny, nanomaterials have significant surface vicinity and high-adsorption properties.Furthermore, nanoparticles are highly mobile in solutions [14].The exclusion of organic pollutants [15], heavy metals [16], inorganic anions [17], and microorganisms [18] by dissimilar nanomaterial kinds has been reported.Lu et al. [19] presented a complete review of the significance of nanomaterial for purifying as well as treating waste water.Yaqoob et al. [20] discussed the nanomaterial role in the treatment of water.Recently, Zahmatkesh et al. [21] discussed the capability of nanomaterials to eliminate viruses, fungus, and organic matter from wastewater.Numerous studies have shown that nanoparticles hold considerable promise for use in water and wastewater treatment.
The study of flow by incorporating non-Newtonian liquids has received a lot of interest recently.These liquids are used often in the production of coated sheets, plastic polymers, food, drilling muds, optical fibers, etc.These fluids have very complex connections between the flow field and shear stress, which presents researchers with fascinating difficulties.Despite these obstacles, the field's experts are even contributing significantly to the study of non-Newtonian fluids [22][23][24][25].The research into Powell-Eyring liquid across a constantly movable porous plane was scrutinized by Jalil et al. [26].They observed that the examination of the Powell-Eyring fluid model has a benefit, as it is developed from the kinetic theory of gases rather than from empirical relations.Roşca and Pop [27] looked into flow features and heat transfer through a shrinkable surface.To identify a stable solution, they carried out a stability analysis and observed that the first solution is stable and acceptable whilst the second solution is unstable and physically not acceptable.Rahimi et al. [28] utilized the collocation approach to achieve the outcome of the fluid flow of Eyring-Powell liquid across a stretchy sheet.The velocity is seen to decrease when the fluid substance parameter is increased, and to increase when the Eyring-Powell fluid substance parameter is expanded.The magneto influence on the radiative flow of Eyring-Powell fluid from a stretchable defaming surface with slip effects and chemical reaction was inspected by Reddy et al. [29].Larger radiative heat transfer values have been found to affect temperature more significantly.Larger radiative heat transport values have been found to affect temperature more significantly.Khan et al. [30] scrutinized the Joule heating influence on the MHD flow and heat transfer induced by Eyring-Powell fluid with irregular viscosity.Recently, Hassan et al. [31] examined the characteristics of fluid flow past an irregular thicker sheet with Eyring-Powell fluid, with an entropy effect.It was discovered that both temperature and velocity are affected differently by the time-scale factor.
The phenomena of an irregular heat sink or source have applications in both medical and engineering fields, such as the metallic sheet cooling, the design of thrust bearings, the revitalization of unrefined oil, etc. Tawade et al. [32] discussed the motion of unsteady thin film and heat transfer from a stretchable sheet with magnetic and irregular heat sink/source effects.The performance of heat transport is impacted by irregular heat parameters.Das et al. [33] inspected the impact of heat source/sink on the fluid flow and heat transfer, incorporating a water-based nanofluid moving past a shrinkable sheet.Thumma et al. [34] inspected the magnetic impact on the convective motion of nanoliquid with an unpredictable heat sink/source caused by a stretching sheet.The outcome is obtained by utilizing the well-known Keller box numerical technique.Acharya et al. [35] examined the impact of a non-uniform heat source on the magneto flow across a timedependent stretchable cylinder.Ramadevi et al. [36] inspected the consequences of an inconsistent thermal sink/source on the flow of an electrically conductive fluid through a varied thickness sheet.Lately, Khan et al. [37] have examined the characteristics of a jet flow incorporated H 2 O/Al 2 O 3 nanofluid moving past a slippery stretchable surface in a porous medium with an erratic heat source/sink.It has been concluded that the temperature is enhanced due to the heat source and declined owing to the heat sink factor.
According to the literature, no study has looked at the effects of pollutant concentration and irregular heat source/sink on the movement of water-based silver nanoparticles induced by Eyring-Powell fluid across a stretching/shrinking sheet.The aforementioned consequences have been explored in the present research work.The distinguished numerical bvp4c solver is employed to construct and work out the leading equations for current flow scenarios.There is a visual impact on many flow characteristics with pertinent physical significance.Notably, the present research is the initial attempt to examine the flow of heat and mass transfer in a non-Newtonian Eyring-Powell fluid incorporating the discharge concentration.The modeling of the nanofluid flow led to the development of the following research questions:

•
What effect does the heat sink/source have on the temperature disparity within the fluid layers near the boundary and the temperature distribution on the behavior of nanofluid?

•
What are the effects of skin friction, Nusselt number, and Sherwood number on the behavior of Eyring-Powell fluid in the presence of silver nanoparticles?

•
How do the silver nanoparticles affect the reduction in the pollutant dispersal?

Problem Formulation
Consider a 3D flow of water-based silver nanoparticles with heat and mass transport phenomenon in a non-Newtonian (Eyring-Powell) nanofluid past a deformable (stretching/shrinking) horizontal plane surface.The impacts of pollutant concentration and irregular heat source/sink are incorporated.Figure 1 shows the Eyring-Powell fluid flow structure when silver nanoparticles are included, where x q , y q and z q are Cartesian coordinates with velocities u q , v q , and w q , respectively.Here, the z q − axis is normal to the deformable sheet, and the x q − and y q − axes are placed along the streamwise and spanwise directions, respectively.Assume that the streamwise and spanwise velocities at the surface of the horizontal plane are described by u wq x q = Ax n q and v wq x q = Bx n q , where A, and B refer to the constant stretching rates in the x q − and y q − directions, respectively (see Weidman [38]).In addition, w q x q = w 0 communicates the mass transpiration velocity with w 0 > 0, and w 0 < 0 refer to the suction and injection, respectively, while the impermeable surface of the plane is symbolized by w 0 = 0. Additionally, T wq and T ∞ are thought to stand for the far-field temperature and the constant surface temperature, respectively.Furthermore, C wq and C ∞ illustrate, respectively, the consistent concentration at the wall surface and the outside (far-field) flow condition.The Eyring-Powell constitutive model has the following form and is in accordance with Narsimhan [39] and Ara [40].
The first expression describes the viscosity effect, and the second expression describes the elastic component.Taking into account the Maclaurin series from the second term, this has the form as follows: or By ignoring extraneous components and only making use of the first term and second term from the approximations, one can derive As a result, the fundamental model for Eyring-Powell in Equation ( 1) changes to the form: where β q and δ q indicate the Eyring-Powell fluid or material parameters and µ n f refers to the nanofluid viscosity.Considering these presuppositions, the leading governing equations of the nanofluid in terms of PDEs can be constructed as follows, with the assistance of the aforementioned premises and the boundary layer approximations [8,[26][27][28][29]: ρ n f u q ∂u q ∂x q + v q ∂u q ∂y q + w q ∂u q ∂z q = µ n f + 1 The physical boundary conditions (BCs) listed below are u q = γ b u wq x q , v q = γ b v wq x q , w q = w 0 , T q = T wq , C q = C wq at z q = 0, Here, β q and δ q indicate the Eyring-Powell fluid or material parameters, T q the nanofluid temperature, C q the nanofluid concentration, and D f the concentration molecu-lar diffusivity; Q q and b * correspond to the external pollutant strengths, A * b refers to the exponentially decaying space coefficient, and B * b signifies the time-dependent heat absorption/generation.In addition, the phenomenon of heat source and heat sink is created by the positive and negative consistent values of A * b and B * b .Also, γ b states the deformable (stretching and shrinking) horizontal plane surface with γ b < 0 for shrinking case and γ b > 0 for stretching case, whereas, γ b = 0 is the case of a static plane surface. ( q b wq q q b wq q q q wq q wq q Here, q β and q δ indicate the Eyring-Powell fluid or material parameters, q T the nanofluid temperature, q C the nanofluid concentration, and f D the concentration mo-  In addition, the rest of the mathematical symbols utilized in the aforesaid equations for the physical correlations of the silver nanofluid are nf ρ the density, nf k the thermal In addition, the rest of the mathematical symbols utilized in the aforesaid equations for the physical correlations of the silver nanofluid are ρ n f the density, k n f the thermal conductivity, and ρc p n f the specific heat capacity.The correlations are placed in Table 1.

Properties Correlations
Dynamic viscosity The notation ϕ s c suggests the solid-fraction volume of nanoparticles, whereas the particular case ϕ s c = 0 shrinks the above correlations into a regular water-based fluid.Also, the subscripts n f , s c , and f indicate the corresponding nanofluid, nanoparticles, and regular fluid.However, c p signifies the heat capacitance at constant or uniform pressure.Table 2 is constructed to show the data of the water liquid and the silver nanoparticles obtained from the experiment.The similarity variables are presented which obey the equation of the continuity Equation ( 6) as: where ξ is the pseudo-similarity variable and prime signifies the derivative with respect to ξ, n the power-law index.Also, C = C x q and C x q = dC(x q ) dx q .Using Equation ( 12), the momentum Equations ( 7) and ( 8) in the streamwise and spanwise directions are altered into the following form: The following concept was chosen as follows, in order to achieve the similarity equations (see Weidman [38]): Applying integration by parts, one now arrives at where K b and D b are the requisite arbitrary constants.Additionally, δ q = δ b x 3n−1 2 q and β q = β b x −( 3n−1 2 ) q are the Eyring-Powell fluid or material parameters with constants δ b and β b .Using Equations ( 15) and ( 16) and the above working variable expressions in the streamwise and spanwise momentum Equations ( 13) and ( 14) yield: Hence, we set the coefficients of FF and FG as equal to unity, so that and Equations ( 17) and ( 18) are simplified into the following form: 1 1 Here, the dimensionless influential parameters are Powell fluid parameters, and E b = B A the stretching-rate ratio parameter.Additionally, employing Equation ( 12), along with S(ξ) = the energy Equation ( 9) and concentration Equation ( 10) are transmuted into the following posited form: 1 Pr Using the strength of the variable external pollutant and Equations ( 15), ( 16) and ( 19), , the same procedure was used to update Equations ( 22) and ( 23), resulting in the following form 1 Additionally, the mass transpiration parameter for the surface of a horizontal plane is provided by and these are the new boundary conditions: where Pr = υ f /α f represents the Prandtl number and f wb the constant mass flux parameter, with f wb > 0 and f wb < 0 indicating the suction and injection, respectively, while f wb = 0, subject to the impermeable surface of the plane, In this investigation, it should be noted that if Σ b = Σ c = 0 (Newtonian fluid), f wb = 0 (impermeable surface of the plane), γ b = 1 (stretched surface of the plane), and ϕ s c = 0 (viscous fluid) are investigated, the above model reduces to Weidman [38] (see the third extension of the Bank's problem).Weidman [38] did not incorporate the heat transfer features because of the absence of the energy equation.Furthermore, Equations ( 20) and ( 27) can be simplified to align with work by Miklavčič and Wang [42] for linear shrinkable and Fang [43] for nonlinear shrinkable when viscous liquid, Newtonian fluid, and a porous flat plate are taken into account.Similar to previous papers, none of them examined the necessary energy equation or the concentration equation.
The local Nusselt number Nu x q , the Sherwood number Sh x q , and the coefficients of the skin friction in the streamwise C f and spanwise C g directions are the gradients of engineering physical significance.
Hence, Re x q x q = u wq (xq)xq , and Re y q x q = v wq (xq)xq υ f are the local Reynolds number.

Graphical Results and Discussion
The symmetrical heat source/sink and suction impact in a three-dimensional Eyring-Powell nanofluid induced by a deformable (expanding/contracting) horizontal plane surface with waste discharge concentrations was considered.The set of similarity equations involved various influential dimensionless parameters, namely, the nanoparticle volume fractions ϕ s c , the material or Eyring-Powell fluid parameters Σ b and Σ c , the Prandtl number Pr, the Schmidt number Sc b , the mass suction/injection parameter f wb , the deformable (expanding/contracting) parameter γ b , the local external pollutant source parameter δ a , and the local external pollutant source variation parameter δ b .The impression of these parameters on shear stress in the streamwise (SMW) and spanwise (SNW) directions, heat transfer rate (HTR), and mass transfer rate (MTR) are bounded via several graphs.In addition, this portion of the work exemplifies the explanation of the acquired graphs in terms of physical viewpoints as well as the validation of the considered approach for a limiting situation.For dual results, the prominent bvp4c was employed to work out the retrieved set of dimensionless similarity equations.The three-stage Lobatto formula is specifically executed by this built-in code, which is part of the MATLAB package and is based on a finite difference procedure (refer to Shampine et al. [44]).
The authentication/validation and accuracy of the considered study scheme are quantitatively shown in Table 3, for a limiting example.This computational table shows a comparison of the up-to-date HTR values with prior research works of Devi and Devi [45], Wang [46], and Khan and Pop [47], owing to the sundry values of Pr when f wb = 0, Σ b = 0.0, Σ c = 0.0, ϕ s c = 0, Sc b = 0.0, δ b = 0.0, δ a = 0.0, A * b = 0.0, B * b = 0.0, n = 1.0 and γ b = 1.0.From these tabular values, it is proved that the existing work is accurately aligned with prior research works.Further, this exceptional matching gives us excellent confidence that the located solutions for both (stable and unstable) branches are perfectly acceptable.Figures 2-7 are exclusively constructed to sightsee the influence of several controlling parameters on gradients for the stable branch solution (SBS) and unstable branch solution (UBS).In general, the SBS and UBS exist for the cases of suction and shrinking horizontal plane surface.The current problem raises dual (stable branch (SB) and unstable branch (UB)) solutions for the variance in the different distinct parameters, where the SB is indicated by the solid black lines and the UB solution is indicated by the dashed red lines.However, the point where the SB and UB merge is essentially referred to as a bifurcation point (BP).
Table 3.Comparison values of heat transfer rate for regular base fluid (ϕ s c = 0) with several values of Pr when Figure 2a-d display the consequence of the power-law index n on shear stress coefficients (SSCs) in both the SMW and SNW directions, and HTR and MTR versus f wb of the (water/Ag) nanofluid for the SBS and UBS, respectively.From Figure 2a,b, it is acknowledged that the shear stress in both the SMW and SNW directions increases for the SBS owing to the superior impacts of n, but a contradictory pattern is perceived near the points of bifurcation.However, the SSCs in both the SMW and SNW directions near the BP decline and then incline, with superior impacts of n for the UBS.In addition, the SSC in the SNW direction is higher for the UBS as compared to the outcome obtained in the direction of SMW because of the mounting values of the power-law index parameter.On the other hand, the HTR and MTR profiles were boosted up in both SBS and UBS, owing to the amplified values of n, see Figure 2c,d.Additionally, the SBS and UBS merged at a single position which is called the bifurcation or critical point (BP or CP).This BP is represented in all the graphs via a solid blue ball.The BP is not unique in these diagrams, because it is changed with the changing values of n.Meanwhile, the BP is mathematically represented by f wb C; where the outcome is unique, then it is expressed as f wb = f wb C. Also, the outcomes are not unique for the situation when f wb > f wb C, and no single and multiple outcome possibility arises with f wb < f wb C. The magnitude of the bifurcation values (BVs) escalates due to the superior impressions of n, and hence the flow separation of the boundary layer declines from the horizontal plane surface.The effects of the material or EPF parameter Σ c on SSCs in both (SMW and SNW) directions, and HTR and MTR of the (water/Ag) nanofluid for the SBS and UBS are depicted in Figure 3a-d, respectively.Advancing values of Σ c lead to a remarkable enhancement in the SSCs in both (SMW and SNW) directions and MTR for the UBS, but the HTR substantially reduces.However, the heat transfer rate upsurges for the SBS with the higher impacts of the non-Newtonian parameter Σ c , while the SSCs in both directions and the MTR fall for the SBS.In general, the material or non-Newtonian fluid parameter Σ c is inversely proportional to µ f of the fluid.The significant rise in the constraint Σ c generates a reduction in the dynamic viscosity of the fluid, and, as a result, the velocity enhances.Due to this circumstance, a higher velocity of the fluid strongly recommends that the friction factor coefficient in both paths for the Ag nanoparticles should decline.Moreover, the absolute properties of the BVs are uplifted with larger values of Σ c .This behavior graphically illustrates that the flow separation of the boundary layer from the surface of the plane decays successively, due to the larger consequences of Σ c .

Pr
Figure 4a-d show the effect of the solid nanoparticle volume fractions ϕ s c on the SSCs in both (SMW and SNW) directions, and HTR and MTR for the SBS and UBS, respectively.It is professed that for a fixed ϕ s c , two distinct (SB and UB) outcomes of each SSCs in both (SMW and SNW) directions, HTR and MTR, originated.It is seen in Figure 4a,b that the SSCs in both (SMW and SNW) directions are lower for the case of shrinking and higher for the case of stretching, due to the superior impacts of ϕ s c .Meanwhile, they both perform distinctly for the UBS with higher effects of ϕ s c .In addition, the SSC in the SNW direction is slightly better and stronger as compared to the value obtained in the case of the SMW direction.Also, the dual SBS and UBS exist in the case of shrinking, while for the case of stretching a single solution has been calculated.On the other hand, the MTR inclines for the UBS and declines for the SBS, due to the higher impacts of ϕ s c while the heat transfer rate was remarkably boosted up in both solution branches (see Figure 4c,d).Physically, the TCN of the non-Newtonian fluid enriches with superior consequences of ϕ s c , and thus the rate of heat transfer enlarges.Furthermore, the outcome of two distinct (SB and UB) branches is met at a single point called bifurcation.At this position, the fallouts are exceptional b > 0 (the heat-source parameter) leads to a noteworthy reduction in both the SBS and UBS.However, a contradictory behavior is seen in both outcomes when we upsurge the values of A * b , B * b < 0 (heat-sink parameter).Physically, this propensity is due to the incremental choices of A * a < 0, B * a < 0 which generates the heat in the familiar form of energy, and thus the rate of heat transfer rises.However, the heat produced is much less or negligible because of the augmentation in the sequential values of the heat generation parameter A * a > 0, B * a > 0, and as a response, the rate of heat transfer decelerates.In addition, the critical values are the same if we change the value of A * b and B * b to simultaneously positive and negative, see Figure 5a,b.Hence, the gap between the LB and UB curves is better for the circumstance of the heat source as compared to the heat sink parameter.
Figure 6a,b describe the influence of the external pollutant parameters δ a , and δ b on the MTR of the (water/Ag) nanofluid against the expanding/contracting parameter γ b for the SBS and UBS, respectively.From the graphs, it is seen that dual solutions are only possible in the case of contracting, while only a single solution is noted for the circumstance of expanding.In addition, the MTR declined for both solution branches due to the continuous improvement in the external pollutant parameters δ a , and δ b .Generally, by increasing the external pollutant parameters, a more significant influx of pollutants into the deformable horizontal plane surface increases pollutant concentrations.The fluid flow patterns may be affected, as the higher pollutant concentration alters the velocity distribution and may cause flow instabilities and increase temperature and concentration but reduce the rates of HT and MT.More significant improvement in the pollutant parameter will improve the pollutant load, which challenges the efficiency of pollution control techniques.Also, the critical point (CP) will remain the same if we change the values of δ a and δ b .Meanwhile, the gap between the solution curves for the MTR is better for the variations in the value of δ b as compared to δ a .
The consequences of the external pollutant parameters δ a and δ b on the MTR of the (water/Ag) nanofluid against the suction parameter f wb for the SBS and UBS are presented in Figure 7a,b, respectively.Here, the complete dual (SB and UB) solutions are found against the suction parameter for varying the influence of the parameters δ a and δ b .It is seen that the MTR declined continuously in both SBS and UBS with the sophisticated external pollutant parameters δ a and δ b .By intensifying the values of external pollutant parameters δ a and δ b , the bifurcation values | f wb C| will remain fixed and constant.Moreover, the mass transfer rate solution gap between the two curves (SB and UB) for the varying values of δ b is higher than the outcomes obtained for the changing values of δ a .Hence, these parameters can significantly impact the system's behavior and the effectiveness of waste treatment and pollution control.

Conclusions
In this manuscript, the separations of the boundary layer and multiple (SB and UB) solutions were explored for a 3D Eyring-Powell nanofluid flow past a porous deformable (expanding/contracting) horizontal plane surface.The impacts of erratic heat source/sink and pollutant concentration for the heat and mass transfer analysis were also taken into account.The new set of variable ansatz was used to transform the problem from PDEs to ODEs.The substantial outcomes of this model are listed below.

•
The dual (SB and UB) outcomes occur only for the limiting cases like shrinking surface and mass suction parameters, with variation of the other influential parameters.

•
The boundary layer separations lessened with the higher influence of n and Σ c but were augmented with superior impacts of ϕ s c .

•
The inclusion of the silver nanoparticles ϕ s c improves the heat transfer rate for both (SB and UB) consequences, whereas the rate of mass transfer declines for the SBS and rises for the UBS.

•
The mass transfer rate declined for both (SB and UB) solutions with the superior impact of the local external pollutant concentration parameters.

•
The shear stress coefficients in both (SMW and SNW) directions initially incline and then decline for the SBS owing to the higher impacts of ϕ s c , while the contrary behaviors are found for the UBS.

•
The mass transfer rate values increases with the significant impact of pollutant concentration parameters δ a , and δ b .
Further, this effort can be expanded by investigating an unsteady flow or buoyancy flow, due to their significant applications like nuclear reactors, flows in the atmosphere, and the ocean due to seasonal change, and varying densities along the vertical direction of a lake causing circulation.Heat transfer rate Sh x q Mass transfer rate C f , C g Skin friction coefficients in streamwise and spanwise direction Re x q x q , Re y q x q Local Reynolds numbers D f Concentration molecular diffusivity (m 2 /s) n Power-law index k Thermal conductivity (W/m K) Q q , b * External pollutant strengths w 0 Mass transpiration(m/s) velocity u q , v q , w q Velocity components along x q −, y q and z q − axes (m/s) x q , y q , z q Coordinates (m) Greek symbols β q , δ q Eyring-Powell fluid characteristics lecular diffusivity; q Q and * b correspond to the external pollutant strengths, * b A refers to the exponentially decaying space coefficient, and * b B signifies the time-dependent heat absorption/generation.In addition, the phenomenon of heat source and heat sink is created by the positive and negative consistent values of * b A and * b B .Also, b γ states the deformable (stretching and shrinking) horizontal plane surface with 0 b γ < for shrinking case and 0 b γ > for stretching case, whereas, 0 b γ = is the case of a static plane surface.

Figure 1 .
Figure 1.Physical model of the problem.

Figure 1 .
Figure 1.Physical model of the problem.

δ
a = b * C wq − C ∞ corresponds to the local pollutant external source parameter, andδ b = Q b A(C wq −C ∞ )denotes the local external pollution source variation parameter.

Figure 2 .Figure 3 .
Figure 2. (a) Shear stress in the streamwise direction, (b) shear stress in the spanwise direction, (c) heat transfer rate, and (d) mass transfer rate versus f wb for various values of n.

Figure 3 .
Figure 3. (a) Shear stress in the streamwise direction, (b) shear stress in the spanwise direction, (c) heat transfer rate, and (d) mass transfer rate versus f wb for various values of Σ c .

Figure 4 .Figure 5 .
Figure 4. (a) Shear stress in the streamwise direction, (b) shear stress in the spanwise direction, (c) heat transfer rate, and (d) mass transfer rate versus γ b for various values of ϕ s c .

Figure 5 .
(a) Heat transfer rate versus γ b for the various values of the internal heat source parameter and (b) for various values of the internal heat sink parameter.

Figure 6 .
Figure 6.(a) Mass transfer rate versus b γ for the different values of a δ and (b) for various values

Figure 6 .Figure 7 .
Figure 6.(a) Mass transfer rate versus γ b for the different values of δ a and (b) for various values of δ b .

Figure
Figure 2a-d display the consequence of the power-law index n on shear stress co- efficients (SSCs) in both the SMW and SNW directions, and HTR and MTR versus wb f of the (water/Ag) nanofluid for the SBS and UBS, respectively.From Figure 2a,b, it is

Figure 7 .
Figure 7. (a) Mass transfer rate versus f wb for the different values of δ a and (b) for various values of δ b .
(or unique) and scientifically expressed as γ b = γ b C, where γ b C corresponds to the BP or CP.For the specific domain areas such as γ b C < γ b < ∞, and −∞ < γ b < γ b C, the dual outcomes and no outcomes existed, respectively.Intensifying values of ϕ s c decreases |γ b C|, that is, it raises the flow separation of the boundary layer from the HPS.The impression of the internal heat source parameter A * b , B * b > 0, and the internal heat sink parameter, A * b , B * b < 0, on the rate of heat transfer of the H 2 O-based Ag particles for the dual (SB and UB) outcomes are depicted in Figure 5a,b, respectively.Growing values of A * b , B * heat absorption/generation c p Specific heat (J/Kg.K) E b Stretching-rate ratio parameter K b , D b Arbitrary constants u wq x q Velocity of the deformable surface (m/s) C wq Wall Concentration T wqTemperature of wall (K)T ∞ Ambient temperature (K) F(ξ), G(ξ)Dimensionless velocities in streamwise and spanwise direction S Specific heat capacity (J kg −1 K −1 ) ρ Density (kg m −3 ) ψ Stream function Σ b , Σ cEyring-Powell fluid

Table 2 .
[41]thermophysical properties of the base fluid and the silver nanoparticles[41].