Impact of Electroosmosis and Wall Properties in Modelling Peristaltic Mechanism of a Jeffrey Liquid through a Microchannel with Variable Fluid Properties

The current work emphasizes the modelling of the electroosmosis-modulated peristaltic flow of Jeffery liquid. Such flows emerge in understanding the movement of biological fluids in a microchannel, such as in targeted drug delivery and blood flow through micro arteries. The non-Newtonian fluid flows inside a non-uniform cross-section and an inclined microchannel. The effects of wall properties and variable fluid properties are considered. The long wavelength and small Re number approximations are assumed to simplify the governing equations. Debye-Hückel linearization is also utilized. The nonlinear governing equations are solved by utilizing the perturbation technique. MATLAB is used for the solution, velocity, temperature, skin friction, coefficient heat transport, concentration, shear wood number, and streamlines expressions. The obtained result in optimal electroosmotic velocity (or Helmholtz-Smoluchowski velocity) increases from −1 to 6; the axial circulation has substantial momentum. For larger optimal electroosmotic velocity, a subsequent boost in an axial electric field causes a significant deceleration. Further, the study helps biomedical engineers to create biomicrofluidics devices that may aid in carrying biological fluids.


Introduction
In the human body, peristalsis is the reflex response initiated when lumen contents stretch the gut wall. It mainly occurs in all parts of the gastrointestinal tract, from the esophagus to the rectum. The primary role of the mechanism is to transport the fluids via a distensible tube facilitated by utilizing a progressive wave of expansion and contraction. The movement of food through the throat, chyme transport, urine conveyed through the ureter, ova transportation from the ovary to the uterus, liquid stream in ductus efferent, blood development from tiny blood conduits, and lymph transportation through the lymphatic vessels are some examples of peristalsis. Additionally, it is utilized for planning different contemporary clinical supplies, such as roller-finger siphons and siphoning of blood in the heart-lung machine. The earlier studies on peristaltic flow were carried out by taking Newtonian fluid. However, recent studies have been carried out by taking non-Newtonian fluids because of its broader application in complex structures. The work on Newtonian fluid was initiated by Latham [1] to investigate the peristaltic mechanism in the ureter for urine flow. After that, numerous studies have been conducted recently with taking non-Newtonian liquids in different geometrical assumptions. Inspired by this, the first investigation on peristalsis of a non-Newtonian fluid was carried out by Raju and Devanathan [2] using the power-law model. Manjunatha et al. [3] studied the slip effects on peristaltic transport of Casson fluid in an inclined elastic tube with porous walls. Furthermore, Vaidya et al. [4] investigated the rheological properties and peristalsis of Rabinowitsch fluid through compliant porous walls in an inclined channel. Their study shows that the Biot number increases the temperature for dilatant fluid, whereas it decreases for Newtonian and pseudoplastic fluids. In continuation of [4], Manjunatha et al. [5] studied the effects of wall properties on the peristaltic transport of Rabinowitsch fluid through an inclined non-uniform slippery tube. Their results show that an increase in the value of velocity and thermal slip parameters enhances the velocity and temperature profiles for viscous and shear-thinning fluids. Oudina et al. [6] investigated the entropy and convection effect on magnetized hybrid nano-liquid flow inside a trapezoidal cavity with a zigzagged wall. Their studies claim that magnetic field application strongly influences the flow of nanofluid, and the thermal performance of the cavity is improved with the augmentation of the Reynolds and Hartmann numbers. Furthermore, the statistical approach has been conducted on dual stratification on stagnation point Walter's B nanofluid via radiative Riga plate by Shafiq et al. [7].
In modelling biological liquids, the study of heat transfer characteristics plays an enormous role in understanding fluid behavior. The heat transfer in the body uses four heat exchange systems to maintain homeostasis: conduction, radiation, convection, and evaporation. These mechanisms rely on the property of heat to stream from a greater concentration to a lower concentration. Therefore, heat exchange mechanisms vary in rate according to the temperature level and conditions of the environment. Driven by the application, the heat transfer characteristics of the peristaltic mechanism through porous annulus were examined by Vajravelu et al. [8]. Srinivas and Kothandapani [9] scrutinized the heat transport characteristics of a peristaltic mechanism through an asymmetric channel. Ali et al. [10] scrutinized the impact of heat transport over peristaltic flow through a curved channel. Tanveer et al. [11] examined the heat transport characteristics of peristaltic flow through modified Darcy's law. Further, mass transfer is another essential factor in the scattering interaction, such as the dispersion of supplements from the blood to adjoining cells. In the more significant part of the modern applications, mass transfer assumes a huge role in understanding the cycles engaged with invert assimilation, film division measure, dissemination of synthetic pollutions, and refining measure. Different analysts have considered the impact of heat and mass transport by thinking about these cut off points with various numerical models [12]. Iftikhar et al. [13] studied the peristaltic flow of an Eyring Prandtl fluid in a diverging tube with heat and mass transfer. Their investigations have observed that the rate of heat transfers and mass transfers decreases with an increasing Grashof number. Magesh et al. [14] investigated the heat and mass transfer analysis on non-Newtonian fluid motion driven by peristaltic pumping in an asymmetrically curved channel. Numerical investigation of nonlinear radiation with magneto-micropolar stagnation point flow past a heated stretching sheet was done by Warke et al. [15]. Motivated by the above research, Rajashekhar et al. [16] investigated the effect of mass and heat transport impact on the peristaltic flow of Ree-Eyring liquid with variable properties for hemodynamic flow.
The elastic wall properties strongly affect the wavy fluid flow behavior regarding mean pressure, axial velocity, and net flow rate. These wall properties were articulated as a flexible wall model preserved through elastic spring to disperse the vertical wall displacement and wall damper to damp the average wall velocity. Consideration of wall properties in peristalsis is of unique value in the scrutiny of blood flow in arteries and veins, urine flow in the urethras, and airflow in the lungs. Moreover, investigations on the influence of wall properties and Poiseuille flow in peristalsis were done by Mittra and Prasad [17]. Further, Srinivas et al. [18] studied the effect of heat and mass transfer on MHD peristaltic flow through a porous space with compliant walls. The series solution for the axial velocity, temperature, and concentration distribution has been computed. Hayat et al. [19] investigated the peristaltic transport of viscous fluid in a curved channel with compliant walls. Eldesoky et al. [20] investigated the combined effect of heat transfer, complaint wall properties, and slip conditions on the peristaltic flow through a tube. Recently, Javed et al. [21] studied the peristaltic flow of a realistic fluid in a complaint channel.
The variable fluid properties (viscosity and variable conductivity) play an exceptional role in understanding blood flow in micro arteries. Most of the investigations carried out on peristalsis have considered the viscosity and thermal conductivity of the fluid to be constant. However, recent studies have shown the importance of considering the variable properties in investigating biological liquids. It can be noticed from the literature that the viscosity of the blood varies concerning the thickness and temperature. Nadeem et al. [22] modelled the peristaltic flow by considering the variation in viscosity concerning the thickness of the tube. Further, the impact of slip and variation in viscosity with a thickness on the peristaltic mechanism was investigated by Sinha et al. [23]. Further, Nadeem et al. [24] investigated the effects of temperature-dependent viscosity on the peristaltic flow of a Jeffrey-six constant fluid in a non-uniform vertical tube. The expressions for pressure rise, temperature, pressure gradient, and stream functions are sketched for various embedded parameters and interpreted in their study. However, recently, the variation of viscosity concerning thickness and thermal conductivity concerning temperature variations on the peristaltic flow of non-Newtonian fluids with different geometrical configurations was carried out by [25]. Vaidya et al. [26] studied the peristaltic motion of non-Newtonian fluid with variable liquid properties in a convectively heated non-uniform tube utilizing the Rabinowitsch fluid model. Chu et al. [27] analyzed the entropy analysis in the Rabinowitsch fluid model through the inclined wavy channel by considering both constant and variable properties. Their study shows that a small amount of entropy is observed at the bottom of the channel, and maximum entropy is noted at the ciliated walls under the effect of the Brinkman number. Rajashekhar et al. [28] studied the unsteady flow of Rabinowitsch fluid on peristaltic transport in a non-uniform channel with temperature-dependent properties. Vaidya et al. [29] studied the effect of variable liquid properties on the peristaltic flow of a Rabinowitsch fluid in an inclined convective porous channel with variable thermal conductivity. The study reveals that an increase in the value of porous parameters diminishes the occurrence of the trapping phenomenon for Newtonian and dilatant fluid models. Recently, Rajashekhar et al. [30] studied the peristalsis of Ree-Eyring fluid with homogeneous and heterogeneous reactions addressing hemodynamic applications.
Most of the above studies carried out on peristaltic transport were carried out in the absence of an electric field. However, it has considerable benefits in microfluidic pumping. It permits an effective policy of microchannel circulation areas employing electrical fields because it does not need the standard relocating components featured in traditional micropumps. Electroosmotic pumping has been observed to create a continual pulse-free circulation via applying an electrical field to an electrolyte in contact through an area. The viscous drag triggers the fluid to flow tangentially to the area and creates a consistent net migration of ions. The demand for quick, reputable, and reliable microfluidic systems has generated the strict need for tiny, easy-to-handle, and inexpensive analytic devices. Such applications consist of drug delivery, deoxyribonucleic acid (DNA) sequencing/analysis systems, and biological/chemical agent discovery sensors over integrated circuits. One of the most significant challenges in microfluidic transfer systems in a few of the above applications hinges on increasing the price of fluid circulation to the optimum possible limit through organizing for ideal mixes of flow activating mechanisms. On the other hand, electro-osmosis has ended up being a more straightforward mechanism to develop virtually plug-shaped velocity profiles in microchannels through manipulating the communication of an applied axial electrical field via the local fee density gradients close to the liquid-solid user interface. It has considerable benefits in microfluidic pumping and permits very reliable microchannel flow areas through electric fields. It does not call for the typical relocating elements included in conventional micro-pumps. Electroosmotic pumping has been examined to create a continual pulse-free flow. These pumps are also a lot more open to construction at the microscale. They are significantly being released in biomedical systems, such as medicinal distribution, plasma splitting, electro-osmotically activated bio-microfluidic systems, etc. Flow actuation is realized using an electric field with an electrolyte in contact through a surface. The connection of the surface area through the electrolyte causes a web charge thickness. The viscous drag causes the fluid to move tangentially to the area and creates a constant web movement of ions. Driven by the application of electroosmotic transport on the peristaltic device, Chakraborty [31] designed the peristaltic circulation by taking an electric field. Furthermore, detailed information about electroosmosis and nanopores has been mentioned in a textbook written by Bruus [32]. Ai et al. [33] analyzed and tested the porous membrane with platinum electrodes coated on both sides of a low-voltage nano-porous electroosmotic pump. Wu et al. [34] studied an alternating current electroosmotic pump based on conical nanopore membranes. Asandei et al. [35] investigated the electroosmotic trap against the electrophoretic force near a protein nanopore. Bengtsson and Robinson [36] presented an all-plastic, flexible electroosmotic pump constructed from commercially available materials with the help of an adjustable electroosmotic pump. Huang et al. [37] analyzed the electroosmotic capture and ionic discrimination of peptide and protein biomarkers with FraC nanopores. Chinappi et al. [38] investigated the analytical model for particle capture in nanopores, which elucidates competition among electrophoresis, electroosmosis, and dielectrophoresis. Talapatra and Chakraborty [39] investigated squeeze flow between charged parallel plates by taking an electric field. Later, Bandopadhyay [40] investigated the electroosmosis-modulated peristaltic transport in microfluidic channels. Waheed et al. [41,42] studied the heat and mass transfer in the electroosmosis flow of third-order fluid in the absence and presence of slip effects. Analysis of heat stream in electroosmotic bio-fluid flow in a straight microchannel via peristalsis was done by Noreen et al. [43]. Recently, Tripathi et al. [44] investigated the numerical simulation of double-diffusive convection and electroosmosis during peristaltic transport of a micropolar nanofluid on an asymmetric microchannel.
The mathematical model is developed to examine the electroosmosis modulated peristaltic flow of Jeffrey liquid via a non-uniform channel. The impact of the inclination angle with the horizontal under the effects of wall properties is considered. Further, the variation in viscosity and thermal conductivity are also taken into consideration. The present investigation is the first of its kind in considering the electroosmotic peristalsis of Jeffrey liquid in the presence of changeable liquid properties via a non-uniform inclined channel. The semi-analytical technique (perturbation method) is used to solve the governing nonlinear equations. Outcomes are presented graphically to explore the impact of the parameters emerging in the problem. The current work incorporates flow characteristics and thermal characteristics in the model, which will help stimulate drug delivery in circulatory systems.

