MHD flow study of viscous fluid through a complex wavy curved surface due to biomimetic propulsion under porosity and second-order slip effects

Pei-Ying Xiong,Khurram Javid,Mohsin Raza,Sami Ullah Khan,M Ijaz Khan and Yu-Ming Chu

1 College of Science,Hunan City University,Yiyang 413000,China

2 Department of Mathematics,Northern University,Wattar-Walai Ziarat,Kaka Saib Road,24110,Nowshera,KPK,Pakistan

3 Department of Mathematics,COMSATS University Islamabad,Sahiwal 57000,Pakistan

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

5 Department of Mathematics,Huzhou University,Huzhou 313000,China

Abstract The steady laminar flow of viscous fluid from a curved porous domain under a radial magnetic field is considered.The fluid flow by a curved domain is due to peristaltic waves present at the boundary walls.The whole analysis is based on porosity(Darcy number)effects.Moreover,the effects of second-order slip on the rheology analysis are also discussed.Due to the complex nature of the flow regime,we have governed the rheological equations by using curvilinear coordinates in the fixed frame.The physical influence of magnetic(Hartmann number)and porosity(Darcy number)parameters on the rheological features of peristaltic transportation are argued in detailed(in the wave frame).Additionally,in the current study,the complex wavy pattern on both boundary walls of the channel is used.The whole rheological study is based on ancient,but medically valid,assumptions of creeping phenomena and long wavelength assumptions.Analytical solutions of the governing equations are obtained by using the simple integration technique in Mathematica software 11.0.The core motivation of the present analysis is to perceive the physical influence of embedded parameters,such as the dimensionless radius of the curvature parameter,magnetic parameter,porosity parameter,different amplitude ratios of complex peristaltic waves,first-and second-order slip parameters,on the axial velocity,pressure gradient,local wall shear stress,tangential component of the extra-stress tensor,pumping and trapping phenomena.

Keywords:peristaltic waves,magnetic field,porous medium,pumping phenomena,secondorder slip condition

Nomenclature

1.Introduction

Peristaltic propulsion is a biological phenomenon in which the motion of physiological fluids takes place through physiological vessels and arteries due to the continuous contraction and relaxation of these small vessels and arteries,such as the movement of a small food bolus through the oesophagus,urine transportation through the ureter,the motion of blood through blood vessels,air and gas movement through small arteries of the lung system,the movement of food from the mouth to the stomach,transportation of bile in bile ducts,intestinal chyme movement,ovum movement in the fallopian tube,and other hollow tubes.In these phenomena,physiological fluids are moved from one part of the vessel where the pressure is higher to another part of the vessel where the pressure is smaller by a continuous process of muscle contraction and relaxation.The periodic transportations of these vessels are actuated due to electro–chemical responses occurring in the body.Many industrial pumps,such as micro-pumps,mini-pumps,handy-pumps,integrating filling systems,roller pumps,compact peristaltic pumps,and standard peristaltic pumps,are based on the peristaltic pumping phenomenon.Applications in engineering fluid mechanics utilize caustic fluid(atomic industries),destructive chemicals,high solid slurries and numerous materials which are transported by peristaltic pumps.In the mid 1960s,Latham[1]first studied the mathematical modelling of peristaltic pumping with a seminal study(both analytically and experimentally)on biomechanical pumping.Two years after Latham’s work[1],in 1968,Fung and Yih[2]theoretically performed the same study as Latham[1].After these two studies,a new horizon was developed in the research community and researchers performed dynamic research in this domain.In 1969,Shapiro et al[3]studied the mathematical modelling of the peristaltic motion of a viscous fluid and verified the fundamental findings of Latham[1].Their analysis is based upon bio assumptions and arguments over the reflux and trapping phenomena during peristalsis.In this era,some fundamental contributions were performed,as given in[4–6].

Many industrial instruments and systems,including smart lubrication systems[7],nanomaterial fabrication[8],biomimetic pumps[9],separation devices,and magnetic blood pumps[10],are based on the study of magneto-hydrodynamics(MHD)to control the flow system,to acquire the performance augmentation,material influence,etc.A massive range of phenomena emerges in flows,such as the Hall current,Ohmic dissipation,ferromagnetism,oscillation,Landau damping,and oscillation.The utilization of magnetic fields in electrically conducting streams are produced by Lorentz powers(Lorentz forces),which can control rheological features.

