Impact of activation energy and variable properties on peristaltic flow through porous wall channel

Aug 21, 2023

Scientific Reports volume 13, Article number: 3219 (2023) Cite this article

661 Accesses

1 Citations

Metrics details

The current study discusses the peristaltic flow of Jeffrey fluid through a porous wall channel. Magnetohydrodynamic (MHD) effects are also considered while formulating the problem. Heat and mass transfers are discussed in the presence of activation energy and constant heat source/sink effects. A chemical reaction is also part of the analysis. The Lubrication approach is adopted for the simplification of resulting non-linear equations. MATHEMATICA command, NDSolve, is used to discuss the results graphically for various flow parameters like Hartman number \((M)\), porosity parameter \((k)\), slip parameters (\(\gamma ,{\gamma }_{1},{\gamma }_{2}\)), Schmidt \((Sc)\), Soret \((Sr)\) and Prandtl \((Pr)\) numbers, and many others. Parabolic behavior for velocity and sinusoidal nature for heat transfer and pressure gradient is noticed. Results indicate that the velocity is greatly affected by varying values of slip parameters (γ′s) and Hartman number \((H)\). Enhancing the viscoelastic nature of fluid causes an increase in velocity. Similar behavior is noticed for velocity and temperature profiles. The decreasing trend is shown by concentration when the value of the chemical reaction and temperature ratio parameters is enhanced. Thus, the study presented in the current analysis can be used to study many human physiological systems especially, the blood flow. Since Jeffrey's fluid exhibits the same characteristics as observed for blood.

Mathematical modeling is employed in biomechanics to investigate the physiological systems. Bio-fluid mechanics is a section of biomechanics which reveals the kinematics and dynamics of body fluids in living beings. Advancement in bio-fluid mechanics enable scientists to study the blood vessels liquid stream, the respiratory tract, the lymphatic system, the gastrointestinal tract, the urinary tract and various other. Recent investigations disclose the clinical applications like artificial organs, vascular vessel advancement, designing medical instruments, creation of material membranes for orthopedics, and many more. Analogous processes of bio-liquids flow can be seen in variety of situations within human body among which the prominent is peristalsis and can be seen as basis for the current study. The main purpose of peristalsis is to move the fluids through the tubular structure without requiring an overall pressure difference. The term peristalsis originates from the Greek word peristaltilkos, meaning "compressing and clasping". According to Merriam-Webster's1 peristalsis is consecutive waves of involuntary contraction passing along the walls of a hollow muscular structure and pushing the contents forward. The human body's peristalsis mechanism started working after food was chewed, swallowed as a bolus, and passed through the oesophagus. To prevent the bolus from voyaging back toward the mouth, smooth muscles behind it constrict. It was first described as a kind of motility where there is contraction above and relaxation below a transported by Bayliss and Starling2. Industrial application of peristaltic pumping is exploited in different applications which include sterile fluids exchange, blood pump in heart lung machines, transportation of internecine and dangerous liquids to prevent their involvement in surrounding environment etc. A noteworthy modern use of peristaltic pumping can be seen in designing the roller pumps which are used to avoid the contact of fluid with the pumping equipment. Peristaltic transport on viscous liquids was first introduced by Latham3, in 1966. This study was further extended by Shapiro et al.4.

In fact, not every fluid possesses the characteristics of a Newtonian fluid. Therefore, we include non-Newtonian fluids in our discussion. However, realistically speaking, complicated fluids like food boluses travelling through the oesophagus, urine passing through the ureter, or chyme traversing the gastrointestinal track would not adhere to Newton's principles of viscosity. As a result, no single constitutive relation can predict the characteristics of all fluids. In response to this issue, a number of constitutive models have been proposed to identify the properties of non-Newtonian fluids. The equations for the flow of pigment oil suspensions used in printing ink-type fluid were expanded by Casson5. The theory of the micro-polar fluid was covered by Eringen6, who also thoroughly investigated features such couple stresses, body couples, micro-rotation, and micro inertial effects. The theory of micro-fluids, of which micro-polar is a specific instance, was first presented by Eringen6. Among all of these models, the Jeffrey fluid, which possesses both relaxing and retarding qualities, is comparatively one of the simplest types of viscoelastic fluid. Because Jeffrey fluid can forecast relaxation/retardation time effects, which are crucial for the analysis of viscoelastic properties in the polymer industries and human physiological systems. Ramanamurthy et al.7 investigated the peristaltic flow of viscous fluid through two-dimensional curved channel. Their main aim is to analyze unsteady nature of the flow. Nadeem et al.8 performed the endoscope analysis of Prandtl liquid in both fixed and wave frame of references. They have addressed the effect of various wave shapes on endoscope. Power-law fluid flow through cylindrical tube is examined by Sadeghi and Talab9. The results indicate enhancement in the fluid flow due to high values of power-law index. Tripathi et al.10 developed mathematical model to discuss intestinal flow by taking bi-viscosity non-Newtonian fluid. Their study helps in better understanding of gastric-hydrodynamics. Axisymmetric flow of Bingham fluid through cylindrical geometry is investigated by Fusi and Farina11. Ramesh and Devaker12 modeled endoscope problem to discuss its application to biomedicine. They have taken couple stress fluid to model physiological fluid. Application of peristalsis to chyme movement in gastrointestinal tract can be seen through Vaidya et al.13. Their study discloses the increasing impact of variable viscosity on bolus size. In the culmination of above mentioned literature, we can observe the real life applications of non-Newtonian fluid flow in medical and industry.

