7.19多孔介质边界条件
多孔介质模型适用的范围非常广泛,包括填充床,过滤纸,多孔板,流量分配器,还有管群,管束系统。当使用这个模型的时候,多孔介质将运用于网格区域,流场中的压降将由输入的条件有关,见Section 7.19.2.同样也可以计算热传导,基于介质和流场热量守恒的假设,见Section 7.19.3.
通过一个薄膜后的已知速度/压力降低特性可以简化为一维多孔介质模型,简称为“多孔跳跃”。多孔跳跃模型被运用于一个面区域而不是网格区域,而且也可以代替完全多孔介质模型在任何可能的时候,因为它更加稳定而且能够很好地收敛。见Section 7.22.
7.19.1 多孔介质模型的限制和假设
多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流动阻力。本质上,多孔介质模型仅仅是在动量方程上叠加了一个动量源项。这种情况下,以下模型方面的假设和限制就可以很容易得到:
因为没有表示多孔介质区域的实际存在的体,所以fluent默认是计算基于连续性方程的虚假速度。做为一个做精确的选项,你可以适用fluent中的真是速度,见section7.19.7。
? 多孔介质对湍流流场的影响,是近似的,见7.19.4。
? 当在移动坐标系中使用多孔介质模型的时候,fluent既有相对坐标系也可以使用绝对坐标系,当激活相对速度阻力方程。这将得到更精确的源项。相关信息见section7.19.5和7.19.6。
? ?
当需要定义比热容的时候,必须是常数。
7.19.2 多孔介质模型动量方程
多孔介质模型的动量方程是在标准动量方程的后面加上动量方程源项。源项包含两个部分:粘性损失项(达西公式项,方程7.19-1右边第一项),和惯性损失项(方程7.19-1右边第二项)
(7.19-1)
式中,si是i(x,y,z)动量方程的源项,是速度大小,D和C是矩阵。动量源项对多孔介质区域的压力梯度有影响,生成一个与速度大小(速度平方)成正比的压降。
对于各向同性多孔介质简单情况下:
(7.19-2)
式中是渗透性系数,对角上的系数分别为
是惯性阻力系数,也就是将D,C矩阵简化为对角矩阵,和
,其它元素都是0.
同样fluent也可以将源项设定为速度的幂函数型:
(7.19-3)
式中and 是用户自定义的经验系数。.
在幂函数型模型中,压降是均匀的,多孔介质中的达西定律
的单位是国际单位制。
通过多孔介质的层流,典型的压降是与速度大学成正比,常数C2可以认为是0。
忽略对流加速度和扩散,多孔介质的动量方程源项就可以化简为达西定律:
(7.19-4)
坐标轴三个不同方向的压降fluent计算如下:
(7.19-5)
式中
,
是方程 7.19-1中矩阵, and
的项,
是 , , and 方向速度,
是 , , and 方向的多孔介质厚度.
这里, 多孔介质的厚度( , , or ) 是模型中的实际厚度. 如果
进行调整。
模型中的厚度不是实际厚度,就需要对输入参数
多孔介质材料惯性损失
在高速流动时,方程7.19-1中常数C2是对多孔介质中关于惯性损失的修正。这个常数被认为是流动方向单位长度的损失系数,压降定义为动水头的函数。 如果模拟多孔板或者管束系统,有时候可以忽略渗透项,只使用惯性损失项,就得到如下的多孔介质压降方程:
(7.19-6)
如下x,y,z方向的压力损失项: