Next Article in Journal
Matrix Factorization Techniques in Machine Learning, Signal Processing, and Statistics
Next Article in Special Issue
Optimal Dirichlet Boundary Control for the Corotational Oldroyd Model
Previous Article in Journal
On the Structure of the Mislin Genus of a Pullback
Previous Article in Special Issue
Analysis of Control Problems for Stationary Magnetohydrodynamics Equations under the Mixed Boundary Conditions for a Magnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid

1
Department of Basic Sciences, Preparatory Year Deanship, King Faisal University, Al-Ahsa 31982, Saudi Arabia
2
Department of Computer Science, National University of Sciences and Technology (NUST), Balochistan Campus (NBC), Quetta 87300, Pakistan
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2673; https://doi.org/10.3390/math11122673
Submission received: 17 May 2023 / Revised: 9 June 2023 / Accepted: 9 June 2023 / Published: 12 June 2023
(This article belongs to the Special Issue Theoretical Research and Computational Applications in Fluid Dynamics)

Abstract

:
Nanofluids are extremely beneficial to scientists because of their excellent heat transfer rates, which have numerous medical and industrial applications. The current study deals with the peristaltic flow of nanofluid (i.e., Casson nanofluid) in a symmetric elastic/compliant channel. Buongiorno’s framework of nanofluids was utilized to create the equations for flow and thermal/mass transfer along with the features of Brownian motion and thermophoresis. Slip conditions were applied to the compliant channel walls. The thermal field incorporated the attributes of viscous dissipation, ohmic heating, and thermal radiation. First-order chemical-reaction impacts were inserted in the mass transport. The influences of the Hall current and mixed convection were also presented within the momentum equations. Lubricant approximations were exploited to make the system of equations more simplified for the proposed framework. The solution of a nonlinear system of ODEs was accomplished via a numerical method. The influence of pertinent variables was examined by constructing graphs of fluid velocity, temperature profile, and rate of heat transfer. The concentration field was scrutinized via table. The velocity of the fluid declined with the increment of the Hartman number. The effects of thermal radiation and thermal Grashof number on temperature showed opposite behavior. Heat transfer rate was improved by raising the Casson fluid parameter and the Brownian motion parameter.

1. Introduction

