Dimensionless Temperature Profiles Research Articles

Investigation of Soret and Dufour Effects on Chemically Reacting Free Convective Fluid Flowing Over a Vertical Plate Along with Viscous Dissipation Using the Laplace Adomian Decomposition Method

This study comprehensively analyses heat and mass transfer phenomena in a chemically reacting free convective fluid flow along a vertically moving plate. The flow is influenced by thermo-diffusion, diffusion-thermo, and viscous dissipation effects. To simplify the analysis, scaling group analysis and appropriate similarity transformations are used to transform the governing equations into nonlinear ordinary differential equations. These equations are then solved using a combination of Laplace transform and the Adomian decomposition method. The study conducts a parametric investigation to explore the impact of various control parameters on the dimensionless velocity, temperature, and concentration profiles. The parameters considered include the Prandtl number, Schmidt number, Eckert number, chemical reaction parameter, Soret parameter, Dufour parameter, solutal Grashof number, and thermal Grashof number. These parameters are depicted graphically and analysed quantitatively. The results reveal that an increase in the Schmidt number leads to a decrease in velocity and concentration profiles while temperature varies monotonically. Elevating the Eckert number enhances velocity and temperature profiles, with a slight decrease in concentration profiles. A rise in the Prandtl number decreases the temperature profile, with minimal effects on velocity and concentration profiles. Increasing the solutal Grashof number decreases temperature and concentration profiles, whereas the thermal Grashof number is directly proportional to the velocity profile. An increase in the Dufour parameter boosts velocity and temperature profiles while reducing the concentration profile. The presence of the Soret parameter increases velocity and concentration profiles but decreases the temperature profile. This study aims to enhance comprehension of the complex interactions within flow characteristics, providing valuable insights into the fundamental mechanisms of such systems. It also highlights their potential applications in various engineering and industrial processes.

Deduction of the Dimensionless Groups and Type Curves of Temperature Profiles in Two-Layer Soils with Water Flow at Depth

In the common hydrogeologic scenarios of horizontal groundwater flow and a water table below the surface, the steady-state 2D thermal field resulting from the coupling between water flow and heat flow and transport gives rise to a vertical temperature profile that develops progressively over a finite extent of the domain. Beyond this region, the temperature profiles are linear and independent of horizontal position. Such profiles are related to the groundwater velocity so they can be usefully used to estimate this velocity in the form of an inverse problem. By non-dimensionalization of the governing equations and boundary conditions, this manuscript formally derives the precise dimensionless groups governing the main unknowns of the problem, namely, (i) extent of the profile development region, (ii) time required for the steady-state temperature profile solution to be reached and (iii) the temperature–depth profiles themselves at each horizontal position of the development region. After verifying the mathematical dependencies of these unknowns on the deduced dimensionless groups, and by means of a large number of accurate numerical simulations, the type curves related to the horizontal extension of the development of the steady-state profiles, the characteristic time to develop such profiles and the dimensionless vertical temperature profiles inside the characteristic region are derived. These universal graphs can be used for the estimation of groundwater horizontal velocities from temperature–depth measurements.

Semianalytical investigation on heat transfer in porous fins with temperature‐dependent thermal conductivity via the homotopy perturbation Sumudu transform approach

AbstractA unique investigation has been undertaken to analyze the heat transmission by convective and radiative mechanisms in a fully saturated penetrable fin of a longitudinal structure positioned on a leaning surface. This study introduces the fusion of the realms of Homotopy perturbation and Sumudu transform techniques to address a previously unexplored problem involving a moving fin with temperature‐dependent thermal conductivity. In prior research papers, the Homotopy Perturbation Sumudu Transform Method (HPSTM) was utilized to obtain analytical solutions for fins featuring temperature‐dependent thermal conductivity. However, in our current study, we employ the HPSTM to tackle a novel problem involving a moving porous fin. This fin exhibits temperature‐dependent thermal conductivity and is subjected to convection and radiation effects. Through a comparison with numerical results, the present study has validated the dependability of its findings. The dimensionless temperature profile has been investigated by studying its relationship with several parameters. Here we observed that when the Peclet number is augmented by 400%, there is a corresponding 1.11% increase in thermal outline at the fin's extremity. Enhancing the value of radiation parameter by 400% declines the temperature of the fin tip by 14.079%. This study encourages the application of the Homotopy perturbation Sumudu transform technique in more complex fin problems.

Analytical approach for micropolar fluid flow in a channel with porous walls

