基于多孔介质模型的膜式氧合器内部流场分析

doi:10.16085/j.issn.1000-6613.2019-1068

摘要/Abstract

摘要：

膜式氧合器内部流体运动特性对其性能有重要影响，利用计算流体力学（CFD）对氧合器模型进行流体动力学分析是预测其性能的重要方法之一。本文基于压降实验计算氧合器纤维束的黏性阻力系数，建立了各向同性多孔介质模型。采用RNGk-ε湍流模型对不同流量下氧合器内部流场进行计算，得到了血液速度、压力和壁面剪切应力分布云图。发现随着流量的增加，氧合器内部速度梯度分布形式基本保持不变，压力分布呈倾斜状态且逐渐减小，大部分压力损失位于纤维束内，其中53.3%位于氧合室，42.6%位于变温室。氧合器血液的入口及出口位置为血液损伤的高发区域。采用溶血数值预估模型计算得到了氧合器的标准溶血指数NIH。结果表明：在低流量1.65～3.00L/min下，各向同性多孔介质模型的模拟结果与实验结果基本一致，模拟数值与实验数值的偏差会随着液体流量的增加而变大；流量为1.65～6.00L/min时，标准溶血指数NIH为0.0084～0.0835g/100L，满足人体生理允许的使用范围。

关键词: 膜, 氧合器, 计算流体力学, 多孔介质, 溶血

Abstract:

The fluid motion characteristics of the membrane oxygenator have an important influence on its performance, the computational fluid dynamics (CFD) hydrodynamic analysis of the oxygenator model by computational fluid dynamics is one of the important methods to predict its performance. Based on the pressure drop experiment, the viscous drag coefficient of the oxygenator fiber bundle was calculated and an isotropic porous media model was established. The RNGk-ε turbulence model was used to calculate the internal flow field of the oxygenator under different flow rates, and the cloud velocity map of blood velocity, pressure and wall shear stress was obtained. It is found that with the increase of flow rate, the internal velocity gradient distribution of the oxygenator is basically unchanged, the pressure distribution is inclined and gradually decreases, and most of the pressure loss is located in the fiber bundle, of which 53.3% is located in the oxygenation chamber, and 42.6% is located in the temperature chamber. The inlet and outlet locations of the oxygenator blood are high-risk areas of blood damage. The standard hemolysis index NIH of the oxygenator was calculated using the hemolysis numerical prediction model. The results show that the simulation results of the isotropic porous medium model are basically consistent with the experimental results at 1.65—3.00L/min. The deviation between the simulated value and the experimental value will increase with the increase of the liquid flow rate; at 1.65—6.00L/min, the standard hemolysis index NIH is 0.0084—0.0835g/100L, which meets the physiologically acceptable use range of human body.

Key words: membrane, oxygenator, computational fluid dynamics, porous media, hemolysis

中图分类号:

TQ 021.1

叶非华,廖虎,易国斌. 基于多孔介质模型的膜式氧合器内部流场分析[J]. 化工进展, 2020, 39(3): 898-905.

Feihua YE,Hu LIAO,Guobin YI. Internal flow field analysis of membrane oxygenator based onporous media model[J]. Chemical Industry and Engineering Progress, 2020, 39(3): 898-905.

图/表 18

图1 膜式氧合器工作原理

图2 实验流程

图3 氧合室和变温室压降-流量拟合曲线

图4 膜式氧合器几何模型示意图

图5 氧合器结构尺寸（单位：mm）

图6 网格无关性验证

图7 管道轴面速度云图

图8 全局速度流线图

图9 中间管道轴切面平均速度变化

图10 出口管道轴切面平均速度变化

图11 氧合器径向中间截面的速度云图

图12 出口管道速度矢量图

图13 氧合器压力分布

图14 壁面剪切应力分布云图

图15 湍流强度分布云图

图16 流量1.65～6.00L/min时红细胞停留时间

图17 较大流量时红细胞停留时间分布

图18 标准溶血指标NIH随流量的变化

参考文献 21

