Entropy generation minimization of higher

Scientific Reports volume 12, Article number: 17688 (2022) Cite this article

1130 Accesses

5 Citations

Metrics details

The present investigation aims to analyze higher-order endothermic/exothermic chemical reactions with activation energy by considering thermophoresis and Brownian motion effects on MHD mixed convective flow across a vertical stretching surface. The influence of velocity slip, thermal slip, and concentration slip along with an inclined external magnetic field is also considered. The governing coupled non-linear partial differential equations are transformed into ordinary differential equations using similarity transformation. The resulting system of non-linear ODEs is solved by the Newton Raphson shooting technique using the RK-4 algorithm. The impact of various physical parameters discovered in the problem viz. endothermic/exothermic reaction variable, thermophoresis parameter, activation energy parameter, Brownian motion parameter, chemical reaction parameter have been analyzed on velocity profile, temperature profile, and concentration profile. The effects of these parameters on skin-friction coefficient, Nusselt number, and Sherwood number are displayed in tabular form as well as surface plots. The impact of various physical parameters that appeared in the entropy generation is shown using surface and contour plots. The numerical findings are in good agreement with the previously published results. It is observed that an increment in thermophoresis and Brownian motion parameters results in a declination of entropy profiles, whereas an increment in Bejan number profiles is observed. A small region near the surface exhibits an inclination in concentration profiles with an increase in the order of the chemical reaction. In contrast, the opposite effect is analyzed near the boundary layer. Also, the contour and surface plots are displayed to portray real-world applications in industrial and technical processes and the physical depiction of flow characteristics that arise in the current study.

The mixed convection flows with simultaneous heat mass transfer involving Arrhenius activation energy with chemical reactions have been studied in recent years due to their vast applications. Most industrial end products’ quality is determined by cooling rates and chemical reactions, either the reaction rate or the type of chemical reaction. The present model includes activation energy, which most researchers have not included in earlier studies. Activation energy is widely considered while studying various physical phenomena, such as oil storage and engineering. A few theoretical publications on the role of activation energy in fluid dynamics are available. In 1889, Arrhenius made a groundbreaking attempt to present the concept of activation energy. The activation energy is the smallest amount of energy required by the reactants for a chemical reaction to occur. This phenomenon is used extensively in nuclear reactors, oil and water emulsion mechanics, and food processing. Menzinger and Wolfgang1 explained the detailed meaning of Arrhenius energy. Bestman2 was the first to develop and study this phenomenon in boundary layer transport. Makinde et al.3 investigated unsteady natural convection flow numerically under the impact of nth-order chemical reaction and activation energy. In the presence of thermal radiation, Maleque4 analyzed the effects of endothermic/exothermic chemical reactions having Arrhenius activation energy on MHD free convection flow. Shafique et al.5 used a numerical technique to quantitatively report a rotating viscoelastic flow with activation energy. Tripathi et al.6,7 discussed the influence of chemical reaction on blood flow considering variable viscosity model. Dhalmini et al.8 addressed the entropy generation and activation energy in viscous nanofluids containing higher-order chemically reactive species. Ullah9 investigated the activation energy associated with exothermic/endothermic reactions on magnetized nanomaterials flowing through a Darcy–Forchheimer porous medium. Dawar et al.10 studied mixed convective MHD flow of magnetite ferroparticles (Fe\(_3\)O\(_4\)) with blood as a base fluid past a non-isothermal vertical flat plate. Dawar et al.11 performed a mixed convective MHD flow of a water-based Al\(_2\)O\(_3\) nanofluid towards the stagnation region of an angularly revolving sphere.

In the intervening years, a considerable amount of attention has been received by combining heat and mass transfer problems with chemical reactions. The temperature differences lead to heat transfer. Therefore an enormous variety of heat transfer equipment was created to address these differences, such as boilers, condensers, radiators, furnaces, refrigerators, solar collectors, compact heat exchangers, and many more. The impact of a drop of dye in water is an example of mass transfer. Mass transport has many industrial applications. Air and water pollution processes are also diffusion controlled. Heat and mass transport co-occur in various processes, such as flow in a desert cooler, evaporation from the top of any water body, and energy transfer in a wet cooling tower. Kandasamy et al.12 researched the impact of heat and mass transportation with the thermal stratification effects on MHD flow through a stretching surface. Sharma et al.13 discussed the transfer of mass and heat on 3D flow through porous medium. Rajeswari et al.14 studied the influence of suction on heat and mass transport through a porous vertical surface. Ahmad and Khan15 studied heat and mass transmission with viscous dissipation over a moving wedge. The influence of heat and mass transport on nanofluid flow over a stationary/moving vertical plate embedded in a porous media was investigated by Madhura et al.16. Sharma and Kumawat17 studied heat and mass transportation considering Ohmic heating and variable viscosity effects via a stretching surface. The heat transfer analysis of a water-based hybrid nanofluid, including ferrous and graphene oxide nanoparticles, over a flat plate was studied by Dawar et al.18 using magnetohydrodynamics.

MHD is important to cover the liquid's mechanical properties in fluid dynamics and deals with the cooperation among the electrically conductive and electromagnetic fluids. The current may be generated whenever the conductive liquid particles move under the concurrent action of the electric and magnetic fields. The interaction with the magnetic field contributes to a body force on the liquid. MHD flow takes place in the sun as well as in the interior of the earth. Many new devices in the laboratory make full use of the MHD interaction, such as propulsion units and power generators, including interactions between fluid and electromagnetic fields, such as dynamics of electron beams and traveling wave tubes. As many critical engineering and industrial applications exist, combining mass and heat transport by MHD flow with chemical reaction inclusion has also received ample attention. Mansour et al.19 considered the presumptions of Soret and Dufour effects on the transportation of heat and mass through MHD free convective flow. Samad et al.20 explored heat and mass transport free convective flow with heat generation considering the magnetic field. Rajesh21 studied the impact on MHD free flow of thin gray fluid. Jafar et al.22 studied the transportation of heat and MHD flow over shrinking or stretching sheets. Under the presumption of varying permeability, Sharma et al.23,24 reported the chemical reaction's influence on the micropolar fluid flow that replicates the microscopic effects owing to the local behavior and micro-motion of the liquid particles. Waqas et al.25 examined the micro-level impact of micropolar liquid on MHD flow through a non-linear sheet under the convective condition. Srinivasulu et al.26 investigates the effect of an aligned magnetic field with convective boundary conditions using numerical methods through a stretching surface. The impact of an induced magnetic field on the flow of Maxwell nanofluid toward a vertically permeable and stretchy sheet was explored by Walelign et al.27. Dawar et al.28 investigated the effect of an inclined magnetic field and inter-particle spacing on the two-dimensional flow of an electrically conducting water-based copper nanofluid over a stretching surface using a porous medium. Gandhi and Sharma29 studied MHD two-dimensional pulsatile blood flow through a vertical artery with irregular stenosis. A comparison of the magnetohydrodynamic flow of copper and copper oxide hybrid nanofluids based on water and kerosene oil across a bi-directional stretching surface was conducted by Dawar et al.30.