Heat transfer research is extremely useful in a wide range of engineering disciplines, such as evaporative cooling thrust bearing, drag reduction, and circular heat exchanger prototypes. Fluids are frequently used as thermal carriers in applications such as automotive cooling and heating and are a central role in the industrial process. Heat transfer enhancement can enhance the effectiveness of a thermoelectric generator. Heat transfer improvement can result in more compact cooling equipment, conserving fuel, space, and wealth. Comparing nanofluids to their base fluids reveals that they are more effective at transmitting heat and have better thermal conductivity. Consequently, it is widely acknowledged and appreciated from both a philosophical and empirical standpoint that distributing nanoparticles in a biofluid may enhance its thermophysical qualities. Nanofluids reduce energy use, improve thermal efficiency, speed up operations, and lengthen the useful life of the equipment. It was Choi [1] who initially suggested using nanofluids. Once these fluids were accepted by the scientific world, they immediately assumed a key function. Chon et al. [2] securitized how temperature and particle size are significant for enhancing the thermal conductivity of the Al2O3 nanomaterial. Thermophoresis and random motion spreading were named by Buongiorno [3] as the most crucial factors in enhancing common liquids’ capacity to transmit heat. Tripathi and Bég [4] elaborated on the peristalsis of nanofluid in a drug delivery system in biomedical science. The peristaltic motion of a nanoliquid with wall and slip characteristics was explored by Hayat et al. [5]. The hydrothermal characteristics of magnetohydrodynamic (MHD) nanofluid flow and entropy formation in a trapezoidal crevice were examined by Atashafrooz et al. [6]. The peristaltic motion of Carreau–Yasuda nanomaterial with entropy was stated by Ahmed et al. [7]. The peristaltic activity of nanofluid with the Hall effect and mixed convection was considered by Alsaedi et al. [8]. Abbasi et al. [9] analyzed the peristaltic flow of nanofluid, considering the aspects of variable thermal conductivity and entropy optimization. Considering the peristaltic motion of nanofluid having activation energy and gyrotactic motile microbial organisms, Akbar et al. [10] conducted a bioconvection investigation. Hina et al. [11] modeled the peristalsis of the Carreau–Yasuda nanomaterial by incorporating electro-osmotic transport. A few related studies are cited through [12,13,14,15,16,17,18].
The contraction and expansion of an expansible tube or channel characterizes the pumping movement known as peristalsis. This is the movement of material in the wave’s propagation direction. Many technical, industrial, and biological applications make use of this feature. Several scientists are interested in peristaltic nanofluid flow due to its wide array of potential uses. The peristaltic deploying process has been the subject of several theoretical and empirical research studies. These investigations include ovum migration in the fallopian tube, capillary roller pumps, chyme movement in the intestine, and flow from the kidney to the bladder. Initially, Latham [19] examined the peristaltic activity inside a pump. A peristaltic framework with a low Reynolds number and long wavelength was described by Shapiro et al. [20]. Shugan and Smirnov [21] discussed the peristaltic flow with mass transfer in a channel, considering their understanding of wall vibrations. Ellahi et al. [22] reported the aspects of mass and thermal transport with peristalsis in a duct of non-uniform geometry. The consequences of a magnetic field for peristaltic activity under the influence of a double electrical layer were discussed by Tripathi et al. [23]. Aspects of radiative heat (called thermal radiation, in other words) and activation energy for the peristaltic movement of Eyring–Powell nanoliquid were scrutinized by Nisar et al. [24]. Peristalsis with a Maxwell fluid pattern incorporating convective properties was studied by Iqbal et al. [25]. Khushi and Abbasi [26] examined the peristaltic motion of TiO2–Ag/EG (hybrid nanofluids) with the Hall current. Nisar et al. [27] studied aspects of the peristalsis of a non-Newtonian nanomaterial. Javed et al. [28] modeled the peristaltic activity of a wavy micro-channel. Yasin et al. [29] inquired about the numerical inquiry for the peristaltic mechanism of Eyring–Powell fluid with Joule heating and slip aspects. Zhang et al. [30] explored the wall characteristics of the peristalsis of a bionic reactor. The compliant walls and channel also have an impact on the sinusoidal wave’s shape. This led various studies to look into how wall characteristics affect heat-exchange peristaltic transport. Few studies in this direction have been seen through [31,32,33,34].
Due to numerous uses in sectors, including petroleum drilling, polymer production, and others, non-Newtonian fluids are important. Various materials, such as blood, caramel, sauces/ketchup, honey, blood, jellies, porridge, etc., exhibit this model’s shear-thinning tendency. The most popular type of model for non-Newtonian properties is the Casson fluid, which represents the polymeric material. Yield stress is demonstrated by this model. When the shear stress is lower than the yield stress, Casson liquid acts like a solid. After yield stress exceeds shear stress, it begins to flow like a melted substance. Casson [35] first proposed the rheological model for the Casson fluid. Divya et al. [36] reported the aspects of the peristaltic flow of Casson fluid with thermal properties. Abbas et al. [37] analyzed the peristalsis of Casson liquid with wall characteristics. Priam et al. [38] accounted for the numerical inquiry of the peristaltic flow of Casson fluid. Hafez et al. [39] assessed the mass and thermal transport with the peristaltic flow of Casson material via the inclined channel.
The features of a chemical reaction from a peristaltic flow are important in biochemical engineering. The process of mass transfer is impacted by the concentration disparity among species that are chemically interacting. Chemical species travel from a highly concentrated location to a less concentrated one in certain circumstances. The manufacturing of food, the initiation and diffusion of fog, the creation of ceramics and polymers, the freezing of crops to cause damage, the hydrometallurgical industry, geothermal reservoirs, nuclear reactor cooling, and the recovery of thermal oil are all examples of practical uses for chemical reactions. A few findings in this field are quoted in the refs [40,41,42,43,44,45,46].
In this study, we examined the Casson fluid flow with the peristaltic phenomenon in a symmetric channel. A first-order chemical reaction was taken in this study. Conditions of slip were applied to the compliant/elastic channel. The Buongiorno version was used, and it included the novel characteristics of thermophoretic and Brownian motion. Several innovative features set apart the current study. It may have useful applications in disciplines such as biomedical engineering, microfluidics, and heat transfer and has the potential to advance our understanding of fluid dynamics in complicated systems. In addition, viscous dissipations are an important aspect of engineering analysis and design, enabling engineers to optimize fluid flow, heat transfer, and energy efficiency in a wide range of applications, such as turbomachinery, aerodynamics, and mighty planets. To solve the resulting problem, numerical techniques were employed. Graphs were plotted for the velocity and thermal fields.

2. Statement

