Electromagnetohydrodynamic Electroosmotic Flow and Entropy Generation of Third-Grade Fluids in a Parallel Microchannel.

The present paper discusses the electromagnetohydrodynamic (EMHD) electroosmotic flow (EOF) and entropy generation of incompressible third-grade fluids in a parallel microchannel. Numerical solutions of the non-homogeneous partial differential equations of velocity and temperature are obtained by the Chebyshev spectral collocation method. The effects of non-Newtonian parameter Λ, Hartman number Ha and Brinkman number Br on the velocity, temperature, Nusselt number and entropy generation are analyzed in detail and shown graphically. The main results show that both temperature and Nusselt number decrease with the non-Newtonian physical parameter, while the local and total entropy generation rates exhibit an adverse trend, which means that non-Newtonian parameter can provoke the local entropy generation rate. In addition, we also find that the increase of non-Newtonian parameter can lead to the increase of the critical Hartman number Hac.


Introduction
Microfluidic devices are widely demonstrated in areas of biomedical and biochemical analysis, and have been one of the powerful tools for studying basic physical processes [1][2][3]. In these processes, pressure gradients, electrical fields, magnetic fields or their suitable combinations are the popular actuation mechanisms. Compared with the previous single pattern of pressure-driven flow, increasing attention has been attached to electroosmotic and electromagnetic actuation mechanisms in recent years. With the rapid development of lab-on-a-chip technologies, electroosmosis has been widely utilized with advantages of high reliability and simple operation [4]. A variety of theoretical and experimental studies of electroosmotic flow (EOF) in a microchannel have been extensively performed for both Newtonian fluid [5,6] and non-Newtonian fluid [7][8][9][10][11]. The hydrodynamic dispersion-combined electroosmosis and magnetohydrodynamic effect has been analyzed in microchannels with slowly varying wall zeta potentials [12]. Heat-transfer phenomena that are associated with electroosmotic and pressure-driven flows in microchannels have also been studied for thermally fully-developed flows [13][14][15] and thermally developing flows [16,17].
Aside from electroosmosis mentioned above, the magnetohydrodynamic (MHD) flow has also attracted much attention due to its prospective applications in engineering and medical fields [18][19][20]. Meanwhile, in order to achieve more efficient flow control, electromagnetohydrodynamic (EMHD) flow has also received widespread attention, i.e., both the external electrical and magnetic field are applied to the conducting fluid. The interaction of electric field and transverse magnetic field can produce the Lorentz force, which is a non-intrusive way to influence the motion of EMHD flow. 3 of 18 are numerically solved by the Chebyshev spectral collocation method. Heat transfer characteristics, represented by the temperature distribution and Nusselt number, have been sketched in this study. Moreover, considering heat diffusion and heat convection, Joule heating effect, coupling electromagnetic diffusion, magnetic field and viscous friction, the entropy generation is discussed for different values of several pertinent non-dimensional parameters.

Physical Model and Explanation of the Problem
Here we consider thermally fully developed flow of third grade fluids in a microchannel under the magnetohydrodynamic electroosmotic influence. The flow is assumed to be steady, incompressible, viscous and electrically conducting. The physical model and the coordinate system are shown in Figure 1, Two external electrical field E x (in x * -axis direction) and E z (in z * -axis direction) are tangential to the charged surface, and a uniform magnetic field of strength B (in y * -axis direction) is perpendicular to the charged plates. A constant pressure gradient is imposed along the direction of flow, i.e., x * -axis direction. We assume the channel length L in x * -direction is much larger than both the channel height 2h in the y * -direction and the channel width W in the z * -direction, i.e., 2h, W << L; and the ratio of height 2h to width W of the channel is small enough, i.e., δ = 2h/W << 1. Therefore, the rectangular pipe flow is transformed into parallel plates flow.
Micromachines 2020, 11, x 3 of 18 energy equations are numerically solved by the Chebyshev spectral collocation method. Heat transfer characteristics, represented by the temperature distribution and Nusselt number, have been sketched in this study. Moreover, considering heat diffusion and heat convection, Joule heating effect, coupling electromagnetic diffusion, magnetic field and viscous friction, the entropy generation is discussed for different values of several pertinent non-dimensional parameters.

