Insightful Facts on Peristalsis Flow of Water Conveying Multi-Walled Carbon Nanoparticles Through Elliptical Ducts With Ciliated Walls

In this research, a mathematical model is disclosed that elucidates the peristaltic flow of carbon nanotubes in an elliptic duct with ciliated walls. This novel topic of nanofluid flow is addressed for an elliptic domain for the very first time. The practical applications of current analysis include the customization of the mechanical peristaltic pumps, artificial cilia and their role in flow control, drug delivery and prime biological applications etc. The dimensional mathematical problem is transformed into its non-dimensional form by utilizing appropriate transformations and dimensionless parameters. Exact mathematical solutions are computed over the elliptic domain for the partial differential equations appearing in this convection heat transfer problem. A thorough graphical assessment is performed to discuss the prime results. The graphical visualization of the flow in this elliptic duct is obtained by plotting streamlines. The viscous effects are playing a vital role in the heat enhancement as compared to the molecular conduction. Since the incrementing Brinkman number results in a declined conduction due to viscous dissipation that eventually results in an enhanced temperature profile. This research first time elucidates the impacts of nanofluid flow on the peristaltic pumping through an elliptic domain having ciliated walls. Considering water as base fluid with multi-wall Carbon nanotubes for this ciliated elliptic domain having sinusoidal boundaries.


INTRODUCTION
Peristalsis is concerned with fluid dynamics problems involving flow inside ducts with wavelike deformable walls. The flexible walls' sinusoidal movement eventually aids in the development of flow along the duct's axis. This fascinating subject has a wide range of functional, engineering, and medical applications, making it a hot research topic. Slurries, various heated mixtures, and chemicals are among its primary industrial applications; see Jaffrin and Shapiro [1]. These are some finest mathematical models furnished by different precise researchers that elucidate the peristaltic flow in distinct geometries like cylinders, Barton [2] had numerically examined the peristaltic flow inside a cylindrical conduit, Siddiqui and Schwarz [3] had conveyed the perturbation solution analysis on Peristaltic pumping of a non-Newtonian fluid, Nadeem et al. [4] had interpreted the Homotopy Analysis solutions on the peristaltic pumping of non-Newtonian Sisko fluid, curved channels, Nadeem and Maraj [5] had disclosed the Homotopy Perturbation solutions on the Peristaltic pumping of Hyperbolic Tangent non-Newtonian fluid thorough a curved domain, Rashid et al. [6] had also utilized the Homotopy Perturbation analysis to interpret the magnetohydrodynamics effects on the peristaltic pumping of Williamson fluid inside a curved domain, rectangular ducts, Nadeem and Akram [7] had provided exact solutions on Jeffrey fluid flow within a rectangular domain having sinusoidal wall fluctuations, Ellahi et al. [8] provided the Homotopy Perturbation analysis on the three dimensional flow inside a rectangular conduit working as a peristaltic pump, Akram and Saleem [9] had interpreted the Carreau fluid transportation inside a rectangular domain operating as a peristaltic pump, asymmetric geometry, Nadeem and Akram [10] had disclosed the peristaltic pumping activity for asymmetric domain by using the regular perturbation expansion method, Akbar et al. [11] had highlighted the impact of nanofluids on peristaltic pumping through asymmetric domain, Akbar [12] had also considered the magnetohydrodynamic effects on the peristaltic pumping through asymmetric domain and the peristaltic pumping through an elliptic duct via exact solutions was interpreted by Saleem et al. [13].
Nanofluids consist of small-scale (1-100 nm in size) nano particles that are scattered in the base fluid. Scientifically, experts prefer multi-walled carbon nanotubes (MWCNTs) been a special form of carbon nanoparticles because they are nested inside one another and possess elongated hollow cylindrical tiny molecules of sp2 carbon [14,15]. These are useful in industrial problems where thermal conductivity is a prime factor [16]. Akbar and Nadeem [17] had reported a mathematical study that elucidates the combine analysis of peristaltic flow and nanofluids. Akbar and Butt [18] had mathematically elucidated the peristaltic flow of water as a base fluid with nano particles having distinct structures. Akbar [19] had numerically examined the nanofluid flow inside asymmetric domain via Runge-Kutta technique. Nadeem et al. [20,21] had evaluated the nanofluids impact on the peristaltic pumping through an annular in addition to rectangular domain. Akbar et al. [22] evaluated the magnetohydrodynamic impacts on the nanofluid pumping through a vertically held asymmetric domain.
The consequence of ciliated walls on the peristalsis driven flow is also a subject of keen interest. Cilia are tiny structures (like human hair) that play a prime role in the flow development. The peristaltic flow mechanism is also studied for distinct geometries having ciliated boundaries like ciliated cylindrical tube, Butt et al. [23] had provided exact solutions on the Phan Thien Tanner non-Newtonian flow inside a ciliated cylindrical domain, Saleem et al. [24] had interpreted the effects of hybrid nanofluid on the peristaltic pumping through a ciliated curved domain, Akbar and Khan [25] highlighted the impacts of MHD on the peristaltic pumping via a ciliated asymmetric domain. Sadaf and Nadeem [26] had presented a mathematical model that relates the combine effects of ciliated walls and peristalsis with nanofluid flow. Akbar [27] had reported the mathematical investigation of peristaltic flow of CNT in a symmetric geometry having ciliated boundaries. Akbar and Butt [28] elucidated the impacts of MHD and entropy interpretation on the peristaltic pumping of water based nanofluids through a cylindrical ciliated domain. Nadeem and Sadaf [29] had investigated the trapping phenomenon of nanofluids flow via a ciliated annular domain. Vaidya et al. [30,31] had highlighted the impacts of chemical reactions on the peristaltic pumping of Ree-Eyring non-Newtonian fluid passing through a permeable domain in addition to an asymmetric domain under MHD effects. Ashraf et al. [32] had examined the MHD impacts on nanofluid flow of Casson fluid through a cylindrical domain with sinusoidal wall fluctuations. Samuel et al. [33][34][35] had interpreted the effects of MHD, nanofluidic volume fraction with convection analysis through numerical techniques. Saleem et al. [36] highlighted the impacts of activation energy and magnetohydrodynamics on ciliary propelled flow of nanofluid. Ramesh et al. [37] modelled the cilia driven couple stress fluid flow under MHD effects. Farooq et al. [38] had analysed the micro polar flow developed due to cilia propelled rhythmic wave. Akbar et al. [39] modelled the pressure driven flow generated due to rhythmic activity of cilia.
A thorough and careful investigation of the available literature has disclosed that the peristaltic flow of carbon nanotubes in an elliptic duct with ciliated walls is not mathematically reported yet. Thus, the intent of the present mathematical computation is to elucidate this important topic for an elliptic domain for the first time. The base fluid considered in this mathematical study is water and multi-wall carbon nanotubes are taken into account. The viscous dissipation effect is also incorporated for a detailed study of convection heat transfer. Exact mathematical solutions are computed and the results are verified through graphical assessment. Streamlines depict the visualization of flow pattern.