Physical Model
Let us consider a 2D electrokinetic modulated flow of Jeffrey liquid through a nonuniform inclined microfluidic channel having half width l(x). The propagation of sinusoidal waveforms thus flow at a continuous speed c along the compliant channel walls (See Figure 1). Let Y = ±h(X, t) be the upper and lower wall, respectively. Here, propagation of waves is along the axial direction (x − axis). The heat and mass transport scrutiny are maintained through giving T 0 , C 0 and T 1 , C 1 to the inferior channel wall and the superior channel wall, respectively. It is assumed that the electric field E is required axially for the liquid circulation. The microchannel wall equation is given by [14]: where l(x) = a 0 + Kx. The governing equation in the laboratory frame of reference for the circulation of Jeffrey liquid in electrohydrodynamics (EHD) are

∂U ∂X
The constitutive equation for Jeffrey liquid is given through [3]: where τ denotes the Cauchy stress tensor, I is the identity tensor, λ 2 denotes the retardation time, S is the extra stress tensor, λ 1 presents the ratio of relaxation to retardation times, p is the pressure, • γ represents the shear rate, and dots on the quantities specify differentiation with respect to time.
The Poisson Equation in the microchannel is [41] where ρ e denotes the total charge density, ξ is the dielectric permittivity, and φ is the electric potential.
The net charge density ρ e follows the Boltzmann variation, given as [42] ρ e = −z v e n − − n + .
The anions n − and cations n + are defined by ρ e of the Boltzmann equation: where n 0 represents bulk concentration, z v the charge balance, T av the average temperature, e the electronic charge, and K B the Boltzmann constant. Using Debye-Hückel linearization approximation, Equation (8) turns into where m e represents the electroosmotic term. The analytical solution of above Equation (11) subject to boundary conditions ∂φ ∂y = 0, at φ = 1 and y = 0, at y = h(x) is attained as [42] φ(y) = cos h(m e y) cos h(m e h) .
The conversion from laboratory frame X, Y to wave frame of references (x, y) is as follows: The above conversions are utilized in Equations (2)-(6) and then introduce the dimensionless variables: Here, λ D is the Debye length, Sc the Schmidt number, y is the dimensionless transverse coordinate, U hs is the Helmholtz-Smoluchowski velocity, m e is the electroosmotic term, p is dimensionless pressure, Ω is the concentration, α 3 is the peristaltic wave number, θ is the temperature, P r is the Prandtl number, ψ is a dimensionless stream function, S is the dimensionless shear stress, R e is the Reynolds number, Br is the Brinkman number, is the amplitude ratio, Sr is the dimensionless Soret number, x is the dimensionless axial coordinate, and E c is the Eckert number.
The governing equation of movement of the flexible wall may be articulated as where R is an operator that is utilized to represent the movement of stretched membrane through viscous damping forces, such that Continuity of stress at y = h and utilizing x-momentum equation yield where τ represents the elastic tension in the membrane, n 1 denotes the mass per unit area, n 2 is the coefficient of viscous damping forces, h represents the dimensional slip term, and p 0 denotes the pressure over the outside area of the wall due to the tension in the muscles. We supposed p 0 = 0. Equation (2) is identically satisfied and, from Equations (3)-(6), subject to longwavelength assumption and ignoring superior powers of α 3 , we acquire ∂p ∂x The non-dimensional boundary conditions are The expression for variable viscosity is The thermal conductivity term is where α 1 is the coefficient of viscosity and α 2 is the coefficient of thermal conductivity.