Over the past two decades,the mathematical formulations of the MHD flow of numerous fluids due to peristaltic phenomena have gained considerable attention.The magnetic force has diverse applications,especially in MRI,that is used nowadays in hospitals to control blood motion(haematological diseases)during complex surgeries.Tripathi and Beg[11]have studied the magnetic and porosity effects,respectively,on the peristaltic motion of couple stress bio-fluid through coaxial channels.Ali et al[12]discussed the peristaltic propulsion of Newtonian fluid through a two-dimensional channel with variable viscosity.Furthermore,they have argued the rheological features under the influence of magnetic field and slip effects in detail under bio assumptions.Hayat and Ali[13]have deliberated the rheological study of third grade fluid through a deformable tube under a magnetic field.They have obtained the analytical solution by using the perturbation technique.In another research study,Hayat and Ali[14]have discussed the rheological analysis related to the peristaltic transportation of a viscoelastic Jeffrey fluid through a tube with a sinusoidal wave travelling down its boundary walls.The whole rheological analysis is based on a uniform magnetic field.Further productive research studies in this area are given in[15–17].

The internal motion of physiological fluids through a confined porous channel with contraction and relaxation walls has significant applications in biological flows,for example,blood flow through blood vessels,artificial dialysis,and air circulation in the lung system.Further applications include pulsating diaphragms,cerebral hydrodynamics,oxygen diffusion in the capillaries,the mechanics of the cochlea in the human ear,filtration in tissue,and cosmetics material manufacturing.Hayat et al[18]have performed dynamic research and highlighted the influence of partial slip and porosity effects on the biomimetic motion of Newtonian fluid.They obtained the solution for the stream function and axial velocity analytically by using the Adomian decomposition method.Alsaedi et al[19]examined the transportation of an incompressible couple stress bio-fluid through a porous channel under creeping phenomena.They attained the exact solutions of a system of rheological equations.Hayat et al[20]have argued over the peristaltic flow of Maxwell fluid in a planar channel under uniform porosity and Hall effects,respectively.They modelled the governing equations with the help of modified Darcy law(MDL).An impressive amount of research in this regard has been performed by some researchers and cited in[21–23].

The motion of Maxwell fluid through a moving plate has been studied by Lin and Guo[24]under magnetic effects.They have considered the second-order slip between the walls and the fluid.An analytical solution of a system of rheological equations is acquired by using Laplace transformation.They also obtained the solution for the simple slip(first-order slip)condition.Aly and Ebaid[25]have performed a flow study which deals with the peristaltic motion of nanofluid through an asymmetric channel under the influence of second-order velocity slip.They have obtained the exact solutions of the flow equations in terms of the stream function.Many mathematicians have performed dynamic research concerning the simple slip,and the second-order slip is cited in[26–28].

To the best of our knowledge,in 2000,Sato et al[29]were the first to perform productive research by using the curved channel.They have discussed the motion of twodimensional Newtonian fluid through a curved domain in the fixed frame.They obtained the rheological equations by using curvilinear coordinates because of the complex shape of flow geometry.After the productive research work of Sato et al[29],in 2008,Ali et al[30]prolonged the study of Sato et al[29]in the wave frame.They solved the Navier–Stokes equations analytically with the help of an integration technique.In these studies[29,30],their main concern was to focus the sound effects of the curvature parameter on a Newtonian fluid flow.Their whole research studies are based on bio approximations,respectively.These two research studies generated a new method of research for mathematicians.Researchers performed good research studies by using the different nature of fluids through the curved channel under numerous physical effects,such as,magnetic effects,electric effects,porosity effects,complaint wall effects,electro–magneto–hydrodynamic effects,heat transfer effects,heat and mass transfer effects,etc,as fluid flows because of the sinusoidal contraction and relaxation of boundary walls.The physical influence of heat transfer on the peristaltic transportation of a viscous fluid through a curved domain has been studied by Ali et al[31].They obtained the numerical solution by using the Runge–Kutta algorithm.They observed the effects of the curvature parameter(physical parameter related to flow geometry)and Brinkman number(the heat flux parameter)on numerous rheological features.They also provided a comparison between curved and straight channels.They noticed that the rate of heat transfer in the curved domain is smaller than the straighter channel.Tanveer et al[32]presented an investigation related to the numerical simulation for the transportation of a Carreau–Yasuda nanofluid model through a curved channel.They also discussed the porosity effects by using MDL.After these productive studies,some researchers[33–38]discussed the rheological characteristics of numerous fluids through curved geometries under the lubrication approach and bio assumption,respectively,where they highlighted the influence of various physical forces,such as electric forces,magnetic forces,electro–magneto–hydrodynamics forces,heat transfer effects and mass transfer effects by using the different nature of fluids that flow through the curved geometry.Their main intention was to highlight the physical influence of the curvature parameter on numerous rheological features.They solved the rheological equations both analytically(the integration method,perturbation technique,singular perturbation technique)and numerically(the shooting method,Runge–Kutta technique,finite difference technique,finite element technique,BVP4C(boundary value problem)technique).

To the best of our knowledge,no studies have so far been performed with regard to the peristaltic transport of a viscous fluid from a curved channel in porous medium.Complex sinusoidal waves are present at the boundary walls of the flow regime.The main objective of the present manuscript is therefore to scrutinize the influence of the dimensionless curvature parameter,porous media,second-order slip parameter and magnetic field on the rheological characteristics of Newtonian fluid under creeping phenomena.An analytical solution is developed with an integration technique in Mathematica software 11.0.This is the novelty of the current study.