The study of magneto-hydro-dynamics is concerned with the motion of highly conducting fluids in the presence of a magnetic field. The velocity of the conducting fluid through magnetic field generates electric current that change the magnetic field and generates mechanical forces that change the fluid flow14. Due to its substantial applications, such as in materials processing, magneto-hydro-dynamic (MHD) energy generators15, cancer therapy16, and biomedical flow control and separation devices17, the impact of bio-magnetic fluid dynamics has attained a prominent position. Its use in biomedical engineering includes magneto-fluid rotary blood pumps, MHD medication targeting, and controlling hyperthermia in the cardiovascular system18,19. A device that uses a magnetic field and a highly sensitive sensor to detect minute movements of an object within the magnetic field is known as a giant magneto resistive (GMR) device. Research into peristaltic activity inside tubular organs like the colon, fallopian tubes, and even the vas deferens has improved because to this technique. Satyanarayana et al.20 explores the impact of magnetic field on peristaltic flow of micropolar fluid through asymmetric channel. Their study reveals that velocity increases with the enhancement in micro-rotaion parameter. Peristaltic flow with variable thermal conductivity and viscosity is addressed by Latif et al.21. Their study concluded that Newtonian fluid has low heat transfer rate as compared to third order fluid. Prakash et al.22 considered Williamson fluid to model blood flow in the presence of magnetic field. Selvi and Sirinivas23 addresses the peristaltic transport of Herschel-Bulkley fluid through non-uniform tube. They verify their results by comparing them to Vajravelu et al. Shera et al.24 worked on mathematical analysis of thermal therapy used for the treatment of cancer.

Heat transfer plays vital role in industrial and medical applications. Especially, within human body heat transfer is an essential area of research. Bio-heat transfer in tissues has attracted the attention of biomedical engineers for thermotherapy25 and the human thermoregulation system26. The heat transfer inside human body occurs as conduction in tissues, perfusion of the arterial-venous blood through the pores of the tissue, metabolic heat generation whereas the annihilation of cancerous cells, dilution technique of blood flow and vasodilation. In connection with peristalsis heat transfer become significant in oxygenation and hemodialysis. Several researchers investigated about heat transfer in peristaltically induced flows. In general, viscous dissipation effect is ignored while performing theoretical analysis of fluid flow problems. But under certain situations this supposition may lead to doubtful results. The need of considering viscous dissipation effects is felt while dealing with strong temperature-dependent viscosity, high viscosity fluids and high-speed gas dynamics. Heat produced due to viscous dissipation may increase the temperature of tube wall consequently decreasing viscosity which results in increased velocity and temperature. Thus, instead of viscous dissipation effects variable thermal conductivity is considered to make analysis more realistic. While discussing heat transfer, we cannot ignore mass transfer phenomenon as simultaneous occurrence of both can be seen through many applications such as drying, energy transfer in wet cooling tower, evaporation at the surface of a water body and the ow in dessert cooler. Additionally, mass transfer manifests as a result of the different species' concentrations in a combination liquid. Such characteristics change as a mixture is transported from areas of higher to lower concentration. Additionally, activation energy, which is referred to as the least obligatory energy that chemical reactants must gain before the chemical reaction takes place, is one of the most crucial characteristics for chemical reactants. Chemical engineering, geothermal reservoirs, oil emulsions, and water mechanics, among other fields, all heavily depend on the consideration of mass transfer with both chemical reaction and activation energy. Tanveer et al. presented the analysis for electroosmotic peristaltic flow of nanofluid with non-Newtonian base fluid (see refs.27,28). The study focuses the attention on thermal aspects of the flow. Thus, finding promising applications in micro-fabrication and chemical industry. Maryam et al.29 presented the meta-analysis of homogeneous–heterogeneous chemical reaction on peristaltic flow through wavy curved geometry. Ahmed et al.30 studies thermal radiation effects on peristaltic flow of nanofluid with mixed convection. The study unveils that the magnetic field tends to increase the thermal energy of the flow. Microfluidic peristaltic flow is examined by Noreen et al.31 in view of heat transfer and electroosmotic effects. Convective heat and mass transfer for Prandtl nanofluid through non-uniform channel is analyzed by Akram et al.32. Khazayinejad et al.33 presented mathematical model for the treatment of cancer by injecting graphene particles in to blood stream and then, applying external magnetic field. Nanoparticle volume fraction enhancement results in greater heat transfer. Hence, shows great impact on the destruction of cancerous cells. Imran et al.34,35,36 discusses effect of various types of chemical reactions on peristaltic flow through various geometries. Analysis of physiological liquid flow between absorbent barriers, such as blood, became crucial, particularly in the lungs. According to Fung and Tang37 and Gopalan38, the lung can be seen as a duct encircled by two thin, porous media layers. Therefore, other researchers like Naveed et al.39,40 presented mathematical formulation of different non-Newtonian fluids through porous medium as real life application of peristalsis. Some other important and meaningful research on miscellaneous disciplines i.e., nanoribbons, non-uniform wall thickness, granular thermodynamics framework41,42,43, material analysis44,45,46, fluid flow and mass transport in multi-scale heterogeneous media and porosity effects47,48,49,50. Saima et al.51 focus their attention on the heat and mass transfer analysis of micro-polar nanofluid through a lid driven cavity. They adopted finite element method to solve the obtained nonlinear system of equations. The results indicate that great mass diffusion occurs inside the cavity for small Schmidt number. Rasool et al.52 presented the study for thermal analysis of Maxwell nanofluid over isothermally heated surface. Their study is based on the comparison on convective and non-convective boundary conditions. In addition, entropy generation analysis of multi-walled carbon nanoparticles (MWCN) inside vertical Cleaveland Z-staggered cavity is also performed by Rassol et al.53. The outcome of the study shows that for higher Reynolds number Bejan number decreases creating great impact on entropy generation. Investigation of elrctro-magneto-hydrodynamic nanofluid flow over riga plate is investigated in54. The study unveils that the heat transfer across the wall can be controlled by adjusting convective conditions.

