多孔介质参考资料

《多孔介质参考资料》由会员分享，可在线阅读，更多相关《多孔介质参考资料（13页珍藏版）》请在金锄头文库上搜索。

1、多孔介质条件多孔介质模型可以应用于很多问题，如通过充满介质的流动、通过过滤纸、穿孔圆盘、 流量分配器以及管道堆的流动。当你使用这一模型时，你就定义了一个具有多孔介质的单元 区域，而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定。通过介质的热传 导问题也可以得到描述，它服从介质和流体流动之间的热平衡假设，具体内容可以参考多孔 介质中能量方程的处理一节。多孔介质的一维化简模型，被称为多孔跳跃，可用于模拟具有已知速度/压降特征的薄 膜。多孔跳跃模型应用于表面区域而不是单元区域，并且在尽可能的情况下被使用(而不是 完全的多孔介质模型)，这是因为它具有更好的鲁棒性，并具有更好的收敛性。详细内容

2、请 参阅多孔跳跃边界条件。1、多孔介质模型的限制如下面各节所述，多孔介质模型结合模型区域所具有的阻力的经验公式被定义为“多 孔”。事实上多孔介质不过是在动量方程中具有了附加的动量损失而已。因此，下面模型的 限制就可以很容易的理解了。 流体通过介质时不会加速，因为事实上出现的体积的阻塞并没有在模型中出现。这对于过渡流是有很大的影响的，因为它意味着 FLUENT 不会正确的描述通过介质的过渡时 间。 多孔介质对于湍流的影响只是近似的。详细内容可以参阅湍流多孔介质的处理一节。2、多孔介质的动量方程 多孔介质的动量方程具有附加的动量源项。源项由两部分组成， 一部分是粘性损失项 (Darcy)，另一个是

3、内部损失项：S 二工 D pv + 工 C p v viij jij 2 j jj=1j=1其中S_i是i向(x, y, or z)动量源项，D和C是规定的矩阵。在多孔介质单元中，动量损 失对于压力梯度有贡献，压降和流体速度(或速度方阵)成比例。 对于简单的均匀多孔介质：S = v + C p v vi a i 2 2 j j其中a是渗透性，C2是内部阻力因子，简单的指定D和C分别为对角阵1/a和C2，其它 项为零。FLUENT 还允许模拟的源项为速度的幂率：其中C_0和C_1为自定义经验系数。注意：在幕律模型中，压降是各向同性的，C_0的单位为国际标准单位。3、多孔介质的 Darcy 定律

4、通过多孔介质的层流流动中，压降和速度成比例，常数 C_2 可以考虑为零。忽略对流加速 以及扩散，多孔介质模型简化为Darcy定律：Vp =-巴va在多孔介质区域三个坐标方向的压降为：ApxApy=工巴v Ana j x j =1 xj=工旦v Ana j y j =1 yjAp =乞v Anz a j z j =1 zj其中l/a_ij为多孔介质动量方程1中矩阵D的元素，v_j为三个方向上的分速度，D n_x、D n_y、以及D n_z为三个方向上的介质厚度。在这里介质厚度其实就是模型区域内的多孔区域的厚度。因此如果模型的厚度和实际厚度不同，你必须调节1/a_ij的输入。.4、多孔介质的内部损

5、失 在高速流动中，多孔介质动量方程1 中的常数 C_2 提供了多孔介质内部损失的矫正。这一常数可以看成沿着流动方向每一单位长度的损失系数，因此允许压降指定为动压头的函 数。如果你模拟的是穿孔板或者管道堆，有时你可以消除渗透项而只是用内部损失项，从而得到 下面的多孔介质简化方程：並=工 C 1 p v|v|dx2ij 2 j ji j=1写成坐标形式为：Apx=工 C An2 xj xj=1ApyApz=工 CAn2 yj xj=1=工 CAn2 zj xj=15、多孔介质中能量方程的处理对于多孔介质流动， FLUENT 仍然解标准能量输运方程，只是修改了传导流量和过度项 在多孔介质中，传导流量

