多孔介质中二维对流传热的分叉α.pdfVIP

第 17 卷第 3 期 计 算 力 学 学 报 V o l. 17 N o. 3 2000 年 8 月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS A ugust 2000 文章编号:(2000) 多孔介质中二维对流传热的分叉 孔祥言, 　鹿蓬勃, 　吴建兵 ( 中国科学技术大学 力学系, 合肥 230027) 摘　要: 本文利用分叉理论研究了流体饱和的二维多孔介质从底部加热所引起的自然对流。用有限 差分方法确定对流的分叉进程; 揭示其模式转换机理及分叉对非正常流动图象形成的影响; 同时确 定了矩形截面宽高比与临界瑞利数的关系。还提出了一个判别分支稳定性的简明方法。 关键词: 多孔介质; 自然对流; 分叉; 稳定性; 非线性 中图分类号: O 346　　　文献标识码:A 0　引　　言 [ 1 ] 多孔介质中对流问题不仅有重要的理论意义, 而且有广泛的实际应用价值 。关于矩形截 面内多孔介质中自然对流的分叉问题, 近几年的研究取得了一些进展, R ileg W in ters[ 2, 3 ] 利 用有限元方法研究了底部加热问题。文献[ 4 ]用有限差分方法研究了瑞利数R a ≤100 的分叉 结构; 文献[ 5, 6 ]研究了介质内部有热源的自然对流情形; 文献[ 7 ]研究了这种对流的混沌结 构。 本文在文献[4 ] 的基础上作了进一步研究。主要是基于分叉理论, 即在代数方程和微分方 程中确定奇异点。基本思路是补充分叉点需满足的条件以扩展控制方程, 然后利用差分方法解 这些方程。把问题转换成利用线性稳定性分析来确定奇异点, 再利用延拓方法跟踪由奇异点引 出的非线性分支, 得出分支的路径。本文通过计算N esselt 数对R ayleigh 数的变化曲线, 进一 步刻画这些分支的路径。同时还对流动状态, 各种对流模式解分支的稳定性作了进一步研究, 建立了稳定性的充要判据, 避免了不确定性。 1　物理模型和控制方程 考虑一个二维矩形腔体, 内部有为单相牛顿流体所饱和的多孔介质, 矩形宽高比为W H = 。选取坐标系的原点位于矩形中点, x 和y 轴分别向右和向上。介质的孔隙度、渗透率、热容 量和热导率分别为 , K , (c) 和k , 流体的热膨胀系数、运动粘度、热容量和热导率分别为, s s ( ) ( ) ( ) ( ) , c f 和kf , 总的有效热导率和有效热容量分别为km = kf + 1 - k s 和 c m = c f + 1 1 ( ) ( c) 。设矩形四边不透水; 上下边界等温, 温度分别为 T - T , T + T ; 而侧壁 1 - s 0 0 2 2 绝热。并设固体与液体处于热平衡状态, 惯性项可以忽略,Bou ssin e sq 近似成立, 故对非