基于光滑耗散粒子动力学的血流模拟

来源 :吉林大学 | 被引量 : 0次 | 上传用户:pdswzjhxr
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
血液作为运输物质和氧气的主要载体,对生命维持发挥着至关重要的作用.它也储存着人体的健康信息,许多疾病都与血液有着密切的联系,如疟疾、贫血、血栓等,因此验血成为了很多疾病诊治的先导手段.对血液的研究不仅有助于深入理解血液循环的机理,而且也有利于研发与血液相关疾病的诊治手段.本论文采用数值方法对血流进行模拟,主要研究了血浆的微流控分离技术以及红细胞压积的时空不均匀性.考虑到血管网络的复杂性以及细胞的大小尺度,我们采用基于粒子的介观数值方法—光滑耗散粒子动力学,来描述微血管网络内流体的行为;同时采用浸入边界法来处理细胞与流体间的相互作用.这两者结合便是本论文所采用的主要数值方法(SDPD-IBM),它不仅便于解决复杂计算区域内的流固耦合问题,而且还能在一定程度上描述介观尺度的随机振荡现象.然而,该方法也有一些不足之处,主要包括进出口边界条件的实现和计算效率偏低两个问题.基于此,本论文主要贡献有:1)开发了精确的进出口边界条件,2)构建高效的并行框架,3)数值设计血浆分离的微流控芯片,4)数值研究红细胞压积的时空不均匀性.不同于基于网格的数值方法,基于粒子的方法一般采用Lagrange描述,在进口处需要有粒子不断流入,而在出口处需要不断删除流出的粒子.因此进出口边界条件是基于粒子类方法的重要挑战之一.为此,我们在计算区域两端分别增加一个进口区域和一个出口区域.进口区域是一个可以不断产生新的粒子和细胞的生成区域.在该区域内,我们采用周期边界条件;当一个粒子或细胞流过边界时,该粒子或细胞被放回该区域.为了保证粒子和细胞能够不断地流入,该粒子或细胞同时也被复制到计算区域.出口区域是一个处理流出粒子和细胞的删除区域.当一个粒子或细胞流出该区域时,若直接被删除,会导致系统质量和动量不守恒.因此,为了保证质量和动量守恒,我们对出口区域内的每个粒子施加保守力和耗散力自适应地补偿质量和动量的损失.模拟结果显示,我们所提出的进出口边界条件可以准确地实现多出口的血流模拟.为了保证计算精度,基于粒子的方法往往采用大量的粒子来离散计算区域,这样计算效率便成为其另一个重要的挑战.为此,我们构建了一个具有主从结构的MPI并行框架.主线程负责划分任务并且接收计算结果,而从线程分为三种类型:细胞线程、进口线程和混合线程.一个细胞分配一个细胞线程,用来计算细胞的力学行为,如细胞的变形和细胞间的聚集,这样大大减少了细胞线程间的信息传递量.进口线程和混合线程采用区域分解进行划分,分别负责更新进口区域和计算区域的细胞和粒子的位置和速度等物理量.对于一个典型的血流模拟算例,约60000个粒子和40个细胞,该并行框架在计算时间上比串行方式提升了近30倍.血浆分离是实验分析和制药过程中常见的操作之一.我们通过对实验设计的一款血浆分离的微流控芯片进行分析研究,数值设计了一个更加高效的血浆分离芯片.首先,我们将实验芯片拆分为单个分叉管道,研究管道流量及分叉角度对分离纯度和效率的影响;确定了一组最优参数,即分叉角度为2π/3,细胞与血浆出口的流量比为2.5:1.然后,以最优参数组合设计了新的血浆分离芯片,并与实验芯片分离效果进行比较.结果表明,在两种芯片都没有细胞堵塞的情况下,新设计芯片的分离纯度为100%,效率为64%,而实验芯片的分离纯度为100%,效率为25%.最后,我们进一步确定使用新设计芯片时的血液注入的流量以及稀释程度.结果显示,随着进口处红细胞压积的增大,分离效率先增大后减小,在红细胞压积为10.4%时达到最大值.此外,当进口流量小于13.3μl/h,分离纯度可达到100%;而当进口流量大于该值时,纯度将降低.因此,新芯片进口处的最佳注入流量为13.3μl/h,最佳红细胞压积为10.4%,此时分离效率为64%,分离纯度为100%.这体现了数值模拟在生物医学设备设计方面的辅助作用.红细胞压积是人体健康状况的重要指标之一,对疾病的诊治具有重要意义.我们针对三种不同的微血管网络,研究不同时刻和不同位置的红细胞压积分布,展示细胞压积的时空不均匀性.从时间上,细胞压积有的变化比较平稳,而有的一直处于周期振荡.从空间上,不同直径血管的细胞压积分布不均匀;而且不同类型血管,如直管、弯管以及分叉管,也体现出不同的细胞压积分布.研究结果显示,这种时空不均匀性主要归因于网络中存在大量的分叉管道,如果血管为直管道,那么细胞压积将在时空上比较平稳.基于此,我们详细研究了分叉血管中血液流量、管道直径、管道曲率对细胞压积的影响.相比于当前流行的细胞压积测量手段,如抽血化验(只能给出总体平均测量数据)和核磁共振成像(只能得到二维测量数据),我们的研究结果(可以提供不同时刻不同位置的三维数据)对相关疾病的诊治提供了更加准确的信息.
其他文献
湖南省湘西自治州降水丰富,其地貌类型多样,以中低山为主,地质条件复杂,是地质灾害易发区。根据《湘西自治州地质灾害防治规划(2011—2020)》等资料统计,湘西自治州总共包括708处地质灾害隐患点,其中主要包含滑坡、泥石流、崩塌等,其中实有滑坡隐患点407处,占地质灾害隐患点总数的57.49%。本文利用湘西自治州实有滑坡灾害隐患点数据作为参考值,并结合该区域的高程、坡度等地形地貌数据,断层、地层年
学位
学位
学位
随着智能制造工业4.0时代的到来,虚拟数控仿真技术得到了迅速的发展。在智能制造加工领域内,整体叶轮叶片之类的复杂曲面零件由于其材料去除率大,叶片薄、叶展长等特点造成了其切削性能较差,易发生变形、颤振以及共振等现象。在对此类零件的实际加工过程中,由于其加工过程相关切削物理量的未知性,因此极易造成切削材料以及人力的浪费。故针对此类零件的铣削加工的特点,构建虚拟数控铣削仿真系统以实现高效高精加工是当今制
跨境并购作为快速进军海外市场的重要手段得到了很多企业的青睐,尤其是在我国从2015年开始正式实施“一带一路”的建设工作之后,对外投资的环境逐渐地得到了改善,经济要素的流动性也随之逐渐提高,再加上供给侧改革等一系列政策的逐步落实和深入推进,从而引导我国企业积极响应“走出去”的战略,使得跨境并购进入了快速发展时期。但是跨境并购相对于境内并购来说会存在着额外的风险,比如交易双方所在国家或地区的政治关系可
冶金高炉用铸钢冷却壁,具有延伸率高、抗拉强度高、熔点高、抗热冲击性好等特点,应用于高炉生产能显著提高高炉寿命。但铸钢冷却壁铸件的浇注温度高,浇注凝固时预置的冷却水管极易发生变形和熔穿,凝固后冷却壁外表面粘砂也较严重。本文以解决上述问题为目标,研究了铸钢冷却壁的防熔穿工艺、开发了铸钢件防粘砂涂料,对铸钢冷却壁取代铸铁冷却壁在高炉上推广应用具有重要意义。本文首先采用铸造有限元模拟软件,对铸钢冷却壁的充
在能源问题以及环境问题愈加严重的现在,汽车轻量化变得更加迫切。使用铝-钢异种结构同传统钢结构相比重量更轻,且比全铝结构更加安全。而铝钢之间物理性能差异较大,难以相互连接,同时铝钢又极易在焊接过程中发生冶金反应,生成具有高脆硬相的铁铝金属间化合物,导致构件的力学性能的下降,限制其应用。激光熔钎焊具有热输入可控性高、焊接效率高、便于实现自动化等优点,在汽车板焊接中应用广泛。同时大量研究表明,磁场对于导
汽车车身选材作为汽车生产中最关键的环节之一,对汽车的安全性、操稳性、耐久性等各方面的性能起着决定性的作用。而由于车身零件种类繁多,结构复杂,性能要求千变万化,一套完整可行的车身选材方案在生产环境中能够起到巨大的推动作用。本文针对现存的多属性决策选材方案在备选材料评估属性、属性赋权方法、整体选材流程等方面的不足,基于专业化的车身材料数据库,提出了综合考虑经验主观性和物理客观性的多属性决策选材方案,同
切削过程中,随着切削参数的连续变化,切屑的流动方向即流屑角会在某一时刻发生突跳性的改变,这就是流屑角突变现象。实验发现,流屑角突变发生瞬间,切削力发生最高可达75%的变化,其对工艺效果的影响和在节能降耗方面的应用潜力不容小觑。通过理论分析和实验研究,探明流屑角突变的内在机理并建立能准确预报流屑角突变的数学模型,是合理利用该突变现象的前提,具有重要的意义。本文首先在分析不等分剪切区模型的参数敏感性的