Darcy–Forchheimer Magnetized Nanofluid flow along with Heating and Dissipation Effects over a Shrinking Exponential Sheet with Stability Analysis

Nanoparticles have presented various hurdles to the scientific community during the past decade. The nanoparticles dispersed in diverse base fluids can alter the properties of fluid flow and heat transmission. In the current examination, a mathematical model for the 2D magnetohydrodynamic (MHD) Darcy–Forchheimer nanofluid flow across an exponentially contracting sheet is presented. In this mathematical model, the effects of viscous dissipation, joule heating, first-order velocity, and thermal slip conditions are also examined. Using similarity transformations, a system of partial differential equations (PDEs) is converted into a set of ordinary differential equations (ODEs). The problem is quantitatively solved using the three-step Lobatto-three formula. This research studied the effects of the dimensionlessness, magnetic field, ratio of rates, porosity, Eckert number, Prandtl number, and coefficient of inertia characteristics on fluid flow. Multiple solutions were observed. In the first solution, the increased magnetic field, porosity parameter, slip effect, and volume percentage of the copper parameters reduce the velocity field along the η-direction. In the second solution, the magnetic field, porosity parameter, slip effect, and volume percentage of the copper parameters increase the η-direction velocity field. For engineering purposes, the graphs show the impacts of factors on the Nusselt number and skin friction. Finally, the stability analysis was performed to determine which solution was the more stable of the two.


Introduction
Flow across a Darcy-medium has extensive uses and considerable value in medicinal, chemical, and contemporary ecological frameworks. Numerous processes, including oil applications in different processes, thermal management in geothermal exchange formats, atomic head waste processes, water development, and water purification, are included. The classic model of Darcy is comprise the non-Darcian porous medium, which is a constrained form of this medium for the aforementioned applications and operations. The term Darcy's Law clarifies the fluid movement through a porous medium. This rule is applicable in scenarios with low porosity and low velocity. Forchheimer showed the velocity square factor in Darcy's velocity equation to analyze the boundary and inertia features. Muskat and Meres [1] referred to this as the Forchheimer word using Maxwell nanofluid flow through an isothermally heated stretching sheet. Rasool et al. [2] explained the effect of of The influence of slip conditions on fluid flow have not received a lot of attention, particularly when it comes to nanofluids, as shown by a comprehensive review of the research that has been published. A great number of important fluid applications exhibit slip boundary circumstances, such as the improvement of valves of the heart and cavity interiors, as well as the cleaning of prosthetic valves of the heart [35]. It is important to note that the condition of there being no slippage does not always hold in actual practice. It is possible to provide a straightforward explanation for the slip velocity condition by pointing out that flowing liquids do not have nil velocity with respect to their interaction with the barrier of a solid. Andersson [36] was likely the first person to present the idea of slip impact having an impact on the flow of a boundary layer. Wang et al. [37] looked at what happened to a Maxwell nanofluid when it slipped over a stretch that grew longer and longer at an exponential rate. The researchers discovered that as the slip parameter velocity component was increased, a diminishing behavior was seen in both directions. According to the findings of Saleem and Abd El-Aziz [38], the effect of fluid friction irreversibility decreases near to a surface as the slip parameter increases; nevertheless, farther away in the flow system, the irreversibility of heat transfer has been shown to be a more prominent factor. The unsteady flow on an exponential sheet was examined by Haider et al. [39], who discovered that the slip condition had a diminishing influence on the skin friction coefficient. Imran et al. [40] examined the behavior of an MHD generalized Maxwell fluid when it was applied to an exponentially accelerating infinite vertical surface along with a slip condition. According to Reddy et al. [41], a rise in the slip condition causes the velocity profiles to increase as well.
Joule heating has also been one of the most intriguing impacts to be implemented since joule heating has a significant impact on the MHD flow of fluids. Ohmic or joule heating is the process of converting electrical energy into thermal energy via material resistive losses. In addition, the heating of joule impact is used extensively and experimentally in the majority of electrical and electronic equipment. Reddy et al. [42] described the effect of ohmic heating on the flow of nanofluids along elastic barriers. They explained that the existence of ohmic heating might result in a temperature rise. Maskeen et al. [43] investigated the MHD flow in a Cu-Al 2 O 3 /water hybrid nanofluid via a vertically extending cylinder. Sajid et al. [44] reported a computational solution for the MHD flow of a ferrofluid in a curved channel with semi-porosity and ohmic heating. Patel and Singh [45] considered joule heating in their study of Walters-B liquid flow. Furthermore, Kamran et al. [46], Gholinia et al. [47], and Khan et al. [48] investigated the joule heating effect numerically in the analysis of micropolar fluid, Williamson fluid, and Casson fluid. Khan et al. [49] and Hussain et al. [50] independently examined the application of joule heating to a stretching cylinder for MHD Carreau and MHD Sisko nanofluid flows. Yan et al. [51] examined the MHD Cu-Al 2 O 3 /water flow of a hybrid nanofluid via an exponentially reducing permeable sheet and discovered binary solutions. Khan et al. [52] recently investigated the Eyring-Powell flow of fluids along with ohmic heating and the variable viscosity impacts for the material of the wire coating application.
The current work theoretically studies the MHD Darcy-Forchheimer flow of a waterbased nanofluid across an exponentially diminishing surface in the presence of ohmic heating, viscous dissipation, and first-order slip conditions, as inspired by the aforementioned literature. Nanofluids have been thoroughly discussed in the aforementioned literature; thus, a Cu-water-based nanofluid was selected as an example to explore the influences of the different physically applied factors. Four aspects of the heat transfer concert of nanofluids are called into question. (i) If no suction is applied, are double solutions of similarity conceivable for flow caused by an exponentially contracting sheet? (ii) Does an increase in the Forchheimer parameter and Eckert (ohmic heating) number cause an interruption in the separation of the boundary layer and a decrease in the rate of heat transfer? (iii) Does the use of a Cu-water nanofluid increase the rate of heat transfer compared to a conventional fluid? (iv) Does the increasing volume fraction of copper nanoparticles in a fluid based on water increase the rate of heat transfer? Consequently, this paper will provide answers to each of the preceding questions. Using exponential transformations, the governing model is transformed into ODEs. All computations are performed with a bvp4c solver, and, in a few instances, it is attempted to compare the numerical outcomes with the previously published literature. The authors are assured that the current piece of research is novel and will have a substantial impact on future researchers in the field of fluid dynamics.

