弹性动脉血管中血液流动特性的模拟和分析

研究了两个不同的非牛顿血液流动模型:低粘性剪切简单幂律模型和低粘性剪切及粘弹性振荡流的广义Maxwell模型．同时利用这两个非牛顿模型和牛顿模型,研究了磁场中刚性和弹性直血管中血液的正弦型脉动．在生理学条件下,大动脉中血液的弹性对其流动性态似乎并不产生影响,单纯低粘性剪切模型可以逼真地模拟这种血液流动．利用高剪切幂律模型模拟弹性血管中的正弦型脉动流,发现在同一压力梯度下,与牛顿流体相比较,幂律流体的平均流率和流率变化幅度都更小．控制方程用Crank-Niclson方法求解．弹性动脉中血液受磁场作用是产生此结果的直观原因．在主动脉生物流的模拟中,与牛顿流体模型比较,发现在匹配流率曲线上,幂律模型的平均壁面剪切应力增大,峰值壁面剪切应力减小．讨论了弹性血管横切磁场时的血液流动,评估了血管形状和表面不规则等因素的影响．     

A sphere in a second degree polynomial creeping flow parallel to a wall

Comprehensive results are provided for the creeping flow arounda spherical particle in a viscous fluid close to a plane wall,when the external velocity is parallel to the wall and variesas a second degree polynomial in the coordinates. By linearityof Stokes equations, the solution is a sum of flows for typicalunperturbed flows: a pure shear flow, a ‘modulated shearflow’, for which the rate of shear varies linearly inthe direction normal to the wall, and a quadratic flow. Solutionsconsidered here use the bipolar coordinates technique. Theycomplement the accurate results of Chaoui and Feuillebois (2003)for the pure shear flow. The solution of Goren and O'Neill (1971)for the quadratic flow is reconsidered and a new analyticalsolution is derived for the ambient modulated shear flow. Theperturbed flow fields for these two cases are presented in detailand discussed. Results for the force and torque friction factorsare provided with a 5 x 10^–17 accuracy as a reference.For the quadratic flow, there is a force and a torque on a fixedsphere. A minimum value of the torque is found for a gap ofabout 0·18a, where a is the sphere radius. This minimumis interpreted in term of the corresponding flow structure.For the modulated shear flow, there is only a torque. The freemotion of a sphere in an ambient quadratic flow is also determined.     

Numerical computation of pulsatile flow through a locally constricted channel

This paper deals with the numerical solution of a pulsatile laminar flow through a locally constricted channel. A finite difference technique has been employed to solve the governing equations. The effects of the flow parameters such as Reynolds number, flow pulsation in terms of Strouhal number, constriction height and length on the flow behaviour have been studied. It is found that the peak value of the wall shear stress has significantly changed with the variation of Reynolds numbers and constriction heights. It is also noted that the Strouhal number and constriction length have little effect on the peak value of the wall shear stress. The flow computation reveals that the peak value of the wall shear stress at maximum flow rate time in pulsatile flow situation is much larger than that due to steady flow. The constriction and the flow pulsation produce flow disturbances at the vicinity of the constriction of the channel in the downstream direction.     

Dynamic deformation of a stationary rubber flow

The effect of dynamic deformation on the stationary flow of a rubber composition has been experimentally investigated for comparable values of the stationary and dynamic strain rates. The dependence of the effective viscosity on the stationary shear rate is not equivalent to its dependence on the periodic shear rate amplitude. An expression is given for calculating the effective viscosity in the case of combined stationary and dynamic shear deformation. The effectiveness of the dynamic deformation, estimated in terms of the effective viscosity, depends on whether it is superimposed on the stationary flow at constant stationary shear rate or at constant stress. It is proposed to estimate the effectiveness of dynamic deformation of a stationary non-Newtonian flow in terms of the change in the power of the stationary forces. When the effective viscosity is reduced by dynamic deformation of the stationary flow, the power of the stationary forces increases at constant shear stress and falls at constant stationary shear rate.Institute of Polymer Mechanics, Academy of Sciences of the Latvian SSR, Riga. Translated from Mekhanika Polimerov, No. 3, pp. 489–496, May–June, 1971.     

Experimental study of simple shear flow of monodisperse fibrous composites

