This research deals with the mathematical model for the development of the peristaltic principle of the combination of the pressure and electroosmotic flow (EOF) of ionic liquid across microchannels with electrokinetic effects. For thermomechanical dynamics, the convective conditions on the boundary for mass and heat transfer at the walls of the channel are quantified. For the microchannel, a porous structure is presumed. Soret, Dufour, and Joule heating are also listed in the scope of the problem addressed. The corresponding equations for the ionic fluid flow, mass, and heat transfer along with the Poisson–Boltzmann equation within the electrical double layer (EDL) are studied. The exact solution has been obtained based on lubrication theory (i.e., low Reynolds number and long wavelength approximations). The channel height is therefore believed to be much higher than the electrical double layer (EDL) thickness. Various dimensionless pertinent parameters illustrate the important aspects of electroosmotically controlled flow and subsequent convective mass/heat transfer attributes in a microchannel. A linear dependency on the fluid flow rate is exhibited by the pressure drop. The analysis shows that the electroosmotic parameter gives a reducing effect on the channel permeability. The distribution of temperature and concentration is greatly affected by convective heat and mass parameters, respectively. In biomedical engineering, the application areas of the study proposed are for the design of the devices such as a microfluidic pump to pump a small amount of ionic liquids by regulating the variation in temperature and concentration.

1. Introduction

Electroosmosis has been seen to a considerable extent in recent decades. Microelectromechanical System (MEMS) was the first device with microfluidic technology assembled in the 1980s. The installations of these devices are widespread in a range of contexts comprising biological and biomedical industry sectors, where microfluidic products are termed as LOC (lab-on-chip devices). These devices can also be used for a variety of purposes like biomedical therapeutics and microbial or toxic contamination. Nowadays, a strong demand is raised on compact, simple, and inexpensive standardized equipment for pharmaceutical dissemination and DNA-based analysis because of the need for effective, rapid, and considerable standardized reliability of microfluidics. Microfluidic devices frequently drive fluid throughout the whole network employing electroosmosis. Electroosmosis represents the electrolyte flow across a channel with an applied voltage charge boundary. In several biological, clinical, and manufacturing systems, electroosmosis develops, including dialysis of fluids, porous membranes, zoological processes, transportation in human skin, tube/canal flow, and techniques of separation. The simulation of EOF was projected by Patankar and Hu [1]. In capillary assessments, Huang et al. [2] demonstrated the electroosmotic flow (EOF). The composition characteristics of microflows were examined on an EOF basis by Haswell [3]. The applications of microfluidics have been addressed by Gravesen et al. [4]. In the annulus of capillary walls, Kang et al. [5] also analyzed the EOF.

Investigations on controlled microchannel flow continued to draw interest, leading to its widespread use in various fields of science, especially genetics, medicine, and medical technology. Studies such as this one have shown, on the one hand, that they are capable of revealing interesting characteristics, especially if such channels reflect a longitudinal wave that travels across the length of walls, in this scenario called the peristalsis. The peristalsis is therefore characterized as a wavelike movement inside flexible tubes triggered and operated by a physiological flexion and extension of adjacent locations, which leads to a movement of the within contents. We also note that since this mechanism is becoming increasingly relevant in many ways, it has recently undergone comprehensive study with an even wider impact. Considering its relevance to various industrial and physiological fields, many observations were taken in the studies [68]. Furthermore, several biomedical electronic types of equipment, for example, the cardiovascular device and plasma pumps during dialysis, have found this spontaneous phenomenon. However, after one of Latham’s first major investigations [9], several researchers had already been motivated to explore the dynamics of peristaltic motion. At that time, there have been numerical [68, 10] and experimental [1113] approaches used, where there have been those who have used analytical [14, 15] approach to examine and understand the effects of such oscillatory movement’s mechanics and physiology. The way that this type of action takes place is by the squeezed displacement of the walls of the channel and their contents. A good example is seen as food is swallowed via the esophagus. Thus, peristaltic pumps are usually demanded in pharmaceutical companies due to their systemic potential for transportation, thus generating considerable interest in transportation studies as they have been done in [13].