Mathematical Description of the Problem
Let us take into account the effects of viscous dissipations, joule heating, and slip conditions in the context of the steady two-dimensional incompressible flow of an electrically conducting nanofluid across a diminishing surface in a porous medium along with Darcy-Forchheimer. The physical model of the considered problem can be seen in Figure 1.
interruption in the separation of the boundary layer and a decrease in the rate of hea transfer? (iii) Does the use of a -water nanofluid increase the rate of heat transfe compared to a conventional fluid? (iv) Does the increasing volume fraction of coppe nanoparticles in a fluid based on water increase the rate of heat transfer? Consequently this paper will provide answers to each of the preceding questions. Using exponentia transformations, the governing model is transformed into ODEs. All computations are performed with a bvp4c solver, and, in a few instances, it is attempted to compare the numerical outcomes with the previously published literature. The authors are assured that the current piece of research is novel and will have a substantial impact on future researchers in the field of fluid dynamics.

Mathematical Description of the Problem
Let us take into account the effects of viscous dissipations, joule heating, and slip conditions in the context of the steady two-dimensional incompressible flow of an elec trically conducting nanofluid across a diminishing surface in a porous medium along with Darcy-Forchheimer. The physical model of the considered problem can be seen in Figure 1. The assumption of a system of cartesian coordinates, in which the -axis is paralle to the contracting surface and the -axis is perpendicular to it, is also taken into account To the perpendicular of the shrinking sheet, a uniform magnetic field ( ) = i applied with strength, where is a constant magnetic resilience. The properties of the nanofluid and solid along with water are given in Tables 1-2. Table 1. Thermophysical properties of nanofluid [44].

Properties Nanofluid
Dynamic viscosity = (1 − ) . The assumption of a system of cartesian coordinates, in which the x-axis is parallel to the contracting surface and the y-axis is perpendicular to it, is also taken into account. To the perpendicular of the shrinking sheet, a uniform magnetic field B(x) = B 0 e x 2l is applied with strength, where B 0 is a constant magnetic resilience. The properties of the nanofluid and solid along with water are given in Tables 1 and 2. Table 1. Thermophysical properties of nanofluid [44].