Methodology
The transformed governing equation and boundary conditions need to be solved using suitable analytical or semi-analytical methods to obtain the solution for velocity, temperature concentration, and streamlines. The semi-analytical technique, the perturbation method, is employed to solve the system of non-linear transformed equations.

Perturbation Solution
Equations (18)-(21) specify a governing differential system with boundary conditions (22) and (23) that are extremely nonlinear. It is impossible to find an exact solution to this system of equations. As a result, the system can be solved analytically using the perturbation method, which uses variable viscosity α 1 and variable thermal conductivity α 2 as perturbation terms. Thus, ψ and θ can be expanded as Substituting Equations (26) and (27) in Equations (18) and (20) and collecting the same power terms and equating to zero yields the following systems.

Zeroth Order System
Subject to the following boundary conditions

First Order System
Subject to the following boundary conditions The solution of zeroth and first order systems along with boundary conditions are obtained using the MATLAB code. The separate code has been written to acquire the solutions for streamline, velocity, and temperature. Further, the solution for Equation (21) is acquired by twice integrating the temperature term through the aid of boundary conditions (22) and (23).

Results and Discussion
The present investigations emphasize the impact of electroosmotic velocity (U hs ), the electroosmotic parameter (m e ), the Jeffrey parameter (λ 1 ), variable viscosity (α 1 ), variable thermal conductivity (α 2 ), angle of inclination (α), the non-uniform parameter (m 1 ), the wall tension parameter (E 1 ), the mass characterization parameter (E 2 ), the wall damping term (E 3 ), the wall rigidity term (E 4 ), the elastic term of the wall (E 5 ), the Brinkmann number (Br), the Schmidt number (Sc), and the Soret number (Sr) on velocity (u), temperature (θ), concentration (φ), and streamlines (ψ). The outcomes of these parameters are shown through Figures 2-13. Figure 2 depicts the effect of relevant terms on velocity. The velocity profiles have a parabolic shape due to an optimum velocity in the conduit's center. Figure 2a shows that as U hs , or the optimal electroosmotic velocity (or Helmholtz-Smoluchowski velocity), increases from −1 to 6, the axial circulation has substantial momentum. For larger U hs , a subsequent boost in an axial electric field causes a significant deceleration. The circulation is supported and obstructed by an essential electrical field strength. The more critical values of m e in Figure 2b enhance the velocity profiles. Furthermore, there is a slight variation in velocity profiles for higher values. The mutual of the Debye length is the Debye-Hückel specification. Because increased ion migration occurs as we move away from the indicted area, the decreasing Debye length, i.e., the rising Debye-Hückel term, is thought to raise electrical capacity. As a result, the Debye length is a vital style criterion in managing potential electrical variation, which significantly impacts the axial velocity field. Figure 2c has been drawn for multiple values of the Jeffrey parameter. The velocity profile improves as the Jeffrey criterion rises, as shown in this graph. It is worth noting that by taking λ 1 = 0, the outcome decreases to Newtonian liquid. Figure 2d depicts the effect of changeable viscosity over velocity profiles. The velocity profile enhances for more significant values of variable viscosity because a boost in variable viscosity parameters decreases the fluid's viscosity and thus improves the fluid circulation. Figure 2e reveals an augment in velocity profiles for more significant values of inclination parameters. The increase in the inclination angle gives rise to a gravitational force, and thus the velocity of the fluid increases for the higher value of the inclination angle. The effect of the nonuniform term over the velocity profile is depicted in Figure 2f. Here, a rise in the value of the non-uniform term enhances the flow region and thereby increases the velocity profile. The impact of the wall tension term, the mass characterizing parameter, the wall damping term, the wall rigidity parameter, and the elastic parameter of the wall on the velocity profile is sketched in Figure 2g. The velocity can be seen to be decreasing as we move from the center of the channel towards the walls for each of these elastic parameters. Moreover, As E 1 and E 2 increase, the velocity rises, whereas with higher values of E 3 , E 4 , and E 5 , the fluid decelerates. Clinically, when blood walls are flexible or the mass per unit surface rises (as in blood capillaries), the exchange of oxygen, water, and other nutrients becomes simple. The damping tends to react oppositely.