The flow field can be dramatically altered by the suction or injection of fluid through the bounding surfaces, as in mass transfer cooling, which then affects the rate of heat transfer from the bounding surfaces. In general, injection functions in the opposite ways of suction, which tends to enhance skin-friction and heat transfer coefficients. Fluid injection or withdrawal via porous heated or cooled surfaces is generally of importance in real-world issues like blood flow through arteries, fluid flow through urinary tract etc. This may result in the system being heated (or cooled) more effectively. Moreover, viscoelastic nature of blood can be exactly exhibited by the considered Jeffrey fluid model. Also, different types of chemical reactions within human physiological systems, particularly during blood flow, greatly affect the fluid flows. Therefore, giving consideration to such effects provide us with more realistic mathematical model of various physiological systems. Therefore, goal of this study is to develop a mathematical model for Jeffrey fluid that incorporates activation energy and peristaltic flow via a porous walled channel. In order to adjust the MHD flow brought on by peristalsis, the presence of porous walled channels and the effects of activation energy on non-Newtonian fluid will be crucial. The construction of nonlinear coupled equations via mathematical modelling will be discussed, and the lubrication approach will then be used to simplify them. The plotting of graphs will be used to explain the effects of applicable factors and parameters on the flow via built-in NDSolve command of MATHEMATICA.

Through a flexible walled conduit, we examine the incompressible 2D peristaltic flow of Jeffrey fluid. The walls of the channels are permeable. A magnetic field is taken that is parallel to the flow. Sinusoidal waves with a wavelength of λ are moving at a speed of \(c\) along the walls. Let \(u\) and \(v\) represent the axial and transverse velocity components, respectively. The Cartesian coordinate system is used to discuss the flow (see Fig. 1).

Geometry of the problem.

Wall surface geometry mathematical expression is expressed as under3,4:

where \(\overline{t }\) denotes time, \(\lambda\) wavelength, \(b\) wave amplitude, \(a\) width of channel, \(c\) is the wave speed and \(\overline{X }\) is the horizontal direction. The governing equations for the flow are55,57,58,59:

where shear stress components are given as60:

where variable thermal conductivity is taken as56:

The boundary conditions for the analysis are considered as26,27:

In above equations, \(\overline{U }\) and \(\overline{V }\) stands for velocity components in \(\overline{X }\) and \(\overline{Y }\) directions respectively whereas ρ denotes density. \(\overline{{V }_{0}}\) is the velocity through porous walls, \(\overline{P }\) is the pressure and shows electrical conductivity. \(\overline{{B }_{0}}\) is magnetic field strength and specific heat capacity is shown by \({C}_{p}\). \(D\) is mass diffusion coefficient and \({Q}_{0}\) is heat source/sink parameter. \({K}_{1}\) is chemical reaction parameter and \({K}_{2}\) is the reaction rate. \({E}_{a}\) stands for activation energy, \(n\) is the fitted rate constant and \({K}^{*}\) is Boltzmann constant, \({K}_{m}\) be the thermal conductivity at constant temperature and α be the constant. \(\overline{\gamma }\), \(\overline{{\gamma }_{1}}\) and \(\overline{{\gamma }_{2}}\) are velocity, temperature and concentration slip parameters respectively.

The transformation between steady and unsteady reference frames is

where \(\overline{{V_{0} }}\). is the suction/injection parameter at the walls. The dimensionless constants/variables are defined as:

After using the above transformation, dimensional equations take the form

and shear stress components are given as:

where \(\delta\) is the wave number, \(Re\) the Reynolds number, \(H\) stands for Hartman number, \(Pr\) is the Prandtl number whereas \(B\) shows constant heat source/sink parameter. \(Sc\) exhibits Schmidt number, \({K}_{c}\) is chemical reaction parameter, \(\Omega\) is temperature ratio parameter, \(E\) the activation energy parameter whereas \(n\) is the fitted rate constant. Employing long wavelength approximation to above equations, we get:

and the shear stress components become:

whereas, variable thermal conductivity takes the form:

Wall surface geometry expressed as:

The non-dimensional boundary conditions are:

where \(\gamma , {\gamma }_{1}, {\gamma }_{2}\) are velocity, temperature and concentration slip parameters respectively. The heat transfer at the walls is expressed as:

In the fixed frame, the instantaneous ow rate is given by58:

In a wave frame Eq. (42) takes the form58:

Over one period of the peristaltic wave, the average non-dimensional volume flow rate \(\overline{Q }\) is defined as58:

Since the equations established in the previous section are coupled and nonlinear, they cannot be precisely solved. As a result, exact solutions to these equations cannot be found. But as technology has advanced, a variety of built-in software programs have emerged that can offer the best numerical approximation to such system of equations especially, in confined domains where a maximum to minimum range is attainable. In this method, the entire system is added to obtain a direct graphical illustration of the problem under consideration. Such a method has the advantage of offering a level of accuracy with the least amount of CPU time (5–25 min) every evaluation. One such example is the problem under consideration, which may be solved using the built-in command NDSolve in the computer program MATHEMATICA which provide graphical output with the necessary boundary conditions. Since, the technique is based on shooting method which shows good efficiency for boundary value problem. However, the technique automatically choose starting initial condition which may create problem. As having good starting point become crucial for finding the best solution. Since, this problem can be fixed after choosing appropriate initial point for starting the calculations. Therefore, various researchers have adopted the technique to solve and justify the obtained solutions (see refs.28,30).