Physical Model and Explanation of the Problem
Here we consider thermally fully developed flow of third grade fluids in a microchannel under the magnetohydrodynamic electroosmotic influence. The flow is assumed to be steady, incompressible, viscous and electrically conducting. The physical model and the coordinate system are shown in Figure 1, Two external electrical field Ex (in x * -axis direction) and Ez (in z * -axis direction) are tangential to the charged surface, and a uniform magnetic field of strength B (in y * -axis direction) is perpendicular to the charged plates. A constant pressure gradient is imposed along the direction of flow, i.e., x * -axis direction. We assume the channel length L in x * -direction is much larger than both the channel height 2h in the y * -direction and the channel width W in the z * -direction, i.e., 2h, W << L; and the ratio of height 2h to width W of the channel is small enough, i.e., δ = 2h/W<<1. Therefore, the rectangular pipe flow is transformed into parallel plates flow.

Electrical Potential Distribution
The electric double layer (EDL) originates from the chemical contact interaction between the electrolyte solution and the channel wall. For a symmetric electrolyte solution, we assume that the ionic species concentration obeys the Boltzmann distribution. Therefore, the Poisson-Boltzmann equations are applied to describe the electric potential ψ* of the EDL , * 2 ε ρ ψ e − = ∇ (1) where ρe = −2n0zvesinh(zve0ψ/kbTa) is the local volumetric net charge density, ε is the dielectric constant of the medium, n0 is the bulk ionic concentration, zv is the ion valence, e is the electron charge, kb is the Boltzmann constant, and Ta is the absolute temperature. The parallel plates are assumed to be charged and bear a uniform wall zeta potential of ψ0, which is small compared to the thermal potential, i.e., ( ). Due to Debye-Hückel linearization approximation (sinh(zve0ψ/kbTa) ≈ zve0ψ/kbTa), the Poisson-Boltzmann equation and corresponding boundary conditions are simplified as:

Electrical Potential Distribution
The electric double layer (EDL) originates from the chemical contact interaction between the electrolyte solution and the channel wall. For a symmetric electrolyte solution, we assume that the ionic species concentration obeys the Boltzmann distribution. Therefore, the Poisson-Boltzmann equations are applied to describe the electric potential ψ* of the EDL where ρ e = −2n 0 z v esinh(z v e 0 ψ/k b T a ) is the local volumetric net charge density, ε is the dielectric constant of the medium, n 0 is the bulk ionic concentration, z v is the ion valence, e is the electron charge, k b is the Boltzmann constant, and T a is the absolute temperature. The parallel plates are assumed to be charged and bear a uniform wall zeta potential of ψ 0 , which is small compared to the thermal potential, i.e., ( z v e 0 ψ < k b T a ). Due to Debye-Hückel linearization approximation (sinh(z v e 0 ψ/k b T a ) ≈ z v e 0 ψ/k b T a ), the Poisson-Boltzmann equation and corresponding boundary conditions are simplified as: where κ is the Debye-Hückel parameter and 1/κ denotes the characteristic thickness of the EDL. The local volumetric net charge density ρ e can be easily evaluated

