多孔介质中粘弹性流体热对流的线性稳定性分析.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：093135 多孔介质中粘弹性流体热对流的线性 稳定性分析 牛骏1，符策基2 ，谭文长3 （北京大学工学院力学与空天技术系，北京 100871） (Tel: 010 Email: tanwch@pku.edu.cn) 摘要：对底部受来流加热、上端等温多孔介质中粘弹性流体进行了线性稳定性分析。引入了毕奥数Bi 描述底部加热条件。使用了特征松弛时间λ和迟滞时间ε表征粘弹性流体性质。我们选择了三组不同 的 −4 4 λ和ε ，对三种情况中的临界瑞利数 进行了数值查找，查找范围为边长0 到4， 从 到 ， R B i 10 10 c 底部边条件从接近等温情况变到接近等热流情况。结果显示随着B i 增加，震荡对流被推迟。松弛时间 与迟滞时间比值越大，对流产生越早。流动模式仅与迟滞时间有关。扰动温度等值线图显示，其样式 与B i ，边长a 和弹性参数有关。 关键词：粘弹性流体，多孔介质，线性稳定性，传热 0 前言 虽然由于粘性和非线性的存在，对多孔介质中流体热对流的研究十分复杂，但这方 面的研究从来没有停止过，因为它在很多领域都有广泛的用途。在现有的成果中，对多 孔介质中牛顿流体的稳定性已经做了很广泛的研究[1-3]。但至今为止对于非牛顿流体稳 定性的研究还十分匮乏，仅仅只有对第一类、第二类边条件情况的研究，但这些完美的 条件在现实中是很难找到的。Kubitschek 等通过引入毕奥数对多孔介质中的牛顿流体底 部受来流加热的情况做了研究，它们发现，通过改变毕奥数的取值，整个讨论过程能够 包含第一、第二、第三类边条件。但是对于非牛顿流体，还没有这方面的研究工作。在 本文中我们对上端等温底部被来流加热多孔介质中的非牛顿流体进行了线性稳定性分 析，通过引入毕奥数研究了不同的底部加热模式，得出了震荡对流产生时的临界瑞利数 值及流动模式，并对影响对流产生的因素做了分析。 1 问题分析 * 考虑一个高度为 H,长度为 的二维正方形多孔介质层。垂直的两侧面完全绝热， a 并且液体不能透过。上端等温，下端被温度为T * ，对流换热系数为h 的无穷来流加热。 ∞ 假设多孔介质的渗透率为 K，并且其中充满了不可压粘弹性流体，该流体的动力学粘度 β 为 ，热膨胀系数为 ，密度为 ，整个系统的热扩散系数为κ 。通过运用修正的Darcy µ ρ 定律和Boussinesq 近似，可以得到粘弹性流体在多孔介质中运动的控制方程组[4]: * * 1 ∇⋅v 0 (1) 基金项目：国家自然科学基金资助项目（No10572006）  ∂  * * * ˆ µ  ∂  * (2) 1 p g 1