Impact of Nonlinear Thermal Radiation and Entropy Optimization Coatings with Hybrid Nanoliquid Flow Past a Curved Stretched Surface

The present investigation addresses the flow of hybrid (nickel–zinc ferrite and ethylene glycol) nanoliquid with entropy optimization and nonlinear thermal radiation coatings past a curved stretching surface. Analysis was carried out in the presence of magnetohydrodynamic, heat generation/absorption, and convective heat and mass flux conditions. Solution of the modeled problem was attained numerically using MATLAB built-in function bvp4c. Impacts of prominent parameters on betrothed distributions were depicted through graphs and were well supported by requisite discussions. Numerically calculated values of Sherwood number were established in a tabulated form and were scrutinized critically. An excellent concurrence was achieved when results of the presented model were compared with previously published result; hence, dependable results are being presented. It was observed that concentration field diminished with increasing values of curvature parameter, though the opposite trend was noticed for velocity and temperature distributions. Further, it was detected that Nusselt number decreased with augmented values of radiation and curvature parameters.


Introduction
Numerous applications of heat transfer liquids or coolants can be found in a variety of fields, such as automobiles, industry, electronics, and cooling processes.In all such industrial applications, cooling by liquids has been used for years.The process of cooling by fluids may be the single phase (where there is no phase change in the coolant) or two-phase (where coolant liquid will experience a phase change).In the latter, latent heat influences the cooling efficiency [1].Several coolants, such as water, ethylene glycol, blend of water and glycol, propylene glycol and amalgamation of water, and propylene glycol, are used as coolants in automobiles and industrial cooling processes.In the last two decades, several researchers have devoted their efforts to increasing thermal conductivity of coolants, thereby improving heat transfer capabilities.The pioneering work of Choi et al. [2] introduced nanofluids by insertion of solid nanoparticles into liquids, thus enhancing the thermal properties of these liquids.
This pioneering work has remarkably revolutionized modern engineering and the industrial world.Nanofluids are an amalgamation of suspended solid material particles and customary liquids (ethylene glycol, water).This new type of advanced material possesses amazing capabilities that trigger the process of heat transfer and augments the thermal conductivity of the base fluid.Enhancement in the thermal conductivity and heat transfer is visualized once ferrite nanoparticles are added into the base liquid.Several examples featuring heat transfer can be quoted, including chemicals, cooling and heating system of buildings, and avionics cooling systems.Nanofluids exhibit potentially exceptional features in comparison to macrometer-sized particles.This is because nanoparticles have sufficiently larger surface area compared to micrometer-sized particles; this is the reason nanofluids possess incomparable capabilities of heat transfer [3].
In several electromagnetic applications with high permeability, e.g., electromagnetic wave absorbers and inductors, usage of nickel-zinc ferrite can be noticed.To minimize energy losses related to bulk powders, usage of nickel-zinc nanoparticles has been recommended by a number of researchers [4][5][6].In addition, a majority of electronic gadgets require such materials to be compressed into outsized shapes with the required thickness, which is reasonably challenging if the size of these particles is large enough.Several methods have been proposed to get nickel-zinc ferrite, including ball milling, precipitation, and hydrothermal.Ferrofluids are colloidal fluids comprising ferromagnetic or ferrimagnetic nanoparticles suspended in an electrically insulated hauler fluid.In the current examination, ethylene glycol (C 2 H 6 O 2 ) was taken as a carrier fluid.The assumed ferrite nanoparticle was nickel-zinc ferrite (NiZnFe 2 O 4 ) crystallize in the normal spinal structure.Typically, at room temperature, the inverted spinals are ferromagnetic and normal spinals are paramagnetic.Moreover, zinc ferrites act like antiferromagnetic in nature at low temperature.This feature makes ferromagnetic nanofluids more relevant in different real-world applications [7,8].The ferrofluid's flow with the effect of thermal gradients and the magnetic field was discussed by Neuringer [9].Majeed et al. [10] demonstrated the heat transfer investigation in a ferromagnetic fluid flow.
The subject of fluid flow past stretched surfaces has diverse engineering and industrial applications, including paper production, glass blowing, crystal growing, hot rolling, manufacturing of rubber sheets, annealing of copper wires, etc.The coined work of Crane [11] discussing the flow past a linearly stretching surface urged fellow researchers to discover more avenues in this exciting and interesting subject.This was followed by the remarkable work of Gupta and Gupta [12] who pondered on the flow past a spongy surface.Then, Chakrabarti and Gupta [13] examined the hydromagnetic flow past a stretched surface.Andersson et al. [14] considered the flow of power-law fluid past a surface, which was linearly stretched under the influence of magnetic forces.The flow of an Oldroyd-B fluid with the impact of generation/absorption was deliberated by Hayat et al. [15].Muhammad et al. [16] discussed the effect of thermal stratification in the ferromagnetic fluid on a stretching sheet.Ramzan and Yousaf [17] demonstrated that the elastic viscous nanofluid finished a bi-directional stretching surface in view of Newtonian heating.Hussain et al. [18] utilized the exponentially stretching sheet to scrutinize the flow of Jeffrey nanofluid with radiation effects.Some recent explorations highlighting various fluid flows past stretched surfaces with coatings can be found in references [19][20][21][22].
In today's cutting-edge engineering technology, curved stretching has a broad relevance because of its different uses in industry, for example, in transportation and electronics.Sanni et al. [23] attained a numerical solution for the viscous fluid flow on a curved stretched channel.Sajid et al. [24] inspected the ferrofluid (Fe 3 O 4 ) flow on a curved sheet with effects of Joule heating and magnetic forces.Rosca and Pop [25] studied time-dependent flow along a spongy curved surface.Imtiaz et al. [26] introduced the effect of homogeneous/heterogeneous reactions in ferrofluid embedded in a stretching surface.Naveed et al. [27] calculated heat transfer and used the micropolar fluid to analyze the effects over a curved surface with thermal radiation.
A literature review has specified that copious studies are available relating to nanofluids with linear/nonlinear/exponential stretching surfaces but comparatively less research work is available highlighting curved stretched surfaces.This gets even narrower when we talk about the study of hybrid nanoliquid with entropy optimization past curved surfaces.Therefore, our task here is to discuss hybrid nanoliquid flow comprising ferromagnetic nanoparticle, i.e., nickel-zinc ferrite (NiZnFe 2 O 4 ), and the base fluid, i.e., ethylene glycol (C 2 H 6 O 2 ), over a curved surface with entropy optimization coating.The whole analysis was performed with added impressions of nonlinear thermal radiation with entropy optimization coatings.The analysis was supported by the convective heat and mass flux boundary conditions.Numerical solution of the envisioned model was obtained by utilizing bvp4c from MATLAB.The traits of the sundry parameters on involved distributions were thoroughly discussed keeping their physical justification in mind.