Flow Analysis and Mathematical Formulation
The velocity of the incompressible flow is governed by the continuity equation and the momentum governing equation: where ρ is the density of the fluid, V* = (u*, v*, w*) is the velocity vector, τ is stress tensor, and f is the body force vector acting on the flow, which is mainly composed of Lorentz force and the electrical force resulting from the electrokinetic effects. It can be written as where E = E x e x − E z e z is the applied electrical field vector along x * -axis and z * -axis directions, and both components in two directions are supposed to be invariant. B = Be y is the applied constant magnetic field along y * -direction. Here, we mention in particular that the induced magnetic field is neglected due to the small magnetic Reynolds number. J is the local ion current density vector and obeys the Ohm's law: where σ is the electrical conductivity. In general, the Cauchy stress tensor τ for incompressible third grade fluids is given by [37]: where p is the pressure, I is the identity tensor. µ denotes the dynamic viscosity and α i (i = 1, 2) and β i (i = 1, 2, 3) are the material constants. Ai (i = 1, 2, 3) are kinematic tensors with the following expressions: Due to continuity Equation (5) and the hypotheses 2h, W << L, only the axial velocity component u * (y * ) along x * -axis is taken into account. This assumption has been proven to correct by making a comparison with the analytical solution for velocity in three directions in reference [25].
Substituting the velocity field Equation (12) and Equation (7) in Equation (6), and the pressure gradient is assumed to be a constant C 0 , the momentum governing equation along x * -axis direction can be converted into the following form: Equation (13) (See Brief Steps in Appendix A) is subjected to the following no-slip boundary conditions: To obtain the dimensionless form of Equation (13), non-dimensional parameters are defined as follows: where U eo is the Helmholtz-Smoluchowski electroosmotic velocity, Λ is the dimensionless parameter related to the non-Newtonian behavior, Ha is the Hartman number, which represents the relative strength of the electromagnetic forces and the viscous forces, S is a non-dimensional quantity denoting the strength of the lateral electric field. Ω represents the estimate of applied pressure compared to electroosmotic force. K is so-called electrokinetic width. So the dimensionless form of velocity in Equations (13) and (14) are: u(y) = 0 ,at y = ±1.

Thermal Transport for Thermally Fully Developed Flow
It is well known that viscous dissipation becomes significant in the microscale conduits. Sadeghi and Saidi concluded that viscous dissipation effects in combined pressure and electroosmotically driven flow had an important infection to the microscale thermal transport characteristics [57]. Considering the viscous dissipation, electromagnetic coupling heat and volumetric joule heating, the energy governing equation can be given as [15]: where T * is the local temperature of the liquid, c p is the specific heat at constant pressure, k th is the thermal conductivity of the liquid. Considering a steady state and thermally fully developed flow, Equation (18) can be written as: The second and forth terms on the right side of Equation (19) represent the volumetric energy generation caused by the viscous dissipation and Joule heat, which induced by Joule heating effect together with the contribution from electromagnetic effect, especially for the case of a large strength of magnetic field [58].
Furthermore, in a thermally fully developed case, we have: where T s * and T m * are the local wall and mainstream temperature, respectively. Under the imposed constant heat flux (q s ) boundary condition, i.e., q s =h(T s * − T m * ) = const., whereh is the convective heat transfer coefficient, we obtain: Based on the above assumptions, the energy Equation (19) and relevant boundary conditions are: An overall energy balance for an elemental control volume on a length of duct dx * was considered: where u m * denotes axial mean velocity, and can be written as: Then, the constant mean temperature gradient dT * m/ dx * can be expressed from Equation (24) as: where, The coefficients A and D of Equation (27) are given by: It is convenient to express Equation (22) in the non-dimensional form by introducing the nondimensional parameters and variables as follows: Physically, the parameters S x and S z stand for the relative strength of Joule heating to wall heat flux, which can be viewed as the dimensionless Joule heat parameters, and Br is Brinkman number, which represents the ratio of heat generated by viscous dissipation to the applied wall heat flux. Then the non-dimensional Equation (22) is expressed by: The corresponding boundary conditions of the dimensionless energy equation are: Based on the dimensionless variables defined earlier, the bulk mean temperature θ m can be defined as: An important heat transfer parameter expressed as Nusselt number Nu can be written as: where D h denotes hydrodynamic diameter and D h = h for a half of microchannel height. From Equations (32) and (33), the finial local Nusselt number (at the upper wall) can be expressed as:

Entropy Generation Rate
According to the entropy generation minimization concept [46,59], the local volumetric rate of entropy generation based on the above obtained velocity and temperature field can be expressed as: where S G * is the volumetric entropy generation rate, which consists of five parts: heat diffusion irreversibility, Joule heating effect, coupling electromagnetic diffusion, magnetic field and viscous friction of the fluids. In this work, they are written respectively as: By use of the characteristic entropy transfer rate (k th /h 2 ), the dimensionless form of entropy generation rate can be given as: where, where Θ = (k th T s * )/(q s h) is a constant determined by the unaltered wall temperature and the heat flux, Pe = (ρc p U eo h)/k th is Peclet number, and the variable F in Equation (38) is: In addition, the total non-dimensional entropy generation can also be obtained:

Numerical Solution
Among numerous numerical methods, the Chebyshev spectral method has higher accuracy and wider application. In this paper, by utilizing the Chebyshev spectral collocation method, we study the non-dimensional EMHD velocity, temperature, Nusselt number and entropy generation of third-grade fluids between two parallel micro-plates, owing to the fact that the analytical solutions of these physical quantities are difficult to obtain for third-grade fluids when EMHD electroosmotic effects are all taken into condition.
The physical domain is [−1, 1] in the present analysis, and the Chebyshev points y j = cos(jπ/N), j = 0, 1, . . . , N are chosen to discretize the interval. Let u = [u(y 0 ), u(y 1 ), . . . , u(y N )] be the undetermined vector at the Chebyshev points, then we obtain a Chebyshev polynomial P of degree at most equal to N, i.e., P(y i ) = u(y i ), i = 0, 1, . . . , N. By differentiating P and evaluating at the grid points, we can transform the differential equation into linear algebraic equations, and the numerical solution of Equation (16) under Equation (17) can be obtained. By use of the velocity values obtained, the value of dimensionless temperature θ can be easily computed at each Chebyshev point by Equation (30) which is a second order differential equation, and then the value of the Nusselt number and entropy generation can be easily calculated by Equations (34) and (37).

Results and Discussion
In the following sections, the distributions for dimensionless velocity, temperature, Nusselt number and entropy generation rate will be discussed. The effects of non-Newtonian parameter, Joule heating, magnetic field intensity and electro-kinetic parameters on the above physical quantities will be shown graphically. Before proceeding, the permissible ranges of relevant physical parameters should be given firstly.
For typical microscale fluid flow, half-height of the channel h is about 100 µm, viscosity µ~10 −3 kg/(ms), electrical conductivity σ~2.2 × 10 −4 -10 6 S/m. The applied magnetic field B is 0.018-0.44T [21], so the range of Hartmann number Ha varies from 0 to 3 calculated from Equation (15) [32,60]. The strength of the applied electric field varies from 0 to 20 V/m and the electroosmotic velocity is U eõ 100 µm/s. Generally, S x and S z are both positive, and for simplicity we assume S x + S z = 1 [51]. We suppose the electrokinetic width K = 10 unless there is a special announcement [43], which is a typical thickness value for non-overlapping EDLs. Brinkman number (Br) is 0-0.04. What is worth noting, according to the theoretical analysis given by Sarkar et al. [25], is that the value of lateral electric field should not be too large, otherwise the induced transverse flow will not be neglected, which will contradict the assumption of unidirectional flow. Thus, if there is no special statement, the value of S is set to 1. In addition, Péclet number Pe is restricted to 0.5 [61].
Firstly, we have conducted comparisons including two special cases. In Figure 2, the result of the present numerical velocity (for Λ = 0) is compared to the result obtained by Chakraborty et al. [24], who have discussed the Newtonian fluid which is subjected to the combined action of electroosmotic and electromagnetic forces. Secondly, in Figure 3 the present results of velocity distributions and temperature distributions (for S = 0 and Ha ≈ 0) for different non-Newtonian parameters Λ are compared to the results of Akgül et al. [33], who have obtained approximate analytical solutions of the electroosmotic flow of the third grade fluids by the perturbation techniques. We can see that the present results are in good agreement with the earlier conclusion when the non-Newtonian parameter Λ is small enough, however, for larger Λ (0.025), perturbation method is no longer appropriate to solve the present problem. Therefore, the spectral collocation method in our analysis has a broader application.
which will contradict the assumption of unidirectional flow. Thus, if there is no special statement, the value of S is set to 1. In addition, Péclet number Pe is restricted to 0.5 [61].
Firstly, we have conducted comparisons including two special cases. In Figure 2, the result of the present numerical velocity (for Λ = 0) is compared to the result obtained by Chakraborty et al. [24], who have discussed the Newtonian fluid which is subjected to the combined action of electroosmotic and electromagnetic forces. Secondly, in Figure 3 the present results of velocity distributions and temperature distributions (for S = 0 and Ha ≈ 0) for different non-Newtonian parameters Λ are compared to the results of Akgül et al. [33], who have obtained approximate analytical solutions of the electroosmotic flow of the third grade fluids by the perturbation techniques. We can see that the present results are in good agreement with the earlier conclusion when the non-Newtonian parameter Λ is small enough, however, for larger Λ (0.025), perturbation method is no longer appropriate to solve the present problem. Therefore, the spectral collocation method in our analysis has a broader application.