Thermophoresis is a temperature gradient-induced migration of suspended particles from higher to lower locations. Nuclear power plants, micro contamination, and aerosol collector are examples of their use. The suspended particles or water droplets in the air are known as aerosols. The thermophoretic force is the force that causes these aerosol particles to migrate due to a temperature gradient. Brownian motion is the ‘indecisive’ random movement of particles suspended in a fluid due to collisions with the fluid's fast-moving molecules. Thermophoresis and Brownian motion are essential in heat and mass transmission in fluids. Hayat et al.31 examined MHD squeezing flow across a penetrable stretching surface in the presence of nanoparticles under the effect of Brownian movement. Sulochana et al.32 investigated the impact of stagnated point flow of a Carreau fluid past a shrinking/stretching plate with Brownian movement and MHD. Reddy et al.33 studied flow through a stretching geometry considering the effects of varying drag force and thermophoresis. Shah et al.34 studied thermophoresis's variation on 3D nanofluid flow using a rotating system between parallel plates. Considering the non-Newtonian Prandtl fluid approach, Soomro et al.35 examined the effects of Brownian motion and thermophoresis via a vertical stretching surface on the MHD stagnation-point flow of nanofluid. The investigation of heat and mass transport in MHD Williamson nanofluid flow over a vertical Riga plate with nonlinear thermal radiation was carried out by Rooman et al.36. Dawar et al.37 addressed thermophoresis and Brownian motion effects considering Cu and CuO nanoparticles of different shapes.

The irreversible processes of viscous dissipation and Joule heating show how electrical and kinetic energy is transformed into thermal energy. Viscous dissipation is the work done via the fluid on nearby layers due to shear forces. In contrast, Joule heating is a mechanism in which conduction electrons transfer into the atoms of a conductor due to a collision procedure. Beg et al.38 investigated Joule heating effects on MHD Hartmann-Couette flow with Hall current. To introduce the outcome of viscous dissipation, Sahoo39 presented the second-grade MHD fluid flow past a transversally stretching surface. Many researchers40,41,42 discussed the influence of Joule heating with various additional conditions such as MHD flow, solar radiations, convective boundary conditions, and partial slip. Hsaio43,44 discussed the impact of viscous dissipation on coupled electrical MHD heat transfer and micropolar nanofluid flow. Gayatri et al.45 explored the Joule heating effects and viscous dissipation across a stretching sheet with varying thickness, utilizing slip parameters. Seethamahalskshmi et al.46 investigated the MHD mixed convective flow through a semi-infinite vertical plate under the influence of Joule heating and viscous dissipation using the two-term perturbation technique. Gandhi et al.47 analyzed the simultaneous effects of viscous dissipation and Joule heating through a stenosed artery considering the variable viscosity model.

Different forms of thermal systems are linked to the irreversibility process, which may be characterized using entropy generation, and are significant to viscous dissipation, magnetic fields, heat and mass transport, etc. Various studies employed the first law of thermodynamics to improve this irreversibility process, but the results were insufficient. Later, several researchers used the second law of thermodynamics to optimize these irreversibilities, concluding that the 2nd law is more efficient than the 1st. Rashidi et al.48 explored entropy generation in a rotating disc MHD fluid flow with variable properties. Dalir et al.49 used the Keller-box scheme to investigate the entropy formation for MHD heat and mass transfer flow of Jeffrey nanofluid across a stretched sheet. Baag et al.50 computed entropy generation by applying the second law of thermodynamics to MHD heat and mass transfer of an electrically conducting viscoelastic fluid past a stretching surface considering Darcy dissipation besides viscous and Joule dissipation. Bhatti et al.51 investigated entropy generation on the boundary layer flow with chemical reaction effects. Khan et al.52 performed entropy generation analysis in a mixed convective flow of nanomaterials with thermophoresis and Brownian movement considering the Buongiorno nanofluid model. Sohail et al.53 performed entropy generation computation for Casson fluid past a bi-directional stretched surface with heat and mass conveyance having a variable thermal conductivity. Hayat et al.54 studied irreversibility in the Darcy–Forchheimer flow of nanofluid by a curved stretching sheet with a spiral shape. Sharma et al.55 performed entropy analysis through a multi-stenosed artery in the prsence of hybrid nanoparticles (Au-Al\(_2\)O\(_3\)/Blood). Gandhi et al.56 conducted an entropy analysis of the MHD blood flow of hybrid nanoparticles of various shapes through an irregularly stenosed permeable walled artery under periodic body acceleration. Sharma et al.57 analyzed entropy generation effects on EMHD Jeffrey fluid flow over a vertical stretching surface.

The current physics of flow across a vertical stretching surface is expected to serve as the foundation for various medical science, engineering, and technology applications. The combined effect of physical characteristics may aid scientists in comprehending their findings. To the best of our knowledge, no effort has been yet made to perform the entropy generation minimization of higher-order endothermic/exothermic chemical reactions with activation energy on MHD mixed convective flow over a vertical stretching surface in the presence of thermophoresis and Brownian motion. Hence, the motivation from the above studies inspired us to perform this analysis. The following are some significant novel aspects included in this study: (1) to perform entropy generation minimization of higher-order endothermic/exothermic chemical reactions with activation energy, (2) to incorporate thermophoresis and Brownian motion effects along with the imposition of a time-dependent inclined magnetic field, (3) to add velocity, thermal, and concentration slips along with injection/suction effects. The current problem could help researchers use this method to solidify liquid metal from the mushy zone, build a metallic layer around a thermonuclear fusion-fission hybrid reactor, and produce drug delivery systems and gene therapy. The present study initiative is organized into six sections, as follows:

The first section is an introduction, which describes the various physical quantities in this and other relevant studies.

The model's geometry and the flow's governing equations are included in the second part, namely Mathematical Formulation.

The third section comprises similarity transformation. Also, the given PDEs are converted to nonlinear coupled ODEs using these similarity transformations. This section introduces non-dimensional variables used in the present study to generate governing equation solutions.

The fourth section is the numerical solution, which explains the numerical procedure used to solve the nonlinear coupled ODEs in which RK-4, along with the Newton Raphson shooting technique, are employed.

The fifth section describes the minimization of entropy generation and Bejan number effects.

Finally, there's the section on Results and Graphical Analysis. The results are shown graphically in MATLAB, and the graphical results are then elaborated. The surface and contour plots are drawn to precisely analyze the flow parameter's behavior.

An unsteady, incompressible, laminar, viscous, electrically conducting MHD boundary layer flow across a stretching vertical sheet with viscous dissipation, thermophoresis, Brownian motion, Joule heating, and higher-order endothermic/exothermic chemical reaction is under consideration. The cartesian coordinate system is used with the \(x_1^*\) and the \(y_1^*\) axis. Within the fluid medium, the origin is considered fixed with ambient temperature \(T_\infty ^*\), and the surface is kept at uniform temperature \({\tilde{T}}_w\). The surface concentration is kept uniform at \({\tilde{C}}_w\) whereas the ambient fluid concentration is \(C_\infty ^*\). The stretching sheet velocity is

along \(x_1^*\) axis at time \(t_1^*=0\), where \({\tilde{p}}\) and \({\tilde{r}}\) are constants. Here, \({\tilde{p}}\) represents the initial stretching rate, while \(\frac{{\tilde{p}}}{(1-{\tilde{r}}t_1^*)}\) represents the effective stretching rate over time. An inclined magnetic field \({\tilde{B}}(t_1^*)\) with an acute angle \(\xi\) is applied along the \(x_1^*\)-direction. The magnetic Reynold's number is assumed to be very small \((Re \ll 1)\) in this study, so the induced magnetic field effect can be neglected. Figure 1 depicts the pictorial representation of the model. Based on the above assumptions and using the order of magnitude approach along with the usual Boussinesq's approximation for the boundary layer, the governing equations are58,59,60:

A pictorial representation of the model.

The boundary conditions subject to the flow are58,61:

and

\({\tilde{V}}_w^*\) is specified by

where, \({\tilde{V}}_w^* > 0\) demonstrates injection and \({\tilde{V}}_w^* < 0\) demonstrates suction.

