Entropy Analysis on Electro-Kinetically Modulated Peristaltic Propulsion of Magnetized Nanofluid Flow through a Microchannel

A theoretical and a mathematical model is presented to determine the entropy generation on electro-kinetically modulated peristaltic propulsion on the magnetized nanofluid flow through a microchannel with joule heating. The mathematical modeling is based on the energy, momentum, continuity, and entropy equation in the Cartesian coordinate system. The effects of viscous dissipation, heat absorption, magnetic field, and electrokinetic body force are also taken into account. The electric field terms are helpful to model the electrical potential terms by means of Poisson–Boltzmann equations, ionic Nernst–Planck equation, and Debye length approximation. A perturbation method has been applied to solve the coupled nonlinear partial differential equations and a series solution is obtained up to second order. The physical behavior of all the governing parameters is discussed for pressure rise, velocity profile, entropy profile, and temperature profile.


Introduction
During the recent few decades, Nanotechnology has received a remarkable attention due to its numerous applications in industry, since material containing nano-sized particles exhibits various chemical and unique physical features.Major applications include hybrid power engines, fuel cells, microelectronics, domestic refrigerators, pharmaceutical processes, chillers, cooling of nuclear reactors, space technology, boiler flue gas temperature control, and grinding.Nanofluid is the interaction of Nano-sized particles within the fluid.Nanofluid heat transfer is very much favorable to enhance the heat transfer and thermal conductivity of the fluid up to 40%.Nanofluid consists of nano-sized particles (diameter: 1-100 nm) in different Newtonian and non-Newtonian fluid such as water, ethylene glycol, and oil.In nanofluid, nanoparticles are made up of carbon or carbide nanotubes, oxides, and metals.Most of the fluids have low thermal conductivity and thus they are unable to match today's need of cooling rate.For this purpose, a new way to enhance the thermal conductivity of such type of fluids is by adding nano-sized particles in the base fluids (water) to form slurries.Nanofluid has become now the major attention of various researchers for the new manufacturing of automotive and plant cooling systems, and for transfer of heat in different heat exchanger devices.Choi [1] was the first who described experimentally the thermal conductivity of fluid with nanoparticles.Khanafer et al. [2] described the heat transfer enhancement in a two-dimensional enclosure using nanofluid under the influence of buoyancy force.Oztop et al. [3]  natural convection in a partially heated rectangular enclosure with nanofluid models.Sheremet and Pop [4] discussed the natural convection in the presence of sinusoidal temperature distributions on nanofluid flow propagating through a porous medium.He modeled the governing flow with the help of Darcy-Boussinesq approximation and the Buongiorno's nanofluid model.They used a second order finite difference method to obtain the numerical solutions.In this study, they found that the Sherwood number and Nusselt number behave as an increasing function due to greater influence of the buoyancy-ratio parameter, thermophoresis parameter, and Rayleigh number, and behave conversely for the Brownian motion parameter and Lewis number.
Xiao et al. [5] presented an analytical model for effective thermal conductivity of nanofluid with Brownian motion effects of nanoparticles.They found that the greatest influence of the concentration of nanoparticles is to enhance the thermal conductivity of nanofluids.Moreover, they also observed that the effective thermal conductivity of nanofluids for the smaller magnitude of nanoparticles is higher than the larger size at a given concentration.Sheremet et al. [6] discussed numerically the transient natural convection through a porous wave wall cavity filled with nanofluid.He used the mathematical model proposed by Tiwari and Das and utilized the Darcy-Boussinesq approximation with thermal dispersion.Cai et al. [7] highlighted the recent developments on fractal-based approaches to nanoparticle aggregation and nanofluids.Recently, Xia et al. [8] presented the optimization of a fractal-like architecture of porous medium associated with diffusivity, thermal conductivity, and permeability.They observed that the ratio of dimensionless permeability over dimensionless effective diffusivity diminished due to the decrement of tortuosity and porosity fractal dimensions.Hence, the lower tortuosity and porosity fractal dimensions were favorable to water/wind resistant fabric, as it is helpful to diminish the ratio of dimensionless permeability over dimensionless effective diffusivity.Further, they also found that the ratio of dimensionless total effective thermal conductivity over effective diffusivity enhances with tortuosity fractal dimension which reveals that the lower tortuosity fractal dimensional was helpful to clothing insulations, because it diminishes the ratio of dimensionless effective thermal conductivity over effective diffusivity.Some more pertinent studies on the said topic can be found from the references [9][10][11][12][13].
In a human body, symmetrical contraction and suspension of smooth muscles create a sinusoidal wave which is known as a peristaltic wave.It is a very well-known mechanism that helps in transporting various biological fluids, such as the transport of urine, in the human body.Chyme movement and vasomotion of small blood vessels are the common examples of peristaltic flow.Due to an extensive range of applications, peristaltic nanofluid flow grabbed the attention of different researchers.For instance, Bhatti et al. [14] discussed the peristaltic nanofluid flow in the presence of titanium nanoparticles with magnetic and endoscopic effects.He further obtained the series solutions using a homotopy perturbation method.Tripathi and Beg [15] presented a mathematical and theoretical study on nanofluids with peristalsis.They used the viscous fluid model and presented the exact solutions for the velocity, temperature, and nanoparticle concentration profiles.Akbar and Nadeem [16] considered PTT fluid model with nanoparticles propagating peristaltically in a diverging tube.Bhatti et al. [17] presented the two-phase flow simulation with heat and mass transfer in the presence of EDL ("Electric double layer") and MHD ("Magnetohydrodynamics") effects.A few more studies on the said topic can be found from the references [18][19][20].
Entropy generation can be defined as the phenomena of irreversibility processes of different thermal systems that are related to viscous dissipation, heat transfer, magnetic field and mass transfer etc.To enhance the performance of these irreversibilities, the second law of thermodynamics has been used [21][22][23][24][25][26] by various scientists.However, the first law of thermodynamics has also been used by various researchers and it is found that the second law of thermodynamic is more efficient and convenient as compared to the first law of thermodynamics.Rashidi et al. [27] investigated the entropy generation in steady magnetohydrodynamic flow of nanofluid over a rotating porous disk.Ranjit and Shit [28] recently investigated the entropy generation with a peristaltic flow in the presence of magnetic field.Akbar et al. [29] also discussed the entropy generation on peristaltic flow of water Entropy 2017, 19, 481 3 of 15 with copper nanoparticles.Sheikholeslami and Ganji [30] applied Lattice Boltzmann method (LBM) to determine the entropy generation on nanofluid flow with magnetic effects.Bhatti et al. [31,32] studied numerically the entropy generation with non-Newtonian fluid models through a shrinking and stretching sheet.
The purpose of the present study is to analyze the entropy generation on electro-kinetically modulated peristaltic propulsion on the magnetized nanofluid flow through a microchannel with joule heating.The effects of viscous dissipation, heat absorption, magnetic field, and electrokinetic body force are also taken into account.According to the best of our knowledge, no such analysis has been presented before.

Mathematical Formulation
The geometry model of the electro-osmotically peristaltic nanofluid flow through a microchannel is presented in Figure 1.The wavy wall surface of the microchannel is described as [33]: where h, a, c, t, φ, λ are the transverse vibration of the wall, half width of the channel, velocity of the wave, time period, amplitude of the wave, and wave length, respectively.The temperature at the wall and the center line are presented as The governing equations of motion, concentration, temperature, and continuity are expressed as [17]: With the help of Poisson Boltzmann equation, the electric potential with in a microchannel can be written as [17,33]: where ε is the permittivity, ρ e = ez[ n + − n − ] is the electrical charge density, e is the elementary charge, n + and n − are positive and negative charge ions consists of bulk concentration ("number density") n 0 and the valence of z + and z − respectively.Moreover, in order to analyze the potential distribution, it is compulsory to define the charge number density.For this purpose, Nernst-Planck equation is helpful to describe the number distributions of an individual species as [34].
In the above equation, D is the diffusivity of an ionic species, k B is the Boltzmann constant.The approximation of equal ionic diffusion coefficients for both species and the mobility of a species is described with the help of the Einstein formula.
Entropy 2017, 19, 481 4 of 15 In the above equation, D is the diffusivity of an ionic species, B k is the Boltzmann constant.
The approximation of equal ionic diffusion coefficients for both species and the mobility of a species is described with the help of the Einstein formula.To move from the fixed frame to the wave frame, we define the following variables as , , , .
Let us define the following non-dimensional quantities The nonlinear terms in Equation ( 9) are of ( ) O Pe δ , whereas