The purpose of the current part is to thoroughly examine how embedded parameters affect flow quantities. For some parameter values, the study that was conducted is useful. Results are discussed for fixed values of \(x = 0.3\) and \(\varepsilon = 0.1\). Whereas, the other parameters range is taken as, \(0\le k\le 0.5\), \(0\le H\le 2.0\), \(0\le {\lambda }_{1}\le 1.5\), \(0\le \alpha \le 0.5\), \(0\le Pr\le 2.0\), \(0\le B\le 1.0\), \(0\le Sc\le 0.5\), \(0\le \xi \le 3.0\), \(0\le \Omega \le 2.0\), \(0\le {\gamma ,\gamma }_{1} ,{\gamma }_{2} \le 0.3\) and \(0\le E\le 2.0\).

This subsection is prepared to analyze the impact of various parameters on temperature \(\theta\). Figure 2 discloses the effect of variable thermal conductivity parameter \(\alpha\) on \(\theta\). It can be seen that temperature decreases for higher values of \(\alpha\). This is due to the fact that material ability to absorb or disperse heat is enhanced for greater values of \(\alpha\). To study impact of \(B\) on \(\theta\), Fig. 3 is plotted. The figure shows the enhanced temperature for increasing the values of \(B\). Such enhancement in temperature is caused by the heat generation resulted from the friction between fluid layers. Figure 4 depicts the behavior of thermal slip parameter \({\gamma }_{1}\) on \(\theta\). Temperature rises for greater values of \({\gamma }_{1}\). The reason behind this increase can be linked with velocity. As temperature is the average kinetic energy of particles. Increasing slip will result in increased velocity hence, accelerating fluid flow. Therefore, temperature rises. Similar behavior is shown by suction/injection parameter \(k\) (see Fig. 5). Increase in the values of \(k\) results in large pores. In turn, more fluid will flow through them resulting in higher velocity. Since, temperature is the average kinetic energy of particles. Therefore, temperature enhances. Figure 6 shows the impact of \(Pr\) on \(\theta\). Rise in temperature is noticed for higher values of \(Pr\). Given that Prandtl is the ratio of momentum to thermal diffusivity. Enhancing Prandtl number causes increase in internal friction among fluid layers. Therefore, a higher \(Pr\) indicates a higher temperature.

Temperature presentation of α.

Temperature presentation of B.

Temperature presentation of \({\gamma }_{1}\).

Temperature presentation of \(k\).

Temperature presentation of \(Pr\).

This subsection analysis the behavior of involved parameters on velocity \(u\). Figure 7 shows the effect of slip parameter on \(u\). Velocity declines for higher values of \(\gamma\). This is due to the fact that as the slippage occurs on the boundary, velocity become different from the boundary. Which, in turn, decreases the effect of boundary on the motion of fluid, resulting in decreased velocity. Figure 8 portrays the significance of Hartman number \(H\) on \(u\). Since, magnetic field is the resistive kind of force. Therefore, it decreases velocity. Figure 9 discusses the impact of \(k\) on \(u\). Greater the value of suction/injection parameter, larger the number of pores, thus more fluid flow results. Fluid parameter \({\lambda }_{1}\) effects the velocity in increasing manner. Since it is viscoelastic parameter. Therefore increase in \({\lambda }_{1}\) always results in increased elastane causing fluid velocity to increase. (see Fig. 10).

Velocity presentation of \(\gamma\).

Velocity presentation of \(H\).

Velocity presentation of \(k\).

Velocity presentation of \({\lambda }_{1}\).

Through the energy equation, the heat transfer phenomena has a significant impact on the peristalsis process. Therefore, heat transmission considerations are necessary for physiological and industrial applications. Thus, Figs. 11, 12, 13, 14 and 15 study the effect of various parameters on heat transfer profile. It is noticed from the figures that heat transfer shows sinusoidal behavior. This is due to fact that waves at the boundary are of sinusoidal nature due to the consideration of peristaltic waves. Figures 11 and 12 depict the behavior of variable thermal conductivity parameter α and thermal slip parameter \({\gamma }_{1}\). It is noticed from the figures that magnitude of heat transfer decreases for enhancing values of α and \({\gamma }_{1}\). Since increasing the value of α reduces the ability of fluid to absorb heat. Thus, resulting in decrease heat transfer rate. Also, the presence of pores at the boundaries greatly effects the fluid velocity which is directly related with temperature. Therefore, decrease in temperature results in less heat transfer. Figure 13 is sketched to study the effect of suction/injection parameter on \(Z\). It is noticed from the figure that rate of heat transfer decreases as we increase the value of \(k\). Because higher number of pores results in increased velocity causing temperature to rise. Thus, resulting in higher rate of heat transfer. As we consider \(B\) less than \(0\), excessive heat is absorbed resulting in decreased rate of heat transfer. On the other hand, taking \(B\) greater than \(0\) generates more heat causing higher rate of heat transfer (see Fig. 14). Impact of Prandtl number \(Pr\) on heat transfer rate can be depicted through Fig. 15. Since, temperature shows increasing trend for \(Pr\), thus have higher rate of heat transfer.

Heat transfer presentation of \(\alpha\).

Heat transfer presentation \({\gamma }_{1}\).

Heat transfer presentation \(k\).

Heat transfer presentation \(B\).

Heat transfer presentation \(Pr\).

The axial pressure gradient expressed in terms of independent variable \(x\) is plotted for various corresponding parameters. Fluctuating behavior is shown by pressure gradient attaining its minimum at (0.0…, 0.95…), whereas approaching maximum at 0.5 …. This shows the existence of high level flow through the tube without the need of greater pressure gradient. The reason behind fluctuating behavior of pressure gradient is the peristaltic wave traveling along the boundary. For dealing with the results for pressure gradient and pressure rise Figs. 16, 17, 18, 19, 20, 21, 22 and 23 are plotted. For velocity slip, \(dp/dx\) enhances, whereas, similar behavior is noticed for Jeffrey fluid parameter \({\lambda }_{1}\) (see Figs. 16 and 18). The higher values of Hartman number \(H\) results in greater pressure gradient (see Fig. 17). Figure 19 studies the impact of suction/injection parameter on \(dp/dx\). It is observed that pressure gradient rises for \(k\).

