Significance of Chemical Reaction and Lorentz Force on Third-Grade Fluid Flow and Heat Transfer with Darcy–Forchheimer Law over an Inclined Exponentially Stretching Sheet Embedded in a Porous Medium

Abstract

The combined impact of a linear chemical reaction and Lorentz force on heat and mass transfer in a third-grade fluid with the Darcy–Forchheimer relation over an inclined, exponentially stretching surface embedded in a porous medium is investigated. The proposed process is mathematically expressed in terms of nonlinear and coupled partial differential equations, with the symmetry of the conditions normal to the surface. To solve the mathematical model of the proposed phenomenon, the partial differential equations are first reduced to ordinary differential equations; then, MATLAB built-in Numerical Solver bvp4c is used to obtain the numerical results of these equations. The influence of all the pertinent parameters that appeared in the flow model on the unknown material properties of interest is depicted in the forms of tables and graphs. The physical attitude of the unknown variables is discussed with physical reasoning. From the numerical solutions, it is inferred that, as Lorentz force parameter $M$ is increased, the velocity of the fluid decreases, but fluid temperature and mass concentration increase. This is due to the fact that Lorentz force retards the motion of fluid, and the increasing resistive force causes the rise in the temperature of the fluid. It is also noted that, owing to an increase in the magnitude of chemical reaction parameter $R$, the velocity profile and the mass concentration decline as well, but the fluid temperature increases in a reasonable manner. It is noted that, by augmenting the values of the local inertial coefficient $Fr$ and the permeability parameter $K^{*}$, the velocity field decreases, the temperature field increases, and mass concentration also increases with reasonable difference. Increasing values of Prandtl number $Pr$ results in a decrease in the profiles of velocity and temperature. All the numerical results are computed at the angle of inclination $\alpha =\pi /6$. The current results are compared with the available results in the existing literature for this special case, and there is good agreement between them that shows the validation of the present study. All the numerical results show asymptotic behavior by satisfying the given boundary conditions.

1. Introduction

The study of non-Newtonian fluid flows has received a great deal of interest in recent decades because of its importance in engineering and industry. Food processing, papermaking, and lubricating processes all rely heavily on non-Newtonian fluids with heat and mass transfer. Because of the widespread use of non-Newtonian fluids in industry, the research community has developed a variety of non-Newtonian fluid models. The third-grade fluid model is one of the most comprehensive fluid models, displaying all shear-thinning and shear-thickening fluid properties. The researchers focused a great deal of work on the numerical and analytical studies of Newtonian and rheological fluid flows, as well as purely theoretical studies in theoretical fluid dynamics with an emphasis on symmetry concepts arising from group studies. The symmetricity breaks under the conditions of a high Reynolds number. In light of the above, the outlined physical applications of such flow problems with different flow characteristics have been investigated by many researchers. The research work performed by the research community concentrating on the above-mentioned issues is highlighted in the following paragraphs.

Rashidi et al. [1] investigated the entropy generation on a magnetohydrodynamic third-grade fluid over a stretching sheet. Riaz et al. [2] investigated simultaneous heat and mass transfer in a non-Newtonian fluid under the effects of a magnetic field over a vertical surface. Ebrahimi et al. [3] took into study the mechanism of the fully developed flow of a third-grade fluid in the plane duct, with the assumption that there is convection on the walls. The heat transportation through the generalized Fourier’s law in a third-grade fluid under the impact of a chemical reaction over the stretching sheet has been proposed by Imtiaz et al. [4]. Alazwri et al. [5] explored the study of entropy optimization of the third-grade viscoelastic nanofluid flow over a stretching surface by using the Keller-Box Scheme. Sahoo and Do [6] numerically investigated the phenomenon of heat transfer in third-grade fluid flow via a linearly stretching sheet with the effects of slip condition. Jamshed et al. [7] encountered the study of an unsteady Casson nanofluid over the surface of a stretching sheet under the effects of thermal radiation. A fluid flow and heat transfer analysis on non-Newtonian nanofluid over a porous surface was performed by Maleki et al. [8]. Maleki et al. [9] paid attention to the heat transfer and fluid flow of a pseudo-plastic nanofluid over a moving porous surface, with viscous dissipation and heat generation and absorption effects. Abu-Hamdeh et al. [10] highlighted the physical significance of solar energy by considering the non-Newtonian Powell–Eyring nanofluid with the impact of the thermal jump. The study of hybridized Casson nanofluid on porous moving and fixed plates under the effects of suction and magnetic field effects was conducted by El-Zahar et al. [11].

By considering the magnetic field and convective boundary conditions, Jawanmard et al. [12] concentrated on the fully developed flow of a third-grade fluid in a pipe and solved the model using the Runge–Kutta method in order four. The investigation of third-grade fluid flow and heat transfer with a chemical reaction and magnetic field effects was undertaken by Hayat et al. [13]. Bioconvective magnetohydrodynamic nanofluid flow due to an exponentially stretching sheet in the presence of a chemical reaction effects was considered by Asjad et al. [14]. The entropy generation analysis of magneto-third-grade fluid flow through a micro-channel was performed in [15]. Zhang et al. [16] investigated the electro magnetohydrodynamic behavior of third-grade fluid flow between parallel plates, taking into account the effects of Joule heating and dissipation. Shehzad et al. [17] researched heat transportation and flow analysis in a third-grade fluid using generalized transfer principles and an exponentially stretched sheet. Reddy et al. [18] investigated the unsteady flow of third-grade fluid down an evenly heated cylinder when vertically viewed using the Bejan function. Ref. [19] focuses on the bioconvection heat transfer via a third-grade nanofluid over a stretched surface, with significant effects on the activation energy and magnetic field. The analytical study of third-grade fluid flow over a rotating cone in the presence of nanoparticles was the focus of Nadeem and Saleem [20].