2.Mathematical formulation

In the current study,we have considered a curved channel that is illustrated in figure 1.When the boundary walls of the curved channel are brought under the influence of peristaltic waves,initially,one section of the wall starts to contract at the inlet and the component positioned immediately ahead of this starts to relax,and this propagation of contraction and relaxation is continuous throughout the length of boundary walls.Finally,this propagation of the wave moves towards the outlet.This phenomenon continues until the whole transportation takes place.The fluid which is present inside the confined domains of the curved channel is initially at rest during the boundary walls when they are in a stationary state.But as the peristaltic waves are generated at the boundary walls,then fluid starts to move with the same speed and direction as the peristaltic waves.Additionally,the sinusoidal peristaltic wave existing at the boundary walls is complex in shape.This type of pattern is used to increase the efficiency and accuracy of roller and peristaltic pumps,respectively.For five years,this complex pattern of waves has been attracting the attention of researchers.Some dynamic research has been performed by various mathematicians and researchers by using the complex wavy pattern of peristaltic waves and obtaining productive outcomes that are mentioned in[39,40].Additionally,a radial magnetic field is applied to the rheostat of the fluid transportation.

The mathematical expression of the complex peristaltic wave that propagates along the geometric walls is defined as:

MHD is an outlet of applied physics that deals with the relation between the magnetic field and electrically conducting fluids,for example,electrolytes,liquid metals,magnetofluids,including plasmas and briny water.Physically,the essential thought behind MHD is that when the magnetic field is applied to a moving conductive fluid,then the fluid is polarized and,reciprocally,it deviates from the magnetic field.The external force exerted on a moving charge particleQwith a velocity fieldthrough both electricand magneticfields,respectively,is called the Lorentz force.The mathematical expression of the Lorentz force is defined below:

whereσis the fluid electrical conductivity.In the current manuscript,we will consider the weaker strength of an electric field,say=0.Then,the above equation will reduce to the following form:

Using equation(5)in(4),we get

According to Maxwell law,we have

The mathematical expression of the Cauchy-stress tensorfor a viscous fluid is defined as:

wherePis the pressure term,is the identity tensor,andis the extra-stress tensor.The vector form of the extra-stress tensor for the Newtonian fluid model is defined as:

Fluid mechanics has three fundamental equations:(i)continuity equation,(ii)momentum equation,and(iii)energy equation,respectively.These three dynamic equations provide the physics of dynamic systems.Their mathematical statement is based upon the three fundamental laws of physics:the law of conservation of mass,momentum and energy.These equations are expressed in a differential form of time and space coordinates,respectively.

The vector representations of rheological equations are:

where,

The mathematical forms of the rheological equations in terms of curvilinear coordinates are expressed as:

In the above equations,the mathematical termsare the components of the extra-stress tensor andυis the dynamic viscosity.The dimensional form of the components of the extra-stress tensor are given below:

In this section,we will introduce a linear transformation between fixed and wave frames,respectively.The wave frame is moving with the wave speed‘a(chǎn)’.The main cause of linear transformation is to shift the flow system from an unsteady state to a steady state.The mathematical expressions of these linear relations are defined as:

The non-dimensional or dimensionless numbers and physical parameters are those which are significant in order to explain the numerous rheological features of fluids.These dimensionless numbers are a combination of numerous physical characteristics,such as force,density,speed,length parameter,etc.The Hartmann number(magnetic parameter)and Reynolds number are two famous dimensionless numbers.Another main advantage of these dimensionless numbers is that the quantity of the physical embedded variables is condensed.

Introducing the dimensionless variables:

The stream function is a scalar function.Mathematically,these scalar functions are expressed in terms of the derivatives of time and space coordinates,respectively.Each derivative term represents a velocity component orthogonal to any direction.These stream functions give information about the rotational or irrotational flow of a fluid.If a fluid satisfies Laplace’s equation,it is said to be an irrotational fluid:

wherexis the axial component of geometry,ris the radial component of geometry,vis the radial velocity,uis the axial velocity,Reis the Reynolds number,his the upper half of the regime,δis the wave number,Ω is the dimensionless radius of curvature,Hais the dimensionless magnetic parameter(Hartmann number),Dais the dimensionless porosity parameter(Darcy number),(α,γ)are the physical parameters,(?1,?2)are the dimensionless symbols of the different amplitudes of complex peristaltic waves,Pis the pressure term,τijis the dimensionless form of extra-stress tensor,(β1,β2)are the first-and secondorder slip parameters,respectively.After using these dimensionless variables and stream functions,we apply bio assumptions into a system of equations in a wave frame.These assumptions are called low Reynolds numbers(where the viscous force is much larger than the inertial force,and mathematically it is denoted asRe→0)and long wavelength(where the wavelength of the sinusoidal waves is much larger than the half width of the channel,and mathematically it is denoted asλ?A)approximations.Then,we get

