肺泡内颗粒运动轨迹的检测方法技术

一种肺泡内颗粒运动轨迹的检测方法，包括：根据肺泡的几何数据对肺泡进行网格建模，建立肺泡动网格模型；对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹。本发明专利技术方案提高了获取颗粒在肺泡内的运动轨迹的准确率，进而可以根据不同受力条件不同颗粒大小下颗粒的运动轨迹确定颗粒沉积规律，从而可以提高气溶胶的药性。

【技术实现步骤摘要】
肺泡内颗粒运动轨迹的检测方法
本专利技术涉及检测
，特别是涉及一种肺泡内颗粒运动轨迹的检测方法。
技术介绍
人体呼吸的主要功能是为身体的各个组织提供氧气和排除二氧化碳废气，人体的呼吸过程可以分为两个阶段：由外界环境向血液输送气体；气体经由血液进入到各个组织。随着社会、经济的不断进步，人类对生活环境的质量要求不断提高，对生存环境的保护意识逐渐增强。工业生产和环境恶化所带来的颗粒物污染已经成为评价生活质量和大气质量的一个重要指标之一。大气气溶胶颗粒物污染中，其中的部分微小的气溶胶颗粒，尤其是可吸入颗粒对人类的健康的影响更是深远，在进入人体呼吸道后，没有沉积在呼吸道的传导气管上，而是深入到人体呼吸道终末处的气体交换区沉积，很多研究表明这些颗粒对人体健康的危害最大。研究指出，人类许多疾病都和吸入颗粒物污染有着直接或间接的联系。因此，研究可吸入颗粒物在呼吸道中的运动特性，对于帮助了解可吸入颗粒的致病机理以及气溶胶治疗有着非常重要的意义，对于保护人类的健康，保护环境，提高生活质量具有积极意义。传统技术中通过建立简化的气管支气管模型，借助计算流体动力学对其中的空气流动进行数值模拟，获得颗粒运动轨迹，根据运动轨迹确定颗粒沉积规律。然而，目前技术中的模型与人体相关性较小，获得的运动轨迹与实际情况相差较大。
技术实现思路
基于此，有必要针对获取颗粒运动轨迹准确率低的问题，提供一种肺泡内颗粒运动轨迹的检测方法。一种肺泡内颗粒运动轨迹的检测方法，包括：根据肺泡的几何数据对肺泡进行网格建模，建立肺泡动网格模型；对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹。上述肺泡内颗粒运动轨迹的检测方法，通过根据实际的肺泡几何数据对肺泡进行网格建模，大大提高了模型与人体肺泡的相似度。通过对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；从而可以根据颗粒数值模拟模型确定稳态和非稳态气流下、不同受力条件下、不同颗粒直径下颗粒在肺泡内的运动轨迹。提高了获取颗粒在肺泡内的运动轨迹的准确率，进而可以根据不同受力条件不同颗粒直径下颗粒的运动轨迹确定颗粒沉积规律，从而提高气溶胶的药性。附图说明图1为本专利技术肺泡内颗粒运动轨迹的检测方法实施例的流程示意图；图2为本专利技术实施例中肺泡动网格模型的正面示意图；图3为本专利技术实施例中肺泡动网格模型的状态示意图；图4为呼吸过程中非稳态时气流流量示意图；图5为肺泡囊的二维几何模型示意图；图6为本专利技术应用实例中不同受力条件下0.1μm颗粒在肺泡模型内稳态气流下的运动轨迹示意图；图7为本专利技术应用实例中不同受力条件下1.0μm颗粒在肺泡模型内稳态气流下的运动轨迹示意图；图8为本专利技术应用实例中不同受力条件下10μm颗粒在肺泡模型内稳态气流下的运动轨迹示意图；图9为本专利技术应用实例中在完整力作用下各个时间间隔中0.1μm颗粒所在位置示意图；图10为本专利技术应用实例中在重力作用下各个时间间隔中0.1μm颗粒所在位置示意图；图11为本专利技术应用实例中在失重作用下各个时间间隔中0.1μm颗粒所在位置示意图；图12为本专利技术应用实例中在完整力作用下各个时间间隔中1.0μm颗粒所在位置示意图；图13为本专利技术应用实例中在重力作用下各个时间间隔中1.0μm颗粒所在位置示意图；图14为本专利技术应用实例中在失重作用下各个时间间隔中1.0μm颗粒所在位置示意图；图15为本专利技术应用实例中在完整力作用下各个时间间隔中10μm颗粒所在位置示意图；图16为本专利技术应用实例中在重力作用下各个时间间隔中10μm颗粒所在位置示意图；图17为本专利技术应用实例中在失重作用下各个时间间隔中10μm颗粒所在位置示意图。具体实施方式下面结合实施例及附图对本专利技术作进一步详细说明，但本专利技术的实施方式不限于此。如图1所示，为本专利技术肺泡内颗粒运动轨迹的检测方法实施例的流程示意图，包括步骤：步骤S101：根据肺泡的几何数据对肺泡进行网格建模，建立肺泡动网格模型；肺泡的几何数据可以是指某一节肺泡的几何数据。所述几何数据可以包括：管长参数，内腔直径，肺泡直径，两肺泡间的中心距离，肺泡开口角，呼吸参数。比如，可以是第18节肺泡。几何数据可以包括：管长为740μm，内腔直径为320μm，肺泡直径为320μm，两肺泡间的中心距离为340μm，肺泡开口角为120°，呼吸参数为0.283。本实施例还采用动网格模型模拟肺泡，提高模拟的真实性和相关性。肺泡动网格模型可以为3D模型。步骤S102：对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；进行空气流场模拟的方法有很多种，比如可以采用Fluent软件进行空气流场模拟。可以采用SIMPLE(Semi-ImplicitMethodforPressure-LinkedEquations)方法、SIMPLEC方法或PISO(Pressure-ImplicitwithSplittingofOperators)等方法计算空气流场的速度。步骤S103：根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；颗粒数值模拟模型是一种数学模型，是可以用于计算颗粒在肺泡内运动轨迹的模型。步骤S104：根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹。不同受力条件可以包括完整作用、失重作用和重力作用等情况。不同颗粒可以是指直径不同的颗粒。比如，直径分别为0.1μm，1.0μm和10μm的颗粒。在稳态气流测试时，可以将这些颗粒均匀的从肺泡的入口随气流一起射入肺泡内。其中，颗粒的入口速度可以与空气流场速度一致。本实施例通过根据实际的肺泡几何数据对肺泡进行网格建模，大大提高了模型与人体肺泡的相似度。通过对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；从而可以根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹。提高了获取颗粒在肺泡内运动轨迹的准确率，进而可以根据不同受力条件不同大小颗粒的运动轨迹确定颗粒沉积规律，从而提高气溶胶的药性。在其中一个实施例中，所述颗粒数值模拟模型包括：其中，其中，Xp表示t时刻颗粒在空气流场中的位置，up表示颗粒的速度，u表示空气流场的速度，FD表示作用在颗粒上的拽力，ρp表示颗粒的密度，ρ表示空气流场的密度，gi表示重力矢量，FB表示作用于颗粒上的布朗力，μ表示空气流场的粘度，dp表示颗粒的直径，CC表示滑移修正系数，λ表示分子的平均自由程，表示一个具有零平均值和单位方差的高斯随机矢量，S0表示白噪声过程的幅度，v表示空气的动力粘度，kB＝1.38×10-16，T表示空气流场的温度。本实施例中的FD充分考虑了微小颗粒下，分子间的滑移作用。通过本实施例中的颗粒数值模拟模型，可以获得各个时刻、不同颗粒大小、不同受力条件下颗粒的运动轨迹。在其中一个实施例中，具体介绍了一种获得空气流场的速度的方法。即步骤S102，包括：A1：设置肺泡动网格模型初始化的空气流场本文档来自技高网...

