两个自然对流换热问题的LB方法模拟.pdfVIP

半吲，??半舟，。＼?‘。／＼。’彳??△???△，?会?????蝗齀?／：?，?一彳。??瑀??，驰桃㈣川刈卜如∽沪矿∽??两个自然对流换热问题的?方法模拟【??!#瑂?——，?赵婷杨茉??档パ骞?余敏选用??取緊?岢龅腡脱?模型，密度分布函数?拔露确植己?鼺演化方程为：其中?豪?为粒子迁移速度，血，出分别是格子宽度和时间步长。前言?方法计算选用模型???簔???????格子玻尔兹曼????．?????琇?方法作为模拟流体流动与换热的一种新的数值方法发展得十分迅速。?方法从微观层次描述流体内部流动，将连续介质看作大量位于网格点上的离散流体质点粒子，粒子按碰撞和迁移规则在网格间运动，通过对各网格点上粒子运动特征的统计以获得流体宏观运动规律。由于?方法将微观和宏观规律联系起来，提供了一种不同于基于连续介质假设的???．???方程的描述和求解流动和传热问题的新思路，并被认为可能解决目前还较难处理的跨尺度流动和传热问题，因此引起了学术界的极大兴趣。已有许多?方法的成功算例，包括二维?????流?А⑾罡乔?动流【?、绕流圆柱【?、有自由表面的流动?】等。但是，?方法作为一种较新的方法，目前尚处于发展和被考核阶段，特别是对于有换热问题的模拟。人们正在用它模拟多种问题，提供新的算例，考察其适应能力和发展新的模型及数值处理技术。方腔内竖直平板的自然对流换热问题由于在计算区域内部存在一个孤立体，数值处理上较复杂，需要特殊考虑。底部加热的方腔内自然对流换热问题由于其特殊的加热和流动条件，能够表现出分岔等各种非线性特性。文献??肔?椒ǔ晒Φ囟哉庵治侍饨?行模拟，并获得了振荡解。但目前还没有采用?方法预测这种问题存在的多解性问题。本文采用?方法模拟方腔内竖直平板的自然对流换热和底部加热的方腔内自然对流换热问题，对孤立体进行了数值处理，预测出了自然对流问题的多个解和振荡解，并给出了振荡解的速度相图。二维空间中微观粒子速度??头植己?齔离散为图??镜?个方向：?，??—?万．．?—?万、??“??讯作者：．?汉＞黅路??牛????琓?：???????珽??：???????畁?中国???程熟物理学会学术会议论文传热传质学编号：???上海理工大学动力工程学院上海?．???????图?肷⒎较?国家自然科学基金重点项目????????????△?．， 吼炉哆怿??孚。???海桓?矿∽忙和????兰?萑???扣一??，等。?垡籵．????蟶?硓?琭?以功?以曩。?硓?珼，“五??黴?√???莩?琾㈨???盯尝??材???五号?，??刈????甜???蜀?弧ぁ隆鞫雚·??丁一瓦??渲械眆??时，％?；否则％?。江。数值计算与结果????齇???，??，???封闭腔内竖直平板的自然对流换热方形封闭腔内有一宽???、高??且温度为死?的孤立竖直平板，方腔的左右两侧壁温为乙?，上下壁面均为绝热。缸???琍????，特征流速选取为图??给出了???、?肚襌噝萨??????盠?椒?Ｄ獾?哆，吼无?悄Ｐ筒问?憾?宏观的密度、速度、温度及压力定义如下：式??械穆辣硎居芍亓σ?鸬牧鞫?牖蝗锐詈舷睿?泶锸轿#?通过?印??瓻???喑叨日箍7椒n?梢酝频汲龊旯鄣腘．?匠獭?】，且粘性系数及导温系数表示式如下：图?境隽朔獗涨荒谑?逼桨宓淖匀欢粤骰蝗任侍狻Ｈ缤?所示，边长??的正矿?∥．??鞫。????，△丁?咭籞，本文处理计算域中孤立体时采用如下方法：孤立体表面宏观速度均为???龉铝⑻迩?蚰诘姆植己??唤?醒莼??且在流场与孤立体界面上采用无滑移边界条件，边界格点上不进行碰撞，但入射粒子同转??沿原路返同。剑的流线图。为了比较，图??和??直鸶?隽吮疚挠镁哂?????址桨傅腟???椒?Ｄ馔??侍獾玫降牧飨咄????盯????／??．?×??????，??图?校患乐薄?板的封闭腔声瓦．??“；?。”， ￡??：?瑈?骸危甮．△丁．日：???，?裕后纬浮?底部加热的方腔内自然对流换热问题如图??尽７角豢砣???叩???闹芪???底部加热的方腔内自然对流砌?×??那樾谓?屑扑悖???愎怀さ氖奔洌?钪盏慕峁?鞒【?锏轿榷ㄗ刺?比?砌?×???，?秩??，??砌???，且奇数表示?醭。?际?硎?初场。如图?????，??一定态解。??和??分别表示当尺萨???和尼??时以?醭『?初场计算所得和文献?】给出的同一问题流场的烟可视化实验结果。由图?杉?街质?捣椒?笆笛榻?果一致，表明?方法和本文对孤立体的数值处理技术能够获得正确结果。图??隽耸?逼桨逶诓煌?恢檬钡任孪叩氖?到峁?ＮA薊匕较，本文还用???方法???袷浇?辛耸?的Ｄ狻ＮA私谑∑7?瑂???椒ǖ慕峁?饫镂戳谐觥１?较表明，两种方法的模拟的等温线结果非常一致。固体壁面，取无滑移条件。两侧壁面绝热，下壁面温度为死?，上壁面温度为疋?。水平方