It is obvious from equation(22)that the pressure term is independent of the radial coordinate.In equation(23),two dependent variables are present:one is the pressure term and the second is the stream function.In order to eliminate the pressure from equation(23)by using cross derivatives between equations(22)and(23),we obtain:

An alternative form of equation(24)is

Equation(26)is the required fourth order ODE that describes the rheological topographies of a Newtonian fluid by the curved geometry in terms of the stream function.

The second-order slip boundary conditions(see[24–28])in terms of the stream function are:

where F is the dimensionless flow rate,β1is the first-order slip parameter andβ2is the second-order slip parameter.

The pressure rise can be expressed as the mathematical form defined as[30–35]:

The mathematical term is related to the pressure gradient,which is present in the above equation(29)as dp/dx.The pressure gradient defines the rate of change(derivative or slope)of the pressure term w.r.t the axial coordinate,while ΔPis the pressure rise per wavelength in each cross-section.

The analytical solution of equation(26)is obtained by using an integration technique in the Mathematica software 11.0:

The axial velocity is obtained by taking the partial differential of equation(30)w.r.t r:

whereAi(i=1 -4)are the constant of integrations.These constants can be obtained by utilizing Mathematica software 11.0.

The local wall shear stress at the upper wall of the channel is:

Taking partial differential of equation(32)w.r.t r:

atr=h

Special Case:

The rheological results of[30]are retrieved by neglecting the magnetic effects,porosity effects and complex wavy scenario of peristaltic waves from the rheological equations.

All these results give authenticity to the correctness and validation of obtaining governing equations.The range of physical parameters are given in table 1.

Table 1.Range of physical parameters.

3.Graphical consequences and interpretations

The transport problem amounts to a 4th-order system of linear and ODE defined by equations(22)–(23)with BC’s in equations(27)and(28).This BVP is solved analytically using the four-time integration in Mathematica.Validation of the current flow analysis has been conducted with existing results in the literature.The results of[30]are acquired by eliminating some physical flow features.

In this section,we shall examine the physical influence of involved parameters such as the curvature parameter Ω,different amplitude ratios of complex peristaltic waves(?1,?2),the Hartmann numberHa,Darcy numberDa,first-and second-order slip parameters(β1,β2)on the axial velocity,pressure gradient,pumping phenomena,extra shear stress tensor,the local shear stress at the upper walls,and trapping phenomena,respectively with the help of graphical illustration.Figures 2–4 illustrate selected flow distributions of the physical variables,which is,axial velocity,pressure gradient,pressure rise,component of the extra stress tensor,local shear stress,and trapping phenomena with radial and axial coordinates,respectively.

Figure 1.Flow geometry.

Figure 2.The variation of the velocity proflie for the first-order slip parameter(β1) atΩ=3.5,?1 =0.3,?2 =0.4,x=π,A=1,Ha=1,Da=0.1,Q =1.

Figure 3.The variation of the velocity profile for the first-order slip parameter(β2)atΩ=3.5,?1 =0.3,?2 =0.4,x=π,A=1,Ha=1,Da=0.1,Q =1.

3.1.Axial velocity profile

Figures 2–5 show the impacts of the Hartmann numberHa,Darcy numberDa,first-and second-order slip parameters(β1,β2),i.e.,with the radial coordinate.The magnitude of the axial velocity of Newtonian fluid through a curve is decreased by increasing the length of both slip parameters(see figures 2 and 3).These results are obtained under the influence of two external forces:one is the magnetic force and the other is the porous medium.The difference between these two figures is that the magnitude of the axial velocity sharply decreases by enhancing the second-order slip parameter from 0.1 to 2 and fixing the first-order slip parameter,sayβ1=0.1.Physically,these two results show that the magnitude of the velocity profile at both the walls is larger in the absence of both slip parameters as compared with the non-zero values of the slip parameters.Figure 4 depicts the impact of the Hartmann number on the velocity profile under larger porosity effects.At smaller strengths of magnetic field,say Ha=0.1,the magnitude of the velocity profile is larger.But as the numeric value of the Hartmann number changes from 0.1 to 2,not only is the magnitude of the velocity profile larger but it also shifts from the lower half to the upper half of the channel.Additionally,when the magnetic strength is enhanced from 2 to 10,then the magnetic effects are dominant over the porosity effects and the magnitude of the velocity profile reduces.These results show that a larger strength of magnetic field has a dynamic impact on the velocity profile and sharp changes are predicted under a larger magnetic strength.Here,the magnetic effects are dominant over the curvature effects.The physical influence of the Darcy number on the axial velocity is displayed in figure 5 under the larger strength of the magnetic field.Here,we have noticed two distinct results for smaller and larger values of the Darcy number.When Da changes from 0.1 to 0.5,then the amplitude of the axial velocity is enhanced near the boundary walls of the curved channel.Just like the previous result of Ha,it is also observed that the velocity profile shifts toward the lower half of the curved channel as the numeric value of Da changes from 0.5 to 2.Here,it can be noticed that the magnitude of the velocity profile is decreased at the lower half as Da increases from 2 to∞(no porosity).Physically,it means that under larger porosity effects,the magnitude of the velocity profile is smaller as compared with the weaker strength of the porosity effects.