dp/dx for \(\gamma\).

dp/dx for \(H\).

dp/dx for \({\lambda }_{1}\).

dp/dx for \(k\).

Δp for \(H\).

Δp for \(\gamma\).

Δp for \(k\).

Δp for \({\lambda }_{1}\).

Additionally, pressure rise for one wavelength is used to examine the maximum pressure rise at which peristalsis acts as a pump. Non-Newtonian fluid characteristics can easily be described by nonlinear nature of these curves. All these figures are composed of four primary sections: the peristaltic pumping zone \((\Delta P > 0)\), the free pumping region \((\Delta P=0)\), and the co-pumping region \((\Delta P< 0)\). Peristalsis, which occurred as a result of pressure difference, causes flow rate to be positive in the zone of peristaltic pumping, whereas peristalsis of the channel boundaries produces a free-pumping region. Negative pressure difference helps the peristalsis-related flow in the co-pumping zone. The effect of \(H\) on \(\Delta p\) is studied through Fig. 20. It is noticed from the figure that making increment in \(H\) results in decrease of \(\Delta p\) in co-pumping region. Figure 21 depicts the effect of on pressure rise. Increase in \(\Delta p\) can be seen in co-pumping region whereas shows decrease in \(\Delta p\) in free pumping region. Figure 22 depicts the impact of suction/injection parameter \(k\) on \(\Delta p\). It is noticed that pumping curves meet at point \(Q\approx 0.1\). For \(Q>0.1\), pressure rise enhances whereas opposite behavior is noticed for \(Q<0.1\). Figure 23 shows the effect \({\lambda }_{1}\) of on \(\Delta p\). Pumping rate is found to be enhancing in co-pumping region whereas shows opposite behavior in free pumping region.

This subsection addresses the effect of various parameters on concentration profile \(\phi\). The concentration profiles exhibit the opposite behavior from the temperature profiles. Physically speaking, this makes sense because heat and mass are known to be opposites. Additionally, the patterns show that the fluid particles are more concentrated close to the channel walls. For this purpose Figs. 24, 25, 26, 27, 28 and 29 are plotted. It can be analyzed through Fig. 24 that concentration increases for higher values of variable thermal conductivity parameter \(\alpha\). Since variation in thermal conductivity of material greatly influence the temperature, thus effecting concentration. Figure 25 depicts the impact of \(\xi\) (rate of chemical reaction) on ϕ. The concentration is decreasing function of \(\xi\). To elucidate the effect of concentration slip parameter \({\gamma }_{2}\), Fig. 26 is plotted. Concentration increases for higher values of \({\gamma }_{2}\). To demonstrate the effect of Schmidt number \(Sc\) on \(\phi\), Fig. 27 is sketched. Enhancement in the values of \(Sc\) decreases \(\phi\). Increase in \(Sc\) causes concentration diffusion to decrease, hence resulting in less concentration at the point. Impact of activation energy \(E\) can be depicted through Fig. 28. It shows that concentration enhances for higher values of \(E\). Figure 29 elucidates the effect of temperature ratio parameter . Decrease is noticed for concentration when values of Ω are enhanced.

ϕ for \(\alpha\).

ϕ for \(\xi\).

ϕ for \({\gamma }_{1}\).

ϕ for \(Sc\).

ϕ for \(E\).

ϕ for \(\Omega\).

The current study is performed to analyze the 2-dimensional peristaltic flow through porous walled channel in the presence of magnetic field. Energy equation is modelled by considering heat source/sink parameter. Whereas, activation energy is considered for the analysis of concentration profile. In addition, slip conditions are imposed to study flow quantities. The main outcomes of the analysis are:

Velocity shows parabolic behavior near the center of channel.

Decreasing behavior is noticed for slip parameter \(\gamma\) whereas enhancement is noticed for \(k\) and \({\lambda }_{1}\).

Similar profile is exhibited by temperature and velocity.

Temperature is enhanced for \(B\) and \({\gamma }_{1}\) whereas decreases for \(\alpha\).

Heat transfer and pressure gradient are of sinusoidal nature.

Pressure rise shows increasing behavior for \(\gamma , k\) and \({\lambda }_{1}\) in co-pumping region.

Concentration depicts increase for higher values of \({\gamma }_{2}\) and α whereas show decrease for \(Sc\) and \(\xi\).

The study presented in this article finds interesting application in various human physiological systems. Majority organs within human body show porous nature. In addition, presence of fluids make boundaries slippery in nature. Thus, we can say that the mathematical model developed in current analysis can be used to predict functioning of variety of systems. From future perspective, addition of nanoparticles, consideration of viscous dissipation with thermal radiations can be taken as mathematical model to study the treatment of cancer within any physiological system. Also, taking convective conditions at the boundaries in addition to Soret/Dufour effects may help to develop model for thermal analysis of digestive system. Thus, in the culmination of above analysis peristalsis finds numerous applications in multiple disciplines.

The data used to support the findings of this study are available from the corresponding author upon request.

Dimensional/non-dimensional wall surface

Time

Wavelength

Width of the channel, wave amplitude

Speed of wave

Cartesian coordinates

Dimensional/non-dimensional velocity components along horizontal and vertical directions

Dimensional/non-dimensional suction/injection velocity

Fluid density

Pressure

Electrical conductivity/magnetic field strength

Stress tensor

Specific heat capacity

Dimensional/non-dimensional Temperature

