1. 前言
多孔介质具有较大的比表面积、良好的均匀透过性和耐高温、耐气候性、抗腐蚀等性能,使多孔材料具有其它材料难以取代的优异性质,它可以用于冶金、化工、环保、能源、生物等行业中的气体分离、液体分离、液相色谱柱填充剂、离子选择性电极、放射性废弃物的处理、催化剂或酶的载体等方面 [1]。多孔碳材料种类繁多包括活性炭、介孔碳、活性炭纤维、碳纳米管、多孔碳分子筛等,多孔碳材料具有大量纳米孔道,巨大的比表面积,在吸附、催化、储能、电容等方面潜在巨大应用价值 [2]。基于碳材料的制备方法有两种:活化法和模板法。它的制法简单,与此同时,它也具有比表面积大、导电传热性好、质量轻、化学稳定性好等特点 [3]。而单位体积的多孔介质中最大孔隙团体积所占的比率反映了多孔介质中流体的运动与扩散现象 [4] [5],表示孔隙的连通状态可以决定多孔介质的主要性质。
通过研究流体在多孔介质内的渗流问题进而可以提高多孔介质的性能,国内外学者取得了较多成果。十九世纪五六十年代,Darcy研究城市地下水问题得出达西定律,渗流力学从此蓬勃发展;Pascal研究发现当流体处于低速流动时压力梯度克服流体主要的力是粘性阻力 [6];Brinkman将Darcy流与Stokes流结合起来,对达西定律引入有限黏度并用于高孔隙率多孔介质;十九世纪八十年代初,Vafai和Tien假设流体连续流动,获得了应用较广的广义达西定律。这些在渗流模拟方面目前常用的方法有:格子法 [7]、有限元法 [8] 与逾渗法 [9] 等。
诸多研究表明 [10] [11] [12],多孔介质的渗透效果由孔隙结构参数共同决定,并不是某一孔隙结构参数的单值函数,也不随某一结构参数发生简单的递增或递减变化,定量描述这种多值函数和对多孔介质孔隙精确取值是当前需要解决的一大难题,随着计算机断层扫描(CT)技术分辨率和精度的提高,CT技术逐渐应用于分析散体孔隙特征 [13] [14] [15] [16] [17]。在国内,利用CT扫描方法获取多孔介质孔隙参数工作开展较晚,Lin等 [18] 学者借助CT技术揭示了矿岩散体在浸出过程中孔隙率随浸柱高度的空间变化规律。Miller等人 [19] 通过µ-CT扫描微型浸柱,分析了矿物揭露程度与矿石颗粒径的关系。杨保华 [20] 通过医学CT扫描浸柱,分析了二维孔隙率与浸柱断层的分形关系,同时也研究了矿石间孔隙的空间分布特征。CT技术在模拟的研究取得了重大的成果 [21] - [27],具有丰富的理论指导。它和三维重构模型方法相结合,建立微观孔隙结构模型,为其渗流模拟奠定基础。
本文以煤基多孔碳为研究对象,利用CT扫描技术获得多孔碳二维切片照片,通过三维图像处理技术处理二位照片重构得到多孔碳微观模型结构,再用FLUENT进行模拟。探究多孔碳渗流不同介质(H2O/CO2)的模拟过程和影响因素,得出有利于多孔碳渗流的规律,最大化提高多孔碳的适用范围和性能。为后续多孔介质的研究提供基础。
2. 模型与数值分析
2.1. 模型构建
在通过CT成像技术得到二维切片照片时由于操作环境和技术缺陷等的影响,使得在CT成像过程中存在噪音,得到的二维切片照片总体较暗、前后存在差异以及等问题,因此就需要对图像进行处理,包括:CT图片的选取、亮度调节、滤波处理等步骤。
2.1.1. CT图片的选取和调节
本文选取尺寸为直径8 mm长度为20 mm的圆柱试样进行CT扫描,最终获得了1000张像素为986 × 1014体素为2.033 μm的16位二维灰度切片图片,体素即就是所能识别的最小尺寸,这些CT图片总体比较暗,如图1(a)展示了CT照片中第400张图处理前后照片。考虑计算机计算能力等的问题,图片过大需要千万级的网格数量,求解困难。在不改变灰度值的条件下,对这1000张照片进行亮度调节和选择,如图1(b)。最终选取了300~699张像素为200 × 200的连续照片。
(a) 第400张未处理前CT图 (b) 第400张处理后的照片
Figure 1. The 400th CT image before and after processing
图1. 第400张处理前后CT图
2.1.2. 三维重构
三维重构时多孔介质有限元模型构建至关重要的一步,它包括图片导入、体渲染、阈值分割、孔隙提取、骨架化等步骤。把经过处理的CT图片在三维可视化软件AVIZO中,设置X、Y、Z三方向的体素为2.033 μm,由于设备的局限性造成了图片有锯齿状的条纹产生。如图2(a)所示,亮的为基质暗的孔隙,略有模糊,为了进一步看清基质和孔隙的分布,随之对图2(a)所示多孔碳进行体渲染如图2(b)所示,图中黄色为基质绿色的孔隙。
(a) 多孔碳渲染前YZ面
(b) 多孔碳渲染后YZ 面
Figure 2. Three-dimensional reconstruction of porous carbon before rendering and after rendering
图2. 渲染前后的多孔碳三维重构
2.1.3. 网格划分
![](//html.hanspub.org/file/1-2190081x13_hanspub.png)
Figure 3. Three-dimensional grid division of porous carbon
图3. 多孔碳三维孔隙网格划分
利用ICEM对三维重构的模型进行网格修复,采用非结构化网格,手动去除重复面、补孔、去除相交面等,修补好的网格如图3所示。在经过一系列网格处理后,总的网格数量为45万个,网格质量均大于零,可以进行模拟。
2.2. 数值分析
2.2.1. 数值计算
本文选取多孔碳建立有限元模型,假定流体为连续流体,其流动规律可用无滑移边界条件的N-S方程进行阐述。流体遵循质量、能量以及动量守恒。
1) 质量守恒方程
质量守恒方程是单位时间内流体微元的质量变化率为0,其关系式如式1所示:
(1)
式中,ρ为流体的密度;ux、vy、wZ表示流体速度分量。
2) 动量守恒方程
动量守恒方程(N-S)流体微元动量对时间的变化率与外界作用力矢量和。如式2示:
(2)
式中,ρ为流体微元受到的压力;τij为由于分子之间的粘性力分量;Fi为作用在流体微元上的质量力。
3) 能量守恒方程
能量守恒方程也称为能量守恒定律,它的基本表述关系如式3:
(3)
式中,λ代表传热系数;Jj代表组分j的扩散流量;Sh指代的是流体粘性热耗散项;E代表流体的内能、动能、势能的总和。E表达关系如式4:
(4)
2.2.2. 边界条件设置
压力进口、压力出口并采用不同流体介质进行渗流模拟,模拟不同介质CO2和H2O来对多孔碳渗流进行分析,进口压力范围为1.01 × 102~1.01 × 107 Pa (小于多孔碳的最大抗压强度12.17 MPa),出口压力为0 Pa,再选取不同进口流速分别模拟。这里假设介质和多孔碳模型孔隙都是稳态情况下进行。
3. 不同流体结果分析
本文研究了以H2O、CO2为流体介质,在不同压力作用下通过多孔材料时的渗流情况。图4分别为以H2O(a)和CO2(b)为流体介质,进口压力101 kPa下多孔介质整体压力分布云图。图5分别为以H2O(a)和CO2(b)为流体介质,进口压力101 kPa下多孔介质整体速度分布云图。
图4可得出:H2O、CO2的压力云图在整体规律上呈现出一致性,但在出口面以及压力最小值方面存在差异,以水为流体介质时,红色部分(高压)较少,浅绿色(低压)较多。图5可看出:H2O、CO2的速度云图在整体规律上也呈现出相似性,但是在局部存在差别,其中以水为流体介质时,蓝色部分(低速)最多,黄色部分(高速)最少,以二氧化碳为流体介质时较少,这反映了其在多孔碳的渗流速度。
![](//html.hanspub.org/file/1-2190081x19_hanspub.png)
Figure 4. Overall pressure distribution of porous carbon at 101 kPa
图4. 101 kPa下多孔碳整体压力分布
![](//html.hanspub.org/file/1-2190081x20_hanspub.png)
Figure 5. Overall velocity distribution of porous carbon at 101 kPa
图5. 101 kPa下多孔碳整体速度分布
![](//html.hanspub.org/file/1-2190081x21_hanspub.png)
Figure 6. Velocity distribution (H2O) of each section of porous carbon at 101 kPa
图6. 101 kPa下多孔碳各截面速度分布(H2O)
![](//html.hanspub.org/file/1-2190081x22_hanspub.png)
Figure 7. Velocity distribution (CO2) of each section of porous carbon at 101 kPa
图7. 101 kPa下多孔碳各截面速度分布(CO2)
![](//html.hanspub.org/file/1-2190081x23_hanspub.png)
Figure 8. XOZ and YOZ surface velocity distribution (H2O) at 101 kPa
图8. 101 kPa下XOZ与YOZ面速度分布(H2O)
![](//html.hanspub.org/file/1-2190081x24_hanspub.png)
Figure 9. XOZ and YOZ surface velocity distribution (CO2) at 101 kPa
图9. 101 kPa下XOZ与YOZ面速度分布(CO2)
为了研究多孔碳内部在不同渗流介质下的渗流规律,依次模拟了在101 kPa下h = 100、80、60、40、20、0 mm处各XOY面以及XOZ面与YOZ面速度分布。图6与图7分别为以H2O和CO2为流体介质在101 kPa下截面高度为100、80、60、40、20、0 mm处的各截面速度分布云图,图8与图9分别为以H2O和CO2为流体介质在101 kPa下在XOZ面、YOZ面以及XOZ与YOZ正交面的渗流速度分布云图。
分析图6~9的各截面速度云图可看出:在这两种介质下,平均速度整体呈现从进口先增大再平缓变化再减小再增大的规律,流体局部出现回流,从图6与图7可看出,进口面(a)速度变化范围最小,出口面(f)速度变化范围最大,截面在80(b)、60(c)、40(d) mm的3个截面处速度整体大于进口,比较平稳,这说明该多孔碳中部位置孔隙较多分布相对规律。这说明该多孔材料边缘孔隙较多,截面中间位置孔隙较少,根据速度云图整体来看具有平稳的渗流效果。
进一步对H2O和CO2在不同压力下进行模拟的结果如图10和图11。两幅图展示如下的特点:两种流体介质的在不同压力下,进口压力与出口流速呈现非线性递增的关系,随着进口压力的增加,出口流速先快速再缓慢增加,然后逐渐趋于平缓。其中,液态水的增长较为平缓,气体的介质二氧化碳的增长曲线较陡峭。而进口压力与稳态时间呈现反比例函数关系,液态水所需的稳定时间较长,下降比较缓慢,气态二氧化碳所需的稳定时间较短,曲线急速下降后,稳定时间不随进口压力而变化。此外,H2O和CO2两种介质在该多孔介质内的存在启动压力与稳定压力,并且水的启动压力在3 MPa左右,二氧化碳的启动压力在0.5 MPa左右,水的稳定压力在9 MPa左右,二氧化碳的稳定压力在6 MPa左右,从而得出液
![](//html.hanspub.org/file/1-2190081x25_hanspub.png)
Figure 10. Relationship between inlet pressure and outlet flow rate of H2O and CO2
图10. H2O和CO2进口压力与出口流速关系图 [28]
<img src=//html.hanspub.org/file/1-2190081x26_hanspub.png width=992.125015258789 height=462.374992370605 class=100 />
图11. H2O和CO2进口压力与稳定时间关系图 [28]
态的介质的启动压力与稳定压力均大于气态介质的稳定压力与启动压力,并且气态介质之间二者的数值也与气态介质的性质有关。
4. 结论
1) 多孔碳的渗流过程是个复杂的过程,并非某一特定结构参数所决定,不同的渗流介质会有不同的渗流效果,多孔碳对某些渗流介质展现出了特定的吸附性能。
2) 多孔碳渗流不同介质总体上呈现相似规律,局部存在差异,相同条件下,气体在多孔介质内的流速远大于液体,并且不同的流体介质在多孔介质中渗流时存在不一样的启动压力与稳定压力,水、二氧化碳启动压力3 MPa和0.5 MPa;稳定压力9 MPa和6 MPa。
基金项目
徐州市重点研发计划(KC19222),中央高校基本科研业务费专项资金资助(2019GF02)。