Properties Nanofluid
Dynamic viscosity Density where subscript s indicates the solid properties of the copper. Thermal conductivity Table 2. Water and copper thermo-physical characteristics [11]. We disregard the induced magnetic field, and, subsequently, the Reynolds number is so very low. The governing equations for a steady nanofluid flow under these conditions are as follows [40,51]: ∂u ∂x The associated BCs related to (1)(2)(3) are In these equations v and u represent components of velocity in the yand x-directions, correspondingly, ρ denotes the density of the fluid, σ is the electrical conductivity of the fluid, b is the local inertia coefficient, K is the porous medium permeability, T is the fluid temperature, and subscript nf shows the properties of the nanofluid explained in Table 1.
is the wall temperature, and T ∞ is the ambient temperature. In addition, 2l is the velocity slip factor, and D = D 1 e − x /2l is the thermal slip factor.
These transformations of similarity are used to obtain the corresponding similarity solutions: The stream function ψ is expressed in components of velocity as K = 2K 0 e − x /l is used to represent the porous media permeability. By applying Equations (5)-(6) to Equations (1)-(3), the equation of continuity is fulfilled, and the equations of momentum and energy can be expressed as along with the BCs: is the velocity slip, and 2ϑ f l is the thermal slip parameter. The coefficient of skin friction and the local Nusselt number are all relevant physical quantities, and they can be written as By applying Equations (5)- (6) in Equation (10), we obtain

Linear Stability Analysis
Recently, numerous scholars have examined multiple fluid solutions for various fluid kinds and flow circumstances. From an experimental standpoint, it is worthwhile to determine which solution is physically dependable and applicable. Therefore, linear stability is necessary for validating the dependability of solutions. Merkin [53] and Weidman et al. [54] suggested reducing the controlling boundary layer of Equations (3)-(5) to the following unsteady forms for the stability analysis: where t represents time. Additionally, a novel similarity transformation is presented as Using Equation (14), Equations (12)-(13) can be written as follows: To test the stability of the steady flow solutions solving the boundary value problem (7-9), where θ(η) = θ 0 (η) and f (η) = f 0 (η), one can express where f 0 (η) and θ 0 (η) are minor relatives of F(η) and G(η), respectively. Moreover, ε represents the unknown eigenvalues. When the eigenvalue problem (16)(17)(18) is solved, an endless collection of eigenvalues is obtained. We must choose the least eigenvalue from this set. If the smallest eigenvalue (ε) is negative, this shows that the flow is unstable and demonstrates the growth of disruptions and physical impossibility. If the smallest eigenvalue is positive, on the other hand, this indicates that the solution is physically reliable and stable. The following equations result from applying Equation (18) in (15)(16): Subject to the BCs, Conferring with Haris et al. [55], in order to find the stability of (19-21), one boundary condition on F 0 (η) and G 0 (η) must be relaxed. In this case, it is important to note that we relaxed F 0 (η) → 0 as η → ∞ into F 0 (0) = 1. In addition, we fixed all parameters, including F S = 0.5, K 1 = 0.5, Pr = 6.2, Ec = 0.2, δ = δ T = 0.1, and varied values of S and ∅.

Result and Discussion
The system of PDEs (2-3) is transformed into ODEs (7-8) using exponential transformations (5). The converted ODEs (7)(8) are then resolved by employing the bvp4c function along with the BCs (9) in MATLAB. A water-based nanofluid along with Darcy-Forchheimer phenomena over an exponentially shrunk sheet under the influence of different physical effects, parameters, the Hartmann number, the Prandtl number, the Eckert number, thermal slip, suction, Forchheimer, velocity slip, and porosity were studied. Before beginning to analyze the outcomes of the problem at hand, it is necessary to validate that the numerical coding is functioning effectively; in this regard, Table 3 is built for comparison with the previously published data by Waini et al. [56] when Pr = 6.2, φ = 0 (regular fluid), K 1 = F S = M = Ec = δ = δ T = 0, for a stretching surface (i.e., f (0) = 1), and a great agreement is observed. Therefore, the existing numerical method and its MATLAB coding are used with full sureness to solve the problem at hand. In addition, we considered fluids at a temperature of 25 • C; therefore, the researchers recommended that the Prandtl number for water at 25 • C be equal to 6.2. Table 3. Magnitude of f (0) and −θ(0) under various estimates of S.