This article aims to analyze and reach new physical achievements in a well-known nonlinear problem of two-dimensional, steady, and steady micropolar fluid flow between two parallel permeable porous walls. Micropolar fluids are considered a class of non-Newtonian fluids. The nonlinear governing equations of this problem were analyzed using two analytical methods. The achievements and novelty of this article, apart from showing the application of two mathematical methods in solving non-linear equations, are more about achieving new results that are obtained from solving the governing equations of this article. From the results, it is obvious that with the raise of the coupling parameter, other dimensionless parameters increase, except for the dimensionless microrotation. Also, when the spin gradient viscosity parameter increases, all the average values of the mentioned dimensionless parameters increase. From the problem's solution, it is clear that when the micro-inertia density ratio increases across the channel, all the average values of the dimensionless parameters decrease except for the dimensionless flow function characteristics. Also, when the Peclet number for the heat diffusion parameter is changed, it only affects the dimensionless temperature profile. When the Peclet number for the mass diffusion parameter is changed, it only affects the dimensionless concentration profile.

Application of Bernoulli wavelet method on triple‐diffusive convection in Jeffery–Hamel flow

AbstractIn this article, triple‐diffusive convection in the Jeffery–Hamel flow of viscous fluid is studied. The Jeffery–Hamel flow depends upon the radial component of velocity, while the peripheral velocity is zero. The problem has been articulated as nonlinear coupled partial differential equations (PDEs) together with the pertinent boundary conditions. The reduction of the nonlinear coupled PDEs into new nonlinear coupled ordinary differential equations is achieved via a collection of appropriate transformations, which is solved using the Bernoulli wavelet method. The obtained results were compared with the numerical results, confirming the good agreement between the present and previous numerical results. Dimensionless velocity, temperature, and concentration profiles are discussed for the relevant factors involved. Furthermore, skin friction, heat transfer, and solute diffusions are calculated.

Numerical investigation of MHD flow of hyperbolic tangent nanofluid over a non-linear stretching sheet

This research investigates two-dimensional MHD incompressible boundary layer Hyperbolic Tangent nanofluid flow across a non-linear stretching plate. Similarity transformations are employed to convert the governing non-linear partial differential equations (PDEs) into coupled non-linear ordinary differential equations (ODEs). The MATLAB built-in routine bvp4c has been used for finding the numerical solutions of the dimensionless velocity, temperature, and concentration profiles. The current findings are validated with already published results. The influence of some important parameters on the velocity, temperature, and concentration profiles are displayed through graphs and tables. It is observed that for increasing values of magnetic parameter M and hyperbolic Tangent parameter We, the boundary layer thickness of the velocity profile decreases while it increases for the temperature profile.

Investigation of radiative MHD Carreau nanofluid over an inclined stretching cylinder with activating energy, Wu’s slip, and convective heating

The influence of radiation, convective heating, second-order velocity slip condition, and activating energy on the magnetohydrodynamic (MHD) Carreau nanofluid and heat/mass transfer are depicted in this study. A mathematical model of the Carreau nanofluid is carried out with the effects of thermophoresis and Brownian motion. The governing equations are transformed into a system of nonlinear coupled differential equations. Solutions of the transformed equations are derived by the Matlab BVP4c solver. Efficiency of the solver is verified by the comparison with the previous literature. Moreover, the impacts of the dimensionless velocity, temperature, and nanoparticle concentration profile are elucidated graphically.

Heat transfer in MHD flow of Carreau ternary-hybrid nanofluid over a curved surface stretched exponentially

This investigation aims to study Magnetohydrodynamics (MHD)two-dimensional incompressible boundary layer performing non-Newtonian Carreau ternary-hybrid nanofluid flow with heat transfer through an exponential stretching curved surface. The ternary-hybrid nanofluid has been synthesized with titanium oxide, aluminum oxide, and silver dispersionin the base fluid water. TheNavier Stokes equation and Carreau ternary-hybrid nanofluid model govern the partial differential equations (PDEs), and appropriate similarity transformations are utilized to transfer these PDEs into ordinary differential equations (ODEs). The effects of the pertinent parameters on the dimensionless velocity and temperature profiles are analyzed withfigures. This study provides new insights and solutions to previously unsolved problems related to heat transfer in the MHD flow of a Carreau Ternary-Hybrid Nanofluid over a curved surface stretched exponentially, or it could contribute to the existing knowledge and literature by refining existing models or methods. The surface drag force and Nusselt numbers are studied for the different values of the governing parameters throughgraphs. It is demonstrated that the heat transfer rate and skin friction increase from base fluid to mono, hybrid, and ternary nanofluids. Both heat transfer rate and skin friction increase with the addition of nanoparticles.