Figure 4.The variation of the velocity profile for the magnetic parameter(Ha)at Ω=3.5,?1 =0.3,?2 =0.4,β1=0.01β2=0.1x=π,A=1,Da=50,Q =1.

Figure 5.The variation of the velocity profile for the porosity parameter(Da)atΩ=3.5,?1 =0.3,?2 =0.4,β1=0.01,β2=0.1x=π,A=1,Ha=1,Q =1.

3.2.Pressure gradient

An important rheological characteristic of peristaltic motion that is concerned with the rate of change of the pressure function along the axial direction is called the pressure gradient.The impacts of the Hartmann numberHa,Darcy numberDa,first-and second-order slip parameters(β1,β2),i.e.,with the axial coordinate,are visualized in figures 6–9.Due to the complex nature of peristaltic waves,the wavy behavior is predicted in the graphs of the pressure gradient at each flow rate.In this subsection,the physical impacts of the physical parameters on the pressure gradient are highlighted under magnetic,curvature,porosity,Hartmann number,firstand second-order slip effects,respectively.The effects of a simple slip parameter under both larger curvature and porosity effects are displayed in figure 6.According to the literature survey,slip parameters reduce the magnitude of the pressure gradient.A comparative kind of conduct is anticipated in figure 6.As the numeric value ofβ1is enhanced,the magnitude of the pressure gradient is reduced under larger porosity effects.Physically,the larger values of the first-order slip parameter are dominant over the curvature and porosity effects which is why the magnitude of the pressure gradient is reduced.Figure 7 deals with the effects of the second-order slip parameter on the pressure gradient under curvature,porosity and magnetic effects,respectively.Once again,the magnitude of the pressure gradient is reduced as the length of the second-order slip parameter is increased.Physically,the larger values of the second-order slip parameter are dominant over the curvature and porosity effects which is why the magnitude of the pressure gradient is reduced.Furthermore,it is perceived that the magnitude of the pressure gradient sharply decreases in figure 7 by enhancing the length ofβ2as compared with the results of figure 6.The physical effects of porosity(Darcy number)and magnetic(Hartmann number)parameters on the pressure gradient are plotted in figures 8 and 9,respectively.According to the available literature,the magnetic force plays a role in reducing the magnitude of the pressure gradient.Similar results are obtained in figure 8;the magnitude of the pressure gradient is strongly affected by increasing the magnetic force,while the opposite is obtained in figure 9.The magnitude of the pressure gradient is enhanced as the Darcy number is increased.Physically,it means that under larger porosity effects or smaller values of Darcy number,the magnitude of the pressure gradient is smaller,while the reverse pattern is obtained for smaller porosity effects.Furthermore,it is noted that under the larger strength of both magnetic and porosity parameters,the magnitude of the pressure gradient is smaller as compared with the weaker strength of the magnetic and porosity effects.

Figure 6.The variation of pressure gradient for the first-order slip parameter(β1) at Ω=2,?1 =0.3,?2 =0.4,A=1,Ha=1,Da=0.1,Q=1.

Figure 7.The variation of pressure gradient for the second-order slip parameter(β2)at Ω=2,?1 =0.3,?2 =0.4,A=1,Ha=1,Da=0.1,Q=1.

Figure 8.The variation of pressure gradient for the Hartmann number(Ha)atΩ=2,β1=β2=0.1,?1 =0.3,?2 =0.4,A=1,Da=50,Q =1.

Figure 9.The variation of pressure gradient for the Darcy number(Da)atΩ=2,β1=β2=0.1,?1 =0.3,?2 =0.4,A=1,Ha=1,Q=1.

Figure 10.The variation of pressure rise for the first-order slip parameter(β1) atΩ=2,?1 =0.3,?2 =0.4,A =1,β2=0.1,Ha=1,Da=0.1.

Figure 11.The variation of pressure rise for the second-order slip parameter(β2)at Ω=2,?1 =0.3,?2 =0.4,A=1,Ha=1,β1=0.1,Da=0.1.

Figure 12.The variation of pressure rise for the Hartmann number(Ha)atΩ=2,?1 =0.3,?2 =0.4,A =1,β1=β2=0.1,Da=0.1.