The study of flows saturating porous space has piqued the interest of researchers due to their importance in engineering and industry applications, such as oil reservoirs, resin transfer models, porous insulation, packed beds, geothermal energy, fossil fuel beds, and nuclear waste disposal. The difficulties of permeable media that were developed using Darcy’s relation, in which the pressure gradient and volume average velocity are directly associated, have received a great deal of attention in the literature. Darcy’s idea is only valid when a smaller porosity and lesser velocity are present. Because non-uniform porosity and greater velocities occur in different industrial and technical systems, Darcy’s relation is ineffectual. To attain better results, it is highly desirable to include non-Darcian properties in the model for permeable materials. Inertial effects and high-speed flow are mostly described by Forchheimer’s law. Forchheimer’s law is unique in that it combines the monotonicity of the nonlinear term with the non-degeneracy of Darcy’s component. To examine the properties of inertia and boundary, Forchheimer [21] included the factor of square velocity in the calculation of Darcy’s velocity. This word was coined by Muskat and Wyckoff [22] as the Forchheimer term, which always applies to issues with a high Reynolds number. The combined effect of the heat source, sink and Ohmic heating on the three-dimensional Darcy–Forchheimer flow of a CNTs nanofluid over a stretching sheet are examined by Upreti et al. [23].

The chemical and suction/injection effects on magnetic micropolar fluid through a porous wedge with Hall and ion slip conditions have been numerically solved using the Keller-Box Scheme by Sing et al. [24]. Ramzan et al. [25] perform a theoretical analysis of a Williamson nanofluid with the Darcy–Forchheimer relation in a stratified medium, taking into account the impacts of extended Fourier’s and Fick’s laws, the magnetic field, and the chemical reaction. By encountering generalized Fourier’s and Fick’s laws, Hayat et al. [26] addressed the ideal research of the Darcy–Forchheimer law in anon-Newtonian fluid. Upreti et al. [27] proposed the phenomenon of a Siskonanofluid flow with the Darcy–Forchheimer relation over a stretching sheet fixed in a porous medium under the influence of thermal radiation. The study of the Darcy–Forchheimer relation’s impact on chemically reacting fluid flow over a convectively heated surface with the inclusion of thermal radiation effects was documented by Sadiq et al. [28]. The study on the impact of Darcy–Forchheimer magnetohydrodynamic Carreau–Yasuda nanofluid flow was documented by Khan et al. [29]. The Darcy–Forchheimer effects were used to analyze the Jeffery nanofluid flow through a curved stretching surface under the influence of thermal radiation, a chemical reaction, and convective boundary conditions by Nagaraja et al. [30].

The heat and mass transmission properties of stretched sheets have important engineering and industrial applications. This comprises plastic sheet extrusion, blood pumping, metal processing, pharmaceuticals, fiber production processes, wire coating, and so forth. Due to their widespread use in recent years, researchers have turned their attention to various forms of flow. Kumar et al. [31] investigated heat and mass transmission on exponentially angled stretching sheets embedded in porous media using Soret, Dufour, a magnetic field, slip effects, Joule heating, and a chemical reaction. A heat and mass transfer analysis in the boundary layer flow over an exponentially stretching sheet was considered by Magyari and Keller [32]. The researchers in refs. [33,34,35,36] investigated the numerical study of several processes with varying flows over an exponentially stretched sheet.

From a technological standpoint, magnetohydrodynamic flows in experimental and theoretical research are essential because they have various applications, such as magnetohydrodynamic electrical power generation, geophysics, and so on. There is much physical significance of the magnetic field in the study of the heat transportation process; the idea of the magnetic field in the boundary layer flow has attracted a great deal of attention from researchers. Therefore, an enormous amount the research focusing on the impact of the magnetic field on the different fluid characteristics, along with several geometries with different flow conditions, has been encountered. The magnetohydrodynamic Casson fluid flow over a stretching cylinder was analyzed by Tamoor et al. [37]. The influence of the magnetic field and a chemical reaction on convective heat and mass transfer over an exponentially stretching sheet was studied by Pattnaik et al. [38]. Magnetohydrodynamic flow with thermal radiation effects over an exponentially stretching sheet was tackled numerically by using the Homotopy Analysis Method in [39]. Salahuddin et al. [40] documented the research study of the topic of Williamson fluid flow via the generalized heat transportation law over the stretching sheet with a variable thickness in the presence of the magnetic field. Kumar et al. [41] explored the process of magnetohydrodynamic flow past a cone and wedge with a variable heat source and sink. The effects of the magnetic field of the fluid flow and heat transfer with different flow conditions and fluid characteristics on the diverse flow geometries are discussed in [42,43,44,45,46].