In addition, peristaltic pumping is not infrequent to be subject to other external consequences which can seriously impair it. This mechanism has been extensively investigated and is often used with effects such as heat transfer [8], electroosmosis [16], and magnetic fields [7]. While this type of analysis usually results in nonlinear procedures, a relatively small number has been designed to give rise to linear forms. Consequently, several mathematical models have thus been explored governing peristaltic flow. The potential effects of peristaltic pumping and common external effects (electromagnetic, electroosmotic) [1723], which have been documented in the literature, influence particularly physiological flows. The narrowness of the peristaltic conduits has made it highly appealing for the lubrication assumption. Despite its success, other researchers have prevented this approach, but they have given the solutions for these dynamics [68]. In addition, the studied results alongside peristalsis highlight the non-Newtonian parameters, the thickness of the electrical double layer (abbreviated as EDL), the external electromagnetic field, slip parameter, Joule heating, flow characteristics, and radiation. Surprisingly, there was also an overview of entropy generation analysis [18]. Indeed how this impact alters peristaltic transports significantly is interesting.

In terms of tissue degradation, portable disease testing kits, analysis of dilution strategies in blood flow, and microfabrication technology, heat transfer in microfluidic systems is extremely important. The effect of heat exchange on the third-order fluid in peristaltically driven EOF is investigated by Waheed et al. [24] in EOF. The heat flow impact of the power-law fluid on the microchannel is reported by Babaie et al. [25]. Besides, Chen et al. [26] have studied the thermal exchange of the non-Newtonian suspension in a microchannel. Sinha and Shit [27] were engaged in thermal blood exchange across capillaries. EOF for heat exchange and MHD also was shown by Shit et al. [28]. Also, the tapered walls in the EOF which improve transfer of mass are regarded in medical operations. Therefore, the heat and mass distribution through EDL effects have been represented by Bhatti et al. [29]. Furthermore, the MHD electroosmotic flow of power-law PVC nanofluid in a channel was discovered by Shehzad et al. [30]. Yadav et al. [31] also analyzed peristaltic propelling with EDL. The influence of heat in the EOF was investigated by Narla et al. [32] through time. Furthermore, the transfer of heat in MHD EOF was also portrayed by Yang et al. [33].

Both mass and heat transfer impacts play a vital role in the distribution of organic compounds in saturated soil, radioactive waste disposal, thermal insulation, increased oil recovery, and so on. In combustion, fermentation, and distillation and many other manufacturing processes, the impact of the mass transfer appeared. In the case of devices like steam engines, nuclear plants, heat storage, and so forth, convective heat transfer via boundary conditions is used. The results of the heat/mass transfer on the fluid triggering peristaltic motion in fully compliant curved canals were investigated, in this regard, by Hina and Hayat [34]. The characteristics of heat flux in peristaltic transportation with/without compliant walls were studied by Riaz et al. [35], Hayat et al. [36], and Mehmood et al. [37]. In the presence of convective constraints, different non-Newtonian fluid models with/without peristalsis are analyzed in [3843] and many references therein. Yasmin et al. [44] addressed the impact of convective conditions in the peristalsis of Johnson–Segalman fluid in an asymmetric channel.

With this development, we notice that the solutions may lead to complicated mathematical procedures when the aforesaid nontrivial hypotheses are implemented into the peristaltic microchannels. Analytical solutions have been achieved [16, 45] for basic situations in which the leading governing equations are compatible. The rheological effects of fluid were also studied on the electroosmotic flow in a non-Darcy porous medium [46, 47]. In other cases, the complexity involved the semianalytical approach, which included the use of numerical schemes [4850] and perturbation theory [51] to find a complete solution to the problem. However, some effective research without the use of perturbation theory [68] was observed in other cases.

We realize that no study on convective heat and mass transfer in EOF altered by peristalsis with Soret, Dufour, Joule heating, and porous medium currently exists. This research covers the void. The main objective of this paper is to study an ionic liquid with viscous dissipation in a microchannel with a porous structure, impact of convective heat transfer, and convective mass transfer for EOF. The equations governing the fluid velocity, mass, and heat transfer in the EDL are taken with the famous Poisson–Boltzmann equation. Approximation of lubrication and Debye-Hückel linearization are implemented for this research. The exact solution for the governing equations is obtained. Schematic findings are drawn and discussed thoroughly for physical quantities.

2. Flow Analysis