In Eq. (5), \({\tilde{E}}={\tilde{E}}_0(1-{\tilde{r}}t_1^*)^{1/2}\), \({\tilde{F}}={\tilde{F}}_0(1-{\tilde{r}}t_1^*)^{1/2}\), \({\tilde{G}}={\tilde{G}}_0(1-{\tilde{r}}t_1^*)^{1/2}\) represent velocity, thermal and concentrations respectively with \(E_0\), \(F_0\), \(G_0\), being their initial values. Also, \(E_0=0\), \(F_0=0\), \(G_0=0\) corresponds to the no-slip boundary condition.

Further, it is assumed that

where \({\tilde{p}}>0,{\tilde{q}} \ge 0,{\tilde{r}} \ge 0,{\tilde{s}} \ge 0\) are constants and \({\tilde{r}}t_1^*<1\).

We consider \({\tilde{B}}={\tilde{B}}_0^*(1-{\tilde{r}}t_1^*)^{-1/2}~and~\Gamma (t_1^*)=\Gamma _0(1-{\tilde{r}}t_1^*)^{-1}\) where \({\tilde{B}}_0^*\) represents the magnetic field strength at \(t_1^*=0\) and \(\Gamma _0\) is a constant.

The following similarity transformation is used to obtain the solution to the governing equations:

where \(\psi ^*\) is the stream function that satisfies the continuity equation (1).

The velocity components are: \({\tilde{u}}_1^*=\frac{\partial \psi ^*}{\partial y_1^*}\) and \({\tilde{v}}_1^*=-\frac{\partial \psi ^*}{\partial x_1^*}\).