The study of heat transfer and mass transfer with chemical reactions is of most relevant significance to engineers and scientists because of its universal incidence in many branches of science and engineering. This phenomenon plays a significant role in the chemical industry, the power and cooling industry for drying, evaporation, energy transfer in a cooling tower, and the flow in a desert cooler, etc. Examining transportation techniques and the relationship between these tactics and chemical reactions can be difficult and is closely linked to basic fluid dynamics. In creating applications, such as food processing and atomized liquid fuel ignition, such a thorough inspection and chemical techniques are useful. All industrial and chemical processes are designed to transform low-cost raw materials into high-value products (for the most part, through chemical reactions). Because of its importance, researchers took into account the chemical reactions of various orders in various fluid flow models with various model circumstances and fluid features, taking into account the various flow geometries presented in [47,48,49,50,51]. Misra and Adhkary [52] focused their attention on the magnetohydrodynamic oscillatory channel flow in physiological fluid in the presence of a chemical reaction. Mishra et al. [53] discussed the chemical reaction effects on viscoelastic fluid over a stretching sheet, along with the magnetic field effects. Sivakami et al. [54] shed light on the mechanism of unsteady, free convective immiscible fluid flow in a horizontal channel under the impact of a magnetic field and a chemical reaction. Lv et al. [55] investigated the effects of a chemical reaction and thermal radiation on convective heat and mass transfer over a rotating channel with Hall currents. Fetecau and Vieru [56] gave atheoretical analysis on the symmetric and non-symmetric flows of Burger’s fluid over a sheet embedded in a porous medium. Zeshan et al. [57] studied the electromagnetic flow of SWCN/MWCNT suspensions in two immiscible Newtonian fluids through a porous medium. An entropy generation analysis on magnetohydrodynamic nanofluid flow with the Darcy–Forchheimer relation due to a nonlinear stretching sheet was performed by Rasool et al. [58]. Prakash et al. [59] conducted the discussion on the topic of Peristaltic nanofluid flow through a tapered channel embedded in a porous medium. Ullah et al. [60] conducted the analysis on magnetohydrodynamic Casson fluid flow over a porous stretching cylinder fixed in a porous medium with a chemical reaction, viscous dissipation, heat generation, and absorption effects.

In the already-known literature, due to the widespread use of non-Newtonian fluids in engineering and industry (food processing, papermaking, and lubricating processes), the research community has paid a great deal of attention to the study of third-grade fluid flow and heat transfer on different flow geometries by considering the different fluid characteristics separately. Due to the practical applications of non-Newtonian fluids in industrial processes (food processing, papermaking, and lubricating processes), along with chemical reactions that are used in the chemical industry, the power and cooling industry for drying, evaporation, energy transfer in a cooling tower, and the flow in a desert cooler, etc. in the presence of a magnetic field that has physical applications in electrical power generation and geophysics, the combined study of the effects of the applied magnetic field (Lorentz force effects) and a linear chemical reaction with the Darcy–Forchheimer relation on the heat and mass transfer in non-Newtonian third-grade fluid along the surface of an inclined exponentially stretching sheet embedded in a porous medium is discussed numerically, which has not been encountered before the present work. In Section 2, a mathematical model for the proposed mechanism is developed. In Section 3, the proposed model, in terms of partial differential equations, is then reduced to a set of ordinary equations by sung similarity transformation and is then solved with the help of the MATLAB built-in numerical solver bvp4c. The obtained results regarding physical properties are presented and discussed in Section 4 and Section 5. At the end, the entire study is concluded.

2. Mathematical Modeling

Consider the following assumptions to make the model for the current study:

• Steady flow;
• Two-dimensional third-grade fluid flow;
• Incompressible flow;
• Viscous flow;
• Inclusion of Darcy–Forchheimer relation;
• Inclined stretching sheet embedded in a porous medium;
• Inclusion of Lorentz force effects (applied magnetic field);
• Chemical reactions are taking place in the flow field, and the constant $K_1$ is the first-order chemical reaction rate is taken as $>0$ (destructive reactant);
• The axes $x, y$ are horizontal and normal at the surface of an exponentially inclined stretching sheet. The schematic diagram to show the flow coordinate system and boundary layer is given in Figure 1.

By keeping the above-said assumptions in mind and following refs. [6,31], the flow equations are given below:

$$\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}=0$$

(1)

$$\begin{gathered} u\frac{\partial u}{\partial x}+v\frac{\partial u}{\partial y}=\nu \frac{{\partial }^{2}u}{\partial y^2}+\frac{{\alpha }_1}{\rho }\left(u\frac{{\partial }^{3}u}{\partial x\partial y^2}+\frac{\partial u}{\partial x}\frac{{\partial }^{2}u}{\partial y^2}+3\frac{\partial u}{\partial y}\frac{{\partial }^{2}v}{\partial y^2}+\nu \frac{{\partial }^{3}u}{\partial y^3}\right) \\ +2\frac{{\alpha }_2}{\rho }\frac{\partial u}{\partial y}\frac{{\partial }^{2}v}{\partial y^2} \\ +6\frac{{\beta }_3}{\rho }\left(\frac{\partial u}{\partial y}\right){\frac{{\partial }^{2}u}{\partial y^2}}^2+g{\beta }_T\left(T-T_{\infty }\right)cos\alpha +g{\beta }_C\left(C-C_{\infty }\right)cos\alpha \\ -\frac{\sigma B{\left(x\right)}^2}{\rho }-\frac{\nu u}{K_o}-F{u}^2 \end{gathered}$$

