Numerical Simulation of Seepage Flow in Porous Bioelastic Materials

doi:10.12422/j.issn.1006-396X.2024.02.010

Abstract

Abstract:

The Biological seepage studies the seepage of biofluids in living organisms and fluids containing microorganisms in non-biological porous media.The mass transfer diffusion osmosis phenomenon of a porous bioelastomer material,poly(glycidyl sebacate) (PGS),implanted into human soft tissues was simulated.The pore structure of the PGS material was characterised by N₂ adsorption-desorption and the scaffold model with different pore numbers and pore diameters was designed using the multi-physics field simulation software-COMSOL to investigate the effects of pore and pore diameter parameters on blood osmosis when the pore size was constant.COMSOL was used to design the adapted PGS stent models,numerically simulate the characteristics of the blood flow when blood flowed through the PGS stent,and analyse the kinetic viscosity and the shear rate of the blood field with the theory of fluid dynamics.Comparative analyses of seepage pressure and diffusion at different blood inlet velocities were carried out. The results show that the PGS material is a material that tends to be mesoporous.When the blood flows inside the porous scaffold, the kinetic viscosity varies with the shear rate,indicating that the blood seepage inside the scaffold is a kind of non-Newtonian fluid seepage.The diffusion speed of the blood inside the porous scaffold is different under different inlet velocities and the larger the inlet velocity is,the higher the pressure,and the faster the mass transfer diffusion speed will be.

Key words: Bio?seepage, Porous bioelastomers, Mass transfer diffusion, Non-Newtonian fluids, COMSOL numerical simulation

摘要：

模拟了一种多孔生物弹性体材料—聚癸二酸甘油酯（PGS)被植入人体软组织的传质扩散渗流现象；用N₂吸附?脱附法表征了PGS材料的孔隙结构；用多物理场仿真软件COMSOL设计不同孔道数目和孔径的支架模型，探究了孔隙度不变时孔道和孔径参数对血液渗流的影响；用COMSOL设计适配的PGS支架模型，结合流体动力学理论，对血液流过PGS支架时的血流特性进行数值模拟，对血流场的动力黏度和剪切速率进行了分析；对不同血液入口速度下的渗流压力和扩散情况进行了对比分析。结果表明，PGS为一种趋于介孔的材料；当血液在多孔支架内流动时，动力黏度随剪切速率的变化而变化，说明血液在支架内的渗流是一种非牛顿流体渗流；在不同入口速度下，血液在多孔支架内的扩散速度不同，入口速度越大，压力越大，传质扩散速度越快。

关键词: 生物渗流, 多孔生物弹性体, 传质扩散, 非牛顿流体, COMSOL数值模拟

CLC Number:

TQ021

Miaochao CHEN, Sheng YANG, Kaixuan GUO, Jinbao FENG, Jiao YU. Numerical Simulation of Seepage Flow in Porous Bioelastic Materials[J]. Journal of Petrochemical Universities, 2024, 37(2): 73-80.

陈妙超, 杨盛, 郭凯旋, 冯金宝, 玉姣. 多孔生物弹性材料的渗流数值模拟[J]. 石油化工高等学校学报, 2024, 37(2): 73-80.

Figures/Tables 12

Fig.1 Physical image of PGS porous bioelastomer material

Fig.2 Pore size distribution and cumulative pore volume map of PGS material

Fig.3 Cuboid pore canal model

Fig.4 Seepage velocity of axial section in the model of pore canal

Fig.5 MATLAB seepage velocity fitting curves

Table 1 MATLAB solution result

孔径/cm	a	b	c	d×10⁴	Q/（m³·s^-1）
0.2	52.708 3	-25.737 2	3.454 8	-1.754 9	0.003 8
0.4	7.955 1	-6.881 3	1.759 4	-1. 569 9	0.015 2

Fig.6 Simulation result of dynamic viscosity and velocity diagram of PGS stent

Fig.7 Dynamic viscosity and shear rate curves of blood on vertical and horizontal lines

Fig.8 Simulation model

Table 2 PGS scaffold and blood basic parameters

参数	数值
支架密度/（kg·m^-3）	$1 130$
支架孔隙率	$0.5$
渗透率/m²	$2.38 × 10 - 13$
血液密度/（kg·m^-3）	$1 060$
血液动力黏度/（Pa·s）	$0.003 5$ ^[19]
血液浓度/(mol·m^-3)	$4.0$
扩散系数/（m²·s^-1）	$1.0 × 10 - 9$

Table 2 PGS scaffold and blood basic parameters