6、使用有效传导系数，过渡项包括了介质固体区域的热惯量：k 竺“2h j +Dp x ( eff x 丿ii(bP h(L-Q)p h)+2 uh )= dtf fs sx/ i fiT du +Sh +(L _0)Sh ik dxfsk其中： 匸”,Dp+ xj j Dti j人_彳=流体的焓 h_s=固体介质的焓 f=介质的多孔性 k_eff=介质的有效热传导系数SAh_f=流体焓的源项SAh_s=固体焓的源项 6、多孔介质的有效传导率多孔区域的有效热传导率k_eff是由流体的热传导率和固体的热传导率的体积平均值计算得 到：k = k +（L -keff fs其中：f=介质的多孔性k_f=流体

7、状态热传导率（包括湍流的贡献k_t）k_s=固体介质热传导率如果得不到简单的体积平均，可能是因为介质几何外形的影响。有效传导率可以用自定义函 数来计算。然而，在所有的算例中，有效传导率被看成介质的各向同性性质。7、多孔介质中的湍流处理在多孔介质中，默认的情况下 FLUENT 会解湍流量的标准守恒方程。因此，在这种默 认的方法中，介质中的湍流被这样处理：固体介质对湍流的生成和耗散速度没有影响。如果 介质的渗透性足够大，而且介质的几何尺度和湍流涡的尺度没有相互作用，这样的假设是合 情合理的。但是在其它的一些例子中，你会压制了介质中湍流的影响。如果你使用 k-e 模型或者 Spalart-Allma

8、ras 模型，你如果设定湍流对粘性的贡献 m_t 为零， 你可能会压制了湍流对介质的影响。当你选择这一选项时，FLUENT会将入口湍流的性质 传输到介质中，但是它对流动混合和动量的影响被忽略了。除此之外，在介质中湍流的生成 也被设定为零。要实现这一解策略，请在流体面板中打开层流选项 。激活这个选项就意味 着多孔介质中的 m_t 为零，湍流的生成也为零。如果去掉该选项（默认）则意味着多孔介 质中的湍流会像大体积流体流动一样被计算。8、概述模拟多孔介质流动时，对于问题设定需要的附加输入如下：L. 定义多孔区域2. 确定流过多孔区域的流体材料3. 设定粘性系数（多孔介质动量方程3中的l/a_ij）以

9、及内部阻力系数（多孔介质动量方 程3中的C_2_ij）,并定义应用它们的方向矢量。幕率模型的系数也可以选择指定。4. 定义多孔介质包含的材料属性和多孔性5. 设定多孔区域的固体部分的体积热生成速度（或任何其它源项，如质量、动量）（此项 可选）。6. 如果合适的话，限制多孔区域的湍流粘性。7. 如果相关的话，指定旋转轴和/或区域运动。在定义粘性和内部阻力系数中描述了决定阻力系数和/或渗透性的方法。如果你使用多孔动 量源项的幕律近似，你需要输入多孔介质动量方程5中的C_0和C_1来取代阻力系数和流 动方向。在流体面板中（下图）你需要设定多孔介质的所有参数，该面板是从边界条件菜单中打开的 （详细内容

10、请参阅边界条件的设定一节）Figure 1:多孔区域的流体面板定义多孔区域正如定义边界条件概述中所提到的，多孔区域是作为特定类型的流体区域来模拟的。 这表明流体区域是多孔区域，请在流体面板中激活多孔区域选项。面板会自动扩展到多孔介 质输入状态。定义穿越多孔介质的流体在材料名字下拉菜单中选择适当的流体就可以定义通过多孔介质的流体了。如果你模 拟组分输运或者多相流，流体面板中就不会出现材料名字下拉菜单了。对于组分计算，所有 流体和/或多孔区域的混合材料就是你在组分模型面板中指定的材料。对于多相流模型，所 有流体和/或多孔区域的混合材料就是你在多相流模型面板中指定的材料。定义粘性和内部阻力系数粘性和

11、内部阻力系数以相同的方式定义。使用笛卡尔坐标系定义系数的基本方法是在 二维问题中定义一个方向矢量，在三维问题中定义两个方向矢量，然后在每个方向上指定粘 性和/或阻力系数。在二维问题中第二个方向没有明确定义，它是垂直于指定的方向矢量和z 向矢量所在的平面的。在三维问题中，第三个方向矢量是垂直于所指定的两个方向矢量所在 平面的。对于三维问题，第二个方向矢量必须垂直于第一个方向矢量。如果第二个方向矢量 指定失败，解算器会确保它们垂直而忽略在第一个方向上的第二个矢量的任何分量。所以你 应该确保第一个方向指定正确。在三维问题中也可能会使用圆锥（或圆柱）坐标系来定义系数，具体如下： 定义阻力系数的过程如下

