2模型建立与网格划分-中国力学学会.docVIP

第八届全国流体力学学术会议 2014年9月18～21日 甘肃兰州 文章编号： CSTAM2014-B01-0087 标题：基于CFD血管模型的模拟研究 作者：周俊杰，李雪丽，房全国 单位：郑州大学化工与能源学院Copyright ? 2014 版权所有 中国力学学会地址: 北京市北四环西路15号　邮政编码:100190 　Address: No.15 Beisihuanxi Road, Beijing 10化工与能源学院，河南省郑州市高新技术开发区科学大道）化工与能源学院河南省郑州市高新技术开发区科学大道运用流牛顿流体模型pa·s、5mpa·s、6mpa·s、7mpa·s)的全血粘度血管进行了数值模拟，模拟得出血管的流场分布。全血液粘度增高会引起血流阻力增加，使血流速度减慢，最后导致血流停滞，直接影响脏器血液供应，导致疾病。通过对肿瘤血管和正常血管的模拟对比，由图1可知肿瘤血管的剪切力较正常血管大和流速的不稳定对血液流动造成了影响，同时还会引发其他疾病，肿瘤的治疗有积极地意义。肿瘤 的的根本因素大多难以预测。血管内压力场及的影响因素具有积极地意义方程[6]主要是连续方程和动量方程（N-S方程）： （1）连续性方程 式中，u、v、w表示速度在x、y、z方向上的分量。 （2）动量方程 式中，xi（i=1,2,3）为坐标系坐标，xj(j=1,2,3)为坐标系坐标，ui(i=1,2,3)为沿i方向的速度分量，fi为沿i方向的质量力，P是静压，ρ为空气密度。 建立与网格划分 直径为mm的球体。 3边界条件及求解参数的设置Velocity-inlet，出口为Outflow，其他边界默认为Wall，计算区域设为fluid。 Velocity-inlet中速度选择0.03m/s，湍流指定方法选择水利直径和湍流强度。入口处的直径为人造血管的内径，所以设定水力直径为8mm。湍流强度为湍流脉动速度与平均速度的比值，等于0.16倍的雷诺数的负八分之一次方，计算得湍流强度为4.5%。 3.2求解参数设置 （1）控制方程采用连续方程和动量方程，而不考虑能量方程。（2），密度恒定，密度ρ=kg/m3。 （3）求解器选择COUPLED IMPLICIT SOLVER。 （4）为了提高计算精度，差分格式采用二阶迎风格式。 （5）速度场和压力场的求解基于Simple算法，松弛因子默认。 （6）中选用层流非牛顿流体结构，壁面区域采用标准壁面函数。 （7）残差设为e-。 图2网格划分 4模拟结果与分析 独立性考核 s的直血管进行模拟计算，得到X=0，Y=0，Z从0变化到0.06m（其中坐标系见图1）时的速度分布及X=0，Y=2mm，Z从0变化到0.06m时的速度分布，如图3、图4所示。 由图3，图4可以看出，17万网格与20万网格的计算结果相差很小，我们认为17万网格已满足计算精度的要求，可以获得与网格无关的解。本文在以后的计算中都选择17万网格。  图3 X=0，Y=0，Z从0到0.06m的速度分布 图4 X=0，Y=2mm，Z从0到0.06m的速度分布  4.2结果分析模拟计算s，5mpa·s，6mpa·s， 7mpa·s。其他参数和图1中的模型参数一样，网格采用考核过的17万网格系统，将模型导入fluent中求解得到血管的流场分布。 流场分析 （a）粘度为4 mpa·s ）粘度为5 mpa·s （c）粘度为6 mpa·s （d）粘度为7 mpa·s 图5不同粘度下X=0处的速度云图 血管内沿Z轴方向速度逐渐增大，管中心到壁面速度逐渐减小。粘度不同的血管内血液的最大流速基本相同，但速度发展的位置，粘度高的血液要迟于粘度低的血液。血液粘度主要来源于血液内分子或颗粒间的摩擦，当温度不变时，全血的粘度主要决定于血细胞的高低，血液的粘度主要决定于血浆蛋白的含量，越低，易形成湍流。7]。4种粘度下管内湍流脉动动能最大值见表1。表1不同粘度管内湍流脉动动能最大值 粘度（mpa·s） 4 5 6 7 湍流脉动动能最大值Kmax(m2/s2) 0.537 0.376 0.255 0.301 由表1可以看出，粘度为6mpa·s、 7mpa·s 的湍流脉动动能较小，此时流动，流动属于非牛顿流体流动，血液的惯性力大于粘滞力的作用，就