分叉血管中动脉狭窄对血液流动的影响

利用血液流动的粘性不可压Navier-Stokes方程，采用有限元方法研究了分叉血管中动脉局部狭窄对血液流动的影响，分别对母血管壁有动脉狭窄病变和无病变时的血液流动进行了数值模拟。针对不同程度的动脉狭窄，得到了当狭窄度S分别为0.05、0.25、0.5时血液流动的流线图和压力图。一方面，从计算的压力图中可以看出当血液流经狭窄区域时，压力会迅速发生变化；另一方面，观察流线图可知，在血液流经狭窄的区域时会出现流动分离现象，产生涡的结构。根据上述结果，可以看出动脉局部狭窄会对血液流动产生很大的影响。动脉粥样硬化等某些血管疾病的发病机制和病变发展与血液流动的力学特征有密切关系。在动脉局部狭窄的初期，动脉狭窄邻近区域的血液几乎光滑地流过狭窄区，流动产生的分离小；随着狭窄区域半径的逐渐增大，血流流动产生比较明显的分离，在分离区内存在层流和回流；而当狭窄区域的半径进一步增大时，血管内血流流动分离区域扩大，而且在分离区内的血液流动会出现局部湍流。这种情况可能会使血液中的血小板、纤维蛋白、脂肪颗粒沉积，最终导致血流不畅或动脉粥样硬化等血管疾病。     

弯曲动脉的血流动力学数值分析

利用计算流体力学的理论和方法对弯曲动脉中的血流动力学进行数值分析，是研究心血管疾病流体动力学机理的一种行之有效的方法。本文将升主动脉、主动脉弓和降主动脉联系起来作为弯曲动脉几何模型，给出了血液流动的边界条件以及计算条件。根据生理脉动流条件，对狗的弯曲动脉几何模型内发展中的血液流动进行了有限元数值模拟，并利用可视化方法对血液流动的轴向速度、二次流、壁面切应力等计算结果进行了分析。研究结果表明，在弯管内侧壁处，同时存在主流方向和二次流方向的回流，此处容易形成涡流。弯管内侧壁比外侧壁的壁面切应力具有更强的脉动性。     

弯曲狭窄管内血液流的分析

与血管狭窄有关的异常血液动力学特征在血管疾病的发生和发展过程中起着重要的作用,由于血管狭窄和弯曲的综合影响,将会出现一系列有趣的流体力学现象,本文研究具有对称狭窄的弯曲小动脉内定常血液流动,在一定的假设条件下,直接从支配血液流动的Navier-Stokes方程求出问题的摄动解,由此求得弯曲狭窄管內血液流动的轴向速度、二次流速度及压力梯度等分析表达式,并进一步求得轴向和周向血管壁切应力。本文的结果是先前有关狭窄直管和弯曲均匀管流动研究的拓广。     

棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

Willis环瞬态血液流体动力特性数值分析

针对左颈内动脉在无狭窄和完全阻塞情况下交通动脉中的血液瞬态流动采用颅内Willis环二维模型进行计算分析.计算得到完整环、交通动脉大脑前动脉缺失情况下各交通动脉血液流量变化趋势.结果表明交通动脉中血液流量变化最大的情况为:右侧大脑前动脉缺失同时左颈内动脉阻塞;在左(右)大脑前动脉近端缺失时前交通动脉中血液流量变化趋势与左(右)大脑前动脉远端相同;左侧后交通动脉在左颈内动脉狭窄比较严重时血液流动方向就会发生改变;在完整环、左颈内动脉阻塞时前交通动脉中的血液流量大于后交通动脉.     

弯曲动脉壁非线性弹性力学性质的理论分析

