中等严重程度冠状动脉病变模型的血流动力学参数分析

基于中等严重程度冠状动脉病变模型,应用流固耦合方法数值研究了中等严重程度面积狭窄率(AS=50%,65%,75%)和病变长度(LL= 0 mm,15 mm,20 mm) 对血流动力学参数的影响.研究发现:随着AS与LL的增大,病变血管分支的壁面剪应力变化愈加剧烈,狭窄段下游的壁面剪应力值逐渐降低,狭窄段下游回流区的长度呈"S"型增长,模型最大剪切速率呈抛物线型增长, 压力分布曲线显著下降.血流动力学参数结果表明, 中等严重程度面积狭窄率和病变长度均是可能引发血栓的因素,临床上应予以重视.      

基于医学影像的血流动力学分析

随着计算机成像技术的发展,计算机断层扫描(CT)和核磁共振(MRI)已经广泛应用于临床,特别是血管类疾病的诊断,比如动脉狭窄、血管瘤、血管畸形等.除了可以提供高分辨率的静态图像,先进的MRI技术还可以通过时间序列直接反映血流动力学的变化.而基于计算机成像和三维重建技术,血流动力学的参数又可以通过计算流体力学的方法进行详细的分析.如何将血流动力学参数和临床诊断相结合是近年来在转化医学领域研究的热点.文中将结合文献调研和作者自己的研究工作对基于医学图像的血流动力学分析进行综述,并探讨未来的研究方向.     

颈动脉分支的血流动力学数值模拟

采用有限元法数值模拟颈动脉分支的血流动力学。根据在体测量的实际尺寸来构造颈动脉分支的几何模型,以保持模型的解剖精确度;利用在体测量的颈内动脉和颈外动脉流量波形以及主颈动脉的压力波形来确定数值计算的边界条件,以保持数值计算的生理真实性。关注的重点是颈动脉窦内的局部血流形态、二次流和壁面剪应力。在心脏收缩的减速期和舒张期的某些时刻,颈动脉窦中部外侧壁面附近产生了流动分离,形成了一个低速回流区。该流动分离是瞬态的,导致了壁面剪应力的振荡,其振荡范围在-2～6dyn／cm^2之间。同时,颈动脉窦中部横截面内的二次流存在于整个心动周期,最大的二次流速度为同时刻轴向速度平均值的1／3左右。     

颅内Willis环二维血液动力学计算分析

采用颅内二维Willis环模型,对血管在不同程度狭窄情况下的血液流动进行稳态计算与分析,计算中采用了进出口处的压力边界条件和无滑移固壁条件.计算得到完整环和交通动脉缺失情况下各进出口血管流量变化曲线.所获结果表明,完整Willis环中,前交通动脉的代偿能力在左颈内动脉狭窄率达到25%以后才能得到比较明显的体现;后交通动脉缺失情况下前交通动脉的代偿能力变弱;Willis环的前部和后部有相对独立性,只有在主要进口动脉严重阻塞的时候,单侧的独立性才丧失.     

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

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

主动脉瓣倾斜角度血流动力学的 PIV 实验研究

瓣叶血栓是主动脉瓣置换术后典型的继发性瓣膜疾病,血流动力学特征异常在其发展过程中至关重要.本文利用粒子图像测速 (particle image velocimetry,PIV) 系统,实验研究了主动脉瓣开口纵向轴线与升主动脉纵向轴线之间倾斜角度 ($\alpha =0^\circ$, $\alpha=5^\circ$,$\alpha =10^\circ$ 和 $\alpha =15^\circ$) 对速度、涡度和黏性剪应力分布等血流动力学特性的影响.研究结果表明:当 $\alpha =0^\circ$ 时,主动脉根部跨瓣血液流动为中心对称流动,而 $\alpha =5^\circ$,$\alpha=10^\circ$ 和 $\alpha =15^\circ$ 时跨瓣血液流动向升主动脉的左冠状动脉一侧倾斜.随着倾斜角度增大,跨瓣血液流动方向倾斜程度增加,血液流动冲击升主动脉壁,损伤内皮细胞导致血栓形成.主动脉瓣倾斜时主动脉窦血液流动速度增大,涡旋也更向主动脉窦底部运动,不利于血液从冠状动脉口流出向心肌供血.同时,主动脉根部的高涡度和高黏性剪应力区域也向升主动脉的左冠状动脉一侧倾斜,主动脉窦的高涡度区域位于主动脉窦底部、高黏性剪应力区域分布于主动脉窦壁面处.主动脉瓣存在倾斜角度时,涡度和黏性剪应力较大,特别是 $\alpha =10^\circ$ 和 $\alpha=15^\circ$,为血栓形成提供了有利环境.研究结果可为临床主动脉瓣置换术参数选择以及继发性瓣膜疾病的避免提供理论依据和技术参考.      

