论文部分内容阅读
气道是人类机体与外界环境之间气体交换所必须的通道,理解生理和病理状态下的气道的通气功能具有重要科学与临床意义。目前,利用多排CT,一次屏气扫描即可获得全肺高分辨率薄层解剖结构图像。然而,多排CT获取的气道结构数据无法明确反映高级生理(通气)功能。计算流体动力学(Computational fluid dynamics, CFD)是一种利用计算机和离散化求解流体力学问题的方法,可用于研究空气在人体气道中的流动情况。本文利用CFD方法研究了人体气道内的空气流动。首先,构建了理想、三维、二分叉、对称气道树模型,研究了模型中的速度和压力分布,考察了相似比、分叉角度和支气管平面翻转角度对空气流动的影响。然后,针对4名左肺动脉吊带(Left pulmonary artery sling, LPAS)患者,从CT图像中提取气管和主支气管,建立结构模型,根据体重和入口面积确定个性化的入口流速,从而计算出LPAS患者气管和支气管中的空气流速、壁面压力、壁面剪切应力的分布。另外,利用单向流固耦合方法,研究了流动作用下LPAS患者气道壁面的变形和应力分布。最后,基于CT图像提取的多级气道(128个出口)结构模型,研究了瞬态情况下的空气流动。主要结果包括:(1)在理想气道树模型中速度和压力分布一般是对称的,但从第四级支气管开始不均匀;内下支气管内的流速小于外上支气管;改变分叉角度和支气管平面翻转角将提高压降。(2)对于LPAS患者,由于气管和(或)主支气管的狭窄,压降可达78.9-914.5 Pa,远远高于正常控制受检者(0.7 Pa);左右主支气管的质量流量比与出口面积不相关,可能的原因是C形气管便于空气流入左主支气管(在惯性力的作用下);与正常对照相比,LPAS患者气道中的速度、壁面压力和壁面剪切应力分布更加不均匀,在狭窄处可观察到高流速、低壁面压力和高壁面剪切应力。(3)LPAS患者气道的壁面变形和等效应力是正常对照者的50-900和90-100倍,气管移动的方向与气管形状有关,与净支反力相反;速度和壁面厚度对变形和应力的影响是非线性的。(4)对于本研究中的健康对象,气道中各时相流线都比较光滑顺畅,没有发现明显涡流;气道通气效率非常高。吸气相的最大压降等于22.85 Pa,呼气相的最大压降等于7.59 Pa;壁面压力分布相对均匀,但吸气相和呼气相的规律不同;壁面剪切应力非常小,在吸气相和呼气相的峰值基本相等,约为3.02 Pa;吸气相和呼气相中气体混合的模式不同;吸气相流量百分数与呼气相流量百分数相差不大,与各肺叶体积比成正比。可见,基于CT影像提取的气道结构模型结合个体化的入口边界条件,CFD可以提供更多关于气道内空气流动的定量化信息,包括速度、压降、壁面压力、壁面剪切应力、气道壁面变形等。这些定量化信息将能帮助人们认识生理和病理情况下的气道通气功能,通过与结构信息相关联,从而揭示不同呼吸道疾病的发病机理。