The impact of a chemical reaction on the heat and mass transfer mechanisms in a dissipative and radiative nanofluid flow over a nonlinear stretching sheet

This paper presents a numerical investigation of the flow of a non-Newtonian tangent hyperbolic nanofluid over a nonlinearly stretched surface, taking into account factors such as thermal radiation, prescribed surface temperature, and a chemical reaction mechanism. Furthermore, the analysis includes the consideration of both viscous dissipation and the influence of a magnetic field within a Darcy porous medium. A mathematical framework for addressing the issue, rooted in the principles of conserving momentum, energy, and mass. The MATHEMATICA tools were employed to apply the shooting technique in order to solve the modeled equations describing the temperature, velocity, and concentration fields of the proposed physical system. Graphs are used to illustrate how certain key parameters affect the profiles of concentration, velocity, and temperature. Data tables are utilized to display information pertaining to the local Nusselt number, local Sherwood number, and local skin friction coefficient. The present results have been confirmed through a comparison with previously published findings. This research holds significant importance as it focuses on the extensive utilization of tangent hyperbolic nanofluids in cooling electronic components that produce substantial heat during their operation. The observed pattern indicates that as the local Weisbsenberg number, magnetic number, local porous parameter, and power law index increase, there is a reduction in the boundary layer thickness. Conversely, in the instances of concentration and temperature distributions, an escalation in these parameters leads to an expansion of the boundary layer thickness.