Figure 13.The variation of pressure rise for the Hartmann number(Da)atΩ=2,?1 =0.3,?2 =0.4,A =1,β1=β2=0.1,Ha=1.

3.3.Pumping phenomena

An important rheological phenomena related to peristalsis is known as the pumping phenomena.The mathematical definition of these phenomena is integrated with the pressure gradient term across a wavelength at each cross-sectional area.Graphically,a linear relation(graph)is obtained between the pressure rise and dimensionless flow rate parameter.The physical effects of the involved parameters,such as the dimensionless radius of curvature Ω,different amplitude ratios of complex peristaltic waves(?1,?2),Hartmann numberHa,Darcy numberDa,first-and second-order slip parameters(β1,β2),on the pumping phenomena are visualized in figures 10–13.The effects ofβ1(first-order slip parameter)on the pumping phenomena are plotted in figure 10 for the complex peristaltic transportation of viscous fluid through a curved channel.In the absence of the firstorder slip parameter,the magnitude of the pressure rise over the flow rate is maximum as compared with the presence ofβ1under porosity,magnetic and curvature effects.The magnitude of the pressure rise decreases by enhancing the first-order slip parameter from 0 to 0.3.A comparative kind of conduct is anticipated in figure 11 as observed in figure 10 for various values of the second-order slip parameter.Once again,the magnitude of the pressure rise decreases by enhancing the length ofβ2from 0 to 0.3.The main difference between these two figures is that the rate of change in the pumping phenomena is larger in figure 11(forβ2)as compared with the outcomes of figure 10(forβ1).Physically,the outcomes of both figures show that the slip effects overcome the magnetic and porosity effects.Furthermore,it is said that the pumping phenomena can be controlled with the help of the slip parameter.The magnetic effects are displayed in figure 12 on the pumping phenomena.In the absence of the Hartmann number,sayHa=0,the porosity and curvature effects are dominant.Due to these effects,the pressure rise is larger in magnitude for a viscous fluid.When the numeric value of the Hartmann number is increased from 0 to 100,then the magnitude of the rise in pressure is sharply reduced.Physically,the larger values of the magnetic field overcome the curvature and porosity effects,respectively.The physical effects of the porosity parameter are visualized in figure 13 on the pumping phenomena for the complex peristaltic motion of a viscous fluid under curvature,slip and magnetic effects,respectively.In this figure,we have observed the reverse behavior by enhancing the numeric values of the porosity parameter.For the smaller values of the Darcy number(when the effects of the porosity medium on the pumping phenomena are larger),the pressure rise has a smaller magnitude.Physically,the porosity effects are dominant over the curvature and magnetic effects.AsDaincreases from 0 to∞,the magnitude of the pumping phenomena is enhanced because of the weaker effects of the porous medium under larger curvature and magnetic effects.Additionally,the larger strength of the magnetic field and porosity medium play a dynamic role in the reduction of peristaltic pumping under curvature effects.

3.4.Shearing stress tensor

In this subsection,the graphs of the shearing component of the stress tensor are plotted through figures 14–17.The shearing stress tensor is also known as the tangential stress tensor.Figures 14–17 display the variation of the shearing stress tensor for different values of the Hartmann numberHa,Darcy numberDa,and first-and second-order slip parameters(β1,β2).The magnitude of the shearing stress tensor at both walls is decreased by increasing the first-order slip(figure 14)effects for the biomimetic transportation of a viscous fluid under curvature,porosity and magnetic effects,respectively.Figure 15 is concerned with the effects of the second-order slip parameter on the shear stress of a viscous fluid.The magnitude of the shear stress tensor at both boundary walls is reduced by enhancing the length of the second-order slip parameter(β2).Physically,in both figures 14 and 15,the magnitude of the shearing stress tensor sharply reduces by enhancing the strength of both slip parameters.The magnetic effects are visualized in figure 16 on the shearing stress tensor.The magnitude of the shearing stress tensor is larger in the absence of a magnetic field in the presence of both porosity and curvature effects,respectively.It is noticed that the magnitude of the shearing stress is reduced by increasing the magnetic force at larger numeric values of the Hartmann number.The boundary layer phenomena is observed for larger magnetic effects.The physical influence of the porous medium is displayed in figure 17 on the shear stress tensor.Here,the reverse behavior is predicted for the larger Darcy number.The magnitude of the shearing stress is increased at both the boundary walls by increasing the Darcy number.Additionally,the boundary layer phenomena is observed near the boundary walls for a smaller Darcy number.Physically,these two figures show that the larger strength of the magnetic field and porous medium play a dynamic role in producing the boundary layer phenomena immediate to the boundary walls.Additionally,these two forces play a dynamic role in controlling the behavior of the shearing stress tensor and larger values of both parameters overcome the curvature effects.

