热肋片宽度对复杂方腔自然对流影响的数值分析.pdfVIP

第28卷第3期 计 算 力 学 学 报 vO1．28，NO．3 2011年 6月 ChineseJournalofComputationalMechanics June201l 文章编号 ：1007—4708(2011)03—0383—05 热肋片宽度对复杂方腔 自然对流影响的数值分析 王小华 ， 魏英杰 ， 朱文芳。 (1．浙江大学 航空航天学院 流体工程研究所 。杭州 310027；2哈尔滨工业大学 航空航天学院，哈尔滨 150090； 3．浙江科技学院 建工学 院，杭州 310012) 摘 要 ：采用笔者等发展 的二阶全展开 Euler—Taylor—Galerkin分裂步有限元方法 ，经过对算法精度验证 ．在确定不 同网格密度对计算结果影响的基础上，确定 出计算所 需网格密度 的必要条件 ，并在指定的网格密度条件下 ，分析 了热肋片宽度对封闭方腔 内典型 自然对流流动的影响，流动对应 的Pr=0．7l，Ra=10。计算结果表 明，随热肋片 宽度的增加，方腔平均 Nu数呈现 出先减小后增加的趋势，转折点肋宽约为 0．7，对应涡核从一个变为两个 。 关键词 ：热肋片；肋片宽度 ；复杂方腔 ；自然对流 ；网格密度 中图分类号 ：0357 文献标志码 ：A 1 引 言 作为从热环境控制问题抽象出来 的模型，封闭 方腔内的自然对流 以其简单 的几何形状和丰富的 对流流动信息，在被用作检验算法典型算例的同 时，其流动特征的分析 ，对揭示 自然对流这一复杂 非线性问题的基本规律具有重要的研究价值 ，同时 (a)不加肋片和具有热肋的计算域 封闭方腔内的 自然对流在工程领域中有着较广 的 应用背景，其研究一直是数值传热学和传热学 中的 经典课题[1]。关于方腔内各种流体 自然对流的数 霪 值模拟 ，已经取得了很多成果 ，其 中相当部分 已被 用于 自然对流换热的控制和优化的研究 。实 际过 程 中对方腔内自然对流的有效控制涉及到诸多 的 研究领域 ，如食物保存、结晶控制等等，确定相应构 件几何特性对温度场、速度场以及对应换热效应的 (b)网格剖分实例 图 1 计算域和 网格剖分实例 影响，具有重要的应用价值 。 本文应用笔者等发展起来的二阶全展开 ETG 2 计算方法 (Euler-Taylor—Galerkin)分裂步有限元方法 [1，分 2．1 控制方程 析了具有热肋的、二维的复杂方腔 中，肋片宽度对 由于耦合了浮力项的作用 ，对应控制方程的数 方腔中自然对流特征的影响，计算网格采用均匀的 值求解难度有所增加 。对应方程 的无量刚化整理 正方形单元网格 ，计算域和网格划分实例如图 1所 如下 ，取无 因次长度 、时间、速度和温度为 示。研究表明，热肋片宽度对方腔 内的自然对流特 一 性具较强的影响 ，随着热肋宽度 的增加 ，方腔 内对 一≠，’ 一一 ，’p户一一 流流动结构及传热效应会发生较大的变化 。 ￡一 ， 一 收稿 日期 ：2009—05—21；修改稿收到 日期 ：20lO一0709