血管壁非线性力学特性是探索心血管中的血液流动规律及脉搏波传播现象的一个重要前提．在血管力学研究中,直管动脉壁的本构方程已有相当深人的研究,唯独象主动脉弓那样的弯管动脉壁本构方程;至今还没有建立一个理论模型．文中在已有的直管动脉壁本构方程研究基础上,提出一个理论方法来分析弯曲动脉壁的非线性弹性力学性质,在弯曲动脉壁被模拟为均质、正交各向异性、不可压缩材料的前提下,作者从理论上建立了一个表达弯管动脉壁非线性弹性性质的三维ｅ指数型本构方程文中还探讨了弯曲动脉壁内的残余应变分析．     

动脉中脉动流的准二维分析与数值计算

本文考虑了血液的菲牛顿流动特性,对血液在动脉管系中的脉动流建立了准二维流动模型。利用有限差分方法得到了动脉管系内非牛顿流体的准二维不定常流动的数值解。并以人体五根主要动脉所组成的动脉管系为例,进行了详细的数值计算。计算结果表明,在动脉的某些典型位置上。计算所得的理论波形与实测波形是相似的。     

工程湍流模式理论综述及展望

本文讨论了国内外湍流模拟的现状和发展趋势．指出湍流模式的建立除了应遵循理性力学原则外,还必须密切结合工程流体的复杂流动现象,如对具有浮力的回流、分离流及强旋流的模拟,对逆梯度的输运模拟,对单相流、多相流,单流体、多流体的湍流牛顿流体及湍流的非牛顿流体的模拟．应加深对湍流机理的认识,改进湍流模拟手段,结合工程实际,提出较为通用的工程湍流模式．      

脑Willis环的一维血流动力学及氧输运特性的数值研究

Willis环是大脑侧枝循环的重要组成部分, 研究其血流动力学特性以及氧输运规律对脑缺血疾病的认知和预防有着非常重要的作用. 该文旨在利用一维血流动力学模型模拟整个Willis环的流量变化和压力分布, 并建立动脉内氧输运的一维模型以模拟Willis 环内氧分压的变化规律, 为研究脑组织内血液流动和氧输运打下基础. 首先, 基于弹性圆管内的一维非线性流动方程和状态方程建立血流动力学模型, 在一维对流扩散方程的基础上, 考虑由管腔向壁面的扩散和壁面细胞的新陈代谢消耗推导出氧输运特性方程. 通过 Lax-Wendroff两步法对血流动力学方程进行离散, 而在进行对流扩散方程的离散时, 则运用迎风格式. 通过数值计算得到了正常情况下Willis环各个血管任意位置的流量、压力和氧分压的变化曲线, 正常情况下各个位置的氧分压处于稳定的平衡状态. 最后, 还通过此模型进一步模拟了右侧颈内动脉狭窄对各个血管内流动的影响. 当狭窄程度达到80%时, 中脑动脉的流量和压力会明显下降, 造成其供应区域的血流减少. 同时, Willis环右侧血管内的氧分压会大大降低, 而左侧血管的氧分压会出现上升趋势, 但幅度要小于右侧血管降低的幅度.     

多层流体的Marangoni对流

介绍微重力环境下矩形液池中多层不相混液体的Marangoni对流及其不稳定性.Marangoni对流流动的形成是由于在该系统施加一个与液体交界面相垂直的外加温度梯度.Marangoni对流的线性不稳定性分析是基于无限延长矩形液池内的两层流体系统.应用数值模拟研究了Marangoni对流的流体动力学和热传输特性,将多层流体的Marangoni对流同两层流体中外加温度梯度与流体交界面平行时引起的热毛细对流的主要特征进行了比较.     

Dynamic scaling of unsteady shear-thinning non-Newtonian fluid flows in a large-scale model of a distal anastomosis