Steady-state simple shear flow of epoxide resin and its composites containing 1, 3, or 5 vol. % fiber material was studied at temperatures ranging from 20 to 70°C. Normal stresses were measured. The decrease in viscosity with increasing shear rate, during the formation of a continuous fiber mesh, is interpreted according to the molecular-kinetic theory of flow of continuous media, on the assumption of viscoelastic flow elements with the relaxation time a function of the shear rate. The conditions are established for a transition from the slip mechanism of flow, which involves fiber linkages, to the cluster mechanism of flow.     

Analysis of blood flow through a modelled artery with an aneurysm

The intention of the present work is to carry out a systematic analysis of flow features in a tube, modelled as artery, having a local aneurysm in presence of haematocrit. The arterial model is treated to be axi-symmetric and rigid. The blood, flowing through the modelled artery, is treated to be Newtonian and non-homogeneous. For a thorough quantitative analysis of the flow characteristics such as wall pressure, flow velocity, wall shear stress, the unsteady incompressible Navier-Stokes equations in cylindrical polar co-ordinates under the laminar flow conditions are solved by using the finite-difference method. Finally, the numerical illustrations presented in this paper provide an effective measure to estimate the combined influence of haematocrit and aneurysm on flow characteristics. It is found that the magnitude of wall shear stress and also the length of separation increase with increasing values of the haematocrit parameter. The length of flow separation increases but the peak value of wall shear stress decreases with the increasing length of aneurysm. The peak value of wall shear stress as well as the length of separation increases with the increasing height of the aneurysm.     

An improved dynamic subgrid-scale model and its application to large eddy simulation of stratified channel flows

In the present paper, a new dynamic subgrid-scale (SGS) model of turbulent stress and heat flux for stratified shear flow is proposed. Based on our calculated results of stratified channel flow, the dynamic subgrid-scale model developed in this paper is shown to be effective for large eddy simulation (LES) of stratified turbulent shear flows. The new SGS model is then applied to the LES of the stratified turbulent channel flow to investigate the coupled shear and buoyancy effects on the behavior of turbulent statistics, turbulent heat transfer and flow structures at different Richardson numbers.     

Instability of power-law fluid flows down an incline subjected to wind stress

In this paper we investigate the effect of a prescribed superficial shear stress on the generation and structure of roll waves developing from infinitesimal disturbances on the surface of a power-law fluid layer flowing down an incline. The unsteady equations of motion are depth integrated according to the von Kármán momentum integral method to obtain a non-homogeneous system of nonlinear hyperbolic conservation laws governing the average flow rate and the thickness of the fluid layer. By conducting a linear stability analysis we obtain an analytical formula for the critical conditions for the onset of instability of a uniform and steady flow in terms of the prescribed surface shear stress. A nonlinear analysis is performed by numerically calculating the nonlinear evolution of a perturbed flow. The calculation is carried out using a high-resolution finite volume scheme. The source term is handled by implementing the quasi-steady wave propagation algorithm. Conclusions are drawn regarding the effect of the applied surface shear stress parameter and flow conditions on the development and characteristics of the roll waves arising from the instability. For a Newtonian flow subjected to a prescribed superficial shear stress, using an analytical theory, we show that the nonlinear governing equations do not admit roll waves solutions under conditions when the uniform and steady flow is linearly stable. For the case of a general power-law fluid flow with zero shear stress applied at the surface, the analytical investigation leads to a procedure for calculating the characteristics of a roll waves flow. These results are compared with those yielded by the numerical procedure.     

Long Eddies in Sheared Flows