Sc
Pe Re = is the ionic Peclet number and Sc denotes the Schmidt number.Therefore, the nonlinear terms dropped in the limit i.e.,

, , 1
Re Pe δ << .The Poisson equation reduced into the following form where is debye length or characteristic thickness EDL ("electrical double layer").After some simplification Equation ( 10), the ionic distribution can be written as Along with following bulk conditions 0,and 0.
Then the Boltzmann distributions for the ions are To move from the fixed frame to the wave frame, we define the following variables as The nonlinear terms in Equation ( 9) are of O Pe δ 2 , whereas Pe = Sc Re is the ionic Peclet number and Sc denotes the Schmidt number.Therefore, the nonlinear terms dropped in the limit i.e., Re, Pe, δ << 1.The Poisson equation reduced into the following form where m = aez 2n 0 /k B ε T = a/λ d is an electro-osmotic parameter and λ d = 1/m is debye length or characteristic thickness EDL ("electrical double layer").After some simplification Equation ( 10), the ionic distribution can be written as Along with following bulk conditions ∂n ± ∂y = 0, and ∂Φ ∂y = 0. ( Then the Boltzmann distributions for the ions are Entropy 2017, 19, 481 5 of 15 Using Equations ( 11) and ( 14), we get To proceed further, Equation ( 15) must be modified by means of low-zeta potential approximation.This approximation is not beneficial due to a broad range of pH and the magnitude of the zeta potential is very small (less than 25 mV).Therefore, Equation (15) reduces to the following form Subject to whereas the potential function is defined as Using Equation (10) in Equations ( 1)-( 7), and applying the approximation long wavelength and creeping flow regime, ∂p ∂x subject to the following boundary conditions ∂u ∂y = 0, θ = 0; at y = 0 and u = 0, θ = 1; at y = h.

Entropy Generation Analysis
The volumetric rate of a local entropy generation is described in the following form as: Heat Transfer Irreversibility (HTI) The above equation describes the dimensional form of entropy generation due to fluid friction irreversibility, heat transfer, and magnetic field.A characteristics entropy generation is defined as Then the dimensionless form of entropy generation is written as Entropy 2017, 19, 481 where M, Br, Ω T are defined above.Furthermore, the integration of Equation ( 23) in the whole computational domains gives the dimensionless entropy generation, N avg G , which can be written as Moreover, the Bejan number B e is a significant parameter that describes the importance of irreversibility process in a domain and can be described as The relative global supremacy of heat transfer irreversibility can be analyzed by average Bejan number B avg e , which can be defined as It is noticeable here that B avg e < 1/2 reveals that the irreversibility process due to fluid friction dominates the flow whereas B avg e > 1/2 reveals that heat transfer irreversibility dominates and B avg e = 1/2 shows that the fluid friction and heat transfer entropy generation are equal.

Solution of the Problem
The exact solutions of Equations ( 19) and ( 21) are difficult to obtain due to a coupled nonlinearity.Therefore, we have applied the homotopy perturbation method [17,25] to obtain the solutions.Therefore, the homotopy for Equations ( 19) and ( 21) are defined as: where ξ is an embedding parameter.The Linear operator L is considered as and the initial guesses for the above linear operators are defined as Defining the following expansions Applying the series expansion in Equations ( 33)- (34) to the homotopy sets in Equations ( 29) and ( 30), then comparing the powers of ξ, the system of linear differential equations can be obtained.Using the property of homotopy perturbation scheme [17,25] i.e., ξ → 1 , we have The final form of solutions for the velocity profile and temperature can be written in simplified form as u(y) = u 0 + u 1 y 2 + u 2 y 3 + u 3 y 4 + u 4 y 5 + u 5 y 6 + u 6 y 8 + u 7 cosh my, (37) whereas the constant appearing i.e., (u i , θ i ; i = 1, 2, 3....) in the above equations can be found using routine calculations.The instantaneous volume flow rate Q can be calculated as The mean flow rate over one period of wavelength can be written as The pressure rise is evaluated numerically with the help of the following expression

Results and Discussion
The main idea of this section is to provide the physical behavior of all the governing parameters that are included in the present flow problem.To determine graphically, a symbolic computational software "Mathematica" has been utilized.For this purpose, Figure 2 to Figure 12 are sketched against velocity profile, temperature profile, Bejan number.and entropy profile for the following parameters i.e., Hartmann number M, thermal Grashof number Gr, electro-osmotic parameter m, Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") U HS , Joule heating parameter S, heat source/sink β, amplitude ratio φ, and Brinkmann number Br. Pumping features have also been evaluated numerically with the help of Equation ( 33) and are discussed.Figure 2 shows the numerical comparison with previously published results [35] by ignoring the magnetic field (M → 0) and the viscous dissipation effects in the governing equations.From this figure, we can see that the present results are in excellent agreement with previously published data.Therefore, it is concluded that the methodology used to obtain the results are correct and the present results converge.In the all the upcoming figures, the red, blue, and green colors reveal the variation of the mentioned variable in ascending order.However, the solid and dashed lines indicate the lower to the higher value for the given parameter.Figure 3 represents the variation of the velocity profile against different values of thermal Grashof number Gr and Hartmann number M. In this figure we can see that the thermal Grashof number significantly enhances the velocity profile in the region y ∈ [0.5, 1], however, the results are completely converse when y < 0.5.On the other hand, the effects of Hartmann number M is completely opposite as compared to thermal Grashof number.The velocity of the fluid tends to diminish when y > 0.5, because when the magnetic field is applied to any conducting fluid, a converse force occurs which tends to oppose the velocity and provide a marked resistance in the flow field.Figure 3 is plotted for electro-osmotic parameter m, and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") U HS .In Figure 4, we can see that the electro-osmotic parameter accelerates the fluid velocity in the region y ∈ [0.5, 1], whereas its behavior becomes the opposite for y < 0.5.Electro-osmotic parameter is the ratio of channel length and a Debye length m(= a/λ d ) which reveals that the increment in the height of the channel produces a significant enhancement in velocity profile and consequently and increment in electrical Debye length.Further, Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") U HS (= −E x ζε/µ c) also boosts the velocity profile in the similar region as compared to the electro-osmotic parameter.S significantly raises the temperature profile.The Joule heating parameter S is a normalized generation term that shows the ratio of Joule heating to the surface heat flux.The result is cooling in the electrokinetic flow through a microchannel with a negative S and heating with a positive S .Further, Joule heating is related to the electric field induced by the internal heating process in a fluid in the presence of an applied potential gradient.Higher values of an electric field cause a higher heat generation which results in heating of the fluid.In Figure 5 we can further notice that higher values of heat source/sink , β  parameter S is a normalized generation term that shows the ratio of Joule heating to the surface heat flux.The result is cooling in the electrokinetic flow through a microchannel with a negative S and heating with a positive S. Further, Joule heating is related to the electric field induced by the internal heating process in a fluid in the presence of an applied potential gradient.Higher values of an electric field cause a higher heat generation which results in heating of the fluid.In Figure 5 we can further notice that higher values of heat source/sink β, produce a positive effect on the temperature profile and significantly boosts the temperature profile.Figure 6 shows that the effects of magnetic field and Brinkmann number accelerate the temperature profile.In fact, an increment in Brinkmann number tends to diminish the conduction of heat that occurs due to viscous dissipation and, as a result, the temperature rises significantly.S significantly raises the temperature profile.The Joule heating parameter S is a normalized generation term that shows the ratio of Joule heating to the surface heat flux.The result is cooling in the electrokinetic flow through a microchannel with a negative S and heating with a positive S .Further, Joule heating is related to the electric field induced by the internal heating process in a fluid in the presence of an applied potential gradient.Higher values of an electric field cause a higher heat generation which results in heating of the fluid.In Figure 5 we can further notice that higher values of heat source/sink , β produce a positive effect on the temperature profile and significantly boosts the temperature profile.Figure 6 shows that the effects of magnetic field and Brinkmann number accelerate the temperature profile.In fact, an increment in Brinkmann number tends to diminish the conduction of heat that occurs due to viscous dissipation and, as a result, the temperature rises significantly. .It can be viewed from Figure 8 that an enhancement in the electro-osmotic parameter , m creates a noticeable increment in the peristaltic pumping region and retrograde pumping, whereas its attitude is completely opposite in the co-pumping region.Further we can see that amplitude ratio φ acts similarly to an electro-osmotic parameter in all the regions but it has a greater impact on pressure rise.m, creates a noticeable increment in the peristaltic pumping region and retrograde pumping, whereas its attitude is completely opposite in the co-pumping region.Further we can see that amplitude ratio φ acts similarly to an electro-osmotic parameter in all the regions but it has a greater impact on pressure rise.
pumping regions.It is a very important and substantial part of the peristaltic flow that is beneficial in the propagation of multiple fluids.It can be observed from Figure 7 that the Hartmann number , M and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") HS U behaves in a similar way on pressure rise.For instance, with the increment of both the parameters , HS U M, the pumping rate increases in the retrograde pumping region 0, 0 p Q Δ > <    , whereas the behavior of pumping rate changes in peristaltic pumping 0, 0 .
 and becomes opposite in co-pumping region 0, 0 p Q Δ < >    .It can be viewed from Figure 8 that an enhancement in the electro-osmotic parameter , m creates a noticeable increment in the peristaltic pumping region and retrograde pumping, whereas its attitude is completely opposite in the co-pumping region.Further we can see that amplitude ratio φ acts similarly to an electro-osmotic parameter in all the regions but it has a greater impact on pressure rise.