We consider the Casson nanofluid movement in a 2 D elastic channel with a width of 2 d 1 . The channel walls are subjected to the magnetic field ( B 0 ) in a perpendicular direction. With low magnetic Reynolds numbers, the induced magnetic field disappears. Ohmic heating, thermal radiation, and viscous dissipation are present in the thermal equation. Further, slip conditions are applied on the channel walls (see Figure 1). The governing expressions of an incompressible Casson nanofluid in terms of Cartesian coordinates are given by:
y = ± η ( x , t ) = ± [ d 1 + a sin 2 π λ ( x c t ) ] ,
where t and c are the time and speed of the wave, respectively. The current flow problem’s rheological state equation is stated as
τ i j = { 2 ( μ β + p y / 2 π ) e i j ,             π > π c 2 ( μ β + p y / 2 π c ) e i j ,             π < π c    
where p y is the fluid yield stress, μ β is the dynamic viscosity, π = e i j e i j , and e i j represents the ( i , j ) t h factors of the rate of deformation while π c implies the significant value of this product-based non-Newtonian material. The Lorentz force (denoted by F) is calculated using Ohms law:
F = ( σ B 0 2 ( 1 + m 2 ) ( u m v ) , σ B 0 2 ( 1 + m 2 ) ( v + m u ) , 0 )
Below are the equations that govern the particular flow problem.
u x + v y = 0 ,
ρ f ( u t + u u x + v u y ) = p x + [ μ { 1 + 1 β } { 2 u x 2 + 2 u y 2 } ] σ B 0 2 ( 1 + m 2 ) ( u m v ) + g ρ f β T ( T T 0 ) + g ρ f β C ( C C 0 ) ,
ρ f ( v t + u v x + v v y ) = p y + μ ( 1 + 1 β ) ( 2 v x 2 + 2 v y 2 ) σ B 0 2 ( 1 + m 2 ) ( v + m u ) ,
ρ f c f ( T t + u T x + v T y ) = k ( 2 T x 2 + 2 T y 2 ) + [ ( u y + v x ) 2 + 4 ( u x ) 2 ] μ ( 1 + 1 β ) + [ D T T m { ( T y ) 2 + ( T x ) 2 } + D B { C y T y + C x T x } ] ρ p c p q r y + σ B 0 2 ( 1 + m 2 ) ( u 2 + v 2 ) ,
C t + u C x + v C y = ( 2 C x 2 + 2 C y 2 ) D B + ( 2 T x 2 + 2 T y 2 ) D T T m k 1 ( C C 0 ) .
The corresponding boundary conditions are displayed as
u ± β 1 ( 1 + 1 β ) ( u x + u y ) = 0 ,   T ± β 2 T y = { T 1 T 0 } , C ± β 3 C y = { C 1 C 0 }   at   y = ± η ,
[ τ 1 3 x 3 + m 1 3 x t 2 + d 2 t x ] η = ρ f ( u t + u u x + v u y ) + μ ( 1 + 1 β ) ( 2 u x 2 + 2 u y 2 ) σ B 0 2 ( 1 + m 2 ) ( u m v ) + g ρ f β T ( T T 0 ) + g ρ f β C ( C C 0 ) , at   y = ± η .
Here, ( u , v ) represent the elements of liquid velocity in ( x , y ) directions, ρ f is the density of nanoliquid, p is the pressure, β is the Casson fluid variable, ν is the kinematic viscosity, σ is the electric conductivity of liquid, and k is the thermal conductivity, Furthermore, τ 1 is the elastic tension, D B is the Brownian movement coefficient, D T is the thermophoretic diffusion coefficient, m 1 is the area per mass unit, T m is the mean temperature, τ 1 is the elastic tension, d is the damping coefficient, and ( T 1 , T 0 ) represents the fluid temperature, whose concentration is denoted by ( C 1 , C 0 ) on the upper and lower channel walls. The mathematical illustration of radiative heat flux is demarcated as [26]:
q r = 16 σ ¯ T 0 3 3 k ¯ T y ,
Now we define stream function u = ψ y and v = δ ψ x . Non-dimensional variables are defined as
u * = u c ,   v * = v c   ,   x * = x λ ,   y * = y d 1 ,   t * = c t λ , β i * = β i d 1   ( i = 1 3 ) , η * = η d 1 ,   p * = d 1 2 p c λ μ ,   ϕ = C C 0 C 1 C 0 , θ = T T 0 T 1 T 0 ,
using Equations (4)–(10). The resulting problem later utilizes ( δ 0 ) , which is the long wavelength. In addition, ( Re 0 ) , which describes small Reynolds number suppositions, is given by
( 1 + 1 β ) 4 ψ y 4 ( M 2 1 + m 2 ) 2 ψ y 2 + G r θ y + G c ϕ y = 0 ,
( 1 + P r R n ) 2 θ y 2 + P r N b θ y ϕ y + P r N t ( θ y ) 2 + B r [ ( 1 + 1 β ) ( 2 ψ y 2 ) 2 + M 2 ( ψ y ) 2 ] = 0 ,
N t 2 θ y 2 + N b 2 ϕ y 2 N b S c ζ ϕ = 0 .
Boundary conditions are expressed as
ψ y ± β 1 [ ( 1 + 1 β ) ( 2 ψ y 2 ) ] = 0   at   y = ± η ,
[ E 1 3 x 3 + E 2 3 x t 2 + E 3 2 x t ] η = ( 1 + 1 β ) 3 ψ y 3 M 2 1 + m 2 ψ y + G r θ + G c ϕ   at   y = ± η ,
θ ± β 2 θ y = { 1 0 } ,   ϕ ± β 3 ϕ y = { 1 0 }   at   y = ± η .
The continuation of Equation (3) is automatically satisfied where dimensionless parameters are defined by
δ = d 1 λ , ε = a d 1 , G c = g ρ f β C ( C C 0 ) d 1 2 μ c ,   E c = c 2 c f ( T 1 T 0 ) , P r = ν α ,   R e = c ρ d 1 μ , M = σ μ B 0 d 1 ,   S c = ν D B ,   N t = τ D T ( T 1 T 0 ) T m ν ,   N b = τ D B ( C 1 C 0 ) ν ,   G r = g ρ f β T ( T T 0 ) d 1 2 μ c , β = μ β 2 π c p y ,   B r = P r E c , E 1 = τ d 1 3 λ 3 μ c , E 2 = m 1 c d 1 3 λ 3 μ , E 3 = d d 1 3 λ 2 μ , ζ = k 1 d 2 1 ν .

3. Numerical Results and Analysis

The system of Equations (13)–(15) with subjective boundary conditions (16)–(18) was numerically solved via the in-built algorithm NDSolve [5,7,8,14,24,27,29,44] in MATHEMATICA. This scheme is very advanced and useful for finding solutions to differential equations. It handles a wide range of differential equations, including ordinary differential equations (ODEs) and boundary value problems. It employs high-precision numerical algorithms to obtain accurate solutions. It uses adaptive step-size control, error estimation, and sophisticated numerical techniques to ensure reliable results. This segment was made to analyze the effects of several embedded variables on the thermal field, velocity, heat transfer rate, and concentration.

3.1. Velocity

Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8 were plotted to see aspects of pertinent variables for the velocity profile. Figure 2 describes the temperature curves for the Casson fluid variable β . The velocity of the fluid was enhanced via β .  Figure 3 designates the consequence of velocity slip parameter β 1 for velocity. Enhancing behavior was observed against velocity for velocity slip parameter β 1 . This is because fluid particles which have the slip velocity in direct contact with the surface can successfully transmit some of their momentum to the fluid particles close to them. The fluid velocity increases as a result. The influence of Hall parameter m is displayed in Figure 4. In this particular graph, we noted that fluid velocity was enhanced via Hall current variable m . The significance of the thermal Grashof number G r for the fluid velocity profile is shown in Figure 5. The velocity of the liquid was enhanced via a larger thermal Grashof number G r . When the thermal Grashof parameter was high, the buoyant forces in a fluid flow dominated the viscous forces. Impressions of concentration Grashof number G c are exhibited in Figure 6 . It indicates that higher values of the mass Grashof number G c give rise to an enhancement in velocity. The impacts of E 1 , E 2 , and E 3 on the velocity profile are portrayed in Figure 7. The velocity of the fluid heightened for larger values of the viscous damping force and stiffness while decreasing for rigidity variables. The behavior of the velocity profile in relation to the Hartmann number M is seen in Figure 8. Near the middle of the channel, a clear decline in the velocity profile was shown along with an enhancement in M . Subsequently, the imposed magnetic field created a potent opposite force (the Lorentz force), and the fluid could not move. Additionally, with high Hartmann numbers, the utmost velocity shifted to the lower boundary, causing the fluid velocity to lose its symmetry close to the conduit’s midpoint.

3.2. Temperature

Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15 and Figure 16 depict how different variables affect the thermal field. The graphical representation of velocity for a few rising values of Casson liquid parameter β is shown via Figure 9. The acclivitous β shows a growth in the temperature. The influence of the temperature slip variable β 2 is displayed via Figure 10. As we can see in this figure, we observed an enhancement in the thermal profile of the fluid. Figure 11 depicts the consequence of the Brinkman variable B r for temperature. The temperature of the liquid increased with larger values of B r . This is because B r caused the flow to form a resistance due to its shear, which caused increased generation of heat due to the effects of viscous dissolution. As a result, the fluid’s temperature rose. Effects of thermal Grashof number G r are presented via Figure 12 . From this curve, we noticed that the temperature was enhanced. As the thermal Grashof parameter increased, the buoyancy forces became stronger, resulting in faster fluid motion and increased temperature. Figure 13 demonstrates the effect of Harman parameter M variation on θ . The graph clearly indicates that a rise in Harman number M produces an increase in temperature. Figure 14 illustrates the impressions of the radiation parameter R n on the thermal field. The variable θ declines with radiation parameter R n . The combined consequences of N t and N b for the thermal field are displayed via Figure 15. We observed that the temperature of the fluid was enhanced by increasing both parameters simultaneously. When the Brownian motion rose, so did the kinetic energy of the fluid molecules. As a result, the frequency and strength of collisions between fluid molecules and suspended particles increased. This raised the temperature. The wall parameters E 1 , E 2 , and E 3 are displayed in Figure 16 to show how they affect temperature. We observed that the fluid thermal profile was an ever-increasing function of E 1 and E 2 , and it decayed for E 3 .

3.3. Heat Transfer Rate

Figure 17, Figure 18, Figure 19, Figure 20, Figure 21, Figure 22 and Figure 23 are designed to see the behavior of Brownian motion variable N b , Casson fluid variable β , thermal slip parameter β 2 , Brinkman number B r , radiation parameter R n , Hall current variable m , and thermal Grashof number G r with regards to the rate of heat transfer at upper wall θ ( η ) . When the thermal Grashof parameter is high, it indicates that the buoyancy forces dominate over the viscous forces, leading to enhanced fluid motion and heat transfer. Effects of Casson fluid variable β are represented via Figure 17 . An increasing trend against Casson fluid variable β   was   witnessed . Impacts of thermal slip parameter β 2 are represented via Figure 18 . In this figure, we noticed that the heat transfer rate increases against thermal slip parameter β 2 .  Figure 19 displays the features of Brinkman number B r on θ ( η ) . An increasing trend was noticed for θ ( η ) via B r .  Figure 20 designates the influence of thermal Grashof number G r against heat transfer rate θ ( η ) . An increase in the Brinkman number implies a greater dominance of viscous dissipation over conductive heat transfer, leading to an overall higher heat transfer rate. It can be seen that larger levels of G r increased the rate of heat transfer of nanoliquid acclivities. The behavior of Brownian variable N b is exhibited in Figure 21. The finding of this graph shows that heat transfer rate was boosted. Figure 22 describes the aspects of the Hall current variable for θ ( η ) . A decreasing trend was detected for higher values of the Hall current parameter. Moreover, the consequences of radiation parameter R n for θ ( η ) are shown in Figure 23 . The heat transfer rate of the fluid declined via R n .

3.4. Concentration Field

Numerical values for concentration field ϕ ( 0 ) are represented via Table 1 for different physical variables. Variations of thermophoresis N t can be seen in Table 1. As we can notice, the magnitude value of the concentration declined via thermophoresis variable N t . When the thermophoresis parameter increased, the particles migrated from regions of low temperature to regions of high temperature. In this case, the concentration of the fluid decreased. Impacts of concentration Grashof number G c on concentration field ϕ ( 0 ) are seen in Table 1 . Higher values of concentration Grashof number G c decreased the nanoparticle concentration. In fact, as this parameter increased, the buoyancy force enhanced, and the fluid motion became more turbulent. Turbulence caused the fluid to mix and disperse the concentration more equally, minimizing concentration gradients in the liquid. Impressions of Casson liquid variable β and mass slip parameter β 3 on concentration field ϕ ( 0 ) are presented via Table 1. A decreasing trend was noticed when increasing the values of these parameters. Larger values of chemical reaction variable ζ gave a fall to the concentration field ϕ ( 0 ) (see Table 1). As the reaction progressed, the concentration of the reactants decreased, leading to a decrease in the overall concentration of the fluid. Aspects of N b for concentration field ϕ ( 0 ) can be witnessed in Table 1. In this, we can observe that the concentration field ϕ ( 0 ) intensified because Brownian motion influenced the dispersion and stability of nanoparticles within the fluid.

4. Results Validation

We produced Table 2 to create a validation check for our numerical result. This comparison makes it clear that the current findings show excellent agreement with [5] for Newtonian fluid and G r = G c = R n = m = ζ = 0 .

5. Conclusions

The aim of the current analysis was to see the effects of mixed convection and Hall current on the peristaltic flow of a Casson nanomaterial. Slip conditions were imposed on an elastic channel. Further, the impacts of viscous dissipation and thermal radiation were also present in thermal transport. The nonlinear system of equations was tackled numerically. The key findings of the study are listed below.
  • The velocity of the flowing fluid increased with velocity slip β 1 and Casson fluid parameters β .
  • By enhancing the thermal Grashof number G r , velocity was enhanced.
  • Temperature raised for larger values of nanofluid parameters known as thermophoresis N t and Brownian motion N b .
  • The opposite trend was noticed for radiation parameter R n and Brinkman number B r against temperature.
  • The rate of heat transfer was enhanced for Casson fluid β and thermal Grashof number G r .
  • The role of Hall current number m and thermal slip parameter β 2 in heat transfer rate was inverse.
  • Concentration declined via chemical reaction ζ and thermophoresis N t .
  • Temperature and concentration heightened for elastic parameters E 1 and E 2 and lessened for E 3 .

Author Contributions

Conceptualization, H.Y. and Z.N.; methodology, Z.N.; software, Z.N.; validation, H.Y. and Z.N.; formal analysis, H.Y.; investigation, Z.N.; resources, H.Y.; writing—original draft preparation, Z.N.; writing—review and editing, H.Y.; visualization, H.Y.; funding acquisition, H.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Deanship of Scientific Research, the Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 3609).