Mathematical Formulation
We considered a 2D steady, incompressible nanoliquid flow over a curved stretching channel looped in the form of a circle with a radius R about the curvilinear directions r and x, as shown in Figure 1.Here, a higher value of R corresponds to a marginally curved surface.The stretching velocity is taken as u = u w along the x-direction.A magnetic field is applied normal to the fluid flow and along the r-direction.The electric and induced magnetic fields were overlooked owing to our assumption of small Reynolds number.hybrid nanoliquid with entropy optimization past curved surfaces.Therefore, our task here is to discuss hybrid nanoliquid flow comprising ferromagnetic nanoparticle, i.e., nickel-zinc ferrite (NiZnFe2O4), and the base fluid, i.e., ethylene glycol (C2H6O2), over a curved surface with entropy optimization coating.The whole analysis was performed with added impressions of nonlinear thermal radiation with entropy optimization coatings.The analysis was supported by the convective heat and mass flux boundary conditions.Numerical solution of the envisioned model was obtained by utilizing bvp4c from MATLAB.The traits of the sundry parameters on involved distributions were thoroughly discussed keeping their physical justification in mind.

Mathematical Formulation
We considered a 2D steady, incompressible nanoliquid flow over a curved stretching channel looped in the form of a circle with a radius R about the curvilinear directions r and x, as shown in Figure 1.Here, a higher value of R corresponds to a marginally curved surface.The stretching velocity is taken as u = uw along the x-direction.A magnetic field is applied normal to the fluid flow and along the r-direction.The electric and induced magnetic fields were overlooked owing to our assumption of small Reynolds number.The assumed system is governed by the following equations: (1) ( ) ( ) ) The system of Equations ( 1)-( 5) is supported by the following boundary conditions: , The thermophysical traits of the hybrid nanoliquid are appended in Table 1.The assumed system is governed by the following equations: v ∂u ∂r v ∂C ∂r

)
(5) The system of Equations ( 1)-( 5) is supported by the following boundary conditions: = j w u| r→∞ → 0, ∂u ∂r r→∞ The thermophysical traits of the hybrid nanoliquid are appended in Table 1.The mathematical form of thermophysical properties are given as follows: In Equation ( 4), the nonlinear radiation heat flux term via Rosseland's approximation is given as follows:

Solution Procedure
Here, we used the following dimensionless transformations: Here, prime denotes the derivative w, r, T, ζ and θ w = T f /T ∞ .The above transformation Equation (11) satisfies Equation (1) identically and Equations ( 2)-( 6) are given by the following: 1 Pr Here, Eliminating pressure term from Equations ( 12) and ( 13) by differentiating Equation ( 13) with respect to ζ and then putting in Equation ( 12), we get the following: The surface drag force (C f ), Sherwood number (Sh x ), and Nu x (Nusselt number) along x-direction are defined as follows: where After putting Equations ( 11) and ( 20), Equation ( 19) becomes the following: Here, Re x = sx 2 v f .

Entropy Generation
The equation of entropy generation in dimensional form is given in reference [34].

∂C ∂r
∂T ∂r (22) Equation ( 22) after applying the requisite transformations takes the following form: where

Results and Discussion
The MATLAB built-in function bvp4c was applied to integrate the numerical solution for the system of Equations ( 14), (15), and ( 17) with initial and boundary conditions, Equations ( 16) and (18), for numerous values of K 1 , M, R d , λ 1 , and S c graphically.For this technique, we first changed differential equations with the higher order to the equations of order one by utilizing new variables.The function bvp4c needs an initial guess for the solution and with the tolerance of 10 −7 .The guess we chose needed to satisfy the boundary conditions (Equations ( 16) and ( 18)) and the solution.The validation of our presented results is depicted in Table 2.An excellent agreement with Sanni et al. [23] was observed when M = 1, ϕ = 0.0, and in the absence of temperature and concentration profile.we chose needed to satisfy the boundary conditions (Equations ( 16) and ( 18)) and the solution.The validation of our presented results is depicted in Table 2.An excellent agreement with Sanni et al. [23] was observed when M = 1, φ = 0.0, and in the absence of temperature and concentration profile.The impact of the curvature parameter K1 on velocity, concentration, and temperature profiles are depicted in Figures 4-6.Increasing values of K1 resulted in an increase in fluid velocity and temperature field, while the concentration profile diminished.This was because of the radius of the surface augment when curvature parameter K1 was increased.As a result, the flow increased but it offered more resistance, therefore the temperature rose.The values of other parameters were fixed as Pr = 10, Sc = 0.5, φ = 0.1, Rd = 0.5, θw = 0.5, λ1 = 0.5 and M = 0.1.we chose needed to satisfy the boundary conditions (Equations ( 16) and ( 18)) and the solution.The validation of our presented results is depicted in Table 2.An excellent agreement with Sanni et al. [23] was observed when M = 1, φ = 0.0, and in the absence of temperature and concentration profile.The impact of the curvature parameter K1 on velocity, concentration, and temperature profiles are depicted in Figures 4-6.Increasing values of K1 resulted in an increase in fluid velocity and temperature field, while the concentration profile diminished.This was because of the radius of the surface augment when curvature parameter K1 was increased.As a result, the flow increased but it offered more resistance, therefore the temperature rose.The values of other parameters were fixed as Pr = 10, Sc = 0.5, φ = 0.1, Rd = 0.5, θw = 0.5, λ1 = 0.5 and M = 0.1.The impact of the curvature parameter K 1 on velocity, concentration, and temperature profiles are depicted in Figures 4-6.Increasing values of K 1 resulted in an increase in fluid velocity and temperature field, while the concentration profile diminished.This was because of the radius of the surface augment when curvature parameter K 1 was increased.As a result, the flow increased but it offered more resistance, therefore the temperature rose.The values of other parameters were fixed as Pr = 10, S c = 0.5, ϕ = 0.1, R d = 0.5, θ w = 0.5, λ 1 = 0.5 and M = 0.1.The characteristics of Biot number Bi and heat generation/absorption parameter λ1 on temperature field are displayed in Figures 8 and 9, respectively.Figure 8 shows that the convective heat transfer coefficient intensified for higher estimates of Bi and the temperature subsequently rose.The characteristics of Biot number Bi and heat generation/absorption parameter λ 1 on temperature field are displayed in Figures 8 and 9, respectively.Figure 8 shows that the convective heat transfer coefficient intensified for higher estimates of Bi and the temperature subsequently rose.The characteristics of Biot number Bi and heat generation/absorption parameter λ1 on temperature field are displayed in Figures 8 and 9, respectively.Figure 8 shows that the convective heat transfer coefficient intensified for higher estimates of Bi and the temperature subsequently rose.The characteristics of Biot number Bi and heat generation/absorption parameter λ1 on temperature field are displayed in Figures 8 and 9, respectively.Figure 8 shows that the convective heat transfer coefficient intensified for higher estimates of Bi and the temperature subsequently rose.The impression of Schmidt number Sc on concentration distribution is portrayed in Figure 12.A decrease in concentration field was detected with increasing values of Sc.As the Schmidt number has a converse proportion with the Brownian diffusion coefficient, an increment in the Sc yielded a decay in Brownian diffusion coefficient that brought about a diminishment in concentration and its interrelated boundary layer thickness.The values allocated to other parameters were K1 = 10, M = 0.3, φ = 0.1, Bi = 0.1, θw = 0.5, and Rd = 0.1.

