分形土壤导热模型建立和传热分析.docVIP

中国工程热物理学会 学科类别：传热学 学术会议论文 编号：123175 分形土壤导热模型建立与传热分析 杜欣，李素芬，东明 （大连理工大学，大连 116023） (Tel Email:yushuangdao@) 摘 要 本文以土壤源热泵为研究背景，采用分形理论，视土壤为多孔介质，从多孔介质结构形成的物理机制出发，详尽阐述了土壤分形导热模型的建立过程，并自行编写Matlab程序计算土壤的有效导热系数，得到了不同元胞下土壤有效导热系数与土壤孔隙率之间的变化规律，以及固气相导热系数比对土壤导热系数的影响规律。 关键词 分形理论；土壤；Sierpinski地毯；有效导热系数 0前言 热泵作为绿色空调技术日益受到普遍青睐，而在热泵众多的热源中，土壤源热泵因其运行的稳定、节能性及环保性而被称为最有效的采暖空调装置。针对土壤源热泵的研究，地下换热是其核心内容。以往对U型埋管在土壤中的换热过程，多采用CFD数值模拟方法。通过建立动量方程和能量方程来模拟饱和土壤的传热与传质过程，模型中有效导热系数采用方程式（1）的形式： （1） 式中，kgks为固相导热系数，Ф为土壤孔隙率[1,2]。采用上述方法计算有效导热系数，很难考察孔隙微观结构对有效导热系数的影响。 20世纪80年代, 分形理论开始被学者们用于研究多孔介质内热质传递现象[3,4] ,国内外学者通过大量的实验和计算证明实际多孔介质孔隙结构在一定尺度范围内具有分形的特征[5,6]。分形概念的引入为多孔介质传热传质研究提供了新的方法和手段, 从而有可能更深刻地揭示多孔介质热质传递过程中的现象和规律。 土壤是由大小不同、形状各异的固体骨架和孔隙以一定的形式连接形成的多孔介质。土壤组分、形状及其结合方式决定着土壤的结构，从而进一步影响其有效导热系数。本文基于分形理论构造了不同类型的分形结构来模拟多孔介质,讨论了孔隙率、固气相导热系数等因素对分形多孔介质有效导热系数的影响规律。 1 分形导热模型的建立 1.1 多孔介质分形自相似理论 分形的一个基本特征是具有标度不变性，也就是说，对于分形，不论将其放大或缩小，它的形态、复杂程度、不规则性等各种特性不会发生变化。研究表明，实际多孔介质在一定局域长度范围内呈现出分形特性。 由于多孔介质是由大小不同的颗粒逐渐堆积而成，这与分形的构造过程十分相似，因此决定了多孔质材料的孔隙和粒子可用分形模型进行定量描述。机理是利用分形的构造过程即对初始元胞反复复制生成的原理，来模拟多孔介质中孔隙和粒子的空间构形以及粒子的扩散和聚集。 1.2 分形导热模型 本文采用随机Sierpinski地毯模型模拟自然界真实的多孔介质，这是一种“人工构造分形体”，自然界真实的多孔介质是复杂多变的，构造随机Sierpinski地毯的方法多种多样，在文献[7]中作了详细的阐述。 本文使用类似文献[8]的方法来构造：首先建立N×N的矩阵（N=3，4，5...），用来储存多孔介质的信息，并以数值来代表整个多孔介质基质和孔隙随机分布的信息。例如在矩阵中，1代表基质，以黑色来表示，0代表孔隙，以白色来表示，此时N=3，孔隙率为0.56。以n=3为例：把一个单位正方形，分成九个边长为1/3的小正方形，其中如果被基质占据，那么就标志为1，如果为孔隙所占据，那么就标志为0，类似地，对于上述九个正方形中的每一个，如其为l，那么重复以上操作继续细分；如其为0，则停止细分。这样就构成一个含递归关系的随机分形模型，同时利用matlab编程语言生成分形多孔介质的模型，将信息储存在矩阵中，并根据这些数据集合逐级绘出图形。 本文中采用两种方法建立随机Sierpinski地毯模型，分别是基质生成法与孔隙生成法，鉴于随着孔隙维数P（P代表每级生成的孔隙数量）的变化，模型会与土壤的实际情况相偏离，如P过小时模型中孔隙偏少，所以采用基质与孔隙生成相结合的方法可以更加贴近土壤的实际情况，基于以上模型，本文分析了孔隙率，多孔分形结构以及固气相导热系数比对土壤有效导热系数的影响规律。 1.2.1 基质生成法 绘制Sierpinski地毯的matlab算法步骤如下 (l)设置孔隙维数P，矩阵维数n，矩阵总元素数N×N（N=64;n=4代表64×64的矩阵，最小元胞为4×4） (2)一级分形：用黑色随机填充（16-P）个正方形，作为多孔介质的基质，设其为1;其他的正方形为白色，为多孔介质的孔隙，设其为O;这些