三维方腔自然对流的非结构同位网格数值模拟研究.pdfVIP

三维方腔自然对流的非结构 同位网格数值模拟 赖锡军1 汪德耀2 (1．中国科学院南京地理与湖泊研究所，南京210008； 2．河海大学环境学院，南京210098) 摘要 以基于Boussinesq假定的无量纲化的方程为自然对流计算的控制方程，采用非结 构同位网格上的SIMPI。EC算法计算了三维方腔内水和空气两种介质在受侧壁加热情况下的 自然对流问题。计算的对流项同时采用了非结构网格中的二阶逆风格式和二阶中心差分格式 近似。计算得到了稳态的三维流动结构和温度分布。而且，随着Ra数的增加，对流作用增 强，三维流动特征更为显著。 关键词 三维，自然对流，非结构同位网格 1 前言 在重力扬、离心力场或其他力场的作用下，由于流体的温度差或(和)浓度差形成密度差和 浮升力，使流体产生流动的现象称为自然对流。它有别于由外力(如泵、风机或飞行器的运动 等)而引起的流体受迫对流。自然对流及其换热在自然界和工程上有着广泛的应用。例如，大 气或海洋中的循环，变压器中油的循环，由于室内的供热或制冷而导致空气的自然 对流，某些电机或核反应堆中的自然对流冷却，太阳能集热器中的自然循环等等。 由于人们对大气和海洋中的自然对流现象以及它们在工程中的应用越来越重视，因此自 然对流的研究得到迅速的发展[1’4]。为了增强模型模拟复杂区域内流动和换热问题的能力，非 结构同位网格算法受到了很大的重视[5]。本文分别选择水和空气两种介质，采用非结构同位 网格上的SIMPI。EC算法[6]计算了三维方腔内的自然对流换热过程，得到了稳态的三维流动 结构和温度分布。 自由对流的控制方程 自然对流是在重力场或其他力场的作用下由密度差引起的浮升力产生的，因此动量方程 式中的体积力项不允许忽略。而且，密度差常常是由于温度差或者物质混入引起的。因此，自 然对流中的流动和传热、传质相互耦合，动量方程和温度(浓度)方程必须同时求解。由 BOussinesq假定，可写出自然对流的控制方程组： 挚：o (1) dZi 瓦十吩瓦2一石茬十口蔼十p松』gi玺+吩差=一寿差+口丢亳+阳您i (2) 瓦十哟磊2u丁瓦i云 (3) 雾+哟嚣=D丁器 为了便于分析计算，通常对上述方程进一步无量纲化。以流场的特征长度L和特征温度 △丁为基本量。并取参照速度Uo—Dr／L，参照时间％一L2／D丁，可得如下无量纲量：xi=詈， u一番，P=毒每，口=尝秀，虿i=詈代人可得方程的无量纲形式： 型一n (4) aX；一。 (5) a￡‘。’aXi aXi’～aXiaX，‘”1’惜 警+u爰一蓑+R爰‰+鲫rz磨 一十-，一=一～+rr一+，∥Pr。W 雾+u爰一鑫 (6) 其中，Pr—v／D了’，为Prandtl数；Gr=艮g△丁L 数。g为重力加速度。 3数值解法 3．1对流扩散方程的离散 通用变量西的输运方程的积分形式可写表达成： 毫，∥V+J．A(聊～脚例·眦一，vs∥V (7) 以任意形状的凸多边形作为网格单元(如图1)，可得到用有限体积法离散的方程‘6|，即：