Mathematical modeling and analysis of immiscible metallic based nanofluid flow in a microchannel with non-spherical nanoparticles

The non-spherical shaped nanoparticles have grown in popularity due to their ability to alter the nanofluid thermophysical properties. This nanofluid plays a beneficial role in biomedical imaging, cancer and tumor therapy. In the current study the effect of pressure gradient, transverse magnetic field and electroosmosis on the flow of two immiscible metallic based nanofluids (blood-based Au/Ag Casson nanofluid and water-based Au/Ag Newtonian nanofluid) of different viscosities in separate regions in a microchannel is addressed. The combined impact of electromagnetohydrodynamics (EMHD), modified Darcy’s law, Joule heating and Hall current have been considered into account. The mathematical modeling has been developed with the help of modified Navier Stoke’s equations, and simplified the complex problem using appropriate assumptions and non-dimensional parameters. The complementary function and particular integral concepts have been used to obtain the exact solutions of a non-dimensionalized system of coupled differential equations along with non-dimensionalized conditions. The dimensionless velocity, temperature and concentration profiles of nanofluids have been computed. The influence of numerous relevant fluid transport parameters has been graphically depicted. It is noticed that the magnetic strength declines the nanofluid velocity in both regions. Further, the radiation parameter exhibits an inverse effect on the rate of heat transfer in nanofluids. The heat transfer rate enhanced by 56% and 2% with rising volume fraction [0.01 to 0.02] and zeta potential [0.1 to 0.2], respectively. With increasing strength of magnetic field [2 to 3], thereby 31% increment in shear stress distribution. The results of the current study can be used to enhance the performance of various processes such as surface-enhanced Raman spectroscopy, drug delivery, sensors, and electrochemistry.

Dual Solutions and Stability Analysis of Cu-H2O-Casson Nanofluid Convection past a Heated Stretching/Shrinking Slippery Sheet in a Porous Medium

In this study, we examined the impact of Cu-H2O nanoparticles on two-dimensional Casson nanofluid flows past permeable stretching/shrinking sheet embedded in a Darcy-Forchheimer porous medium in the presence of slipperiness of surface, suction/injection, viscous dissipation, and convective heating. Using some realistic assumptions and appropriate similarity transformations, the governing nonlinear partial differential equations were formulated and transformed into a system of nonlinear ordinary differential equations and then numerically solved by using the shooting technique. Numerical results are displayed for dimensionless fluid velocity and temperature profiles, skin friction, and the local Nusselt number. The impacts of different governing physical parameters on these quantities are presented and discussed using graphs, tables, and a chart. For the specific range of shrinking sheet, the result shows that dual solutions exist, and temporal stability analysis is performed by introducing small disturbances to determine the stable solutions. It is detected that the upper branch solution is hydrodynamically stable and substantially realistic; however, the lower branch solution is unstable and physically unachievable. The fluid flow stability is obtained by enhancing the suction, surface slipperiness, and viscous dissipation parameters. However, augmenting the values of the Casson factor, Cu-H2O nanoparticle volume fraction, porous medium, porous medium inertia, and convective heating parameters increases the blow-up stability of the fluid flow. The rate of heat transfer enhances with the increment in the Casson factor, porous medium, porous medium inertia, suction, velocity ratio, nanoparticle volume fraction, and convective heating parameters, whereas it reduces as the slipperiness of the surface and viscous dissipation parameters rise. Increment of Cu-H2O nanoparticle volume fraction into the Casson fluid boosts the heat transfer enhancement rate higher for the shrinking sheet surface.

Convective micropolar fluid over inclined surface with thermal radiation and velocity slip condition effects: Duality and stability

