平行平板流动腔脉动流切应力的计算

高度远小于横向和纵向几何尺寸的矩形平行平板流动腔是人们用以体外研究细胞在切应力作用下力学行为的主要工具之一。大多数研究者主要对定常层流情进行研究。本文通过对矩形平行平板流动腔内的层流脉动流进行详细分析,给出腔内速和腔室底部切应力的准确计算公式。当Ｗｏｍｅｒｓｌｅｙ数较小时,准确公式简化为准定常公式。数值计算结果表明,在脉动流条件下,对于人们常用的流动腔几何尺寸,准定常公式具有相当高的精度。这为在脉     

具有切应力梯度的平行平板流动腔的构造

在近年来关于流体切应力与细胞力学行为之间关系的研究中,流体切应力梯度被认为是诱发动脉粥样硬化和内膜增生的重要因素之一.本文探讨如何利用常用的离体细胞力学实验工具--平行平板流动腔模拟具有梯度的定常流切应力环境.结果表明,根据Hele-Shaw流的原理和常用复势W(Z)=AZn(n＞1)的特性,可构造出具有各种切应力梯度分布的流动腔.与其它模拟切应力梯度的方法比较,本文的方法更加简洁、可行.     

局部狭窄血管中血液振荡流的速度和切应力

本文建立一种分析局部缓慢狭窄血管中血液振荡流的数学模型，给出了血液的轴向流速，径向流速和切应力的包含压力梯度项的解析表达式，并讨论了血管内由局部狭窄引起的压力梯度沿轴向变化的规律。文章以局部余弦狭窄为例进行数值计算，详细讨论上游均匀管段压力梯度的定常部分和不同次谐波对狭窄管段内流速和切应力的影响。数值结果表明，与均匀管情况相比，在狭窄段内，血液振荡流轴向流速无论平均值还是脉动幅值均明显增大，且径向流速不再为零。但径向流速仍远小于轴向流速。同时，切应力也不再仅由轴向流速梯度提供，径向流速梯度也将产生切应力，但是在计算管壁切向上的切应力时，径向流速梯度的贡献仍相当大。与均匀管管壁切应力沿流运方向保持恒定不同。狭窄管管壁切应力（平均值和脉动值）将随着狭窄高度的增大而增大，在狭窄最大高度处达到最大，因而沿流动方向产生了较大的切应力梯度。     

切应力协同下受热过冷层流液膜的破断特性

针对界面切应力协同下受热过冷层流液膜流动的破断过程, 建立了不同气液流向下的临界液膜厚度和最小润湿量的理论模型, 分析了不同驱动力作用下, 接触角、流体温度、界面切应力和壁面热流密度对液膜破断特性的影响. 研究表明: 临界液膜厚度和最小润湿量均随壁面热流密度的增加而增大; 重力驱动下的接触角影响在不同热流密度下有所不同, 流体温度在不同驱动力下对最小润湿量的影响截然相反; 同向切应力驱动下临界液膜厚度和最小润湿量随切应力增加而减小; 在重力和切应力协同驱动下, 同向切应力对最小润湿量的影响与重力和切应力所起作用的相对大小有关, 反向切应力使得临界液膜厚度和最小润湿量有所增大.      

一种力-电协同驱动的细胞微流控培养腔理论模型

细胞培养液在微流控生物反应器中受到外界物理场（如压力梯度或者电场）作用流动而产生流体剪应力,并进一步刺激种子细胞调控其内部基因的表达,从而促进细胞的分化和生长,这个过程在自然生命组织内的微管中亦是如此。考虑到细胞培养微腔隙中液体流动行为很难实验量化测定,理论建模分析是目前可行的研究手段。因此建立了矩形截面的细胞微流控培养腔理论模型,将外部的物理驱动场（压力梯度与电场）与培养腔内液体的流速、切应力和流率联系起来,分别得到了压力梯度驱动（Pressure gradient driven,PGD）、电场驱动（Electric field driven,EFD）及力-电协同驱动（Pressure-electricity synergic driven,P-ESD）三种驱动方式下的液体流动理论模型。结果表明该理论模型与现有的实验结果基本一致,具体地:力-电协同作用下的解答为压力梯度驱动和电场驱动结果的叠加。细胞培养腔内的流体流速、剪应力及流率幅值均正比于外部物理场强幅值,但随着压力梯度驱动载荷频率的增大而减小,随着电场驱动频率的变化不明显。在压力梯度驱动作用下,细胞贴壁处的切应力随着腔高的增大而线性增大,流率则随着腔高的增大而非线性增大,而电场驱动下的结果不受腔高的影响。生理范围内的温度场变化对压力和电场驱动的结果影响不大。另外,在引起细胞响应的流体切应力水平,电场驱动能提供较大的切应力幅值而压力梯度驱动则能提供较大的流率幅值。该理论模型的建立为细胞微流控生物反应器实验系统的设计及参数优化提供理论参考,同时也为力-电刺激细胞生长、分化机理的研究的提供基础。      