On calculating we have \({\tilde{u}}_1^*={\tilde{U}}_w^* {\tilde{\chi }}'(\eta )\) and \({\tilde{v}}_1^*=-\sqrt{\frac{{\tilde{p}}{\tilde{\nu }}^*}{(1-{\tilde{r}}t_1^*)}}{\tilde{\chi }}\).

The substitution of the similarity transformation introduced in Eq. (9) to the governing Eqs. (2)–(4), the following set of ODEs is obtained:

The non-dimensional parameters used in the above equations are mentioned in Table 1.

The dimensional boundary conditions are reduced to the following non-dimensional boundary conditions:

In Eq. (13), the injection is represented by \(S \le 0\) whereas suction is represented by \(S \ge 0\). Also,

.

Flow chart illustrating solution methodology.

Equations (10)–(12) are the higher-order nonlinear coupled ordinary differential equation. To solve the system of coupled non-linear ODEs (10)–(12) with boundary conditions (13), the Newton Raphson shooting technique is employed in combination with the RK-4 algorithm. The boundary value problem of the physical model is first turned into an initial value problem. The system of Eqs. (10)–(12) comprises three differential equations out of which one is third-order, and the other two are second-order equations. Therefore, it cannot be solved until seven initial conditions are specified. However, initially, there are only four conditions defined as given in Eq. (13). For obtaining the solution, the conditions \({{\tilde{\chi }}}'(\infty )=0\), \({{\tilde{\zeta }}}(\infty )=0\), and \({{\tilde{\phi }}}(\infty )=0\) are replaced by \({{\tilde{\chi }}}''(0)=l_1\), \({{\tilde{\zeta }}}'(0)=l_2\), and \({{\tilde{\phi }}}'(0)=l_3\) (initial guesses). Furthermore, the \(\eta _{\infty }\) should have a finite upper bound. The solution is then calculated using the RK fourth-order approach. Finally, the computed solution will converge if the residuals are less than the error tolerance (\(10^{-6}\)). The Newton's method is used to modify the original guesses if the computed solution fails to meet the convergence condition. The solution procedure is depicted with the help of a flowchart in Fig. 2.

The governing Eqs. (10)–(12) are the nonlinear coupled ordinary differential equation. For solving, these are converted into a system of first order differential equations. Let

Therefore, on introducing these new variables the Eqs. (10)–(12) transformed into the following system-

and the boundary conditions transformed from the equation (13) are-

The entropy of a system is a broad attribute that changes as mass and energy are exchanged. The overall entropy of a system made up of numerous processes is equal to the sum of the entropies produced by each process. The entropy generation rate due to the exchange of momentum, energy, and mass explains the irreversibility of MHD mixed convective flow over a vertical stretching surface with Joule effect, thermophoresis, Brownian motion, viscous dissipation and higher-order endothermic/exothermic chemical with activation energy. The volumetric entropy generation is defined as51,52:

where \({\tilde{F}}^*\) represents viscous dissipation.

Therefore, the entropy generation rate due to the exchange of momentum, energy, and mass is given as:

On applying the similarity transformation given in Eq. (9),

The dimensionless entropy generation number can be described as the ratio of the characteristics entropy generation rate to the actual entropy generation rate.

Be depicts the ratio of irreversibility owing to heat transfer and total irreversibility owing to heat transfer and fluid friction. In mathematical form, it is described as:

The present study discusses the entropy generation effects of higher-order endothermic/exothermic chemical reactions on MHD mixed convective flow across a vertical stretching surface with Joule heating, thermophoresis, Brownian motion, and viscous dissipation. The influence of velocity, thermal, and concentration slip is also investigated. The impact of the flow parameters identified in the problem such as Hartmann number (M), Grashof number (Gr), solutal grashof number (Gc), velocity slip (\(S_v\)), suction parameter (S), inclination parameter (\(\xi\)), Prandtl number (Pr), thermophoresis parameter (\(N_t\)), Eckert number (Ec), endothermic/exothermic reaction parameter (\({{\tilde{\lambda }}}_1\)), Brownian motion parameter (\(N_b\)), activation energy parameter (\({\tilde{E}}^*\)), thermal slip (\(S_t\)), order of chemical reaction (N), Schmidt number (Sc), concentration slip (\(S_c\)), and chemical reaction parameter (\({{\tilde{\sigma }}}_1\)) is explored on the entropy generation (\(N_s\)), Bejan number (Be), velocity profile (\({{{\tilde{\chi }}}}'(\eta )\)), temperature profile (\({{{\tilde{\zeta }}}}(\eta )\)), and the concentration profile (\({{{\tilde{\phi }}}}(\eta )\)) to gain a physical understanding of the problem. For numerical results, some default values for the parameters are depicted in Table 2. These values are considered default unless mentioned in the relevant graphs.

The effect of inclination of the surface, permeability, heat source, and radiation parameter in58 and the influence of thermophoresis, Brownian motion, and higher-order endothermic/exothermic reaction with activation energy in the present work is neglected to validate the current results with Reddy et al.58. Figure 3 depict the present work's velocity and temperature profile with the study done by Reddy et al.58. Further, the results are compared with the available results of Sharma and Gandhi62. Table 3 shows the comparison of \({{\tilde{\chi }}}''(0),-{{\tilde{\zeta }}}'(0),-{{\tilde{\phi }}}'(0)\) for62 and the present work. The current results are in good agreement (under specific limiting conditions), clearly demonstrating the accuracy of the computed results.

Comparative analysis of (a) velocity profile \({{\tilde{\chi }}}'(\eta )\) for M = 3, (b) temperature profile \({{\tilde{\zeta }}}(\eta )\) for Pr = 7.

A contour plot is a 2-D representation of the surface in which similar-responding points are linked to generate contour lines with constant responses. The contour plots help determine the desired response values and operating circumstances. Figure 4 represents the contour plots for entropy generation (\(N_s\)) and Bejan number (Be). The influence of Hartmann number (M) on entropy generation (\(N_s\)) and Bejan number (Be) is shown via contours in Fig. 4a,d. It is noticed that as M increases, entropy drops while Be increases. Figure 4b,e illustrate the contour plots depicting the influence of thermophoresis parameter (\(N_t\)) on entropy generation (\(N_s\)) and Bejan number (Be) respectively. With an increase in \(N_t\) values, entropy decreases. However, an increment in Bejan number is investigated. Figure 4c,f highlight the contours depicting the influence of Brownian motion parameter (\(N_b\)) on entropy generation (\(N_s\)) and Bejan number (Be). As the \(N_b\) values rise, so does Be. On the other hand, the \(N_s\) is declining.

(a) \(N_s\) versus M, (b) \(N_s\) versus \(N_t\), (c) \(N_s\) versus \(N_b\), (d) Be versus M, (e) Be versus \(N_t\), and (f) Be versus \(N_b\).

The effect of different flow parameters namely Hartmann number (M), Grashof number (Gr), solutal Grashof number (Gc), velocity slip (\(S_v\)), suction parameter (S), inclination parameter (\(\xi\)) on dimensionless velocity \(\chi '(\eta )\) is shown in Fig. 5. The velocity profile \(\chi '(\eta )\) for various values of M is shown in Fig. 5a. The velocity profiles show deterioration as M increases from 3 to 5. The Lorentz force, thus generated, opposes the flow and retards the fluid's velocity with the increment in values of M. The velocity profile \(\chi '(\eta )\) for Gr and Gc are shown in Fig. 5b,c. The Gr is the proportion of buoyant and viscous forces in a fluid layer. Since the viscous forces become less dominating as the value of Gr enhances, the resistance to flow diminishes, and fluid flow velocity rises. There is an abrupt increase in velocity near the wall and afterward descends uniformly towards zero. The velocity profile for the solutal Grashof number Gc also increases with increasing Gc values. Figure 5d shows the profiles for non-dimensional velocity for \(S_v\). The velocity profile deteriorates with increasing values of \(S_v\). The momentum boundary layer rises as we increase the values of \(S_v\), yet surface velocity exhibits a declining trend. This mechanism occurs because the fluid's velocity is reduced due to the stretching velocity partially transferring the disturbance caused by frictional retardation between the surface and the fluid's particles. Thus, the velocity profile drops. The non-dimensional velocity profile \(\chi '(\eta )\) for S is depicted by Fig. 5e. With rising values of S, the velocity profiles decline slightly. The heated fluid is propelled towards the wall due to the considerable influence of viscosity, where buoyant forces can interfere to retard the fluid. Figure 5f expresses that increment in \(\xi\) results in the decrement in non-dimensional velocity because the magnetic field strength enhances with the increment in the aligned angle. An opposite flow to force is generated due to this enhanced magnetic field, known as the Lorentz force. This generated force thus becomes a barrier in the fluid's path. Hence, a decrease in \(\chi '(\eta )\) is analyzed.

Non-dimensional velocity profiles for (a) Hartmann number (M), (b) Grashof number (Gr), (c) solutal Grashof number (Gc), (d) velocity slip (\(S_v\)), (e) suction parameter (S), and (f) inclination parameter (\(\xi\)).

The different parameters’ impact such as Prandtl number (Pr), Eckert number (Ec), thermophoresis parameter (\(N_t\)), Brownian motion parameter (\(N_b\)), endothermic/exothermic reaction parameter (\({{\tilde{\lambda }}}_1\)), activation energy parameter (\({\tilde{E}}^*\)), thermal slip (\(S_t\)), and order of chemical reaction (N) on temperature profile \({{{\tilde{\zeta }}}}(\eta )\) is highlighted in Fig. 6. Figure 6a highlights the influence of Pr on \({{{\tilde{\zeta }}}}(\eta )\). The temperature profile diminishes as Pr increases because Pr regulates the relative thickness of the thermal and momentum boundary layers. Because thermal conductivity increases as Pr decreases, heat diffusion from the heated surface happens faster for small Pr values than for large Pr values. The non-dimensional temperature profiles \({{{\tilde{\zeta }}}}(\eta )\) for different Ec values are shown in Fig. 6b. The temperature profiles \({{\tilde{\zeta }}}(\eta )\) rise with an increment in Ec as the internal energy increases. Ec is the relationship between the flow's kinetic energy and the heat transfer's enthalpy driving force. Figure 6c,d demonstrate the impact of \(N_t\) and \(N_b\) on dimensionless temperature profile \({{{\tilde{\zeta }}}}(\eta )\). The purpose of Fig. 6c is to show how \(N_t\) affects the thermal field. Higher values of the thermophoretic parameter (\(N_t\)) are shown to result in elevated temperature profiles in the boundary layer region. This results from the fact that particles close to a hot surface produce a thermophoretic force that aids particle disintegration away from the fluid regime, which increases temperature boundary layer thicknesses. Figure 6d shows that a temperature rise is perceived for \(N_b\). The random movement of suspended particles in the base fluid, known as Brownian motion, is more impacted by the fluid's quickly moving atoms or molecules. It is important to note that Brownian motion is connected to particle size and that these particles frequently take the shape of aggregates or agglomerates. Brownian motion is very weak for massive particles, and the parameter (\(N_b\)) will have the lowest values. As the values of the Brownian motion parameter (\(N_b\)) increase, the temperature profiles in the boundary layer region increase. Figure 6e portrays the influence of \({{\tilde{\lambda }}}_1\) on \({{{\tilde{\zeta }}}}(\eta )\). An increment in the temperature profiles is observed with increasing values of \({{\tilde{\lambda }}}_1\), and the results obtained are consistent with that of9. Figure 6f shows the influence of \({\tilde{E}}^*\) on \({{{\tilde{\zeta }}}}(\eta )\). The temperature profiles improve for higher elevations of \({\tilde{E}}^*\) due to the generative reaction. The effect of \(S_t\) on \({{{\tilde{\zeta }}}}(\eta )\) is depicted in Fig. 6g. There is a declination in temperature profiles with increasing values of \(S_t\) since the thickness of the thermal boundary layer reduces. As a result, non-dimensional temperature profiles are diminished. Figure 6h shows the impact of N on \({{{\tilde{\zeta }}}}(\eta )\). A reduction in \({{{\tilde{\zeta }}}}(\eta )\) is seen with rising values of N. The obtained results show a good agreement with that of 8.

Non-dimensional temperature profiles for (a) Prandtl number (Pr), (b) Eckert number (Ec), (c) thermophoresis parameter (\(N_t\)), (d) Brownian motion parameter (\(N_b\)), (e) endothermic/exothermic reaction parameter (\({{\tilde{\lambda }}}_1\)), (f) activation energy parameter (\({\tilde{E}}^*\)), (g) thermal slip (\(S_t\)), and (h) order of chemical reaction (N).

The effect of several parameters on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\), including Schmidt number (Sc), concentration slip (\(S_c\)), thermophoresis parameter (\(N_t\)), Brownian motion parameter (\(N_b\)), endothermic/exothermic reaction parameter (\({{\tilde{\lambda }}}_1\)), chemical reaction parameter (\({{\tilde{\sigma }}}_1\)), activation energy parameter (\({\tilde{E}}^*\)), and order of chemical reaction (N) is depicted in Fig. 7. Figure 7a shows the impact of Sc on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\). The concentration profiles show declination with increment in the values of Sc for both N = 1 and N = 2. Sc refers to fluid flows simultaneously undergoing both momentum and mass diffusion convection processes. When the Schmidt number is large enough, momentum diffusion takes precedence over mass diffusion, and when it is small, mass diffusion takes precedence over momentum diffusion. The mass transfer boundary layer thickness is lesser than the hydrodynamic boundary layer thickness with increasing values of Schmidt number. Therefore, a decrease in non-dimensional concentration profiles \({{{\tilde{\phi }}}}(\eta )\) is investigated. Figure 7b represents the influence of \(S_c\) on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\). The concentration profile decreases as the value of \(S_c\) increases, which is consistent with the results of 58. It is because slip essentially slows down fluid motion, which ultimately manifests as a reduction in net molecular mobility. Therefore, decreased molecular mobility causes mass fraction fields to decrease. The concentration slip parameter can likely control the mass transport phenomenon as the velocity slip parameter, and the thermal slip parameter can regulate the momentum and temperature inside the flow. In conclusion, the concentration boundary layer can be controlled up to the appropriate levels by modifying the concentration slip parameters. Figure 7c,d depicts the influence of \(N_t\) and \(N_b\) on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\). The concentration in a narrow region near the surface is reduced when the thermophoretic parameter is increased. This phenomenon shows particles being ‘pushed’ out of the heated boundary layer into the colder free stream zone. This behavior changes when we proceed away from the surface into the free stream, and rising \(N_t\) raises the concentration. On the other hand, the solute boundary layer diminishes as there is an increment in \(N_b\). Particle mobility is aided by increasing \(N_b\), which causes the boundary layer to warm, causing particles to move away from the surfaces within the inactive fluid. As a result, solute particle deposition away from the surface rises, causing concentration profiles to fall. The effect of \({{\tilde{\lambda }}}_1\) and \({{\tilde{\sigma }}}_1\) on non-dimensional concentration \({{\tilde{\phi }}}(\eta )\) is demonstrated by Fig. 7e,f. It is found that the non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\) enhances with increase in the value of \({{\tilde{\lambda }}}_1\). Increasing \({{\tilde{\sigma }}}_1\), on the other hand, causes the mass transfer boundary layer to thicken. It has been discovered that increasing the reaction rate constant results in an improvement in the factor \({{\tilde{\sigma }}}_1(1+\delta ^*{{\tilde{\zeta }}})^m exp\bigg (\frac{-{\tilde{E}}^*}{1+\delta ^*{{\tilde{\zeta }}}}\bigg )\). As a result, a destructive chemical reaction occurs, and the concentration profiles \({{{\tilde{\phi }}}}(\eta )\) diminishes. Figure 7g highlights the impact of \({\tilde{E}}^*\) on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\). The definition of activation energy is that it is the least amount of energy necessary to initiate a reaction. It is discovered that higher activation energy causes a decrease in the reaction rate constant, which ultimately causes the chemical reaction to slow down. Additionally, the concentration profiles show enhancement. Figure 7h portrays the impact of N on non-dimensional concentration \({{{\tilde{\phi }}}}(\eta )\). There is an inclination in \({{{\tilde{\phi }}}}(\eta )\) with increasing values of N in a narrow region near the surface. However, it changes its behavior in the boundary layer and shows the opposite trend.