(2)

$$u\frac{\partial T}{\partial x}+v\frac{\partial T}{\partial y}=\alpha \frac{{\partial }^{2}T}{\partial y^2}$$

(3)

$$u\frac{\partial C}{\partial x}+v\frac{\partial C}{\partial y}=D_m\frac{{\partial }^{2}C}{\partial y^2}-K_1\left(C-C_{\infty }\right)$$

(4)

The molded conditions are:

$$\begin{gathered} u=U_w, v=0, T=T_w\left(x\right), C=C_w\left(x\right) \mathrm{at} y=0 \\ u\to 0\hspace{1em}T\to T_{\infty },\hspace{1em}C\to C_{\infty }\hspace{1em}as\hspace{1em}y\to \infty \end{gathered}$$

(5)

where $U_w=U_o{e}^{x/L}, T_w=T_{\infty }+C_o{e}^{x/L},C_w=C_{\infty }+C_o{e}^{x/L},$ representingthe stretching velocity, wall temperature, and wall concentration, respectively. The symbols $U_o,T_o, \mathrm{and} C_o$ are the reference velocity, temperature, and concentration, respectively. Here, $\left(u,v\right)$ are the velocity components along $\left(x,y\right)$ directions, respectively. The notations $({\alpha }_1, {\alpha }_2, {\beta }_3),K_o, F=\frac{{\mathrm{C}}_{\mathrm{b}}}{\sqrt{{\mathrm{K}}_{\mathrm{o}}}}$, $\mathsf{\sigma },\mathrm{B}\left(\mathrm{x}\right)={\mathrm{B}}_{\mathrm{o}}{\mathrm{e}}^{\frac{\mathrm{x}}{\mathrm{L}}}$, $K_1$ are the material moduli, permeability of the porous medium, coefficient of inertia with $C_b$ being the drag coefficient, electrical conductivity, magnetic field strength and chemical reaction coefficient, respectively. The designations $\rho , C_P, ,\nu , {\alpha }_m, D_m$ are the density of the fluid, specific heat at constant pressure, dynamic viscosity, kinematic viscosity, thermal diffusivity, and mass diffusivity, respectively. The symbols $\left(T, C\right), \left(T_{\infty }, C_{\infty }\right)$ are the temperature and concentration within the boundary layer and in the free stream region, respectively.

3. Solution Methodology

Here, we elaborate numerically on the solution methodology to tackle the equations given in Equations (1)–(4), along with boundary conditions (5), for the purpose of solving. The entire solution method is discussed in this section.

3.1. Similarity Formulation

The Equations (1)–(5) are nonlinear partial differential equations; first, we reduce these equations into a set of ordinary differential equations by using the following similarity variables used by [31]:

$$\begin{gathered} u=U_o{e}^{\frac{x}{L}}, v=\sqrt{\frac{\nu U_o}{2L}}\left(f\left(\eta \right)+\eta f^{\prime }\left(\eta \right)\right)e^{\frac{x}{L}}, T=T_{\infty }+C_o\theta e^{\frac{x}{L}},C \\ =C_{\infty }+C_o\varphi e^{\frac{x}{L}}, \\ \eta =\sqrt{\frac{U_o}{2L\nu }}y{e}^{x/L} \end{gathered}$$

(6)

By using the above-mentioned similarity variables given in Equation (6), Equation (1) is satisfied automatically, and Equations (2)–(4) with boundary conditions from Equation (5) adopt the following form:

$$\begin{gathered} f^{‴}+f{f}^{″}-2f^{{\prime }^2}+2Ri\left(\theta +N\varphi \right)cos\alpha \\ +K\left(6f^{\prime }{f}^{‴}-f{f}^{\left(iv\right)}-2\eta f^{″}{f}^{‴}-9f^{{″}^2}\right) \\ -L\left(3f^{{″}^2}+\eta f^{″}{f}^{‴}\right)+3\beta Re{f}^{{″}^2}{f}^{‴}-K^{*}{f}^{\prime }-M{f}^{\prime }-Fr{f}^{{\prime }^2}=0 \end{gathered}$$

(7)

$$f^{\prime }\theta -f{\theta }^{\prime }=\frac{1}{Pr}{\theta }^{″}$$

(8)

$$f^{\prime }\varphi -f{\varphi }^{\prime }=\frac{1}{Sc}{\varphi }^{\prime }-R\varphi$$

(9)

Modeled conditions:

$$\begin{gathered} f\left(\eta \right)=0, f^{\prime }\left(\eta \right)=1, \theta \left(\eta \right)=1, \varphi \left(\eta \right)=1 at \eta =0 \\ {f}^{\prime }\left(\eta \right)\to 0,f^{″}\left(\eta \right)\to 0, \theta \left(\eta \right)\to 0, \varphi \left(\eta \right)\to 0 as \eta \to \infty \end{gathered}$$

(10)

