Analysis of mixed convecto-magnetic Buongiorno nanofluid flow in a non-Darcy porous medium

Analysis of nanofluid flow phenomena through a stretching sheet in a non-Darcy porous medium has great importance in enhancing transport processes in energy systems like hybrid fuel cell technology, glass blowing, etc. The current study looks at two-dimensional flow over a convectively heated linear stretching sheet in presence of a magnetic field. The boundary layer flow caused by a sheet that is linearly stretched has been investigated numerically. Employing the appropriate similarity transformation, the governing partial differential equations, which define the flow regime, are converted into a set of ordinary differential equations. The mathematical calculations are carried out using a finite difference algorithm with the aid of Newton's linearization method, which allows us to handle non-linear terms very smoothly. The outcomes of this work have been demonstrated are based on eight parameters, such as Prandtl number , Lewis number , Brownian motion parameter , thermophoresis parameter , magnetic parameter , porosity parameter , Frochemier number , and convection Biot number . The impacts of aforementioned parameters on thermal and concentration boundary layers are depicted graphically. It has been found that the heat transfer in the boundary layer rises as enhances, and deposits aggravating particle away from the fluid region and increases the volume percentage of nanoparticles. Moreover, it has been found that as the number rises, the concentration profiles become steeper and the species border layer becomes thinner. Furthermore, it has been encountered that the reduced Sherwoord number rises as Le number increases for but decreases as number increases for .


Nomenclature b
favourable constant related to linear stretching Bi = h(v 1 /β)

D c
Lewis number Nusselt number Nur 1 = Re