Waini et al. [56]
Current Results  Figure 2 depicts the variation of the velocity profile f (η) for changed values of the volume fraction parameter: φ = 0.01, 0.05, and 0.1 along different physical parameters. It is observed that the first solution is decreasing; this is because a nanofluid with a high-volume fraction has a lower velocity because a higher number of particles per unit volume means the fluid is more viscous or heavier, requiring more energy to move. On the other hand, the second solution is increasing near the surface for the high-volume fraction. Figure 3 illustrates the results of temperature variations θ(η) for different values of φ = 0.01, 0.05, and 0.1. It is observed from the figure that the temperature variation increases for both solutions. In practice, the nanoparticles disperse energy as heat. More nanoparticles produce more energy, raising the temperature and thickness of the thermal boundary layer. Figure 4 presents the variation of f (η) for changed estimates of the Hartmann number = 0, 0.3, and 0.5. It is observed that the first solution decreases; it is believed that the application of a transverse magnetic field results in a Lorentz force similar to a drag force, which tends to resist the fluid flow and thus reduces its velocity in the profiles of the first solution. Meanwhile, the second solution has an increasing nature with an increase in the values of M. In Figure 5, the temperature variation for different values of M = 0, 0.3, and 0.5 is shown. It is noted that the estimates of θ(η) increase for both solutions with enhancing estimates of M. Due to the fact that a Lorentz force can enhance mass transport, which easily induces convective motions and energy transport, the thermal boundary layer thickness increases in both solutions.
Micromachines 2022, 13, x FOR PEER REVIEW 9 of 21 tion increases for both solutions. In practice, the nanoparticles disperse energy as heat.
More nanoparticles produce more energy, raising the temperature and thickness of the thermal boundary layer. Figure 4 presents the variation of ( ) for changed estimates of the Hartmann number = 0, 0.3, and 0.5. It is observed that the first solution decreases; it is believed that the application of a transverse magnetic field results in a Lorentz force similar to a drag force, which tends to resist the fluid flow and thus reduces its velocity in the profiles of the first solution. Meanwhile, the second solution has an increasing nature with an increase in the values of . In Figure 5, the temperature variation for different values of = 0, 0.3, and 0.5 is shown. It is noted that the estimates of ( ) increase for both solutions with enhancing estimates of . Due to the fact that a Lorentz force can enhance mass transport, which easily induces convective motions and energy transport, the thermal boundary layer thickness increases in both solutions.         . In Figure 6, estimates of ( ) are reduced in the first solution; this is because effective density and permeability are directly proportional to each other, and, as shown in Equation (2), high permeability causes the velocity and momentum thickness of the boundary layer to decrease. On the other hand, profiles of velocity are increased in the second solution as the values of = 0.1, 0.5, and 0.7 rise. In Figure 7, the variation of θ(η) increases in the first solution and reduces in the second solution. The variation of ( ) is presented in Figure 8 for different estimates of the slip parameter: = 0, 0.2, and 0.5. Here, it shows that the behavior of the first solution is decreasing, but the second solution is increasing initially and then decreasing as the value of rises. This decreasing behavior of the velocity profiles is due to the velocity slip, which is the difference between the particle velocity and the undisturbed nanofluid velocity at the particle location. The variation of ( ) is illustrated in Figure 9 and shows an increase in temperature for both the first and second solutions as = 0.1, 0.2, and 0.3.  Figures 6 and 7 depict the variation of f (η) and θ(η) for various values of permeability K 1 . In Figure 6, estimates of f (η) are reduced in the first solution; this is because effective density and permeability are directly proportional to each other, and, as shown in Equation (2), high permeability causes the velocity and momentum thickness of the boundary layer to decrease. On the other hand, profiles of velocity are increased in the second solution as the values of K 1 = 0.1, 0.5, and 0.7 rise. In Figure 7, the variation of θ(η) increases in the first solution and reduces in the second solution. The variation of f (η) is presented in Figure 8 for different estimates of the slip parameter: δ = 0, 0.2, and 0.5. Here, it shows that the behavior of the first solution is decreasing, but the second solution is increasing initially and then decreasing as the value of δ rises. This decreasing behavior of the velocity profiles is due to the velocity slip, which is the difference between the particle velocity and the undisturbed nanofluid velocity at the particle location. The variation of θ(η) is illustrated in Figure 9 and shows an increase in temperature for both the first and second solutions as δ T = 0.1, 0.2, and 0.3.
The variation of θ(η) is shown in Figure 10 for various estimates of the Eckert number: Ec = 0, 0.3, and 0.7. An increase in temperature is noticeable for both the second and first solutions, as the Ec values increase. Physically, Ec describes the enthalpy and kinetic energy connection. Consequently, it minimizes viscous fluid stress by converting kinetic energy into internal energy. A rise in internal energy causes fluid enhancement. The effect of the Prandtl number Pr is discussed in Figure 11. It is noticed that as the values of Pr = 1, 3, and 6.2 increase, the variations of θ(η) decrease for both the first and second solutions. The Prandtl parameter has an opposite relationship with the diffusivity of thermal energy. Greater values of Pr cause a lower diffusivity of thermal energy. The lower diffusivity of thermal energy brings a decrement in the temperature profiles and is connected with the thermal layer thickness.        . Changes of ( ) for . Figure 9. Changes of θ(η) for δ T .
hancement. The effect of the Prandtl number is discussed in Figure 11. It is noticed that as the values of = 1, 3, and 6.2 increase, the variations of ( ) decrease for both the first and second solutions. The Prandtl parameter has an opposite relationship with the diffusivity of thermal energy. Greater values of cause a lower diffusivity of thermal energy. The lower diffusivity of thermal energy brings a decrement in the temperature profiles and is connected with the thermal layer thickness. Figure 10. Changes of ( ) for . Figure 10. Changes of θ(η) for Ec.
Micromachines 2022, 13, x FOR PEER REVIEW 15 of 21 Figure 11. Changes of ( ) for . Figure 12 depicts the variation of ( ) for various values of Forchheimer . It is noted that the increase in = 0, 0.5, and 1 causes a decrease/increase in the variation of ( ) for both the first and second solutions. Physically, an initial drag force is produced due to a rise in , which acts as a barrier for ( ), which causes the velocity profile to dissipate. Figure 13 depicts the effect of the skin friction (0) coefficient Figure 11. Changes of θ(η) for Pr. Figure 12 depicts the variation of f (η) for various values of Forchheimer F s . It is noted that the increase in F s = 0, 0.5, and 1 causes a decrease/increase in the variation of f (η) for both the first and second solutions. Physically, an initial drag force is produced due to a rise in F s , which acts as a barrier for f (η) which causes the velocity profile to dissipate. Figure 13   The fluctuations of the smallest eigenvalues ε 1 against S and ∅ are displayed in Table 4 for F S = 0.5, K 1 = 0.5, Pr = 6.2, δ = δ T = 0.1, and Ec = 0.2. As explained by Equation (18), the flow is steady when the initial disturbance declines with time. This will be the case for ε 1 > 0. In the meantime, the flow for ε 1 < 0 is unstable due to the beginning rise of disturbance as time progresses, e −ετ → ∞ as ε 1 < 0 and τ → ∞ . According to Table 4, the values of ε 1 for the first solution are positive, whereas they are negative for the second solution. Consequently, this result demonstrates that the first solutions are physically stable and dependable while the second solutions are not.   The fluctuations of the smallest eigenvalues against and ∅ are displayed in Table 4 for = 0.5, = 0.5, = 6.2, = = 0.1, and = 0.2 . As explained by Equation (18), the flow is steady when the initial disturbance declines with time. This will be the case for > 0. In the meantime, the flow for < 0 is unstable due to the beginning rise of disturbance as time progresses, → ∞ as < 0 and → ∞. According to Table 4, the values of for the first solution are positive, whereas they are negative for the second solution. Consequently, this result demonstrates that the first solutions are physically stable and dependable while the second solutions are not.

Conclusions
The present study analyzed the effect of Darcy-Forchheimer on nanofluid flow through an exponentially shrinking surface. The governing boundary layer system of the partial equations was reduced into ODEs before solving them by the bvp4c technique in MATLAB software. The presence of duality is shown by the numerical results. The sta-

Conclusions
The present study analyzed the effect of Darcy-Forchheimer on nanofluid flow through an exponentially shrinking surface. The governing boundary layer system of the partial equations was reduced into ODEs before solving them by the bvp4c technique in MATLAB software. The presence of duality is shown by the numerical results. The stability analysis is conducted to check the stability of the first and second solutions. The numerical value of the smallest eigenvalue in the stability table showed that only the first solution is stable.
The coefficient of skin friction increases for the first solution while reducing for the second solution as the values of volume fraction enhance. The rate of heat transfer decreases for the first and second solutions with the increase in the volume fraction values. The variation of the velocity profile decreases/increases in both solutions as the magnitude of the Forchheimer parameter values increases.

Conflicts of Interest:
The authors declare no conflict of interest.