Non-dimensional concentration profiles for (a) Schmidt number (Sc), (b) concentration slip (\(S_c\)), (c) thermophoresis parameter (\(N_t\)), (d) Brownian motion parameter (\(N_b\)), (e) endothermic/exothermic reaction parameter (\({{\tilde{\lambda }}}_1\)), (f) chemical reaction parameter (\({{\tilde{\sigma }}}_1\)), (g) activation energy parameter (\({\tilde{E}}^*\)), and (h) order of chemical reaction (N).

Flow resistance directly correlates with the fluid flow rate and determines the physiological features of the flow. One of the physical quantities that influence the flow of fluid is shear-stress. The shear-stress expression in mathematical form is

Therefore,

To estimate and understand the heat transfer, the \(Nu_x^*\), which is the ratio of convective to the conductive thermal transfer in the fluid across the boundary is calculated whose general expression is

where \({\tilde{q}}_w^*=-{\tilde{k}}^*\bigg (\frac{\partial {\tilde{T}}_f}{\partial y_1^*}\bigg )_{y_1^*=0}\).

The ratio of the mass transfer due to convection to the diffusive mass rate (called Sherwood number) is

where \({\tilde{m}}_w^*=-{{\tilde{\rho }}}^* {\tilde{D}}\bigg (\frac{\partial {\tilde{C}}_f}{\partial y_1^*}\bigg )_{y_1^*=0}\).

Expressing in dimensionless form, we have

Tables 4 and 5 describe the values of \({{\tilde{\chi }}}''(0)\), \(-{{\tilde{\zeta }}}'(0)\), \(-{{\tilde{\phi }}}'(0)\) for contrasting values of velocity slip (\(S_v\)), thermal slip (\(S_t\)) and concentration slip (\(S_c\)) for A=0 and A=0.5, respectively. A decrement in skin-friction coefficient (\(C_f^*\)) is observed with increment in thermal slip (\(S_t\)) and concentration slip (\(S_c\)) whereas it increases with velocity slip (\(S_v\)). With an increment in thermal slip (\(S_t\)), the Nusselt number (\(Nu_x^{*}\)) declines whereas with increment in velocity slip (\(S_v\)) and concentration slip (\(S_c\)) it enhances. As the value of velocity slip (\(S_v\)), thermal slip (\(S_t\)) and concentration slip (\(S_c\)) is increased, a decrement in Sherwood number (\(Sh_x^*\)) is analyzed. The Nusselt number (\(Nu_x^*\)), and the Sherwood number (\(Sh_x^*\)) increase for velocity slip (\(S_v\)), thermal slip (\(S_t\)) and concentration slip (\(S_c\)) as the value of unsteadiness parameter A varies from 0 to 0.5. In contrast, the skin-friction coefficient (\(C_f^*\)) decreases for for velocity slip (\(S_v\)), thermal slip (\(S_t\)) and concentration slip (\(S_c\)) as the value of unsteadiness parameter A varies from 0 to 0.5.

Surface plots are three-dimensional data visualizations. Surface plots demonstrate a functional relationship between a dependent variable and two independent variables rather than individual data points. Figure 8 represents surface plots for skin-friction coefficient (\(C_f^*\)), Nusselt number (\(Nu_x^*\)), and Sherwood number (\(Sh_x^*\)) for different flow parameters. The influence of Hartmann number (M) on skin-friction coefficient (\(C_f^*\)) is depicted in Fig. 8a. Skin-friction coefficient (\(C_f^*\)) decreases with increment in Hartmann number (M) as well unsteadiness parameter (A). Figure 8b shows the effect of Prandtl number (Pr) on Nusselt number (\(Nu_x^*\)). It is analyzed that there is an increase in the value of Nusselt number (\(Nu_x^*\)) with increment in both Prandtl number (Pr) as well as unsteadiness parameter (A). Figure 8c,d highlights the influence of thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters on Nusselt number (\(Nu_x^*\)). It can be viewed that there is an increment in Nusselt number (\(Nu_x^*\)) with increasing values of both thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters. The impact of Schmidt number (Sc) on Sherwood number (\(Sh_x^*\)) is illustrated in Fig. 8e. The Sherwood number (\(Sh_x^*\)) increases with increment in both Schmidt number (Sc) and unsteadiness parameter (A). Figure 8f portrays the effect of chemical reaction parameter (\({{\tilde{\sigma }}}_1\)) on Sherwood number (\(Sh_x^*\)). There is a decrement in Sherwood number (\(Sh_x^*\)) with increasing values of chemical reaction parameter (\({{\tilde{\sigma }}}_1\)) since raising the chemical reaction rate increases the thickness of mass transfer boundary layer. Figure 8g,h illustrates the impact of thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters on Sherwood number (\(Sh_x^*\)). There is slight increment in Sherwood number (\(Sh_x^*\)) values with increasing values of thermophoresis parameter (\(N_t\)) whereas a declination in Sherwood number (\(Sh_x^*\)) values is observed with Brownian motion parameter (\(N_b\)).