Data Availability Statement

All data are available within the manuscript.

Acknowledgments

This work was supported by the Deanship of Scientific Research, the Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia (Grant No. 3609).

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

( x , y ) Cartesian coordinates g gravitational acceleration
( u , v ) velocity components d coefficient of viscous damping
a wave amplitude τ 1 elastic tension
t time m 1 mass per unit area
c wave speed R n radiation parameter
d 1 damping half channel width β Casson fluid parameter
ε amplitude ratio β 1 , β 2 , β 3   Slip parameters
p pressure ψ stream function
ρ f Nanofluid density T 1 Fluid temperature at the upper wall
α thermal diffusivity T 0 Fluid temperature at the lower wall
k thermal conductivity C 1 Fluid concentration at the upper wall
λ wavelength C 0 Fluid concentration at the lower wall
σ electrical conductivity T m mean temperature
B 0 applied magnetic field Re Reynolds number
ν kinematic viscosity δ wave number
β T thermal expansion coefficients Ec Eckert number
β C concentration expansion coefficients P r Prandtl number
D B Brownian motion coefficient M Hartman number
D T thermophoretic diffusion coefficient G r thermal Grashof number
N b Brownian motion parameter G c concentration Grashof number
E 1 , E 2 , E 3 wall parameters ζ e t a     p a r a m e t e r .   s   a g i n a s t   h o r e s i s   p a r a m e t e r .   Chemical reaction parameter
N t thermophoresis parameter θ fluid temperature
ϕ fluid concentration

