Statistical Analysis of the Mathematical Model of Entropy Generation of Magnetized Nanofluid

This investigation introduces a mathematical model of entropy generation for Magnetohydrodynamic (MHD) peristaltic wave of nanofluid. The governing equations have been created by the supposition of low Reynolds number and long wavelength estimation. The scientific arrangement has been procured with the help of perturbation technique. The concentration profile, temperature profile, pressure distribution and friction forces are shown graphically for some important parameters. Further, the eventual outcomes of connection between the entropy generation and some various parameters have been plotted by means of correlation and regression. It is fundamental to find the affectability of each parameter on entropy generation.


Introduction
In fluid mechanics, we think about the conduct of particles at each point inside a space under different physical conditions. Numerical models have been utilized for various sorts of fluids, such as Newtonian fluid [1,2] and non-Newtonian fluid [3,4] to depict the physical marvels in fluid mechanics. There are numerous special instances of non-Newtonian fluid (e.g., nanofluids, micropolar fluids) and numerous studies [5][6][7][8][9][10][11][12] have utilized a certain logical technique to illuminate distinctive kinds of fluids by creating models. Nanotechnology is considered to be an ideal innovative answer for tackling the worldwide energy crisis. Indeed, nanofluid is the fluid suspension of nanostructures, which predominates and guarantees the fundamental amplification of heat transfer properties of the fluid. This helps us to encounter the prospective complexities of fluids in various fluid configuration.
The use of heat transfer fluids is one of the technological applications of nanoparticles, and it possesses an enormous capacity to suspend nanoparticles and confront cooling problems in thermal systems. Due to the great demands placed upon heat transfer fluids in terms of decreasing or increasing energy release to systems, a significant research work was undertaken by Choi and Eastman [13,14] in which a mixture of nanoparticles and base fluid were designated as "nanofluid". They defined a liquid of ultra-fine particles with sizes less than 100 nm. In the field of thermal engineering and heat transfer, nanofluid has always been an engrossing term. Peristalsis has various applications in connection with nanofluid, such as in engineering, bio-sciences and industries. Several theoretical and experimental attempts in this area have been conducted in the past. Specifically, the works of Latham [15] play a very important role in this connection. Similarly, because of its multiple advantages, research findings on peristaltic flow have received wider applications in industries, and numerous attempts have been made in literature to explore this direction of research. Abbas et al. [16] discussed the application of drug delivery systems under the influence of MHD with peristaltic motion. Bhatti et al. [17] examined the combined effects of MHD and partial slip on the peristaltic flow of nanofluid. Similarly, Salleh et al [18]

Mathematical Formulation
We present a demonstration of the peristaltic movement with thick, electrically leading and incompressible nanofluid properties through a two-dimensional, non-uniform channel with sinusoidal wave engendering towards down its wall. As shown in Figure 1, a Cartesian coordinate framework is used so that the x axis is considered alongside the middle line, toward the wave propagation, with y axis traversing to it. The B 0 , a uniform, outer attractive field, is forced along the y pivot and the initiated attractive field is thought to be irrelevant. The geometry of the divider surface is characterized as The governing equation of motion, continuity, thermal energy and nanoparticle fraction for peristaltic nanofluid can be written as [16]  The governing equation of motion, continuity, thermal energy and nanoparticle fraction for peristaltic nanofluid can be written as [16] Now let us consider the assumptions of a long wavelength number and low Reynolds approximations in the sense of creeping flow. By using the dimensionless quantities in Equations (2)-(6), we get the resulting equations in a simplified form as subject to the respective boundary conditions: In the presence of a magnetic field, the entropy generation can be derived from energy and entropy balance for the case of heat and mass transfer as [43]  Now let us consider the assumptions of a long wavelength number and low Reynolds approximations in the sense of creeping flow. By using the dimensionless quantities in Equations (2)-(6), we get the resulting equations in a simplified form as 1 + R n Pr subject to the respective boundary conditions: In the presence of a magnetic field, the entropy generation can be derived from energy and entropy balance for the case of heat and mass transfer as [43] The dimensionless form of the entropy generation number can be expressed as follows: where Ω, B r , Λ , Γ, ζ are the dimensionless temperature difference, Brinkman number, concentration difference, diffusive coefficient and constant parameter, represented as For nanofluid, the viscosity model and thermal conductivity can be defined as [50] where µ f κ f and κ p , are the viscosity f base fluid, thermal conductivities of the nanofluid and nanoparticle, respectively.