Variable thermal conductivity

Constant heat source/sink parameter

Dimensional/non-dimensional concentration

Mass diffusion coefficient

Chemical reaction parameter

Reaction rate

Dimensional/non-dimensional activation energy

Fitted rate constant

Boltzman constant

Thermal conductivity at constant temperature

Constant

Dimensional/non-dimensional velocity, temperature and concentration slip

Jeffrey fluid parameters

Wave number

Reynolds number

Hartman number

Prandtl number

Heat source/sink parameter

Schmidt number

Chemical reaction parameter

Temperature ratio parameter

Dimensionless wave amplitude

Heat transfer rate

Instantaneous flow rate in fixed frame

Flow rate in wave frame

Merriam-Webster's. Medical Dictionary (Merriam-Webster's, Inc, 1995).

Google Scholar

Bayliss, W. M. & Starling, E. H. The movements and innervation of the small intestine. J. Physiol. 24, 99–143 (1899).

Article CAS PubMed PubMed Central Google Scholar

Latham, T. W. Fluid Motion in a Peristaltic Pump (MIT, 1966).

Google Scholar

Shapiro, A. H. Pumping and retrograde diffusion in peristaltic waves, In Proceedings of the Workshop on Uretal Reflux in Children 109–126 (National Academy of Sciences, 1967).

Casson, N. A flow equation for pigment oil-suspensions of the printing ink type. In Rheology of Disperse System (ed. Mill, C. C.) 84–105 (Pergamon Press, 1959).

Google Scholar

Eringen, A. C. Nonlinear theory of simple micro-elastic solids-I. Int. J. Eng. Sci. 2, 189–203 (1964).

Article MathSciNet MATH Google Scholar

Ramamurthy, J. V., Prasad, K. M. & Narla, V. K. Unsteady peristaltic transport in curved channel. Phys. Fluids 25, 091903 (2013).

Article ADS MATH Google Scholar

Nadeem, S., Sadaf, H. & Akbar, N. S. Analysis of peristaltic flow for a Prandtl fluid model in an endoscope. J. Power Technol. 94, 1–11 (2014).

Google Scholar

Sadeghi, K. & Talab, H. J. Analytical investigation of peristaltic transport of power law fluid through a tube. Appl. Mech. Eng. https://doi.org/10.4172/2168-9873.1000136 (2014).

Article Google Scholar

Tripathi, D., Akbar, N. S., Khan, Z. H. & Anwar Beg, O. Peristaltic transport of bi-viscosity fluids through a curved tube: A mathematical model for intestinal flow. Proc. Inst. Mech. Eng. H 230, 817–827 (2016).

Article PubMed Google Scholar

Fusi, L. & Farina,. Peristaltic axisymmetric flow of a Bingham fluid. Appl. Math. Comput. 320, 1–15 (2018).

MathSciNet MATH Google Scholar

Ramesh, K. & Devakar, M. Effect of endoscope on the peristaltic transport of a couple stress fluid with heat transfer: Application to biomedicine. Nonlinear Eng. 8, 619–629 (2019).

Article ADS Google Scholar

Vaidya, H. et al. Peristaltic flow of non-Newtonian fluid through an inclined complaint nonlinear tube: Application to chyme transport in the gastrointestinal tract. Eur. Phys. J. Plus 135, 934–948 (2020).

Article CAS Google Scholar

Ferraro, V. C. A. & Plumpton, C. An Introduction to Magneto-Fluid Mechanics (Clarendon Press, 1966).

MATH Google Scholar

Alexiou, C. et al. Loco regional cancer treatment with magnetic drug targeting. Cancer Res. 60, 6641–6648 (2000).

CAS PubMed Google Scholar

Carlton, J. M. R., Yowell, C. A., Sturrock, K. A. & Dame, J. B. Bio magnetic separation of contaminating host leukocytes from plasmodium infected erythrocytes. Exp. Parasitol. 97, 111–115 (2001).

Article CAS PubMed Google Scholar

Oleson, J. R. Hyperthermia by magnetic induction: I. Physical characteristics of the technique. Int. J. Radiat. Oncol. Biol. Phys. 8, 1747–1756 (2001).

Article Google Scholar

Voltairas, P. A., Fotiadis, D. I. & Michalis, M. K. Hydrodynamics of magnetic drug targeting. J. Biomech. 35, 813–820 (2002).

Article CAS PubMed Google Scholar

Hummady, L. & Abdulhadi, A. Inuence of MHD on peristaltic flow of couple-stress fluid through a porous medium with slip effect. Adv. Phys. Theor. Appl. 30, 34–44 (2014).

Google Scholar

Satyanarayana, K. V. V., Sreenadh, S., Lakshminarayana, P. & Sucharitha, G. MHD peristaltic transport of a micropolar fluid in an asymmetric channel with porous medium. Adv. Appl. Sci. Res. 7, 105–114 (2016).

Google Scholar

Latif, T., Alvi, N., Hussain, Q. & Asghar, S. Variable properties of MHD third order fluid with peristalsis. Res. Phys. 6, 963–972 (2016).

Google Scholar

Prakash, J., Balaji, N., Siva, E. P. & Chanrasekaran, A. D. Non-linear blood flow analysis on MHD peristaltic motion of a Williamson fluid in a micro channel. AIP Conf. Proc. 2112, 020048. https://doi.org/10.1063/1.5112233 (2019).

Article CAS Google Scholar

Selvi, C. K. & Srinivas, A. N. S. Peristaltic transport of Herschel–Bulkley fluid in a non-uniform elastic tube. Propul. Power Res. 8, 253–262 (2019).

Article Google Scholar