(a) \(C_f^*\) versus M and A, (b) \(Nu_x^*\) versus Pr and A, (c) \(Nu_x^*\) versus \(N_t\) and A, (d) \(Nu_x^*\) versus \(N_b\) and A, (e) \(Sh_x^*\) versus Sc and A, (f) \(Sh_x^*\) versus \({{\tilde{\sigma }}}_1\) and A, (g) \(Sh_x^*\) versus \(N_t\) and A, and (h) \(Sh_x^*\) versus \(N_b\) and A.

The present study performs entropy generation minimization on mixed convective flow across a vertical stretching sheet with an inclined magnetic field, thermophoresis, Brownian motion, viscous dissipation, higher-order endothermic/exothermic chemical reaction with activation energy, and Joule heating. The RK-4 method, in combination with the shooting method, has been used to solve the resulting set of ODEs. The numerically obtained results have been compared to those published in the literature, and the results were found to be in good agreement. The following are some of the most important findings of the research:

A declination in entropy profiles is observed with an increase in thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters, while Bejan number profiles show an increment.

An increment in inclination parameter (\(\xi\)) and velocity slip (\(S_v\)) exhibits declination in the velocity profiles.

The dimensionless temperature profile declines with an enhancement in the values of Prandtl number (Pr), activation energy parameter (\({\tilde{E}}^*\)) and order of chemical reaction (N), whereas an increment is observed with Eckert number (Ec), thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters.

There is reduction in dimensionless concentration profiles with increasing values of Schmidt number (Sc), chemical reaction parameter (\({{\tilde{\sigma }}}_1\)) and concentration slip (\(S_c\)).

Nusselt number (\(Nu_x^*\)) enhances with increment in Prandtl number (Pr), thermophoresis (\(N_t\)) and Brownian motion (\(N_b\)) parameters.

A decrement in the skin-friction coefficient (\(C_f^*\)) is observed with increment in thermal slip (\(S_t\)) and concentration slip (\(S_c\)) whereas it increases with increment in velocity slip (\(S_v\)).

Sherwood number (\(Sh_x^*\)) decreases with an increment in the values of thermal slip (\(S_t\)), concentration slip (\(S_c\)), and velocity slip (\(S_v\)), respectively.

Entropy generation minimization addressed in the present problem is helpful in several sectors of mainstream thermal engineering and science: cryogenics, heat transfer, education, storage systems, solar power plants, nuclear and fossil power plants, and refrigerators. Further, the applications for heat and mass transfer of the boundary layer flow across a stretching sheet are numerous and diverse, including manufacturing artificial film and fibers and some uses for diluted polymer solutions in the polymer processing sector. Multiple industries, including food processing, water mechanics, oil storage and geothermal energy production, base liquid mechanics, and oil emulsification, use the mass transfer process along with endothermic/exothermic chemical reactions, activation energy, and other related phenomena. The results of this problem can be used in various systems susceptible to considerable fluctuations in gravitational force, heat exchanger designs, wire and glass fiber drafting, and nuclear engineering concerning reactor cooling.

The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request.

Unsteadiness parameter

Uniform magnetic field

Brinkmann number

Concentration of the corresponding fluid

Coefficient of skin-friction

Specific heat at constant pressure

Concentration at the wall

Ambient concentration of the fluid

Brownian diffusion coefficient

Thermophoretic diffusion coefficient

Slip corresponding to velocity

Dimensional activation energy

Dimensionless activation energy

Eckert number

Slip corresponding to temperature

Acceleration by virtue of gravity

Slip corresponding to concentration

Grashof number

Solutal Grashof number

Thermal conductivity

Boltzmann constant

Fitted rate constant

Mass flux

Hartmann number

Order of chemical reaction

Thermophoretic motion parameter

Brownian motion parameter

Local Nusselt number

Prandtl number

Heat flux at the surface

Universal gas constant

Reynolds number

Suction parameter

Schmidt number

Local Sherwood number

Dimensionless concentration slip

Dimensionless thermal slip

Dimensionless velocity slip

Temperature of the corresponding fluid

Temperature at the wall

Ambient temperature of the fluid

Velocity component in \(x_1^*\)-direction

Stretching velocity of the sheet

Velocity component in \(y_1^*\)-direction

Velocity of suction or injection at the wall of the surface

Coefficient of thermal expansion

Coefficient of concentration expansion

Endothermic/exothermic coefficient

Dimensionless velocity

Temperature ratio parameter

Similarity variable

Endothermic/exothermic reaction variable

Fluid's dynamic viscosity

Fluid's kinematic viscosity

Temperature difference parameter

Dimensionless concentration

Stream function

Concentration difference parameter

Fluid's density

Electrical conductivity

Chemical reaction parameter

Heat capacity ratio

Shear stress at the wall

Inclination of magnetic field

Dimensionless temperature

Constant parameter

Menzinger, M. & Wolfgang, R. The meaning and use of the Arrhenius activation energy. Angew. Chem. Int. Ed. Engl. 8(6), 438–444 (1969).

Article CAS Google Scholar

Bestman, A. Natural convection boundary layer with suction and mass transfer in a porous medium. Int. J. Energy Res. 14(4), 389–396 (1990).

Article CAS Google Scholar

Makinde, O. D., Olanrewaju, P. O. & Charles, W. M. Unsteady convection with chemical reaction and radiative heat transfer past a flat porous plate moving through a binary mixture. Afr. Mat. 22(1), 65–78 (2011).

Article MathSciNet MATH Google Scholar

Maleque, K. et al. Effects of exothermic/endothermic chemical reactions with Arrhenius activation energy on MHD free convection and mass transfer flow in presence of thermal radiation. J. Thermodyn. 20, 13 (2013).

Google Scholar

Shafique, Z., Mustafa, M. & Mushtaq, A. Boundary layer flow of Maxwell fluid in rotating frame with binary chemical reaction and activation energy. Res. Phys. 6, 627–633 (2016).

Google Scholar

Tripathi, B. & Sharma, B. Effect of variable viscosity on MHD inclined arterial blood flow with chemical reaction. Int. J. Appl. Mech. Eng. 23(3), 767–785 (2018).

Article CAS Google Scholar

Tripathi, B. & Sharma, B. K. Influence of heat and mass transfer on two-phase blood flow with Joule heating and variable viscosity in the presence of variable magnetic field. Int. J. Comput. Methods 20, 1850139 (2018).

MathSciNet MATH Google Scholar

Dhlamini, M., Mondal, H., Sibanda, P. & Motsa, S. Activation energy and entropy generation in viscous nanofluid with higher order chemically reacting species. Int. J. Ambient Energy 20, 1–13 (2020).

Google Scholar

Ullah, I. Activation energy with exothermic/endothermic reaction and Coriolis force effects on magnetized nanomaterials flow through Darcy–Forchheimer porous space with variable features. Waves Random Complex Media 20, 1–14 (2022).

Google Scholar

Dawar, A. et al. Mixed convective flow of blood biofluids containing magnetite ferroparticles past a vertical flat plate: Shapes-based analysis. Wavesin and Complex Media 20, 1–25 (2022).

Google Scholar

Dawar, A. & Acharya, N. Unsteady mixed convective radiative nanofluid flow in the stagnation point region of a revolving sphere considering the influence of nanoparticles diameter and nanolayer. J. Indian Chem. Soc. 20, 100716 (2022).

Article Google Scholar

Kandasamy, R., Periasamy, K. & Prabhu, K. S. Chemical reaction, heat and mass transfer on MHD flow over a vertical stretching surface with heat source and thermal stratification effects. Int. J. Heat Mass Transf. 48(21–22), 4557–4561 (2005).

Article CAS MATH Google Scholar

Sharma, B., Chaudhary, R. & Sharma, P. Fluctuating mass transfer on three dimensional flow through a porous medium with variable permeability. Adv. Theor. Appl. Math. 2(3), 257–267 (2007).

Google Scholar