Figure 12. Upshot of Sc on concentration profile h(ζ).
The influence of curvature parameter K1 and magnetic parameter M on skin friction coefficient   The impression of Schmidt number Sc on concentration distribution is portrayed in Figure 12.A decrease in concentration field was detected with increasing values of Sc.As the Schmidt number has a converse proportion with the Brownian diffusion coefficient, an increment in the Sc yielded a decay in Brownian diffusion coefficient that brought about a diminishment in concentration and its interrelated boundary layer thickness.The values allocated to other parameters were K1 = 10, M = 0.3, φ = 0.1, Bi = 0.1, θw = 0.5, and Rd = 0.1.

Figure 12. Upshot of Sc on concentration profile h(ζ).
The influence of curvature parameter K1 and magnetic parameter M on skin friction coefficient is depicted in Figure 13.It can be noticed that the surface drag force diminished with increasing value of K1.A contradictory trend was demonstrated in case of M. In Figure 14, the consequence of magnetic parameter M and solid volume fraction φ on shear wall stress is demonstrated.The skin friction profile rose with increase in magnetic parameter M and solid volume fraction φ for fixed values of parameters Pr = 10, Sc = 0.5, Rd = 0.5, θw = 0.5, λ1 = 0.5 and Bi = 0.1.

Upshot of S c on concentration profile h(ζ).
The influence of curvature parameter K 1 and magnetic parameter M on skin friction coefficient x is depicted in Figure 13.It can be noticed that the surface drag force diminished with increasing value of K 1 .A contradictory trend was demonstrated in case of M. In Figure 14, the consequence of magnetic parameter M and solid volume fraction ϕ on shear wall stress is demonstrated.The skin friction profile rose with increase in magnetic parameter M and solid volume fraction ϕ for fixed values of parameters Pr = 10, S c = 0.5, R d = 0.5, θ w = 0.5, λ 1 = 0.5 and Bi = 0.1.The impact of the Brinkman number (Br) on entropy generation is portrayed in Figure 17.From the illustration, it can be seen that entropy optimization was boosted with increasing values of (Br).The reason behind this was that more heat was generated between the layers of the fluid because of augmented values of (Br). Figure 18 displays the relationship between the magnetic parameter (M) and the entropy generation.Again, the same trait as depicted in case of (Br) was witnessed here.Higher values of (M) meant stronger Lorentz force and ultimate strengthening of the dissipation energy, and this was the main cause of irreversibility.Table 3 shows the behavior of Sherwood number  x .
The impact of the Brinkman number (Br) on entropy generation is portrayed in Figure 17.From the illustration, it can be seen that entropy optimization was with increasing values of (Br).The reason behind this was that more heat was generated between the layers of the fluid because of augmented values of (Br). Figure 18 displays the relationship between the magnetic parameter (M) and the entropy generation.Again, the same trait as depicted in case of (Br) was witnessed here.Higher values of (M) meant stronger Lorentz force and ultimate strengthening of the dissipation energy, and this was the main cause of irreversibility.The impact of the Brinkman number (Br) on entropy generation is portrayed in Figure 17.From the illustration, it can be seen that entropy optimization was boosted with increasing values of (Br).The reason behind this was that more heat was generated between the layers of the fluid because of augmented values of (Br). Figure 18 displays the relationship between the magnetic parameter (M) and the entropy generation.Again, the same trait as depicted in case of (Br) was witnessed here.Higher values of (M) meant stronger Lorentz force and ultimate strengthening of the dissipation energy, and this was the main cause of irreversibility.Table 3 shows the behavior of Sherwood number  The impact of the Brinkman number (Br) on entropy generation is portrayed in Figure 17.From the illustration, it can be seen that entropy optimization was boosted with increasing values of (Br).The reason behind this was that more heat was generated between the layers of the fluid because of augmented values of (Br). Figure 18 displays the relationship between the magnetic parameter (M) and the entropy generation.Again, the same trait as depicted in case of (Br) was witnessed here.Higher values of (M) meant stronger Lorentz force and ultimate strengthening of the dissipation energy, and this was the main cause of irreversibility.Table 3 shows the behavior of Sherwood number  Table 3 shows the behavior of Sherwood number Sh x (Re x ) − 1 2 for varied values of S c (Schmidt number), K 1 (curvature parameter), and M (magnetic parameter).It can be seen that for snowballing values of S c , the Sherwood number Sh x (Re x ) − 1 2 increased; however, it diminished for increasing value of K 1 and M.