Velocity Analysis
In order to highlight the electromagnetohydrodynamic electroosmotic effects, in the following sections we will discuss the flow and the heat-transfer characteristics without regard for the pressure which will contradict the assumption of unidirectional flow. Thus, if there is no special statement, the value of S is set to 1. In addition, Péclet number Pe is restricted to 0.5 [61].
Firstly, we have conducted comparisons including two special cases. In Figure 2, the result of the present numerical velocity (for Λ = 0) is compared to the result obtained by Chakraborty et al. [24], who have discussed the Newtonian fluid which is subjected to the combined action of electroosmotic and electromagnetic forces. Secondly, in Figure 3 the present results of velocity distributions and temperature distributions (for S = 0 and Ha ≈ 0) for different non-Newtonian parameters Λ are compared to the results of Akgül et al. [33], who have obtained approximate analytical solutions of the electroosmotic flow of the third grade fluids by the perturbation techniques. We can see that the present results are in good agreement with the earlier conclusion when the non-Newtonian parameter Λ is small enough, however, for larger Λ (0.025), perturbation method is no longer appropriate to solve the present problem. Therefore, the spectral collocation method in our analysis has a broader application.

Velocity Analysis
In order to highlight the electromagnetohydrodynamic electroosmotic effects, in the following sections we will discuss the flow and the heat-transfer characteristics without regard for the pressure