Rajeswari, R., Jothiram, B. & Nelson, V. Chemical reaction, heat and mass transfer on nonlinear MHD boundary layer flow through a vertical porous surface in the presence of suction. Appl. Math. Sci. 3(49–52), 2469–2480 (2009).

MathSciNet MATH Google Scholar

Ahmad, R. & Khan, W. A. Numerical study of heat and mass transfer MHD viscous flow over a moving wedge in the presence of viscous dissipation and heat source/sink with convective boundary condition. Heat Transfer-Asian Res. 43(1), 17–38 (2014).

Article Google Scholar

Madhura, K. et al. Impact of heat and mass transfer on mixed convective flow of nanofluid through porous medium. Int. J. Appl. Comput. Math. 3(1), 1361–1384 (2017).

Article MathSciNet Google Scholar

Sharma, B. K. & Kumawat, C. Impact of temperature dependent viscosity and thermal conductivity on MHD blood flow through a stretching surface with ohmic effect and chemical reaction. Nonlinear Eng. 10(1), 255–271 (2021).

Article ADS Google Scholar

Dawar, A., Islam, S., Alshehri, A., Bonyah, E. & Shah, Z. Heat transfer analysis of the MHD stagnation point flow of a non-Newtonian tangent hyperbolic hybrid nanofluid past a non-isothermal flat plate with thermal radiation effect. J. Nanomater. 2022, 2556 (2022).

Article Google Scholar

Mansour, M., El-Anssary, N. & Aly, A. Effects of chemical reaction and thermal stratification on MHD free convective heat and mass transfer over a vertical stretching surface embedded in a porous media considering soret and dufour numbers. Chem. Eng. J. 145(2), 340–345 (2008).

Article CAS Google Scholar

Samad, M. & Mohebujjaman, M. MHD heat and mass transfer free convection flow along a vertical stretching sheet in presence of magnetic field with heat generation. Res. J. Appl. Sci. Eng. Technol. 1(3), 98–106 (2009).

CAS Google Scholar

Rajesh, V. Radiation effects on MHD free convection flow near a vertical plate with ramped wall temperature. Int. J. Appl. Math. Mech. 6(21), 60–677 (2010).

MATH Google Scholar

Jafar, K., Nazar, R., Ishak, A. & Pop, I. MHD flow and heat transfer over stretching/shrinking sheets with external magnetic field, viscous dissipation and joule effects. Can. J. Chem. Eng. 90(5), 1336–1346 (2012).

Article CAS MATH Google Scholar

Sharma, B., Singh, A., Yadav, K. & Chaudhary, R. Effects of chemical reaction on magneto-micropolar fluid flow from a radiative surface with variable permeability. Int. J. Appl. Mech. Eng. 18(3), 833–851 (2013).

Article Google Scholar

Sharma, B. K., Tailor, V. & Goyal, M. Heat source and soret effects on megneto-micropolar fluid flow with variable permeability and chemical reaction. Glob. J. Pure Appl. Math. 13(9), 5195–5212 (2017).

Google Scholar

Waqas, M. et al. Magnetohydrodynamic (MHD) mixed convection flow of micropolar liquid due to nonlinear stretched sheet with convective condition. Int. J. Heat Mass Transf. 102, 766–772 (2016).

Article Google Scholar

Srinivasulu, T. & Goud, B. S. Effect of inclined magnetic field on flow, heat and mass transfer of Williamson nanofluid over a stretching sheet. Case Stud. Thermal Eng. 23, 100819 (2021).

Article Google Scholar

Walelign, T., Haile, E., Kebede, T. & Awgichew, G. Heat and mass transfer in stagnation point flow of Maxwell nanofluid towards a vertical stretching sheet with effect of induced magnetic field. Math. Probl. Eng. 20, 21 (2021).

MathSciNet Google Scholar

Dawar, A., Islam, S., Shah, Z., Mahmuod, S. & Lone, S. A. Dynamics of inter-particle spacing, nanoparticle radius, inclined magnetic field and nonlinear thermal radiation on the water-based copper nanofluid flow past a convectively heated stretching surface with mass flux condition: A strong suction case. Int. Commun. Heat Mass Transf. 137, 106286 (2022).

Article CAS Google Scholar

Gandhi, R. & Sharma, B. K. Unsteady mhd hybrid nanoparticle (Au-Al\(_2\)O\(_3\)/Blood) mediated blood flow through a vertical irregular stenosed artery: Drug delivery applications. In Nonlinear Dynamics and Applications 325–337 (Springer, 2022).

Chapter Google Scholar

Dawar, A., Islam, S. & Shah, Z. A comparative analysis of the performance of magnetised copper-copper oxide/water and copper-copper oxide/kerosene oil hybrid nanofluids flowing through an extending surface with velocity slips and thermal convective conditions. Int. J. Ambient Energy 20, 1–19 (2022).

Article Google Scholar

Hayat, T., Muhammad, T., Qayyum, A., Alsaedi, A. & Mustafa, M. On squeezing flow of nanofluid in the presence of magnetic field effects. J. Mol. Liq. 213, 179–185 (2016).

Article CAS Google Scholar

Sulochana, C., Ashwinkumar, G. & Sandeep, N. Transpiration effect on stagnation-point flow of a Carreau nanofluid in the presence of thermophoresis and brownian motion. Alex. Eng. J. 55(2), 1151–1157 (2016).

Article Google Scholar

Reddy, J. R., Sugunamma, V. & Sandeep, N. Thermophoresis and Brownian motion effects on unsteady MHD nanofluid flow over a slendering stretching surface with slip effects. Alex. Eng. J. 57(4), 2465–2473 (2018).

Article Google Scholar

Shah, Z. et al. Three dimensional third grade nanofluid flow in a rotating system between parallel plates with Brownian motion and thermophoresis effects. Res. Phys. 10, 36–45 (2018).

CAS Google Scholar

Soomro, F. A., Haq, R. U. & Hamid, M. Brownian motion and thermophoretic effects on non-Newtonian nanofluid flow via Crank-Nicolson scheme. Arch. Appl. Mech. 91(7), 3303–3313 (2021).

Article ADS Google Scholar

Rooman, M., Jan, M. A., Shah, Z., Kumam, P. & Alshehri, A. Entropy optimization and heat transfer analysis in mhd williamson nanofluid flow over a vertical riga plate with nonlinear thermal radiation. Sci. Rep. 11(1), 1–14 (2021).

Article Google Scholar

Dawar, A., Saeed, A. & Kumam, P. Magneto-hydrothermal analysis of copper and copper oxide nanoparticles between two parallel plates with Brownian motion and thermophoresis effects. Int. Commun. Heat Mass Transf. 133, 105982 (2022).

Article CAS Google Scholar

Bég, O. A., Zueco, J. & Takhar, H. S. Unsteady magnetohydrodynamic Hartmann–Couette flow and heat transfer in a Darcian channel with Hall current, ionslip, viscous and Joule heating effects: Network numerical solutions. Commun. Nonlinear Sci. Numer. Simul. 14(4), 1082–1097 (2009).

Article ADS Google Scholar

Sahoo, B. Effects of slip, viscous dissipation and Joule heating on the MHD flow and heat transfer of a second grade fluid past a radially stretching sheet. Appl. Math. Mech. 31(2), 159–173 (2010).

Article MathSciNet MATH Google Scholar

Sharma, B., Sharma, M., Gaur, R. & Mishra, A. Mathematical modeling of magneto pulsatile blood flow through a porous medium with a heat source. Int. J. Appl. Mech. Eng. 20(2), 385–396 (2015).

Article Google Scholar

Hayat, T., Waqas, M., Shehzad, S. & Alsaedi, A. A model of solar radiation and Joule heating in magnetohydrodynamic (MHD) convective flow of thixotropic nanofluid. J. Mol. Liq. 215, 704–710 (2016).