Concluding Remarks
In this paper, the flow of nanoliquid comprising nickel-zinc ferrite-ethylene glycol (NiZnFe 2 O 4 -C 2 H 6 O 2 ) accompanied by entropy optimization coating past a curved stretching surface with convective heat and mass flux boundary was examined.The solution of the envisioned system of equations was found numerically by applying MATLAB built-in function bvp4c.The impact of numerous arising parameters on involved profiles was depicted via graphical illustrations with requisite discussions.
The conclusions of the current study are as follows: • An increase in curvature parameter accounted for increasing velocity and temperature fields and diminishing concentration distribution.

•
Under the considerable influence of magnetic parameter, an increased axial velocity field was attained.

•
For the increasing estimates of the solid volume fraction, the temperature and velocity profiles showed increasing behavior.

•
The temperature profile improved with increasing values of Biot number and heat generation/absorption parameter.

•
The value of friction factor profile augmented for larger values of M and ϕ, but it decreased for K 1 and M.

•
The Nusselt number NuRe −1/2 x declined with increasing values of K 1 and R d .

Figures 2
Figures 2 and 3 show the impression of solid volume fraction ϕ on velocity and temperature profiles.Both velocity fields increased with increasing values of solid volume fraction ϕ.Further, the momentum and thickness of the thermal boundary layers were boosted with a larger value of ϕ.The values given to other parameters were Pr = 10, S c = 0.5, K 1 = 10, R d = 0.5, θ w = 0.5, λ 1 = 0.5, and M = 0.1.

Figures 2
Figures 2 and 3 show the impression of solid volume fraction φ on velocity and temperature profiles.Both velocity fields increased with increasing values of solid volume fraction φ.Further, the momentum and thickness of the thermal boundary layers were boosted with a larger value of φ.The values given to other parameters were Pr = 10, Sc = 0.5, K1 = 10, Rd = 0.5, θw = 0.5, λ1 = 0.5, and M = 0.1.

Figures 2
Figures 2 and 3 show the impression of solid volume fraction φ on velocity and temperature profiles.Both velocity fields increased with increasing values of solid volume fraction φ.Further, the momentum and thickness of the thermal boundary layers were boosted with a larger value of φ.The values given to other parameters were Pr = 10, Sc = 0.5, K1 = 10, Rd = 0.5, θw = 0.5, λ1 = 0.5, and M = 0.1.

Figure 7
Figure 7 demonstrates the variation in the velocity field for numerous estimates of magnetic parameter M. Here, increments in M led to a decline in the magnitude of fluid's velocity.This was because of the resistive force (called Lorentz force) triggered by the magnetic field, which lowered the velocity of the fluid's velocity flow.The values of the other parameters were fixed as K1 = 10, Sc = 0.5, φ = 0.1, Rd = 0.5, θw = 0.5, λ1 = 0.5, and Bi = 0.1.