Conclusions
The present study deals with the theoretical study of entropy generation on electro-kinetically modulated peristaltic propulsion on magnetized nanofluid flow through a microchannel with joule heating.The effects of heat absorption, viscous dissipation, electrokinetic body force, and magnetic field are also considered in the present flow.Series solutions have been obtained for the velocity profile and temperature profile, whereas numerical integration has been used to evaluate the expression for pressure.The major findings for the governing flow problem are as follows: i.
The velocity of the fluid behaves in a similar manner due to the increment in thermal Grashof number and magnetic field.ii.
The velocity profile increases when y > 0.5 for higher values of electro-osmotic parameter and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity").iii.
The temperature profile reveals a significant increment due to the higher presence of Brinkmann number and magnetic field.iv.
An increment in the Joule heating parameter and heat source/sink significantly accelerates the temperature profile.v.
Pressure rise exhibits similar behavior to higher values of an electro-osmotic parameter, magnetic field, and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity").vi.
The entropy profile also shows a positive response due to the greater impact of Brinkmann number, magnetic field, and the Joule heating parameter.

Figure 1 .
Figure 1.Flow structure of the electro-osmotically modulated peristaltic flow.

Figure 1 .
Figure 1.Flow structure of the electro-osmotically modulated peristaltic flow.

