Abstract
Cilia-induced flow of fractional Burgers fluid is studied in an inclined tube for both symplectic and antiplectic wave patterns. The solution of the problem is persued under the long wave length limitation. The fractional Adomian decomposition method is employed to evaluate the pressure gradient. Mathematical expressions for the axial velocity, frictional force, pressure gradient, and stream function are obtained and the influence of the main operating parameters is discussed in detail. It is noted that the velocity profile is more dominant in the case of antiplectic metachronal waves compared to symplectic ones, which confirms former results on the better capability of antiplectic waves to transport mucus, obtained with more complex numerical solvers.
1. Introduction
Dutch light microscopist Antoni was the first to discover cilia in 1675 and Sharpey was the first to discuss cilia in English language in 1835. Exhaustive studies on the ciliary structures have been carried out during the nineteenth century (Sleigh et al., ). Cilia and flagella oscillate in a waving fashion during motion to transport fluids and propel cells (Vélez-Cordero and Lauga, ). Cilia motion constitutes a pivotal role in a wide variety of physiological processes, such as alimentation, circulation, locomotion, respiration, and reproduction (Maiti and Pandey, ). Ciliated cells can indeed be found in many human organs e.g., photoreceptor cells in retina, in hair bundles on ear, epithelial cells in the respiratory tract, in Fallopian tubes (Ashraf et al., ), in kidney (Guirao and Joanny, ), in the ependymal cells of the brain that generate cerebrospinal flow among few examples (Lieberstein, ). Malfunctioning of the ciliary activity may be responsible of many respiratory diseases (Rubin, ), like severe asthma.
As the amalgamated motion of cilia occurs, the cilia's upper layer can be seen to have a metachronal wave generated as a result of a small phase lag between neighboring cilia. This collective motion of cilia supports many physiological processes (Lieberstein, ; Murakami and Takahashi, ; Takahashi and Shingyoji, ). Different types of metachronal waves are classified depending upon the dynamics and strokes of the moving cilia. Symplectic beat patterns are produced if the directions of the propagative metachronal waves and the main flow are the same and antiplectic patterns are recognized when the directions of the wave propagation and the flow are opposite (Knight-Jones, ; Blake, ).
To date, many efforts have been made mathematically to understand the involvement of ciliary motion related to different biofluids, such as bronchial mucus, semen …in humans which are represented by Newtonian and non-Newtonian fluid models. Maiti and Pandey () applied a numerical approach to study the flow of a power-law fluid representing semen in an axisymmetric tube (representing the efferent ducts of the male human reproductive channel) due to cilia motion and concluded that, more than only the ciliary activity, there are several other factors like the smoothness of muscles, the constant fluid secretion or the vacuum created during ejaculation, which may be responsible for the semen flow. Siddiqui et al. (, ) obtained exact solutions of cilia induced flow problems of viscous fluid and power-law fluid models representing semen in a cylindrical tube and infinite channel under the long wavelength approximation.
Mucus plays a significant role in individual's health, therefore many researchers discussed the bronchial mucus transportation due to ciliary activity (see for examples in Barton and Raynor, ; Ross and Corrsin, ; Fulford and Blake, ; Maqbool et al., ). Norton et al. () developed a transport model, where the mucus, considered successively as a Doi-Edwards, Jeffrey, and Maxwell fluid, is transported as a rigid body and the metachronal wave exhibits a symplectic behavior. Vélez-Cordero and Lauga () applied the regular perturbation method to solve a problem in which the tracheobronchial mucus is considered as a Carreau fluid. This solution only portrays the Newtonian effects when second-order perturbations are considered and non-Newtonian effects are captured when the perturbation analysis is pushed up to the fourth-order. Maqbool et al. () considered the geometry of Siddiqui et al. () and studied the mucus flow treating mucus as a Jeffrey fluid. Smith et al. () noted that the most advanced model to investigate the mucociliary clearance process is the Maxwell viscoelastic model and proposed a fluid-structure interaction model to examine the complex fluid flow problem arising due to ciliary activity.
More recently, the immersed boundary (IB) method (Hao and Zhu, ) has been extensively applied to study biofluid flows representative of humans/animals. Dillon et al. () used the IB method to examine the two dimensional flow of three cilia in a mucus layer such that the mucus layer is treated as an elastic solid instead of a viscoelastic fluid. Dauptain et al. () used the IB method to examine the motion of fluid due to one row of cilia on a ctenophore Pleurobrachia pileus, which is commonly known as a sea gooseberry for Reynolds numbers Re within the range [50−200]. They found that as the beating of cilia increases, it spreads more power to the interacting fluid and this work may be considered as a guideline for solving the fluid-structure interaction problem. Very recently, Chatelin and Poncet () investigated by 3D simulations the influence of the mucus viscosity, fluid height, cilia length, and beating frequency on the mucociliary process in a two-phase environment. Chateau et al. () performed 3D simulations of the transport and mixing induced by beating cilia at Re up to 20 in a two-phase environment composed of Newtonian fluids using a coupled IB/lattice Boltzmann method.
All reported studies (Barton and Raynor, ; Ross and Corrsin, ; Fulford and Blake, ; Dillon et al., ; Smith et al., ; Vélez-Cordero and Lauga, ; Maqbool et al., ) confirmed the fact that more efforts are still required in terms of scientific research to better understand the internal flow structure due to cilia motion and their interaction with the surrounding fluid. In the present work, one will focus on the mucociliary clearance process. It is well-known now that bronchial mucus exhibits complex rheological properties: stress relaxation, tensile stresses, shear thinning, yielding stress, and thixotropic behavior (see in Lafforgue et al., , ). Though the advanced numerical solvers developed by Chatelin and Poncet () and Chateau et al. () considered the two-phase character of the problem and the behavior of each single cilium, such approaches do not account for the rheology of mucus. Moreover, simulating the mucociliary clearance process in human airways remains very challenging for such methods due to the multiscale character of the problem: from the micrometer scale when considering each individual cilium to the decimeter scale when looking at the main air flow within the trachea. So analytical approaches or simplified models like the envelope model are still deemed necessary if one wants to simulate the complete problem. The present paper is an attempt to demonstrate that analytical solutions obtained by the Adomian decomposition method can provide useful informations regarding the mucociliary clearance process. Momani and Odibat () used successfully the Adomian decomposition method to solve a time-fractional Navier–Stokes equation in a tube and demonstrated both the efficiency and simplicity of its method. In this paper, bronchial mucus is considered as a fractional Burgers fluid. Motion is generated by linear pressure produced by the tips of the moving cilia under the long wavelength and low Reynolds number approximations (Shapiro et al., ). Various illustrations highlighting the effects of the most important parameters are also sketched. Another motivation of the present paper is that the literature is scarce on fractional fluid models (see the monograph of Oldham and Spanier, for example).
2. Mathematical Model
The fluid motion characteristics of an incompressible fractional Burgers fluid in a ciliary tube having an inclination angle θ are considered. The metachronal wave and inclined tube move with the same speed c to the right as shown in Figure 1.
Figure 1
The cilia tips follow an elliptical path as suggested by Maqbool et al. () and Siddiqui et al. (, ), which can be represented by
where Equations (1, 2) are the parametric equations representing the cilia motion in which a is the mean tube radius, c is the wave speed, is time, Z0 is the reference position of cilia, α* is the eccentricity of ellipse, ϵ is a dimensionless varying parameter and η is the wavelength of the metachronal wave. One assumes that the fluid and adjacent cilia tips have the same velocity resulting in a no-slip condition. Therefore, the axial and radial velocities can be written as
and
Introducing Equations (1, 2) into Equations (3, 4) yields
and
where and are the radial and axial velocity components. The wave and fixed frames are related through the transformations
and
where , and are the pressure, radius, and axial position in the fixed (wave) frame. The extra stress tensor for fractional Burgers' fluid is given as
where Γ is the shear stress tensor, is the strain rate, λi(i = 1, 2, 3) are constitutive parameters, α and β are the fractional derivative and integral defined as Hilfer ()
and
with 0 ≤ α ≤ 1. The fractional Burgers fluid model reduces to fractional Oldroyd B model for λ2 = 0 and the classical viscous model can be obtained if λi = 0(i = 1, 2, 3).
The velocity components for the fractional Burgers' fluid model for an inclined tubular flow should satisfy the following equations
where g is the gravitational acceleration, and are the radial and axial velocity components, τij are the shear stress tensor components and ρ is the fluid density. To solve the problem, the following non-dimensional parameters are introduced
where η, a, c, and μ denote the wavelength, tube radius, wave speed, and dynamic viscosity, respectively, Re and Fr are the Reynolds and Froude numbers while β* is the wave number. Under the long wavelength and low Reynolds number approximations, the flow may be considered as a Stokes flow. Thus, during the non-dimensionalizing of Equations (12, 13) with the help of Equation (15), we have ignored the terms involving β*, β*2, β*3, and Re, Re2, Re3 …but terms involving Re/Fr are retained as the orders of Re and Fr numbers are the same. Equations (8–14) (after dropping hats) simplify to the following non-dimensional forms
with appropriate boundary conditions given as
where
The following two equations are used to determine the initial guesses required by the Adomian decomposition method
Integrating Equation (17) with respect to r and using the boundary condition (Equation 21), one gets
Further integrating Equation (24) and applying the boundary conditions (Equations 18, 19) will yield
The results for a Newtonian fluid in a horizontal tube can be deduced using the limits λi → 0 for i = 1, 2, 3 and θ → 0, therefore one gets
The volume flow rate is defined as
which, in light of Equation (24), becomes
so that q the dimensional volume rate and Q the dimensionless volume flow rate in the fixed frame are related as
The mean volume flow rate can be calculated using the time period T in Equation (30)
Equation (26) in view of Equation (29) gives
The stream function ψ in the wave frame is computed with the help of Equations (25, 28, 31) as
3. Solution Methodology
Equation (31) can be simplified as
where with the initial conditions
and
so that Equation (32) will take the form
where
An infinite series solution for l(z, t) by using the Adomian decomposition method (Adomian, ; Babolian and Biazar, ) is given by
where l0, l1, l2, l3, …ln+1 are determined as
From ln(z, t) (n ≽ 0), the other components can also be obtained. Finally an approximate solution of Equation (37) by truncating the series can be written as
where
The pressure difference Δp and friction force F across one wavelength are given by
4. Results and Discussion
In order to evaluate the pressure rise Δp, pressure gradient dp/dz, frictional force F, and streamlines ψ, the software Mathematica 8.0 has been used. The flow characteristics of Burgers fluids in an inclined ciliated tube are presented by controlling the fractional parameters α and β, the cilia length ϵ, and the angle of inclination θ. The Reynolds number Re and the Froude number Fr have been both fixed to 0.1. The Reynolds number based on the cilia tip speed is usually around 10−5 in human airways but Chateau et al. () recently demonstrated that there is no significant influence of Re as long as it remains lower than 1.
Table 1 depicts that the flow of the fractional Burgers fluid model gives the larger magnitude of the frictional force: typically F is 16.5 and 38.9 times higher for the fractional Burgers fluid model compared to the Burgers and the fractional Oldroyd-B models, respectively. The pressure difference Δp for the Oldroyd-B model is greater than for the fractional Burgers model, meaning that this last fluid is more ease to be transported by the mucociliary clearance process.
Table 1
| Fluid models | Rheological properties | Δp | dp/dz | F |
|---|---|---|---|---|
| Fractional burgers | λ1 = λ2 = λ3 = 1, α = β = 0.5 | 3.62293 | −1.9542 | 3.46156 |
| Burgers | λ1 = λ2 = λ3 = 1, α = β = 1 | 7.90151 | −0.624569 | 0.209561 |
| Fractional Oldroyd-B | λ1 = λ3 = 1, λ2 = 0, α = β = 0.5 | 7.75994 | 0.251982 | 0.0889968 |
| Oldroyd-B | λ1 = λ3 = 1, λ2 = 0, α = β = 1 | −7.27648 | 0.794008 | −0.814099 |
Pressure rise, pressure gradient, and frictional force for the four rheological models considered here.
Results obtained for r = 0.1 and z = 0.5 at t = 0.1.
The variations of Δp with are examined in Figure 2 for different values of α, β, ϵ, and θ. It is found that the pressure rise Δp decreases with an increase of the fractional parameter α, cilia length ϵ in the pumping expanse region (Δp > 0) of the tube and an opposite trend is noted in the copumping expanse region (Δp < 0) of the tube. Also in the pumping (resp. copumping) expanse region, the pressure rise increases (resp. decreases) for increasing values of the fractional parameter β. As expected increasing the cilia length or decreasing the tube radius provide similar results in terms of pressure variations but the effect is more pronounced by changing the tube radius. The ratio β/ϵ may be the relevant geometric parameter governing the transport of fluid by the mucociliary clearance process. The influence of the inclination angle θ on Δp is more straightforward as shown in Figure 2D. Δp increases by increasing θ in the pumping expanse region and an opposite trend is reported in the copumping expanse region. Larger pressure differences are observed for a vertical tube.
Figure 2
The variations of the pressure gradient are examined in Figure 3 for different values of α, β, ϵ, and θ. It is noted that remains small and the fluid can flow smoothly without the application of a large pressure gradient in the expanse regions 0 ≤ z ≤ 0.2 and 0.8 ≤ z ≤ 1. On the other hand, for 0.2 < z < 0.8, a large amount of pressure gradient is required to maintain the flow. The magnitude of the pressure gradient increases by increasing the parameters β, ϵ, and θ. These parameters provide the resistive force to the flow thus a larger value of is required to maintain the fluid flow, whereas the parameter α provides the deriving force to the flow thus a smaller value of is required to maintain the fluid flow.
Figure 3
It is observed through Figure 4 that the frictional force F varies linearly by increasing α, β, and ϵ and it tends to increase in magnitude by increasing β and ϵ. The magnitude of the resistive force decreases by increasing α, in the expanse region . On the contrary, its magnitude increases in the same region when increasing β and/or ϵ. The magnitude of the frictional force increases in the expanse by increasing the tube inclination θ but an opposite trend is noted for . The inclination of the tube provides then a large magnitude of F to oppose the flow.
Figure 4
Figure 5 displays the streamline patterns and trapping for α* = β* = 0.7 and by increasing the dimensionless cilia length ϵ. The center line symmetry bifurcates the boluses of the fluid particles circulating along the closed stream lines. The boluses are confined and move with the velocity of the metachronal waves. In Figures 5A–C, keeping all the other parameters constant, the number and size of the boluses increase by increasing ϵ. Thus, the cilia length significantly affects the flow dynamics by generating boluses in the inclined tube.
Figure 5
Figure 6 displays a 3D view of the axial and radial velocity components for ϵ = 0.2, α = 0.7, β = 0.7, and . It is observed that the waves are elliptic waves and the fluid velocity is in both forward and backward directions with a symmetric behavior due to metachronism of the ciliary motion.
Figure 6
Figure 7 displays the radial distributions of the axial and radial velocity components for four flow rates and both antiplectic and symplectic waves formed by the cilia tips. As expected, both components increase by increasing the flow rate. More interestingly, at any value of , the axial velocity produced by antiplectic waves is larger than the one due to symplectic waves. The increase in the velocity profiles is also faster for antiplectic waves compared to symplectic waves. Regarding the radial velocity component, the two types of metachronal waves provide symmetric profiles in the radial direction. All in all, it confirms the former results of Chateau et al. () obtained using a more advanced numerical model but for a two Newtonian fluid configuration. For a single layer of fractional Burgers fluid, antiplectic waves are also more efficient to transport fluid than symplectic waves.
Figure 7
5. Conclusions
In this paper, fractional Burgers' fluid flow in an inclined ciliated tube is examined. Using the long wavelength approximation, a semi-analytic solution is developed. Frictional force, pressure rise, pressure gradient and streamlines are plotted for different values of the main operating parameters and the main results can be summarized as follows:
Δp decreases by increasing α and ϵ in the pumping region and an opposite trend is noted in the copumping region of the tube.
Δp increases in the pumping region and decreases in the copumping region when increasing β.
is insignificant for 0 ≤ z ≤ 0.2 and 0.8 ≤ z ≤ 1. On the contrary, for 0.2 < z < 0.8, a higher value of is required to maintain the flux.
The magnitude of increases by increasing β, ϵ, and θ.
It is noticed that the parameters β, ϵ, and θ greatly influence the pressure gradient whereas the parameter α provides a smaller amount of the pressure gradient to the fluid flow.
The frictional force varies linearly when increasing α, β, and ϵ and its magnitude increases by increasing β and ϵ but it decreases in magnitude by increasing α and θ in the pumping region.
The magnitude of the frictional force increases in the region by increasing the tube inclination θ but an opposite trend is noted in the region .
The magnitude of the pressure difference is larger for fractional generalized Burgers model in comparison to the generalized Burgers model for and a reverse trend is observed for .
A greater magnitude of the pressure gradient is needed to the fluid flow for fractional generalized Burgers fluid as compared with the generalized Burgers fluid to pump the same amount of fluid.
The fractional generalized Burgers fluid provides a greater amount of frictional force in comparison to the generalized Burgers fluid.
The number and size of boluses increase by increasing ϵ.
The magnitude of the velocity profile increases significantly for antiplectic waves as compared to symplectic waves confirming that they are more efficient to transport mucus.
The fractional Adomian decomposition method appears then as a valuable tool to get useful and realistic results for the mucociliary transport process. Such method provides much faster results compared to more complex solvers (Chatelin and Poncet, ; Chateau et al., ) and could be associated to 3D flow solvers to solve the multiscale problem of the mucus clearance in the airways.
Statements
Author contributions
This work is the product of intellectual effort of all authors. They read and approved the final draft of the manuscript. KM and AM suggested and formulated the problem. KM, SS, and AM solved the problem jointly. SS wrote the paper. SP helped in data analysis and reviewed the manuscript.
Funding
SP acknowledges the Natural Sciences and Engineering Research Council of Canada for its financial support through a discovery grant (RGPIN-2015-06512).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AdomianG. (1994). Solving Frontier Problems of Physics: The Decomposition Method. Dordrecht: Kluwer Academic Publishers.
2
AshrafH.SiddiquiA. M.RanaM. A. (2018). Fallopian tube assessment of the peristaltic-ciliary flow of a linearly viscous fluid in a finite narrow tube. Appl. Math. Mech.39, 437–454. 10.1007/s10483-018-2305-9
3
BabolianE.BiazarJ. (2002). Solution of nonlinear equations by modified Adomian decomposition method. Appl. Math. Comput.132, 167–172. 10.1016/S0096-3003(01)00184-9
4
BartonC.RaynorS. (1967). Analytical investigation of cilia-induced mucous flow. Bull. Math. Biophys.29, 419–428.
5
BlakeJ. R. (1972). A model for the micro-structure in ciliated organisms. J. Fluid Mech.55, 1–23.
6
ChateauS.D'OrtonaU.PoncetS.FavierJ. (2018). Transport and mixing induced by beating cilia in human airways. Front. Physiol.9:161. 10.3389/fphys.2018.00161
7
ChatelinR.PoncetP. (2016). A parametric study of mucociliary transport by numerical simulations of 3D non-homogeneous mucus. J. Biomech.49, 1772–1780. 10.1016/j.jbiomech.2016.04.009
8
DauptainA.FavierJ.BottaroA. (2008). Hydrodynamics of ciliary propulsion. J. Fluids Struct.24, 1156–1165. 10.1016/j.jfluidstructs.2008.06.007
9
DillonR. H.FauciL. J.OmotoC.YangX. (2007). Fluid dynamic models of flagellar and ciliary beating. Ann. N. Y. Acad. Sci.1101, 494–505. 10.1196/annals.1389.016
10
FulfordG. R.BlakeJ. R. (1986). Muco-ciliary transport in the lung. J. Theor. Biol.121, 381–402.
11
GuiraoB.JoannyJ. F. (2007). Spontaneous creation of macroscopic flow and metachronal waves in an array of cilia. Biophys. J.92, 1900–1917. 10.1529/biophysj.106.084897
12
HaoJ.ZhuL. (2010). A lattice Boltzmann based implicit immersed boundary method for fluid–structure interaction. Comput. Math. Appl.59, 185–193. 10.1016/j.camwa.2009.06.055
13
HilferR. (2000). Applications of Fractional Calculus in Physics. River Edge, NJ: World Sci. Publishing.
14
Knight-JonesE. W. (1954). Relations between metachornism and the direction of ciliary beat in metazoa. Quart. J. Microsc. Sci.95, 503–521.
15
LafforgueO.BouguerraN.PoncetS.SeyssiecqI.FavierJ.ElkounS. (2017). Thermo-physical properties of synthetic mucus for the study of airway clearance. J. Biomed. Mater. Res. Part A105, 3025–3033. 10.1002/jbm.a.36161
16
LafforgueO.SeyssiecqI.PoncetS.FavierJ. (2018). Rheological properties of synthetic mucus for airway clearance. J. Biomed. Mater. Res. Part A106, 386–396. 10.1002/jbm.a.36251
17
LiebersteinH. M. (1975). Mathematical treatment of some problems in physiology. Comput. Math. Appl.1, 43–67.
18
MaitiS.PandeyS. K. (2017). Rheological fluid motion in tube by metachronal waves of cilia. Appl. Math. Mech.38, 393–410. 10.1007/s10483-017-2179-8
19
MaqboolK.ShaheenS.MannA. B. (2016). Exact solution of cilia induced flow of a Jeffrey fluid in an inclined tube. Springerplus5, 1–16. 10.1186/s40064-016-3021-8
20
MomaniS.OdibatZ. (2006). Analytical solution of a time-fractional Navier-Stokes equation by Adomian decomposition method. Appl. Math. Comput.177, 488–494. 10.1016/j.amc.2005.11.025
21
MurakamiA.TakahashiK. (1975). Correlation of electrical and mechanical responses in nervous control of cilia. Nature break 257, 48–49.
22
NortonM. M.RobinsonR. J.WeinsteinS. J. (2011). Model of ciliary clearance and the role of mucus rheology. Phys. Rev. E83:011921. 10.1103/PhysRevE.83.011921
23
OldhamK. B.SpanierJ. (1974). The Fractional Calculus. New York, NY: Academic Press.
24
RossS. M.CorrsinS. (1974). Results of an analytical model of mucociliary pumping. J. Appl. Physiol.37, 333–340.
25
RubinB. K. (2014). Secretion properties, clearance, and therapy in airway disease. Transl. Respir. Med.2, 1–7. 10.1186/2213-0802-2-6
26
ShapiroA. H.JaffrinM. Y.WeinbergS. L. (1969). Peristaltic pumping with long wavelength at low Reynolds number. J. Fluid Mech.37, 799–825.
27
SiddiquiA. M.FarooqA. A.RanaM. A. (2014). Study of MHD effects on the cilia-induced flow of a Newtonian fluid through a cylindrical tube. Magnetohydrodynamics50, 249–261.
28
SiddiquiA. M.HaroonT.RaniM.AnsariA. R. (2010). An analysis of the flow of a power law fluid due to ciliary motion in an infinite channel. J. Biorheol.24, 56–69. 10.1007/s12573-011-0026-3
29
SleighM. A.BlakeJ. R.LironN. (1988). The propulsion of mucus by cilia. Am. Rev. Respir. Dis.137, 726–741.
30
SmithD. J.GaffneyE. A.BlakeJ. R. (2008). Modelling mucociliary clearance. Respir. Physiol. Neurobiol.163, 178–188. 10.1016/j.resp.2008.03.006
31
TakahashiK.ShingyojiC. (2002). Control of eukaryotic flagella and cilia. Zool. Sci.19, 1393–1394. 10.2108/zsj.19.1393
32
Vélez-CorderoJ. R.LaugaE. (2013). Waving transport and propulsion in a generalized Newtonian fluid. J. Non Newton. Fluid Mech.199, 37–50. 10.1016/j.jnnfm.2013.05.006
Summary
Keywords
burgers fluid, cilia, fractional Adomian decomposition method, inclined tube, long wave length approximation
Citation
Mann AB, Shaheen S, Maqbool K and Poncet S (2019) Fractional Burgers Fluid Flow Due to Metachronal Ciliary Motion in an Inclined Tube. Front. Physiol. 10:588. doi: 10.3389/fphys.2019.00588
Received
21 September 2018
Accepted
26 April 2019
Published
16 May 2019
Volume
10 - 2019
Edited by
Francisco Monroy, Complutense University of Madrid, Spain
Reviewed by
Claudia Tanja Mierke, Leipzig University, Germany; Zakia Hammouch, Moulay Ismail University, Morocco
Updates
Copyright
© 2019 Mann, Shaheen, Maqbool and Poncet.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Amer Bilal Mann a_mann27@yahoo.co.uk
This article was submitted to Biophysics, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.