1	D'ONOFRIO C,VAN L R,ROLLAND S,et al.Three-dimensional computational model of a blood oxygenator reconstructed from micro-CT scans[J].Medical Engineering & Physics,2017,47(190):1-8.
2	FUNAKUBO A,HIGAMI T,SAKUMA I,et al.Development of a membrane oxygenator for ECMO using a novel fine silicone hollow fiber[J].ASAIO Journal,1996,42(5):837-840.
3	ASAKAWA Y,FUNAKUBO A,FUKUNAGA K,et al.Development of an implantable oxygenator with cross-flow pump[J].ASAIO Journal,2006,51(3):291-295.
4	SCHLANSTEIN P C,HESSELMANN F,JANSEN S V,et al.Particle image velocimetry used to qualitatively validate computational fluid dynamic simulations in an oxygenator: a proof of concept[J].Cardiovascular Engineering and Technology,2015,6(3):340-351.
5	FILL B,GARTNER M,JOHNSON G,et al.Computational fluid flow and mass transfer of a functionally integrated pediatric pump-oxygenator configuration[J].ASAIO Journal,2008,54(2):214-219.
6	ROLAND G,RALF B,JUTTA A,et al.Improving oxygenator performance using computational simulation and flow field-based parameters[J].Artificial Organs,2010,34(11):930-936.
7	ZHANG J,NOLAN T D C,TAO Z,et al.Characterization of membrane blood oxygenation devices using computational fluid dynamics[J].Journal of Membrane Science,2007,288(1):268-279.
8	JIAFENG Z,XIAOBING C,JUN D,et al.Computational study of the blood flow in three types of 3D hollow fiber membrane bundles[J]. J. Biomech. Eng.,2013,135(12):121009.
9	PASCAL V.The role of computational fluid dynamics for artificial organ design[J].Artificial Organs,2002,26(7):569–570.
10	MADHANI S P,D’ALOISO B D,FRANKOWSKI B,et al.Darcy permeability of hollow fiber membrane bundles made from membrana polymethylpentene fibers used in respiratory assist devices[J].ASAIO Journal,2016,62(3):329-331.
11	GAGE K L,GARTNER M J,BURGREEN G W,et al.Predicting membrane oxygenator pressure drop using computational fluid dynamics[J].Artificial Organs,2015,26(7):600-607.
12	BHAVSAR S S,THOMAS S R,ULRICH S J A O.Numerical modeling of anisotropic fiber bundle behavior in oxygenators[J].Artificial Organs,2011,35(11):1095-1102.
13	JONES C C,MCDONOUGH J M,PATRIZIO C,et al.Improved computational fluid dynamic simulations of blood flow in membrane oxygenators from X-ray imaging[J].Annals of Biomedical Engineering,2013,41(10):2088-2098.
14	ZHANG J,NOLAN T D C,TAO Z,et al.Characterization of membrane blood oxygenation devices using computational fluid dynamics[J].Journal of Membrane Science,2007,288(1):268-279.
15	CONSOLO F,FIORE G B,PELOSI A,et al.A numerical performance assessment of a commercial cardiopulmonary by-pass blood heat exchanger[J].Medical Engineering & Physics,2015,37(6):584-592.
16	PELOSI A,SHERIFF J,STEVANELLA M,et al.Computational evaluation of the thrombogenic potential of a hollow-fiber oxygenator with integrated heat exchanger during extracorporeal circulation[J].Biomechanics and Modeling in Mechanobiology,2014,13(2):349-361.
17	GIERSIEPEN M,WURZINGER L J,OPITZ R,et al.Estimation of shear stress-related blood damage in heart valve prostheses:in vitro comparison of 25 aortic valves[J].International Journal of Artificial Organs,1990,13(5):300-306.
18	武悦,朱良凡,罗云.计算流体力学方法分析一例喷射悬浮血泵的液力、悬浮及溶血特性[J].机械工程学报,2018,54(20):52-58.
	WU Yue,ZHU Liangfan,LUO Yun.Computational fluid dynamics analysis of an injection suspension blood pump on the hydraulic, suspension and hemolysis property[J].Journal of Mechanical Engineering,2018,54(20):52-58.