参数	数值
支架密度/（kg·m^-3）	$1 130$
支架孔隙率	$0.5$
渗透率/m²	$2.38 × 10 - 13$
血液密度/（kg·m^-3）	$1 060$
血液动力黏度/（Pa·s）	$0.003 5$ ^[19]
血液浓度/(mol·m^-3)	$4.0$
扩散系数/（m²·s^-1）	$1.0 × 10 - 9$

Fig.9 The process of seepage pressure distribution and blood concentration diffusion when inlet velocity is 0.1 m/s

Fig.10 The process of seepage pressure distribution and blood concentration diffusion when inlet velocity is 0.2 m/s

References 27

1	邓英尔.高等渗流理论与方法［M］.北京：科学出版社，2004.
2	周鑫，范方政.孔隙介质中渗流运动与扩散过程的关系研究［C］//北京力学会第二十五届学术年会会议论文集.北京：北京力学会，2019：6.
3	郭尚平，于大森，吴万娣.脏器渗流多孔介质的物理特征［J］.力学学报，1982，1（1）：26⁃32.
	GUO S P，YU D S，WU W D.Physical characteristics of porous media in organ seepage［J］.Chinese Journal of Mechanics，1982，1（1）：26⁃32.
4	DASKALAKIS E，LIU F，COOPER G，et al.Bioglasses for bone tissue engineering：Bio⁃materials and prototyping applications in medicine［M］.2nd ed.Switzerland：Springer International Publishing AG，2021.
5	李培超，孔祥言，卢德唐.饱和多孔介质流固耦合渗流的数学模型［J］.水动力学研究与进展，2003，18（4）：419⁃426.
	LI P C，KONG X Y，LU D T.Mathematical modeling of flow in saturated porous media on account of fluid⁃structure coupling effect［J］.Journal of Hydrodynamics，2003，18（4）：419⁃426.
6	赵静一，李侃，吴厚正，等.组织工程支架内流场的数值模拟［J］.液压与气动，2006（7）：19⁃21.
	ZHAO J Y，LI K，WU H Z，et al.Numerical simulation of flow field for tissue engineering scaffold［J］.Chinese Hydraulics & Pneumatics，2006（7）：19⁃21.
7	ALI D，OZALP M，BLANQUER S B G，et al.Permeability and fluid flow⁃induced wall shear stress in bone scaffolds with TPMS and lattice architectures：A CFD analysis［J］.European Journal of Mechanics⁃B/Fluids，2020，79：376⁃385.
8	OMAR A M，HASSAN M H，DASKALAKIS E，et al.Geometry⁃based computational fluid dynamic model for predicting the biological behavior of bone tissue engineering scaffolds［J］.Journal of Functional Biomaterials，2022，13（3）：104.
9	LOH X J，KARIMA A A，OWHA C.Poly（glycerol sebacate） biomaterial：Synthesis and biomedical applications［J］.Journal of Materials Chemistry B，2015，3（39）：7641⁃7652.
10	CHHABRA R P.Rheology of complex fluids［M］.Springer：Springer Science Business Media，LLC，2010.
11	SHIBESHI S，COLLINS W E.The rheology of blood flow in a branched arterial system［J］.Applied Rheology，2005，15（6）：398⁃405.
12	BERNSDORF J，WANG D.Non⁃Newtonian blood flow simulation in cerebral aneurysms［J］.Computers & Mathematics with Applications，2009，58（5）：1024⁃1029.
13	刘海龙，吴淑红.低渗透油藏渗流压力差分方法数值模拟［J］.辽宁石油化工大学学报，2016，36（4）：38⁃42.
	LIU H L，WU S H.Numerical simulation of fluid flow in low permeability reservoir using finite difference method［J］.Journal of Liaoning Shihua University，2016，36（4）：38⁃42.
14	张文娟，王媛，倪小东.Forchheimer型非达西渗流参数特征分析［J］.水电能源科学，2014，32（1）：52⁃54.
	ZHANG W J，WANG Y，NI X D.Analysis of Forchheimer type non⁃Darcy seepage parameters［J］.Hydropower Energy Science，2014，32（1）：52⁃54.
15	贾双珠，吴林媛，陈炷霖，等.微介孔材料的孔结构分析表征［J］.分析测试技术与仪器，2019，25（3）：141⁃147.
	JIA S Z，WU L Y，CHEN Z L，et al.Analysis and characterization of pore structure of micro⁃mesoporous materials［J］.Analysis and Testing Technology and Instruments，2019，25（3）：141⁃147.