Consider a fully developed unsteady ionic liquid flow via a microchannel bounded by two walls with rhythmic contraction, , in accordance with an alternating external electric field and the enacted pressure gradient. Due to wave propagation, the constituent equation for deforming of the wall (see Figure 1) takes the formwhere the axial channel coordinate, time, half-width of the channel, the amplitude of the wave, the wavelength, and the wave velocity are represented by , and respectively. It is presumed that the flow is symmetrical along the channel center , which means that the flow domain in the area is studied. At the walls of the channel, the zeta potential is taken as constant. produces both an applied external electric field alongside the microchannel axis and the required electroosmotic fluid flow driving force as well as the peristaltic wave-induced flow. The no-slip conditions for velocity have been taken into consideration, whereas convective heat along with convective mass conditions is applied in this analysis.

Because of the existence of electrical double layer (EDL) in the microchannel as defined by the Poisson–Boltzmann equation [52], for a symmetric binary electrolyte solution, the electric potential distribution developed in the Cartesian coordinate can be written asin which is the dielectric constant of the medium, represents the density of the ionic change, and is the vacuum permittivity which is equal to (density of the ionic energy) for a symmetric electrolyte is written asin which is the number of anions’ density and is the number of cations’ density and these can be given by Boltzmann distribution with the consideration of no overlapped EDL:where is the bulk ions concentration, which is independent of electrochemistry surface, is the Boltzmann constant, represents the electric charge, and denotes the average electrolytic solution temperature.

This distribution of the ionic concentration tends to be true if the ionic concentration in the microchannel is not axial and, for fluid flow, the substantially small Peclet number is supposed. From equations (2)–(4), one can get the Poisson–Boltzmann equation as follows:

Introduce now a standardized electroosmotic potential function with medium potential and other dimensionless variables:

Equation (5) is now read in view of the nondimensionally defined variables from equation (6) asin which the ionic energy parameter is defined as and the Debye length is denoted by . As the potential function is symmetrical, the following condition is taken into account:and the potential function at the borderline should be

Further, it is assumed that the wall zeta potential is insufficiently small to enable the approximation of the Debye-Hückel linearization.

The linearized Poisson-Boltzmann equation (7) can be resolved in order to achieve the potential distribution by using the boundary conditions (8) and (9) aswhere the electroosmotic parameter is represented by . The study of the presented model is conducted under the impact of all the above assumptions; the biofluid flow equations and the heat transfer via permeable microchannel are considered as [52, 53]where the velocity components along x-axis and y-axis are denoted as and , respectively, the fluid density is indicated by is the viscosity, is the pressure, is the specific heat with constant pressure, T is the ionic fluid temperature, with the thermal conductivity being k1, the permeability of the medium is k, Ex is the electrical field in the x direction, σ is the electrical conductivity, C is the fluid concentration, kT is thermal diffusion ratio, D is the coefficient of concentration diffusivity, cs is the susceptibility of concentration, and Tm is the mean temperature. The first term on the left-hand side of each equation from (12) to (14) shows the transient state. The terms and on the right-hand side of equations (12) and (13) show the permeability of the microchannel, whereas the first term on the right-hand side of equations (12) and (13) represents the applied pressure gradient. The term of equation (14) provides the effects of Joule heating. In this study, the Joule heating with viscous dissipation in a microchannel is considered. The last term in equation (14) represents the Soret effects, whereas the last term in equation (15) represents the Dufour effects.

The conditions for the particular problem are defined aswhere ht and hm represent the convective heat and mass transfer coefficients, respectively, Ta is the ambient temperature at the walls, and Ca is the ambient concentration at the walls. The conditions on heat and mass transfer are called convective boundary conditions. The source of these boundary conditions is given in [40, 41, 44] and many references therein.

Let us now establish the following dimensionless variable quantities:where p is the dimensionless pressure, Re is the Reynolds number, δ is the wave number, k is the permeability parameter, Pr is the Prandtl number, γ is defined as the Joule heating parameter, Br is the Brinkman number, Du, Sc, and Sr are the Dufour, Schmidt, and Soret numbers, respectively, and Bit and Bim indicate the Biot numbers for heat and mass transfer, respectively. Introducing these dimensionless variables into equations (11)–(15), the nondimensional governing equations may be put aswhere denotes the medium mobility and represents the maximum electroosmotic or Helmholtz–Smoluchowski velocity defined by . By following the approximations of long wavelength and low Reynolds number, equations (18)–(22) correspond to