Introduction
The manipulation of NF's involves immersing of NP's of an average size of less than 10 2 nm in conventional HT fluids like ethylene glycol, oil, and water.Thermal transmission of water, oil, and ethylene glycol with NP's simulates a significant portion of the HT between the HMT medium and surface.The fluids' heat transmission behaviour has been significantly improved due to the increase in the TC.Many engineering applications require high heat transfer efficiency.The TC of NF varies linearly with temperature and variable in many engineering and industrial uses.The mixture of water and nanomaterials is known as a nanofluid.Thermal conductivity and convective heat movement are improved when a surfactant is added to a base fluid.A modern and efficient method of heat transmission is nanofluid technology.Nanofluids have applications in cancer therapy, surgery based on heat treatment applications, as well as solar energy applications.Better oils and lubricants can be made using nanofluid technology for practical uses [1].Over the last few decades, several researchers have been working on fluid to increase efficiency for various thermal uses.By utilizing nanotechnology for HT, the word "NF" is used to discuss about new HT issues [2][3][4][5][6][7].Numerous authors have projected that NF's will be more effective at convective HT than basic fluids.Nanofluids' ability to transfer heat and mass has incorporated many practical uses for thermal and subatomic stratification [8][9][10][11][12][13].Patil et al. [14] studied an unsteady mixed convection NF flow in an exponentially stretching surface with a transverse magnetic field.Using homotopy analysis, [15] showed the entropy analysis and gyrotactic microorganisms of Buongiorno's NF within a pair of stretchable rotating discs.Nainaru et al. [16] showed that how HT characteristics impact on 3D MHD flows of a NF on stretching a surface using thermal radiation.They noticed that the variable TC increased along with the fluid temperature and movement.[17] studied the forced convection NF flow through a horizontal microchannel with cross-flow injection and slip velocity on sidewalls of the microchannel.Increasing the angle of injection up to 94 0 , they noted that both velocity slip and HT rate are increased.Suhail and Siddiqui [18] studied the natural convection NF flows through a vertical annular region numerically.They discussed the simulations for different heat flux values and nanoparticle concentrations.The convective boundary layer is crucial for supplying heat to the fluid via a particle surface with a limited HE because of its numerous and practical applications, particularly in various manufacturing processes like transpiration cooling, laser pulse heating etc. [19][20][21].Convective HT is crucial in a variety of devices, including thermal insulation, gas turbines, and nuclear power plants [22,23].Magnetohydrodynamics is the study of electrically conducting fluids, such as seawater, antioxidants, plasma, and metal compounds (MHD).The acronym MHD was created by Alfven [24].MHD is influenced by magnetic induction intensity.Numerous technical and commercial uses for the MHD fluid flow include crystal growth, nuclear freezing, MHD detectors, and energy generation, which has numerous uses in the medical and biopharmaceutical industries, including hyperglycemia, emergency treatment, radiation therapy, and many others.Using the Keller box method, [25] evaluated the melting HT and nonuniform heat source on magnetic NF flow past a porous cylinder.Jha and Malgwi [26] examined the effects of buoyant forces, pressure gradients, Hall and Ionslip currents on the mixed MHD convection flow of a viscous fluid through a microchannel placed vertically.Using explicit finite difference method, Haque [27] investigated the impacts of magnetic field and a heat sink on the micropolar fluid flows through a semi-infinite vertical porous plate.They discovered that the species concentration is influenced by both thermal diffusion and heat sink and that lighter particles concentrated in the fluid highly than heavier particles Alam et al. [28] investigated the biomagnetic flows of blood with gold nanoparticles over a stretched surface with suction, and velocity slip using bvp4c method.They took into consideration the principles of magneto hydrodynamics and ferrohydrodynamics.
Porous media is a surface with fluid-filled holes (voids).Reverse osmosis, heterogeneous catalysts, and evaporative freezing are among the industrial applications for porous materials [29][30][31].Researchers are engaged in analyzing the fluid flow processes in porous media for having numerous applications in production of scientific, biological, and industrial goods, including processing of crude oil, casting etc. [32][33][34][35][36]. Darcy established a theory in 1856 that asserts the pressure gradient is directly related to the volumetric flux of fluid across a medium.Practically, the Darcy term [37] is frequently applied to modelling and analysis issues involving flow-saturating porous regions.Only when the velocity is low and the porosity is low does Darcy's rule apply.This rule does not take into account inertia or boundary effects.The NDP space that integrates inertial and boundary impacts is the outcome of a modified version of the traditional Darcy's principle.As a consequence, Forchheimer [38] used a U 2 term in the momentum term to account for inertia.Non-Darcian variants are variations on the basic Darcy idea, which also takes into account inertial drag, vorticity dispersion, and their combined effects.It is claimed that the Darcy law fails in the presence of high velocities, highly porous media, and extremely high Re.This technique incorporates the Forchheimer expression into the square velocity term of the momentum equation to account for the inertia effects of pressure drop.Temperature shifts in HMT flows of DP and NDP media were studied by Pop, and Ingham [39] and Vafai [40].Using the ND solve method, [41] performed the DF 3D flows of NF over a rotating disc considering activation energy effects and increase and decrease in heat.Using the shooting method, Asma et al. [42] investigated the convective HT of DF 3D flows of NF through a rotating disc with Arrhenius activation energy.[43] used the bvp4c technique to study the impacts of DF permeable media, CC heat flux, and homogeneous-heterogeneous reactions on the melting HT between two parallel discs.Using the successive over-relaxation method, [44] performed the MHD convection flows of a hybrid NF past a stretching porous surface with viscous dissipation impacts.By homotopy analysis method, Kareem and Abdulhadi [45] looked into the axisymmetric MHD DF flows of 3 rd grade fluid over a stretching cylinder considering the CC effects.Use of homotopy analysis scheme has been made, Jawad et al. [46] to provide the entropy analysis of MHD radiative DF 3D Casson NF flow over a rotating disc taking into account the Arrhenius activation energy.Buongiorno's model was employed by Sheremet and Pop [47] to quantitatively analyze the NF flow in the cavity.In recent years, Bhatti et al. [48] studied the flow of an incompressible, 3-D, unsteady, axisymmetric squeezed Williamson NF between circular and rotating parallel plates with magnetic field, porosity, and suspension of motile gyrotactic microorganisms.They concluded that the axial and tangential velocity curves fall as the squeezing Re rises.In addition, they stated that the magnetic field weakens both axially and tangentially as the magnetic Re rises.
Despite the highly anticipated future of Buongiorno NF flow over a stretching surface saturated with a NDP medium, no study has been done till now.The majority of the work has been performed on DP medium for different types of NF flow.The aforementioned reasons encourage us to accomplish the current work.In the current manuscript, the non-Darcian NF flow and HT across a stretched surface in a NF-saturated porous media has been addressed.The aim is to study the influence of D%, Pr, Le, N b , N t , Fr 1 , K, and M of non-Darcy NF over a stretching surface with magnetohydrodynamics. Using the finite difference approach and Newton's linearization technique, the system of coupled nonlinear differential equations has been numerically solved.Graphical representations have been made of how embedded factors affect fluid velocity, temperature, and particle concentration.There is a good agreement between the findings of the current research and those of [49][50][51][52].This work is an extension of the study of Makinde and Aziz [51], Khan and Pop [52], Kuznetsov and Nield [53], and Nield and Kuznetsov [54].It is anticipated that the findings will complement earlier research in addition to offering useful information for applications.The findings in this study are particularly beneficial to many technical communities, including those that work with transpiration cooling, fabric cleaning, laser pulse heating, etc.The authors claim that no previous integrated effort of this nature has been made so far.