12、：1. 定义方向矢量。 使用笛卡尔坐标系，简单指定方向1矢量，如果是三维问题，指定方向2矢量。每 一个方向都应该是从（0,0）或者（0,0,0）到指定的（X,Y）或（X,Y,Z）矢量。（如果方向不正 确请按上面的方法解决） 对于有些问题，多孔介质的主轴和区域的坐标轴不在一条直线上，你不必知道多孔 介质先前的方向矢量。在这种情况下，三维中的平面工具或者二维中的线工具可以 帮你确定这些方向矢量。1. 捕捉Snap平面工具（或者线工具）到多孔区域的边界。（请遵循使用面工具和线工具中的说明，它在已存在的表面上为工具初始化了位置）。2. 适当的旋转坐标轴直到它们和多孔介质区域成一条线。3. 当成一条线之

13、后，在流体面板中点击从平面工具更新或者从线工具更新按钮。FLUENT 会自动将方向 1 矢量指向为工具的红（三维）或绿（二维）箭头所 指的方向。 要使用圆锥坐标系（比方说环状、锥状顾虑单元），请遵循下面步骤（这一选项只 用于三维问题）：1. 打开圆锥选项2. 指定圆锥轴矢量和在锥轴上的点。圆锥轴矢量的方向将会是从（0,0,0）到指定的 （X,Y,Z）方向的矢量。FLUENT将会使用圆锥轴上的点将阻力转换到笛卡尔坐 标系。3. 设定锥半角（锥轴和锥表面之间的角度，如下图），使用柱坐标系，锥半角为0.Figure 1:锥半角对于有些问题，锥形过滤单元的主轴和区域的坐标轴不在一条直线上，你不必知道

14、锥轴先前的方向矢量以及锥轴上的点。在这种情况下，三维中的平面工具或者二维 中的线工具可以帮你确定这些方向矢量。一种方法如下：1. 在点击捕捉到区域按钮之前，你可以在下拉菜单中选择垂直于锥轴矢量的轴过 滤单元的边界区域。2. 点击捕捉到区域按钮，FLUENT会自动将平面工具捕捉到边界。它也会设定锥 轴矢量和锥轴上的点（需注意的是你还要自己设定锥半角）。 另一种方法为：1. 捕捉Snap平面工具到多孔区域的边界。（请遵循使用面工具和线工具中的说明，它在已存在的表面上为工具初始化了位置）。2. 旋转和平移工具坐标轴，直到工具的红箭头指向锥的轴向。工具的起点在轴上。3. 当轴和工具的起点成一条线时，在

15、流体面板中点击从平面工具更新按钮。FLUENT会自动设定轴向矢量以及在轴上的点（注意：你还是要自己设定锥 的半角）。2. 在粘性阻力中指定每个方向的粘性阻力系数1/a,在内部阻力中指定每一个方向上的内 部阻力系数C_2 （你可能需要将滚动条向下滚动来查看这些输入）。如果你使用锥指定 方法，方向1为锥轴方向，方向2为垂直于锥表面（对于圆柱就是径向）方向，方向3 圆周（q）方向。在三维问题中可能有三种可能的系数，在二维问题中有两种： 在各向同性算例中，所有方向上的阻力系数都是相等的（如海绵。在各向同性算 例中你必须将每个方向上的阻力系数设定为相等。 在三维问题中只有两个方向上的系数相等，第三个方向上的阻力系数和前两个不 等，或者在二维问题中两个方向上的系数不等，你必须准确的指定每一个方向上的 系数。例如，如果你得多孔区域是由具有小洞的细管组成，细管平行于流动方向， 流动会很容易的通过细管，但是流动在其它两个方向上（通过小洞）会很小。如果 你有一个平的盘子垂直于流动方向，流动根本就不会穿过它而只在其它两个方向 上。 在三维问题中还有一种可能就是三