Temperature Profiles
In this section, we examine the impact of various specifications on temperature profiles. The temperature distributions are parabolic through a maximum temperature level in the middle of the conduit. The ramifications of viscous dissipation, which raises the temperature in the center of the channel, may affect these habits. The viscosity of the liquid causes this phenomenon, which changes kinetic energy in the liquid into thermal energy inside it. Figure 3a reveals that the temperature of the liquid rises through a boost in the maximum electroosmotic velocity. On the other hand, the temperature profile decreases for the increasing Debye-Hückel parameter (See Figure 3b). Figure 3c is drawn to see the impact of the Jeffrey specification on the temperature profile. The figure reveals that temperature is a rising function of the Jeffrey specification. Figure 3d is sketched to see the influence of variable thermal conductivity over the temperature profile. Thermal conductivity is the fluid's ability to protect or release heat to its settings. As a result, when the thermal conductivity of the liquid within the channel is higher than the wall temperature, the fluid's thermal evaluation increases. The impact of variable viscosity over temperature profile is illustrated in Figure 3e. Variable viscosity raises the temperature profile, and when variable viscosity amounts to absolutely zero, consistent viscosity can be derived. The reality is that, as the fluid's variable viscosity increases, the heat transfer prospective declines and raises the temperature profile. Figure 4a illustrates that the higher values of the inclination angle boost the temperature profiles. Figure 4b shows that higher values of the Brinkmann number improve the temperature distribution. It is an argument that less heat conducts the viscous dissipation. The Brinkman number shows dissipation intensity. Given that the wall temperature is convectively regulated, heat transfers are considered because of the impacts of heat dissipation. A similar trend is noticed to that of Brinkmann number for an increase in the non-uniform parameter. (See Figure 4c). The effects of the elastic parameters E 1 , E 2 , E 3 , E 4 , and E 5 are sketched in Figure 4d. From the figure, it has been determined that the temperature profiles increase for larger values E 1 and E 2 in the non-uniform channel, whereas the opposite behaviour is noticed for E 3 , E 4 , and E 5     Figure 5a-e shows the effect of electroosmotic velocity, the electroosmotic parameter, the Jeffrey parameter, changeable thermal conductivity, and variable viscosity on the concentration profiles. The concentration profile behaves opposite to temperature profiles in Figure 5a-e. This behavior is predicted physically because mass and heat are comprehended to be opposite in direct proportion. In addition, the profiles reveal that particle matter in the liquid is more at the periphery than in the central location of the conduit. From a biological point of view, this behavior assists in spreading essential nutrients to nearby cells and tissues in the blood and other fluxes. The effect of the inclination angle on the concentration is drawn in Figure 5f. It can be noticed from the figure that enhancing the angle of inclination boosts the concentration profiles. Figure 6a,b is sketched to see the impact of the Schmidt and Soret numbers on the concentration profiles. From the figures, the concentration profiles notably decrease through an augmented value of the Schmidt and Soret numbers. The impact of the Brinkman number and the non-uniform term on the concentration profiles are illustrated in Figure 6c,d. It can be noticed that the concentration profiles decrease through a rise in the value of the Brinkman number and the non-uniform parameter. Figure 6e represents the impact of the elastic parameters E 1 , E 2 , E 3 , E 4 , and E 5 on concentration profiles. It has been noted that the rise in the values of E 1 and E 2 enhances the concentration, whereas the opposite behavior is observed for an increase in the values of E 3 , E 4 , and E 5 .