Figure 7
Figure 7 demonstrates the variation in the velocity field for numerous estimates of magnetic parameter M. Here, increments in M led to a decline in the magnitude of fluid's velocity.This was because of the resistive force (called Lorentz force) triggered by the magnetic field, which lowered the velocity of the fluid's velocity flow.The values of the other parameters were fixed as K1 = 10, Sc = 0.5, φ = 0.1, Rd = 0.5, θw = 0.5, λ1 = 0.5, and Bi = 0.1.

Figure 7
Figure 7 demonstrates the variation in the velocity field for numerous estimates of magnetic parameter M. Here, increments in M led to a decline in the magnitude of fluid's velocity.This was because of the resistive force (called Lorentz force) triggered by the magnetic field, which lowered the velocity of the fluid's velocity flow.The values of the other parameters were fixed as K1 = 10, Sc = 0.5, φ = 0.1, Rd = 0.5, θw = 0.5, λ1 = 0.5, and Bi = 0.1.

Figure 7 15 Figure 7 .
Figure 7 demonstrates the variation in the velocity field for numerous estimates of magnetic parameter M. Here, increments in M led to a decline in the magnitude of fluid's velocity.This was because of the resistive force (called Lorentz force) triggered by the magnetic field, which lowered the velocity of the fluid's velocity flow.The values of the other parameters were fixed as K 1 = 10, S c = 0.5, ϕ = 0.1, R d = 0.5, θ w = 0.5, λ 1 = 0.5, and Bi = 0.1.

Figure 9
illustrates the behavior of λ1.To increase the estimation values of heat absorption/generation parameter, the temperature profile and thermal boundary layer thickness were increased.The values of other parameters were fixed as K1 = 10, M = 0.3, φ = 0.5, Rd = 0.1, θw = 0.5, and Pr = 10.

Figures 10 and 11
Figures 10 and 11 show the impacts of nonlinear radiation parameter Rd and Prandtl number Pr on temperature distribution, respectively.It can be seen that the temperature field fell with increasing Prandtl number Pr.As Prandtl number is linked in a reciprocal way to the thermal diffusivity, a quick augmentation in the Prandtl number Pr lessened the temperature and thickness of the thermal boundary layer.The temperature profile increased with increment in nonlinear radiation parameter Rd.Physically, the radiative heat flux increased with increasing values of Rd which ultimately boosted the temperature of the fluid.The values assigned to other parameters were K1 = 10, M = 0.3, φ = 0.5, Bi = 0.1, θw = 0.5, and Sc = 5.

Figure 9
illustrates the behavior of λ1.To increase the estimation values of heat absorption/generation parameter, the temperature profile and thermal boundary layer thickness were increased.The values of other parameters were fixed as K1 = 10, M = 0.3, φ = 0.5, Rd = 0.1, θw = 0.5, and Pr = 10.

Figures 10 and 11
Figures 10 and 11 show the impacts of nonlinear radiation parameter Rd and Prandtl number Pr on temperature distribution, respectively.It can be seen that the temperature field fell with increasing Prandtl number Pr.As Prandtl number is linked in a reciprocal way to the thermal diffusivity, a quick augmentation in the Prandtl number Pr lessened the temperature and thickness of the thermal boundary layer.The temperature profile increased with increment in nonlinear radiation parameter Rd.Physically, the radiative heat flux increased with increasing values of Rd which ultimately boosted the temperature of the fluid.The values assigned to other parameters were K1 = 10, M = 0.3, φ = 0.5, Bi = 0.1, θw = 0.5, and Sc = 5.

Figure 9
illustrates the behavior of λ1.To increase the estimation values of heat absorption/generation parameter, the temperature profile and thermal boundary layer thickness were increased.The values of other parameters were fixed as K1 = 10, M = 0.3, φ = 0.5, Rd = 0.1, θw = 0.5, and Pr = 10.