计算流体动力学和计算水动力学的进展和展望

从计算流体动力学与计算水动力学的三个主要方面──数学模型研究、数值模拟方法研究和应用计算机软件研制谈其发展状况和进展。     

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

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

颅内动脉瘤血液动力学和破裂机理的数值模拟进展

通过有限元法数值模拟来研究颅内动脉瘤形成、发展和破裂 与其血液动力学参数的关系,对于临床医学特别是介入疗法具有重 要的指导意义。本文对于近年国外出现的各种计算模型及其主要结 果做一综述。     

不同心排出量下主动脉瓣血流动力学的PIV实验研究

主动脉瓣发生病变时导致心排出量(cardiac output, CO)减少,而心排出量减少与主动脉瓣血流动力学耦合作用, 引发瓣膜继发性疾病.本文基于医学影像数据三维重构带有冠状动脉的主动脉根部,制备高度光滑和透明的主动脉根部实验模型, 构建体外脉动循环模拟系统,利用粒子图像测速技术(particle image velocimetry,PIV)研究冠状动脉存在时心排出量对主动脉瓣速度分布、黏性剪应力(viscous shear stress, VSS)和雷诺剪应力(Reynolds shear stress, RSS)等血流动力学的影响.研究结果表明: 冠状动脉的存在改变了主动脉窦中的涡旋运动和涡度,冠状动脉存在时流体经由冠状动脉流出, 主动脉窦中的涡旋运动逐渐消失,涡度较早开始减小. 峰值期, 中心对称流动两侧区域存在正、负高黏性剪切区域,存在冠状动脉一侧的升主动脉下游存在高雷诺剪应力区域.心排出量显著影响主动脉瓣的速度分布、VSS和RSS等血液流动和受力状况.随着心排出量增大, 冠状动脉存在时峰值期的最大速度、VSS和RSS增大, 即$CO=2.1$, 2.8, 3.5和4.2 l/min时, 最大速度分别为0.98, 1.13, 1.21和1.37 m/s, 最大VSS分别为0.87, 0.95, 0.96和1.02 N/m$^{2}$, 最大RSS分别为103.76, 116.25, 138.68和146.55 N/m$^{2}$. 心排出量较低时,主动脉瓣较低的跨瓣流动速度和黏性剪应力易导致血栓形成,研究结果可为主动脉瓣置换术提供理论参考.      

Medical image based hemodynamic analysis

随着计算机成像技术的发展,计算机断层扫描(CT)和核磁共振(MRI)已经广泛应用于临床,特别是血管类疾病的诊断,比如动脉狭窄、血管瘤、血管畸形等.除了可以提供高分辨率的静态图像,先进的MRI技术还可以通过时间序列直接反映血流动力学的变化.而基于计算机成像和三维重建技术,血流动力学的参数又可以通过计算流体力学的方法进行详细的分析.如何将血流动力学参数和临床诊断相结合是近年来在转化医学领域研究的热点.结合文献调研和作者自己的研究工作对基于医学图像的血流动力学分析进行综述,并探讨未来的研究方向.     

Hemodynamic analysis of intracranial aneurysms using phase-contrast magnetic resonance imaging and computational fluid dynamics

