Abstract
The main focus of this article is to mathematically formulate the microfluidics-based mechanical system for nanofluids. A 50:50 mixture of propylene glycol (PG) and water is used as a heat transfer fluid because of its tremendous anti-freezing properties, and nontoxicity and it is safe to be utilized at the domestic level. Titanium dioxide (titania) nanoparticles are suspended in the working fluid to enhance its heat transfer ability. The fluid flow is induced by electroosmosis in a microtube, which is further assisted by cilia beating. The impacts of Joule heating and non-linear thermal radiation are also considered. The simplification of the dimensionless system is done under lubrication theory and the Debye-Hückel linearization principle. The nonlinear system of equations is executed for a numerical solution by adopting the symbolic mathematical software Maple 17 using the command “dsolve” along with the additional command “numeric” to get the numerical solution. This command utilizes a low-ordered method along with accuracy-enhancing schemes such as the deferred correction technique and Richardson extrapolation to get a numerical answer of desired accuracy, where we can choose the accuracy level and mesh points according to our requirements. The detailed analysis of results obtained from the numerical treatment of the considered problem indicates that the efficiency of the PG + water enhances due to the suspension of the nanoparticles and heat is rapidly removed from the system. Further, the velocity of the fluid is augmented by decreasing the thickness of the electric double layer and raising the strength of the electric field in the forwarding direction.
Similar content being viewed by others
Introduction
Most substances acquire a net surface charge when they come across an aqueous polar medium. Due to the interaction between a charged surface and ionic species, the diffuse layer is generated. The ions in the diffuse layer experience an axial body force in the presence of an external force. Such a kind of motion is termed as electroosmosis which is considered as a very suitable mechanism for transportation purposes in microfluidic devices. In recent years, microfluidic devices have been triggered by the growing needs of biomedical engineering as well as energy systems. For the fabrication of microfluidics devices, there is a need to develop theoretical as well as experimental models of electroosmotic flow with heat transfer assisted by the various pumping mechanisms. Xuan et al.1 have presented a basic model of electroosmotic flow where they have discussed the joule heating effects. Recently, Khan and Sasmal2 investigated the electroosmotic flow phenomenon for both Newtonian and viscoelastic fluids through a porous environment containing long micropores. An investigation on the coupled impact of pressure-driven flow and electrokinetics is carried out by Yuan et al.3 where they did not ignore the adverse pressure impact on electroosmotic flow. A mathematical model for electroosmotic transport of two immiscible fluids is developed and analyzed by Alyousef et al.4 in which they utilized the Ellis fluid model for the viscoelastic behavior of the considered fluid. Cilia beating is a natural mechanism of physiology that generates effective stroke and recovery stroke. Due to an effective stroke, water/fluid moves backward and cilia move forward. Due to recovery strokes, cilia further come back to their original position. These effective recovery strokes are rhythmic processes and generate metachronal waves. Cilia beating is a key mechanism of the respiratory system that cleans dust particles. In some physiological systems, cilia work as sensors and actuators. Some of the basic studies5,6,7,8,9 on cilia transport and its roles have been reported. From a natural perspective, it is an essential mechanism of the physiological system. Based on the extensive role of this mechanism, nowadays, there are challenging demands for artificial cilia and artificial magnetic cilia for the sensors and actuators10, for rotating machines to compute the frequency11, for fluid pumping12, and many more. Inspired by the natural and artificial roles of the cilia beating in fluid dynamics and biomedical sciences, most recently researchers have developed some mathematical models on the cilia-assisted rheological fluid flow13,14,15,16,17,18,19 and MHD fluid flow20,21,22,23 where they have discussed the effects of key parameters such as cilia length, the eccentricity of the metachronal waves, curvature of the flow geometry, magnetic field and rheological parameters on fluid pumping. In order to enhance the efficiency of microfluidic devices, different flow-causing mechanisms can be combined to achieve improved flow and heat transfer rates. The electroosmotic flows, triggered by ciliary motion can be helpful in the development of electromechanical devices at the microscale such as in micro and nanofluidic pumps, hemodialysis, cooling chips, microfabricated liquid devices, etc. Owing to the stunning efficiency of combined ciliary-actuated electroosmotic flows, recently few researchers have put forward their efforts in investigating the properties of flow driven by the combined phenomenon. Abdelsalam and Zaher24 studied the impact of electroosmotic forces on the movement of sperm in a ciliated cervical canal in which they employed the hyperbolic-tangent fluid to model the cervical fluid. Imran et al.25 investigated the heat transfer properties of Williamson fluid flow where the fluid flow is generated by electroosmosis and cilia beating in the micro-ciliated channel. Some more significant studies are reported by26, Ijaz et al.27, and Javed et al.28.
Solar thermal technology consists of setups that collect the energy from the sun and convert it into heat energy. This heat energy is then stored by using water as a working fluid and utilized in various applications. The components of a solar heating system contain a source for collecting the solar energy, i.e. collector, a means of exchanging the heat using the flow of the fluid, namely the heat transfer circuit, the heat transferring fluid, and an energy storage system. In a climate, where there occur extreme temperature conditions in different parts of the year, the usage of anti-boiling and anti-freezing fluids is recommended in order to avoid damage to the solar system due to corrosion, overheating, and freezing. For this purpose, multiple types of anti-freezing agents are present, such as ethylene glycol (EG), tri-ethylene glycol, propylene glycol (PG), etc., which have the capability of protecting the solar system from freezing and rust. But these antifreeze agents in their pure form do not possess good thermal properties, which is the main requirement of a solar system. Therefore, a mixture of these anti-freezing liquids with water is utilized for this purpose. Out of different antifreeze/water solutions, the most commonly used solution is propylene glycol/water solution. Due to its low cost, high performance, low corrosivity, least toxicity, and environmentally friendly nature, as found by Heinonen et al.29 in a comparative study, it is the most suitable anti-freezing fluid, especially for use in domestic-level devices. Shojaeizadeh et al.30 performed an experiment to investigate the performance of a flat plate solar cell using different concentrations of PG in water. They observed that upon increasing the concentration of PG in water from 25 to 75%, an improvement in the efficiency of the solar cell occurs. Jugar and Crook31 compared the performance of PG and EG with their different concentrations in water and concluded that it is necessary to maintain a 50/50 concentration of water and antifreeze in order to keep a balance between heat transfer properties and antifreeze protection.
Titanium dioxide (TiO2), which is one of the oxides of titanium, was first discovered from ilmenite in 1791. Due to its high refractive index, and low chances of decolorization, it is mainly used in pigments32. However, other than its large-scale usage in pigments, TiO2 is also used in enamels, water purification, cosmetics, food items33, and in energy storage devices. Titanium dioxide is most beneficial when utilized in the form of small particles (usually nanosized) due to its enhanced light absorption tendency, photocatalytic reduction, and high surface photoactivity. Titanium dioxide nanomaterials are employed as photocatalysts in water purification as they can effectively degrade the organic and inorganic pollutants present in wastewater34,35. Due to the semiconducting nature of titanium dioxide nanoparticles, they are widely utilized in electronic components, catalytical electrodes, and solar cells. Titanium dioxide nanoparticles are very preferable to be used in domestic-level devices due to their non-toxic nature36 and easy and cost-effective availability.
The heat transfer mechanism offered by ordinary fluids such as ethylene glycol, water, engine oils, etc. is highly inefficient due to their innate poor conductivity. As already mentioned the heat transfer process is one of the major factors affecting the efficiency of solar cells but the addition of antifreeze in water lowers its heat transfer ability. In order to overcome the demands of energy-efficient devices in this modern era, a new class of fluid offering potentially improved heat transfer rates is required. It was Choi37, who contributed his efforts to the generation of a new class of fluids named as “nanofluids” offering exciting applications in improving the heat transfer ability of the fluids by the insertion of nanosized particles (1 nm–100 nm) in the conventional heat transfer fluids. The extraordinarily high heat transfer ability of these new-generation fluids triggers their usage in pharmaceuticals, nano-scaled lubrication in rocket engines, cooling/heating of electronic devices, vehicle thermal management, solar cells, and many more. Such promising applications of nanofluids attracted the fluid dynamists to investigate such kinds of fluid and even more innovative fluids having more than one type of solid nanoparticles which are usually named as “hybrid nanofluids”. In an experimental study performed by Javidan and Moghadam38 an effective cooling for photovoltaic cells is tested by using SiC/water nanofluid as a cooling agent. Saleem et al.39 investigated the importance of thermophoresis and Brownian motion on the motion of water consisting of three different types of nanoparticles. Elnaqeeb et al.40 presented the importance of dual stretching and suction in the flow of water-based ternary-hybrid nanofluids by assuming multiple shapes and densities of nanoparticles. Animasaun et al.41 studied the motion of water, containing carbon nanotubes, copper, and graphene nanoparticles of different shapes, over stagnant moveable walls. They concluded that a minimum value of local skin friction coefficient can be achieved when a larger fraction of nanoparticles are used. Many other investigations focusing on the applications of multiple types of nanofluids are contributed by Rasool et al.42, Adnan and Ashraf43, Xiu et al.44, and Shah et al.45.
Due to the above-mentioned advantages of the propylene-glycol/water mixture and titanium dioxide nanoparticles, a combination of these two as a nanofluid can be adapted to work as heat transfer fluid in domestic-level devices. Propylene-glycol/water-based TiO2 nanofluid can provide a sufficient rate of heat transfer along with the most important constraint of non-toxicity. So the aim here is to investigate the heat and mass transfer capability of propylene-glycol/water-based TiO2 nanofluid which is driven by the combined influence of cilia beating and the electroosmotic forces. Nevertheless, none of the studies have focused on the mathematical model of the cilia-assisted electroosmotic flow of nanofluids which could be applicable in various biomedical and energy fields. Based on research gaps, herein, a mathematical model is presented to study the propylene glycol + water-based TiO2 nanofluids flow driven by the electroosmosis and cilia beating in a microtube. The present results suggest that this model can be applied to the household solar system. This model is structured in the form of an introduction followed by mathematical formulation and numerical solutions. Thereafter, a thorough discussion of simulated results has been presented. In last, the concluding remarks of the present analysis have been listed.
Mathematical formulation
Ciliary movement
Here we have considered the characteristics of fluid motion of 50:50 ionic water-propylene glycol (PG) solution-based titanium dioxide nanoparticles through a ciliated tube. A metachronal wave produced by the cumulative beating of cilia travels with the speed c towards the right (see Fig. 1). Further, the electroosmotic motion in the fluid is generated by the implementation of an external electric field across the tube in the z-direction. The envelope for the elliptical path followed by cilia tips is represented mathematically by choosing the polar coordinate system (\(\tilde{r},\tilde{z},\tilde{t}\)) and given as46:
in which d designates the tube radius, \(\epsilon_{1}\) the dimensionless length of the cilia, \(\tilde{z}_{0}\) the cilia reference point and \(\tilde{\alpha }\) the eccentricity for the elliptic path followed by cilia.
Imposing no-slip conditions across the tube walls, the velocity of the fluid adjacent to the cilia is the same as the velocity of the cilia, therefore, the axial \(\tilde{w}\) and radial velocity \(\tilde{u}\) of the cilia at the reference point can be expressed as:
Using the relations given in Eqs. (1)–(2) in Eqs. (3)–(4) result in the following expressions:
Governing equations
This analysis focuses on the flow phenomenon of PG/water solution-based titanium dioxide (TiO2) nanoparticles driven by the combined metachronal wave and the electroosmotic body forces. The considered nanofluid exhibits the Newtonian characteristics which are governed through the Navier–Stokes equation. The conservation equations for the nanofluid are formulated by employing the Buongiorno flow model along with the relation of thermal conductivity and viscosity proposed by Corcione’s model. The electric potential distribution within the fluid medium is governed by the linearized Poisson-Boltzmann equation and simplified under the Debye-Hückel approximation of low zeta potential. The heat transfer analysis is performed in the presence of nonlinear thermal radiation and Joule heating. The slip conditions for temperature and zero mass flux for the nanoparticle volume fraction are applied across tube walls.
Subject to above defined physical conditions, the constitutive equations for momentum, energy, and concentration equations are formulated as:
The nonlinear radiative heat flux computed by Roseland’s approximation is given by:
where \(\sigma^{*}\) and \(\kappa^{*}\) symbolize the Stefan-Boltzmann constant and mean absorption coefficient respectively.
In the above equations, \(\tilde{p}\) stands for pressure, \(\rho_{nf}\) for the nanofluid density, \(\rho_{e}\) for the electric current density, \(U_{{E_{{\tilde{r}}} }}\) and \(U_{{E_{{\tilde{z}}} }}\) for the radial and the axial electric forces, \(\left( {\rho \gamma } \right)_{nf}\) for the thermal expansion coefficient, \(\beta_{{\tilde{\Phi }}}\) for the coefficient of mass expansion, \(\tilde{T}\) and \(\tilde{\Phi }\) for temperature and the nanoparticle volume fraction respectively, \(\sigma_{nf}\) for the nanofluid electric conductivity, \(D_{T}\) and \(D_{B}\) for the coefficient of thermophoretic and diffusion motion respectively and \(\left( {\rho C} \right)_{nf}\) for the specific heat capacity of nanofluid.
The relations for nanofluid properties are expressed as40:
with \(\rho_{bf}\) and \(\rho_{p}\) representing the densities of the base fluid and titanium oxide nanoparticles and \(\Phi_{0}\) the fraction of nanoparticles being scattered in the base fluid.
Corcione’s model for thermal conductivity is given by47:
In the above relations, \(T_{fr}\) specifies the freezing temperature of the base fluid, \(R\) the nanoparticles’ Reynolds number, Pr the Prandtl number and \(d_{bf}\) and \(d_{p}\) the diameters of base fluid molecules and nanoparticles respectively.
The Maxwell–Garnett model for electrical conductivity is expressed as48:
The Poisson equation for electric potential distribution in the fluid medium, due to by motion of ions generated by the external electric field, derived in polar coordinates is given by:
where εr and \(\varepsilon_{0}\) designate the relative permittivity for the base fluid and the permittivity of vacuum respectively.
The electric charge density in terms of the number densities of anions \(n^{ - }\) and cations \(n^{ + }\) having valency \(z\) is defined as:
Defining the transformation relations between a fixed frame of reference and wave frame
Scaling analysis
The following dimensionless parameters are defined as:
In which, \({\text{Re}}\) designates the Reynolds number, \(U_{0}\) the electroosmotic velocity parameter, \(k\) the Debye length parameter, \(N_{t}\) the dimensionless thermophoretic parameter, \({\text{Gr}}_{t}\) the thermal Grashof number, \(N_{b}\) the Brownian motion parameter, \(R_{d}\) the radiation parameter, \(S\) the Joule heating parameter, \({\text{Pr}}\) the Prandtl number, \(\theta\) the dimensionless parameter, \({\text{Gr}}_{\Phi }\) the mass transfer Grashof number, \(\delta\) the wavenumber, \(\eta_{w}\) the temperature ratio parameter, \(\Psi\) the stream function, and \(\Phi\) the dimensionless nano fraction parameter.
Using Eqs. (20) and (21) in the dimensionless analysis of Eqs. (7)–(11) and (18)–(19) and adopting the assumption of small wavenumber and the domination of viscous forces over inertial forces, we get the following simplified set of equations,
The Boltzmann distribution for the local ionic density of each ionic species is defined as:
Substituting the Boltzmann relations in Eq. (27), we get,
In order to linearize Eq. (29), we can make use of Debye–Hückel linearization principle which is basically assumption of smaller zeta potential in the diffuse layer. Use of this linearization principle is quite justified as for a wide range of PH of electrolyte solution, zeta potenetial remains lower than 25 mV so we can make use of assumption \(\sinh \left( {U_{E} } \right) \approx U_{E}\) in Eq. (29) which reduces to
Subject to the appropriate boundary conditions defined below
Equation (30) can be solved directly to get an analytical expression for the potential field as:
Here \(I_{0}\) represents the modified Bessel functions of 1st kind.
The suitable boundary conditions in the dimensionless form are prescribed as:
where \(\beta\) represents the temperature slip parameter.
The instantaneous volume flow rate in the laboratory frame is given as
The dimensionless time-averaged mean flow rate in the laboratory and moving frame is defined as
So we get
in which \(Q\) is a time-mean flow rate calculated for a single period of the wave.
Solution methodology
The linearized Poisson-Boltzmann solution is solved analytically and the resulting expression for electric potential is utilized in the momentum equation. The simplified set of governing equations is defined in Eqs. (24)–(26) is highly nonlinear and cannot be executed for the analytical solutions. Therefore, the system of governing equations along with boundary conditions given in Eq. (33) is numerically simulated using suitable symbolic computer software i.e. Maple 17 and Mathematica. Furthermore, the graphical results for the flow properties such as temperature, velocity, Nusselt number, and streamlines are computed through Maple 17 to highlight the impact of various physical parameters of the interest.
Analysis and discussion of results
The main purpose of this section is to outline the impact of the various physical factors on the flow phenomenon of propylene-based titanium dioxide nanofluid through graphical results. The nanoparticles of an average diameter of about 10 nm are dispersed at an initial temperature of 303 K. The fluid flow pattern is visualized by plotting the contour graphs for stream function. Graphical results are prepared by choosing the direction of the applied electric field in such a way that assists the flow of the fluid generated by the beating of the cilia. However, the effect of removing the electric field and reversing the direction of the electric field is also observed. From the mathematical point of view, negative values of the electroosmotic velocity parameter represent the forwarding direction electric field and positive values characterize the opposing electric field. A volume fraction of up to 0.06 vol% of titanium dioxide nanoparticles is scattered in the PG + water solution. The thermophysical attributes of the base solution and nanoparticles are calculated at 303 K and are listed in Table 1. Based on the numerical values given in Table 1, the value of the Prandtl number is found to be 45.078.
Analysis of results
This subsection is devoted to facilitating an insight into the influence of various involved on the velocity distribution of PG + water solution base titanium dioxide nanofluid. The effect of the cilia length parameter on the fluid velocity is characterized in Fig. 2. It can be found that for increasing the length of the cilia in the range of \(0.5 \le \epsilon_{1} \le 0.65\), the velocity of the fluid decreases for \(0 \le r \le 0.8\) and it increases in the region defined by \(0.8 \le r \le 1.5\). The effect of an increase in the eccentricity of the path followed by moving cilia the on velocity profile is analyzed in Fig. 3. Here a rise in the velocity is observed in the central region of the pipe. A noticeable uplift in the fluid motion is produced with a rise in the wave number in the range of \(0.05 \le \delta \le 0.20\) as seen in Fig. 4. Figure 5 delineates the impression of the Debye length parameter on the flow phenomenon. A significant rise in the velocity of the fluid in the central region is observed. Another important physical parameter involved in the investigation of the electroosmotic flow phenomenon is the electroosmotic velocity parameter whose effect on fluid velocity is analyzed in Fig. 6. Here the graphical results are being drawn for \(U_{0} > 0\), \(U_{0,} < 0\) and \(U_{0} = 0\) and it is found that velocity is maximum for \(U_{0} < 0\) and it is minimal for \(U_{0} > 0\). The impact of variation in the fraction of dispersed titanium dioxide nanoparticles on velocity under the same physical conditions is depicted in Fig. 7 and a decline in velocity profile is observed in the resulting panel.
The trapping phenomenon is one of the key features associated with the flows following the pattern of peristaltic pumping. In this process, few streamlines under special circumstances enclose a quantity of the fluid known as trapped bolus which is carried along with the metachronal wave. The area enclosed by the trapping bolus and its shape is mainly affected by involved physical parameters. Therefore, in Figs. 8, 9, 10, 11, contour graphs are plotted for multiple values of the electroosmotic velocity parameter, Debye length parameter, eccentricity, and cilia length parameter to visualize their impact on the circulatory flow pattern. Figure 8a–c illustrate the effect of the electroosmotic velocity parameter on the trapping process. Evidently, with the inclusion of the external electric field in the direction of the main flow direction, less volume of the fluid is trapped when compared with the case of the absence of electric field and the reversal of the electric field direction. Furthermore, the number of closed streamlines is maximum for opposing electric fields and minimum for assisting the electric field. The evolution in the circulatory flow pattern for the rise in the Debye length parameter is depicted in Fig. 9a–c. It can be clearly seen from the resulting panels that there is an increase in the volume occupied by the closed streamlines retaining the same number of closed streamlines in each case. The variation in the size of the trapping bolus via a larger eccentricity parameter is revealed through Fig. 10a–c. Minor growth in the size of the trapped bolus is observed for a rise in the eccentricity of the path followed by the beating of the cilia. Clearly, longer cilia trace a bigger ellipse and the fluid particles attached to these cilia enclose a comparatively larger fluid bolus. A prominent growth in the circulatory flow pattern is observed for enhancement in the characteristic length of the cilia as manifested in Fig. 11a–c. Thus it can be concluded that the length of the cilia generating the metachronal wave is the most significant parameter affecting the fluid flow dynamics.
The alteration in the temperature subject to variation in different embedded parameters is displayed in Figs. 12, 13, 14, 15, 16, 17. The modification in the temperature profile for rising values of the radiation parameter from 2.0 to 3.5 is indicated in Fig. 12. A remarkable reduction in temperature profiles is observed for increment in the thermal radiation parameter. The development in the temperature distribution for multiple values Joule heating parameter in the range of \(4.5 \le S \le 6.0\) is analyzed in Fig. 13. It is clear from the resulting sketch that there is a rise in temperature of the nanofluid subject to a rise in S. The impact of the thermal slip parameter on temperature distribution is elaborated in Fig. 14. Here a consistent augmentation in the temperature of the fluid is observed. Figure 15 is plotted to analyze the variation in the thermal distribution for higher values of the temperature ratio parameter from \(1.1 \le \eta_{w} \le 1.4\). A substantial suppression in the temperature profile is depicted in the resulting graph. Figure 16 portends a significant drop in the temperature of the fluid for increment in the nanoparticle volume fraction from 0.03 to 0.06 in the base fluid. Figure 17 is prepared to examine the influence of the cilia length parameter on temperature. The temperature of the nanofluid is elevated when longer cilia are used for generating the metachronal wave.
The non-dimensional Nusselt number is expressed as:
Figures 18, 19, 20, 21, 22 are sketched to observe the characteristics of the Nusselt number for variation in different involved parameters. Nusselt number is a parameter that facilitates measuring the heat transfer ability of the working fluid. The variation in the distribution of the Nusselt number along the micro tube walls for rising values of the temperature ratio parameter is assessed in Fig. 18. One may note that there is a depression in the magnitude of the Nusselt number when the temperature ratio is increased. Nusselt number distribution for multiple values of the Joule hearing parameter can be visualized in Fig. 19. The change in Nusselt number for rising values of the thermal radiation parameter is depicted in Fig. 20. The resulting graph reveals a decline in Nusselt number. It is noticed from Fig. 21 that for increasing the length of the cilia involved in generating the fluid flow, Nusselt number declines. The impact of the temperature slip parameter on the Nusselt number along the tube wall is elucidated in Fig. 22. An intensification is observed in the Nusselt number when more thermal slip is experienced by the fluid at microtube walls.
Discussion of the results
Physically a rise in the length of the cilia corresponds to an increase in the amplitude of the metachronal wave as well as in its wavelength, therefore, the fluid velocity is boosted for larger \(\epsilon_{1}\) as observed through Fig. 2. The reason behind the increment in velocity profile for increase in eccentricity parameter can be related to the fact that eccentricity parameter is directly proportional to the cilia length which causes a rise in the driving force generated by the beating of the cilia. The result obtained in Fig. 4 for larger wave number can be validated by the fact that when the wavenumber is increased, a dominance of inertial forces over viscous forces is produced which facilitates the fluid motion. Basically, the Debye length parameter is inversely related to the thickness of the electric double layer and a thin EDL corresponds to the uneven distribution of electric potential within the fluid medium. Consequently, the electroosmotic fluid flow is augmented and a rise in axial velocity is observed as depicted in Fig. 5. The Debye length parameter is one of the most important parameters in the electro-osmotically modulated flows and it can be visualized that it is a very influential parameter that is very helpful in controlling the speed of the fluid. As already mentioned that negative values of electroosmotic velocity represent that the axial electric field is oriented in the positive z-direction which physically means that in this case, the electroosmotic velocity is in the direction of the fluid flow generated by the motion of cilia and a positive value means that electroosmotic process is occurring is the backward direction opposite to main flow direction, therefore, it can be observed from Fig. 6 that velocity of the nanofluid is maximum for a negative value of \(U_{0}\) and it is minimum for positive values. However, \(U_{0} = 0\) physically corresponds to the absence of the external electric field and in this case, fluid motion is occurring only due to ciliary motion, therefore, the velocity profile lies in between the above-mentioned cases of assisting and resisting electric field. As the viscosity of the base fluid rises when the concentration of nanoparticles is enhanced which in turn resists the fluid flow, and a decline in velocity profile is observed in Fig. 7. In all of the Figs. 2, 3, 4, 5, 6, 7, an opposite response of the velocity profile is observed near the wall of the channel, which can be justified by the fact that in order to maintain the fixed flow rate, velocity of the fluid depicts the opposite trends.
For growth in the radiation parameter, the process of conductive heat transfer from the fluid medium is boosted which in turn reduces the temperature profile (See Fig. 12). The Joule heating phenomenon measure the resistance offered by electrolyte solution to the passage of current. A larger Joule heating parameter quantifies the larger amount of resistance experienced by electric current. As a result, more electric energy is converted into the heat energy and temperature of the fluid increases as observed in Fig. 13. As \(\eta_{w}\) is the ratio of temperature at the upper wall and at the lower wall, so larger \(\eta_{w}\) corresponds to an elevation in temperature difference which boosts the heat transfer rate within the fluid medium and the overall temperature of the fluid drops as noticed through Fig. 15. A keen examination of thermal conductivity relation for nanofluid given in Eq. (15) clarifies that the thermal conductivity of the nanofluid is directly influenced by the fraction of nanoparticles being suspended in the base fluid which justifies its effect of boosting the cooling phenomenon and decreasing the fluid temperature in Fig. 16. These enhanced thermal properties of the fluid boost the efficiency of the working fluid by helping in controlling the fluid temperature. When the cilia inside the tube have a relatively larger length, their movement causes a stronger pushing force on the fluid which causes the temperature of the fluid to increase as shown in Fig. 17.
Nusselt number quantifies the relative measure of heat transfer by convection as compared to conductive heat transfer. It can be concluded from Figs. 18 and 20 that for a larger temperature difference which is related with temperature ratio can boost the other modes of heat transfer such as heat transfer by radiation, heat transfer associated with Brownian and thermophoretic diffusion of nanoparticles, therefore a reduction is convective rate of heat transfer occurs and Nusselt number declines. Moreover, as Joule heating parameter \({\text{S}} = \frac{{\sigma_{bf} U_{{E_{z} }}^{2} d^{2} }}{{k_{bf} (T_{1} - T_{0} )}}\) is strongly dependent on the strength of applied electric field, which is the source of fluid motion in the forwarding direction in case of negative values of \(U_{0}\), therefore larger values of S raises the convective heat transfer and the Nusselt number increases significantly (See Fig. 19). The result in Fig. 21 can be justified by the fact that using relatively larger cilia facilitates the movement of fluid which raises the heat transfer due to convection. A larger temperature jump parameter quantifies larger temperature difference between the fluid and the solid surface, which corresponds to larger kinetic energy of the fluid particles. As a result, momentum diffusivity increases and Nusselt number also raises due to rapid movement of heated fluid as observed in Fig. 22.
Concluding remarks
Herein, a theoretical investigation has been conducted to study the heat transfer characteristics of 50:50 PG + water -based ionic solution whose motion is driven by the combined effects of electroosmotic body forces and the beating of the cilia. In presence of nonlinear thermal radiation and Joule heating, a nonlinear coupled system of governing equations is simulated by mathematical software Maple 17. The significant findings of the current analysis are listed as:
-
Propylene glycol and water-based Titania nanofluid is the most suitable heat transfer fluid in the household solar device which enhances heat transfer capability with the reduction of fluid temperature.
-
The application of the electric field along the positive axial direction boosts the fluid velocity which results in the rapid removal of the heat from the system.
-
The length of cilia involved in generating the fluid motion strongly influence the velocity, temperature, and trapping phenomenon and it also boosts the heat transfer tendency of the fluid
-
Increasing the thermal radiation parameter results in the reduction of fluid temperature.
Data availability
The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.
Abbreviations
- \(d\) :
-
The half-width of the channel \({\text{m}}\)
- \({\upeta }\) :
-
Wavelength \({\text{m}}\)
- \(\tilde{r}, \tilde{z}\) :
-
Radial and axial coodinates \({\text{m}}\)
- \(\tilde{u} ,\tilde{w}\) :
-
Radial and axial velocity component \({\text{m}}\,{\text{s}}^{ - 1}\)
- \(\rho_{nf} ,\rho_{fb} ,\rho_{p}\) :
-
The nanofluid, base fluid, and nanoparticles density respectively \({\text{Kg}}\,{\text{m}}^{ - 3}\)
- \(\tilde{t}\) :
-
Time scale \({\text{s}}\)
- \(\tilde{p}\) :
-
Pressure force \({\text{N}}\,{\text{m}}^{ - 2}\)
- \(\rho_{e}\) :
-
Local electric charge density \({\text{C}}\,{\text{m}}^{ - 3}\)
- \(\gamma_{nf} ,\gamma_{s} ,\gamma_{bf}\) :
-
Thermal expansion coefficient for nanofluid, solid particles and basefluid \({\text{K}}^{ - 1}\)
- \(U_{{E_{{\tilde{r}}} }} ,U_{{E_{{\tilde{z}}} }}\) :
-
The radial and axial electric body field \({\text{Kg}}\,{\text{ms}}^{ - 3} \,{\text{A}}^{ - 1}\)
- \(C_{nf} , C_{p} ,C_{bf}\) :
-
Specific heat of nanofluid and nanoparticles \({\text{J}}\,{\text{Kg}}^{ - 1} \,{\text{K}}^{ - 1}\)
- \(\tilde{T}\) :
-
Temperature field \({\text{K}}\)
- \(D_{B}\) :
-
Brownian diffusion parameter \({\text{m}}^{2} \,{\text{s}}^{ - 1}\)
- \(\tilde{\Phi }\) :
-
Concentration field \({\text{Kg}}\,{\text{m}}^{ - 3}\)
- \(\sigma_{nf} ,\sigma_{bf} ,\sigma_{p}\) :
-
Electrical conductivity of the nanofluid, base fluid, and nanoparticles. \({\text{S}}\,{\text{m}}^{ - 1}\)
- \(D_{{\tilde{T}}}\) :
-
Thermophoretic parameter \({\text{m}}^{2} \,{\text{s}}^{ - 1}\)
- \(k_{nf} ,k_{bf} ,k_{p}\) :
-
Nanofluid, basefluid, and nanoparticle thermal conductivity. \({\text{W}}\,{\text{m}}^{ - 1} \,{\text{K}}^{ - 1}\)
- \(d_{p}\) :
-
The diameter of nanoparticles \({\text{m}}\)
- \(\sigma^{*}\) :
-
Stefan-Boltzmann constant \({\text{W}}\,{\text{m}}^{ - 2} \,{\text{K}}^{ - 4}\)
- \(\kappa^{*}\) :
-
Mean absorption coefficient \({\text{m}}^{ - 1}\)
- \(T_{fr}\) :
-
The freezing temperature of the base fluid \({\text{K}}\)
- \(U_{{\tilde{E}}}\) :
-
Electric potential \({\text{V}}\)
- \(\varepsilon_{0}\) :
-
Vacuum’s dielectric constant \({\text{Fm}}^{ - 1}\)
- \(n_{0}\) :
-
The bulk concentration of ionic species \({\text{Kg}}\,{\text{m}}^{ - 3}\)
- \(e\) :
-
Charge of electron \({\text{C}}\)
- \(\Phi_{0}\) :
-
Nanoparticle volume fraction [–]
- \(\varepsilon_{r}\) :
-
Fluid’s relative permittivity [–]
- \(z\) :
-
Valence [–]
- \(p\) :
-
Dimensionless pressure parameter [–]
- \(u,w\) :
-
Dimensionless velocity components [–]
- \(k\) :
-
Debye length parameter [–]
- \(U_{0}\) :
-
Electroosmotic velocity [–]
- \(N_{t}\) :
-
Dimensionless thermophoretic parameter [–]
- \({\text{Gr}}_{\Phi }\) :
-
Mass Grashof number [–]
- \(R_{d}\) :
-
Dimensionless radiation parameter [–]
- \(\eta_{w}\) :
-
Temperature ratio parameter [–]
- \(N_{b}\) :
-
Dimensionless Brownian diffusion parameter [–]
- \(\Psi\) :
-
Dimensionless stream function [–]
- \(S\) :
-
Joule heating parameter [–]
- \({\text{Gr}}_{{\text{t}}}\) :
-
Temprature Grashof number [–]
- \({\text{Pr}}\) :
-
Prandtl number [–]
- \(\delta\) :
-
Wave number [–]
- \({\text{Re}}\) :
-
The Reynolds number [–]
- \(\beta\) :
-
Thermal slip parameter [–]
- \(\xi\) :
-
Dimensionless zeta potential [–]
- \({\text{Nu}}\) :
-
Nusselt number [–]
- \(Q\) :
-
The dimensionless averaged flux [–]
- θ:
-
Dimensionless temperature [–]
- Φ:
-
The dimensionless nanoparticle volume fraction [–]
References
Xuan, X., Xu, B., Sinton, D. & Li, D. Electroosmotic flow with Joule heating effects. Lab Chip 4(3), 230–236 (2004).
Yuan, S., Zhou, M., Liu, X. & Jiang, B. Electro-elastic instability in electroosmotic flows of viscoelastic fluids through a model porous system. Eur. J. Mech. B. Fluids 97, 173–186 (2023).
Yuan, S., Zhou, M., Liu, X. & Jiang, B. Effect of pressure-driven flow on electroosmotic flow and electrokinetic mass transport in microchannels. Int. J. Heat Mass Transf. 206, 123925 (2023).
Alyousef, H. A. et al. Mathematical modeling and analysis of the steady electro-osmotic flow of two immiscible fluids: A biomedical application. Coatings 13, 115 (2023).
Sleigh, M. A., Blake, J. R. & Liron, N. The propulsion of mucus by cilia. Am. Rev. Respir. Dis. 137(3), 726–741 (1988).
Brennen, C. & Winet, H. Fluid mechanics of propulsion by cilia and flagella. Annu. Rev. Fluid Mech. 9(1), 339–398 (1977).
Sade, J., Eliezer, N., Silberberg, A. & Nevo, A. C. The role of mucus in transport by cilia. Am. Rev. Respir. Dis. 102(1), 48–52 (1970).
Vilfan, A. & Jülicher, F. Hydrodynamic flow patterns and synchronization of beating cilia. Phys. Rev. Lett. 96(5), 058102 (2006).
Lardner, T. J. & Shack, W. J. Cilia transport. Bull. Math. Biophys. 34(3), 325–335 (1972).
Asadnia, M. et al. From biological cilia to artificial flow sensors: Biomimetic soft polymer nanosensors with high sensing performance. Sci. Rep. 6, 32955 (2016).
Mendes, R. U., de Castro, H. F., Cavalca, K. L. & Ferreira, L. O. S. Magnetic actuator modelling for rotating machinery analysis. In Vibration Problems ICOVP 2011 537–542 (Springer, Dordrecht, 2011).
Wang, Y., Gao, Y., Wyss, H. M., Anderson, P. D. & den Toonder, J. M. J. Artificial cilia fabricated using magnetic fiber drawing generate substantial fluid flow. Microfluid. Nanofluid. 18(2), 167–174 (2015).
Bhatti, M. M., Ishtiaq, F., Ellahi, R. & Sait, S. M. Novel aspects of cilia-driven flow of viscoelastic fluid through a Non-Darcy medium under the influence of an induced magnetic field and heat transfer. Mathematics 11(10), 2284 (2023).
Akbar, N. S., Tripathi, D., Khan, Z. H. & Bég, O. A. Mathematical modeling of pressure-driven micropolar biological flow due to metachronal wave propulsion of beating cilia. Math. Biosci. 301, 121–128 (2018).
Hall, J. & Clarke, N. The mechanics of cilium beating: quantifying the relationship between metachronal wavelength and fluid flow rate. J. Fluid Mech. 891, A20 (2020).
Ramesh, K., Tripathi, D. & Bég, O. A. Cilia-assisted hydromagnetic pumping of biorheological couple stress fluids. Propuls. Power Res. 8(3), 221–233 (2019).
Ishtiaq, F., Ellahi, R., Bhatti, M. M. & Alamri, S. Z. Insight in thermally radiative cilia-driven flow of electrically conducting non-Newtonian jeffrey fluid under the influence of induced magnetic field. Mathematics 10(12), 2007 (2022).
Asghar, Z., Javid, K., Waqas, M., Ghaffari, A. & Khan, W. A. Cilia-driven fluid flow in a curved channel: Effects of complex wave and porous medium. Fluid Dyn. Res. 52(1), 015514 (2020).
Wu, A. et al. A shear-rate-dependent flow generated via magnetically controlled metachronal motion of artificial cilia. Biomech. Model. Mechanobiol. 19, 1713–1724 (2020).
Akbar, N. S., Tripathi, D., Bég, O. A. & Khan, Z. H. MHD dissipative flow and heat transfer of Casson fluids due to metachronal wave propulsion of beating cilia with thermal and velocity slip effects under an oblique magnetic field. Acta Astronaut. 128, 1–12 (2016).
Sadaf, H. & Nadeem, S. Fluid flow analysis of cilia beating in a curved channel in the presence of magnetic field and heat transfer. Can. J. Phys. 98(2), 191–197 (2020).
Mekheimer, Kh. S. & Abd Elmaboud, Y. The influence of heat transfer and magnetic field on peristaltic transport of a Newtonian fluid in a vertical annulus: Application of an endoscope. Phys. Lett. A 372, 1657–1665 (2008).
Abdelsalam, S. I., Bhatti, M. M., Zeeshan, A., Riaz, A. & Bég, O. A. Metachronal propulsion of a magnetized particle-fluid suspension in a ciliated channel with heat and mass transfer. Phys. Scripta 94(11), 115301 (2019).
Abdelsalam, S. I. & Zaher, A. Z. On behavioral response of ciliated cervical canal on the development of electroosmotic forces in spermatic fluid. Math. Model. Nat. Phenom. 17, 27 (2022).
Imran, A., Raja, M. A. Z., Shoaib, M., Zeb, M. & Nisar, K. S. Electro-osmotic transport of a Williamson fluid within a ciliated microchannel with heat transfer analysis. Case Stud. Therm. Eng. 45, 102904 (2023).
Gul, F., Maqbool, K. & Mann, A. B. Thermal analysis of electroosmotic flow in a vertical ciliated tube with viscous dissipation and heat source effects. J. Therm. Anal. Calorim. 143, 2111–2123 (2021).
Ijaz, S., Abdullah, M., Sadaf, H. & Nadeem, S. Generalized complex cilia tip modeled flow through an electroosmotic region. J. Cent. South Univ. 30, 1217–1230 (2023).
Javid, K., Khan, R., Al-Khaled, K., Khan, M. I., Khan, S. U., Guedri, K., & Galal, A. M. EMHD ciliated pumping of viscoelastic liquid in a complex convergent/divergent channel. Waves Random Complex Media 1–23 (2022)
Heinonen, E. W., Wildin, M. W., Beall, A. N. & Tapscott, R. E. Assessment of antifreeze solutions for ground-source heat pump systems. ASHRAE Trans. 103(2), 747–756 (1997).
Shojaeizadeh, E., Veysi, F., Yousefi, T. & Davodi, F. An experimental investigation on the efficiency of a Flat-plate solar collector with binary working fluid: A case study of propylene glycol (PG)–water. Exp. Therm. Fluid Sci. 53, 218–226 (2014).
JuGer, J. J. & Crook, R. F. Heat transfer performance of propylene glycol versus ethylene glycol coolant solutions in laboratory testing. J. Eng. 108, 71–81 (1999).
Winkler, J. Titanium Dioxide (Vincentz, Hannover, 2003).
Weir, A., Westerhoff, P. & Fabricius, L. Kiril Hristovski and natalie von Goetz titanium dioxide nanoparticles in food and personal care products. Environ. Sci. Technol. 46, 2242–2250 (2012).
Carey, J. H., Lawrence, J. & Tosine, H. M. Photodechlorination of PCB’s in the presence of titanium dioxide in aqueous suspensions. Bull. Environ. Contam. Toxicol. 16, 697–701. https://doi.org/10.1007/BF01685575 (1976).
Kiser, M. A. et al. Titanium nanomaterial removal, and release from wastewater treatment plants. Environ. Sci. Technol. 43, 6757–6763 (2009).
Suh, W. H., Suslick, K. S., Stucky, G. D. & Suh, Y. H. Nanotechnology, nanotoxicology, and neuroscience. Prog. Neurobiol. 87(3), 133–170 (2009).
Choi, S. U. S. Enhancing thermal conductivity of fluids with nanoparticles, developments and application of non-Newtonian flows. ASME J. Heat Transf. 66, 99–105 (1995).
Javidan, M. & Moghadam, A. J. Effective cooling of a photovoltaic module using jet-impingement array and nanofluid coolant. Int. Commun. Heat Mass Transfer 137, 106310 (2022).
Saleem, S. et al. Insight into the motion of water conveying three kinds of nanoparticles shapes on a horizontal surface: Significance of thermo-migration and Brownian motion. Surf. Interfaces 30, 101854 (2022).
Elnaqeeb, T., Animasaun, I. L. & Shah, N. A. Ternary-hybrid nanofluids: Significance of suction and dual-stretching on three-dimensional flow of water conveying nanoparticles with various shapes and densities. Z. Naturforsch. 76, 231–243 (2021).
Animasaun, I. L., Oke, A. S., Al-Mdallal, Q. M. & Zidan, A. M. Exploration of water conveying carbon nanotubes, graphene, and copper nanoparticles on impermeable stagnant and moveable walls experiencing variable temperature: Thermal analysis. J. Therm. Anal. Calorim. 148, 4513–4522 (2023).
Rasool, G., Wakif, A., Wang, X., Shafiq, A. & Chamkha, A. J. Numerical passive control of alumina nanoparticles in purely aquatic medium featuring EMHD driven non-Darcian nanofluid flow over convective Riga surface. Alex. Eng. J. 68, 747–762 (2023).
Adnan, W. Ashraf, analysis of heat transfer performance for ternary nanofluid flow in radiated channel under different physical parameters using GFEM. J. Taiwan Inst. Chem. Eng. 146, 104887 (2023).
Xiu, W., Animasaun, I. L., Al-Mdallal, Q. M., Alzahrani, A. K. & Muhammad, T. Dynamics of ternary-hybrid nanofluids due to dual stretching on wedge surfaces when volume of nanoparticles is small and large: Forced convection of water at different temperatures. Int. Commun. Heat Mass Transfer 137, 106241 (2022).
Shah, Z., Roomana, M. & Shutaywi, M. Computational analysis of radiative engine oil-based Prandtl-Eyring hybrid nanofluid flow with variable heat transfer using the Cattaneo-Christov heat flux model. RSC Adv. 13, 3552–3560 (2023).
Akbar, N. S., Tripathi, D., Khan, Z. H. & Bég, O. A. Mathematical model for ciliary-induced transport in MHD flow of Cu–H2O nanofluids with magnetic induction. Chin. J. Phys. 55(3), 947–962 (2017).
Corcione, M. Empirical correlating equations for predicting the effective thermal conductivity and dynamic viscosity of nanofluids. Energy Convers. Manag. 52, 789–793. https://doi.org/10.1016/j.enconman.2010.06.072 (2011).
Wakif, A., Boulahia, Z., Ali, F., Eid, M. R. & Sehaqui, R. Numerical analysis of the unsteady natural convection MHD couette nanofluid flow in the presence of thermal radiation using single and two-phase nanofluid models for Cu–water nanofluids. Int. J. Appl. Comput. Math. 4, 81 (2018).
Tippa, S., Narahari, M. & Pendyala, R. Unsteady natural convection flow of nanofluids past a semi-infinite isothermal vertical plate. AIP Conf. Proc. 1787, 020014 (2016).
Author information
Authors and Affiliations
Contributions
J.A. contributed in Methodology, Software, Investigation, Solutions, Analysis and Writing- Original draft preparation. N.S.A. supervised, validated and reviewed the work.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Akram, J., Akbar, N.S. Electroosmotically actuated peristaltic-ciliary flow of propylene glycol + water conveying titania nanoparticles. Sci Rep 13, 11801 (2023). https://doi.org/10.1038/s41598-023-38820-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-38820-4
This article is cited by
-
Non-Newtonian rheology of blood in elliptical cross-section artery affected by several stenosis: Prandtl fluid model
Archive of Applied Mechanics (2024)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.