Figures 10 and 11
Figures 10 and 11 show the impacts of nonlinear radiation parameter Rd and Prandtl number Pr on temperature distribution, respectively.It can be seen that the temperature field fell with increasing Prandtl number Pr.As Prandtl number is linked in a reciprocal way to the thermal diffusivity, a quick augmentation in the Prandtl number Pr lessened the temperature and thickness of the thermal boundary layer.The temperature profile increased with increment in nonlinear radiation parameter Rd.Physically, the radiative heat flux increased with increasing values of Rd which ultimately boosted the temperature of the fluid.The values assigned to other parameters were K1 = 10, M = 0.3, φ = 0.5, Bi = 0.1, θw = 0.5, and Sc = 5.

Figures 10 and 11 15 Figure 10 .
Figures 10 and 11 show the impacts of nonlinear radiation parameter R d and Prandtl number Pr on temperature distribution, respectively.It can be seen that the temperature field fell with increasing Prandtl number Pr.As Prandtl number is linked in a reciprocal way to the thermal diffusivity, a quick augmentation in the Prandtl number Pr lessened the temperature and thickness of the thermal boundary layer.The temperature profile increased with increment in nonlinear radiation parameter R d .Physically, the radiative heat flux increased with increasing values of R d which ultimately boosted
depicted in Figure 13.It can be noticed that the surface drag force diminished with increasing value of K1.A contradictory trend was demonstrated in case of M. In Figure 14, the consequence of magnetic parameter M and solid volume fraction φ on shear wall stress is demonstrated.The skin friction profile rose with increase in magnetic parameter M and solid volume fraction φ for fixed values of parameters Pr = 10, Sc = 0.5, Rd = 0.5, θw = 0.5, λ1 = 0.5 and Bi = 0.1.

Figure 11 . 15 Figure 10 .
Figure 11.Upshot of Pr on temperature field θ(ζ).The impression of Schmidt number S c on concentration distribution is portrayed in Figure12.A decrease in concentration field was detected with increasing values of S .As the Schmidt number has a converse proportion with the Brownian diffusion coefficient, an increment in the S c yielded a decay in Brownian diffusion coefficient that brought about a diminishment in concentration and its interrelated boundary layer thickness.The values allocated to other parameters were K 1 = 10, M = 0.3, ϕ = 0.1, Bi = 0.1, θ w = 0.5, and R d = 0.1.

Figure 14 .Figure 15 15 Figure 13 .
Figure 14.Upshot of M and ϕ on wall shear stress − 1 2 C f Re 1/2 x .Figure15shows the effect of Biot number Bi and solid volume fraction ϕ on Nusselt number Nu x (Re x ) − 1 2 .It was detected that for higher value of Bi and ϕ, the surface heat transfer rate upsurged when values of parameters were given as K 1 = 10, M = 0.3, S c = 5.0, R d = 0.1, θ w = 0.5, λ 1 = 0.5 and Pr = 10.

15 .
Upshot of Bi and on ϕ Nusselt number NuRe −1/2 x .The outcome of curvature parameter K 1 and nonlinear radiation parameter R d on Nusselt number Nu x (Re x ) − 1 2 is examined in Figure 16.Here, a reduction in Nusselt number was noted for increasing values of curvature parameter K 1 and the opposite trend was seen for nonlinear radiation parameter R d for fixed values of ϕ = 0.1, M = 0.3, S c = 5.0, θ w = 0.5, λ 1 = 0.5 and Pr = 10.

Figure 16 .
Figure 16.Upshot of K1 and Rd on Nusselt number

−
for varied values of Sc (Schmidt number), K1 (curvature parameter), and M (magnetic parameter).It can seen that for snowballing values of Sc, the Sherwood number 1 2 Sh (Re ) x x − increased; however, it diminished for increasing value of K1 and M.

−
for varied values of Sc (Schmidt number), K1 (curvature parameter), and M (magnetic parameter).It can be seen that for snowballing values of Sc, the Sherwood number 1 2 Sh (Re ) x x − increased; however, it diminished for increasing value of K1 and M.

−
for varied values of Sc (Schmidt number), K1 (curvature parameter), and M (magnetic parameter).It can be seen that for snowballing values of Sc, the Sherwood number 1 2 Sh (Re ) x x − increased; however, it diminished for increasing value of K1 and M.

Figure 18 .
Figure 18.Upshot of M on N G

Table 2 .
Comparison of presented results for skin friction coefficient 1 2 C f (Re x )

Table . 2
. Comparison of presented results for skin friction coefficient

Table . 2
. Comparison of presented results for skin friction coefficient