MATHEMATICAL MODEL
The peristaltic flow of carbon nanotubes in an elliptic duct with ciliated walls is mathematically investigated. A Cartesian coordinate system approach is utilized to interpret this mathematical study. The geometry of the problem is given by Figure 1 providing the model of elliptic domain having ciliated walls. The sinusoidal wall fluctuations are evident in Figure 1 and the ciliated wall are also mentioned. Figure 1 contains the front together with a side view of this geometry, the front view highlights the elliptic domain of this duct.
The envelope model of cilia motion that considers the elliptical shape movements is taken into account. The equations considered for this model incorporates the sinusoidal metachronal wave movements generated due to rhythmic cilia motion. The computational equations that narrate this envelope model of cilia tips are referred as [40].
The following computational equations are used to express the velocities of cilia tips (i.e. axial and radial velocity respectively). (2) The simplification of Eqs 1, 2 provides the following equations The above mentioned velocities W and V mainly differentiate between the two strokes of cilia motion, (i.e., effective and recovery strokes respectively). The cilia tips are assumed to circulate in an elliptical course, as revealed in Figure 2. Cilium remains inflexible and rigid during its effective stroke. On contrary, it shows flexibility during its recovery stroke and retreats loosely.
The mathematical modelling for the incompressible flow with carbon nanotubes and heat convection is provided in dimensional form as follows Density Thermal expansion coefficient and density relation The unsteady and steady frames of reference are related by using the following equations The non-dimensional variables that are incorporated in present analysis are given as Whenever we deal with a non-circular domain then the concept of hydraulic diameter needs to be considered. It is considered as four times the cross-sectional area divided by the wetted perimeter of cross-section. The mathematical expression for the hydraulic diameter of elliptic duct is narrated as The computational equation that narrates the eccentricity of ellipse is e 1 − δ 2 , where δ is the aspect ratio and the second kind elliptical integral E(e) is referred as Yang et al. [41].
The useful transformations provided in Eq. 10, the dimensionless variables given in Eq. 11 are used in Eqs 5-9 and the approximation (λ → ∞ ) is used to avail the following set of dimensionless equations zp zy 0, Eq. 17 contains the dimensionless brinkman number B r that is the product of two dimensionless numbers Prandtl (P r ) and Eckert (E r ), i.e., B r P r E r . The dimensionless boundary conditions over the elliptic domain are , for θ 0, for Where a E(e) π [ 1 δ + ϕSin(2πz)], b E(e) π [1 + ϕSin(2πz)]. Where a and b provide the dimensionless form of wall functions that tackle the sinusoidal movement of travelling boundaries. Table 1 Akbar [42], provides the thermophysical properties of base fluid water and nanofluid. Table 2 Akbar [27], shows the mathematical properties of nanofluid.

EXACT SOLUTION
The dimensionless form of momentum equation is a nonhomogeneous partial differential Eq. 16 subject to nonhomogeneous boundary conditions (18). In order to solve this PDE, a polynomial solution method to solve PDE's over elliptic domain is taken into account. Consider the solution in a six constants polynomial form for velocity profile as narrated below Substituting Eq. 20 into Eq. 16 and equating the coefficients of x 2 , y 2 , x 0 , y 0 on both sides of equation, we get Now using Eq. 20 in boundary condition (18) and equating the coefficients of x 4 , x 2 , x 0 , we have following equations The simultaneous solution of above equations gives the following values of constants The values of these constants is inserted in Eq. 20 and exact mathematical solution is obtained for velocity as follows w x, y dp dz The non-dimensional mathematical expression for flow rate q(z) is availed by integrating Eq. 21 over the elliptic domain. The pressure gradient is calculated by using Eq. 22 and given as dp dz The pressure rise for a single wave-length is computed by using The technique that was utilized to get an exact velocity solution is employed again to get an exact temperature function given as

ANALYSIS OF RESULTS
The results obtained in the exact solution section, that is exact solution of momentum equation given by Eq. 21 exactly satisfies the corresponding momentum Eq. 16 and boundary condition (18). Further, the exact solution of temperature profile given by Eq. 25 also satisfies the corresponding energy Eq. 17 with boundary condition (19). The validation of these solutions is also confirmed by using Mathematica software. According to the modelled problem and relevant boundary conditions, the flow profile in addition to the temperature profile should have a maximum value in the centre of elliptic domain and both the profiles should decline towards the walls of channel. Thus, it is evident from the graphical results that flow and temperature are maximum in the centre and reduces towards the boundaries. The temperature graphs show that temperature profile declines toward the boundaries and becomes zero as it was considered in the boundary condition (19). If x 2 a 2 + y 2 b 2 1 is inserted in Eq. 21 then the reduced form gives the boundary condition (18). Further, if x 2 a 2 + y 2 b 2 1 is inserted in Eq. 25 then the reduced form provides the boundary condition (19). This also validates the considered boundary conditions. FIGURE 11 | (A) dp dz plot against axial coordinate for aspect ratio. (B) dp dz plot against axial coordinate for occulsion. (C) dp dz plot against axial coordinate for volume fraction of nanofluid. (D) dp dz plot against axial coordinate for flow rate.

DISCUSSION OF RESULTS
The graphical results that are incorporated in this section provide a comprehensive evaluation of the mathematical computations given in the exact solution segment. A thorough graphical evaluation is considered for velocity, temperature and various physical characteristics of peristalsis. Figures 3-5 unfold the effect of some important physical parameters on the velocity profile. Figure 3 (i.e., Figures 3A,B) provides the graphical outcome of velocity profile for increasing flow rate Q. Figure 3A reveals that velocity depicts an increasing behaviour for the escalating numerical values of parameter Q. Figure 3B shows the dependence of velocity on both the independent variables x and y and it is 3D-graphical plot of velocity for the escalating flow rate. An axial symmetry flow profile is noted here in this graph and maximum flow velocity is observed in the centre while it is diminishing towards duct walls. Figure 4 (i.e., Figures 4A,B) depicts the graphical result of velocity for the distinct escalating numerical values of wave number. Figure 4A unfolds that velocity shows an escalating behaviour for the distinct increasing numerical values of β in the centre of elliptic duct while a converse behaviour is observed near walls. Figure 4B gives a 3D-graphical plot of velocity for distinct escalating values of β. Figure 5 (i.e., Figures 5A,B) gives the graphical assessment of velocity for the escalating numerical values of α. Figure 5A discloses that velocity increases with the escalating numerical values of α in the centre of duct while a converse behaviour near the walls is noted also for escalating numerical values of α. This converse behaviour of velocity at the ciliated walls is due to increasing eccentricity of elliptical motion of cilia movement. Figure 5B discloses the 3Dgraphical result of velocity for distinct escalating values of α. All these graphical outcomes of velocity reveal an axially symmetry flow that is a fully developed flow and velocity is diminishing towards walls of duct. Figures 6-10 unfold the graphical outcomes of temperature profile for distinct parameters of physical importance. Figure 6 (i.e., Figures 6A,B) reveals the graphical assessment of temperature for distinct escalating numerical values of Q. Figure 6A unfolds that the temperature depicts an escalating behaviour for distinct incrementing numerical values of Q. Figure 6B depicts the 3D-graphical plot of temperature for different escalating numerical values of Q. An axially symmetric behaviour is observed for temperature profile. Figure 7 (i.e., Figures 7A,B) reveals the graphical result of temperature profile for various escalating numerical values of ψ. The parameter ψ is representing the concentration of carbon nanotubes in the base fluid. Figure 7A depicts a diminishing temperature profile with escalating concentration of carbon nanotubes in the base fluid (i.e., water). Since the increasing concentration of carbon nanotubes in the base fluid escalates the thermal conductivity of fluid and this eventually results in decline of temperature of base fluid. Further, it also verifies the inverse relation between the thermal conductivity of fluid and the temperature of fluid. Figure 7B is the 3D-graphical plot of temperature profile for escalating numerical values of ψ. Figure 8 (i.e., Figures 8A,B) is plotted to represent the effect of β on temperature. Figure 8A discloses that the increment in the numerical value of wave number results in the escalation of temperature. Figure 8B is the 3D-graphical plot of temperature for distinct escalating numerical values of β. Figure 9 (i.e., Figures 9A,9B) is the graphical assessment of temperature profile for distinct escalating numerical values of B r Figure 9A unfolds that the temperature escalates rapidly with even a small increase in the numerical value of B r It reveals that the viscous effects are playing a vital role in the heat enhancement as compared to the molecular conduction. Figure 9B unfolds the 3D-graphical plot of temperature profile for escalating numerical values of B r This graph beautifully shows the dependence of temperature profile on the independent variables x and y. Figure 10 (i.e., Figures  10A,B) is the graphical representation of temperature for escalating numerical values of α. Figure 10A discloses an escalation in the temperature profile for distinct incrementing numerical values of α. Figure 10B provides a 3D-graphical plot of temperature for enhancing numerical values of α. All these temperature graphs reveal an axial symmetry profile of temperature. Further, temperature gains maximum escalation in the centre of duct and drops to zero at the duct walls. Figures 11A-D provide a graphical assessment of dp dz for different escalating numerical values of distinct parameters. Figure 11A reveals an escalation in the value of dp dz for different enhancing numerical values of δ. Figure 11B depicts an escalation in the value of dp dz in the first half length of axial coordinate while a drop in the value of dp dz is noted in the next half length of axial coordinate for escalating numerical values of ϕ. Basically, the escalation is noted here for crest of peristaltic wave while the drop is noted for the trough. This drop in the pressure gradient for the trough of peristaltic wave eventually accelerates the flow in the axial direction. Figure 11C depicts a drop in the value of dp dz for distinct accelerating numerical values of ψ. Figure 11D reveals that the value of dp dz diminishes for the escalating numerical values of Q. Figures 12A-C provide the graphical results of ΔP plotted against the dimensionless flow rate Q. Figure 12A unfolds an escalation in the value of ΔP in the region ΔP > 0 while a converse behaviour is observed for the region ΔP < 0 with enhancing numerical values of δ. Figure 12B reveals that ΔP escalates in the region ΔP > 0 and drops in the region ΔP < 0 with enhancing numerical values of ϕ. Figure 12C unfolds a decline in the value of ΔP with escalating concentration of carbon nanotubes ψ. Figures 13A-D provide the graphical visualization of flow in the elliptic duct. These graphical results are obtained for distinct escalating numerical values of flow rate. The escalation in the trapping size can be clearly seen in these graphs for incrementing Q. Streamlines disclose the behaviour of flow regime inside this ciliated elliptic duct.

CONCLUSIONS
The peristaltic flow of carbon nanotubes in an elliptic duct with ciliated walls is mathematically investigated. A Cartesian coordinate system approach is utilized to interpret this mathematical study. The prime outcomes are given below.
• All the graphical outcomes of velocity reveal an axially symmetry flow that is a fully developed flow and velocity is diminishing towards walls of duct. • The increasing concentration of carbon nanotubes in the base fluid escalates the thermal conductivity of fluid and this eventually results in decline of temperature of base fluid. • The inverse relation between the thermal conductivity of fluid and the temperature of fluid is noted. • All the temperature graphs reveal an axial symmetry profile of temperature. Further, temperature gains maximum escalation in the centre of duct and drops to zero at the duct walls. • The escalation in the trapping size can be clearly seen in these graphs for incrementing Q.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.