封闭方腔内空气和水自然对流每一次分岔数值预报 - 应用力学学报.pdf

第27 卷 第1 期 应 用 力 学 学 报 Vol.27 No.1 2010 年3 月 CHINESE JOURNAL OF APPLIED MECHANICS Mar.2010 文章编号：1000- 4939 （2010 ）01-00101-07 封闭方腔内空气和水自然对流 每一次分岔数值预报 1 2 3 王小华 朱文芳 魏英杰 （浙江大学 310027 杭州）1 （浙江科技学院 310023 杭州）2 （哈尔滨工业大学 150090 哈尔滨）3 摘要：采用二阶全展开ETG 分裂步有限元方法，通过对流动拓扑的详细分析，在排除网格密度影 响的基础上，结合二分法给出封闭方腔内空气和水两种典型流体自然对流发生第一次分岔时的临 界Rayleigh 数。计算结果表明，该方法可用于进行不同Prandtl 数条件下方腔内自然对流流动第一 次分岔的数值预报，可作为后续各阶分岔及转捩数值预报研究的基础。在相应的条件下，封闭方 腔内空气比水更容易发生分岔，且空气的流动结构相对于水表现出一定的倾斜性。 关键词：封闭方腔；自然对流；第一次分岔；数值预报；流动拓扑 中图分类号：O351 文献标识码：A 绝大多数的研究都是针对给定的几何边界条件 1 引 言 下，某个给定的Ra 范围内的流动特征展开分析， 或者分析几何以及物理边界条件对流动拓扑的影 近年来，由于封闭腔体的几何形状简单，且 响，并加入相应的流动控制手段整合流动结构以 有广泛的应用背景，但其流动机制复杂，大量学 达到某种实际的要求；文献[2]研究了不同边界条 者从物理实验和数值模拟的角度对内部的自然对 件和温度场对封闭腔体内自然对流的影响；文献[3] 流问题进行了深入的研究。这些研究成果从不同 数值分析了对应于Davis给定的物理边界条件下 角度不断加深人们对这一复杂非线性问题的认 的矩形封闭方腔内的自然对流问题。其结果均表 识，并且相关研究结果已在生产实际中得到应用。 明在发生流动向湍流的转捩后，对应Nusselt数的 现有的研究结果表明，随着 Ra(Rayleigh 数) 的不 变化会强烈依赖于相应的Ra ，且这一发生转捩过 断增加，封闭腔体内的流动会发生一系列的分岔， 程所对应的临界RaCr 的取值并不依赖于Pr(Prandtl 先是从一种定常层流转化为另外一种更为复杂的 数)[4-6]，但关于转捩的研究所给出的RaCr 的取值范 6 8 定常层流，再转化为按照某种周期变化的层流流 围却相当宽泛，在 10 ~ 10 之间变化；文献[7]研 动。最终随着 Ra 的进一步增加，发生转捩形成 究了底部边界加热、两侧边界和上部边界冷却条 复杂的湍流流动。流体流动经历了从稳定到局部 件下不同长宽比对矩形腔体内自然对流的影响， 失稳再到全局失稳的复杂过程。 其研究结果表明流动对应的Nu 是矩形腔体长宽 从文献[1]的数值分析开始，尽管通过大量的 比值的函数；文献[8]则分析了不同倾斜角度条件 数值和物理研究对上述过程有了定性的了解。但 下，矩形腔体内的自然对流流动，在分析中考虑