Javid Khurram, Riaz Muhammad, Chu Yu-Ming, Ijaz Khan M, Ullah Khan Sami, Kadry S
Department of Mathematics, Northern University, Nowshera, 24100, KPK, Pakistan.
Department of Mathematics, Huzhou University, Huzhou 313000, PR China; Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering, Changsha University of Science & Technology, Changsha 410114, PR China.
Comput Methods Programs Biomed. 2021 Mar;200:105926. doi: 10.1016/j.cmpb.2020.105926. Epub 2021 Jan 4.
Now-a-days in medical science, the transport study of biological fluids through non-uniform vessels are going to increase due to their close relation to the reality. Motivated through such type of complex transportation, the current study is presented of cilia hydro-dynamics of an aqueous electrolytic viscous fluid through a non-uniform channel under an applied axial electric field. Mathematical Formulations: Because of the complexity shape and nature of flow channel, we have used curvilinear coordinates in the derivation of continuity and momentum equationsin a fixed frame of reference. A linear transformation is used to renovate the flow system of equations from fixed (laboratory) to moving (wave) frame. For further simplification, the dimensionless variables are introduced to make the flow system of equations into the dimensionless form and at last convert these equations in term of stream function by using the mathematical terminologies of streamlines. The whole analysis is performed under (low Reynolds number) creeping phenomena and long wavelength approximation, respectively. Additionally, small ionic Peclet number and Debye-Huckel linearization are used to simplify the Nernst-Planck and Poisson-Boltzmann equations. The BVP4C technique is used to obtain the numerical solution for velocity distribution, pressure gradient, pressure rise and stream function through MATLAB.
The amplitude of velocity distribution is increased (decreased) at larger values of non-uniform parameter (cilia length). The non-uniform parameter played a vital role not only in the enhancement of circulation at the upper half of the channel but also the length of bolus increased. Results of straight channel are gained for larger value of the dimensionless radius of curvature parameter as well as cilia length.
如今在医学领域,由于生物流体通过非均匀血管的输运研究与实际情况密切相关,此类研究正日益增多。受这种复杂输运的启发,本文开展了在轴向电场作用下,水性电解粘性流体通过非均匀通道的纤毛流体动力学研究。
由于流动通道形状和性质复杂,我们在固定参考系中推导连续性方程和动量方程时采用了曲线坐标。通过线性变换将流动方程组从固定(实验室)参考系转换到运动(波动)参考系。为进一步简化,引入无量纲变量使流动方程组化为无量纲形式,最后利用流线的数学术语将这些方程转化为流函数形式。整个分析分别在(低雷诺数)蠕动现象和长波长近似条件下进行。此外,采用小离子佩克莱数和德拜 - 休克尔线性化来简化能斯特 - 普朗克方程和泊松 - 玻尔兹曼方程。利用BVP4C技术通过MATLAB获得速度分布、压力梯度、压力升和流函数的数值解。
在非均匀参数(纤毛长度)较大时,速度分布的幅值增大(减小)。非均匀参数不仅对通道上半部分的循环增强起着至关重要的作用,而且团块长度也增加了。对于较大的无量纲曲率半径参数值以及纤毛长度,可得到直通道的结果。