底部加热长方体腔内自然对流的非线性特性研究.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083449 底部加热长方体腔内自然对流的非线性 特性 1,2 *1 战乃岩 杨 茉 ( 1 上海理工大学动力学院，上海 200093 ；2 吉林建筑工程学院，吉林 130021) 021 Email：velon0066@ 摘要：本文采用具有QUICK差分格式的SIMPLE算法对底部加热三维封闭长方体腔内空气的自然对流 进行了数值计算，根据计算结果，探讨了方腔内流体流动与换热的静态分岔与振荡情况。计算结果表 明，在固定的几何尺寸和不同Ra 的情况下，当初始场不同时，会出现若干不同的解，即存在解的静态 分岔；在固定的几何尺寸和相同的初始场情况下，较小Ra 时流动与换热处于稳态，当Ra超过某一临界 值时，流动与换热就会发生非稳态振荡。并且随着Ra增加，经历从唯一稳定解到倍周期解到混沌解的 变化过程；当方腔的几何尺寸不同时，从唯一稳定解到倍周期解到混沌解的变化的临界Ra也发生变化。 关键词：自然对流；数值模拟；静态分岔；振荡 0 前言 本文关注了底部加热的长方体腔内空气的自然对流换热问题。它是Rayleigh-Bénard 对流的一种形式。Rayleigh-Bénard对流在流体机械，传热和技术领域有着广泛的应用 [1][2] 。而且由于Bénard对流存在振荡、分岔等各种非线性现象，也是动力学系统稳定性 和非线性领域的经典问题之一[3]，关于此问题已有许多相关的研究报道。本文主要关注 的是底部加热长方体腔内的自然对流换热中的非线性特性。Krishmurti R[4]研究了从底部 加热的水平水银层的自然对流换热，发现当Ra2400 时流动和换热由稳态转变为非稳态。 J.Prat和J.M.Massaguer对固体边界条件的Bénard对流进行模拟，也证明有混沌的存在 [5][6] 。Leong[7]也通过激光干涉实验发现：随着Ra 的增加，流场从轴对称转变成非轴对称。 本文的主要工作是通过对流动和换热的数值模拟探讨底部加热长方体腔内的自然对流 换热的振荡、分岔以及混沌等非线性特性。 1 物理问题及数值方法 本文研究的物理模型如图1所示，长方形封闭腔底面加热，顶面冷却，流体为空气 Pr ＝0.701，底面温度为T =1 ，顶面温度为T =0 ，T T 。宽高比L （＝W/H ）固定为4 ， h c h c z 长高比Lx （＝L/H ）分别为2 ，4 ，8。四周为固体壁面，取无滑移条件，U =V = W =0 。两 侧壁面绝热，Ra 的参考尺度取H 。 图 1 方腔示意图 本文引入Boussinesq假设后，其无量纲微分方程如下： uv ∇ U 0 （1） ∂U uv ∂P P γ 2 （2 ） + ∇ U U − + ∇ U ( ) 1 ∂F ∂X ( P γ • R a ) 2 资助项目：国家自然科学基金项目(No.；上海市重点学科建设项目资助(No. J50501) 通讯作者：上海市军工路 516 号，200093 ，Tel:021Email:yangm66@ ∂V uv ∂P P γ 2 Θ （3 ） +∇ =− + ∇ +