Velocity Analysis
In order to highlight the electromagnetohydrodynamic electroosmotic effects, in the following sections we will discuss the flow and the heat-transfer characteristics without regard for the pressure gradient effect. Figure 4 delineates the profiles of flow velocity with different fluid physical parameters. It can be seen from Figure 4a that the velocity decreases with the increase of non-Newtonian parameter Λ. The reason is that the increase of Λ is actually related to the increase of the viscosity of third-grade fluids. The effects of Hartman number Ha on the dimensionless velocity are shown in Figure 4b,c. The results are the same as those discussed in reference 24, for small Ha (Ha ≤ 1), the aiding force (σBE z ) is greater than the retarding one (−σB 2 u) which can be seen in Equation (16), therefore, the velocity increases with Ha (as depicted in Figure 4b). With the further increase of Ha, the equilibrium between aiding force and retarding force arrives, the corresponding value of Ha is so-called critical Hac. For Ha beyond the critical Hac (as depicted in Figure 4c), the retarding magnetic force becomes the leading factor and triggers a progressive reduction in the flow velocity. Finally, we can conclude from Figure 4d that the velocity increases with the augment of the strength of the lateral electric field S. Clearly by increasing the magnitude of the lateral electric field S, the aiding force (HaS) is dominant which promotes the increase of velocity. In order to gain a better understanding of the critical Hac, the profiles of average velocity have been pictorially depicted for S = 1 in Figure 5. It is easy to observe that the average velocity profile has been divided into two regions by the critical Hac, and shows an increasing-decreasing trend with the increase of Ha. It is worth mentioning that the value of critical Hac shows an increasing trend with the increase of non-Newtonian parameter Λ. the viscosity of third-grade fluids. The effects of Hartman number Ha on the dimensionless velocity are shown in Figure 4b,c. The results are the same as those discussed in reference 24, for small Ha (Ha ≤ 1), the aiding force (σBEz) is greater than the retarding one (−σB 2 u) which can be seen in Equation (16), therefore, the velocity increases with Ha (as depicted in Figure 4b). With the further increase of Ha, the equilibrium between aiding force and retarding force arrives, the corresponding value of Ha is so-called critical Hac. For Ha beyond the critical Hac (as depicted in Figure 4c), the retarding magnetic force becomes the leading factor and triggers a progressive reduction in the flow velocity. Finally, we can conclude from Figure 4d that the velocity increases with the augment of the strength of the lateral electric field S. Clearly by increasing the magnitude of the lateral electric field S, the aiding force (HaS) is dominant which promotes the increase of velocity. In order to gain a better understanding of the critical Hac, the profiles of average velocity have been pictorially depicted for S = 1 in Figure 5. It is easy to observe that the average velocity profile has been divided into two regions by the critical Hac, and shows an increasing-decreasing trend with the increase of Ha. It is worth mentioning that the value of critical Hac shows an increasing trend with the increase of non-Newtonian parameter Λ.

Temperature Analysis
The effects of various dimensionless parameters, including non-Newtonian parameter Λ, the magnetic field Ha and viscous dissipation Br, on the distribution of fluid temperature are discussed in Figure 6. Firstly, a decreasing trend in temperature with non-Newtonian parameter Λ is observed in Figure 6a. This variation tendency is consistent with the previous result in velocity because the decrease of velocity leads to the decrease of heat exchange in non-Newtonian fluids. From Figure 6b

Temperature Analysis
The effects of various dimensionless parameters, including non-Newtonian parameter Λ, the magnetic field Ha and viscous dissipation Br, on the distribution of fluid temperature are discussed in Figure 6. Firstly, a decreasing trend in temperature with non-Newtonian parameter Λ is observed in Figure 6a. This variation tendency is consistent with the previous result in velocity because the decrease of velocity leads to the decrease of heat exchange in non-Newtonian fluids. From Figure 6b,c, it can be found that the magnitude of dimensionless temperature increases with small Ha (Ha < Hac), the variation tendency is opposite for big Ha (Ha > Hac), and the maximum temperature has been obtained in the center of channel. Finally, the effects of viscous dissipation Br are discussed in Figure 6d and we can see that the increase of Br can trigger a tiny increase in temperature. The reason is that the viscous dissipation can be viewed as an energy source to increase the temperature of the fluid.