Additional hemodynamic parameters are highly desirable in the clinical management of intracranial aneurysm rupture as static medical images cannot demonstrate the blood flow within aneurysms. There are two ways of obtaining the hemodynamic information—by phase-contrast magnetic resonance imaging(PCMRI) and computational fluid dynamics(CFD). In this paper, we compared PCMRI and CFD in the analysis of a stable patient's specific aneurysm.The results showed that PCMRI and CFD are in good agreement with each other. An additional CFD study of two stable and two ruptured aneurysms revealed that ruptured aneurysms have a higher statistical average blood velocity,wall shear stress, and oscillatory shear index(OSI) within the aneurysm sac compared to those of stable aneurysms.Furthermore, for ruptured aneurysms, the OSI divides the positive and negative wall shear stress divergence at the aneurysm sac.     

Flow effects of blood constitutive equations in 3D models of vascular anomalies

The effects of different blood rheological models are investigated numerically utilizing two three‐ dimensional (3D) models of vascular anomalies, namely a stenosis and an abdominal aortic aneurysm model. The employed CFD code incorporates the SIMPLE scheme in conjunction with the finite‐volume method with collocated arrangement of variables. The approximation of the convection terms is carried out using the QUICK differencing scheme, whereas the code enables also multi‐block computations, which are useful in order to cope with the two‐block grid structure of the current computational domain. Three non‐Newtonian models are employed, namely the Casson, Power‐Law and Quemada models, which have been introduced in the past for modelling the rheological behaviour of blood and cover both the viscous as well as the two‐phase character of blood. In view of the haemodynamical mechanisms related to abnormalities in the vascular network and the role of the wall shear stress in initiating and further developing of arterial diseases, the present study focuses on the 3D flow field and in particular on the distribution as well as on both low and high values of the wall shear stress in the vicinity of the anomaly. Finally, a comparison is made between the effects of each rheological model on the aforementioned parameters. Results show marked differences between simulating blood as Newtonian and non‐Newtonian fluid and furthermore the Power‐Law model exhibits different behaviour in all cases compared to the other models whereas Quemada and Casson models exhibit similar behaviour in the case of the stenosis but different behaviour in the case of the aneurysm. Copyright © 2005 John Wiley & Sons, Ltd.     

浅谈动脉流体动力学分析

就我们所熟知,绝大部分正常动脉流,其血液的流动特性是属于层流范围,但随着弯曲和分支部分会产生血液流之二次回流区,进而形成所谓近似非稳态流及紊流. 因此动脉流体的特性会随动脉外形及条件的改变而改变. 在某些情形下,异常动脉的血液动力特性会造成动脉的病变. 因此,近年来动脉血液流体的特性的研究,常着重于异常动脉的血液动力特性所形成剪应力和病变部位动脉粥状硬化关系的探讨.动脉血液流动经常包含分离流或二次回流运动,而这是流体力学的分析或数值模拟最困难的部分. 有关分离流或二次回流的研究包括正常血管流和异常血管流,藉由二次回流的模拟与测量可以观察血管病变的形成与演变,其中最受注目探讨题目是窄化血管如粥状斑块相关的血液流动分析.将回顾二维和三维、稳态、非稳态之动脉血流与窄化血管相关的几何外形作模拟研究和实验. 并提供对血液动力学的研究方向,以作为未来医疗诊断与发展相关器材之参考.     

The effect of aneurismal-wall mechanical properties on patient-specific hemodynamic simulations： two clinical case reports

Hemodynamic factors such as the wall shear stress play an important role in the pathogenesis and treatment of cerebral aneurysms. In present study, we apply computational fluid-structure interaction analyses on cerebral aneurysms with two different constitutive relations for aneurismal wall in order to investigate the effect of the aneurismal wall mechanical properties on the simulation results. We carry out these analyses by using two patient-specific models of cerebral aneurysms of different sizes located in different branches of the circle of Willis. The models are constructed from 3D rotational angiography image data and blood flow dynamics is studied under physiologically representative waveform of inflow. From the patient models analyzed in this investigation, we find that the deformations of cerebral aneurysms are very small. But due to the nonlinear character of the Navier-Stokes equations, these small deformations could have significant influences on the flow characteristics. In addition, we find that the aneurismal-wall mechanical properties have great effects on the deformation distribution of the aneurysm, which also affects the wall shear stress distribution and flow patterns. Therefore, how to define a proper constitutive relation for aneurismal wall should be considered carefully in the hemodynamic simulation.     