19	WU Z J,TASKIN M E,ZHANG T,et al.Computational model-based design of a wearable artificial pump-lung for cardiopulmonary/respiratory support[J].Artificial Organs,2012,36(4):387-399.
20	André GARON,Marie‐Isabelle FARINAS.Fast three-dimensional numerical hemolysis approximation[J].Artificial Organs,2004,28(11):1016-1025.
21	寿宸,郭勇君,苏磊,等.基于快速溶血预估模型的离心血泵叶轮特性数值分析[J].生物医学工程学杂志,2014,31(6):1260-1264.
	SHOU Chen,GUO Yongjun,SU Lei,et al.Numerical assessment of impeller features of centrifugal blood pump based on fast hemolysis approximation model[J].Journal of Biomedical Engineering,2014,31(6):1260-1264.

[1]	张祚群, 高扬, 白超杰, 薛立新. 二次界面聚合同步反扩散原位生长ZIF-8纳米粒子制备聚酰胺混合基质反渗透（RO）膜[J]. 化工进展, 2023, 42(S1): 364-373.
[2]	李宁, 李金科, 董金善. 乙烯裂解炉多孔介质燃烧器的研究与开发[J]. 化工进展, 2023, 42(S1): 73-83.
[3]	王云飞, 秦蕊, 郑利军, 李焱, 李清平. 旋转填充床CFD模拟研究进展[J]. 化工进展, 2023, 42(S1): 1-9.
[4]	孙继鹏, 韩靖, 唐杨超, 闫汉博, 张杰瑶, 肖苹, 吴峰. 硫黄湿法成型过程数值模拟与操作参数优化[J]. 化工进展, 2023, 42(S1): 189-196.
[5]	徐茂淯, 陶帅, 齐聪, 梁林. 圆盘式环路热管的启动特性及温度波动[J]. 化工进展, 2023, 42(9): 4531-4537.
[6]	罗成, 范晓勇, 朱永红, 田丰, 崔楼伟, 杜崇鹏, 王飞利, 李冬, 郑化安. 中低温煤焦油加氢反应器不同分配器中液体分布的CFD模拟[J]. 化工进展, 2023, 42(9): 4538-4549.
[7]	徐杰, 夏隆博, 罗平, 邹栋, 仲兆祥. 面向膜蒸馏过程的全疏膜制备及其应用进展[J]. 化工进展, 2023, 42(8): 3943-3955.
[8]	潘宜昌, 周荣飞, 邢卫红. 高效分离同碳数烃的先进微孔膜：现状与挑战[J]. 化工进展, 2023, 42(8): 3926-3942.
[9]	李雪佳, 李鹏, 李志霞, 晋墩尚, 郭强, 宋旭锋, 宋芃, 彭跃莲. 亲水和疏水改性膜的抗结垢和润湿能力的对比[J]. 化工进展, 2023, 42(8): 4458-4464.
[10]	王报英, 王皝莹, 闫军营, 汪耀明, 徐铜文. 聚合物包覆膜在金属分离回收中的研究进展[J]. 化工进展, 2023, 42(8): 3990-4004.
[11]	卜治丞, 焦波, 林海花, 孙洪源. 脉动热管计算流体力学模型与研究进展[J]. 化工进展, 2023, 42(8): 4167-4181.
[12]	汪健生, 张辉鹏, 刘雪玲, 傅煜郭, 朱剑啸. 多孔介质结构对储层内流动和换热特性的影响[J]. 化工进展, 2023, 42(8): 4212-4220.
[13]	陆诗建, 刘苗苗, 杨菲, 张俊杰, 陈思铭, 刘玲, 康国俊, 李清方. 改良型CO₂湿壁塔内气液两相流动规律及传质特性[J]. 化工进展, 2023, 42(7): 3457-3467.
[14]	冯江涵, 宋钫. 阴离子交换膜电解池的研究进展[J]. 化工进展, 2023, 42(7): 3501-3509.
[15]	陈香李, 李倩倩, 张甜, 李彪, 李康康. 自愈合油水分离膜的研究进展[J]. 化工进展, 2023, 42(7): 3600-3610.