Geometry configuration
We have considered the steady and incompressible MHD nano fluid flow induced by linear movement of the stretching sheet with a linear velocity variation with distance x i.e.
u w = bx (1) In the normal direction, B 0 (strong magnetic field) is employed.The effects of viscous and ohmic discharge are both ignored.By selecting a small Re, the aspects of the induced magnetic field are disregarded.

Mathematical formulation
The work of Makinde and Aziz [51], Khan and Pop [52], Kuznetsov and Nield [53], and Nield and Kuznetsov [54] has been extended in the present work to take into account Buogiorno NF flow through a NDP medium, which is governed by the following conservation of mass, momentum, thermal energy, and volume fraction equations in 2-D Cartesian coordinates.
As per the studies of Makinde and Aziz [51], Khan and Pop [52], Kuznetsov and Nield [53], and Nield and Kuznetsov [54], the following boundary conditions have opted for.
The following dimensionless quantities are introduced, where The equation ( 2) can satisfied by the equation (10).Following is an analysis of the y direction momentum equation (normal to the sheet) by order of magnitude employing the standard BL approximations: Using the concept of equations (11)(12) in equation (4), we have Now, the equations (2-8) are thus reduced to the following set of equations (14-16) when the pressure gradient in the y direction is omitted (equation 13).
The boundary conditions (7-8) become, Due to the buoyancy effects caused by the gradients in nanoparticle concentration, there is no heat transfer when N b = 0. Here, we perceive that equation 14 has a closed-form solution, which is given by f = 1 − e η, along with the corresponding boundary conditions on f provided by expression (17).

Numerical calculations and validation
Using the boundary conditions (17)(18), the equations (14-16) have been solved by FDM.At the outset, the finite difference equation corresponding to the 3 rd order ordinary differential equation ( 14) is derived as follows: Using the substitution of in equation 14, we have Using the concept of finite difference form in equation 20, we obtain Since the preceding equation is non-linear, the accompanying generalization of Newton's linearization technique is employed to linearize the equation, which is then followed by the Thomas algorithm.In this procedure, the variable h 1 's solution at the (n + 1) th step is written as The linearized forms of equations (21-22) generate a system of tri-diagonal equations as follows: Where i = 0 and i = m are at the boundaries since the solution at i = 0 and i = m as The linear system of equation 23 is solved to determine Δh 1 n+1 i .Subsequently, the solution of h n+1 i is modified through equations (20)(21), where f is obtained by Tri diagonal matrix algorithm (TDMA) from the following equation Similarly non-linear 2 nd order ordinary differential equations (15)(16) are expressed as below: Equations (27)(28) are also solved numerically in the similar way.The numerical outcomes of Nur 1 are compared with studies of [49][50][51][52] in Table 1 to evaluate the accuracy of our numerical calculations.These asserted to circumstances where the temperature of the stretching sheet is fixed and thermophoresis effects are not present.The findings presented here to agree to three decimal places with those of [49][50][51][52] for various values of Pr listed in table 1.The findings in table 1 show a pattern resembling natural convection from a vertical plate in a regular fluid, where the Nur 1 rises as the Pr rises.
The constant temperature findings have been retrieved by putting Bi = ∞ in expression 17, yielding the condition θ 1 (0) = 1 (isothermal condition).Figures 2(a-b) compare our findings for the Nur 1 and Shr 1 to those reported by [52] at Le = 10, Pr = 10, Bi = ∞.The outcomes showed a flawless agreement between our findings and the results reported in [52] for every combination value of N t , N b used in our computations.From the figure 2(a), it has been demonstrated that the profile of Nur 1 sharply declines as N b is increased from 0.1-0.5 for fixed N t = 0.1.Even so, the Shr 1 is improved by a rise in N t (figure 2(b)), which is consistent with the findings of [52].The thermal boundary layer thickens as the Brownian motion expands, affecting a greater area of the fluid, which reduces the Nur 1 .