Here, the parameters $Ri=\frac{Gr}{R{e}^2}, N=\frac{{\beta }_C{C}_o}{{\beta }_T{T}_o}, K=\frac{{\alpha }_1{U}_o e^{\frac{x}{L}}}{2\rho \nu L}, L=\frac{{\alpha }_2{U}_o e^{\frac{x}{L}}}{\rho \nu L}, \beta =\frac{{\beta }_3{U}_o e^{\frac{x}{L}}}{\rho \nu L}, K^{*}=\frac{2\nu L}{K_1{U}_o{e}^{\frac{x}{L}}}$, $M=\frac{2\sigma B_o^2}{\rho U_o}, Fr=\frac{C_{b}2L}{\sqrt{K_o}},$ Pr = $\frac{\nu }{\alpha }$, $Sc$=$\frac{\nu }{D_m}$, $R=\frac{2L{K}_1}{U_o{e}^{\frac{x}{L}}}$ are the Richardson number, the buoyancy ratio parameter, the non-dimensional viscoelastic parameter, the cross-viscous parameter, the third-grade fluid parameter, the permeability parameter, the local inertial coefficient, the Prandtl number, the Schmidt number, and the chemical reaction parameter, respectively. Here, $Re=\frac{U_{o}L}{\nu }$ and $Gr=\frac{g{\beta }_T\left(T_w-T_{\infty }\right)L^3}{{\nu }^2}$ are the Reynolds number and the Grash of number, respectively.

The mathematical expressions for the drag force, the Nusselt number, and the Sherwood number are:

$$C_f=\frac{{\tau }_w}{\rho U_w^2},Nu=\frac{q_{w}x}{k\left(T-T_{\infty }\right)},Sh=\frac{q_{m}x}{D_m\left(C-C_{\infty }\right)}$$

(11)

where

$$\begin{gathered} {\tau }_w=\left(\frac{\partial u}{\partial y}+\frac{{\alpha }_1}{\mu }\left(2\frac{\partial u}{\partial x}\frac{\partial u}{\partial y}+v\frac{{\partial }^{2}u}{\partial y^2}+u\frac{{\partial }^{2}u}{\partial x\partial y}\right)+\frac{2{\beta }_3}{\mu }{\left(\frac{\partial u}{\partial y}\right)}^2\right) \\ {q}_w=-k\left(\frac{\partial T}{\partial y}\right), q_m=-D_m\left(\frac{\partial C}{\partial y}\right)\mathrm{at}y=0 \end{gathered}$$

(12)

Are the typical stress tensor, heat, and mass flux at the surfaces, respectively. Finally, the expressions for drag force, local Nusselt, and local Sherwood are given below:

$$\left.\begin{array}{c}{C}_{f_L}=\frac{2}{\sqrt{R{e}_L}}\left[f^{″}\left(0\right)+K\left(3f^{\prime }\left(0\right)f^{″}\left(0\right)-f\left(0\right)f^{″}\left(0\right)\right)+2\beta R{e}_x{f}^{{″}^2}\left(0\right)\right],\\ R{e}_L^{-\frac{1}{2}}N{u}_L=-\theta \left(0\right),\\ R{e}_L^{-\frac{1}{2}}S{h}_L=-\varphi \left(0\right),\end{array}\right\}$$

(13)

where $R{e}_L=\frac{U_{w}L}{\nu }$ is the local Reynolds number parallel to the $x$-direction.

3.2. Solution Technique

The exact solution for highly nonlinear coupled ordinary differential equations is not possible. Therefore, the numerical solutions for the above Equations (7)–(9), along with boundary conditions (10), are obtained for different values of the governing parameters, namely, the Richardson number $Ri$, the buoyancy ratio parameter $N$, the third-grade fluid parameter $\beta$, the viscoelastic parameter $K$, the cross-viscous parameter $L$, the permeability parameter $K^{*}$, the local inertial coefficient $Fr$, the magnetic field parameter $M$, the Prandtl number Pr, the Schmidt number Sc, and the chemical reaction parameter $R$. The numerical solutions for the proposed model are solved by utilizing the MATLAB built-in Numerical Solver bvp4c. In the computation, η∞ = 10.0 is taken, and the axis is defined according to the clear figure visibility. The Numerical Solver bvp4c is a finite difference code that implements the three-stage Lobatto formula. This is a collocation formula, and the collocation polynomial is a C1-continuous solution with fourth-order accuracy uniformly in the interval of integration. Mesh selection and error control are based on the residual of the continuous solution. The collocation technique uses a mesh of the points to divide the interval of integration into subintervals. The solver determines a numerical solution by solving a global system of algebraic equations resulting from the boundary conditions and the collocation condition imposed on all the subintervals. The solver estimates the error of the numerical solution on each subinterval. If the solution does not satisfy the tolerance criteria, the solver adapts the mesh and repeats the process. There is a need to provide the points of the initial mesh, as well as an initial approximation of the solution at the mesh points. The results produced with the help of this numerical method in the current work are compared with those previously published, which indicate the excellent agreement between both results and show the accuracy and validation of the present results. The Equations (7)–(9), along with boundary condition (10), are converted into a system of first-order ordinary differential equations and are then put into the bvp4c numerical solver code for the final solutions. We set as below:

$$\begin{gathered} f=y\left(1\right), f^{\prime }=y\left(2\right),f^{″}=y\left(3\right), f^{‴}=y\left(4\right),\theta =y\left(5\right), {\theta }^{\prime }= \\ y\left(6\right), \varphi =y\left(7\right),{\varphi }^{\prime }=y\left(8\right) \end{gathered}$$

(14)