The characteristic feature of the wide variety of hydraulic shear flows analyzed in this study is that they all contain a critical level where some of the fluid is turned relative to the ambient flow. One example is the flow produced in a thin layer of fluid, contained between lateral boundaries, during the passage of a long eddy. The boundaries of the layer may be rigid, or flexible, or free; the fluid may be either compressible or incompressible. A further example is the flow produced when a shear layer separates from a rigid boundary producing a region of recirculating flow. The equations used in this study are those governing inviscid hydraulic shear flows. They are similar in form to the classical boundary layer equations with the viscous term omitted. The main result of the study is to show that when the hydraulic flow is steady and contained between lateral boundaries, the variation of vorticity ω(ψ) cannot be prescribed at any streamline which crosses the critical level. This variation is, in fact, determined by (1) the vorticity distribution at all streamlines which do not cross the critical level, by (2) the auxiliary conditions which must be satisfied at the boundaries of the fluid layer, and by (3) the dimensions of the region containing the turned flow. If at some instant the vorticity distribution is specified arbitrarily at all streamlines, generally the subsequent flow will be unsteady. In order to emphasize this point, a class of exact solutions describing unsteady hydraulic flows are derived. These are used to describe the flow produced by the passage of a long eddy which distorts as it is convected with the ambient flow. They are also used to describe the unsteady flow that is produced when a shear layer separates from a boundary. Examples are given both of flows in which the shear layer reattaches after separation and of flows in which the shear layer does not reattach. When the shear layer vorticity distribution has the form ωαyⁿ, where y is a distance measure across the layer, the steady flows are of Falkner-Skan type inside, and adjacent to, the separation region. The unsteady flows described in this paper are natural generalizations of these Falkner-Skan flows. One important result of the analysis is to show that if the unsteady flow inside the separation region is strongly sheared, then the boundary of the separation region moves upstream towards the point of separation, forming large transverse currents. Generally, the assumption of hydraulic flow becomes invalid in a finite time. On the other hand, if the flow inside the separation region is weakly sheared, this region is swept downstream and the flow becomes self-similar.     

On the homotopy solution for Poiseuille flow of a fourth grade fluid

This paper develops a mathematical model with an aim to compute the analytic solution for the flow of a fourth grade fluid between two fixed porous walls. The flow is induced under the application of a constant pressure gradient. The arising nonlinear problem is treated analytically yielding a series solution by homotopy analysis method (HAM). Results of velocity and shear stresses at the walls are obtained. The impacts of several flow parameters are examined on the velocity and shear stresses.     

Hele Shaw Convection with Imposed Shear Flows: Boundary-Layer Formulation

The nonlinear convection forced by the boundaries of a Hele Shaw cell to align perpendicular to an imposed shear flow was analytically investigated by the boundary-layer method. The imposed shear flow may be a Couette flow that extends throughout the convecting layer or flow confined to a boundary, depending on the geometry of the Hele Shaw cell. This study examined the case in which the imposed shear flow has a boundary-layer structure and its interaction with the convecting interior. Analytical solutions for both the boundary layer and interior were obtained. The study revealed the following.For large aspect ratio A , the interaction of the imposed shear flow and convection is confined to the boundary layer. The boundary layer is a viscous rather than a thermal layer. The results showed that the range of validity of the Hele Shaw equations used in the literature is of order 1/ A ². For an asymptotically large aspect ratio A up to order 1/ A ², the velocity in the y -direction must be zero. The velocity in the x -direction and the z -direction has a parabolic dependence on y , but the temperature perturbation does not depend on y . These results may have implication for convection in porous media.     

Mathematical modeling of pulsatile flow of non-Newtonian fluid in stenosed arteries

The pulsatile flow of blood through mild stenosed artery is studied. The effects of pulsatility, stenosis and non-Newtonian behavior of blood, treating the blood as Herschel–Bulkley fluid, are simultaneously considered. A perturbation method is used to analyze the flow. The expressions for the shear stress, velocity, flow rate, wall shear stress, longitudinal impedance and the plug core radius have been obtained. The variations of these flow quantities with different parameters of the fluid have been analyzed. It is found that, the plug core radius, pressure drop and wall shear stress increase with the increase of yield stress or the stenosis height. The velocity and the wall shear stress increase considerably with the increase in the amplitude of the pressure drop. It is clear that for a given value of stenosis height and for the increasing values of the stenosis shape parameter from 3 to 6, there is a sharp increase in the impedance of the flow and also the plots are skewed to the right-hand side. It is observed that the estimates of the increase in the longitudinal impedance increase with the increase of the axial distance or with the increase of the stenosis height. The present study also brings out the effects of asymmetric of the stenosis on the flow quantities.     

Nonlinear mathematical analysis for blood flow in a constricted artery under periodic body acceleration