【技术保护点】
一种肺泡内颗粒运动轨迹的检测方法，其特征在于，包括：根据肺泡的几何数据对肺泡进行网格建模，建立肺泡动网格模型；对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹。

【技术特征摘要】
1.一种肺泡内颗粒运动轨迹的检测方法，其特征在于，包括：根据肺泡的几何数据对肺泡进行网格建模，建立肺泡动网格模型；对所述肺泡动网格模型中进行空气流场模拟，获得空气流场的速度；根据肺泡的空气流场的速度对肺泡的颗粒进行模拟，获得颗粒数值模拟模型；根据颗粒数值模拟模型确定稳态和非稳态气流下不同受力条件下不同颗粒在肺泡内的运动轨迹；在建立肺泡动网格模型过程中，还包括：建立肺泡囊的二维几何模型，用于模拟呼吸过程中的肺泡的扩张和收缩，其中，将前后开口的圆柱管形模拟呼吸道中的气管，将连接在圆柱壁面上的圆形模拟呼吸道中的肺泡，肺泡位于该段气管的正中位置；在建立肺泡囊的二维几何模型过程中，还包括设置肺泡和气管的各个特征尺寸为随时间变化的正弦函数：其中，L(t)表示肺泡或气管的特征尺寸的值，L0表示几何模型中的各个特征尺寸在吸气开始时或呼气结束时的值，f表示呼吸的频率，f＝2π/T，T表示呼吸的周期时间，β表示各个特征尺寸的扩张幅度，β＝(C+1)1/3-1，C＝(Vmax-Vmin)/Vmin，Vmax和Vmin表示肺泡动网格模型的最大体积和最小体积，所述最大体积和最小体积分别对应在实际呼吸过程中呼气终了和吸气终了时的体积；在建立肺泡囊的二维几何模型过程中，还包括：根据呼吸过程中进出肺泡的气体流量与流过肺泡所处气管的气体流量的比值确定肺泡在呼吸道上所处位置，根据所述位置建立肺泡囊的二维几何模型，其中：Qa(t)表示呼吸过程中进出肺泡的气体流量，Ra表示肺泡半径，Qd(t)表示流过肺泡所处气管的气体流量，Vt表示将经过该段气管后的所有气管及气管上所有肺泡的体积求和获得的值。2.根据权利要求1所述的肺泡内颗粒运动轨迹的检测方法，其特征在于，所述颗粒数值模拟模型包括：其中，其中，Xp...

【专利技术属性】
技术研发人员：李德波，徐齐胜，刘亚明，
申请(专利权)人：广东电网公司电力科学研究院，
类型：发明
国别省市：广东;44