Figure 14.The variation of the shearing stress tensor for the firstorder slip parameter(β1) atΩ=3.5,β2=0.1,A=1,Ha=1,x=π,Q=1,Da =0.1.

Figure 15.The variation of the shearing stress tensor for the secondorder slip parameter(β2)atΩ=3.5,A =Ha =Q =1,β1=0.1,x=π,Da=0.1.

Figure 16.The variation of the shearing stress tensor for the magnetic parameter(Ha)atΩ=2,A =1,β1=β2=0.1,x=π,Q=1,Da=50.

Figure 17.The variation of the shearing stress tensor for the porosity parameter(Da)atΩ=3.5,A =Q=1,x =π ,β1=β2=0.1,Ha =1.

Figure 18.The variation of local wall shear stress for the first-order slip parameter(β1) atΩ=3.5,β2=0.1,A=1,Ha=1,x=π,Da =0.1.

Figure 19.The variation of local wall shear stress for the secondorder slip parameter(β2)atΩ=3.5,A=1,Ha =1,β1=0.1,x=π,Da=0.1.

Figure 20.The variation of local wall shear stress for the magnetic parameter(Ha)atΩ=3.5,A =1,β1=0.01,Q =1,β2=0.1,Da =1.

Figure 21.The variation of local wall shear stress for the porosity parameter(Da)atΩ=3.5,A =1,β1=0.01,Q =1,β2=0.1,Ha =1.

Figure 22.(a)–(d).The variation of streamlines forHa (0,1,5,50) atΩ=3.5,A =1,β1=β2=0.1,Q=1,Da=50.

Figure 23.(a)–(d).The variation of streamlines forDa(∞,1,0.2,0.02) atΩ=3.5,A =1,β1=β2=0.1,Q=1,Ha=0.02.

Figure 24.(a)–(d):The variation of streamlines forβ1 (0,1,2,5)atΩ=3.5,A =1,β2=0.1,Q=1,Da=0.1,Ha=0.1.

Figure 25.(a)–(d):The variation of streamlines forβ 2(0,0.5,2,5)atΩ=3.5,A =1,β1=0.1,Q=1,Da=0.1,Ha=0.1.

3.5.Local wall shear stress

The variation of local shear stress at the upper wall of the curved channel is depicted in figures 18–21 with a variation of embedded physical parameters.The shape of its graph is the same as the shape of the peristaltic waves present at the upper wall of the channel.It is clear that the wavy graph is predicted because of the complex nature of the peristaltic waves.The components of the first-and second-order slip parameters are displayed in figures 18 and 19 for a viscous fluid.In the absence of slip parameters,the magnitude of local shear stress is maximum.When the length of the slip parameters is increased,then the magnitude of local shear stress is strongly affected.Because of this,the magnitude of the local shear stress decreases as the length of the slip parameters is increased.Physically,in both figures 18 and 19,the magnitude of the wall shear stress is sharply reduced by enhancing the strength of both slip parameters.The effects of magnetic force on the local shear stress under stronger porosity effects are depicted in figure 20 for the peristaltic flow of viscous fluid through a curved domain.In the absence of a magnetic parameter,the magnitude of local shear stress is larger.The magnitude of local shear stress is decreased by enhancing the magnetic effects.The magnitude of local shear stress is increased by increasing the Darcy number(figure 21).These results are in contrast to the magnetic effects.Additionally,in the presence of a smaller Darcy number,the effects of the porous medium are stronger,and due to this,the magnitude of local shear stress is reduced.Physically,under larger magnetic and porosity effects,the magnitude of local shear stress is strongly affected because these two external forces overcome the curvature effects.Further,these two figures(figures 20 and 21)show that the larger strength of the magnetic field and porous medium play a dynamic role in disturbing the sinusoidal behavior of wall shear stress under curvature effects.Additionally,these two forces play a dynamic role in controlling the behavior of wall shear stress and the larger values of both physical parameters overcome the curvature effects.

3.6.Trapping phenomena