Results and discussion
Using graphical representations, illustrative numerical findings are presented in this section.Calculations have been performed for various values of physical factors, like Le, N b , N t , Bi, Pr, M, K, Fr 1 , andK.In the boundary layer region, Le number is a measure of the thermal diffusion rate, which contributes to the species diffusion rate.When Le is equal to 1, heat and species disseminate equally, and when Le is greater than 1.0, heat diffuses more quickly than species.When Le number increases, the species border layer become thin, and the temperature profiles become steeper.At N b = N t = 0.1, Pr = 5, Le = 5, M = 0.5, K = 0.3, Fr 1 = 0.5, the impact of Bi number on the distribution of temperature has been depicted in Figure 4(a).The Bi number links surface convection and solid conduction.Physically, the surface thermal resistance drastically decreases as the Bi number increases.The rise in convection raises the surface temperature.As the Bi number rises, the buoyancy force grows and the fluid transfers heat energy more quickly.The density of fluid changes with temperature due to thermal expandability, and the buoyancy force affects the movement.An inadequate convective thermal scenario is represented taking Bi = 0.1.Considering of high values of Bi number allows to achieve the constant wall temperature case θ 1 (0) = 1.A rise in Bi number is accompanied by a significant increase in the temperature profiles.The impact of Pr number on the distribution of temperature has been demonstrated in figure 4(b) at N b = N t = 0.1, Le = 5, Bi = 0.1, M = 0.5, K = 0.3, Fr 1 = 0.5.In general, as Pr numbers increase, the temperature profiles tend to decrease.In engineering and production procedures, the parameter Pr is plays a crucial role.The Pr number relates force diffusivity to warm diffusivity.Thermal diffusivity varies inversely with Pr.For higher values of Pr, the viscosity of the fluid becomes high, and for lower values of Pr, the fluid becomes less viscous.With an increase in Pr values, a substantial reduction in temperature profiles is seen.This pattern resembles the boundary layer movement of a free convective fluid [55].The physical meaning of the Pr number's findings is thus very consistent.