Temperature Analysis
The effects of various dimensionless parameters, including non-Newtonian parameter Λ, the magnetic field Ha and viscous dissipation Br, on the distribution of fluid temperature are discussed in Figure 6. Firstly, a decreasing trend in temperature with non-Newtonian parameter Λ is observed in Figure 6a. This variation tendency is consistent with the previous result in velocity because the decrease of velocity leads to the decrease of heat exchange in non-Newtonian fluids. From Figure 6b and c, it can be found that the magnitude of dimensionless temperature increases with small Ha (Ha < Hac), the variation tendency is opposite for big Ha (Ha > Hac), and the maximum temperature has been obtained in the center of channel. Finally, the effects of viscous dissipation Br are discussed in Figure 6d and we can see that the increase of Br can trigger a tiny increase in temperature. The reason is that the viscous dissipation can be viewed as an energy source to increase the temperature of the fluid.  Ha leads to a decreasing-increasing variation trend of Nu, regardless of the magnitude of Brinkman number. In Figure 7a, in the case of small Ha (i.e., the aiding effect of lateral electric field is stronger than that of magnetic field), the increase of flow velocity with Ha results in the decrease of the bulk mean temperature and the convective heat-transfer coefficient, and then the Nusselt number presents a downward trend. For different non-Newtonian parameter Λ, the Nusselt number Nu profiles have been graphically depicted in Figure 7c. We can observe that the Nusselt number decreases with increasing non-Newtonian parameter Λ for the case of Ha = 3, which is a natural result corresponding to the previous variations in velocity and temperature.
magnitude of the quantity Ts * −Tm * gradually increasing with the Brinkman number. The augment of values of Ha leads to a decreasing-increasing variation trend of Nu, regardless of the magnitude of Brinkman number. In Figure 7a, in the case of small Ha (i.e., the aiding effect of lateral electric field is stronger than that of magnetic field), the increase of flow velocity with Ha results in the decrease of the bulk mean temperature and the convective heat-transfer coefficient, and then the Nusselt number presents a downward trend. For different non-Newtonian parameter Λ, the Nusselt number Nu profiles have been graphically depicted in Figure 7c. We can observe that the Nusselt number decreases with increasing non-Newtonian parameter Λ for the case of Ha = 3, which is a natural result corresponding to the previous variations in velocity and temperature.

Entropy Generation Analysis
The profiles of local entropy generation of third-grade fluids at different values of Ha have been delineated in Figures 8a and b. The variation trend of local entropy generation is opposite for low values of Ha and the high ones. For low values of Hartman number, we can reduce entropy by increasing the intensity of the magnetic field, and we can see the maximum values of entropy have been reached at the center of parallel plates. This is a reasonable result by noticing the fact that most changes of velocity and temperature distributions in the previous figures occur at the center.

Entropy Generation Analysis
The profiles of local entropy generation of third-grade fluids at different values of Ha have been delineated in Figure 8a,b. The variation trend of local entropy generation is opposite for low values of Ha and the high ones. For low values of Hartman number, we can reduce entropy by increasing the intensity of the magnetic field, and we can see the maximum values of entropy have been reached at the center of parallel plates. This is a reasonable result by noticing the fact that most changes of velocity and temperature distributions in the previous figures occur at the center.
Moreover, the influence of Br on the local entropy generation is illustrated in Figure 8c. It can be observed that the local entropy generation shows an increasing trend with the increase of Br. From Figure 8d, it can be seen that entropy reduction can be achieved by reducing the value of non-Newtonian parameter Λ. Finally, Figure 8e shows the entropy generation number falls with increasing Péclet number Pe. This is physically true due to the augment of Pe meaning a decrease of thermal conductivity of the fluid which can be seen in the definition of Pe and, therefore, a decreasing trend in entropy generation is observed [62]. Moreover, the influence of Br on the local entropy generation is illustrated in Figure 8c. It can be observed that the local entropy generation shows an increasing trend with the increase of Br. From Figure 8d, it can be seen that entropy reduction can be achieved by reducing the value of non-Newtonian parameter Λ. Finally, Figure 8e shows the entropy generation number falls with increasing Péclet number Pe. This is physically true due to the augment of Pe meaning a decrease of thermal conductivity of the fluid which can be seen in the definition of Pe and, therefore, a decreasing trend in entropy generation is observed [62].
In Figure 9, we can observe that the total entropy generation rates both show an increasing trend with the increase of Λ and Br for any fixed value of Ha, similar to the tendency of the local entropy generation SG. In particular, it can be noticed that the total entropy generation rate presents In Figure 9, we can observe that the total entropy generation rates both show an increasing trend with the increase of Λ and Br for any fixed value of Ha, similar to the tendency of the local entropy generation S G . In particular, it can be noticed that the total entropy generation rate presents a decreasing-increasing variation trend with Ha, which is consistent to the former conclusion in Figure 8a,b. a decreasing-increasing variation trend with Ha, which is consistent to the former conclusion in Figure 8a and b. Figure 10 illustrates total entropy generation rate versus magnetic field parameter Ha and the ratio of the viscous dissipation to the applied wall heat flux Br for different values of the lateral electric field S. There is a great decrease in the total entropy generation rate with the increasing value of S for given values of Ha and Br.