Another important rheological phenomena related to peristalsis is known as the trapping phenomena.This phenomenon is concerned with the bolus dynamics near the boundary walls.The shape of the bolus is the same as the shape of the peristaltic walls.The bolus is moved with the same speed and direction as the peristaltic waves.The influence of flow parameters on the trapping phenomena is visualized in figures 22–25.The effects of the magnetic parameter(Hartmann number)on the trapping phenomena under curvature,porosity and second-order slip effects are plotted in figure 22.We have noticed these results under larger porosity effects.In the absence of a magnetic field and the curved nature of the flow regime,the circulation of the bolus is present at the upper half of the channel under porosity effects.When the magnetic strength is increased by increasing the numeric values of the Hartmann number from 0 to 1,the circulation of the bolus is shifted to the lower half(figure 22(b)).As was our expectation,the circulation of the bolus disappears from both halves by increasing the Hartmann number from 5 to 50.Additionally,it is noticed that not only does the circulation of the bolus disappear,but the number of the streamlines is reduced too at a larger strength of the magnetic field.The effects of the porous medium are visualized in figure 23 under the curvature,second-order slip parameter and magnetic effects.Here,we have obtained the reverse behavior as compared with the results of figure 22.Under smaller porosity effects,the circulation of the bolus exists at the upper half of the regime at smaller porosity whenDa→∞.The circulation of the bolus is disturbed by enhancing the physical effects of the porous medium(see figures 23(c)and(d)).Physically,the larger strength of both external forces,the magnetic force and porous medium,on the trapping phenomena have sound effects to reduce the circulation of the bolus and number of streamlines form both halves of the channel.The effects of the second-order slip condition on the trapping phenomena are displayed in figures 24 and 25 under curvature,magnetic and porosity effects,respectively.According to the literature survey,the slip parameter plays a dynamic role in effecting the circulation of the bolus and streamlines.Here,the streamlines are present at the lower half of the channel due to the curvature and porosity effects.As the length of the slip parameters are increased,the number of the streamlines are reduced at both halves of the channel.Due to the complex nature of biomimetic waves,the wavy behavior in both the streamlines and trapping phenomena can be observed.Additionally,these results are obtained under larger curvature effects.

4.Conclusions

A mathematical study has been presented for the laminar and incompressible flow of Newtonian fluid through a curved domain with peristaltic waves at the boundary walls with a complex scenario.Due to the curved shape of flow geometry,curvilinear coordinates are used to govern the rheological equations.These equations have been non-dimensionalized by using similarity transformations.The integration technique has been used to solve the boundary value problem(BVP).The effects of curvature,magnetic,porosity and slip parameters on the axial velocity,pressure gradient,shear component of the stress tensor,pressure rise,local wall shear stress and trapping phenomena have been presented graphically and discussed in detail.The mathematical formulation which is relevant to the peristaltic pumping of a viscous fluid through the curved nature of flow geometry has shown that:

In axial velocity,the magnitude of the velocity profile is reduced near the boundary walls by increasing the length of the slip parameters.The magnitude of velocity is reduced for larger magnetic and porosity effects.The amplitude of the pressure gradient is reduced sharply as the slip parameters are enhanced under larger porosity effects.The larger strength of the magnetic field and porous medium reduced the magnitude of the pressure gradient for a viscous fluid.The magnitude of the pressure rise related to the peristaltic flow of Newtonian fluid through a curved domain is reduced linearly by increasing the slip effects and amplitude ratios of peristaltic waves.The larger strength of the magnetic field and porosity medium have a dynamic influence on reducing the pumping phenomena.The magnitude of the pressure gradient is reduced(increased)at larger values of the Hartmann number(Darcy number).The magnitude of the shearing stress tensor is reduced by increasing the length of the slip parameters.The boundary layer phenomena is obtained for larger magnetic and porosity effects,respectively.The magnitude of local wall shear stress is reduced by increasing the length of the slip parameters under porosity and curvature effects,respectively.Under the larger strength of magnetic and porous medium effects,the amplitude of local wall shear stress is reduced.The circulation of the bolus disappeared under larger magnitude of the magnetic and porosity forces,respectively,in the upper half of the channel.Under larger magnetic effects,a larger number of streamlines are present at the lower half of the channel due to the curved nature of the geometry.The circulation of the bolus appears at the upper half of the straight channel under larger magnetic effects.While no circulation of the bolus is present under larger porosity effects,the number of streamlines are larger at the lower half for the curved channel.The number of streamlines is reduced by increasing the length of both slip parameters.

In the current manuscript,the combined effects of magnetic field,porosity and second-order slip parameters on the peristaltic propulsion of viscous fluid through a curved regime are discussed.The main reason for performing this study is to understand the transportation of biological fluids through a curved channel.This study shows that the curved shape of the channel has dynamic effects on the transport phenomena of physiological fluids.Additionally,the role of the magnetic field on the moving fluid is more useful in the curative domain to control the motion of biological fluids through curved vessels.Furthermore,the porosity effects on transportation phenomena are utilized in chemical engineering.The role of complex sinusoidal waves on the boundary walls is utilized to boost the pumping phenomena.Ultimately,the second-order slip has sound impacts on the rheological features of bio fluids.This investigation can be extended to non-Newtonian liquids.This investigation gives information about the movement of biological fluids through complex peristaltic pumps under a larger porosity medium for engineers and scientists.Additionally,these effects are observed in the presence of the wavy nature of peristaltic propulsion.These are the novelties of the current study.This rheological analysis has various applications in medical,chemical and petroleum engineering domains for the manufacturing of micro-scale devices that are useful in medical therapy,blood pumping,drug delivery systems,chemical cooling systems and geophysics,respectively.Some recent contributions in this domain are cited in[41–48].

Conflict of interest

The authors declare no conflict of interest.