Determination of separated region in a curved tube

The formation of atherosclerosis in a curved aorta is closely related to the existence of separated vortex region. This paper deals with the steady laminar motion of an incompressible Newtonian fluid through a curved tube with circular cross-section whose curvature is small and whose curvature gradient is not too large. Using the momentum integral method and the approximation of quasi-constant curvature, an equation which determines the location of separation and reattachment is derived. From this equation the earliest point of separation and the corresponding critical Reynolds number are obtained, and the relation between the position of separation and reattachment and Reynolds number R_e for different azimuthal angle are revealed. It is concluded that the separation first emerges at the position whose curvature gradient has the maximum absolute value. With increasing R_e, the separation region extends in the direction of mainstream, azimuthal angle and radius vector, and then forms a three-dimensional separated vortex, which gradually enlarges in all three directions with the increase of Reynolds number. The theoretical results also very clearly demonstrate the following striking experimental fact: if a symmetrical curved tube exhibits a separated vortex at the outside of the upstream, then it must have another one symmetrically placed at the inside of its downstream.     

The displacement wave theory of blood vessel

On the basis of the medical and mechanical analysis and explanations in this paper the visco-elastic simply supported beam model is proposed to treat the displacement wave of the blood vessels. The relationships between the displacement wave and blood vessel elasticity as well as the viscous dissipation of the blood and blood vessel are obtained. The corresponding relations of such kinds of pulses in the traditional Chinese medicine as smooth pulse, surface pulse and deep pulse to the displacement waves of blood vessels are also found. The computational results are in good agreement with those acquired in the experiments with ultrasonic wave.The Projects Supported by National Science Fundation of China.     

无限长圆管中任意轴对称体绕中心轴旋转的小Re数问题

本文以Sampson旋转子在圆管中反演得到的无穷级数作为基本奇点,对分布函数采用离散奇点线分布法和分段二次抛物分布法,解决了任意轴对称体在无限长圆管中绕中心轴旋转的Stokes流动问題。通过对长球的数值计算,比较了两种方法,结果表明分段二次抛物分布法各方面都优于离散法。此外,本文还研究球串在无穷长圆管中绕中心轴旋转的小Re数问题,第一次给出了每个球所受的阻力矩系数。     

The effect of the vasomotion on the mass transfer in microcirculation

Detailed structure of the attracting set of the piecewise linear Henon mapping(x, y)→(1- a|x|+by,x)with a=8/5 and b=9/25 is described in this paper using the method of dual line mapping. Let A and B denote the fixed saddles in the first quadrant, and in the third quadrant, respectively. It is claimed that(1)the attracting set is the closure of the unstable manifold of saddle B, which includes the unstable manifold of A as its subset, and(2)the basin of attraction is the closure of the stable manifold of A, bounded by the stable manifold of B, which is in the limiting set of the stable manifold of A.Relations of the manifolds of the periodic saddles with the manifolds of the fixed point are given. Symbolic dynamics notations are adopted which renders possible the study of the dynamical behavior of every piece of the manifolds and of every homoclinic or heteroclinic point.     

The inlet effect on the drag factor of a sphere in a tube

The creeping motion Ground a sphere situated axisymmetrically near the entrance of asemi-infinite circular cylindrical tube is analyzed using infinite series solutions for thevelocity components. pressure and the stream function. Truncating the infinite series. thecorresponding coefficients in the series are determined by a collocation technique. The dragfactor and the stress distribution on the surface of the sphere are calculated for the sphere inmotion in quiescent fluid and for the flow with uniform velocity at the entrance past a rigidlyheld sphere. The results indicate that a sphere near the entrance which has a uniformentrance velocity profile will suffer larger drag than that in infinite tube.Theconvergence of the collocation technique is tested by numerical calculation. It is shown thatthe technique has good convergence properties.     

THE METHOD OF CONTINUOUS DISTRIBUTION OF SINGULARITIES TO TREAT THE STOKES FLOW OF THE ARBITRARY OBLATE AXISYMWETRICAL BODY

This paper deals with the Stokes flow of the arbitrary oblate axisymmetrical body by means of constant density and quadratic distribution function approximation for the method of continuous distribution of singularities. The Sampson spherical infinite series arc chosen as fundamental singularities. The convergence, accuracy and range of application of both two approximations are examined by the unbounded Stokes flow past the oblate spheroid. It is demonstrated that the drag factor and pressure distribution both conform with the exact solution very well. Besides, the properties, accuracy and the range of application are getting belter with the improving of the approximation of the distribution function. As an example of the arbitrary oblate axisymmetrical bodies, the Stokes flow of the oblate Cassini oval are calculated by these two methods and the results are convergent and consistent. Finally, with the quadratic distribution approximation the red blood cell, which has physiologic meaning, is considered and for the first time the (orresponding drag factor and pressure distribution on the surface of the cell are obtained.     

The motion of viscous liquid column with finite length in a vertical straight capillary tube

This paper deals with the motion of viscous liquid column with finite length and two free surfaces in a vertical straight capillary tube. It is assumed that fluid is Newtonian. Linearizing the boundary conditions, analytic expressions in the form of infinite series have been obtained for velocity, piessure and free surface at low Reynolds number. The numerical calculation is carried out for a set of cylinder’s length of water and blood. It has been revealed that there are considerable circulating currents at the upper and lower meniscuses. Its maximum velocity is about 57% of the average velocity of the mainstream. Iner-tial effect is also studied in this paper. Using the time-dependent method in finite difference techniques, numerical solution of the corresponding nonlinear equation at Re<24.5 is computed. Comparing it with analytic exact solution at low Reynolds number shows that inertial effect is negligible provided Re<24.5.     

多项式基函数法

提出一种新型的数值计算方法--基函数法.此方法直接在非结构网格上离散微分算子,采用基函数展开逼近真实函数,构造出了导数的中心格式和迎风格式,取二阶多项式为基函数,并采用通量分裂法及中心格式和迎风格式相结合的技术以消除激波附近的非物理波动,构造出数值求解无粘可压缩流动二阶多项式的基函数格式,通过多个二维无粘超音速和跨音速可压缩流动典型算例的数值计算表明,该方法是一种高精度的、对激波具有高分辨率的无波动新型数值计算方法,与网格自适应技术相结合可得到十分满意的结果.     

Basic function scheme of polynomial type

A new numerical method named as basic function method is proposed. It can directly discretize differential operators on unstructured grids. By expanding the basic function to approach the exact function, the central and upwind schemes of derivative are constructed. By using the second-order polynomial as a basic function and applying the flux splitting method and the combination of central and upwind schemes to suppress non-physical fluctuation near shock waves, a second-order basic function scheme of polynomial type is proposed to solve inviscid compressible flows numerically. Numerical results of typical examples for two-dimensional inviscid compressible transonic and supersonic steady flows indicate that the new scheme has high accuracy and high resolution for shock waves. Combined with the adaptive remeshing technique, satisfactory results can be obtained.     

弯曲圆管内漩涡结构分叉现象的理论研究

利用拓扑结构分析方法,分析了弯曲圆管内定常流在横截面上流线的奇点个数及分布规律,给出了二次流的漩涡数目由2个变为4个,流态结构发生分叉现象的理论判据。进而,利用Galerkin方法,得到了弯曲圆管内定常流的流函数和轴向速度的半解析表达式,给出了流态结构发生分叉现象的临界Dean数,所得结果与理论判则一致。