在以幂律流体饱和的水平多孔层中的垂直通流的线性不稳定.docx

在以幂律流体饱和的水平多孔层中的垂直通流的线性不稳定 摘要 研究了非牛顿力学定律流体的垂直流动对水平多孔层中对流不稳定性的影响。扩展的达西动量扩散模型与Oberbeck-Boussinesq近似一起使用。垂直流动的固定基本解决方案是分析确定的。基本速度和温度场表现为与非牛顿流变学无关。进行线性稳定性分析，导致四阶特征值问题。使用特征值问题的数值解来获得中性稳定性曲线和不稳定性发生的临界瑞利数。转移到不稳定的控制参数是与通流相关的Péclet数和流体的幂律指数。这些参数影响中性稳定性曲线的位置以及临界瑞利数。详细讨论了绝对假塑性，绝对膨胀和小Péclet数的渐近情况。后一种情况导致简单的分析解决方案，当Péclet数小于1时，可以很好地近似数值数据。 关键字 多孔介质 非牛顿流体 垂直流 对流不稳定 瑞利数 介绍 地壳中的垂??直温度梯度，如果足够大，可能产生具有局部自由对流单元的流体流。因此，地下水的多孔流循环可以驱动不同地区的可溶性矿物。这可能导致所述的矿化问题。或化学污染物分散可以由地壳中地下水的自由对流驱动，从而产生含水层的严重污染。包括从下面加热并被牛顿流体或非牛顿流体饱和的多孔层的系统是用于研究这些现象的基本模型。在这种多孔层中产生的自由对流室可以与由压力梯度引起的强制垂直通流相结合。例如，在地球物理学中，净垂直通流可能由下地壳中的岩石静压力梯度引起。 在由Woodt，Sutton和Homsy和Sherwood的开创性论文中开发了从下面加热并受到垂直通流的水平多孔层中的对流不稳定性的建模。 通过采用基于达西定律和基于Oberbeck-Boussinesq近似的方案研究饱和多孔介质中的流动。 木材采用无定形厚度的层，而萨顿和霍姆斯和舍伍德认为有一个有限厚度的层。 通过线性分析以及通过能量方法非线性分析来研究垂直通流的稳定性。 后来，进行进一步研究调查不同的效果，如各向异性，非线性，内部发热，粘滞耗散，异质性和局部热不平衡。 关于这个主题的文献的详尽调查可以在Nield和Bejan的书中找到。 通过流动的主要特征是，如果与没有通过流动时在多孔层中发生的纯Rayleigh-Bénard不稳定性相比，其效果是稳定的。 Jones和Persichetti观察到当多孔层的边界面被认为是同量异位的并因此对于不稳定的对流辊是可渗透的情况时，这种行为的轻微偏离。 地球物理学中的许多问题包括非牛顿流体， 例如涉及碳氢化合物和油藏的那些。 在 这些情况下，分析了对流细胞在多孔中的发生 层意味着偏离古典达西定律，由于 非牛顿流体的流体。 一个典型的非牛顿 行为是一些的剪切变稀或剪切增稠性质 流体。 这种行为可以由Ostwald-de Waele建模 幂律[16]。 Christopher和Middleman [17]制定了一个 修正了当饱和流体时应用的达西定律 具有幂律行为。 几个进一步的研究修改 达西定律流体力学定律是从这以来进行的 开创性的纸，提供实验和 理论结果验证这个模型。 彻底的评论和 关于这一主题的参考书目可以在Shenoy的论文中找到 [18]和Di Federico et al。 [19]。 几个研究已经 最近出版的关于流动和流动的幂律流体的流动 多孔介质[20-30]。 本文的目的是扩展经典分析的 在水平多孔层中开始对流不稳定 垂直通过流入非牛顿流体的幂律 类型。 我们提到已经进行了类似的分析 一种粘弹性流体在Shivakumara的几篇论文中 Sureshkumar [31,32]。 我们的线性稳定性分析是意味着 确定不同幂律的中性稳定性曲线 索引，n和不同的Péclet数，Pe，与 垂直通流。 这种分析将进行数值， 即使进行中性稳定性数据的分析评价 对于绝对假塑性的极限情况en！ 0T，绝对 膨胀en！ 1，或非常小的Péclet数字。 特殊情况 的缺失通过流动?Pe?0π产生完美的协议 Barletta和Nield [21]的结论。 2.数学模型 我们考虑具有厚度H的平面水平多孔层， 饱和的Ostwald-de Waele（幂律）流体。 让 垂直轴z平行于重力加速度g，但是 具有相反的方向（参见图1）。 边界平面，z = 0 和z = H，是可渗透的，并保持等温在T = T0 + DT T = T0。 这里T是温度场，T0是a 参考温度和DT是正温差。 速度场u具有笛卡尔分量（u; v; w） 2.1控制方程 其中u是一致性因子，n是幂律指数。 在牛顿流体的情况下，我们具有n = 1并且一致性因子与动态粘度一致。 稠度系数的SI测量单位为Pa，对于牛顿流体，u的SI单位为Pa。在n 1的情况下获得假塑性流体行为，而由n 1描述膨胀流体。 当幂律流体饱和多孔介质时， 流动由热浮力驱动，通过Oberbeck-Bouss