References

  1. Choi, S.U.S. Enhancing thermal conductivity of fluids with nanoparticles. In Developments and Applications of Non-Newtonian Flows; Siginer, D.A., Wang, H.P., Eds.; Argonne National Lab. (ANL): Argonne, IL, USA, 1995; Volume 231, pp. 99–105. [Google Scholar]
  2. Chon, C.H.; Kihm, K.D.; Lee, S.P.; Choi, S.U.S. Empirical correlation finding the role of temperature and particle size for nanofluid (Al2O3) thermal conductivity enhancement. Appl. Phys. Lett. 2005, 87, 153107. [Google Scholar] [CrossRef]
  3. Buongiorno, J. Convective transport in nanofluids. ASME J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  4. Tripathi, D.; Bég, O.A. A study on peristaltic flow of nanofluids: Application in drug delivery systems. Int. J. Heat Mass Transf. 2014, 70, 61–70. [Google Scholar] [CrossRef]
  5. Hayat, T.; Nisar, Z.; Ahmad, B.; Yasmin, H. Simultaneous effects of slip and wall properties on MHD peristaltic motion of nanofluid with Joule heating. J. Magn. Magn. Mater. 2015, 395, 48–58. [Google Scholar] [CrossRef]
  6. Atashafrooz, M.; Sajjadi, H.; Delouei, A.A. Interacting influences of Lorentz force and bleeding on the hydrothermal behaviors of nanofluid flow in a trapezoidal recess with the second law of thermodynamics analysis. Int. Commun. Heat Mass Transf. 2020, 110, 104411. [Google Scholar] [CrossRef]
  7. Ahmed, B.; Hayat, T.; Alsaedi, A.; Abbasi, F.M. Entropy generation analysis for peristaltic motion of Carreau–Yasuda nanomaterial. Phys. Scr. 2020, 95, 055804. [Google Scholar] [CrossRef]
  8. Alsaedi, A.; Nisar, Z.; Hayat, T.; Ahmad, B. Analysis of mixed convection and hall current for MHD peristaltic transport of nanofluid with compliant wall. Int. Commun. Heat Mass Transf. 2021, 121, 105121. [Google Scholar] [CrossRef]
  9. Abbasi, F.M.; Gul, M.; Shanakhat, I.; Anjum, H.J.; Shehzad, S.A. Entropy generation analysis for magnetized peristaltic movement of nanofluid through a non-uniform asymmetric channel with variable thermal conductivity. Chin. J. Phys. 2022, 78, 111–131. [Google Scholar] [CrossRef]
  10. Akbar, Y.; Alotaibi, H.; Iqbal, J.; Nisar, K.S.; Alharbi, K.A.M. Thermodynamic analysis for bioconvection peristaltic transport of nanofluid with gyrotactic motile microorganisms and Arrhenius activation energy. Case Stud. Therm. Eng. 2022, 34, 102055. [Google Scholar]
  11. Hina, S.; Kayani, S.M.; Mustafa, M. Aiding or opposing electro-osmotic flow of Carreau—Yasuda nanofluid induced by peristaltic waves using Buongiorno model. Waves Random Complex Media 2022, 1–17. [Google Scholar] [CrossRef]
  12. Muhammad, K.; Abdelmohsen, S.A.; Abdelbacki, A.M.; Aziz, A. Cattaneo-Christov (C–C) heat flux in Darcy-Forchheimer (DF) flow of fourth-grade nanomaterial with convective heat and mass conditions. Case Stud. Therm. Eng. 2022, 3, 102152. [Google Scholar] [CrossRef]
  13. Atashafrooz, M.; Sajjadi, H.; Delouei, A.A. Simulation of combined convective-radiative heat transfer of hybrid nanofluid flow inside an open trapezoidal enclosure considering the magnetic force impacts. J. Magn. Magn. Mater. 2023, 567, 170354. [Google Scholar] [CrossRef]
  14. Nisar, Z.; Yasmin, H. Analysis of motile gyrotactic micro-organisms for the bioconvection peristaltic flow of Carreau-Yasuda bionanomaterials. Coatings 2023, 13, 314. [Google Scholar] [CrossRef]
  15. Yasmin, H.; Giwa, S.O.; Noor, S.; Aybar, H.Ş. Reproduction of Nanofluid Synthesis, Thermal Properties and Experiments in Engineering: A Research Paradigm Shift. Energies 2023, 16, 1145. [Google Scholar] [CrossRef]
  16. Akram, S.; Athar, M.; Saeed, K.; Razia, A.; Muhammad, T.; Hussain, A. Hybrid double-diffusivity convection and induced magnetic field effects on peristaltic waves of Oldroyd 4-constant nanofluids in non-uniform channel. Alex. Eng. J. 2023, 65, 785–796. [Google Scholar] [CrossRef]
  17. Yasmin, H.; Giwa, S.O.; Noor, S.; Sharifpur, M. Thermal Conductivity Enhancement of Metal Oxide Nanofluids: A Critical Review. Nanomaterials 2023, 13, 597. [Google Scholar] [CrossRef]
  18. Iqbal, J.; Abbasi, F.M.; Alkinidri, M.; Alahmadi, H. Heat and mass transfer analysis for MHD bioconvection peristaltic motion of Powell-Eyring nanofluid with variable thermal characteristics. Case Stud. Therm. Eng. 2023, 43, 102692. [Google Scholar] [CrossRef]
  19. Latham, T.W. Fluid Motion in a Peristaltic Pump. Master’s Thesis, MIT, Cambridge, MA, USA, 1966. [Google Scholar]
  20. Shapiro, A.H.; Jaffrin, M.Y.; Weinberg, S.L. Peristaltic pumping with long wavelength at low Reynolds number. J. Fluid Mech. 1969, 37, 799–825. [Google Scholar] [CrossRef]
  21. Shugan, I.V.; Smirnov, N.N. Peristaltic mass transfer in a channel under standing walls vibrations. Phys. Syst. Vib. 2001, 9, 71–78. [Google Scholar]
  22. Ellahi, R.; Bhatti, M.M.; Vafai, K. Effects of heat and mass transfer on peristaltic flow in a non-uniform rectangular duct. Int. J. Heat Mass Transf. 2014, 71, 706–719. [Google Scholar] [CrossRef]
  23. Tripathi, D.; Bhushan, S.; Bég, O.A. Transverse magnetic field driven modification in unsteady peristaltic transport with electrical double layer effects. Colloids Surf. A Physicochem. Eng. Asp. 2016, 506, 32–39. [Google Scholar] [CrossRef]
  24. 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. 2020, 116, 104655. [Google Scholar] [CrossRef]
  25. Iqbal, N.; Yasmin, H.; Bibi, A.; Attiya, A.A. Peristaltic motion of Maxwell fluid subject to convective heat and mass conditions. Ain Shams Eng. J. 2021, 12, 3121–3131. [Google Scholar] [CrossRef]
  26. Khushi, S.; Abbasi, F.M. Hall current and Joule heating effects on peristalsis of TiO2–Ag/EG hybrid nanofluids via a curved channel with heat transfer. Waves Random Complex Media 2022, 1–24. [Google Scholar] [CrossRef]
  27. Nisar, Z.; Hayat, T.; Alsaedi, A.; Ahmad, B. Mathematical modeling for peristalsis of couple stress nanofluid. Math. Meth. Appl. Sci. 2022, 1–19. [Google Scholar] [CrossRef]
  28. Javed, M.; Aslam, R.; Ibrahim, N. Peristaltic mechanism in a micro wavy channel. Therm. Sci. Eng. Prog. 2023, 38, 101530. [Google Scholar] [CrossRef]
  29. Yasin, M.; Hina, S.; Naz, R.; Abdeljawad, T.; Sohail, M. Numerical Examination on Impact of Hall Current on Peristaltic Flow of Eyring-Powell Fluid under Ohmic-Thermal Effect with Slip Conditions. Curr. Nanosci. 2023, 19, 49–62. [Google Scholar]
  30. Zhang, S.; Liang, W.; Li, C.; Wu, P.; Chen, X.D.; Dai, B.; Deng, R.; Lei, Z. Study on the effect of wall structures and peristalsis of bionic reactor on mixing. Chem. Eng. Sci. 2023, 267, 118373. [Google Scholar] [CrossRef]
  31. Shugan, I.V.; Smirnov, N.N.; Legros, J.C. Streaming flows in a channel with elastic walls. Phys. Fluids 2002, 14, 3502–3511. [Google Scholar] [CrossRef]
  32. Nisar, Z.; Hayat, T.; Alsaedi, A.; Ahmad, B. Wall properties and convective conditions in MHD radiative peristalsis flow of Eyring–Powell nanofluid. J. Therm. Anal. Calorim. 2021, 144, 1199–1208. [Google Scholar] [CrossRef]
  33. Choudhari, R.; Makinde, O.D.; Mebarek-Oudina, F.; Vaidya, H.; Prasad, K.V.; Devaki, P. Analysis of third-grade liquid under the influence of wall slip and variable fluid properties in an inclined peristaltic channel. Heat Transf. 2022, 51, 6528–6547. [Google Scholar] [CrossRef]
  34. Abd-Alla, A.M.; Abo-Dahab, S.M.; Thabet, E.N.; Bayones, F.S.; Abdelhafez, M.A. Heat and mass transfer in a peristaltic rotating frame Jeffrey fluid via porous medium with chemical reaction and wall properties. Alex. Eng. J. 2023, 66, 405–420. [Google Scholar] [CrossRef]
  35. Casson, N. A flow equation for the pigment oil suspension of the printing ink type. Rheology Disper. Syst. 1959, 84–102. [Google Scholar]
  36. Divya, B.B.; Manjunatha, G.; Rajashekhar, C.; Vaidya, H.; Prasad, K.V. Analysis of temperature dependent properties of a peristaltic MHD flow in a non-uniform channel: A Casson fluid model. Ain Shams Eng. J. 2021, 12, 2181–2191. [Google Scholar] [CrossRef]
  37. Abbas, Z.; Rafiq, M.Y.; Hasnain, J.; Javed, T. Peristaltic transport of a Casson fluid in a non-uniform inclined tube with Rosseland approximation and wall properties. Arab. J. Sci. Eng. 2021, 46, 1997–2007. [Google Scholar] [CrossRef]
  38. Priam, S.S.; Nasrin, R. Numerical appraisal of time-dependent peristaltic duct flow using Casson fluid. Int. J. Mech. Sci. 2022, 233, 107676. [Google Scholar] [CrossRef]
  39. Hafez, N.M.; Abd-Alla, A.M.; Metwaly, T.M.N. Influences of rotation and mass and heat transfer on MHD peristaltic transport of Casson fluid through inclined plane. Alex. Eng. J. 2023, 68, 665–692. [Google Scholar] [CrossRef]
  40. Akram, J.; Akbar, N.S.; Maraj, E. Chemical reaction and heat source/sink effect on magnetonano Prandtl-Eyring fluid peristaltic propulsion in an inclined symmetric channel. Chin. J. Phys. 2020, 65, 300–313. [Google Scholar] [CrossRef]
  41. Nisar, Z.; Hayat, T.; Alsaedi, A.; Momani, S. Peristaltic flow of chemically reactive Carreau-Yasuda nanofluid with modified Darcy’s expression. Mater. Today Commun. 2022, 33, 104532. [Google Scholar] [CrossRef]
  42. Vaidya, H.; Prasad, K.V.; Khan, M.I.; Oudina, F.M.; Tlili, I.; Rajashekhar, C.; Elattar, S.; Khan, M.I.; Al-Gamdi, S.G. Combined effects of chemical reaction and variable thermal conductivity on MHD peristaltic flow of Phan-Thien-Tanner liquid through inclined channel. Case Stud. Therm. Eng. 2022, 36, 102214. [Google Scholar] [CrossRef]
  43. Hayat, T.; Nisar, Z.; Alsaedi, A. Bioconvection and Hall current analysis for peristalsis of nanofluid. Int. Commun. Heat Mass Transf. 2021, 129, 105693. [Google Scholar] [CrossRef]
  44. Hussein, S.A.; Ahmed, S.E.; Arafa, A.A.M. Electrokinetic peristaltic bioconvective Jeffrey nanofluid flow with activation energy for binary chemical reaction, radiation and variable fluid properties. Z. Angew. Math. Mech. 2023, 103, e202200284. [Google Scholar] [CrossRef]
  45. Abbasi, A.; Khan, S.U.; Farooq, W.; Mughal, F.M.; Khan, M.I.; Prasannakumara, B.C.; El-Wakad, M.T.; Guedri, K.; Galal, A.M. Peristaltic flow of chemically reactive Ellis fluid through an asymmetric channel: Heat and mass transfer analysis. Ain Shams Eng. J. 2023, 14, 101832. [Google Scholar] [CrossRef]
  46. Saba, S.; Abbasi, F.M.; Shehzad, S.A. Magnetized peristaltic transportation of Boron-Nitride and Ethylene-Glycol nanofluid through a curved channel. Chem. Phys. Lett. 2022, 803, 139860. [Google Scholar] [CrossRef]