Figure 3
Figure 3 represents the variation of the velocity profile against different values of thermal Grashof number Gr and Hartmann number .M In this figure we can see that the thermal Grashof number significantly enhances the velocity profile in the region 0.5,1 y ∈     , however, the results are completely converse when 0.5.y < On the other hand, the effects of Hartmann number M is completely opposite as compared to thermal Grashof number.The velocity of the fluid tends to diminish when 0.5, y > because when the magnetic field is applied to any conducting fluid, a converse force occurs which tends to oppose the velocity and provide a marked resistance in the

Figure 2 .
Figure 2. Comparison of velocity distribution with previously published results [35].

Figure 2 .
Figure 2. Comparison of velocity distribution with previously published results [35].

Figure 3
Figure 3 represents the variation of the velocity profile against different values of thermal Grashof number Gr and Hartmann number .M In this figure we can see that the thermal Grashof number significantly enhances the velocity profile in the region 0.5,1 y ∈     , however, the results are completely converse when 0.5.y <On the other hand, the effects of Hartmann number M is completely opposite as compared to thermal Grashof number.The velocity of the fluid tends to diminish when 0.5, y > because when the magnetic field is applied to any conducting fluid, a converse force occurs which tends to oppose the velocity and provide a marked resistance in the flow field.Figure3is plotted for electro-osmotic parameter , m and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") .HS UIn Figure4, we can see that the electro-osmotic parameter accelerates the fluid velocity in the region 0.5,1 y ∈     , whereas its behavior becomes the opposite for 0.5.y < Electro-osmotic parameter is the ratio of channel length and a Debye length ( ) / d m a λ =  which reveals that the increment in the height of the channel produces a significant enhancement in velocity profile and consequently and increment in electrical Debye length.Further, Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") ( ) / HS x c E U ζε μ = −  also boosts the velocity profile in the similar region as compared to the electro-osmotic parameter.