一种力–电协同驱动的细胞微流控培养腔理论模型

细胞培养液在微流控生物反应器中受到外界物理场(如压力梯度或者电场)作用流动而产生流体剪应力,并进一步刺激种子细胞调控其内部基因的表达,从而促进细胞的分化和生长,这个过程在自然生命组织内的微管中亦是如此.考虑到细胞培养微腔隙中液体流动行为很难实验量化测定,理论建模分析是目前可行的研究手段.因此建立了矩形截面的细胞微流控培养腔理论模型,将外部的物理驱动场(压力梯度与电场)与培养腔内液体的流速、切应力和流率联系起来,分别得到了压力梯度驱动(pressure gradient driven,PGD)、电场驱动(electric field driven,EFD)及力–电协同驱动(pressure-electricity synergic driven,P-ESD)三种驱动方式下的液体流动理论模型.结果表明该理论模型与现有的实验结果基本一致,即力–电协同作用下的解答为压力梯度驱动和电场驱动结果的叠加.细胞培养腔内的流体流速、剪应力及流率幅值均正比于外部物理场强幅值,但随着压力梯度驱动载荷频率的增大而减小,随着电场驱动频率的变化不明显.在压力梯度驱动作用下,细胞贴壁处的切应力随着腔高的增大而线性增大,流率则随着腔高的增大而非线性增大,而电场驱动下的结果不受腔高的影响.生理范围内的温度场变化对压力和电场驱动的结果影响不大.另外,在引起细胞响应的流体切应力水平,电场驱动能提供较大的切应力幅值而压力梯度驱动则能提供较大的流率幅值.该理论模型的建立为细胞微流控生物反应器实验系统的设计及参数优化提供理论参考,同时也为力–电刺激细胞生长、分化机理的研究的提供基础.     

一种骨小管中液体流动产生的流量及切应力模型

骨组织受力变形后其内部液体就会流动,同时在其微观结构——骨单元壁中扩散,并进一步产生一系列与骨液流动相关的物理效应,如流体剪切应力、流动电位等,这些物理效应被细胞感知并做出破骨或成骨等反应,来使骨适应外部载荷环境.鉴于骨组织产生的内部液体流动很难实验测定,理论模拟是目前的主要研究手段.基于骨单元的多孔弹性性质建立了骨小管内部液体的流动模型,该模型将骨单元所受的外部载荷与骨小管内部液体的压力、流速、流量和切应力联系起来,并进一步可以研究其力传导与力电传导机制.骨小管模型的建立分别基于中空和考虑哈弗液体的骨单元模型,并考虑了骨单元外壁的弹性约束和刚性位移约束两种边界条件.最终得到骨单元在外部轴向载荷作用下,骨小管内部液体的流量及流体切应力的解析解.结果表明:骨小管中的液体流量与流体切应力都正比于应变载荷幅值和频率,并由载荷的应变率决定.因此应变率可以作为控制流量和流体切应力的一种生理载荷因素.流量随着骨小管半径的增大而非线性增大,而流体切应力则随着骨小管半径的增大而线性增大.此外,在相同的载荷下,含哈弗液体的骨单元的模型中,骨小管中液体的流量和切应力均大于中空骨单元模型.     

刚性圆管中血液周期振荡流的切应力分布

本文通过求解圆管内血液振荡流的基本方程，求得圆管内血液流的压力梯度与切应力之间的关系式。在此基础上，详细讲座了圆管中轴向流速和切变率谐波的变化规律，指出流速谐波和切变率谐波的幅值都将随着谐波次数的增大而逐渐减小。为了使所得结果便于应用。文章通过管轴向中心线流速与压力梯度之间的关系式，进一步给出一种利用管轴向中心线流速计算管内切应力分布的简便方法。该方法用于检测活体血管内血液振荡流的切应力分布，具有操作简单，精度较高的优点。最后，以人体颈动脉为例，讨论血液周期振荡流的切应力的分布特性。发现在任意时刻，除了邻近管壁处切应力急剧增大到一定数值之外，沿管截面切应力分布相当均匀且接近于零，呈现出与定常流不同的切应力分布特征。     