Figure 1. Schematic of the problem.
Figure 1. Schematic of the problem.
Mathematics 11 02673 g001
Figure 2. Consequences of β for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 2. Consequences of β for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g002
Figure 3. Consequences of β for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 3. Consequences of β for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g003
Figure 4. Consequences of m for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 4. Consequences of m for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g004
Figure 5. Consequences of G r on u where x = ϵ = 0.2 , t = 0.1 , M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 5. Consequences of G r on u where x = ϵ = 0.2 , t = 0.1 , M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g005
Figure 6. Consequences of G c for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 6. Consequences of G c for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g006
Figure 7. Consequences of E 1 , E 2 ,   and   E 3 for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 7. Consequences of E 1 , E 2 ,   and   E 3 for u where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g007
Figure 8. Consequences of M for u where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 8. Consequences of M for u where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g008
Figure 9. Consequences of β for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 , β 1 = β 2 = β 3 = 0.1 .
Figure 9. Consequences of β for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 , β 1 = β 2 = β 3 = 0.1 .
Mathematics 11 02673 g009
Figure 10. Consequences of β 2 for θ where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 3 = 0.1   and   β = 0.5 .
Figure 10. Consequences of β 2 for θ where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g010
Figure 11. Consequences of Br for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 11. Consequences of Br for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g011
Figure 12. Consequences of G r for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 12. Consequences of G r for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g012
Figure 13. Consequences of M for θ where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 13. Consequences of M for θ where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g013
Figure 14. Consequences of R n for θ where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 14. Consequences of R n for θ where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g014
Figure 15. Consequences of N b   and   N t for θ where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 15. Consequences of N b   and   N t for θ where x = ϵ = 0.2 , t = 0.1 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g015
Figure 16. Consequences of E 1 , E 2 ,   and   E 3 for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 16. Consequences of E 1 , E 2 ,   and   E 3 for θ where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 0.5 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g016
Figure 17. Consequences of β for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1 .
Figure 17. Consequences of β for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   x = ϵ = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1 .
Mathematics 11 02673 g017
Figure 18. Consequences of β 2 for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 3 = 0.1   and   β = 0.5 .
Figure 18. Consequences of β 2 for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g018
Figure 19. Consequences of B r for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 19. Consequences of B r for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g019
Figure 20. Consequences of G r for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 20. Consequences of G r for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g020
Figure 21. Consequences of N b for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 21. Consequences of N b for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g021
Figure 22. Consequences of m for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 22. Consequences of m for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g022
Figure 23. Consequences of R n for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Figure 23. Consequences of R n for θ ( η ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 2 ,   G c = 0.7 ,   N b = N t = 0.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = β 3 = 0.1   and   β = 0.5 .
Mathematics 11 02673 g023
Table 1. Effects of physical parameters on ϕ ( 0 ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 1.5 ,   N b = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = 0.1   a n d   B r = 1 .
Table 1. Effects of physical parameters on ϕ ( 0 ) where x = ϵ = 0.2 , t = 0.1 ,   M = 0.5 ,   m = 1 ,   G r = 1.5 ,   N b = 0.5 ,   R n = 1.5 ,   E 1 = 0.01 ,   E 2 = 0.01 ,   E 3 = 0.02 ,   β 1 = β 2 = 0.1   a n d   B r = 1 .
ParametersConcentration
N t G c β ζ β 3 N b ϕ ( 0 )
0.1 0.70.510.10.50.2484573
0.3 0.70.510.10.50.1952146
0.1 10.510.10.50.2464970
0.1 1.50.510.10.50.2430272
0.1 0.70.910.10.50.2365245
0.1 0.71.510.10.50.2226838
0.1 0.70.520.10.50.1589335
0.1 0.70.530.10.50.1100197
0.1 0.70.510.20.50.2265843
0.1 0.70.510.30.50.2076794
0.1 0.70.510.10.70.2152228
0.1 0.70.510.110.2209559
Table 2. Comparison of numerical results with ref. [5].
Table 2. Comparison of numerical results with ref. [5].
ParametersPresent WorkHayat et al. [5]
N b N t M β 2 β 3 θ y | y = η θ y | y = η
0.100.100.100.100.10−0.062447−0.062441
0.200.100.100.100.10−0.075485−0.075478
0.100.300.100.100.10−0.088061−0.088054
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yasmin, H.; Nisar, Z. Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid. Mathematics 2023, 11, 2673. https://doi.org/10.3390/math11122673

AMA Style

Yasmin H, Nisar Z. Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid. Mathematics. 2023; 11(12):2673. https://doi.org/10.3390/math11122673

Chicago/Turabian Style

Yasmin, Humaira, and Zahid Nisar. 2023. "Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid" Mathematics 11, no. 12: 2673. https://doi.org/10.3390/math11122673

APA Style

Yasmin, H., & Nisar, Z. (2023). Mathematical Analysis of Mixed Convective Peristaltic Flow for Chemically Reactive Casson Nanofluid. Mathematics, 11(12), 2673. https://doi.org/10.3390/math11122673

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop