## Analytical and Numerical Approaches for Complicated Nonlinear Equations

View this Special IssueResearch Article | Open Access

# Effect of the Velocity Second Slip Boundary Condition on the Peristaltic Flow of Nanofluids in an Asymmetric Channel: Exact Solution

**Academic Editor:**Abdon Atangana

#### Abstract

The problem of peristaltic nanofluid flow in an asymmetric channel in the presence of the second-order slip boundary condition was investigated in this paper. To the best of the authors’ knowledge, this parameter was here incorporated for the first time in such field of a peristaltic flow. The system governing the current flow was found as a set of nonlinear partial differential equations in the stream function, pressure gradient, nanoparticle concentration, and temperature distribution. Therefore, this system has been successfully solved exactly via a very effective procedure. These exact solutions were then proved to reduce to well-known results in the absence of second slip which were published very recently in the literature. Effect of the second slip parameter on the present physical parameters was discussed through graphs and it was found that this type of slip is a very important one to predict the investigated physical model. Moreover, the variation of many physical parameters such as amplitudes of the lower and upper waves, phase difference on the temperature distribution, nanoparticle concentration, pressure rise, velocity, and pressure gradient were also discussed. Finally, the present results may be viewed as an optimal choice for their dependence on the exact solutions which are obtained due to the highly complex nonlinear system.

#### 1. Introduction

Peristalsis is an interest subject and has recently attracted much attention due to its importance in engineering and medical applications. In human body, peristalsis is found in the swallowing food through the esophagus, chyme motion in the gastrointestinal tract, vasomotion of small blood vessels such as venules, capillaries and arterioles, urine transport from kidney to bladder, and intrauterine fluid flow within the uterine cavity. Various experimental and theoretical studies have been conducted to understand peristaltic flow in asymmetric channel or an axisymmetric tube; see [1–7]. According to De Vries et al. [8], they observed that myometrial contractions are peristaltic-type motion and therefore these contractions of the uterine wall may occur in both symmetric and asymmetric directions. Due to the importance of this field of research, a great effort was devoted to study this type of flow for Newtonian and non-Newtonian fluids in an asymmetric channel [9–17].

Recently, a few papers have been published in the field of nanofluid flows under the peristaltic action [18–23]. Nanofluids usually lead to the enhancement of the thermal conductivity of the base fluid [24]. The concept of nanofluids, analogous to that of nanoparticles, may be initiated by Choi et al. [25]. The fact that nanofluids have higher thermal conductivity than the other heat transfer fluids because of their nanostructure has attracted many engineers and theoretical scientists to investigate their behavior. This higher thermal conductivity may be very useful in the treatment of tumors by injecting the blood vessel nearest to the tumor with magnetic nanoparticles along with placing a magnet close to the tumor, where these particles act like heat sources, in the presence of the applied magnetic field of alternating nature. Hence, the undesirable tissues (cancer's tissues) are destroyed when the temperature reaches 42–45°C [6]. Besides, the drug may be placed on the magnetic nanoparticles and is injected near the tumor. The drug is then absorbed by the tumor through a high gradient magnetic field, which is concentrated near the tumor center [26]. In fact, the drug absorption due to high concentration of magnetic particles increases and magnetic force prevents uniform drug distribution in circulatory system. This approach reduces the side effect and allows using high dose of anticancer drug [27].

Very recently, Ebaid and Aly [28] investigated the system of partial differential equations describing the peristaltic flow of a nanofluid with slip effect of the velocity, temperature, and concentration has been analytically solved. The obtained exact solutions have been applied to study effects of the slip parameter, thermophoresis, Brownian motion parameters, and many other parameters on the pressure rise, velocity profiles, temperature distribution, nanoparticle concentration, and pressure gradient. These exact solutions were proved to reduce to the results in the literature at special cases when . Moreover, on comparing the present results with those in [22], remarkable differences were noticed for behavior of the included previous physical phenomena.

The main feature of the previous studies on the nanofluid flow under peristaltic action [18–23] is the ignorance of the second slip effect, while this effect has been recently discussed by many authors; see for example [29–33]. These papers showed clearly that the second-order slip flow model is necessary to predict the flow characteristics accurately. Therefore, as an extension of these studies and for the first time, in this paper we aim to investigate effect of the second slip on the peristaltic flow of nanofluids in an asymmetric channel. Hence, we focus here on the problem discussed by Akbar et al. [22] and Ebaid and Aly [28] with slight differences in the boundary conditions.

#### 2. The Mathematical Model

Consider peristaltic transport of an incompressible Newtonian nanofluid in an asymmetric channel with flexible walls. The channel asymmetry is generated by propagation of waves on the channel walls traveling with different amplitudes and phases but with the same constant speed . In the Cartesian coordinates system of the fixed frame, the upper and lower walls and , respectively, are given by (see Figure 1)

where and are amplitude of the waves, is the wave length, and is the width of the channel. The phase difference varies in the range , where and correspond to symmetric channel with waves out of the phase and in the phase, respectively. Further, , , , , and have to satisfy the following condition [34]: with the following nondimensional phenomena [22]: Furthermore, in the moving frame of references , we have which have been used with the following nondimensional phenomena [22]: to obtain the present physical model. On taking into account (i) heat transfer along with nanoparticle phenomena, (ii) long wavelength, and low Reynolds number approximation, recently, Akbar et al. [22] found that the flow is governed by the following system of partial differential equations in nondimensional form: where , , , and are the stream function, temperature distribution, nanoparticle concentration, and pressure gradient, respectively. In addition, , , , and are the Brownian motion parameter, thermophoresis parameter, local temperature Grashof number, and nanoparticle Grashof number, respectively. The above system of PDEs. (6)–(9) has to be solved subject to the following boundary conditions: where , , , and represent the first-order slip parameter, second-order slip parameter, thermal slip parameter, and concentration slip parameter, respectively. It should be noted that the full formulation of the above model can be checked in [22], and there is no need to repeat it again here. However, as mentioned in Section 1 in the present work, we focused on obtaining the general closed form solution in the presence of the second-order slip parameter, as a very important one to predict the investigated physical model.

#### 3. The General Closed Form Solution

As mentioned by Ebaid and Aly [28], we have the following exact solutions for the temperature distribution and nanoparticle concentration: where where can be obtained from the following implicit algebraic equation: It can be seen from (14) that it is of complex structure to be solved exactly in its general form. However, it has an exact solution for when and given by Regarding the stream function, it can be solved to give where Applying the boundary conditions on the -equation, we obtain the following system:

where

On solving the system (18a)–(18d), we get To get the pressure gradient , we obtain from (7) that where The pressure rise in terms of the flow rate is given as follows: Therefore, the exact expression for the pressure rise at the center of the channel is given as

#### 4. Results and Discussion

This section is devoted for investing the exact solutions which have been obtained in the previous section for the temperature distribution, nanoparticle concentration, velocity, pressure gradient, and pressure rise. These solutions are expected to provide us with the correct physical effect of the second slip for the five investigated physical phenomena. It is observed from (11) that there is no effect for the second slip parameter on the solutions of the temperature distribution and nanoparticle concentration. Although (11) have been used by Ebaid and Aly [28] to discuss the effect of many parameters on the temperature distribution and nanoparticle concentration for comparing with the results in [22], many other parameters were ignored. These important parameters are investigated in the next subsections.

##### 4.1. The Temperature Distribution

Very recently in [28], Ebaid and Aly discussed the effect of , , and on the temperature distribution. Here, we aim to discuss the effect of many other parameters on this phenomenon. In Figures 2(a) and 2(b), effect of the amplitude on the temperature distribution is depicted at two different values for the Browning motion parameter and fixed values for the other investigated parameters. From these figures, it can be concluded that there are no remarkable influences for the variation of on in the first part of the channel, while an increase in occurs in the upper part with increasing . The variation of the amplitude on the temperature distribution is presented in Figures 3(a) and 3(b). These figures show that very small decrease in occurs as increases. This means that while the amplitude has remarkable effect on , the influence of the amplitude on may be ignored. In Figures 4(a) and 4(b), the variation of the phase difference on is introduced. Although the increase in between the waves of the channel increases the temperature distribution, this increase is very small. Instead of studying the variation of at separate values of Brownian motion parameter , we present in Figure 5(a) the variation of against a continuous range for the values of . This figure indicates a clear description for the variation of Brownian motion parameter on . In addition, Figure 5(b) shows other descriptions for the variation of the thermophoresis parameter on the temperature distribution. It can be concluded from the last two figures that always decreases with the increase in and .

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

##### 4.2. The Nanoparticle Concentration

The nanoparticles concentration is depicted in Figures 6(a) and 6(b). It is observed from these figures that reaches its highest value at the lower wall of the channel and then it decreases in a certain domain across the channel. A converse of this behavior occurs after that domain. It is also detected that at a higher value of Brownian motion parameter the domain in which decreases becomes wider. In addition, the wide of such domain increases with increasing the amplitude of the upper wave. The variation of the lower wave amplitude on is presented in Figures 7(a) and 7(b). For all values of , it is noticed that decreases until it reaches a certain point. After that point the curves become identical whatever the value of becomes. This refers to that there is no effect for on after this point.

**(a)**

**(b)**

**(a)**

**(b)**

In Figures 8(a) and 8(b), the variation of the phase difference on is displayed. Although the increase in the phase difference between the waves of the channel increases the nanoparticles concentration, this occurs in certain domain of the channel. However, after that domain the curves are identical whatever the value of the phase difference. Moreover, instead of studying the variation of at separate values of Brownian motion parameter , the variation of against a continuous range for the values of is displayed in Figure 9(a). In addition, Figure 9(b) shows other descriptions for the variation of the thermophoresis parameter on the nanoparticle concentration. It can be concluded from Figure 9(a) that always decreases when is in the range from 1 to 10 at fixed values for the other parameters. However, the situation is different for the variation of against at certain values of , where it is noticed that when then decreases within the channel. In addition, at values of higher than 4, we observe that decreases in certain domain of the channel and changes this behavior in the rest of channel.

**(a)**

**(b)**

**(a)**

**(b)**

##### 4.3. The Pressure Rise

The exact expression for the pressure rise is given by (24) which is valid for any set of the physical parameters. On using this equation, we observe that the pressure rise is always a decreasing function in terms of the flow rate. This notice is clarified in Figures 10(a)–10(d). It is also observed from Figure 10(a) that the pressure rise increases with increasing the negativity of the second slip parameter in the whole range of the flow rate . The effect of the thermophoresis parameter on the pressure rise is depicted in Figure 10(b). It can be seen from this figure that the pressure rise decreases with the increase in the values of in the whole region of .

**(a)**

**(b)**

**(c)**

**(d)**

Figure 10(c) shows that the pressure rise decreases with decreasing the amplitude when lies in the region , while a converse behavior occurs for . Figure 10(d) reveals that the pressure rise increases with increasing the local nanoparticle Grashof number . This behavior differs from the results obtained by Akbar et al. [22] in the absence of the second slip effect. Therefore, the present results reveal that the existence of this parameter has a significant effect on the pressure rise and leads to different behavior when compared with its ignorance.

##### 4.4. The Velocity

Regarding the axial velocity, it was plotted in Figures 11(a)–11(d) at several values of the present physical parameters. Variations of the second slip parameter , Brownian motion parameter , flow rate , and local temperature Grashof number on the velocity profile have been displayed in Figures 11(a)–11(d), respectively. It is found from Figure 11(a) that the velocity increases near the lower wall with increasing , while a converse behavior occurs near the upper wall; that is, the velocity decreases with increasing . Unlike the study made by Akbar et al. [22], the effect of Brownian motion parameter leads to that there are three different regions in which alternating behavior for the velocity occurs. As shown in Figure 11(b), in a central region of the channel, that is, , the velocity increases with the increase in , while in the other two regions of the channel, that is, and , a different situation is detected. At the same three regions just mentioned, it can be seen from Figure 11(d) that effect of the local temperature Grashof number on the velocity profiles happens in a converse manner regarding Brownian motion parameter . However, it is clear from Figure 11(c) that an increase in the flow rate leads to an increase in the velocity and this is similar to the variation of the velocity profile against flow rate as in Akbar et al. [22].

**(a)**

**(b)**

**(c)**

**(d)**

##### 4.5. The Pressure Gradient

Figures 12(a)–12(d) indicate the pressure gradient for different values of , , , and , respectively. From these figures, it is shown that magnitude of the pressure gradient increases in view of an increase in , , and , and it decreases with an increase in . In addition, it is observed that our results on the effect of , , and on the pressure gradient are completely different comparing with those obtained in [22] for the effect of the first slip parameter only. This also may refer to the importance of including the second slip parameter in such problems of peristaltic flow. However, the effect of the amplitude on the pressure gradient remains the same even in the presence of , where Figure 12(d) shows that an increase in the amplitude leads to an increase in the pressure gradient. Besides, one of the main features of these graphs is that the maximum pressure gradient occurs at the same point when . In the absence of , this value was shifted slightly to the left; that is, it occurred at as in [22]. The final notice here is that the present numerical discussion is based on the obtained exact analytical solutions for the current physical model, unlike the approximate series solutions obtained in [22] whose disadvantages have been proved by Ebaid and Aly [28].

**(a)**

**(b)**

**(c)**

**(d)**

#### 5. Conclusion

The effect of the second slip condition on the peristaltic transport of a nanofluid in an asymmetric channel was investigated for the first time in the present paper. The system describing the problem was exactly solved via a simple, but very effective, analytical procedure. Based on the obtained exact solutions, some important results were introduced through graphs to indicate effect of the second slip on the velocity profiles, pressure rise, and pressure gradient. Moreover, the variation of many physical parameters such as amplitudes of the lower and upper waves, phase difference on the temperature distribution, and nanoparticle concentration were also discussed for a new set of the presented parameters rather than those in the literature. It was found that the temperature distribution always decreases with the increase in and . In addition, the amplitude has remarkable effect on , while the influence of the amplitude on may be ignored. For all values of , nanoparticles concentration decreases until it reaches a certain point; however, there is no effect for after this point. In a certain domain of the channel, the phase difference between the waves of the channel increases ; however, the curves are identical whatever value of the phase difference after that. Finally, it should be mentioned that the second slip condition had a significant effect on the various phenomena involved in the current model. Accordingly, the present suggested technique may be useful for obtaining exact solutions for many other similar problems. In addition, such type of the second slip should be included for modelling the peristaltic flow.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. (130-023-D1434). The authors, therefore, acknowledge DSR technical and financial support.

#### References

- T. W. Latham,
*Fluid motion in a peristaltic pump, [M.Sc. thesis]*, MIT, Cambridge, Mass, USA, 1966. - A. Ramachandra Rao and S. Usha, “Peristaltic transport of two immiscible viscous fluids in a circular tube,”
*Journal of Fluid Mechanics*, vol. 298, pp. 271–285, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH - J. C. Misra and S. K. Pandey, “Peristaltic transport in a tapered tube,”
*Mathematical and Computer Modelling*, vol. 22, no. 8, pp. 137–151, 1995. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - K. S. Mekheimer, “Peristaltic transport of a couple stress fluid in a uniform and non-uniform channels,”
*Biorheology*, vol. 39, no. 6, pp. 755–765, 2002. View at: Google Scholar - K. Vajravelu, S. Sreenadh, and V. R. Babu, “Peristaltic transport of a Herschel-Bulkley fluid in an inclined tube,”
*International Journal of Non-Linear Mechanics*, vol. 40, no. 1, pp. 83–90, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. S. Mekheimer and Y. Abd elmaboud, “The influence of heat transfer and magnetic field on peristaltic transport of a Newtonian fluid in a vertical annulus: application of an endoscope,”
*Physics Letters A*, vol. 372, no. 10, pp. 1657–1665, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - K. S. Mekheimer and Y. Abd elmaboud, “Peristaltic flow of a couple stress fluid in an annulus: application of an endoscope,”
*Physica A*, vol. 387, no. 11, pp. 2403–2415, 2008. View at: Publisher Site | Google Scholar - K. De Vries, E. A. Lyons, G. Ballard, C. S. Levi, and D. J. Lindsay, “Contractions of the inner third of the myometrium,”
*American Journal of Obstetrics and Gynecology*, vol. 162, no. 3, pp. 679–682, 1990. View at: Publisher Site | Google Scholar - G. C. Shit, M. Roy, and E. Y. K. Ng, “Effect of induced magnetic field on peristaltic flow of a micropolar fluid in an asymmetric channel,”
*International Journal for Numerical Methods in Biomedical Engineering*, vol. 26, no. 11, pp. 1380–1403, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. Yildirim and S. A. Sezer, “Effects of partial slip on the peristaltic flow of a MHD Newtonian fluid in an asymmetric channel,”
*Mathematical and Computer Modelling*, vol. 52, no. 3-4, pp. 618–625, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - K. S. Mekheimer, S. Z. A. Husseny, and Y. Abd Elmaboud, “Effects of heat transfer and space porosity on peristaltic flow in a vertical asymmetric channel,”
*Numerical Methods for Partial Differential Equations*, vol. 26, no. 4, pp. 747–770, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - K. S. Mekheimer, S. Z.-A. Husseny, and A. I. Abd El Lateef, “Effect of lateral walls on peristaltic flow through an asymmetric rectangular duct,”
*Applied Bionics and Biomechanics*, vol. 8, no. 3-4, pp. 295–308, 2011. View at: Publisher Site | Google Scholar - S. Srinivas and R. Muthuraj, “Effects of chemical reaction and space porosity on MHD mixed convective flow in a vertical asymmetric channel with peristalsis,”
*Mathematical and Computer Modelling*, vol. 54, no. 5-6, pp. 1213–1227, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. K. Pandey and M. K. Chaube, “Peristaltic transport of a maxwell fluid in a channel of varying cross section induced by asymmetric waves: application to embryo transport within uterine cavity,”
*Journal of Mechanics in Medicine and Biology*, vol. 11, no. 3, pp. 675–690, 2011. View at: Publisher Site | Google Scholar - S. Srinivas, R. Gayathri, and M. Kothandapani, “Mixed convective heat and mass transfer in an asymmetric channel with peristalsis,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 4, pp. 1845–1862, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - K. Das, “Influence of slip and heat transfer on mhd peristaltic flow of a Jeffrey fluid in an inclined asymmetric porous channel,”
*Indian Journal of Mathematics*, vol. 54, pp. 19–45, 2012. View at: Google Scholar | MathSciNet - Y. Abd Elmaboud, S. Kh. Mekheimer, and A. I. Abdellateef, “Thermal properties of couple-stress fluid flow in an asymmetric channel with peristalsis,”
*Journal of Heat Transfer*, vol. 135, Article ID 044502, 2013. View at: Google Scholar - N. S. Akbar and S. Nadeem, “Endoscopic effects on peristaltic flow of a nanofluid,”
*Communications in Theoretical Physics*, vol. 56, no. 4, pp. 761–768, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH - N. S. Akbar, S. Nadeem, T. Hayat, and A. A. Hendi, “Peristaltic flow of a nanofluid in a non-uniform tube,”
*Heat and Mass Transfer*, vol. 48, no. 3, pp. 451–459, 2012. View at: Publisher Site | Google Scholar - N. S. Akbar and S. Nadeem, “Peristaltic flow of a Phan-Thien-Tanner nanofluid in a diverging tube,”
*Heat Transfer*, vol. 41, no. 1, pp. 10–22, 2012. View at: Publisher Site | Google Scholar - M. Mustafa, S. Hina, T. Hayat, and A. Alsaedi, “Influence of wall properties on the peristaltic flow of a nanofluid: analytic and numerical solutions,”
*International Journal of Heat and Mass Transfer*, vol. 55, no. 17-18, pp. 4871–4877, 2012. View at: Publisher Site | Google Scholar - N. S. Akbar, S. Nadeem, T. Hayat, and A. A. Hendi, “Peristaltic flow of a nanofluid with slip effects,”
*Meccanica*, vol. 47, no. 5, pp. 1283–1294, 2012. View at: Publisher Site | Google Scholar | MathSciNet - O. A. Bég and D. Tripathi, “Mathematica simulation of peristaltic pumping with double-diffusive convection in nanofluids: a bio-nano-engineering model,”
*Journal of Nanoengineering and Nanosystems*, vol. 225, no. 3, pp. 99–114, 2012. View at: Publisher Site | Google Scholar - S. U. S. Choi, “Enhancing thermal conductivity of fluids with nanoparticles,” in
*Proceedings of the ASME International Mechanical Engineering Congress and Exposition*, vol. 66 of*FED 231/MD*, p. 99, ASME, San Francisco, Calif, USA, 1995. View at: Google Scholar - S. U. S. Choi, Z. G. Zhang, W. Yu, F. E. Lockwood, and E. A. Grulke, “Anomalous thermal conductivity enhancement in nanotube suspensions,”
*Applied Physics Letters*, vol. 79, no. 14, pp. 2252–2254, 2001. View at: Publisher Site | Google Scholar - Q. A. Pankhurst, J. Connolly, S. K. Jones, and J. Dobson, “Applications of magnetic nanoparticles in biomedicine,”
*Journal of Physics D*, vol. 36, no. 13, pp. R167–R181, 2003. View at: Publisher Site | Google Scholar - M. R. Habibi, M. Ghassemi, and M. H. Hamedi, “Analysis of high gradient magnetic field effects on distribution of nanoparticles injected into pulsatile blood stream,”
*Journal of Magnetism and Magnetic Materials*, vol. 324, no. 8, pp. 1473–1482, 2012. View at: Publisher Site | Google Scholar - A. Ebaid and E. H. Aly, “Exact analytical solution of the peristaltic nanofluids flow in an asymmetric channel with flexible walls and slip condition: application to the cancer treatment,”
*Computational and Mathematical Methods in Medicine*, vol. 2013, Article ID 825376, 8 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - T. Fang, S. Yao, J. Zhang, and A. Aziz, “Viscous flow over a shrinking sheet with a second order slip flow model,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 15, no. 7, pp. 1831–1842, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. M. Nandeppanavar, K. Vajravelu, M. S. Abel, and M. N. Siddalingappa, “Second order slip flow and heat transfer over a stretching sheet with non-linear Navier boundary condition,”
*International Journal of Thermal Sciences*, vol. 58, pp. 143–150, 2012. View at: Publisher Site | Google Scholar - M. Turkyilmazoglu, “Heat and mass transfer of MHD second order slip flow,”
*Computers & Fluids*, vol. 71, pp. 426–434, 2013. View at: Publisher Site | Google Scholar - A. V. Rosca and I. Pop, “Flow and heat transfer over a vertical permeable stretching/shrinking sheet with a second order slip,”
*International Journal of Heat and Mass Transfer*, vol. 60, pp. 355–364, 2013. View at: Publisher Site | Google Scholar - E. H. Aly and A. Ebaid, “Exact analytical solution for suction and injection flow with thermal enhancement of five nanofluids over an isothermal stretching sheet with effect of the slip model: a comparative study,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 721578, 14 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - M. Mishra and A. Ramachandra Rao, “Peristaltic transport of a Newtonian fluid in an asymmetric channel,”
*Zeitschrift fur Angewandte Mathematik und Physik*, vol. 54, no. 3, pp. 532–550, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2014 Emad H. Aly and Abdelhalim Ebaid. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.