细颗粒在气管二叉管内运动的大涡模拟.pdfVIP

中国工程热物理学会 多相流 学术会议论文 编号j066154 细颗粒在气管二叉管内运动的大涡模拟 林江“2胡桂林1樊建人2 l浙江科技学院，杭州，3100142浙江大学机械与能源工程学院，．杭州．310027 0571 melinj@zju．edu．cn 摘要t本文直接从医学矸扫描得到的原始数据出发，利用图像辨识技术，．重建了一个真实人体气 管支气管的三维几何模型．采用大涡模拟的方法计算了复杂几何结构内的气体流动现象，并在拉格 朗日框架下跟踪颗粒的运动规律。得到7气流场的三维分布，以及颗粒的运动轨迹和沉积分布情况。 研究结果可为进一步分析可吸入颗粒的致病机理和气雾剂治病效果提供可靠的依据。 关键词：气流：颗粒沉积：气管支气管；大涡模拟 1、引言 大气中的悬浮小颗粒是大气污染的一大主要因素，被人体吸入后可能会引起气管 炎，支气管炎和肺气肿等呼吸道疾病；另一方面，。国内外越来越多气雾剂被用来治疗 呼吸道疾病，如小儿哮喘等。因此，深入了解气流和颗粒在呼吸道内的运动和沉积规 律，不仅可阻明晰治病机理。有效防病，而且还可以为准确投掷药物到患病部位，提 供可靠的理论依据，大大提高治疗效果。 迄今为止，国内外已有不少研究人员从事过有关呼吸道内气流和颗粒运动的计算 流体力学模拟研究。其中。大多数模型采用低雷诺数k一∞模型“。1和标准k-e模型“1描 述气体流动，拉格朗日框架下跟踪颗粒运动“～，而没有用到更高精度的模型。在计算 域几何模型方面，’以往研究广泛采用wbibel模型A气管．支气管对称性模型(见图1)“1， 口腔．咽．喉—气管￡1捌和理想化口腔．咽喉模型嘲，而直接矛Ⅱ用真实人体数据重建的几何模 型非常少见饽j。 本研究从筛选真人、医学cT扫描、得到的 原始数据图像开始，并利用图像辨识技术；重 建了一个真实的数据完整的人体气管支气管的 三维几何模型；并采用更为精确的大涡模拟方 法求解气流和颗粒运动，因此研究在几何模型 的真实性和数学物理模型的准确性方面都可望 有较大提高。 圉LWdbcI模型A气管，支气管对称性模型 2、模型描述 乏1几何模型 原始CT数据图像是由16层螺旋CT西门子SOMATOMSensationl6扫描一个肺 部有疾病的37岁男病人得到的，具体扫描参数：120kv,140naAs，每周扫描生成16 幅图像，每秒旋转两周，层厚为O．75咖。病人平躺，扫描方向为从头到脚．扫描时 告诉病人吸气到最大，摒住呼吸．以免在呼吸过程中的吸气和呼气不同阶段由于气道 的张缩运动导致图像的重影和失真。采用图像辨识技术，提取得到的三维气管支气管 几何模型如图2所示。从图中可以看到．管的形状和圆柱形相比还是存在不小的差异， 而且上下级三根管之间均呈一定的立体夹角，左支气管和右支气管是明显不对称的。 右支气管的偏角小。因此，通过扫描数据重建的模型与以往研究文献中大多采用的简 化模型有较大的差异，精确度大大提高，建立在此模型基础上的数值模拟研究也将更 具说服力，更具实际价值。 2．2控制方程 基于大涡模拟方法中的亚网格应力(SGS)模型，气相场的张量形式控制方程为： 塑+塑■。o· (1) at axi c：， 掣+掣=一拿i+oOXyk+“)(著+普]] 式中；i、，=l、2、3，，为流体的密度，t是时间，“是流体豹平均速度矢量，工为坐 标方向，p是流体微元体上的平均压力，掣为流体粘性-湍流粘性“采用Smagorinsky 亚格子模式模拟： 一=心△)2阼√2瓦i” (3) 式中：△为网格特征尺度，S“为滤波后湍流粘性应变变化率张量，Cs是模型常量， 本文中取0．1。