16	DUCK F A.Physical properties of tissue：A comprehensive reference book［M］.London：Academic Press，1990：138.
17	GIJSEN F J，ALLANIC E，VAN DE VOSSE F N，et al.The influence of the non⁃Newtonian properties of blood on the flow in large arteries：Unsteady flow in a 90 degrees curved tube［J］.Journal of Biomechanics，1999，32（7）：705⁃713.
18	POMERANTSEVA I，KREBS N，HART A，et al.Degradation behavior of poly（glycerol sebacate）［J］.Journal of Biomedical Materials Research.Part A，2009，91A（4）：1038⁃1047.
19	刘小月，木合塔尔·克力木，巨刚，等.支架置入前脑动脉瘤血流场的特性分析［J］.中国组织工程研究，2016，20（38）：5722⁃5729.
	LIU X Y，KELIMU M H T E，JU G，et al.Blood flow field characteristics of cerebral aneurysm before stent implantation［J］.Chinese Journal of Tissue Engineering Research，2016，20（38）：5722⁃5729.
20	刘鑫，齐磊，左明星，等.天然生物材料在组织工程和再生医学中的应用［J］.医疗卫生装备，2019，40（3）：98⁃103.
	LIU X，QI L，ZUO M X，et al.Application of natural biomaterials in tissue engineering and regenerative medicine［J］.Chinese Medical Equipment Journal，2019，40（3）：98⁃103.
21	LEE S J，YOO J J，ATALA A.Tissue Engineering［M］.Singapore：Springer Nature Singapore Pte Ltd，2018：17⁃51.
22	王东青，魏静，于春梅，等.MCM⁃41/β介孔—微孔复合分子筛的合成、表征及催化性能［J］.东北石油大学学报，2013，37（3）：100⁃104.
	WANG D Q，WEI J，YU C M，et al.Synthesis，characterization and catalytic performance of MCM⁃41/β mesoporous microporous composite molecular sieves［J］.Journal of Northeast Petroleum University，2013，37（3）：100⁃104.
23	封心领，石泉，顾凯，等.介孔碳片的合成及吸附脱硫性能［J］.石油化工高等学校学报，2011，24（6）：14⁃18.
	FENG X L，SHI Q，GU K，et al.Synthesis and adsorption desulfurization properties of mesoporous carbon platelets［J］.Journal of Petrochemical Universities，2011，24（6）：14⁃18.
24	郭阳，李士凤，姚淑华，等.高比表面积介孔α⁃FeOOH的制备及其Sb（Ⅴ）吸附性能［J］.化工环保，2019，39（5）：562⁃567.
	GUO Y，LI S F，YAO S H，et al.Preparation of mesoporous α⁃FeOOH with high specific surface area and its adsorption capability to Sb（Ⅴ）［J］.Chemical Environmental Protection，2019，39（5）：562⁃567.
25	胡智辉，宫秀坤，张东辉.负载氨基的介孔SiO₂用于CO₂/N₂分离［J］.低碳化学与化工，2009，34（5）：1⁃4.
	HU Z H，GONG X K，ZHANG D H.Amino⁃functionalized silica for CO₂ and N₂ separation［J］.Low⁃Carbon Chemistry and Chemical Engineering，2009，34（5）：1⁃4.
26	谢娟，魏雨，赵新强，等.自组装技术制备Au/SiO₂催化剂及其催化空气氧化环己烷的性能研究［J］.石油炼制与化工，2009，40（2）：31⁃36.
	XIE J，WEI Y，ZHAO X Q，et al.Preparation of Au/SiO₂ catalysts by self⁃assembly technique and its catalytic performance of oxidating cyclohexane with air［J］.Petroleum Processing and Petrochemicals，2009，40（2）：31⁃36.
27	孙伟杰，王素艳，谷桂娜，等.介孔分子筛SBA⁃16固定化木瓜蛋白酶性质的研究［J］.当代化工，2010，39（2）：129⁃131.
	SUN W J，WANG S Y，GU G N，et al.Immobilization of papain on mesoporous molecular sieve SBA⁃16［J］.Contemporary Chemical Industry，2010，39（2）：129⁃131.

[1]	Huan WANG, Xiaoyun XIONG, Yunfeng ZHENG, Xiangbo SUN, Huimin GUAN, Qiang LI, Lijuan SONG. Study on the Correlation between Mass Transfer Performance and Pore Structure and Acidity of FCC Catalysts [J]. Journal of Petrochemical Universities, 2024, 37(3): 58-65.
[2]	Mengyao Tang, Tianyi Sun, Jiashu Li, Yan Ding, Jinlong Li. Molecular Simulation of Hydrogen Adsorption Performance in Graphene and Carbon Nanotubes [J]. Journal of Petrochemical Universities, 2022, 35(4): 10-17.