用晶格玻尔曼方法研究血液在弹性管中的流动.pdf

● 用晶格玻尔兹曼方法研究血液在弹性管中的流动 摘要 Boltzmann 90年代初，基于单弛豫时间近似的晶格玻尔兹曼方法(1attice method，简写为I。BM)被提了出来，该方法具有相互作用的完全局域性一适合于 并行计算、边界条件易于施加、易于编程的优点，同时又去除了诸如统计噪音、 ● 流体压强依赖于速度等缺陷。该模型因而被广泛应用于研究湍流、两相流、反应 扩散问题、颗粒流和悬浮体等物理与化学系统，取得了令人信服的结果。血液流 动系统是由心脏、血液和血管构成的边界特别复杂的生物系统，从机制上讲血流 系统主要受到两个因素控制，即心脏收缩所产生的动脉压力和由血管及血液本身 的物理性质对血流所产生的作用力。本文集中对后一个因素进行研究，用LBM对 血流系统的弹性动脉内血液流动性质及其血管壁运动产生的反作用进行初步探 ● 索。我们的工作主要包括以下几个方面 (1)研究血管壁的弹性运动对血管流量的影Ⅱ向。把血管作为极限纵向强约束 弹性管，对血管壁进行离散化处理，从管壁物性方程出发，研究血管壁在流体压 力、弹力和弯曲力作用下的径向运动对血液流动的影响。我们发现主要有两个因 素影响血流量，一是血液的粘滞性，二是血管本身的弹性，而血管弹性的变化起 ● 着至关重要的作用。在数值模拟中我们发现血管弹性系数存在一个临界值QC 大于该值血流量对血管弹性不敏感，小于该值血流量随弹性系数减小而显著增 加，这为改善血管柔性从而增加身体供血功能提供了理论依据。 (2)研究血管壁发生局部病变时血管瘤形成的机理，血管瘤对血液流的反作 用。通过假设血管壁某一段的弹性系数变小，从而需要增大局部形变以获得与流 体压力的平衡，再现了血管瘤的形成过程。在d、Reynolds数下，利用张量积分法 对形成后的血管瘤壁的剪切应力进行计算，表明在血管瘤内由于回流涡旋的产 生，作用于血管瘤壁的剪切应力大大减小，而在血管瘤的近端、远端与正常血管 结合处剪切应力出现峰值，这一减一增将可能导致血管瘤内血栓的形成以及接口 处血管的破裂。这些结果与用其它数值方法所得结论是一致的。 (3)研究血液流的热性质，建立13速热格子玻尔兹曼模型的宏观方程。在该 ● 模型中定义了内能，可用来研究不同温度下血液的粘滞流动。在热传导方程中弓 104的上边界驱动的空腔流模拟的数值结果表明了该方法的精确性和可靠性。当 空腔的左右边界存在温差时，将出现热对流现象，在该模型中引入简化的浮力作 用，较好地模拟了Rayleigll数介于103～106的空腔热对流。该方法与其它使用多 驰豫时问或多个配分函数的玻尔兹曼方法相比，在计算热流体方面表现出编程简 单、占用更少计算机CPU时间和内存等优点，可望在血液流的热性质研究上有较 广阔的应用前景。 关键词：晶格玻尔兹曼方法，血液流，血管壁，血管瘤，剪切应 力，空腔流，热对流。 ． ● 2 ● Boltzmannmethodsimulationsofblood in Lattice flowing