Shera, M. D., Gladman, A. S., Davidson, S. R., Trachtenberg, J. & Gertner, M. R. Helical antenna arrays for interstitial microwave thermal therapy for prostate cancer: Tissue phantom testing and simulations for treatment. Phys. Med. Biol. 46, 1905–1918 (2001).

Article CAS PubMed Google Scholar

Sanyal, D. C. & Maji, N. K. Thermoregulation through skin under variable atmospheric and physiological conditions. J. Theor. Biol. 208, 451–456 (2001).

Article ADS CAS PubMed Google Scholar

Hayat, T., Nisar, Z. & Alsaedi, A. Impacts of slip in radiative MHD peristaltic flow of fourth grade nanomaterial with chemical reaction. Int. Commun. Heat Mass Transf. 119, 104976 (2020).

Article CAS Google Scholar

Tanveer, A., Mahmood, S., Hayat, T. & Alsaedi, A. On electroosmosis in peristaltic activity of MHD non-Newtonian fluid. Alex. Eng. J. 60, 3369–3377 (2021).

Article Google Scholar

Tanveer, A. & Malik, M. Y. Slip and porosity effects on peristalsis of MHD Ree–Eyring nanofluid in curved geometry. Ain Shams Eng. J. 12, 955–968 (2021).

Article Google Scholar

Javed, M., Imran, N., Arooj, A. & Sohail, M. Meta-analysis on homogeneous–heterogeneous reaction effects in a sinusoidal wavy curved channel. Chem. Phys. Lett. 763, 138200 (2021).

Article CAS Google Scholar

Ahmed, B., Hayat, T., Abbasi, F. M. & Alsaedi, A. Mixed convection and thermal radiation effect on MHD peristaltic motion of Powell–Eyring nanofluid. Int. Commun. Heat Mass Transf. 126, 105320 (2021).

Article CAS Google Scholar

Noreen, S., Waheed, S., Lu, D. C. & Tripathi, D. Heat stream in electroosmotic bio-fluid flow in straight microchannel via peristalsis. Int. Commun. Heat Mass Transf. 123, 105180 (2021).

Article Google Scholar

Akram, S., Athar, M. & Saeed, K. Hybrid impact of thermal and concentration convection on peristaltic pumping of Prandtl nanofluid in non-uniform inclined channel and magnetic field. Case Stud. Therm. Eng. 25, 100965 (2021).

Article Google Scholar

Khazayinejad, M., Hafezi, M. & Dabir, B. Peristaltic transport of biological graphene-blood nanofluid considering inclined magnetic field and thermal radiation in a porous media. Powd. Technol. 384, 452–465 (2021).

Article CAS Google Scholar

Imran, N., Javed, M., Sohail, M. & Tlili, I. Simultaneous effects of heterogeneous–homogeneous reactions in peristaltic flow comprising thermal radiation: Rabinowitsch fluid model. J. Mater. Res. Technol. 9, 3520–3529 (2020).

Article CAS Google Scholar

Imran, N. et al. Influence of chemical reactions and mechanism of peristalsis for the thermal distribution obeying slip constraints: Applications to conductive transportation. J. Mater. Res. Technol. 9, 6533–6543 (2020).

Article CAS Google Scholar

Imran, N., Javed, M., Sohail, M., Thounthong, P. & Abdelmalek, Z. Theoretical exploration of thermal transportation with chemical reactions for Sutterby fluid model obeying peristaltic mechanism. J. Mater. Res. Technol. 9, 7449–7459 (2020).

Article CAS Google Scholar

Fung, Y. C. & Tang, H. T. Longitudinal dispersion of particles in the blood flowing in a pulmonary alveolar sheet. ASME J. Appl. Mech. 42, 536 (1975).

Article ADS Google Scholar

Gopalan, N. P. The Pulsatile blood flow in a rigid pulmonary alveolar sheet with porous walls. Bull. Math. Biol. 43, 563577 (1981).

Article MathSciNet MATH Google Scholar

Imran, N., Javed, M., Sohail, M., Gokul, K. C. & Roy, P. Exploration of thermal transport for Sisko fluid model under peristaltic phenomenon. J. Phys. Commun. 4, 065003. https://doi.org/10.1088/2399-6528/ab9557 (2020).

Article CAS Google Scholar

Imran, N., Javed, M., Sohail, M., Farooq, S. & Qayyum, M. Outcome of slip features on the peristaltic flow of a Rabinowitsch nanofluid in an asymmetric flexible channel. Multidiscip. Model. Mater. Struct. 17, 181–197 (2021).

Article Google Scholar

Li, S., Zhao, M., Xue, J. & Zhao, R. Effects of edge type and reconstruction on the electronic properties and magnetism of 1T-ReS2 nanoribbons: A study based on DFT calculations. J. Mag. Magn. Mater. 567, 170351 (2022).

Article Google Scholar

Zhang, Y., Liu, G., Ye, J. & Lin, Y. Crushing and parametric studies of polygonal substructures based hierarchical cellular honeycombs with non-uniform wall thickness. Compos. Struct. 299, 116087 (2022).

Article Google Scholar

Bai, B., Zhou, R., Yang, G., Zou, W. & Yuan, W. The constitutive behavior and dissociation effect of hydrate-bearing sediment within a granular thermodynamic framework. Ocean Eng. 268, 113408 (2023).

Article Google Scholar

Zhang, L. et al. Effect of bimetallic modification on blast furnace slag and its application in low-temperature selective catalytic reduction. J. Chemical Technol. Biotechnol. https://doi.org/10.1002/jctb.7298 (2022).

Article Google Scholar

Zhang, X. et al. Preparation of PI porous fiber membrane for recovering oil-paper insulation structure. J. Mater. Sci. Mate. Electr. 31, 13344–13351 (2020).

Article CAS Google Scholar