Article CAS Google Scholar

Ramzan, M., Farooq, M., Hayat, T. & Chung, J. D. Radiative and Joule heating effects in the MHD flow of a micropolar fluid with partial slip and convective boundary condition. J. Mol. Liq. 221, 394–400 (2016).

Article CAS Google Scholar

Hsiao, K.-L. Combined electrical MHD heat transfer thermal extrusion system using Maxwell fluid with radiative and viscous dissipation effects. Appl. Therm. Eng. 112, 1281–1288 (2017).

Article CAS Google Scholar

Hsiao, K.-L. Micropolar nanofluid flow with MHD and viscous dissipation effects towards a stretching sheet with multimedia feature. Int. J. Heat Mass Transf. 112, 983–990 (2017).

Article Google Scholar

Gayatri, M., Reddy, K. J. & Babu, M. J. Slip flow of Carreau fluid over a slendering stretching sheet with viscous dissipation and Joule heating. SN Appl. Sci. 2(3), 1–11 (2020).

Article Google Scholar

Seethamahalskshmi, V., Ramana Reddy, G., Sandhya, A. & Sateesh Kumar, D. Study MHD mixed convective flow of a vertical porous surface in the presence viscous dissipation. In AIP Conference Proceedings, vol. 2375, p. 030009 (AIP Publishing LLC, 2021).

Gandhi, R., Sharma, B., Kumawat, C. & Bég, O. A. Modeling and analysis of magnetic hybrid nanoparticle (Au–Al\(_2\)O\(_3\)/blood) based drug delivery through a bell-shaped occluded artery with Joule heating, viscous dissipation and variable viscosity effects. In Proceedings of the Institution of Mechanical Engineers, Part E: Journal of Process Mechanical Engineering, 09544089221080273 (2022).

Rashidi, M., Kavyani, N. & Abelman, S. Investigation of entropy generation in MHD and slip flow over a rotating porous disk with variable properties. Int. J. Heat Mass Transf. 70, 892–917 (2014).

Article Google Scholar

Dalir, N., Dehsara, M. & Nourazar, S. S. Entropy analysis for magnetohydrodynamic flow and heat transfer of a Jeffrey nanofluid over a stretching sheet. Energy 79, 351–362 (2015).

Article Google Scholar

Baag, S., Mishra, S., Dash, G. & Acharya, M. Entropy generation analysis for viscoelastic MHD flow over a stretching sheet embedded in a porous medium. Ain Shams Eng. J. 8(4), 623–632 (2017).

Article Google Scholar

Bhatti, M., Rashidi, M. & Pop, I. Entropy generation with nonlinear heat and mass transfer on MHD boundary layer over a moving surface using SLM. Nonlinear Eng. 6(1), 43–52 (2017).

Article ADS Google Scholar

Khan, M. I., Hayat, T., Waqas, M., Khan, M. I. & Alsaedi, A. Entropy generation minimization (EGM) in nonlinear mixed convective flow of nanomaterial with Joule heating and slip condition. J. Mol. Liq. 256, 108–120 (2018).

Article CAS Google Scholar

Sohail, M., Shah, Z., Tassaddiq, A., Kumam, P. & Roy, P. Entropy generation in MHD Casson fluid flow with variable heat conductance and thermal conductivity over non-linear bi-directional stretching surface. Sci. Rep. 10(1), 1–16 (2020).

Article Google Scholar

Hayat, T., Qayyum, S., Alsaedi, A. & Ahmad, B. Entropy generation minimization: Darcy–Forchheimer nanofluid flow due to curved stretching sheet with partial slip. Int. Commun. Heat Mass Transf. 111, 104445 (2020).

Article CAS Google Scholar

Sharma, B., Gandhi, R. & Bhatti, M. Entropy analysis of thermally radiating MHD slip flow of hybrid nanoparticles (Au-Al\(_2\)O\(_3\)/Blood) through a tapered multi-stenosed artery. Chem. Phys. Lett. 20, 139348 (2022).

Article Google Scholar

Gandhi, R., Sharma, B. & Makinde, O. Entropy analysis for MHD blood flow of hybrid nanoparticles (Au–Al\(_2\)O\(_3\)/blood) of different shapes through an irregular stenosed permeable walled artery under periodic body acceleration: Hemodynamical applications. J. Appl. Math. Mech.https://doi.org/10.1002/zamm.202100532 (2022).

Article Google Scholar

Sharma, B., Kumar, A., Gandhi, R. & Bhatti, M. Exponential space and thermal-dependent heat source effects on electro-magneto-hydrodynamic jeffrey fluid flow over a vertical stretching surface. Int. J. Mod. Phys. B 20, 2250220 (2022).

Article Google Scholar

Anki Reddy, P. B., Reddy, S. & Suneetha, S. Magnetohydrodynamic flow of blood in a permeable inclined stretching surface with viscous dissipation, non-uniform heat source/sink and chemical reaction. Front. Heat Mass Transf. 10, 25 (2018).

Article Google Scholar

Bird, R. B., Stewart, W. E., Lightfoot, E. N. & Meredith, R. E. Transport phenomena. J. Electrochem. Soc. 108(3), 78C (1961).

Article Google Scholar

Khanduri, U. & Sharma, B. K. Entropy analysis for mhd flow subject to temperature-dependent viscosity and thermal conductivity. In Nonlinear Dynamics and Applications 457–471 (Springer, 2022).

Chapter Google Scholar

Awais, M., Hayat, T., Ali, A. & Irum, S. Velocity, thermal and concentration slip effects on a magneto-hydrodynamic nanofluid flow. Alex. Eng. J. 55(3), 2107–2114 (2016).

Article Google Scholar

Sharma, B. & Gandhi, R. Combined effects of Joule heating and non-uniform heat source/sink on unsteady mhd mixed convective flow over a vertical stretching surface embedded in a Darcy-Forchheimer porous medium. Propuls. Power Res. 11(2), 276–292 (2022).

Article Google Scholar

Download references

The authors would like to acknowledge and express their gratitude to the United Arab Emirates University, Al Ain, UAE for providing fnancial support with Grant No. 12S086. Also, we would like to thank the reviewers for taking the time and making effort to review the manuscript. We sincerely appreciate all valuable comments and suggestions, which helped us to improve the quality of the manuscript.

Department of Mathematics, Birla Institute of Technology and Science Pilani, Pilani, Rajasthan, India

B. K. Sharma & Rishu Gandhi

Department of Basic Science, College of Science and Theoretical Studies, Saudi Electronic University, Riyadh, 11673, Saudi Arabia

Nidhish K Mishra

Department of Mathematical Sciences, College of Science, UAE University, P.O. Box 17551, Al-Ain, United Arab Emirates

Qasem M. Al-Mdallal

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

B.K.S.: conceptualization, drafting the manuscript, supervision, validation, review and editing. R.G.: conceptualization, drafting the manuscript, methodology, investigation, software handling, supervision, validation, review and editing. N.K.M.: methodology, investigation, software handling, analysis, interpretation of data. Q.M.A.-M.: analysis, interpretation of data, review and editing.

Correspondence to Qasem M. Al-Mdallal.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Sharma, B.K., Gandhi, R., Mishra, N.K. et al. Entropy generation minimization of higher-order endothermic/exothermic chemical reaction with activation energy on MHD mixed convective flow over a stretching surface. Sci Rep 12, 17688 (2022). https://doi.org/10.1038/s41598-022-22521-5

Download citation

Received: 22 May 2022

Accepted: 17 October 2022

Published: 21 October 2022

DOI: https://doi.org/10.1038/s41598-022-22521-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Journal of Thermal Analysis and Calorimetry (2023)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.