Solution of Problem
With the help of HPM [16], Equations (7)-(9) can be written as: and the initial guess and linear operators for Equations (16)-(18) are defined as which defines the following expansion: Using the expansion series defined in the terms mentioned in Equations (23)- (25) and incorporating them into the Equations (16)-(18) we get a system of linear differential equations and their relevant boundary conditions. Applying the scheme of HPM and comparing the powers ofq, we obtain the solution asq → 1 . We obtained the temperature distribution, velocity profile, and concentration profile.
Utilizing the expensing arrangement characterized in terms of (ϑ(x, y),(Θ(x, y) and (w(x, y)) as referenced in Equations (23)- (25) and using into the Equations (16)-(18), we get an arrangement of direct differential equations with their significant limit conditions. By contrasting the forces of q → 1 , we apply the scheme of HPM to determine the arrangement as q → 1, and obtain the required arrangements of temperature circulation, speed profile, and fixation profile.

Results and Discussion
In this section the obtained results are discussed. As shown in Figure 2, for higher values of N b and N t the temperature profile increases. This is because the Brownian motion creates micro-mixing, which increases thermal conductivity. Figure 3 shows that concentration force has the opposite behavior for various values of N b and N t . It is observed from Figure 4 that pressure distribution has opposite effect for the various values of thermal Grashof parameter Gr T . and basic density Grashof number Gr F . The thermal Grashof parameter proclaims the general impact of thick hydrodynamic power and thermal buoyancy forces. For Gr T < 1, the peristaltic routine is ruled by viscous powers, which is the other way around for Gr T > 1. The Gr F parameter is basically the proportion of species buoyancy forces to the thick hydrodynamic forces. For a situation where the two forces are equivalent (e.g., Gr F = 1), velocity will be minimized. We can conclude from Figure 5a that pressure rise reduces for the larger values of magnetic parameter M, which shows the fact that pressure can be controlled by using the suitable magnetic field. Also, it is concluded from this figure that flow can pass easily without imposing higher pressure inside the channel. From Figures 5b and 6, it is observed that friction force has completely the opposite behavior for the different values of the same physical parameters as compared with pressure rise distribution. Figure 7 shows a comparison of the graphical results of pressure rise and velocity profile. In Equation (7), by taking W e = 0, M = 0, κ → ∞, Gr F = 0, Gr T = 0, the current results can be reduced to the results obtained by Shapiro et al. [39] for the Newtonian case. Further results can also be reduced to the similar results obtained by Gupta and Seshadri [41] and Mekheimer [42] by taking W e = 0, κ → ∞, Gr F = 0, Gr T = 0 . 1, we apply the scheme of HPM to determine the arrangement as q ́ →1, and obtain the required arrangements of temperature circulation, speed profile, and fixation profile.

Results and Discussion
In this section the obtained results are discussed. As shown in Figure 2, for higher values of and the temperature profile increases. This is because the Brownian motion creates micromixing, which increases thermal conductivity. Figure 3 shows that concentration force has the opposite behavior for various values of and . It is observed from Figure 4 that pressure distribution has opposite effect for the various values of thermal Grashof parameter and basic density Grashof number . The thermal Grashof parameter proclaims the general impact of thick hydrodynamic power and thermal buoyancy forces. For < 1, the peristaltic routine is ruled by viscous powers, which is the other way around for > 1. The parameter is basically the proportion of species buoyancy forces to the thick hydrodynamic forces. For a situation where the two forces are equivalent (e.g., = 1), velocity will be minimized. We can conclude from Figure   5a that pressure rise reduces for the larger values of magnetic parameter , which shows the fact that pressure can be controlled by using the suitable magnetic field. Also, it is concluded from this figure that flow can pass easily without imposing higher pressure inside the channel. From Figure 5b and Figure 6, it is observed that friction force has completely the opposite behavior for the different values of the same physical parameters as compared with pressure rise distribution. Figure 7 shows a comparison of the graphical results of pressure rise and velocity profile. In Equation (7), by taking = 0, = 0, → ∞, = 0, = 0, the current results can be reduced to the results obtained by Shapiro et al. [39] for the Newtonian case. Further results can also be reduced to the similar results obtained by Gupta and Seshadri [41] and Mekheimer [42] by taking = 0, → ∞, = 0, = 0.   1, we apply the scheme of HPM to determine the arrangement as q ́ →1, and obtain the required arrangements of temperature circulation, speed profile, and fixation profile.

Results and Discussion
In this section the obtained results are discussed. As shown in Figure 2, for higher values of and the temperature profile increases. This is because the Brownian motion creates micromixing, which increases thermal conductivity. Figure 3 shows that concentration force has the opposite behavior for various values of and . It is observed from Figure 4 that pressure distribution has opposite effect for the various values of thermal Grashof parameter and basic density Grashof number . The thermal Grashof parameter proclaims the general impact of thick hydrodynamic power and thermal buoyancy forces. For < 1, the peristaltic routine is ruled by viscous powers, which is the other way around for > 1. The parameter is basically the proportion of species buoyancy forces to the thick hydrodynamic forces. For a situation where the two forces are equivalent (e.g., = 1), velocity will be minimized. We can conclude from Figure   5a that pressure rise reduces for the larger values of magnetic parameter , which shows the fact that pressure can be controlled by using the suitable magnetic field. Also, it is concluded from this figure that flow can pass easily without imposing higher pressure inside the channel. From Figure 5b and Figure 6, it is observed that friction force has completely the opposite behavior for the different values of the same physical parameters as compared with pressure rise distribution. Figure 7 shows a comparison of the graphical results of pressure rise and velocity profile. In Equation (7), by taking = 0, = 0, → ∞, = 0, = 0, the current results can be reduced to the results obtained by Shapiro et al. [39] for the Newtonian case. Further results can also be reduced to the similar results obtained by Gupta and Seshadri [41] and Mekheimer [42] by taking = 0, → ∞, = 0, = 0.                    From Table 1, the R-square entropy generation for various values of magnetic parameter M is determined to be 0.809, meaning that approximately 81% of the variability of entropy generation is explained by parameter M in the model. On the other hand, the adjusted R-square 0.799 indicates that about 80% of the variability of entropy generation is accounted for Magnetic parameter M by the model. The entropy generation value for Brownian motion parameter N b is 0.998, which indicates that approximately 99% of the variability of entropy is due to parameter N b in the model, while the adjusted R-square 0.999 indicates that about 99% of the variability of entropy is accounted for N b by the model. In the R-square, the value of entropy generation for thermophoresis parameter N t is 0.403, which reveals that approximately 40% of the variability of entropy is explained by parameter N t in the model, while the adjusted R-square 0.370 indicates that about 37% of the variability of entropy is accounted for N t by the model. The entropy value for different values of B r is 1.00 with 100% of the variability of entropy accounted for parameter B r in the model, while the adjusted R-square 1.00 indicates that about 100% of the variability of entropy is accounted for B r by the Model. From Table 2, a decrease of −2.562 in entropy for independent variable M and an increase of 2.029 in Entropy for N b can be concluded. Similarly, an increase of 6.307 in entropy for parameter N t and an increase of 68.492 in entropy for B r scores can be concluded for every single-unit increase in Iv, assuming all other variables in the model are constant. Table 3 is plotted to break down the relationship of entropy generation for some delicate parameters. It is determined from these outcomes that a huge impeccable significant positive connection exists between Brinkman number B r and its entropy generation at 0.01 level. A solid positive relationship has been seen from the connection results between entropy and the parameters N t and N b . An extremely negative relationship was seen between M and its entropy.

Conclusions
The following outcomes demonstrated through this study are: • Temperature profile increases with higher values of N b and N t . • Pressure distribution and friction have an inverse conduct for bigger estimations of the magnetic parameter, Brownian movement parameter and thermophoresis parameter.

•
The variability of entropy generation is 81% for the values of M, while 99% variability for the parameter N b .

•
The variability of entropy generation is 40% for the values of N t , while 100% variability for the parameter B r . Greek Symbols κ nanofluid thermal conductivity (W/m K) µ viscosity of the fluid N s/m 2 Φ nanoparticle volume fraction σ electrical conductivity (S/m) δ wave number m −1 c p effective heat capacity of nanoparticle (J/K) ν nanofluid kinematic viscosity m 2 /s (ρ) p nanoparticle mass density kg/m 3 ρ f fluid density kg/m 3 ρ f 0 fluid density at the reference temperature (T 0 ) kg/m 3 ζ volumetric expansion coefficient of the fluid (ρc) f heat capacity of fluid (J/K) λ wavelength (m) Φ amplitude ratio µ n f viscosity of nanofluid θ Temperature