Figure 3 .
Figure 3. Velocity distribution against multiple values of M and Gr .Figure 3. Velocity distribution against multiple values of M and Gr.

Figure 3 . 15 Figure 4 .
Figure 3. Velocity distribution against multiple values of M and Gr .Figure 3. Velocity distribution against multiple values of M and Gr.Entropy 2017, 19, 481 9 of 15

Figures 5
Figures 5 and 6 represent the temperature profile against different values of Hartmann number , M Joule heating parameter , S heat source/sink , β and Brinkmann number .Br It can be noticed from Figure 4 that Joule heating parameter ,S significantly raises the temperature profile.The Joule heating parameter S is a normalized generation term that shows the ratio of Joule heating to the surface heat flux.The result is cooling in the electrokinetic flow through a microchannel with a negative S and heating with a positive S .Further, Joule heating is related to the electric field induced by the internal heating process in a fluid in the presence of an applied potential gradient.Higher values of an electric field cause a higher heat generation which results in heating of the fluid.In Figure5we can further notice that higher values of heat source/sink , β

Figure 4 .
Figure 4. Velocity distribution against multiple values of U HS and m.

Figures 5
Figures 5 and 6 represent the temperature profile against different values of Hartmann number M, Joule heating parameter S, heat source/sink β, and Brinkmann number Br.It can be noticed from Figure 4 that Joule heating parameter S, significantly raises the temperature profile.The Joule heating

Figure 4 .
Figure 4. Velocity distribution against multiple values of HS U and m .

Figures 7 and 8
Figures 7 and 8 are sketched to analyze the behavior of peristaltic pumping in different pumping regions.It is a very important and substantial part of the peristaltic flow that is beneficial in the propagation of multiple fluids.It can be observed from Figure 7 that the Hartmann number , M and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") HS U behaves in a similar way on pressure rise.For instance, with the increment of both the parameters , HS U M, the pumping rate increases in the retrograde pumping region 0, 0 p Q Δ > <    , whereas the behavior of pumping rate changes in peristaltic pumping 0, 0 .p Q Δ > >    and becomes opposite in co-pumping region 0, 0 p Q Δ < >    .It can be viewed from Figure 8 that an enhancement in the

Figure 6 .
Figure 6.Temperature distribution against multiple values of M and Br.

Figures
Figures 7 and 8 are sketched to analyze the behavior of peristaltic pumping in different pumping regions.It is a very important and substantial part of the peristaltic flow that is beneficial in the propagation of multiple fluids.It can be observed from Figure 7 that the Hartmann number M, and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") U HS behaves in a similar way on pressure rise.For instance, with the increment of both the parameters U HS , M, the pumping rate increases in the retrograde pumping region [∆p > 0, Q < 0], whereas the behavior of pumping rate changes in peristaltic pumping [∆p > 0, Q > 0]. and becomes opposite in co-pumping region [∆p < 0, Q > 0].It can be viewed from Figure 8 that an enhancement in the electro-osmotic parameter Figures 7 and 8 are sketched to analyze the behavior of peristaltic pumping in different pumping regions.It is a very important and substantial part of the peristaltic flow that is beneficial in the propagation of multiple fluids.It can be observed from Figure 7 that the Hartmann number M, and Helmholtz-Smoluchowski velocity (or "maximum electro-osmotic velocity") U HS behaves in a similar way on pressure rise.For instance, with the increment of both the parameters U HS , M, the pumping rate increases in the retrograde pumping region [∆p > 0, Q < 0], whereas the behavior of pumping rate changes in peristaltic pumping [∆p > 0, Q > 0]. and becomes opposite in co-pumping region [∆p < 0, Q > 0].It can be viewed from Figure 8 that an enhancement in the electro-osmotic parameter