$$\begin{array}{l}{f}^{\left(iv\right)}=yy1=\left(1/K\ast y\left(1\right)\right)\ast (2\ast Ri\ast \left(N\ast y\left(7\right)+y\left(5\right)\right)\ast \cos \alpha \\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}+y\left(4\right)+y\left(1\right)\ast y\left(3\right)-2\ast y{\left(2\right)}^2+K\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\ast \left(6\ast y\left(2\right)\ast y\left(4\right)-2\ast \eta \ast y\left(3\right)\ast y\left(4\right)-9\ast y{\left(3\right)}^2\right)\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}-L\ast \left(3\ast y{\left(3\right)}^2+\eta \ast y\left(2\right)\ast y\left(4\right)\right)+3\ast \beta \ast Re\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\ast y{\left(3\right)}^2\ast y\left(4\right)-y\left(2\right)\ast (K^{*}+M+y\left(2\right)+Fr\\ \hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\hspace{1em}\ast y\left(2\right))))\end{array}$$

(15)

$${\theta }^{″}=yy2=Pr*(y(2)*y(5)-y(1)*y(6))$$

(16)

$${\varphi }^{″}=\mathrm{Sc}*(\mathrm{y}(2)*\mathrm{y}(7)-\mathrm{y}(1)*\mathrm{y}(8))-\mathrm{R}*\mathrm{y}(7))$$

(17)

Boundary Conditions:

$$\begin{gathered} \mathrm{y}0\left(1\right)=0,\mathrm{y}0\left(2\right)=1,\mathrm{y}0\left(5\right)=1,\mathrm{y}0\left(7\right)=1,\mathrm{yinf}\left(2\right)=0,\mathrm{yinf}\left(3\right) \\ =0,\mathrm{yinf}\left(5\right)=0,\mathrm{yinf}\left(7\right)=0 \end{gathered}$$

(18)

The numerical solutions for the unknown material properties, such as the velocity profile $f^{\prime }\left(\eta \right)$, the temperature field $\theta \left(\eta \right),$ and the mass concentration $\varphi \left(\eta \right)$ profiles, are calculated at an angle of inclination $\alpha =\pi /6$, presented in graphical form, and discussed with physical reasoning. The gradients of the quantities outline above, such as the skin friction coefficient $f^{″}\left(0\right)$, the heat transfer rate $-{\theta }^{″}\left(0\right)$, and the mass transfer rate $-{\varphi }^{″}\left(0\right)$ at the surface of the geometry, are calculated and tabulated. The forthcoming sections are devoted to an analysis with detailed results and a discussion of the graphed and tabulated physical quantities of interest.

4. Results

The proposed mechanism was converted into the mathematical model given in Equations (1)–(4), along with the boundary conditions (5). This system of partial differential equations was reduced to a system of ordinary differential equations by using the similarity variable formulation given in Equation (6). The reduced form of ordinary differential equations was converted to a system of first-order differential equations. The first-order ordinary differential equations were solved numerically in the MATLAB built-in numerical solver bv4c. The numerical results of the velocity profile $f^{\prime }\left(\eta \right)$, the temperature profile $\theta \left(\eta \right)$, and mass concentration $\varphi \left(\eta \right)$, along with the skin friction coefficient $f^{″}\left(0\right)$, the rate of heat transfer ${\theta }^{\prime }\left(0\right)$, and the rate of mass transfer ${\varphi }^{\prime }\left(0\right)$ were computed. The parameters against which the numerical solutions of the physical properties were named as the Richardson number $Ri$, the buoyancy ratio parameter $N$, the viscoelastic parameter $K$, the cross-viscous parameter $L$, the third-grade fluid parameter $\beta$, the permeability parameter $K^{*}$, the local inertial coefficient $Fr$, the magnetic field parameter $M$, the Prandtl number Pr, the Schmidt number Sc, and the chemical reaction parameter $R$. All the numerical solutions were calculated at an angle of inclination $\alpha =\pi /6$ of a stretching surface from the $x$-axis. The appropriate ranges of the parameters were taken as $1\le Ri\le 7$; $1\le \beta \le 12$; $0.1\le M,K^{*}, Fr,R, Sc,N\le 7$; and $0.72\le Pr\le 7$. The values of ${\eta }_{\infty }$, the numerical infinity values, were kept large enough and were retained at ${\eta }_{\infty }=10.0$. Actually, this value is dependent on the physical parameters of the phenomenon, and its value $=10.0$ was adequate to simulate $\eta =\infty$ for all cases presented in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21 for the graphical representation of the velocity profile $f^{\prime }\left(\eta \right)$, the temperature profile $\theta \left(\eta \right)$, and mass concentration $\varphi \left(\eta \right)$. The numerical results for the skin friction coefficient $f^{″}\left(0\right)$, the rate of heat transfer ${\theta }^{\prime }\left(0\right)$, and the rate of mass transfer ${\varphi }^{\prime }\left(0\right)$ were computed exactly at the surface of the geometry. The current numerical results were compared with existing results for special cases for the validation of the current study.

5. Discussion

This section is devoted to the detailed discussion of the numerical solutions that were plotted and tabulated. The impact of each pertinent parameter on the velocity field $f^{\prime }\left(\eta \right)$, the temperature field $\theta \left(\eta \right)$, and mass concentration $\varphi \left(\eta \right)$, along with the skin friction coefficient $f^{″}\left(0\right)$, the Nusselt number, and the Sherwood number, are deliberately discussed from a physical viewpoint.

5.1. Impact of the Involved Parameters on Velocity Profile, Temperature Profile, and Mass Concentration

The physical significance of the Richardson number $Ri$, also known as the mixed convection parameter, on velocity $f^{\prime }\left(\eta \right)$, temperature field $\theta \left(\eta \right)$, and mass concentration $\varphi \left(\eta \right)$ is shown in Figure 2, Figure 3 and Figure 4. From the graphical presentation of the results, it can be inferred that, owing to an enhanced $Ri$, the property $f^{\prime }\left(\eta \right)$ grows as well, but the temperature and concentration fields decline in a visible manner. The impact of the third-grade fluid parameter $\beta$ on $f^{\prime }\left(\eta \right), \theta \left(\eta \right)$, and $\varphi \left(\eta \right)$ is illustrated in Figure 5, Figure 6 and Figure 7. It is concluded that, against the increasing value of $\beta$, an intensification in the velocity of the fluid and an attenuation in the temperature and concentration profiles are seen. By increasing $\beta$, the reference velocity is increased, and the viscosity is decreased, which causes the reduction in fluid velocity. The physical significance of the Lorentz force parameter, also termed as the magnetic field parameter or Hartmann number $M$, on $f^{\prime }\left(\eta \right), \theta \left(\eta \right)$, and $\varphi \left(\eta \right)$ is displayed in Figure 8, Figure 9 and Figure 10. From the graphical demonstration of the results, it is noted that, by increasing $M$, the velocity of the fluid reduces, but the temperature and concentration profiles are expanded. From a physical point of view, this phenomenon is true because, when $M$ is enhanced, Lorentz force is basically strengthened, which retards the movement of the fluid, causing a decrease in the velocity. By increasing the value of $M$, resistance within the fluid flow domain is generated, which causes the temperature to rise. The effect of the permeable parameter $K^{*}$ on the velocity, temperature, and concentration profiles is presented in Figure 11, Figure 12 and Figure 13. The graphs highlight that $f^{\prime }\left(\eta \right)$ decreases, but $\theta \left(\eta \right)$ and $\varphi \left(\eta \right)$ increase. Physically, this is true because, when $K^{*}$ is enlarged, the viscosity is augmented, and the porosity of the porous medium is decreased, due to which the velocity is slowed down. The $f^{\prime }\left(\eta \right),\theta \left(\eta \right)$, and $\varphi \left(\eta \right)$ are plotted against the local inertial coefficient $Fr$ in Figure 14, Figure 15 and Figure 16. The results show that $f^{\prime }\left(\eta \right)$ falls, but $\theta \left(\eta \right)$ and $\varphi \left(\eta \right)$ grow. This happens, physically, because an increase in $Fr$ basically causes an increment in the drag coefficient and a decrement in the porosity of the porous medium, which compels the fluid velocity to decrease. The data plotted in Figure 17 and Figure 18 reveal the physical impact of $Pr$ on the considered unknown quantities. The sketches show that, by strengthening $Pr$, the velocity and temperature of the fluid fall. From a physical point of view, the occurring phenomenon is true because, as $Pr$ increases, the viscous force is enhanced, and the thermal conductance of the fluid drops, causing are duction in the velocity and temperature; the opposite scenario is observed in the concentration profile. By increasing the values of $Pr$, the velocity boundary thickness and thermal boundary layer thickness gradually decrease. The Figure 19, Figure 20 and Figure 21 illustrate the impact of the chemical reaction parameter $R$ on the velocity profile, the temperature field, and mass concentration when the rest of the parameters are kept constant. It is seen that, as $R$ is augmented, a reduction occurs in $\varphi \left(\eta \right)$. The graphs show that, as $R$ is increased, the velocity decreases, the temperature increases, and mass concentration is decreased. The increase in temperature and the decrease in mass concentration are due to heat generated from the chemical reaction occurring in the fluid flow domain. All the graphs are plotted at theangle $\alpha =\frac{\pi }{6}$ (inclined exponentially stretching plate), and asymptotic behavior is observed in all the graphs satisfying the molded conditions.

5.2. Impact of the Involved Parameters on Skin Friction Coefficient $f^{″}\left(0\right)$, Heat Transfer Rate ${\theta }^{\prime }\left(0\right),$ and Mass Transfer Rate ${\varphi }^{\prime }\left(0\right)$

Table 1. Comparison of the $-{\theta }^{\prime }\left(0\right)$ values for several values of the Prandtl number for Newtonian fluid when $Ri=0, N=0, Sc=0, \beta =0, L=0, K=0, K^{*}=0, Fr=0, M=0, R=0$.

Pr	Magyari and Keller [32]	Present Study
1.0	0.9547	0.9551
3	1.8691	1.8121
5	2.5001	2.5577
10	3.6604	3.6868

shows the comparison of the present numerical results for the rate of heat transfer with the previously published results. From a close observation of the numerical solutions from both the previously documented and the current results, we conclude that there is good agreement between them, which indicates the validation of the present study.

Table 2. Consequences of $N$ on (a) skin friction coefficient $f^{″}\left(0\right)$, (b) heat transfer rate $-{\theta }^{\prime }\left(0\right)$, and (c) mass transfer rate $-{\varphi }^{\prime }\left(0\right)$, when $Ri=5.0, \beta =1.0, K=0.1, L=0.1, \alpha =\frac{\pi }{6}, K^{*}= 8.0, Fr=5.0, M=0.1, R=0.1$, Pr = 7.0, $Sc$ = 0.1.