Analysis of velocity profile
At fixed value of N t , N b , Pr, Le, K, Fr 1 , Figures 5(a   the drag coefficient, higher values of Fr 1 reduce the velocity profiles.As a result, the drag coefficient rises as F r1 rises.As a response, the resistance force of the fluid rises, and the velocity decreases (figure 6).A rise in Fr 1 values of thickens of the thermal boundary layer makes difficult for the fluid to pass.The scientific justification for this fact is that raising the permeability values reduces the resistance of the porous medium, which tends to decrease the fluid velocity and improvements in the inertia coefficient (Fr 1 ) decrease the velocity due to a rise in the quadratic drag.When Le is equal to 1, heat and species disseminate equally, and when Le is greater than 1.0, heat diffuses more quickly than species.As the Le number rises, the concentration profiles become steeper and the species border layer becomes thinner.At Le = 5, Pr = 5, M = 0.5, K = 0.3, Fr 1 = 0.5, the distribution of concentration has been shown in figure 7(b) for various N t and N b .The nanoparticle volume fraction is greatly influenced by the parameter N t .The procedure for transporting particles reduce temperature due to the effect of temperature gradient, which is known as N t .The boundary layer's HT rises as N t rises, aggravating particle deposition away from the fluid area and, as a result, increasing the volume fraction of nanoparticles.The N b parameter gauges how much the suspended nanoparticles in the NF are moving randomly.Physically, an increase in the temperature causes the particles' energy to grow, which in turn increases their random motion and rapid collision, i.e. increases the Brownian motion.Due to spontaneous diffusivity, the nanoparticles reconfigure themselves and form a novel structure, which in turn enhances the NF's TC.The Brownian motion pushes the particles away from the fluid system and heats the fluid in the boundary layer.As the values of N t and N b increases, the strength of the Brownian motion and thermophoresis are found to have lowest impact on the distribution of concentration.4b), while the steepness of the curves increases.It is observed that the profile of Nur 1 increases as Pr increases, and decreases as N t increases.This is because for greater values of Pr, the convection process is dominant compared to conduction, which increases the HT.The higher thermophoretic force that propels the high thermal conductivity nanoparticles from the hot sheet towards the quiescent fluid, and.

Analysis of Nur 1 and Shr 1
causes this decrease as a result of the increase in the N t parameter.Impact of Le number on the distribution of Nur 1 has been shown in figure 8(b).It has been   found that the distribution of Nur 1 decreases with the increase of the values of Le number.Thermal boundary layer thickness rises with Le number, which in turn causes Nur 1 to decrease.The influence of Le number on Shr 1 has been depicted in Figure 9.It has been found that the profile of Shr 1 increases with the increase of Le number.When the surface velocity exceeds the free stream velocity, Shr 1 rises with an increase in Le number for N t < 0.1 but falls with an increase in Le number for N t greater than 0.1.

Conclusion
In the current article, the non-Darcian nanofluid flow and heat transfer across a stretched surface in a nanofluid-saturated porous media with the influence of a magnetic field has been studied.Using suitable similarity transformations, the basic equations of flow are transformed into ODE's.Then, the mathematical calculations are carried out using a finite difference algorithm with the aid of Newton's linearization method.The effects of different flow-controlled variables have been illustrated numerically and graphically for the distributions of velocity, temperature, concentration, Nur 1 , and Shr 1 .The important outcomes of this work are listed below: The obtained results of the current study will become very much helpful in various engineering applications such as glass blowing, spinning of fibers, and continuous casting of metals, particularly in various manufacturing processes like transpiration cooling, fabric cleaning, and laser pulse heating.
0 are known.Δh n+1 i at i = 0 and i = m is zero.The expression for a i , b i , c i and d i in equation 23 are given:

Figure 3 . 5 . 1 .
Figure 3. Influence of (a) N b , N t , and (b) Le on temperature profiles.5.1.Analysis of temperature profileAt various N b , N t and at Le = 5, Pr = 5, Bi = 0.1, M = 0.5, K = 0.3, Fr 1 = 0.5, Figure3(a) displays the temperature distributions in the TBL.The N t parameter has a significant effect on the nanoparticle volume fraction.N t is the process of moving particles in such a way that the particles become cool under the influence of temperature gradient.With a rise in N t values, a significant improvement in the temperature profile is seen due to the enhancement of molecular movement.HT in the boundary layer increases as N t increases, aggravating particle deposition away from the fluid region and increasing the volume percentage of nanoparticles.Due to spontaneous diffusivity, the nanoparticles reconfigure themselves and form a novel structure, which in turn enhances the NF's TC.The particle fury and increase in the fluid temperature in the boundary layer caused by the N b , increases the temperature distributions.With the increase of the values of N b and N t , it has been found that the surface temperature rises, the boundary layer thickens, and the curves become steep, as a result, the profile of Nur 1 decreases.Figure3(b) displays the profiles of temperature distributions for various values

Figure 4 .
Figure 4. Influence of (a) Bi and (b) Pr on temperature profiles.

Figure 5 .
Figure 5. Impacts of (a) M and (b) K on profiles of velocity.
) depicts the influence of M on the distribution of velocity.Velocity exhibits a decreasing tendency as M rises.With higher M values, a strong Lorentz effect significantly reduces the velocity.Applying a perpendicular magnetic field to an electrically conducting fluid, Lorentz force is produced.The velocity distribution diminishes due to the resistive nature of the Lorentz force, which reduces the motion of the particles.Higher values of the magnetic field indicate a larger body force, which can slow down the liquid flow and increases the thickness of the momentum boundary layer, as viscous force and body force are proportional to hydromagnetic force, as described by M. It shows unequivocally how the transverse magnetic field interferes with the transport phenomena.