Trapping Phenomenon
Trapping is the most interesting and common phenomenon in the peristaltic transport of bio fluids. The behavior of the relevant parameters over the bolus is analyzed with the help of contours. Figure 7a-d is drawn for the influence of electroosmotic velocity on the trapped bolus. From the figure, the formation of boluses increases drastically through the higher values of electroosmotic velocity. This is because powerful electroosmotic velocity means a dominant external electrical field; the variety of enclosed boluses considerably diminishes for large values of electroosmotic velocity. Figure 8 explains that the union of enclosed boluses rises through an enhancement in the electroosmotic term. It can be noticed that the larger the value of electroosmotic parameter, the stronger the enclosed bolus appears in EDL. Figure 9 is drawn for the effect of the Jeffery parameter on the trapped bolus. The figure reveals that the volume of the trapped bolus improves with a larger value of the Jeffery parameter. A similar trend has been observed for the rise in the significance of variable viscosity, angle of inclination, and the non-uniform parameter. (See Figures 10-12). Figure 13 illustrates the influence of the elastic parameter E 1 , E 2 , E 3 , E 4 and E 5 on the trapped bolus. The volume of the trapped bolus decreases for larger values of E 5 , and it increases for E 1 , E 2 , E 3 , and E 4 .

Conclusions
This investigation describes the impact of electroosmotic properties, wall properties, and variable liquid properties on the peristaltic flow of the Jeffrey liquid in an inclined channel. The perturbation technique is utilized to obtain the solutions for velocity and temperature expressions. MATLAB programming is used to discusses the impact of critical terms over velocity, concentration, temperature, and streamlines. The outcomes of this paper have potential applications in both the industry and medical areas of research. Besides considering the slip and variable fluid properties in peristaltic circulation, physical conditions, such as an asymmetric/tapered uneven artery through a porous medium, can be applied in future research. The present investigation provides a fair theoretical approximation for future research in this direction. The significant primary outcomes of the current model are summed as follows: • The impact of optimum electroosmotic velocity enhances the velocity profiles. • Electroosmotic specification is a decreasing function of velocity.

•
The increasing values of the Brinkman number and the non-uniform parameter enhance the velocity profiles.

•
The Jeffrey term boosts the temperature and velocity profiles.

•
The changeable viscosity and variable thermal conductivity boost the velocity and temperature profiles.

•
The concentration profiles significantly diminish for more significant values of the Schmidt and Soret numbers.

•
The influence of E 1 and E 2 enhances the velocity and temperature profiles for Newtonian fluids, and it decreases for E 3 , E 4 , and E 5 .

•
The size of the trapped bolus increases through an expansion in the estimation of changeable viscosity. • Accumulation of streamlines decreases with an increase in maximum electroosmotic velocity.

•
The number of enclosed boluses is augmented through a boost in the electroosmotic criterion. Funding: This research received no external funding.

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