This research examines the flow of 2D micropolar fluid across an inclined linear shrinking/stretching surface with suction, convective, slip and thermal radiation impact. The governing partial differential equations (PDEs) are transformed into an ordinary differential equations (ODEs) system using a linear similarity transformation, and the bvp4c technique is then used to solve the ODEs system in MATLAB. Comparisons of both solutions are presented. The effects of regulating flow parameters on dimensionless temperature profiles, angular velocity, velocity, skin friction, wall couple stress and heat transfer are shown graphically. Stable velocity profiles are inversely related to velocity slip and material parameter. For a given suction and stretching/shrinking parameter value, graphs show double solutions. The couple stress and skin friction in the first solution, drop then increase as the material parameter K increases. Critical points of the stretching/shrinking parameter, [Formula: see text] where both solutions occur, are indicated by the symbols [Formula: see text] (i.e., [Formula: see text], and that there is no solution when [Formula: see text]. A table summarizes the stability study’s results and highlights.

Numerical simulation of mass and heat transport phenomena of hydromagnetic flow of Casson fluid with sinusoidal boundary conditions

AbstractA magnetic field might interact with its surroundings, influencing chemical and physical processes in materials processing, heat exchangers, and other scientific study. Therefore, a computational study of non‐Newtonian (Casson) free convective MHD unsteady fluid flow has been highlighted in this article with mass and heat transit property through a vertical infinite porous plate. A sinusoidal boundary conditions as well as chemical reaction and thermal radiation have been considered. Using a collection of nondimensional variables, the flow related equations are also turned into nondimensional form. The EFDM algorithm is employed in order to arrive at a numerical solution via Compaq Visual Fortran. The reliability of the numerical solution has been confirmed using stability testing and convergence analysis. The whole system is convergent when the values of Prandtl number and Lewis number are greater than or equals to 0.075 and 0.025, respectively. A visual depiction of the impact of the pertinent factors on dimensionless velocity, temperature, and concentration profiles are displayed through graphical representation as well as with tabular representation. It has been inspected that when the magnetic component is regarded, it greatly affects the heat transfer factors of Casson nanofluid and the heat also rises when Eckert number, heat source and radiation parameter accelerate. It is also found that the Sherwood number is increased as the impact of chemical reaction parameter and the Lewis number, also the skin friction is decreased as the influence of porosity term got accelerated. The comparison of the current findings with the data that were previously published serves as the final stage in validating the present study.

Assessment of transient thermal distribution in a moving porous plate with temperature-dependent internal heat generation using Levenberg–Marquardt backpropagation neural network

ABSTRACT In the present paper, the transient temperature distribution within a moving porous plate is scrutinized by considering the internal heat generation, convection and radiation condition. Furthermore, an effective artificial intelligence technique is proposed for discussing the transient thermal behavior in a moving plate by employing the Levenberg-Marquardt backpropagation network. Using suitable non-dimensional terms, the governing partial differential equation (PDE) is transformed into its non-dimensional form. The transformed equation is then solved using the finite difference method (FDM) and the obtained thermal values (datasets) are used as target values in the process of neural network mapping. This network uses the training, testing, and validation processes to determine the solution of several scenarios established at various process times. The behavior of dimensionless transient temperature profile has been inspected for diverse values of non-dimensional parameters such as heat generation number, radiative parameter, porosity parameter, and conductive parameter with the assistance of graphs. The significant outcomes elucidate that, the transient temperature profile of the porous plate increases with the change in time. Further, it enriches for rise in the magnitude of Peclet number and heat generation number while it drops for enhanced values of porosity parameter.

Heat transfer analysis of a convective radiative porous moving longitudinal fin exposed to magnetic field by Adomian decomposition sumudu transform method

A massive number of engineering disciplines are concerned with the enhancement of heat exchange. In this regard, the present article focuses on thermal analysis of natural convection and radiation of a fractional ordered moving fin problem exposed to the magnetic field and also for the first time Adomian decomposition sumudu transform method (ADSTM) has been introduced to solve the energy balance equation comprising the Peclet number. The governing equations along with the corresponding boundary conditions were solved with the aid of ADSTM. The accuracy of the present work has been validated through comparison with the numerical results. ADSTM results are in good agreement with the existing numerical results. The effect of various parameters on the dimensionless temperature profile for different fractional values of is discussed in graphical records. We found that the temperature distribution along the tip of the fin increases by 1.53% as the Peclet number increases by 400%, and the thermal gradient along the fin tip decreases by 15.09% as the Hartman number increases by 400%.

Effect of Magnetic Field on the Developing Thermal Field in a Duct Filled with Porous Media under Local Thermal Non-Equilibrium with a Nonlinear Flow Model

In this article, the numerical study of the influence of a magnetic field at the laminar forced convection in a thermally developing region coming under the influence of local thermal non-equilibrium (LTNE) of parallel plate channels completely immersed in the porous material is investigated. Constant wall heat flux boundary conditions are applied to the walls of the channel. In the nonlinear flow model, the Darcy-Brinkman-Forchheimer equation governs the flow field in the porous region, which is assumed to be unidirectional. The system is defined by certain well-known parameters, these being Darcy number (Da), thermal conductivity ratio (κ), Forchheimer number (F), Hartmann number (M), and Biot number (Bi). Numerical solutions have been obtained by applying a successive accelerated replacement (SAR) scheme. Exact solutions for the dimensionless temperature and the fully developed Nusselt number in the absence of the Forchheimer number (F = 0), for the fully developed thermal field, are obtained for the linear flow model, the Darcy-Brinkman model. Plots are given for the dimensionless temperature profiles in the fluid and solid phases, wall temperature, as well as the local Nusselt number at the parallel plate channel, which has been displayed. The effect of the magnetic field and the thermal conductivity ratio has a significant effect on the local Nusselt number.

Dual solutions on three-dimensional nanofluid flow and heat transfer over a permeable non-linearly shrinking surface with second-order velocity slips

PurposeThis paper aims to present the steady dual solutions on three-dimensional flow and heat transfer of nanofluid over a permeable non-linearly shrinking surface with two-order velocity slips conditions. Boundary layer assumption is considered in the mathematical modelling. Authors comprehend from previous studies and papers that the shrinking surfaces are extremely important in current engineering and environmental systems.Design/methodology/approachUsing appropriate similarity variables, the full partial differential equations (PDF) are modified into a specific set of ordinary (similar) differential equations (ODE). The resulting non-linear ordinary differential system is then solved both analytically for some particular cases and numerically for the general case using the function bvp4c from MATLAB for characteristic values of the parameters which govern the equations. The transformed mathematical model is analysed using the bvp4c procedure. Based on the given assumptions, this study is able to produce multiple solutions of the problem.FindingsThe ordinary (similarity) differential equations have two branches solutions, upper and lower branch solutions, given some interval of shrinking and velocity slip parameters. The authors consider here a temporal stability analysis, as they want to establish which of the solutions are stable and which are not. In a distinct paragraph, the authors discuss in detail and present in a graphical manner the effects of shrinking and second-order slip flow model on the skin friction coefficient, surface wall heat flux and dimensionless velocity and temperature profiles. The analysis reveals that the second order slip has a big influence on the flow and heat transfer characteristics.Originality/valueThe present discoveries are unique and truly new for the research of the three-dimensional stretching/shrinking forced convection flow and heat transfer nanofluids. The nanofluid is a water-based nanofluid (H2O), which contains one type of nanoparticles, namely, copper (Cu). Of course, the analysis can be further extended considering other types of nanoparticles such as alumina (Al2O3). The authors assume that the thermal equilibrium is reached for the base fluid together with the suspended nanoparticles and that the nanoparticles are uniform in dimension and form.

Thermal investigation into the Oldroyd-B hybrid nanofluid with the slip and Newtonian heating effect: Atangana–Baleanu fractional simulation

The significance of thermal conductivity, convection, and heat transportation of hybrid nanofluids (HNFs) based on different nanoparticles has enhanced an integral part in numerous industrial and natural processes. In this article, a fractionalized Oldroyd-B HNF along with other significant effects, such as Newtonian heating, constant concentration, and the wall slip condition on temperature close to an infinitely vertical flat plate, is examined. Aluminum oxide (Al2O3) and ferro-ferric oxide (Fe3O4) are the supposed nanoparticles, and water (H2O) and sodium alginate (C6H9NaO7) serve as the base fluids. For generalized memory effects, an innovative fractional model is developed based on the recently proposed Atangana–Baleanu time-fractional (AB) derivative through generalized Fourier and Fick’s law. This Laplace transform technique is used to solve the fractional governing equations of dimensionless temperature, velocity, and concentration profiles. The physical effects of diverse flow parameters are discussed and exhibited graphically by Mathcad software. We have considered 0.15≤α≤0.85,2≤Pr⁡≤9,5≤Gr≤20,0.2≤ϕ1,ϕ2≤0.8,3.5≤Gm≤8, 0.1≤ Sc ≤0.8, and 0.3≤λ1,λ2≤1.7. Moreover, for validation of our present results, some limiting models, such as classical Maxwell and Newtonian fluid models, are recovered from the fractional Oldroyd-B fluid model. Furthermore, comparing the results between Oldroyd-B, Maxwell, and viscous fluid models for both classical and fractional cases, Stehfest and Tzou numerical methods are also employed to secure the validity of our solutions. Moreover, it is visualized that for a short time, temperature and momentum profiles are decayed for larger values of α, and this effect is reversed for a long time. Furthermore, the energy and velocity profiles are higher for water-based HNFs than those for the sodium alginate-based HNF.

Influence of chemical reaction and variable mass diffusivity on non-Newtonian fluid flow due to a rough stretching sheet with magnetic field and Cattaneo-Christov fluxes

A non-Newtonian Williamson fluid flow due to a stretching sheet with radiation, magnetic field, and viscous dissipation effects is described using variable conductivity and variable diffusivity. The Cattaneo-Christov model is used to correctly compute the physical properties of a heat and mass flux model. Both the chemical reaction phenomenon and the slip velocity have an impact on the heat and mass mechanism. The physical problem is represented mathematically as a nonlinear coupled differential system. After that, the shooting method is used to solve the mathematical model numerically. To gain a better understanding of the behavior of governing emergent factors on dimensionless velocity, concentration, and temperature profiles, physical interpretations are created and discussed utilizing graphical and tabular representations. The results show that the Sherwood number and the Nusselt number are both decreased by the magnetic, viscosity, and slip velocity parameters. Also, according to the findings it has been observed that the concentration outlines enhances for the magnetic number, the viscosity parameter, and the slip velocity parameter, but they dwindle for expanding reaction rate values. Finally, after confirmation of our numerical results, the theoretical results show good agreement with previously published work.

Thermal performance of a longitudinal fin under the influence of magnetic field using Sumudu transform method with pade approximant (STM‐PA)

AbstractTemperature distribution, and efficiency of a rectangular profiled longitudinal fin are examined in this investigation with the impact of the magnetic field. By exploiting appropriate non‐dimensional terms, the heat transfer equation incorporating temperature‐dependent thermal conductivity, heat transfer coefficient, and Maxwell expression for the effect of the magnetic field yield a dimensionless nonlinear ordinary differential equation (ODE) with corresponding boundary conditions (BCs). Sumudu transform method with Pade approximant (STM‐PA) has been employed to obtain an analytical solution for the temperature profile of a longitudinal rectangular fin subjected to a uniform magnetic field under multi‐boiling heat transfer. The STM‐PA results are compared to the Runge‐Kutta Fehlberg's fourth‐fifth (RKF‐45) order technique for computational verification and are observed to be in good accordance. The behavior of dimensionless temperature profile has been explicated graphically for diverse values of non‐dimensional parameters such as thermal conductivity parameter, Hartmann number, and thermogeometric parameter. The results of this study show that as the thermal conductivity parameter enriches, the temperature profile of the longitudinal fin improves, whereas it declines for the Hartmann number and thermogeometric parameter. Under multi‐boiling heat transference, fin efficiency varies significantly depending on the impact of pertinent variables.

Numerical Investigation of MWCNT and SWCNT Fluid Flow along with the Activation Energy Effects over Quartic Auto Catalytic Endothermic and Exothermic Chemical Reactions

A mathematical model is created to analyze the impact of Thompson and Troian slip boundaries over a contracting/expanding surface sustaining nanofluid-containing carbon nanotubes along a stagnation point flow. Both multi-wall (MWCNTs) and single-wall (SWCNTs) carbon nanotubes are taken into consideration, with water serving as the base liquid. The flow is obtained due to the stretching or contracting of the surface. The thermal radiation, activation energy, buoyancy impacts, and chemical processes called quartic autocatalysis are additionally added to the original mathematical model. The MATLAB-constructed bvp4c function involving the three-stage Lobatto IIIa formula for the numerical results of dimensionless velocity, concentration, and temperature profiles are used. By contrasting it against a published paper in this limited instance, it is determined whether the suggested mathematical model is legitimate. In this sense, a remarkable consensus is achieved. Graphical representations are used to depict the behavior of many non-dimensional flow variables, such as the slip velocity parameter, the inertia coefficient, the porosity parameter, and the solid volume fraction. Surface drag force computations are reported to examine the effects at the permeable stretching surface. It has been shown that increasing the slip velocity factor increases the fluid streaming velocity while decreasing the surface drag force. If the endothermic/exothermic coefficient increases, the local thermal transfer efficiency falls. For nanofluids, the changing viscosity factor increases axial velocity while decreasing temperature distribution. Additionally, the solid volumetric fraction improves the temperature distributions by lowering the concentration profile and speed.