$\mathit{N}$	${\mathit{f}}^{″}\left(0\right)$	$-{\mathit{\theta }}^{\prime }\left(0\right)$	$-{\mathit{\varphi }}^{\prime }\left(0\right)$
0.1	1.604773055747342	2.893024633389970	0.380437969348776
2.0	0.752124574958732	3.139862909322790	0.435468451438860
5.0	−0.676048985262494	3.402040713252884	0.489783407856660
7.0	−1.460502296892100	3.524921340615097	0.516004890188933

and

Table 3. Consequences of $Sc$ on (a) skin friction coefficient $f^{″}\left(0\right)$, (b) heat transfer rate $-{\theta }^{\prime }\left(0\right)$, and (c) mass transfer rate $-{\varphi }^{\prime }\left(0\right)$, when $Ri=5.0, \beta =1.0, K=0.1, L=0.1, \alpha =\frac{\pi }{6}, K^{*} =8.0, Fr=5.0, M=0.1, R=0.1$, Pr = 7.0, N = 0.1.

$\mathit{S}\mathit{c}$	${\mathit{f}}^{″}\left(0\right)$	$-{\mathit{\theta }}^{\prime }\left(0\right)$	$-{\mathit{\varphi }}^{\prime }\left(0\right)$
0.1	1.604773055747342	2.893024633389970	0.380437969348776
2.0	1.615378497071464	2.886105078686632	1.365657744516334
5.0	1.624426389709029	2.881459769245696	2.383773750441613
7.0	1.628041039070256	2.879970770253318	2.899471644583231

depict the numerical solutions of the skin friction coefficient, the rate of heat transfer, and the rate of mass transfer under the impact of the buoyancy ratio parameter $N$ and the Schmidt number $Sc$, respectively. From the tabular results displayed in Table 2, we can observe that, as $N$ is augmented, the $f^{″}\left(0\right)$ decreases, but there is a visible reduction in $-{\theta }^{\prime }\left(0\right)$ and $-{\varphi }^{\prime }\left(0\right)$. Table 3 illustrates the physical attitude of $f^{″}\left(0\right),-{\theta }^{\prime }\left(0\right)$ and $-{\varphi }^{\prime }\left(0\right)$ under the impact of the Schmidt number. The results indicate that, with the augmentation of Sc, there is an attenuation in $-{\theta }^{\prime }\left(0\right)$ and $-{\varphi }^{\prime }\left(0\right)$, but a rise is noted in $f^{″}\left(0\right),-{\theta }^{\prime }\left(0\right)$. These results are computed exactly at the surface.

6. Conclusions

The purpose of this research was to look at the physical effects of chemical reactions and the Lorentz force on electrically conducting third-grade fluid flow and heat transfer, with a focus on the Darcy–Forchheimer relationship along an inclined exponentially stretched sheet buried in a porous medium. The suggested process was expressed mathematically as partial differential equations, which were then simplified to ordinary differential equations using similarity variables for comparable solutions. The MATLAB built-in numerical solver bvp4c was used to solve the governing flow problem numerically. Graphs and tables were used to display the results of the material properties. The following is a summary of the main findings:

• The velocity $f^{\prime }\left(\eta \right)$ rises with the rising value of $Ri$ and $\beta$, but the reverse trend is noted for increasing values of $M, K^{*}, Fr, R$, and $Pr$.
• The temperature field $\theta$ is increased as $M, K^{*}, R$, and $Fr$ are augmented, but the reverse attitude is observed for an increase in the values of $Ri, \beta$, and $Pr$.
• The concentration profile is enlarged as $K^{*}, M$, and $Fr$ rise, but it falls as $Ri, \beta$, and $R$ are elevated.
• The skin friction coefficient $f^{″}\left(0\right)$ gets stronger owing to a rise in the values of $Sc$; the opposite trend is observed with the increasing values of $N$.
• The rate of heat transfer ${\theta }^{\prime }\left(0\right)$ expands as $N$ rises but reduces with increasing magnitudes of $Sc.$
• The rate of mass transfer $-{\varphi }^{\prime }\left(0\right)$ increases as $N$ and $Sc$ are augmented.
• It is concluded that the plotted results satisfy the given set of boundary conditions, which is evidence of the accuracy of the numerical solutions for $f^{\prime }\left(\eta \right),\theta \left(\eta \right)$, and $\varphi \left(\eta \right)$.
• The tabular results are presented for $f^{″}\left(0\right),-{\theta }^{\prime }\left(0\right)$, and $-{\varphi }^{\prime }\left(0\right)$ exactly at the surface.
• All the numerical results at the inclined exponentially stretching plate fixed at the angle of inclination $\alpha =\pi /6$ were computed.
• All the numerical results presented in the graphs plotted in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21 satisfy the given boundary conditions asymptotically; therefore, the numerical results given in tabular form are accurate.
• The current results are compared with the available results in the existing literature for this special case, and there is good agreement between them that shows validation of the present study.
• In the future, the study will be extended to third-grade nanofluid and hybrid nanofluid, with the inclusion of different flow features and the physical effects of the different fluid characteristics over the exponentially inclined stretching sheet embedded in a porous medium with different flow conditions.