Figure 7 .
Figure 7. Pressure distribution against multiple values of M and HS U .Figure 7. Pressure distribution against multiple values of M and U HS .

Figure 7 . 15 Figure 8 .
Figure 7. Pressure distribution against multiple values of M and HS U .Figure 7. Pressure distribution against multiple values of M and U HS .Entropy 2017, 19, 481 11 of 15

Figures 9 and 10 ΩFigure 8 .
Figures 9 and 10 are sketched to determine the behavior of entropy profile against different values of Hartmann number , M temperature difference , T Ω Joule heating parameter , S and Brinkmann number .Br Figure 9 shows that an increment in Hartmann number M causes a significant increment in the entropy profile and remains positive throughout the region, whereas the parameter T Ω shows opposite attitude on the entropy profile.It can be seen from Figure 10 that when the Joule heating parameter , S and Brinkmann number Br increase, then the entropy profile greatly increases and it remains positive throughout the region.Figures 11 and 12 show the Bejan number variation along multiple values of Hartmann number , M temperature difference , T Ω Joule heating parameter , S and Brinkmann number .Br From both the figures, one can easily observe that an increment in Hartmann number , M and Joule heating parameter , S and Brinkmann number Br tends to diminish the Bejan number significantly, however, converse behavior has been observed against the temperature difference .T Ω

Figures 9 and 10
Figures 9 and 10 are sketched to determine the behavior of entropy profile against different values of Hartmann number M, temperature difference Ω T , Joule heating parameter S, and Brinkmann number Br.Figure9shows that an increment in Hartmann number M causes a significant increment in the entropy profile and remains positive throughout the region, whereas the parameter Ω T shows opposite attitude on the entropy profile.It can be seen from Figure10that when the Joule heating parameter S, and Brinkmann number Br increase, then the entropy profile greatly increases and it remains positive throughout the region.Figures11 and 12show the Bejan number variation along multiple values of Hartmann number M, temperature difference Ω T , Joule heating parameter S, and Brinkmann number Br.From both the figures, one can easily observe that an increment in Hartmann number M, and Joule heating parameter S, and Brinkmann number Br tends to diminish the Bejan number significantly, however, converse behavior has been observed against the temperature difference Ω T .

Figure 9 ΩΩFigure 9 .
Figures 9 and 10 are sketched to determine the behavior of entropy profile against different values of Hartmann number M, temperature difference Ω T , Joule heating parameter S, and Brinkmann number Br.Figure9shows that an increment in Hartmann number M causes a significant increment in the entropy profile and remains positive throughout the region, whereas the parameter Ω T shows opposite attitude on the entropy profile.It can be seen from Figure10that when the Joule heating parameter S, and Brinkmann number Br increase, then the entropy profile greatly increases and it remains positive throughout the region.Figures11 and 12show the Bejan number variation along multiple values of Hartmann number M, temperature difference Ω T , Joule heating parameter S, and Brinkmann number Br.From both the figures, one can easily observe that an increment in Hartmann number M, and Joule heating parameter S, and Brinkmann number Br tends to diminish the Bejan number significantly, however, converse behavior has been observed against the temperature difference Ω T .

Figure 11 .
Figure 11.Variation of Bejan number against multiple values of M and Br .

Figure 10 . 15 Figure 10 .
Figure 10.Entropy profile against multiple values of S and Br.

Figure 11 .
Figure 11.Variation of Bejan number against multiple values of M and Br .Figure 11.Variation of Bejan number against multiple values of M and Br.

Figure 11 .
Figure 11.Variation of Bejan number against multiple values of M and Br .Figure 11.Variation of Bejan number against multiple values of M and Br.

Figure 11 .
Figure 11.Variation of Bejan number against multiple values of M and Br .

Figure 12 .
Figure 12.Variation of Bejan number against multiple values of S and T Ω .Figure 12. Variation of Bejan number against multiple values of S and Ω T .

Figure 12 .
Figure 12.Variation of Bejan number against multiple values of S and T Ω .Figure 12. Variation of Bejan number against multiple values of S and Ω T .