第l卷第1期
2001年南京师大学报(工程技术版)JOURNAL0FNANJINGNORMALUNIVERSITY(ENGINEERINGV01.1No.1ANDTECHNOLoGY)2001
多孔介质传热传质分形理论初析
施明恒,陈永平
(东南大学动力工程系,南京,210096)
[摘要]对分形理论在多孔介质传热传质过程中的应用进行了初步的分析,求出了基于分形理论的多孔介质
固有渗透率和有效导热系数,建立了多孔介质渗流与导热的分形模型.
[关键词]多孔介质;传热传质;分形理论
[中图分类号]TKl24;[文献标识码]A;[文章编号]1008—1925(2001)06一07
O引言
多孔介质是由固体骨架和流体组成的一类复合介质,它构成了地球生物圈的物质基础.多孑L介质传热传质在自然界和人类生产、生活中广泛存在.它对国民经济的发展、科学技术的进步以及人民生活水平的提高具有重要的影响.土壤中水、肥、污染物的吸收、保持和迁移过程的人工控制,节水农业工程的实施,地下岩层中石油、天然气和地下水资源的开采,地热能的开发利用等,都涉及到多孔介质中能量和物质的传输问题;动植物中的生命过程也是在多孔介质中发生的传热传质和生化反应的复杂热物理过程;与人民生活密切相关的农副产品、食品、建材和纺织品的干燥、建筑物的隔热保温也是典型的多孔介质传热传质过程;现代铸造技术、燃烧技术、冷冻技术、催化反应技术和各类轻工技术的发展,都与多孔介质传热传质过程密切相联.因此,研究多孔介质传热传质过程对于改造自然、造福人类都具有重大的经济和社会意义.
从学科发展的角度看,多孔介质传热传质学已经渗透到许多学科和新技术领域,包括能源、材料、环境科学、化学工程、仿生学、生物技术、医学和农业工程,是形成新的交叉和边缘学科的一个潜在生长点.因此,多孔介质传热传质研究,是一项具有重大学术价值、对学科发展和技术创新具有深远影响的研究课题,已成为国内外工程热物理、地球和环境科学中最活跃的前沿研究领域之一.L1矗。
以期以来,人们对多孔介质中的传热传质过程进行了大量的理论和实验研究,在理论模型和热质迁移机理方面已经发展了能量理论、液体扩散理论、毛细流动理论和蒸发冷凝理论等描述多孔介质中热质迁移过程的单一理论模型之后,Philip,DeVries,Luikov又发展了多孔介质热质迁移的热力学理论和综合理论以及相应的数学描述,对多孔介质传热传质的研究起到了重要的推动作用uJ.
但是,由于多孔介质内部结构十分复杂,一般是由大小颗粒、碎片或小组织聚集而成的结构,没有特征尺度且极不规则,其内部发生的热质传递过程与传统的均匀介质中发生的过程有很大的差异,各类迁移参数随着实际多孔介质内部的几何结构的不规律性而出现容积范围内
收稿日期:2000—1z一20
基金项目:国家重点基础研究发展规划资助项目(G200026303)作者简介:施明恒,1939一,东南大学动力工程系教授,博士生导师,主要从事传热传质、汽液两相流方面的研究
万方数据
施明恒,等:多孔介质传热传质分形理论初析
的不均匀性和不确定性.上述各种现有的多孔介质传热传质理论和模型,都是直接或间接地把新研究的多孔介质看作是一种在大尺度上均匀分布的虚拟连续介质,在研究中采用“容积平均”的基本方法,即采用平均物性和空隙的平均几何分布来进行过程的研究.显然,这种“容积平均”的假设与实际多孔介质内部状态存在着很大的差异,因此现有的理论只能近似地在大尺度范围内描述多孑L介质中的传递过程而无法揭示局部和整体之间的本质联系,所得结果与实际测量有较大的偏差.受此限制,多孔介质传热传质研究始终未能有突破性的发展.为此,迫切需要寻找一种描述多孔介质内部结构和迁移参数的新方法,为多孔介质传热传质研究开辟一条新路.
1多孑L介质剖面的分形描述
分形是1975年由美国学者Mandelbrot首先提出的[3].由于分形能反应自然界存在的大量非线性现象和几何形状的客观规律,因此立刻引起了各国科学家的重视,开始了大量的研究,逐步形成了分形几何理论体系.分形几何学是一门以非规则几何形状为研究对象的几何学,它与传统的欧氏几何最大的差别在于空间维数的值域.欧氏几何认为空间的维数是整数,其描述的图形的边界都是规则的且可以用一定的解析式表示,例如直线、平面、球或立方体等等.但是,自然界中大量物体的形状和结构,如土壤、海岸线、多孔物料等,它们在图形上是完全不规则的,使它们的整体与局部都不能用传统的几何语言来描述,其内部发生的过程则不能用简单的线性近似方法来认识和描述.分形是直接从非线性复杂系统本身人手,从未经简化和抽象的研究对象本身去认识其内在的规律性,这是分形理论和线性处理方法的本质区别.分形几何突破了传统几何的局限,认为分形物体的空间维数可以不是整数.如果用数学表达式表示,则有:
N(艿)~∥(1)
式中Ⅳ是分形物体的空间占有积(线、面或体积),艿是度量尺度,d是分形维数.它可以是整数也可以是非整数.两物体只要满足分形维数相等,那么这两个物体是自相似的.
分形可分为两类.一类称为有规分形,它是按一定的数学法则生成的,具有严格的自相似性;另一类是无规分形,其自相似性并不严格,只是在大范围内统计意义下的自相似性.例如多孔介质,其自相似性只有在一定尺度范围内才能成立,且属于统计意义下的自相似性,我们称它为局域分形.
分形理论为描述物体内部复杂结构和空间分布提供了一种新的行之有效的手段,从而为精确研究复杂结构内部发生的各种物理化学过程开辟了一条新路.下面将给出多孔介质剖面的简单分形描述.
图1给出了两种典型的多孑L介质(土壤)剖面图.由图可见,对于每个局部区域来说,多孔介质内部通道呈现出不轨则性,但是从较大范围来看,其剖面骨架面积分布或孔隙分布又具有相似特征.为了证明这一点,下面采用分形理论来计算上述多孔介质剖面的面积分形维数.在剖面图中,取一微元剖面面积S为最小的度量尺度,同时取剖面图中能够反映该多孔介质宏观传递特征的一最小面积元为分形计算的最大度量尺度A对于在S—A面积范围内的不同度量尺度X,如所测得的剖面固体相或孔隙的面积平均值S满足:
S(X)~∥
万方数据 (2)则该多孔介质剖面结构具有统计意义上的自相似性,即剖面面积分布具有分形特征.式(2)中一7——
南京师大学报(工程技术版)第1卷第1期(2001年)d即为土壤剖面面积分形维数.对应于图1的两种多孔介质剖面结构计算出的面积分形维数给出在图2上,由图可见,结构疏松,孔隙度较大的多孔介质,其剖面孔隙面积分形维数较大.
图1两种典型的多孔介质(土壤)的剖面图
坫
M
坞
挖
n加
10.0012.0014.0016.OO
1舭
A.c卢0.8445B.赤0.9173l眦
图2多孔介质剖面孔隙面积分布分形计翼值国
上面给出的是以剖面面积定义的分形维数.描述多孔介质内部构造和通透性还可以采用其它方法定义的分形维数,例如按照组成多孑L介质的颗粒粒度和孔隙分布定义的分形维数,它在描述多孔聚集体形成规律上有特殊的用途.
此外,在研究作随机行走的粒子统计性质及多孔介质内部的输运规律时,需要引入另一个分形维数——谱维数.谱维数与多孔介质中孔隙的通透性或连通性有着密切的关系.谱维数反的定义为:如果在随机行走中粒子经过的时间间隔£以后(或行走步数,z),能够访问的不同孔隙节点或通过的格子数P能够满足:
Pocft/2(3)
则以就定义为谱维数.所谓随机行走是指每经过一个离散的时间f,粒子就在多孔介质划好的网格上行走一步.因此,研究谱维数时的分形结构是由格子组成的分形结构,在空间和时间上是离散的.
图3给出了一个简化的粒子在多孔介质剖面孔隙中随机行走的模型[8].粒子可随机地向周围任一方向移动,每行走一步,运动步数咒加1,同时判断粒子是否是第一次到达该节点(格
一8一万方数据
施明恒,等:多孔介质传热传质分形理论初析
子).若是第一次到达,则访问节点数户
加1;若非第一次到达,则不加.重复多
次行走可得到多组户,,z的值.取对数后
若能进行直线拟合,则认为该多孔通道
存在着谱维数.由拟合后的直线斜率就
可从式(3)求出谱维数以.
2圈丽五瓜面可谪口而丽面万正,+1)囫Tenllinationofwalk多孔介质内部流体渗流的
流体在多孔介质中的流动是一种在分形模型
随机介质中的运动,称为渗流.Darcy早
在1856年的实验中发现[3],在竖立的匀图3一个简化的粒子在多孔介质剖面孔隙中随机行走模型
质砂柱中的饱和水流的质量流速与水头差成正比.石从而导出了著名的达西定律.之后,无论是多孔介质中的饱和渗流还是非饱和渗流,都以Darcy这个宏观经验定律为基础,展开了广泛的研究并取得了很大的进展.但是正如本文一开始就提出的那样,Darcy定律是建立在均匀多孑L介质的基础上的,即使以后有所发展,也是把非均匀多孔介质近似的以当量均匀多孔介质来进行描述,用容积平均的渗透率来计算渗流量.它无法真实的揭示实际多孔介质内部不均匀性和规则的复杂结构对渗流的影响.
多孔介质的固有渗透率,仅受多孔介质结构特性的影响,是一个仅反映流体在多孔介质中渗流运动本质的特性参数.前人对固有渗透率已进行了许多研究与测定.泰勒的研究[43指出:固有渗透率与构成多孔介质的颗粒粒径、孔隙率这两个主要几何特性有关.作者进一步分析认为:固有渗透率还应与多孔介质内部通道的连通性有关.因为即使两类多孔介质的粒径、孔隙率相同,如果其通道的连通性不同(即存在着闭孔和盲道),它们的渗透率也完全不同.因此,本文认为,均匀多孔介质的粒径、孑L隙率和通道的连通性是影响多孔介质固有渗透率的三个主要因素,而对于实际非均匀的多孔介质,影响其固有渗透率的主要因素变为粒径分布、孔隙率分布和通道的连通性.
利用分形理论,可以获得多孔介质粒径分布,孑L隙率和连通性的分形表达式,从而使固有渗透率和多孔介质复杂的几何结构之间建立起有机的联系.
下面首先推导多孔介质颗粒粒径分布的分形表达式.
对于具有分形结构的多孔介质,设M(R)为直径小于R的颗粒质量,M为多孔介质总质量,则直径小于R的颗粒质量百分数为:
Ⅳ(R)一丝婴一口科』Ⅵ(4)
式中6为颗粒粒径分布的分形维数,口为比例系数.
颗粒粒径分布的均方值为:
,PD2+6
露一I—R2dⅣ一口6≤号与JZ十口o(5)
上式为用最大颗粒粒径表示的颗粒粒径均方值的分形表达式.
可以进一步认为具有分形结构的多孔介质其内部孔隙分布也一定具有分形特征.孔隙分万方数据 一9一
南京师大学报(工程技术版)第1卷第1期(2001年)布的分形表达式类似于前述多孔介质剖面面积分形维数表达式的推导过程.假定面孔隙率等于体孔隙率,则剖面孔隙相的面积分形表达式为:
5(X)一C∥,
由此孔隙率分布的分形表达式为:(6)
£,一1一(竺≠)3门一1一(1一c肼叫)吖2
式中X为面积度量尺度,在实际计算中可取最佳度量尺度‘".1
流体在渗流过程中经过时间f能够到达的区域为P(£),其分形的定义式应为:
P(f)一C(r2)‘一凡t72
则(7),.2一(F/C)一锄‘t/%’
由于r2肛是单位时间流体可以占领的渗透区,它表征了多孔介质孔隙的通透能力.其分形表达式为:
尹∥一B2以r2【1一毒J
式中d。一2d,.(8)
Taylor在不考虑土壤通道的通透性的前提下,导出了各相同性均质土壤的固有渗透率与土壤平均颗粒粒径平方以及孔隙率函数£3/(1一£2)成正比的结果[4].本文将在上述分析的基础上,利用分形理论,计及实际多孔介质粒径的分布和孔隙率分布,进一步考虑多孔通道通透性对固有渗透率的影响,提出固有渗透率的下列分形表达式:K一触2高等
式中A为比例系数,F,C为计算谱维数和剖面面积分维数时给出的常数值.
多孔介质单位面积渗流流量应遵循Darcy宏观经验定律,即:
口一K£墨grad^P=郴/c)2/伊一每,筹号籍爿弹㈤(10)
此式即为本文导出的多孔介质渗流流量表达式,与现有对均质多孔介质或容积平均多孔介质渗流流量计算式的不同点在于上式中的固有渗透率是考虑了多孔介质内部的不均匀结构后基于分形理论推导出来的,在多孔介质质量传递与多孔介质内部结构之间建立了有机联系.3多孔介质内部热量传递的分形模型
多孔介质中热量传递是一个复杂的过程.热量既可以通过复杂的固体骨架进行传导,也可以通过孔隙中的流体流动和流体的导热进行传递.如果固体骨架与流体之间存在温差,则它们之间还有导热和对流.一般说来,多孔介质内部的渗流是一个缓慢的过程(特别是对非饱和介质),作为近似分析,可以假设固体骨架和流体之间有热平衡.此外通过流体的热量迁移也可以假定是通过多孔介质孔隙中的气液两相混合物的当量传导进行的.通过上述简化后,多孔介质中的热量传递过程可看作是通过固体骨架和孔隙流体的一组并联热传导过程.
前人曾对多孔介质中的导热过程进行过许多研究,但是由于无法正确的描述内部多孔复杂结构,因此只对各相均匀排列,由相同大小颗粒组成的多孔介质进行过理论推算;或者对不规则结构的多孔介质按容积平均的孔隙率计算平均固体骨架和流体份额,然后进行导热计算,一】O一万方数据
施明恒,等:多孔介质传热传质分形理论初析
造成计算结果误差较大.分形理论的引入,使我们有可能对实际多孔复杂结果进行考虑,即把容积平均的当量导热系数改造成与内部形状直接联系起来的分形导热系数,从而较真实地反映多孔介质内部的导热过程.
为了推导多孔介质的分形导热系数,首先来计算多孔介质固相占有率的分形表达式.
对于一定孔隙率的多孔分形介质,其截面固相面积满足S(x)一似。则单位体积中固相占有份额为:
u一(导)3/2一c3/2x3“_1”2A(11)
作为近似,多孔介质内的导热通道被简化成通过气
隙(气液混合相)和固相的串并联通道,其热阻网络
如图4所示.通过多孔介质的总热阻为:
生!!二生::!±垄生::
“一k九(1一沪73+u)+蟹(p173一c,)P一
圈4多孔介质内的导热通道热阻示意图
最后得多孔介质的有效导热系数为:
~一i一——————永rj夕顷瓦丽F口丽伊砺而丽厂———一】1k九(1一c172x‘4—1’72+c372x3‘。一1’72)+砖(c1/2x‘4—1’/2一C372x3‘扣1’72)u纠,1o、式中砧,凡分别为多孔介质气液混合相和固相导热系数.
式(12)是根据分形理论导出的考虑多孔介质内部几何结构特征的有效导热系数计算公式,由它计算出的导热系数值比通常的容积平均计算公式更接近实际情况.
一般地说,多孔介质中热量的传递在不考虑辐射的状态下,对于未饱和多孔介质和饱和多孔介质来说有较大的差别,前者主要是导热为主.当内部液态介质有相变时(如干燥和冻结),则导热量主要消耗在介质内部的相变上.后者除了导热以外,还有流体在多孔介质孔隙通道中的渗流对流换热.通过单位多孔介质表面的热流密度可分别表示为
g。一九券(纯导热)
gz一九券+咿(导热+内部相变)
口。一九凳+^,(t,一o)(导热+渗流换热)
式中^,为分形渗流换热系数,将在另外的文章中进行分析和研究.
4结束语
本文对分形理论在多孔介质传热传质过程中的应用进行了初步的分析,求出了基于分形理论的多孔介质固有渗透率和有效导热系数.在此基础上,建立了多孔介质渗流与导热的分形模型.虽然本文的模型还是一种建立在分形参数上的当量传热模型,但本文的工作为进一步深人研究复杂的传热传质过程开辟了一条新路.
[参考文献]
[1]林瑞泰.多孔介质传热传质引论[M].北京:科学出版社,1995
[2]施明恒.多孔介质传热传质的进展和展望[M].中国科学基金,1995
[3]张济中.分形[M].北京:清华大学出版社,1995
[4]TaylorDW.FundamentalS0fSoilMechanics[M].JohnWiley&sons,Inc,1948万 方数据—11—
南京师大学报(工程技术版)第l卷第l期(2001年)[5]黑龙江省土壤普查办公室.黑龙江土壤[M].北京:农业出版社,1992
[6]江苏省土壤普查办公室.江苏土壤[M].北京:农业出版社,1995
[7]陈永平,施明恒.基于分形理论的多孔介质导热系数的研究[J].工程热物理学报,1999,20(5):605~612[8]陈永平,施明恒.基于分形理论的多子L介质渗透率的研究[A].中国工程热物理学会传热传质学术会议论
文集[C].苏州:1999
TheFractalTheory0fHeatand
MassTransferinPorousMaterial
ShiMingheng,ChenYOngpjn
(PowerEngineeringDepartmentofSoutheastUniversity,Nanjing,210096,PRC)
Abstract:Thispapermakesapreliminaryanalysisaboutthefractaltheoryusedinheatandmasstransferin
onporousmaterial;thecalculatingequationsofinfiltrationrateandeffectiveheatconductivitycoefficientbased
thetheoryareobtainedandthefractalmodelofinfiltrationandheatconductioninporousmaterialisestablished.Keywords:porousmaterial;heatandmasstransfer;fractaltheory
[责任编辑:刘健]
(上接第5页)
[参考文献]
[1]王树青,元英进.生化过程自动化技术[M].北京:化学工业出版社,1999
[2]王树青.生化反应过程模型化及计算机控制[M].杭州:浙江大学出版社,1998
[3]王树青.发酵过程控制的新进展[J].化工自动化及仪表,1995,22(4):3~10
ComputerIntegratedControlofAntibioticProductionProcess
WangShuq.rTg,ZhangQuan|.ng,ChenQ
(StateKeyLaboratoryofIndustrialCon”olTechnology,ZhejiangUniversity,Hangzhou。310027。PRC)
Abstract:Theresearchonintegratedcontrolofantibioticproductionprocessinvolvesmanipulatingtheinforma—tionaboutproductionprocess,marketinformation,managementinformationandexperimentaldata.Thekeysub—jectsofcomputerintegratedcontrol
toarediscussed,includingdatareconciliation,datamining,optimalcontrolandprocess.optimalmanagementimproveantibioticproductionandfermentation
Keywords:antibioticprocess;computerintegratedcontrol;integratedmanagementcontrol
[责任编辑:刘健]
一12一万方数据
多孔介质传热传质分形理论初析
作者:
作者单位:
刊名:
英文刊名:
年,卷(期):
被引用次数:施明恒, 陈永平东南大学动力工程系,南京,210096南京师大学报(工程技术版)JOURNAL OF NANJING NORMAL UNIVERSITY(ENGINEERING AND TECHNOLOGY)2001,1(1)18次
参考文献(8条)
1. 林瑞泰 多孔介质传热传质引论 1995
2. 施明恒 多孔介质传热传质的进展和展望中国科学基金 1995
3. 张济中 分形 1995
4. Taylor D W Fundamentals of Soil Mechanics 1948
5. 黑龙江省土壤普查办公室 黑龙江土壤 1992
6. 江苏省土壤普查办公室 江苏土壤 1995
7. 陈永平;施明恒 基于分形理论的多孔介质导热系数的研究 1999(05)
8. 陈永平;施明恒 基于分形理论的多孔介质渗透率的研究[会议论文] 1999
本文读者也读过(10条)
1. 陶涛. 陈剑波. 彭晓峰. 王峥. 颜岩 对复杂可变形多孔介质传热传质研究的思考[会议论文]-2001
2. 陈宝明. 王补宣. 张立强 多孔介质传热传质中耦合扩散效应的研究[会议论文]-2003
3. 李小川. 施明恒 多孔介质热导率的数值研究[会议论文]-2007
4. 李小川. 施明恒. LI Xiao-Chuan. SHI Ming-Heng 多孔介质热导率的数值计算[期刊论文]-工程热物理学报2008,29(2)
5. 陈永平. 施明恒 基于分形理论的多孔介质导热系数研究[会议论文]-1998
6. 金仁喜. 淮秀兰. JIN Ren-xi. HUAI Xiu-Lan 多孔介质高强度传热传质的理论研究[期刊论文]-工程热物理学报2005,26(z1)
7. 宫克勤. 孙苗苗 多孔介质中传热传质机理研究[期刊论文]-油气田地面工程2009,28(4)
8. 施明恒. 樊荟 多孔介质导热的分形模型[期刊论文]-热科学与技术2002,1(1)
9. 陈永平. 施明恒. CHEN Yongping. SHI Mingheng 基于分形理论的多孔介质渗透率的研究[期刊论文]-清华大学学报(自然科学版)2000,40(12)
10. 申盛. 刘伟. 陶文铨 多孔介质中自然对流传热的场协同分析[期刊论文]-自然科学进展2003,13(4)
引证文献(18条)
1. 朱群敏. 潘阳. 柏俊 竖直管蓄冰多孔介质局部自然对流的实验研究[期刊论文]-制冷与空调(四川) 2013(5)
2. 司新国. 张玉柱. 张庆军. 莫文玲. 刘丽妹. 马雪刚 低硅烧结矿质量分形维数影响因素[期刊论文]-河北理工大学学报(自然科学版) 2007(3)
3. 张赛. 陈君若. 刘显茜 基于分形和随机扩散描述食品材料的干燥速率[期刊论文]-材料导报 2013(10)
4. 韩博. 张凌. 王童. 张学东. 衣雪峰. 刘向飞 不同剪切作用下水合物钻井液微观结构试验研究[期刊论文]-地质与勘探 2013(4)
5. 韩博. 李承中. 范勇. 姜淼 水合物钻井液微观结构热传导分形模型研究[期刊论文]-地质与勘探 2012(4)
6. 蔡亮. 侯普秀. 虞维平. 张小松 水蒸气在霜层中扩散的分形模型[期刊论文]-中南大学学报(自然科学版)
2010(1)
7. 谢军. 王定标 分形理论在固体吸附式制冷系统中的应用研究[期刊论文]-制冷与空调 2007(6)
8. 李艳松. 张捍民. 孟凡刚. 张兴文. 杨凤林 分形理论在多孔介质渗透性能的研究[期刊论文]-过滤与分离 2006(1)
9. 马德刚. 张书廷. 柴立和 污泥干燥速度曲线的分形维数分析[期刊论文]-天津大学学报 2007(10)
10. 周济福 渗流力学研究的现状和发展趋势[期刊论文]-力学与实践 2007(3)
11. 韩博. 张凌. 蒋国盛. 宁伏龙. 刘力. 范勇. 胡代平. 辜幕丹 不同剪切作用下CMC水溶液的微观结构研究[期刊论文]-地质科技情报 2012(1)
12. 韩博. 张凌. 王童. 张学东. 衣雪峰. 刘向飞 不同剪切作用下水合物钻井液微观结构试验研究[期刊论文]-地质与勘探 2013(4)
13. 杨屹立. 杨玲. 徐武明. 朱留宪. 杨明金 多孔介质物料热风干燥研究进展[期刊论文]-农机化研究 2011(3)
14. 王李斌 多孔介质中的声场和流场及其数值模拟[学位论文]硕士 2006
15. 宋颖韬 多孔介质填充床中流体流动与化学反应过程分析[学位论文]硕士 2005
16. 彭鹏 石墨泡沫的制备及分形研究[学位论文]硕士 2006
17. 刘相东. 杨彬彬 多孔介质干燥理论的回顾与展望[期刊论文]-中国农业大学学报 2005(4)
18. 李艳松 膜污染因子识别与凝胶层结构分析[学位论文]硕士 2006
引用本文格式:施明恒. 陈永平 多孔介质传热传质分形理论初析[期刊论文]-南京师大学报(工程技术版)2001(1)
第l卷第1期
2001年南京师大学报(工程技术版)JOURNAL0FNANJINGNORMALUNIVERSITY(ENGINEERINGV01.1No.1ANDTECHNOLoGY)2001
多孔介质传热传质分形理论初析
施明恒,陈永平
(东南大学动力工程系,南京,210096)
[摘要]对分形理论在多孔介质传热传质过程中的应用进行了初步的分析,求出了基于分形理论的多孔介质
固有渗透率和有效导热系数,建立了多孔介质渗流与导热的分形模型.
[关键词]多孔介质;传热传质;分形理论
[中图分类号]TKl24;[文献标识码]A;[文章编号]1008—1925(2001)06一07
O引言
多孔介质是由固体骨架和流体组成的一类复合介质,它构成了地球生物圈的物质基础.多孑L介质传热传质在自然界和人类生产、生活中广泛存在.它对国民经济的发展、科学技术的进步以及人民生活水平的提高具有重要的影响.土壤中水、肥、污染物的吸收、保持和迁移过程的人工控制,节水农业工程的实施,地下岩层中石油、天然气和地下水资源的开采,地热能的开发利用等,都涉及到多孔介质中能量和物质的传输问题;动植物中的生命过程也是在多孔介质中发生的传热传质和生化反应的复杂热物理过程;与人民生活密切相关的农副产品、食品、建材和纺织品的干燥、建筑物的隔热保温也是典型的多孔介质传热传质过程;现代铸造技术、燃烧技术、冷冻技术、催化反应技术和各类轻工技术的发展,都与多孔介质传热传质过程密切相联.因此,研究多孔介质传热传质过程对于改造自然、造福人类都具有重大的经济和社会意义.
从学科发展的角度看,多孔介质传热传质学已经渗透到许多学科和新技术领域,包括能源、材料、环境科学、化学工程、仿生学、生物技术、医学和农业工程,是形成新的交叉和边缘学科的一个潜在生长点.因此,多孔介质传热传质研究,是一项具有重大学术价值、对学科发展和技术创新具有深远影响的研究课题,已成为国内外工程热物理、地球和环境科学中最活跃的前沿研究领域之一.L1矗。
以期以来,人们对多孔介质中的传热传质过程进行了大量的理论和实验研究,在理论模型和热质迁移机理方面已经发展了能量理论、液体扩散理论、毛细流动理论和蒸发冷凝理论等描述多孔介质中热质迁移过程的单一理论模型之后,Philip,DeVries,Luikov又发展了多孔介质热质迁移的热力学理论和综合理论以及相应的数学描述,对多孔介质传热传质的研究起到了重要的推动作用uJ.
但是,由于多孔介质内部结构十分复杂,一般是由大小颗粒、碎片或小组织聚集而成的结构,没有特征尺度且极不规则,其内部发生的热质传递过程与传统的均匀介质中发生的过程有很大的差异,各类迁移参数随着实际多孔介质内部的几何结构的不规律性而出现容积范围内
收稿日期:2000—1z一20
基金项目:国家重点基础研究发展规划资助项目(G200026303)作者简介:施明恒,1939一,东南大学动力工程系教授,博士生导师,主要从事传热传质、汽液两相流方面的研究
万方数据
施明恒,等:多孔介质传热传质分形理论初析
的不均匀性和不确定性.上述各种现有的多孔介质传热传质理论和模型,都是直接或间接地把新研究的多孔介质看作是一种在大尺度上均匀分布的虚拟连续介质,在研究中采用“容积平均”的基本方法,即采用平均物性和空隙的平均几何分布来进行过程的研究.显然,这种“容积平均”的假设与实际多孔介质内部状态存在着很大的差异,因此现有的理论只能近似地在大尺度范围内描述多孑L介质中的传递过程而无法揭示局部和整体之间的本质联系,所得结果与实际测量有较大的偏差.受此限制,多孔介质传热传质研究始终未能有突破性的发展.为此,迫切需要寻找一种描述多孔介质内部结构和迁移参数的新方法,为多孔介质传热传质研究开辟一条新路.
1多孑L介质剖面的分形描述
分形是1975年由美国学者Mandelbrot首先提出的[3].由于分形能反应自然界存在的大量非线性现象和几何形状的客观规律,因此立刻引起了各国科学家的重视,开始了大量的研究,逐步形成了分形几何理论体系.分形几何学是一门以非规则几何形状为研究对象的几何学,它与传统的欧氏几何最大的差别在于空间维数的值域.欧氏几何认为空间的维数是整数,其描述的图形的边界都是规则的且可以用一定的解析式表示,例如直线、平面、球或立方体等等.但是,自然界中大量物体的形状和结构,如土壤、海岸线、多孔物料等,它们在图形上是完全不规则的,使它们的整体与局部都不能用传统的几何语言来描述,其内部发生的过程则不能用简单的线性近似方法来认识和描述.分形是直接从非线性复杂系统本身人手,从未经简化和抽象的研究对象本身去认识其内在的规律性,这是分形理论和线性处理方法的本质区别.分形几何突破了传统几何的局限,认为分形物体的空间维数可以不是整数.如果用数学表达式表示,则有:
N(艿)~∥(1)
式中Ⅳ是分形物体的空间占有积(线、面或体积),艿是度量尺度,d是分形维数.它可以是整数也可以是非整数.两物体只要满足分形维数相等,那么这两个物体是自相似的.
分形可分为两类.一类称为有规分形,它是按一定的数学法则生成的,具有严格的自相似性;另一类是无规分形,其自相似性并不严格,只是在大范围内统计意义下的自相似性.例如多孔介质,其自相似性只有在一定尺度范围内才能成立,且属于统计意义下的自相似性,我们称它为局域分形.
分形理论为描述物体内部复杂结构和空间分布提供了一种新的行之有效的手段,从而为精确研究复杂结构内部发生的各种物理化学过程开辟了一条新路.下面将给出多孔介质剖面的简单分形描述.
图1给出了两种典型的多孑L介质(土壤)剖面图.由图可见,对于每个局部区域来说,多孔介质内部通道呈现出不轨则性,但是从较大范围来看,其剖面骨架面积分布或孔隙分布又具有相似特征.为了证明这一点,下面采用分形理论来计算上述多孔介质剖面的面积分形维数.在剖面图中,取一微元剖面面积S为最小的度量尺度,同时取剖面图中能够反映该多孔介质宏观传递特征的一最小面积元为分形计算的最大度量尺度A对于在S—A面积范围内的不同度量尺度X,如所测得的剖面固体相或孔隙的面积平均值S满足:
S(X)~∥
万方数据 (2)则该多孔介质剖面结构具有统计意义上的自相似性,即剖面面积分布具有分形特征.式(2)中一7——
南京师大学报(工程技术版)第1卷第1期(2001年)d即为土壤剖面面积分形维数.对应于图1的两种多孔介质剖面结构计算出的面积分形维数给出在图2上,由图可见,结构疏松,孔隙度较大的多孔介质,其剖面孔隙面积分形维数较大.
图1两种典型的多孔介质(土壤)的剖面图
坫
M
坞
挖
n加
10.0012.0014.0016.OO
1舭
A.c卢0.8445B.赤0.9173l眦
图2多孔介质剖面孔隙面积分布分形计翼值国
上面给出的是以剖面面积定义的分形维数.描述多孔介质内部构造和通透性还可以采用其它方法定义的分形维数,例如按照组成多孑L介质的颗粒粒度和孔隙分布定义的分形维数,它在描述多孔聚集体形成规律上有特殊的用途.
此外,在研究作随机行走的粒子统计性质及多孔介质内部的输运规律时,需要引入另一个分形维数——谱维数.谱维数与多孔介质中孔隙的通透性或连通性有着密切的关系.谱维数反的定义为:如果在随机行走中粒子经过的时间间隔£以后(或行走步数,z),能够访问的不同孔隙节点或通过的格子数P能够满足:
Pocft/2(3)
则以就定义为谱维数.所谓随机行走是指每经过一个离散的时间f,粒子就在多孔介质划好的网格上行走一步.因此,研究谱维数时的分形结构是由格子组成的分形结构,在空间和时间上是离散的.
图3给出了一个简化的粒子在多孔介质剖面孔隙中随机行走的模型[8].粒子可随机地向周围任一方向移动,每行走一步,运动步数咒加1,同时判断粒子是否是第一次到达该节点(格
一8一万方数据
施明恒,等:多孔介质传热传质分形理论初析
子).若是第一次到达,则访问节点数户
加1;若非第一次到达,则不加.重复多
次行走可得到多组户,,z的值.取对数后
若能进行直线拟合,则认为该多孔通道
存在着谱维数.由拟合后的直线斜率就
可从式(3)求出谱维数以.
2圈丽五瓜面可谪口而丽面万正,+1)囫Tenllinationofwalk多孔介质内部流体渗流的
流体在多孔介质中的流动是一种在分形模型
随机介质中的运动,称为渗流.Darcy早
在1856年的实验中发现[3],在竖立的匀图3一个简化的粒子在多孔介质剖面孔隙中随机行走模型
质砂柱中的饱和水流的质量流速与水头差成正比.石从而导出了著名的达西定律.之后,无论是多孔介质中的饱和渗流还是非饱和渗流,都以Darcy这个宏观经验定律为基础,展开了广泛的研究并取得了很大的进展.但是正如本文一开始就提出的那样,Darcy定律是建立在均匀多孑L介质的基础上的,即使以后有所发展,也是把非均匀多孔介质近似的以当量均匀多孔介质来进行描述,用容积平均的渗透率来计算渗流量.它无法真实的揭示实际多孔介质内部不均匀性和规则的复杂结构对渗流的影响.
多孔介质的固有渗透率,仅受多孔介质结构特性的影响,是一个仅反映流体在多孔介质中渗流运动本质的特性参数.前人对固有渗透率已进行了许多研究与测定.泰勒的研究[43指出:固有渗透率与构成多孔介质的颗粒粒径、孔隙率这两个主要几何特性有关.作者进一步分析认为:固有渗透率还应与多孔介质内部通道的连通性有关.因为即使两类多孔介质的粒径、孔隙率相同,如果其通道的连通性不同(即存在着闭孔和盲道),它们的渗透率也完全不同.因此,本文认为,均匀多孔介质的粒径、孑L隙率和通道的连通性是影响多孔介质固有渗透率的三个主要因素,而对于实际非均匀的多孔介质,影响其固有渗透率的主要因素变为粒径分布、孔隙率分布和通道的连通性.
利用分形理论,可以获得多孔介质粒径分布,孑L隙率和连通性的分形表达式,从而使固有渗透率和多孔介质复杂的几何结构之间建立起有机的联系.
下面首先推导多孔介质颗粒粒径分布的分形表达式.
对于具有分形结构的多孔介质,设M(R)为直径小于R的颗粒质量,M为多孔介质总质量,则直径小于R的颗粒质量百分数为:
Ⅳ(R)一丝婴一口科』Ⅵ(4)
式中6为颗粒粒径分布的分形维数,口为比例系数.
颗粒粒径分布的均方值为:
,PD2+6
露一I—R2dⅣ一口6≤号与JZ十口o(5)
上式为用最大颗粒粒径表示的颗粒粒径均方值的分形表达式.
可以进一步认为具有分形结构的多孔介质其内部孔隙分布也一定具有分形特征.孔隙分万方数据 一9一
南京师大学报(工程技术版)第1卷第1期(2001年)布的分形表达式类似于前述多孔介质剖面面积分形维数表达式的推导过程.假定面孔隙率等于体孔隙率,则剖面孔隙相的面积分形表达式为:
5(X)一C∥,
由此孔隙率分布的分形表达式为:(6)
£,一1一(竺≠)3门一1一(1一c肼叫)吖2
式中X为面积度量尺度,在实际计算中可取最佳度量尺度‘".1
流体在渗流过程中经过时间f能够到达的区域为P(£),其分形的定义式应为:
P(f)一C(r2)‘一凡t72
则(7),.2一(F/C)一锄‘t/%’
由于r2肛是单位时间流体可以占领的渗透区,它表征了多孔介质孔隙的通透能力.其分形表达式为:
尹∥一B2以r2【1一毒J
式中d。一2d,.(8)
Taylor在不考虑土壤通道的通透性的前提下,导出了各相同性均质土壤的固有渗透率与土壤平均颗粒粒径平方以及孔隙率函数£3/(1一£2)成正比的结果[4].本文将在上述分析的基础上,利用分形理论,计及实际多孔介质粒径的分布和孔隙率分布,进一步考虑多孔通道通透性对固有渗透率的影响,提出固有渗透率的下列分形表达式:K一触2高等
式中A为比例系数,F,C为计算谱维数和剖面面积分维数时给出的常数值.
多孔介质单位面积渗流流量应遵循Darcy宏观经验定律,即:
口一K£墨grad^P=郴/c)2/伊一每,筹号籍爿弹㈤(10)
此式即为本文导出的多孔介质渗流流量表达式,与现有对均质多孔介质或容积平均多孔介质渗流流量计算式的不同点在于上式中的固有渗透率是考虑了多孔介质内部的不均匀结构后基于分形理论推导出来的,在多孔介质质量传递与多孔介质内部结构之间建立了有机联系.3多孔介质内部热量传递的分形模型
多孔介质中热量传递是一个复杂的过程.热量既可以通过复杂的固体骨架进行传导,也可以通过孔隙中的流体流动和流体的导热进行传递.如果固体骨架与流体之间存在温差,则它们之间还有导热和对流.一般说来,多孔介质内部的渗流是一个缓慢的过程(特别是对非饱和介质),作为近似分析,可以假设固体骨架和流体之间有热平衡.此外通过流体的热量迁移也可以假定是通过多孔介质孔隙中的气液两相混合物的当量传导进行的.通过上述简化后,多孔介质中的热量传递过程可看作是通过固体骨架和孔隙流体的一组并联热传导过程.
前人曾对多孔介质中的导热过程进行过许多研究,但是由于无法正确的描述内部多孔复杂结构,因此只对各相均匀排列,由相同大小颗粒组成的多孔介质进行过理论推算;或者对不规则结构的多孔介质按容积平均的孔隙率计算平均固体骨架和流体份额,然后进行导热计算,一】O一万方数据
施明恒,等:多孔介质传热传质分形理论初析
造成计算结果误差较大.分形理论的引入,使我们有可能对实际多孔复杂结果进行考虑,即把容积平均的当量导热系数改造成与内部形状直接联系起来的分形导热系数,从而较真实地反映多孔介质内部的导热过程.
为了推导多孔介质的分形导热系数,首先来计算多孔介质固相占有率的分形表达式.
对于一定孔隙率的多孔分形介质,其截面固相面积满足S(x)一似。则单位体积中固相占有份额为:
u一(导)3/2一c3/2x3“_1”2A(11)
作为近似,多孔介质内的导热通道被简化成通过气
隙(气液混合相)和固相的串并联通道,其热阻网络
如图4所示.通过多孔介质的总热阻为:
生!!二生::!±垄生::
“一k九(1一沪73+u)+蟹(p173一c,)P一
圈4多孔介质内的导热通道热阻示意图
最后得多孔介质的有效导热系数为:
~一i一——————永rj夕顷瓦丽F口丽伊砺而丽厂———一】1k九(1一c172x‘4—1’72+c372x3‘。一1’72)+砖(c1/2x‘4—1’/2一C372x3‘扣1’72)u纠,1o、式中砧,凡分别为多孔介质气液混合相和固相导热系数.
式(12)是根据分形理论导出的考虑多孔介质内部几何结构特征的有效导热系数计算公式,由它计算出的导热系数值比通常的容积平均计算公式更接近实际情况.
一般地说,多孔介质中热量的传递在不考虑辐射的状态下,对于未饱和多孔介质和饱和多孔介质来说有较大的差别,前者主要是导热为主.当内部液态介质有相变时(如干燥和冻结),则导热量主要消耗在介质内部的相变上.后者除了导热以外,还有流体在多孔介质孔隙通道中的渗流对流换热.通过单位多孔介质表面的热流密度可分别表示为
g。一九券(纯导热)
gz一九券+咿(导热+内部相变)
口。一九凳+^,(t,一o)(导热+渗流换热)
式中^,为分形渗流换热系数,将在另外的文章中进行分析和研究.
4结束语
本文对分形理论在多孔介质传热传质过程中的应用进行了初步的分析,求出了基于分形理论的多孔介质固有渗透率和有效导热系数.在此基础上,建立了多孔介质渗流与导热的分形模型.虽然本文的模型还是一种建立在分形参数上的当量传热模型,但本文的工作为进一步深人研究复杂的传热传质过程开辟了一条新路.
[参考文献]
[1]林瑞泰.多孔介质传热传质引论[M].北京:科学出版社,1995
[2]施明恒.多孔介质传热传质的进展和展望[M].中国科学基金,1995
[3]张济中.分形[M].北京:清华大学出版社,1995
[4]TaylorDW.FundamentalS0fSoilMechanics[M].JohnWiley&sons,Inc,1948万 方数据—11—
南京师大学报(工程技术版)第l卷第l期(2001年)[5]黑龙江省土壤普查办公室.黑龙江土壤[M].北京:农业出版社,1992
[6]江苏省土壤普查办公室.江苏土壤[M].北京:农业出版社,1995
[7]陈永平,施明恒.基于分形理论的多孔介质导热系数的研究[J].工程热物理学报,1999,20(5):605~612[8]陈永平,施明恒.基于分形理论的多子L介质渗透率的研究[A].中国工程热物理学会传热传质学术会议论
文集[C].苏州:1999
TheFractalTheory0fHeatand
MassTransferinPorousMaterial
ShiMingheng,ChenYOngpjn
(PowerEngineeringDepartmentofSoutheastUniversity,Nanjing,210096,PRC)
Abstract:Thispapermakesapreliminaryanalysisaboutthefractaltheoryusedinheatandmasstransferin
onporousmaterial;thecalculatingequationsofinfiltrationrateandeffectiveheatconductivitycoefficientbased
thetheoryareobtainedandthefractalmodelofinfiltrationandheatconductioninporousmaterialisestablished.Keywords:porousmaterial;heatandmasstransfer;fractaltheory
[责任编辑:刘健]
(上接第5页)
[参考文献]
[1]王树青,元英进.生化过程自动化技术[M].北京:化学工业出版社,1999
[2]王树青.生化反应过程模型化及计算机控制[M].杭州:浙江大学出版社,1998
[3]王树青.发酵过程控制的新进展[J].化工自动化及仪表,1995,22(4):3~10
ComputerIntegratedControlofAntibioticProductionProcess
WangShuq.rTg,ZhangQuan|.ng,ChenQ
(StateKeyLaboratoryofIndustrialCon”olTechnology,ZhejiangUniversity,Hangzhou。310027。PRC)
Abstract:Theresearchonintegratedcontrolofantibioticproductionprocessinvolvesmanipulatingtheinforma—tionaboutproductionprocess,marketinformation,managementinformationandexperimentaldata.Thekeysub—jectsofcomputerintegratedcontrol
toarediscussed,includingdatareconciliation,datamining,optimalcontrolandprocess.optimalmanagementimproveantibioticproductionandfermentation
Keywords:antibioticprocess;computerintegratedcontrol;integratedmanagementcontrol
[责任编辑:刘健]
一12一万方数据
多孔介质传热传质分形理论初析
作者:
作者单位:
刊名:
英文刊名:
年,卷(期):
被引用次数:施明恒, 陈永平东南大学动力工程系,南京,210096南京师大学报(工程技术版)JOURNAL OF NANJING NORMAL UNIVERSITY(ENGINEERING AND TECHNOLOGY)2001,1(1)18次
参考文献(8条)
1. 林瑞泰 多孔介质传热传质引论 1995
2. 施明恒 多孔介质传热传质的进展和展望中国科学基金 1995
3. 张济中 分形 1995
4. Taylor D W Fundamentals of Soil Mechanics 1948
5. 黑龙江省土壤普查办公室 黑龙江土壤 1992
6. 江苏省土壤普查办公室 江苏土壤 1995
7. 陈永平;施明恒 基于分形理论的多孔介质导热系数的研究 1999(05)
8. 陈永平;施明恒 基于分形理论的多孔介质渗透率的研究[会议论文] 1999
本文读者也读过(10条)
1. 陶涛. 陈剑波. 彭晓峰. 王峥. 颜岩 对复杂可变形多孔介质传热传质研究的思考[会议论文]-2001
2. 陈宝明. 王补宣. 张立强 多孔介质传热传质中耦合扩散效应的研究[会议论文]-2003
3. 李小川. 施明恒 多孔介质热导率的数值研究[会议论文]-2007
4. 李小川. 施明恒. LI Xiao-Chuan. SHI Ming-Heng 多孔介质热导率的数值计算[期刊论文]-工程热物理学报2008,29(2)
5. 陈永平. 施明恒 基于分形理论的多孔介质导热系数研究[会议论文]-1998
6. 金仁喜. 淮秀兰. JIN Ren-xi. HUAI Xiu-Lan 多孔介质高强度传热传质的理论研究[期刊论文]-工程热物理学报2005,26(z1)
7. 宫克勤. 孙苗苗 多孔介质中传热传质机理研究[期刊论文]-油气田地面工程2009,28(4)
8. 施明恒. 樊荟 多孔介质导热的分形模型[期刊论文]-热科学与技术2002,1(1)
9. 陈永平. 施明恒. CHEN Yongping. SHI Mingheng 基于分形理论的多孔介质渗透率的研究[期刊论文]-清华大学学报(自然科学版)2000,40(12)
10. 申盛. 刘伟. 陶文铨 多孔介质中自然对流传热的场协同分析[期刊论文]-自然科学进展2003,13(4)
引证文献(18条)
1. 朱群敏. 潘阳. 柏俊 竖直管蓄冰多孔介质局部自然对流的实验研究[期刊论文]-制冷与空调(四川) 2013(5)
2. 司新国. 张玉柱. 张庆军. 莫文玲. 刘丽妹. 马雪刚 低硅烧结矿质量分形维数影响因素[期刊论文]-河北理工大学学报(自然科学版) 2007(3)
3. 张赛. 陈君若. 刘显茜 基于分形和随机扩散描述食品材料的干燥速率[期刊论文]-材料导报 2013(10)
4. 韩博. 张凌. 王童. 张学东. 衣雪峰. 刘向飞 不同剪切作用下水合物钻井液微观结构试验研究[期刊论文]-地质与勘探 2013(4)
5. 韩博. 李承中. 范勇. 姜淼 水合物钻井液微观结构热传导分形模型研究[期刊论文]-地质与勘探 2012(4)
6. 蔡亮. 侯普秀. 虞维平. 张小松 水蒸气在霜层中扩散的分形模型[期刊论文]-中南大学学报(自然科学版)
2010(1)
7. 谢军. 王定标 分形理论在固体吸附式制冷系统中的应用研究[期刊论文]-制冷与空调 2007(6)
8. 李艳松. 张捍民. 孟凡刚. 张兴文. 杨凤林 分形理论在多孔介质渗透性能的研究[期刊论文]-过滤与分离 2006(1)
9. 马德刚. 张书廷. 柴立和 污泥干燥速度曲线的分形维数分析[期刊论文]-天津大学学报 2007(10)
10. 周济福 渗流力学研究的现状和发展趋势[期刊论文]-力学与实践 2007(3)
11. 韩博. 张凌. 蒋国盛. 宁伏龙. 刘力. 范勇. 胡代平. 辜幕丹 不同剪切作用下CMC水溶液的微观结构研究[期刊论文]-地质科技情报 2012(1)
12. 韩博. 张凌. 王童. 张学东. 衣雪峰. 刘向飞 不同剪切作用下水合物钻井液微观结构试验研究[期刊论文]-地质与勘探 2013(4)
13. 杨屹立. 杨玲. 徐武明. 朱留宪. 杨明金 多孔介质物料热风干燥研究进展[期刊论文]-农机化研究 2011(3)
14. 王李斌 多孔介质中的声场和流场及其数值模拟[学位论文]硕士 2006
15. 宋颖韬 多孔介质填充床中流体流动与化学反应过程分析[学位论文]硕士 2005
16. 彭鹏 石墨泡沫的制备及分形研究[学位论文]硕士 2006
17. 刘相东. 杨彬彬 多孔介质干燥理论的回顾与展望[期刊论文]-中国农业大学学报 2005(4)
18. 李艳松 膜污染因子识别与凝胶层结构分析[学位论文]硕士 2006
引用本文格式:施明恒. 陈永平 多孔介质传热传质分形理论初析[期刊论文]-南京师大学报(工程技术版)2001(1)