Figure 5 (
b) shows how the distribution of velocity is affected by the parameter K.The velocity profiles start to decay as K numbers rise.Porosity physically prevents fluid from moving freely, reducing the fluid's movement.The presence of permeable space improves the safety of the fluid stream by reducing the rate of the fluid.Porousness generally raises fluid flow resistance, which leads to thicker thermal boundary layers and higher temperature profiles.For fixed value of N t , N b , Pr, Le, K, M, the impact of the Fr 1 number on the distribution of velocity has been shown in Figure6.Since the inertia coefficient is directly proportional to

Figure 6 .
Figure 6.Impacts of Fr 1 on profiles of velocity.

Figure 7 (
Figure 7(a) presented the distribution of concentration for various values of Le at N b = N t = 0.1, Pr = 5, M = 0.5, K = 0.3, Fr 1 = 0.5.A higher Le suggests a lower Brownian diffusion coefficient for a base fluid of particular kinematic viscosity.In the boundary layer region, Le number is a measure of how much the thermal diffusion rate contributes to the species diffusion rate.When Le is equal to 1, heat and species disseminate equally, and when Le is greater than 1.0, heat diffuses more quickly than species.As the Le number rises, the concentration profiles become steeper and the species

Figure 8 (
Figure 8(a) has been illustrates the influence of Pr number on the distribution of Nur 1 .Nur 1 measures the dimensionless HT at the sheet.The thickness of TBL's diminishes as the Pr number rises (figure 4b), while the steepness of the curves increases.It is observed that the profile of Nur 1 increases as Pr increases, and decreases as N t increases.This is because for greater values of Pr, the convection process is dominant compared to conduction, which increases the HT.The higher thermophoretic force that propels the high thermal conductivity nanoparticles from the hot sheet towards the quiescent fluid, and.causes this decrease as a result of the increase in the N t parameter.Impact of Le number on the distribution of Nur 1 has been shown in figure 8(b).It has been

Figure 7 .
Figure 7. Influence of (a) Le, and (b) N b , N t on profiles of concentration.

Figure 8 .
Figure 8. Influence of (a) Pr, and (b) Le on profiles of Nur 1 .

Figure 9 .
Figure 9. Influence of Le on profiles of Shr 1 .
i.It has been found that the distribution of fluid temperature and nanoparticle concentration enhances with the increase of the values of N b , and N t .With the increase of the values of N t and N b , the strength of the Brownian motion and thermophoresis has been found to have a minimal impact on the distribution of concentration.ii.It has been demonstrated that the increment of Pr number causes the reduction of temperature distribution.On the other hand, increasing Le number enhances the temperature distribution.iii.It has been found that as the Le number rises, the concentration profiles become steeper and the species border layer becomes thinner.In addition, it has been noticed that a rise in Bi number is accompanied by a significant increase in the temperature profiles.iv.It has been revealed that as the number of M, K, and Fr 1 increases, the distribution of fluid velocity decreases.Moreover, it has been demonstrated that the profile of Nur 1 increases as Pr increases, and decreases as N t increases.v.At the end, it has been encountered that the Shr 1 number raises as Le number increases for N t < 0.1, but decreases as Le number increases for N t > 0.1.

Table 1 .
Numerical outcomes for Nur 1 at N t = N b = 0..