Mathematical modeling
In this model, we consider the steady flow of a nanofluid with a hyperbolic tangent profile, taking into consideration its variable conductivity as it passes over a nonlinearly stretched sheet.The stretched sheet is oriented horizontally, with its length along the x-axis and the y-axis positioned perpendicular to it.The fluid's movement is a result of the stretching of the horizontal sheet.Additionally, we make the assumption that the flow is confined to the region where y is greater than zero.The physical model is presented in a systematic manner in Fig. 1.
A constant magnetic field B in the transverse direction is applied along the y-axis, and the sheet is stretched in the x-axis direction with a velocity represented as u w = ax 1 3 .Also, the constant magnetic field B is assumed to vary with x coordinate in the form B = B 0 x −2 3 .Also, T w , T ∞ and C w represent the temperature at the wall, the ambient temperature, and the concentration at the wall's surface, respectively.We analyze the heat and mass transfer aspects, taking into account the influence of a chemical reaction phenomenon characterized by a conversion rate K and the presence of thermal radiation with a radiative flux denoted as q r .The porous medium that encompasses the stretching sheet follows the fundamental principles of Darcy's law, which govern fluid flow through porous materials.In the context presented here, the discussion takes into account thermophoresis effects and also investigates both viscous dissipation and the Brownian motion of nanoparticles that are suspended within the fluid.In their dimensional form, the governing equations, as described in references 31,32 , are expressed as follows: Here, u, v, c p , µ, D T , D B , C, ρ, Ŵ, n and τ correspond to the longitudinal and perpendicular velocity components, specific heat, kinematic viscosity, thermophoresis diffusion coefficient, Brownian diffusivity, concentration of nanoparticles, fluid density, time constant, power law index and the ratio between the heat capacity of the nanomaterial and that of the fluid.The applicable constraints or limits can be summarized as follows: where v w represents the velocity of suction which can be assumed to vary with x coordinate in the form v w = v 0 x −1 3 .The transformations that follow 33 are introduced with the aim of converting the represented equations into a dimensionless form: In this context, the variable q r is streamlined and expressed as a function of temperature, which aligns with the prior research conducted by Alali and Megahed 25 .In this context, it is important to note that the applicability of linear radiation, based on a linear correlation between temperature and radiation heat transfer, is generally suitable for moderate temperature ranges.Consequently, our assumptions are aligned with scenarios characterized by moderate temperature distributions.Furthermore, we embrace the presumption that the thermal conductivity, denoted as κ , for the Tangent hyperbolic nanofluid varies as a function of dimensionless temperature θ .This connection is elucidated by the subsequent equations 34 : where ε represents the parameter governing thermal conductivity.After implementing the transformations mentioned previously in Eqs. ( 7)-( 8), the governing equations take on the ensuing structure: In addition to the following given boundary conditions: (1) ( where Re x = u w x ν is the local Reynolds number.

Method of solution
Owing to the intricate nonlinearity inherent in Eq. ( 10) through (12), our present physical model does not possess an exact solution.Nevertheless, employing the shooting technique enabled us to convert these equations from boundary value problems to initial value problems.Subsequently, we successfully solved the system numerically by employing an appropriate initial guess.The difficulty encountered revolved around determining suitable initial values and conditions for the modified problem.In order to overcome this challenge, a transformation of the following is necessary: In light of the preceding information, it is possible to deduce the following from the set of Eqs. ( 10)- (12): The iterative shooting method is utilized to seek appropriate estimates for ξ 1 , ξ 2 , and ξ 3 , ensuring the fulfillment of outer boundary conditions χ 2 (∞) , χ 4 (∞) , and χ 6 (∞) .This procedure can be reiterated until the target level of precision, specifically 10 −6 , is attained.

Confirmation of the code's precision
In this study, we utilized the shooting technique in conjunction with the Runge-Kutta scheme to conduct our numerical computations.The shooting method presents numerous advantages, such as its adaptability in addressing boundary value problems related to ordinary differential equations.It is versatile across a broad spectrum of problems and is easy to implement.Notably, the method demonstrates numerical efficiency and doesn't necessitate a detailed formulation of the solution, enhancing its utility for intricate systems.To verify the effectiveness of the provided code, we obtained results for the skin friction coefficient Re 1 2 x Cf x at β = f w = 0 for various values of the power law index n while maintaining the other parameters at constant values, as illustrated in Table 1.The outcomes displayed in Table 1 indicate a strong concurrence between our findings and those acquired in the study conducted by Ahmed et al. 33 .This strong agreement serves to confirm both the precision and the reliability of our code and the numerical outcomes it produces.It is important to note that the originality and novelty of our proposed research, in comparison to the previous work by Ahmed et al. 33 , can be encapsulated by examining crucial phenomena such as viscous dissipation, variable thermal conductivity, and the chemical reaction process within a Darcy porous medium.

Discussion
This portion of the study delves into the numerical analysis of how governing parameters affect the flow behavior of a magnetized Hyperbolic tangent nanofluid over a non-linear stretching sheet.The equations in this model are controlled by various key parameters that have recently emerged in the field.These parameters include the local porous parameter denoted as β , the local suction factor represented as f w , the local Hartmann number indicated as M, the local Weisbsenberg number denoted as Ŵ e , the local Eckert number labeled as Ec, the thermophoresis parameter referred to as δ t , the radiation parameter designated as R, the Brownian motion parameter represented as δ b , the Schmidt number termed as Sc, and the chemical reaction parameter known as γ .Figure 2 illustrates how the temperature θ(η) , concentration φ(η) , and velocity f ′ (η) profiles change with respect to η as the local Weisbsenberg number Ŵ e undergoes variations.The Weissenberg number Ŵ e is a dimensionless quantity that quantifies the relationship between the relaxation time of a material and the time scale associated with fluid flow.As depicted in Fig. 2, as Ŵ e increases, it becomes evident that a longer duration is required for the fluid to undergo its flow, leading to a reduction in velocity.Conversely, both temperature and concentration exhibit an opposite trend in response to the variations in Ŵ e .The influence of Ŵ e on the characteristics of the flow of the nanofluid can be attributed to the fundamental concept behind Ŵ e .Physically, the local Weisbsenberg number fundamentally assesses the speed at which the fluid moves in relation to the material's capacity to deform or respond to external forces.
Figure 3 presents a visualization of the traits associated with the porous parameter β and how it affects the flow properties of the nanofluid.A dimensionless quantity that deals with flows across porous medium is the porous parameter.It describes the porosity or permeability of a substance through which a fluid is flowing.This characteristic is crucial because it affects how easily and with what resistance fluid can pass through porous objects like geological formations, filters, and membranes.Hence, Fig. 3 effectively illustrates that as the porous x Cf x with the previously findings of Ahmed et al. 33 when M = Ŵ e = 0.5 and β = f w = 0. parameter experiences an increase, there is a noteworthy reduction in the velocity of the nanofluid.Additionally, it is evident that the resistance of the fluid to the pores within the porous medium substantially amplifies both the temperature and concentration of the fluid.Physically, the decline in nanofluid velocity due to the porous parameter is mainly linked to the heightened resistance faced by the fluid when moving through the porous medium.This increased resistance acts as an obstacle to the fluid's flow, leading to a decrease in velocity.
Figure 4 provides a visual representation of how the Magnetic parameter M influences the characteristics of flow, heat, and mass within the nanofluid.It is evident that the fluid's velocity experiences a reduction in the presence of a magnetic field.Furthermore, the same graph demonstrates that the presence of a magnetic force leads to an improvement in both the temperature θ(η) and concentration φ(η) distributions.In a physical sense, the Lorentz force exerted on the fluid increases with the rise in the magnetic parameter.Consequently, this causes a reduction in fluid velocity, ultimately slowing down the fluid's movement.Moreover, the preceding research conducted by Yasmin and colleagues 35 provides support and confirmation for the influence of the magnetic field on the heat transfer dynamics within nanofluid flow.
Figure 5 provides a graphical depiction of how changes in the suction parameter impact the characteristics of concentration φ(η) , velocity f ′ (η) , and temperature θ(η) profiles.This graph demonstrates that with an increase in the suction parameter f w , the concentration, temperature, and velocity profiles decline at elevated values.From a physical perspective, the reduction in nanofluid velocity with an increase in the suction parameter can be attributed to the presence of pores on the stretching sheet.These pores introduce resistance to the flow of fluid, leading to the observed decrease in velocity.
Figure 6 illustrates the changes in the graphs of temperature, concentration, and velocity concerning variations in the power law index n.This visual representation showcases how these parameters evolve as the power law index undergoes alterations.It is quite apparent that there is a clear correlation between the increase in the power law index factor and the corresponding elevation of the fluid's temperature, concentration, and the thickness of the thermal boundary layer.Furthermore, based on the information presented in the same figure, it becomes quite evident that when the power law index is intensified, there is a noticeable reduction in the velocity distribution within the boundary layer region.Physically, the reduction in nanofluid velocity associated with the power law index is mainly due to the fluid's response to shear forces.A higher power law index indicates heightened resistance to shear, resulting in a decline in fluid velocity.Figure 7a provides a visual representation of how the thermal radiation parameter R influences the temperature profile θ(η) .Essentially, the radiation parameter assists in assessing how significant thermal radiation is in a specific heat transfer situation, providing insight into whether it plays a prominent or inconsequential role.Consequently, with an increase in the radiation parameter, there is a notable upward trend in θ(η) , resulting in the expansion of the thermal boundary-layer thickness.Additionally, previous studies [36][37][38][39][40][41] provide confirmation of the influence of the thermal radiation parameter on the heat transfer mechanism.Figure 7b provides a depiction of how the temperature distribution, represented by θ(η) , is influenced by changes in the variable thermal conductivity factor ε .As the thermal conductivity factor rises, the temperature of the nanofluid rises.Physically,  www.nature.com/scientificreports/this phenomenon happens as a result of a more effective heat transmission from the sheet to the nanofluid at higher levels of the thermal conductivity factor.Figure 8 visually demonstrates how variations in the thermophoretic parameter δ t influence both temperature θ(η) and concentration φ(η) .To begin, we introduce the concept of the thermophoretic parameter, a non- dimensional value relevant to thermophoresis.Thermophoresis is a phenomenon occurring in fluid dynamics and heat transfer where particles or molecules dispersed within a fluid experience motion due to temperature variations within the fluid.An escalation in the thermophoretic parameter, denoted as δ t , is depicted as leading to an increase in both temperature θ(η) and concentration φ(η) levels.Moreover, the rise in temperature and concentration resulting from alterations in the thermophoretic parameter leads to an increase in the thickness of the boundary layer.Furthermore, as the value of the thermophoretic parameter increases, there is a corresponding decline in the local Nusselt number, while the local Sherwood number exhibits an opposite trend.
Figure 9 portrays the concentration and temperature profiles for different values of Brownian motion parameter δ b .This graph clearly illustrates a direct correlation between the Brownian motion parameter and the profiles of both temperature and concentration.Elevating the Brownian motion parameter leads to a decrease in nanofluid concentration.This occurs because Brownian motion causes the particles to move away from the fluid region as the boundary layer experiences warming.Physically, an elevation in the Brownian motion parameter leads to increased thermal conduction as a result of the enhanced mobility of nanoparticles, thereby causing a rise in the temperature field.
The temperature profile θ(η) is depicted graphically in Fig. 10a in relation to variations in the Eckert number Ec.The correlation between the Eckert number Ec and the temperature distribution is highlighted by this visual representation.The relationship between the Eckert number and temperature differences is inverse; as the Eckert number increases, temperature differences become weaker, resulting in more pronounced temperature profiles and a thicker thermal boundary layer.Figure 10b illustrates the attributes of the chemical reaction parameter γ and how it influences the distribution of concentration φ(η) .A decrease in the concentration of the nanofluid results from increasing the chemical reaction parameter.This effect results from a drop in nanofluid concentration brought on by an enhanced chemical reaction parameter.
Table 2   www.nature.com/scientificreports/summary of how these diverse factors affect these critical physical quantities.As the power law index and the Weissenberg number increase, there is a corresponding decrease in the local Nusselt number, skin friction coefficient, and local Sherwood number.Also, the local Sherwood number, skin friction coefficient, and local Nusselt number all rise in response to a boost in the suction parameter.As the magnetic number is heightened, the skin friction coefficient experiences an increase.However, the local Sherwood number and the local Nusselt number

Concluding remarks
This research offers a quantitative examination of the flow of a non-Newtonian tangent hyperbolic nanofluid over a surface that undergoes nonlinear stretching, while the flow takes place through a porous medium characterized by Darcy's law.This study takes into account several factors, including chemical reactions, variable thermal conductivity, the magnetic field, the dissipation of viscous energy, and the influence of thermal radiation.After using the shooting technique to obtain the solution for our model, the numerical outcomes are presented in both tabular and graphical forms to emphasize the specific details of flow and heat transfer.Below is a compilation of the model's most crucial findings.
1.The velocity profile is diminished by factors such as the porous parameter, magnetic number, suction parameter, power law index, and the Weisbsenberg number.2. The skin friction coefficient exhibits a declining trend as the Weisbsenberg number and power law index rise, while it demonstrates an upward trend with increased values of the porous parameter, magnetic parameter, and suction parameter.3. The nanofluid temperature experiences a decrease when the suction parameter and Prandtl number are increased, whereas it shows an increase as the Eckert number, Brownian motion parameter, radiation parameter, power law index, and magnetic number are increased.4. Heat and mass transfer rates decline with an increase in the Weisbsenberg number, porous parameter, magnetic number, and power law index.Conversely, the suction parameter leads to the opposite effect, causing an increase in these rates.5.The concentration distribution exhibits a decrease with an increase in both the Brownian motion parameter and suction parameter, whereas it shows an increase with higher values of the power law index and magnetic number.

Figure 1 .
Figure 1.Geometry of nanofluid flow over a stretching sheet.
as η → ∞.It's evident that the resulting governing parameters are clearly defined as follows:where β is the local porous parameter, f w is the local suction factor, M is the local Hartmann number, Ŵ e is the local Weisbsenberg number, Ec is the local Eckert number, δ t is the thermophoresis parameter, R is the radiation parameter, δ b is the Brownian motion parameter, Sc is the Schmidt number, γ is the chemical reaction parameter and Pr is the Prandtl Number.Furthermore, we can mathematically express certain distinct physical factors to provide an intriguing insight into mass transfer, flow characteristics, and heat transfer mechanisms.These quantities are denoted as the local Sherwood number Sh x , the local skin-friction coefficient Cf x , and the local Nusselt number Nu x .Their respective definitions are as follows:

Table 1 .
Comparison table of Re presents the influences of various significant factors on key physical parameters, including the skin friction coefficient Cf x Re

Table 2 .
Values of Cf x Refor various values of Ŵ e , β, M, R, n, ε, δ t , δ b , Ec and δ b with Sc = 2.0, γ = 0.2 and Pr = 3.0.because of the Lorentz effect, which hinders fluid movement.As the value of the Brownian motion parameter rises, there is a concurrent decrease in the local Nusselt number, while the local Sherwood number experiences a simultaneous increase.