黏弹性流体在开口多孔介质中的自然对流与传热.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083148 黏弹性流体在开口多孔介质中的自然对流 与传热 牛骏，符策基，谭文长 （北京大学工学院力学与空天技术系，北京 100871） (Tel:010 Email:tanwenchang@pku.edu.cn) 摘要：基于修正的Darcy定律 ，数值模拟了上端开口，从底部等温加热的多孔介质中黏弹性流体的自 然对流与传热，其中使用了特征松弛时间λ和迟滞时间ε 来表征黏弹性流体的性质。我们选择了三组 不同的λ和ε 分别在Rayleigh数Ra≤400的范围内进行了一系列的计算。发现与牛顿流体中对应的情 况相比，黏弹性流体的传热特性有很大不同。另外，上端开口多孔介质中的黏弹性流体比上端闭口多 孔介质中黏弹性流体的自然对流与传热也有很大不同。在不同的Ra数范围，传热规律服从不同的标度 律。 关键词：黏弹性流体，多孔介质，自然对流，传热 0 前言 虽然由于粘性和非线性的存在，对多孔介质中流体热对流的研究十分复杂，但这方 面的研究从来没有停止过，因为它在很多领域都有广泛的用途，而且它是最简单的展现 非线性动力学的系统之一[1,2]。在现有的成果中，对牛顿流体在多孔介质中的对流与传 热已经做了很广泛的研究[3-5]。但至今为止对于非牛顿流体在多孔介质中的对流与传热 的研究还十分匮乏。但是对非牛顿流体的研究是十分必要的，因为非牛顿流体的对流与 传热在石油采收，生物工程等多个领域均有很大的用途[6,7]。在本文中我们通过使用修 正的Darcy定律对上端开口从底部加热多孔介质中的非牛顿流体进行了直接数值模拟， 并将结果与相同情况下牛顿流体以及上端封闭的多孔介质中非牛顿流体的计算结果[8] 进行了对比。 1问题的分析 1.1 控制方程组 考虑一个厚度为 d 的二维正方形多孔介质层，其顶端开口的，与外界相连通。垂直 的两侧面完全绝热，并且液体不能透过。上下两个等温端面存在温度差ΔT ，即整个多 孔介质层可看作是从下面加热的。假设多孔介质的渗透率为 K，并且其中充满了不可压 黏弹性流体，该流体的动力学黏度为μ ，热膨胀系数为β ，密度为ρ ，整个系统的热扩 散系数为κ 。通过运用修正的 Darcy 定律和 Boussinesq 近似，可以得到黏弹性流体在 多孔介质中运动的控制方程组[8]: ∇⋅V 0 (1) 资助项目：国家自然科学基金资助项目(No ⎛ ∂ ⎞ μ ⎛ ∂ ⎞ 1 λ p ρg 1 ε (2) ⎜ + ⎟(−∇ + k ) ⎜ =+ ⎟V ⎝ ∂t ⎠ K ⎝ ∂t ⎠ ∂T +V ⋅∇T κ∇2T (3) ∂t ρ ρ0 (1=−β T −( T 0)) (4) 这里V 是 Darcy 速度，p 是压强，g 是重力加速度，ε和λ分别是表征黏弹性流体的迟 滞时间和松弛时间。k 是沿垂直方向的单位向量，正方向向上，ρ 是在