充满多孔介质的长方形截面通道内充分发展对流换热问题的数值研究.doc

计算传热学课程设计（报告） 题 目  项目 每人任务分配 编程 芦旭红，齐尚超，杨鹏来，朱鹏 问题分析 朱鹏（图像），齐尚超（程序） 图像 杨鹏来 PPT制作 朱鹏 推导 芦旭红，齐尚超，杨鹏来，朱鹏 报告 芦旭红 热工一班 组长：朱鹏 组员：芦旭红，齐尚超，杨鹏来  目录 设计题目................................................3 1.1设计题目.............................................3 1.2已知参数.............................................4 物理与数学模型.. .................................... ..5 2.1物理模型.............................................5 2.2数学模型.............................................5 数值处理与程序设计......................................6 3.1数学模型无量纲化.....................................6 3.2数值求解.............................................8 3.3程序编写. ..........................................11 程序的验证.............................................12 计算结果与分析.........................................14 结论...................................................21 参考文献...............................................21 附录...................................................22 1 设计题目 （多孔介质，矩形a/b,单方程）水在一长方形截面的通道中进行充分发展的层流流动，该通道内充满多孔介质。多孔介质具有良好的强化换热能力，孔隙率ε是其基本结构参数，据此可以计算渗透率K，惯性系数CF，有效导热系数ke,具体表达式见[5]。其内部充满流体时的流动和换热通常采用体积平均法进行建模，即不考虑区域内孔的微结构而假定区域内任意一点处既有流体相又有固体相。由于金属泡沫的固体骨架导热系数较高，因此对于其内部的对流换热，通常采用局部非平衡模型，即考虑区域内流体温度和固体温度的差异。填充孔隙率为ε=0.6的多孔介质，渗透率表示为： 惯性系数表示为： 有效导热系数ke表示为： 沿流动方向的速度方程可以简化为 （1） 截面上的平均流速为wm=0.1m/s,dp/dz的值是恒定的，可以通过下式得到： （2） 其中，w为沿流动方向的速度。换热方程为： （3） 其中z为轴向。假设流动和换热都达到充分发展，外壁面为恒热流边界条件（qw=1000W/）,请基于局部热平衡模型，选取a=0.04m,b=0.02m,以20℃作为水物性的参考温度，参考数值传热学（陶文铨著）4.8节的内容 设计题目 分析孔隙率对渗透率K，惯性系数CF，有效导热系数ke的影响规律，计算随的变化关系； 通过能量守恒将方程（3）化为更简单的形式进而消去z（仅对恒热流条件下实施），对流动和换热方程进行无量纲化处理； 计算轴向无量纲速度的二维数值解，并计算摩擦系数f, 数，分析对f数的影响； 计算无量纲温度的的二维数值解，并计算截面流体平均温度和数； 1.2 已知参数 为了求得数值结果和利用结果进行分析，现给定题目相关已知量，长方形截面:a=0.04m,b=0.02m.填充孔隙率ε=0.6，热流密度qw=1000W/,平均流速，孔隙直径,20℃水物性参数（查表知）:=998.2kg/,=。 渗透率K（在一定压差下，孔隙结构允许流体通过的能力） 惯性系数（用于研究渗流，流体力学中常用的无量纲量） 压降（沿轴向方向由于增速降压时引发的压头损失） 计算 2 物理与数学模型 2.1 物理模型 图1