对空气在有泡沫金属介质管内流动与传热的研究.docVIP

《计算传热学基础》 对空气在有泡沫金属介质管内流动与传热的研究 热能与动力工程系10-1班 张皓威 雒飞 陈诚 白代立 指导教师： 黄善波 巩亮 徐会金 2013年7月 目录 题目 1 一、问题分析 3 二、解题过程 4 （一）对各个模型的流动和换热进行无量纲化 4 1、对各个模型的换热进行无量纲化 4 2、对各个模型的流动进行无量纲化 5 （二）Darcy模型的温度分布 7 （三）Brinkman模型的速度分布和温度分布 9 （四）Forchheimer模型的速度和温度分布 13 （五）Brinkman模型和Forchheimer模型的速度分布和温度分布进行对比。 17 （六）Brinkman模型和Forchheimer模型的f, fRe, Nu值 18 （七）总结 20 附录 21 附录1 计算Darcy模型的温度程序 21 附录2 计算Brinkman模型的速度和温度及Nu程序 23 附录3 计算Forchheimer模型的速度和温度及Nu程序 26 附录4 计算f和fRe的程序 30 参考文献 31 问题三十三（难度：5.0） 一根完全填充多孔介质管外表面为恒热流边界条件（），管内径为，平均流速为的空气在管内流动，其内部层流充分发展流动模型通常有Darcy模型、Brinkman模型和forchheimer模型，管内填充孔隙率为的多孔介质，渗透率表示为： 惯性系数表示为： 有效导热系数表示为： 充分发展的Darcy流动模型： （1） 充分发展的Brinkman模型： （2） 充分发展的forchheimer模型： （3） 质量守恒方程： 动量方程： （4） 式中，J=u/∣u∣是沿坐标轴方向的单位速度矢量；和μf分别为流体的密度和动力粘度；V为速度矢量；K，ε，CF分别为渗透率、孔隙率和惯性系数。动量方程右边的4项分别为：压力梯度、Brinkman项、Darcy项和forchheimer项。 能量守恒方程： （5） 参数： 试通过数值迭代计算： 将该问题的流动和换热进行无量纲化； Darcy模型的温度分布； Brinkman模型的速度分布和温度分布，计算该模型下的f、fRe、Nu数； forchheimer模型的速度分布和温度分布，计算该模型下的f、fRe、Nu数，并与Brinkman模型的速度分布和温度分布进行对比。 组员分工： 1、资料搜集：白代立 2、能量方程无量纲化离散：张皓威 3、Brinkman模型无量纲化离散：雒飞 4、Forchheimer模型无量纲化离散：陈诚、白代立 5、编程：雒飞、陈诚、张皓威 6、word录入：白代立 陈诚 7、ppt总结：张皓威 一 问题分析 （1）分析动量方程 题目的研究对象是充满泡沫金属管内充分发展段的空气，边界条件为恒热流密度，求在不同模型下的速度和温度分布。由于在垂直于轴向的圆形截面内，整个速度和温度分布呈轴对称，所以问题就变成求在圆形横截面内径向方向上的速度和温度分布。 由于流动过程为充分发展段，所以动量方程（4）等号左侧项等于0，可简化为： 根据上式可以得到不同模型的动量方程，即上面给到的式（1）、（2）、（3），所以在求每个模型的速度分布时可以直接对（1）、（2）、（3）式进行处理即可。 分析温度的边界条件 ·0 ·1 ·2 ···· · N-2· N-1· R · N 用方法A在径向从截面圆心到外边界取点如上图所示。 由元体能量法处理第二类边界条件(无内热源) 二 解题过程 对各个模型的流动和换热进行无量纲化 1．对各个模型的换热进行无量纲化： 由于三种模型的能量方程均一样，温度解均用同一能量方程进行编程计算，唯一不同的是：在darcy模型中U=um=1,在Brinkman和Forchheimer模型中计算各点温度需要将已计算出的各点速度带入U即可，均为无量纲量。 对能量守恒方程（5）处理： 其中温度的偏导为轴向的，可有以下得出： 在管子轴向取微元段： T T+dT 2H dx r X 能量守恒得： 所以: 所以（5）式可化