基于地形图像特征与模糊决策的海底地形可导航性分析方法

基于多波束测深数据和图像匹配技术的水下地形辅助导航方法需要对实施区域的地形可导航性进行评估,其中的地形独特性计算需要提取图像特征,并利用各特征参数进行综合评判。地形图像的直方图包含实际地理区域的起伏程度信息,通过定义直方图复杂度对此信息进行了有效提取。均方根对比度和拉普拉斯模方和分别侧重描述了地形整体和局部的变化程度,图像的二维熵从信息论的角度反映了地理区域包含的信息量。此外,依据单一特征参数对于待评价地形的可导航性进行排名,其结果存在不一致性,利用模糊综合评判对各参数的信息进行了有效综合。通过实测多波束地形数据仿真分析,基于图像纹理特征的地形匹配方法,在模糊评判结果最优的地形区域能够将匹配误差限制在2 m以内,结果验证了该可导航性分析方法的有效性。     

Study on hemodynamics in patient-specific thoracic aortic aneurysm

The objective of this study is to investigate the hemodynamics in patient-specific thoracic aortic aneurysm and discuss the reason for formation of aortic plaque. A 3-Dimensional pulsatile blood flow in thoracic aorta with a fusiform aneurysm and 3 main branched vessels was studied numerically with the average Reynolds number of 1399 and the Womersley number of 19.2. Based on the clinical 2-Dimensional CT slice data, the patient-specific geometry model was constructed using medical image process software. Unsteady, incompressible, 3-Dimensional Navier-Stokes equations were employed to solve the flow field. The temporal distributions of hemodynamic variables during the cardiac cycle such as streamlines, wall shear stresses in the arteries and aneurysm were analyzed. Growth and rupture mechanisms of thoracic aortic aneurysm in the patient can be analyzed based on patient-specific model and hemodynamics simulation.     

Soil flow analysis for grouser wheels based on a particle image velocimetry method

This paper presents an analysis, based on a particle image velocimetry method, of soil flow field beneath a grouser wheel traveling over loose soil. Although the grouser wheel is expected to have better traction and mobility over fine, loose soil, its interaction mechanisms with the soil remain to be elucidated. Thus, a particle image velocimetry-based soil flow analysis is conducted to directly observe soil behavior around the grouser wheel. In the experimental analysis, key parameters of the soil flow field, such as general shape, thickness, streamlines of the flow field, soil velocity on the streamlines, and soil failure angle are examined quantitatively. From the results, the soil flow shape periodically changes with wheel rotation, and this change appears, depending on wheel slip varying over time. Furthermore, the experimental result of the soil failure angle differs drastically from its typical theory. These results will contribute to modeling the mechanical interaction between the grouser wheel and soil.     

基于螺旋桨滑流效应的大展弦比机翼气动弹性分析

以高空长航时大展弦比太阳能无人机机翼为研究对象,针对分布式电驱螺旋桨滑流和大展弦比机翼之间耦合的复杂气动干涉问题,采用滑移网格方法、动网格技术、SST k-ω RANS湍流模型和CFD/CSD （Computational Fluid Dynamics/Computational Structural Dynamics）双向流固耦合技术,研究了螺旋桨不同转速、布局方式和气动阻尼对机翼气动弹性响应的影响。数值计算结果表明,螺旋桨滑流会改变机翼表面的压力分布;螺旋桨流场对机翼的扰动频率接近机翼的结构固有频率时,机翼会发生共振;螺旋桨的位置越靠近翼尖,或螺旋桨的数量增多,都将增加机翼气动弹性响应的幅值。