Zhao, G., Shi, L., Yang, G., Zhuang, X. & Cheng, B. 3D fibrous aerogels from 1D polymer nanofibers for energy and environmental applications. J. Mater Chem. A 11, 512–547 (2023).

Article CAS Google Scholar

Bai, B., Wang, Y., Rao, D. & Bai, F. The effective thermal conductivity of unsaturated porous media deduced by pore-sScale SPH simulation. Front. Earth Sci. https://doi.org/10.3389/feart.2022.943853 (2022).

Article Google Scholar

Bai, B., Rao, D., Chang, T. & Guo, Z. A nonlinear attachment-detachment model with adsorption hysteresis for suspension-colloidal transport in porous media. J. Hydrol. (Amsterdam) 578, 124080 (2019).

Article CAS Google Scholar

Zhang, X. et al. Radionuclide transport in multi-scale fractured rocks: A review. J. Hazard. Mater. 424, 127550 (2021C).

Article PubMed Google Scholar

Zhang, X. et al. Application of upscaling methods for fluid flow and mass transport in multi-scale heterogeneous media: A critical review. Appl. Energy 303, 117603 (2021).

Article CAS Google Scholar

Batool, S. et al. Numerical analysis of heat and mass transfer in micropolar nanofluids flow through lid driven cavity: Finite volume approach. Case Stud. Therm. Eng. 37, 102233. https://doi.org/10.1016/j.csite.2022.102233 (2022).

Article Google Scholar

Rasool, G. et al. Significance of Rosseland's radiative process on reactive Maxwell nanofluid flows over an isothermally heated stretching sheet in the presence of Darcy–Forchheimer and Lorentz forces: Towards a new perspective on Buongiorno's model. Micromachines 13, 368. https://doi.org/10.3390/mi13030368 (2022).

Article PubMed PubMed Central Google Scholar

Rasool, G. et al. Darcy–Forchheimer flow of water conveying multi-walled carbon nanoparticles through a vertical cleveland Z-staggered cavity subject to entropy generation. Micromachines 13, 744. https://doi.org/10.3390/mi13050744 (2022).

Article PubMed PubMed Central Google Scholar

Rasool, G., Ali Shah, N., El-Zahar, E. R. & Wakif, A. Numerical investigation of EMHD nanofluid flows over a convectively heated Riga pattern positioned horizontally in a Darcy–Forchheimer porous medium: application of passive control strategy and generalized transfer laws. Waves Random Complex Media https://doi.org/10.1080/17455030.2022.2074571 (2022).

Article Google Scholar

Hayat, T., Farooq, S., Ahmad, B. & Alsaedi, A. Consequences of variable thermal conductivity and activation energy on peristalsis in curved configuration. J. Mol. Liq. 263, 258–267 (2018).

Article CAS Google Scholar

Hayat, T., Abbasi, F. M., Ahmad, B. & Alsaedi, A. MHD mixed convection peristaltic flow with variable viscosity and thermal conductivity. Sains Malays. 43, 1583–1590 (2014).

Google Scholar

Sivaiah, R. & Reddy, R. H. Magnetohydrodynamic peristaltic motion of a Newtonian fluid through porous walls through suction and injection. IOP Conf. Ser. Mater. Sci. Eng. 263, 062007 (2017).

Article Google Scholar

Ramesh Babu, V., Sreenadh, S. & Srinivas, A. N. S. Peristaltic transport of a viscous fluid in a porous channel with suction and injection. Ain Shams Eng. J. 9, 909–915 (2018).

Article Google Scholar

Nisar, Z., Hayat, T., Alsaedi, A. & Ahmad, B. Significance of activation energy in radiative peristaltic transport of Eyring–Powell nanofluid. Int. Commun. Heat Mass Transf. 116, 104655 (2020).

Article CAS Google Scholar

Vaidya, H. et al. Influence of transport properties on the peristaltic MHD Jeffrey fluid flow through a porous asymmetric tapered channel. Res. Phys. 18, 103295 (2020).

Google Scholar

Download references

The authors are thankful to the Deanship of Scientific Research-Research Center at King Khalid University in Saudi Arabia for funding this research (Code Number: RGP 2/23/43).

Department of Mathematics, COMSATS University Islamabad, Attock, 43600, Pakistan

Maimona Rafiq & Asma Shaheen

Department of Physics, College of Science and Arts at Muhayel, King Khalid University, Abha, Saudi Arabia

Youssef Trabelsi

Center of Research, Faculty of Engineering, Future University in Egypt New Cairo, New Cairo, 11835, Egypt

Sayed M. Eldin

Department of Mechanical Engineering, Lebanese American University, Beirut, Lebanon

M. Ijaz Khan

Department of Mathematics and Statistics, Riphah International University, I-14, Islamabad, 44000, Pakistan

M. Ijaz Khan

Mechanical Engineering Department, College of Engineering and Islamic Architecture, Umm Al-Qura University, P. O. Box 5555, Mecca, 21955, Saudi Arabia

Dhia Kadhm Suker

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

M.R. presented the conceptM.R. and A.S. modeled and solved the problemY.T. and D.S. analysis of data and Figure Preparation in the revised manuscript D.S. Reviewed revised conceptsM.K. write-up the manuscriptAll authors reviewed the revised manuscript before submission.

Correspondence to Maimona Rafiq.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Rafiq, M., Shaheen, A., Trabelsi, Y. et al. Impact of activation energy and variable properties on peristaltic flow through porous wall channel. Sci Rep 13, 3219 (2023). https://doi.org/10.1038/s41598-023-30334-3

Download citation

Received: 26 October 2022

Accepted: 21 February 2023

Published: 24 February 2023

DOI: https://doi.org/10.1038/s41598-023-30334-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.