内分液结构调控竖直管间歇流型建模研究

内分液流型调控管依靠微尺度网孔阻气通液的毛细力学特性,调控气液两相间歇流型以实现传热强化.基于Lockhart-Martinelli 分相模型以及Zuber-Findlay 漂移流动模型,建立描述内分液竖直管内流体动力特性的一维数学模型. 采用模型求解实验工况,计算结果与实验结果误差均在20% 以内. 计算发现,液速对流动现象起决定作用,而气速影响通过丝网的渗透程度. 在定性分析基础上,采用三角立方插值与最小二乘B 样条拟合获得了流动特性与气速、液速的定量函数关系. 据此得出结论,当Re_l < 693 7 时,一定出现第1 类工况;当Re_l > 693 7,且Re_g < 67 时,可能会出现第2 类工况,此时较低的气速会促进第2 类工况的出现. 根据建立的模型与拟合关系式可实现内分液调控管的优化设计.      

明流光滑平板边界层特性

本文应用正交偏振差动式激光测速系统,测量了明渠水流中顺流放置的光滑平板上边界层内的流速分布。特别是由于仪器的改进,比较详尽地量测了粘性底层中的流速分布,从而确定了壁面切应力。文中给出了切应力系数与雷诺数Re_x=U_0x/v的关系。 试验研究发现明流光滑平板边界层中,与无限边界绕流平板边界层一样,沿平板有层流边界层,转捩区和紊流边界层三部分。边界层厚度的计算公式也基本相同,明流平板的紊流边界层具有分区性质,沿高度将边界层分为粘性底层,过渡区和紊流区。在过渡区和紊流区中流速分布可用指数方程u/U_0=M(y/δ)~(1/n)表示。过渡区中n值在2—3之间。紊流区中M=1.0,n值随Re_x而变化,n=1.17Re_x~(0.135)紊流区也可用对数公式u~+=5.95logy~++5.48。这些公式中的常数与无限边界绕流平板边界层的情况略有差别,粘性底层中流速为线性分布,u~+=y~+。     

Effects of the side walls on unsteady flow of a second grade fluid over a plane wall

The effects of the side walls on unsteady flow of a second grade fluid over a plan wall are considered. The solution of the governing equation for velocity is obtained by the sine transform method. This gives a correct result for the shear stress at the bottom wall. The shear stress at the bottom wall is minimum at the middle of the plate and it increases near the side walls. It is shown that the mean thickness of the layer of the liquid over the plate increases with time and the ratio of the mean thickness to the distance between the side walls becomes ultimately 0.2714.     

Turbulent separated reattached flow in a two-dimensional curved-wall diffuser

A turbulent separation-reattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter. The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel. At the inlet of the diffuser, Reynolds number based on the diffuser height is 1.2×10⁵ and the velocity is 25.2m/s. The results of experiments are presented and analyzed in new defined streamline-aligned coordinates. The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress. A scale is formed using the maximum Reynolds shear stresses. It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law exists in the forward shear flow. Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model. The length scale is taken from that developed by Schofield and Perry. The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield-Perry velocity scale as well as the edge velocity of the boundary layer. The results of these experiments are presented in this paper.     

磁场和Hall电流对狭窄动脉中血液流动的影响

A micropolar model for blood simulating magnetohydrodynamic flow through a horizontally nonsymmetric but vertically symmetric artery with a mild stenosis is presented. To estimate the effect of the stenosis shape, a suitable geometry has been considered such that the horizontal shape of the stenosis can easily be changed just by varying a parameter referred to as the shape parameter. Flow parameters, such as velocity, the resistance to flow （the resistance impedance）, the wall shear stress distribution in the stenotic region, and its magnitude at the maximum height of the stenosis （stenosis throat）, have been computed for different shape parameters, the Hartmann number and the Hall parameter. This shows that the resistance to flow decreases with the increasing values of the parameter determining the stenosis shape and the Hail parameter, while it increases with the increasing Hartmann number. The wall shear stress and the shearing stress on the wall at the maximum height of the stenosis possess an inverse characteristic to the resistance to flow with respect to any given value of the Hartmann number and the Hall parameter. Finally, the effect of the Hartmann number and the Hall parameter on the horizontal velocity is examined.     

