相变是传质传热中的一种典型情况。而多孔介质中物质的相变在工业生产和日常生活中十分常见。由于相变问题的非线性特点,需要格子Boltzmann方法这种能够模拟流体运动以及为复杂物理现象建模的数值方法来处理它。本文运用格子Boltzmann方法对纯相变以及多孔介质内物质的对流相变过程进行了数值模拟,并对数值结果加以研究和讨论。发现不同的达西数对多孔介质内的对流相变过程中有重要的影响。 Phase change is a basic problem in heat and mass transfer field, especially in porous medi a . In this work, as it’s nonlinear characteristics, Lattice Boltzmann method (LBM) is used to numerical y simulate the h eat and mass trans- fer coupling with phase change in porous media, and result s showed that different Darcy Numbers have some certain effects on phase change process in porous medium.
多孔介质中质量、动量和能量的传递现象[
人们对多孔介质中能量、动量和质量传递现象的关注和研究由来已久。最初的研究领域是地下水勘探与预测。早在1856年,达西(H. Darcy)曾对地下水源进行了研究,提出了著名的适用于一定条件下多孔介质中流体流动的达西定律[
由融化过程中的自然对流是一个缓慢的过程。液相区的自然对流是由重力作用引起的。Boussinesq近似后,多孔介质中对流相变的控制方程写为(Guo[
(1)
(2)
(3)
方程(1)(2)(3)式中,是流体的速度矢量,p是压力,是粘性系数,是多孔介质的孔隙度,是热扩散系数,是热膨胀系数,纯净物与多孔介质的热熔比,c为热容比,它们的下标分别表示固体、流体和恒压,这里等于1。另外,外力项为
(4)
(4)式中的第一项表示流体与多孔介质之间的相互作用力,第二项代表几何结构的影响,最后一项代表重力作用。
(5)
是相变过程中的热源项。
这里有3个描述控制方程的无量纲参数:
(6)
(6)式中为达西数,为普朗特数,为瑞利数。其中l为特征长度。
图1为封闭方腔内多孔介质对流相变示意图。方腔内填充了孔隙度为的多孔介质。的值恒定不变。方腔左壁的温度始终保持在T1。在这个多孔的方腔内充满了温度为的固态纯净物,该纯净物在常压下的相变温度为。上下壁及右壁绝热。图1中的黑色曲线为相变界面,它将方腔内的纯净物分为固相和液相两个区域。
该模型在速度场的LB方程可以写为[
(7)
是粒子的密度分布函数,下标i表示速度方向,x是粒子的位置,t是时间,是时间间隔,是速度场的弛豫时间。是粒子的平衡态密度分布函数。
图1. 多孔介质方腔内对流相变示意图
在D2Q9模型中,离散的粒子速度是给定的:
(8)
考虑到多孔介质对流体流动的影响以及其他外力的作用,公式(7)中的被表示为:
(9)
平衡态密度分布函数为
(10)
权重,,。流体的密度被定义为
(11)
粘滞系数为
(12)
经多尺度展开验证,LB方程(7)式是正确的。
在多孔介质中温度场的LB方程与对流的影响耦合,因而温度场的LB方程为
平衡态温度分布函数改为[6]
(13)
热传导系数为
(14)
这里和分别是格子速度和宏观流体流动速度。为温度场的弛豫时间
以上LB方程经多尺度展开[4,7,8]后便可得到宏观能量方程式(3)式。
在LB方法中,边界条件的选取将直接影响计算的精确性和简洁性,选择一个适当的边界条件将使结果更加可靠。在这里,我们选用了非平衡态外推条件来处理速度场的相变界面和二维方腔的四个壁面。假如格点b是边界上的点,格点f是流场中的点,f点和b点的位置关系是。速度是已知的,密度是未知的,所以b点分布函数就可以从下式得到
(15)
温度场的边界条件与速度场相同,因此b点温度分布函数为
(16)
对于融化过程,能量的传递和流体的流动主要依靠热传导和重力作用。达西数是表征流体流过多孔介质时,其流速与流动方向上的压力梯度成正比趋势的物理量,在流体研究中具有重要意义。为了研究Da对流动的影响,其他参数设定为:,,,,,。该纯净物的相变温度Tm等于0,相变潜热L等于0.25。以下是我们讨论分析了达西数在具体案例中的影响情况。
图2是模型的等温度和流线图,从上到下达西数Da依次是、、、、。五张等温线图中标记的在Y=3L/5处的点的温度等于相变温度。五张流线图中的空白处标明此区域为固相区,没有发生对流。如图2(左1)所示,当时等温线几乎是垂直的,这就意味着在低达西数下,热传导对能量传递起主导作用,而对对流的作用不明显。随着达西数Da的增大,等温线的弯曲形状就越加明显。在高达西数下,等温线弯曲已经十分明显,这意味着在高达西数下对流起主要作用。
相变界面的位置同样是相变研究中最为关心的问题。Huber[
图2. 在= 0.6, Pr = 1.0, Ra =1.0e + 5下的等温线图(左)和流线图(右):从上至下Da = 5.0e − 5, 5.0e − 4, 5.0e − 3, 5.0e − 2, 5.0e − 1
怎样变化呢?图3是达西数分别为,,时不同达西数下,温度沿着位置的横线变化情况。可以看到后随着达西数的增大温度也随着增高,但是在前,可以看到随着达西数的减少温度也增高的现象,而这正是对流对温度的影响,与我们之前在2.2.1中得出的结论相吻合。还可以看到在温度等于“0”位置,也就是相变界面,随着达西数的增大,相变界面的移动速度逐渐加快。
图4(a)、(b)分别是位置的竖直线上速度和温度的变化情况。可以看到达西数的大小不仅会影响热传导,同样也会影响对流。在图4(a)中,速度在方腔顶部和底部附近最大和最小两个峰值。随着达西数Da的增大,两个峰值会变得更加陡峭。对流在顶部和底部最为活跃。在中间区域,速度随着高度h的增加而变大。在图4(b)中,方腔底部的温度是最低的。温度随着高度的上升而上升,在高度时到达最大值。图4(a)和(b)解释了相变液相区中较热的流体上涌,较冷的流体由于重力作用下沉的现象。这个运动导致了速度产生最大和最小两个峰值,以及底部流体最冷,顶部最热。
格子Boltzmann模型适合应用于多孔介质中对流相变问题。考虑到重力和多孔介质的外力的作用,我
图3. 在200 × 200的格子不同达西数下,温度沿着y = 3L/5位置的横线变化情况
图4. 在200 × 200的格子中,速度(a)和温度(b)
沿着x = L/4位置的竖直线的变化情况
们修改了LB方程和平衡态分布函数。根据模拟结果得出如下结论:
1) 随着达西数Da的增大,等温线的弯曲形状将越加明显;
2) 在高达西数Da下,融化过程中的自然对流传热有重要影响;
3) 达西数同样影响了对流和传热。当达西数增大时,相变界面向右移动变快;
4) 融化过程中,较热的流体上涌,较冷的流体由于重力作用下沉。这个运动导致了速度产生最大和最小两个峰值,以及底部流体最冷,顶部最热。
[
[
[
[
[
[
[
[