This study analyses the pulsatile flow of blood through mild stenosed narrow arteries, treating the blood in the core region as a Casson fluid and the plasma in the peripheral layer as a Newtonian fluid. Perturbation method is employed to solve the resulting coupled implicit system of non-linear partial differential equations. The expressions for shear stress, velocity, wall shear stress, plug core radius, flow rate and longitudinal impedance to flow are obtained. The effects of pulsatility, stenosis depth, peripheral layer thickness, body acceleration and non-Newtonian behavior of blood on these flow quantities are discussed. It is noted that the plug core radius, wall shear stress and longitudinal impedance to flow increase as the yield stress and stenosis depth increase and they decrease with the increase of the body acceleration, pressure gradient, width of the peripheral layer thickness. It is observed that the plug flow velocity and flow rate increase with the increase of the pulsatile Reynolds number, body acceleration, pressure gradient and the width of the peripheral layer thickness and the reverse behavior is found when the yield stress, stenosis depth and lead angle increase. It is also recorded that the wall shear stress and longitudinal impedance to flow are considerably lower for the two-fluid Casson model than that of the single-fluid Casson model. It is found that the presence of body acceleration and peripheral layer influences the mean flow rate and mean velocity by increasing their magnitude significantly in the arteries.     

非Newton流体的物质点法模拟研究

准确模拟非Newton流体的运动特性具有重要的工程意义.物质点法作为一种相对新兴的粒子型算法,其结合了Lagrange算法和Euler算法的双重优势,已广泛有效地应用于各个工程领域.基于物质点法,结合人工状态方程,分析了两种非Newton流体(cross流体和幂律流体)在平板Poiseuille流和Couette流情况下的流动特性.结果表明:对Newton流体,物质点模拟结果与理论值一致;对非Newton流体,物质点法可准确模拟其剪切稀化和剪切稠化现象.表明了物质点法在模拟非Newton流体流动问题时的适用性,拓展了物质点法的应用范围.     

Linear Viscoelasticity of Liquid-Crystalline Polymers

A linear (small-amplitude) periodic shear deformation of anisotropic viscoelastic liquids obeying the Akay–Leslie rheological model is considered. The frequency dependences of the real and imaginary components of the complex shear modulus and complex normal-stress coefficient are determined. A comparison between calculation results and test data on the shear flow of poly(-benzylglutamate) in m-cresol is carried out. It is stated that, if the material is characterized by some initial orientation, both components of the complex shear modulus contain a multiplier which depends on the degree of the initial orientation and increases the values of the components compared with those for an initially isotropic material. The model predicts that, in a periodic shear flow, the components of shear and normal stresses are constant and, like the components of shear modulus, are independent of deformation frequency. If the parameter d ₀ of the Akay–Leslie model is equal to zero, the values of its other parameters can be determined from experimental results on periodic shear flow.     

Large shear rate behavior for the Hébraud-Lequeux model

The Hébraud-Lequeux model is a model describing the flow of soft glassy material in a simple shear flow configuration.It is given by a kinetic/Fokker-Planck-type equation whose coefficients depend on the shear rate of the experiment.In this paper we want to study what happens to the stationary solutions of this model when the shear rate is asymptotically large.In order to do that,we expand the solution of the equation using singular perturbation tools.In the end,we rigorously prove the estimate of Hébraud and Lequeux that the material asymptotically behaves as a Newtonian fluid.     

Determination of arterial wall shear stress

The arteries can remodel their structure and function to adapt themselves to the mechanical environment. In various factors that lead to vascular remodeling, the shear stress on the arterial wall induced by the blood flow is of great importance. However, there are many technique difficulties in measuring the wall shear stress directly at present. In this paper, through analyzing the pulsatile blood flow in arteries, a method has been proposed that can determine the wall shear stress quantitatively by measuring the velocity on the arterial axis, and that provides a necessary means to discuss the influence of arterial wall shear stress on vascular remodeling.     

Instability in Parallel Flows Revisited

An example of an unstable inviscid plane parallel shear flow with classical boundary conditions is presented. The complete unstable spectrum is exhibited using techniques of continued fractions for the shear flow with profile U ( y )=cos  m y . For such flows spectral instability implies nonlinear instability. A three-dimensional generalization is discussed.     

Numerical simulation of shear and normal forces in standard and shear flow sequencing batch reactor (SF-SBR)

The purpose of this contribution is to compare the shear and the normal stresses in two different types of bioreactors. In the first one (SBR), the granules are generated by liquid and bubbles flow, whereas in the second one (SF-SBR), the shear rate is achieved by installing a rotating cylinder inside the reactor. Such shear flow can be applied for anaerobic process in the wastewater treatment. The results demonstrate the effective role of the process parameters and the reactor geometry on the shear and normal stresses and consequently on the granulation process. Hereby, the different tendencies of the velocity fields, the particle sedimentation as well as Taylor vortices in (SF-SBR), are observed. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)