Dimensional analysis has been applied to an unsteady pulsatile flow of a shear-thinning power-law non-Newtonian liquid. An experiment was then designed in which both Newtonian and non-Newtonian liquids were used to model blood flow through a large-scale (38.5 mm dia.), simplified, rigid arterial junction (a distal anastomosis of a femorodistal bypass). The flow field within the junction was obtained by Particle Imaging Velocimetry and near-wall velocities were used to calculate the wall shear stresses. Dimensionless wall shear stresses were obtained at different points in the cardiac cycle for two different but dynamically similar non-Newtonian fluids; the good agreement between the measured dimensionless wall shear stresses confirm the validity of the dimensional analysis. However, blood exhibits a constant viscosity at high-shear rates and to obtain complete dynamic similarity between large-scale experiments and life-scale flows, the high-shear viscosity also needs to be included in the analysis. How this might be done is discussed in the paper.     

血流动力学数值模拟与动脉粥样硬化研究进展

血流动力学因素被认为与动脉粥样硬化等病理改变密切相关。目前血流动力学数值模拟的对象,主要集中于分支动脉、弯曲动脉以及因血管内膜增生而导致的局部狭窄动脉,这些都是动脉粥样硬化多发的病灶部位。精确的血流动力学数值模拟,必须依赖于解剖精确的血管几何模型和生理真实的血流与管壁有限变形的非线性瞬态流－固耦合。只有在“虚拟血液流动”的基础上,综合考虑血管内的壁面剪应力、粒子滞留时间和氧气的跨血管壁传输等多种因素,血流动力学的数值模拟才能真正有助于人们理解动脉粥样硬化的血流动力学机理,才有可能应用于有关动脉疾病的外科手术规划中。      

A numerical study of fully developed laminar flows in pipes with two planar curvatures

The spectral element method is applied on unstructured tetrahedral elements to solve the Navier–Stokes equations for fully developed laminar flow in pipes with two planar curvatures. Specific implementations of the spectral element method to double curved pipes and parallelization are described. Previous studies on flows in pipes focused on constant curvature or torsion geometries, as well as pipes with varying curvature. This study focuses on the periodic variation of both the curvature as well as torsion by analysing a pipe having two planar curvatures. The effects of the three parameters defining the pipe are studied to isolate the curvature and torsion effect on the magnitude and angle of the secondary flow. Furthermore, the geometric effects on the wall shear stress are studied, as it is an important fluid flow property, especially in blood flows. Copyright © 2006 John Wiley & Sons, Ltd.     

一种确定均匀动脉壁面切应力的非线性方法

从Ling和Atabek提出的``局部流'理论出发，提出一种利用测量血液黏度、管轴上 的血流速度、压力和管径波形计算均匀动脉管壁切应力的非线性方法. 将这种方法与柳兆荣 等提出的利用测量血液黏度、管轴上的血流速度和平均管径计算切应力的线性方法比较，结 果表明，当管壁脉动幅度较小时，两种方法计算的压力梯度、流速剖面和管壁切应力差别较 小；而当管壁脉动幅度增大时，两种方法计算的压力梯度、流速剖面和管壁切应力差别增大. 对于小幅脉动均匀动脉，用线性方法计算管壁切应力有较高的精度；而对于大变形 均匀动脉，则需要考虑非线性因素对管壁切应力的影响. 由于作为输入量的血液黏度、轴心 血流速度、压力波形和管径波形可在活体上通过无损伤或微损伤的检测方法得到， 所提出的计算切应力的方法为在体或离体研究切应力与动脉重建的关系提供了方法学基础.     

脉动流条件下动脉狭窄血管内脂质浓度极化现象的计算机数值模拟

在脉动流条件下,用计算机数值模拟的方法对低密度脂蛋白(LDL)在动脉狭窄血管段内的质量传输进行了研究。计算结果表明．无论是在定常流还是在脉动流条件下．LDL都将聚积于血管狭窄处峰口附近的流动分离点,LDL壁面浓度在此处最高。在脉动流条件下,LDL在血流受扰动区的聚积高于定常流的值;而且．流动分离点处LDL壁面浓度峰值覆盖的区域也宽于定常流。本文所揭示出的LDL在血管狭窄处的质量传输现象可能在动脉粥样硬化的局部性和动脉狭窄的形成中起着很重要的作用。     

