人体上呼吸道内气流和燃烧污染物颗粒沉积的大涡模拟.pdfVIP

中国工程热物理学会 燃烧学 学术会议论文 编号：064206 人体上呼吸道内气流和燃烧污染物颗粒 沉积的大涡模拟 林江o’2胡桂林1樊建人2 l浙江科技学院，杭州，3100142浙江大学机械与能源工程学院． 杭州， 310027 0571 irfan@zju．edu．cn 摘要：I黪于燃烧污染物控制研究集中在进入人体之前的现状，在分析了现有研究的所用上呼吸道几 何模型的筚础I：．建带r一个更完整的从ILl腔到前三级支气管的二三维几何模型：采用大涡模拟的方 法计算了复杂几何结构内的气体流动和颗粒的运动现象。研究了入【J稳定时，休闲状态、轻微运动 和加大活动量等条件下的气流和颗粒物的运动特性，统计和分析了不同工况下燃烧污染物颗粒在呼 吸道模型内各部位的沉积率，结果表明．颗粒沉积率的大小主要受到呼吸强度与颗粒惯性的影响． 研究结果与文献中的实验数据进行j，比较，结果表明两者吻合很好。 ’ 关键词：气流颗粒沉积人体J：呼吸道模型大涡模拟 1． 引言 燃烧产生的烟气由人量的细微颗粒物组成，是大气污染的一大主要因素。其中， 可吸入颗粒物的污染对人类的影响更大，被人体吸入后可能会引起支气管炎、肺气肿、 尘肺、肺癌等疾病。长期以来对大气中可吸入颗粒物污染控制的研究集中在吸入人体 之前，而对于可吸入颗粒物进入人体后情况的研究是近年来才在国际上兴起的，国还 鲜见研究报道。而深入了解由燃烧等产生的细微颗粒物通过呼吸道进入人体后的运动 特性，了解气流和颗粒在呼吸道内的运动和沉积规律，不仅可以明晰可吸入颗粒物的 污染致病机理，而且可以为治疗准确投掷药物寻找靶点，为提高治疗效果提供可靠的 理论依据，有效地控制燃烧污染。 迄今为止，国内外研究人员有关呼吸道内气流和颗粒运动的计算流体力学模拟研 究，从计算模型看大多数采用低雷诺数k．∞模型0,21和标准k-e模型【3】描述气体流动， 拉格朗日框架下跟踪颗粒的运动【I川，采用更高精度的模型的研究结果未见报道；从采 用的几何模型米看，，“泛采用Weibel模型A气管．支气管对称性模型【41，口腔．咽．喉． 气管【I’2J和理想化口腔．咽喉模型【5J，缺乏关丁．从口腔到二级支气管的更为完整的几何模 型的研究。 因此，本文将基丁理想化口腔．咽喉模型和Weibel模型A气管．支气管模型两者的 合理性和有效性，建立一个更为完整的由口腔到前三节支气管的三维几何模型；选用 能更精确的描述人体上呼吸道流场和颗粒运动的火涡模拟方法求解气流和颗粒运动状 况，为进一步研究由燃烧等产生的颗粒物对人体的影响打F基础。 2．模型描述 2．1几何模型 图l是口腔．前三节支气管三维完整几何模型图。口腔为拱形，咽部是一个前后 668 管壁近似平行的椭圆柱体，喉部为一椭圆管，声门位于喉部突起的部位，气管-支气管 采用Weibel模型A，具体的建模方法和参数可见文献【6】。 圈1人体上呼吸道几何模型 图2关键部位的网格划分 2．2控制方程 由于在呼吸道中存在管道的急剧收缩面和扩展面，以及分叉结构等复杂几何结 构，流动将不仅有层流，而且会有复杂的湍流流动结构。基于大涡模拟方法中的亚网 格应力(SGS)模型，气相场的张量形式控制方程为： ，一、 望+业划：o (1) at axi 掣+掣=一害+毒卜“悟+鲁]]‘．c2， 式中：f、_，=l、2、3，p为流体的密度，t是时间，“是流体的平均速度矢量，工为坐 ●一 标方向，P是流体微元体上的平均压力，∥为流体粘性，湍流粘性鸬采用Smagorinsky 距格子模式模拟：。 “=(e△)2I-i1=厨