Hydrodynamics study in a plane ultrafiltration module using an electrochemical method and particle image velocimetry visualization

 The wall shear stress is determined at the surface of a plane sheet of Plexiglas, taking the place of a membrane, using an electrochemical method. Several microelectrodes are mounted flush to this plane plate, and maps of shear stress are determined for two inlet and outlet configurations and three channel heights. The heterogeneity of the wall shear stress is observed for both configurations. Furthermore, the study of the turbulence features of the flow shows a decreasing fluctuating rate of velocity gradient when the channel height is decreased. The wall velocity gradients and turbulent intensity rates analysis are confirmed by a flow visualization using the particle image velocimetry method. Received: 25 September 2000 / Accepted: 23 April 2001     

Two-phase non-linear model for blood flow in asymmetric and axisymmetric stenosed arteries

The pulsatile flow of a two-phase model for blood flow through axisymmetric and asymmetric stenosed narrow arteries is analyzed, treating blood as a two-phase model with the suspension of all the erythrocytes in the core region as the Herschel-Bulkley material and plasma in the peripheral layer as the Newtonian fluid. The perturbation method is applied to solve the resulting non-linear implicit system of partial differential equations. The expressions for various flow quantities are obtained. It is found that the pressure drop, plug core radius, wall shear stress increase as the yield stress or stenosis height increases. It is noted that the velocity increases, longitudinal impedance decreases as the amplitude increases. For asymmetric stenosis, the wall shear stress increases non-linearly with the increase of the axial distance. The estimates of the increase in longitudinal impedance to flow of the two-phase Herschel-Bulkley material are significantly lower than those of the single-phase Herschel-Bulkley material. The results show the advantages of two-phase flow over single-phase flow in small diameter arteries with stenosis.     

FLOW CHARACTERISTICS AROUND AN INCLINED ELLIPTIC CYLINDER IN A TURBULENT BOUNDARY LAYER

The flow characteristics around an inclined elliptic cylinder located near a flat plate were investigated experimentally. The axis ratio of the elliptic cylinder was AR=2. The pressure distributions along the surface of the cylinder and the flat plate were measured by varying the angle of attack of the elliptic cylinder. The velocity profiles behind the cylinder were measured using hot-wire anemometry. When the angle of attack varies, the peak pressure location on the windward cylinder surface moves towards the rear edge of the cylinder, while that on the leeward surface moves towards the front edge of the cylinder. The vortex-shedding frequency also gradually decreases, defining a critical angle of attack for each gap ratio. The location of the minimum pressure on the flat plate surface moves downstream for positive angles of attack, while it moves upstream for negative angles of attack. Negative angles of attack cause a greater disturbance in the boundary layer near the wall compared to positive angles of attack. This shows that the separated wall shear layer from the boundary layer and the lower shear layer of the cylinder wake are strongly merged compared to other cases.     

Wall pressure and shear stress measurements beneath an impinging jet

This paper presents comprehensive measurements of wall pressure and surface shear stress beneath a plane, two-dimensional, turbulent jet impinging normally onto a flat surface. The results cover a wider range of Reynolds number and ratio of impingement height (H) to nozzle gap (D) than do previous studies. The pressure distributions are nearly Gaussian, independent of Reynolds number, and closely balance the momentum flux from the jet nozzle as H/D varies. Particular attention was paid to probe size in measuring the wall shear stress because this has a significant effect on the results. A range of Preston tubes and Stanton probes were tested from which it was found that a 0.05-mm-high Stanton probe—the smallest that we could make—appeared to give accurate results. As expected, the shape of the wall shear stress distributions depended both on H/D and on Reynolds number. Furthermore, the relation between wall pressure and shear stress from Hiemenz's theoretical solution for stagnation flow is not in agreement with the results. It is postulated that the discrepancy is due to the relatively high free-stream turbulence level in the jet. Future papers will document the mean flow field and turbulence and the time dependence of the surface pressure.     

水平旋转空腔环流的壁面应力

通过对典型的水平旋转内消能泄洪洞空腔环流的试验观测，研究了其壁面应力的变化规律。空腔环流的壁面压强在水平洞的起始段壁面压强急剧减小然后回升，具有过渡段的性质，沿程波状减小，符合对数变化规律，但不同的流态，对数律的参数的变化是不同的：内界面相对压强Po/Pwz。在淹沿流流态时，随(H-h)/h的变化率显著不同，在吸允流流态时却基本相同；壁面切应力沿z的变化规律为先急剧减小，随后缓慢减小至零，主要与环流特性有密切的关系。