Reynolds stress closure applied to axisymmetric,impinging turbulent jets

A second-order, single-point closure model for calculating the transport of momentum in turbulent flows is extended to cover flows that are close to solid surfaces. In such flows the proximity of a solid boundary directly influences the fluctuating pressure field within the main body of the flow and leads to a dampening of velocity fluctuations normal to the wall. These effects are accommodated through the incorporation of an additional contribution in the modelled form of the redistributive fluctuating pressure term used in the Reynolds stress transport equation. Predictions of the extended closure model are compared with available data in configurations where an air jet impinges orthogonally onto a plane surface. The inclusion of the wall reflection model is shown to result in superior predictions of mean velocities, and normal and shear stresses. In particular, normal-to-wall velocity fluctuations and shear stresses are successfully damped resulting in agreement with observations.     

Modification of near-wall turbulence structure in a shear-driven three-dimensional turbulent boundary layer

 Most high Reynolds number flows of engineering interest are three-dimensional in nature. Key features of three-dimensional turbulent boundary layers (3DTBLs) include: non-colateral shear stress and strain rate vectors, and decreasing ratio of the shear stresses to the turbulent kinetic energy with increasing three-dimensionality. These are indicators that the skewing has a significant effect on the structure of turbulence. In order to further investigate the flow physics and turbulence structure of these complex flows, an innovative method for generating a planar shear-driven 3DTBL was developed. A specialized facility incorporating a relatively simple geometry and allowing for varying strengths of crossflow was constructed to facilitate studies where the skewing is decoupled from the confounding effects of streamwise pressure gradient and curvature. On-line planar particle image velocimetry (PIV) measurements and flow visualization results indicate that the experimental configuration generates the desired complex flow, which exhibits typical characteristics associated with 3DTBLs. Furthermore, spanwise shear results in modification of the near-wall turbulence structure. Analysis of near-wall flow visualization photographs revealed a reduction of mean streak length with increasing spanwise shear, while streak spacing remained relatively constant. In the most strongly sheared case, where the belt velocity is twice that of the freestream velocity, the mean streak length was reduced by approximately 50%. Received: 28 October 1997/Accepted: 4 February 1998     

Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels

Eccentric annular pipe flows represent an ideal model for investigating inhomogeneous turbulent shear flows, where conditions of turbulence production and transport vary significantly within the cross-section. Moreover recent works have proven that in geometries characterized by the presence of a narrow gap, large-scale coherent structures are present. The eccentric annular channel represents, in the opinion of the present authors, the prototype of these geometries. The aim of the present work is to verify the capability of a numerical methodology to fully reproduce the main features of the flow field in this geometry, to verify and characterize the presence of large-scale coherent structures, to examine their behavior at different Reynolds numbers and eccentricities and to analyze the anisotropy associated to these structures. The numerical approach is based upon LES, boundary fitted coordinates and a fractional step algorithm. A dynamic Sub Grid Scale (SGS) model suited for this numerical environment has been implemented and tested. An additional interest of this work is therefore in the approach employed itself, considering it as a step into the development of an effective LES methodology for flows in complex channel geometries. Agreement with previous experimental and DNS results has been found good overall for the streamwise velocity, shear stress and the rms of the velocity components. The instantaneous flow field presented large-scale coherent structures in the streamwise direction at low Reynolds numbers, while these are absent or less dominant at higher Reynolds and low eccentricity. After Reynolds averaging is performed over a long integration time the existence of secondary flows in the cross session is proven. Their shape is found to be constant over the Reynolds range surveyed, and dependent on the geometric parameters. The effect of secondary flows on anisotropy is studied over an extensive Reynolds range through invariant analysis. Additional insight on the mechanics of turbulence in this geometry is obtained.