In view of the problem’s physical conditions, the constraints can be mathematically expressed as

Let us introduce the stream function , given by and , which satisfies the continuity equation (18) automatically. Equation (19) can be written as follows:

After eliminating the pressure term, one can obtain

The boundary conditions in terms of stream function can be written as

It is worth mentioning here that, for obtaining the solution of ODE in as it exists in equation (26) after removal of the pressure term, two additional boundary conditions are placed on the stream function. Dimensionless flow rate F is supposed to vary exponentially with constant A and Q0 given by . Constant is accountable for the negative and positive flow rate or according to or , respectively. When the flow is contrary to the peristaltic movement direction, the negative flow rate reveals retrograde or reverse pumping. The positive pumping is, however, the case where both pressure rise and volumetric flow rate are positive. The blood flow rate reduces exponentially over time, experimentally demonstrated by Kikuchi [54]. This variance in flow rates was also proposed to be not dependent on the microchannel’s structural information.

3. Analytical Solutions

Solving the system of equations subject to the relevant boundary conditions yields the following solutions for stream function and temperature. The following solutions exist for stream function, temperature, and concentration profiles for the system according to the required boundary conditions:where

The results presented in the section are in good agreement with the previous study by Ranjit et al. [17] in the absence of convective heat/mass transfer conditions.

4. Results and Discussion

Analytical solutions have been acquired and reported in the previous sections on the axial velocity, stream function, pressure gradient, temperature, and concentration characteristics. We provided a series of figures in this section, describing the effects of different parameters of interest on flow quantities such as the axial velocity , distribution of temperature , and concentration by setting the parametric values , , , , , , , , , , , , , , and . Figures 24 demonstrate a change in axial velocity for several parameters. The consequences on the axial velocity of said electroosmotic parameter m can be seen in Figure 2. It reveals that the velocity profile depends highly on m through the channel. The axial fluid velocity is increased when the values of the electroosmotic parameter m rise in the neighborhood of the channel's wall. Because m is the ratio between channel height y = h and the Debye thickness λD, an increasing trend in channel height results in an increase of the axial fluid velocity. The influence of permeability k on fluid velocity is illustrated in Figure 3. This graph represents that a boost in the value of the permeability parameter k raises the fluid velocity in the middle of the channel, while the fluid speed decreases across the walls of the channel. From Figure 4, we noted that, with decreasing flow rate F, the axial velocity tends to decrease. We noted that the reduced volumetric flow rate F indicates that a phenomenon called reflux flow occurs in reverse flow. A reflux flow (that occurs for the critical flow rate) corresponding to the preference of the physical parameters is therefore very straightforward to quantify. Figure 5 highlights the results of the fluid velocity profile for the mobility of the porous medium β. It indicates that the axial velocity lowers when the value of β grows. As β depends directly on Helmholtz–Smoluchowski velocity , this can be perceived physically as follows: with increased EDL thickness, the fluid's velocity declines, and the fluid flow decreases when EDL is present. In Figures 614, the impact on the fluid temperature of the relevant parameters is shown. Figure 6 shows the influence of the mobility of the medium β. It reveals that the temperature of the liquid decreases with the greater mobility of the medium. As β depends directly on Helmholtz–Smoluchowski velocity , this can be perceived physically as follows: with increased EDL thickness, the fluid’s velocity declines, and the fluid flow decreases when EDL is present; hence, it decreases the fluid temperature. To study the effect of thermal Biot number Bit on fluid temperature, Figure 7 is portrayed. It can be seen that the fluid temperature decays as the value of the thermal Biot number is increased. Figure 8 depicts the effect of the Brinkman number Br on the temperature profile. The fluid temperature rises when the Brinkman number is starting to rise. The principal truth in Br is the resistance creating viscosity. This resistance induces a collision between fluid particles, which causes a rise in temperature in the resulting collision. Increased fluid temperature is observed in Figure 9 due to a rise in the value of the Dufour number Du. For the improved value of the permeability parameter k, Figure 10 indicates an increased fluid temperature. Figure 11 shows the temperature profile difference for various electroosmotic parameter (m) values. The improvement in the value of the permeability parameter increasing the fluid temperature is apparent from this figure. The fluid temperature is increased by an improvement of the value of Schmidt Sc in Figure 12. From Figure 13, it is evident that the fluid temperature is elevated with the rapid increase of the value of Soret number Sr. The effect of the Joule heating parameter γ depicts an increase in the fluid temperature as portrayed in Figure 14. The effect of the relevant parameters on fluid concentration is reported for analysis in Figures 1523. Figure 15 indicates that the distribution of the concentration declines with a rise in the mobility of medium β. In Figure 16, there is an upturn in the concentration profile for increasing the value of the Biot number Bim (mass transfer). For the Brinkman number Br and Dufour number Du, Figures 17 and 18 reflect the concentration distribution . From these figures., it can be noted that the concentration distribution declines as Brinkman number Br and Dufour number Du are increased. As the value of the permeability parameter increases, Figure 19 indicates an increase in fluid concentration. Decay in the concentration profile is seen when the value of an electroosmotic parameter is increased (see Figure 20). For the Schmidt Sc and Soret Sr numbers in Figures 21 and 22, the same behavior of fluid concentration is found, meaning that the fluid concentration decays as the Schmidt or Soret numbers rise. The effect of the Joule heating parameter depicts a decrease in the fluid concentration as shown in Figure 23.

5. Concluding Remarks

This article addresses the electroosmotic flow (EOF) with the impact of Joule heating, Soret, Dufour, and permeability in the peristaltic motion of an ionic fluid. The convective heat and mass conditions are employed at the microchannel walls. The major observations have been summarized as follows:(i)The axial velocity is increased as the values of the electroosmotic parameter m rise in the neighborhood of the channel’s wall.(ii)A boost in the value of the permeability parameter k raises the fluid velocity in the middle of the channel, while the fluid speed decreases across the walls of the channel.(iii)The axial velocity and ionic fluid temperature decay when the value of mobility of medium β grows.(iv)The fluid temperature decays as the value of the thermal Biot number Bim is increased, whereas there is an upturn in the concentration profile for increasing the mass transfer Biot number Bim.(v)An increase in the temperature and concentration is noticed with a rise in the value of the permeability parameter k.(vi)The fluid temperature is elevated with the rapid increase of the value of the Soret number Sr, while reverse behavior is detected for the concentration of the fluid.(vii)In the presence of the Joule heating effect γ, the temperature and concentration distributions are reversed. Hence, it seems important not to neglect the Joule heating effect in the temperature equation to obtain more accurate results.


d:Mean half channel width (m)
c:Wave speed ()
λ:Wavelength (m)
a:Wave amplitude (m)
H, h:Displacement of channel walls (m)
:Velocity vector ()
:Applied electric field ()
:Velocity components ()
:Time (s)
ϕ:Electroosmotic potential function
:Space coordinates (m)
μ:Dynamic viscosity ()
:Specific heat ()
:Dielectric constant of the medium
:Vacuum permittivity ()
:Density of ionic energy ()
:Number of anions’ density
:Number of cations’ density
:Bulk ions concentration ()
:Boltzmann constant
e:Electric charge (C)
:Average electrolyte solution temperature (K)
:Medium potential (V)
:Ionic energy parameter
:Debye length (m)
m:Electroosmotic parameter
k:Permeability of the medium ()
:Electric field along x-axis ()
:Electrical conductivity ()
:Pressure in wave frame ()
:Fluid density ()
k1:Thermal conductivity ()
:Fluid temperature (K)
:Fluid concentration ()
:Mean temperature (K)
D:Mass diffusion coefficient (m2 / s)
:Ambient temperature (K) and concentration at walls
:Heat transfer coefficient (W/m2K)
:Mass transfer coefficient ()
:Thermal diffusion ratio (m2/s)
:Susceptibility of the concentration (m3/kg)
:Dimensionless space coordinates
:Dimensionless velocity components
:Dimensionless wave number
p:Dimensionless pressure
γ:Joule heating parameter
Pr:Prandtl number
Re:Reynolds number
Du:Dufour number
Br:Brinkman number
Sc:Schmidt number
Sr:Soret number
:Heat transfer Biot number
:Mass transfer Biot number
:Mobility of the medium
:Helmholtz–Smoluchowski velocity
:Stream function
F:Dimensionless time mean flow.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This research was financially supported by the Deanship of Scientific Research, King Faisal University, Al-Ahsa, Saudi Arabia, through Nasher Track under grant no. 206135.