Conclusions
In this study, we have discussed a mathematical model for describing the electromagnetohydrodynamic flow and entropy generation of third-grade fluids between two parallel microplates combined with electroosmotic effects. The Chebyshev spectral collocation method has been applied to obtain the numerical solutions of the dimensionless velocity, temperature and entropy generation rate under the unidirectional flow assumption. The influences of dimensionless governing parameters, including non-Newtonian parameter (Λ), magnetic field (Ha) and viscous dissipation (Br) on the above obtained physical quantities are systematically investigated. The following conclusions can be drawn from the above theoretical analysis and numerical simulation. First of all, the results show that effect of non-Newtonian parameters are significant on fluid velocity, temperature and entropy generation rate. The dimensionless flow velocity and temperature are observed to decrease with the increase of non-Newtonian parameters and their maximum values have been reached when Λ = 0 (i.e., the fluid is Newtonian fluid). Under the combined action of electrical field and magnetic field, the variation profiles of velocity and  a decreasing-increasing variation trend with Ha, which is consistent to the former conclusion in Figure 8a and b. Figure 10 illustrates total entropy generation rate versus magnetic field parameter Ha and the ratio of the viscous dissipation to the applied wall heat flux Br for different values of the lateral electric field S. There is a great decrease in the total entropy generation rate with the increasing value of S for given values of Ha and Br.

Conclusions
In this study, we have discussed a mathematical model for describing the electromagnetohydrodynamic flow and entropy generation of third-grade fluids between two parallel microplates combined with electroosmotic effects. The Chebyshev spectral collocation method has been applied to obtain the numerical solutions of the dimensionless velocity, temperature and entropy generation rate under the unidirectional flow assumption. The influences of dimensionless governing parameters, including non-Newtonian parameter (Λ), magnetic field (Ha) and viscous dissipation (Br) on the above obtained physical quantities are systematically investigated. The following conclusions can be drawn from the above theoretical analysis and numerical simulation. First of all, the results show that effect of non-Newtonian parameters are significant on fluid velocity, temperature and entropy generation rate. The dimensionless flow velocity and temperature are observed to decrease with the increase of non-Newtonian parameters and their maximum values have been reached when Λ = 0 (i.e., the fluid is Newtonian fluid). Under the combined action of electrical field and magnetic field, the variation profiles of velocity and

Conclusions
In this study, we have discussed a mathematical model for describing the electromagnetohydrodynamic flow and entropy generation of third-grade fluids between two parallel microplates combined with electroosmotic effects. The Chebyshev spectral collocation method has been applied to obtain the numerical solutions of the dimensionless velocity, temperature and entropy generation rate under the unidirectional flow assumption. The influences of dimensionless governing parameters, including non-Newtonian parameter (Λ), magnetic field (Ha) and viscous dissipation (Br) on the above obtained physical quantities are systematically investigated. The following conclusions can be drawn from the above theoretical analysis and numerical simulation. First of all, the results show that effect of non-Newtonian parameters are significant on fluid velocity, temperature and entropy generation rate. The dimensionless flow velocity and temperature are observed to decrease with the increase of non-Newtonian parameters and their maximum values have been reached when Λ = 0 (i.e., the fluid is Newtonian fluid). Under the combined action of electrical field and magnetic field, the variation profiles of velocity and temperature with increasing Ha has been divided into two regions, separated by the critical Hartmann number Hac, which also increases with non-Newtonian parameter. In addition, we observe that the Nusselt number shows decreasing behavior for an increasing non-Newtonian parameter Λ. Finally, we find that non-Newtonian characteristic can stimulate both the local and the total entropy generation rate.