精华内容
下载资源
问答
  • 研究了一类带有点源和非局部边界条件的多孔介质方程非负解的爆破现象.运用经典的上、下解方法,并通过比较原理,在不同的假设下,构造出了方程合适的上解或下解,得出了该方程解的整体存在性和在有限时间爆破的充分条件.
  • 主要研究层状多孔介质中渗流问题在平衡状态下各种边值问题的平衡解的存在性,揭示了层状介质渗流问题解的一些独特的性质.
  • Washburn方程在分析毛细管流现象中有很好的代表性,其改进方程在预测多孔介质毛细上升高度前期阶段具有较高的吻合度.本文以不同粒径范围的洗砂作为研究对象,用该模型对不同实验条件下的数据进行优化拟合,得出该...
  • 针对多孔介质建模的复杂性,从随机动力学的朗之万方程出发,引入概率的影响,建立构造随机多孔介质的微观生长模型.通过不同的初始固相分布,该模型可以演化出各种不同性质的多孔介质微观结构.由于考虑了颗粒之间的作用力...
  • 利用微扰链-统计缔合流体理论状态方程(PC-SAFT)结合van der Waals-Platteuw模型和毛细管Kelvin模型,建立了用于多孔介质水合物体系的相平衡预测模型。在此模型基础上,针对甲烷水合物和CO2水合物对界面张力作了进一步...
  • 讨论了多孔介质的分形结构和相关的分形维数,利用能量方程,导出了分形维数为D的有限尺度多孔介质中的广义热传导方程。在此基础上,假定热量在多孔介质中的传导路线也是一种分形结构,提出了一个简化的多孔介质并联...
  • SIMPLE算法求解多孔介质的一维流动控制方程 问题介绍 求解思想 压力修正方法的基本思想 两个关键问题 求解步骤及说明 疑惑: a e a_e ae​表示动量方程离散系数?怎么求解? 本案例求解分析过程 本案例计算Python...

    问题介绍

    设流经某多孔介质的一维流动控制方程为:
    C ∣ u ∣ u + d p / d x = 0 C|u|u+dp/dx=0 Cuu+dp/dx=0 d ( u F ) / d x = 0 d(uF)/dx=0 d(uF)/dx=0
    其中,系数 C C C与空间位置有关, F F F为流道的有效截面积。
    对于图1所示的均匀网格,假定:
    C B = 0.25 C_B=0.25 CB=0.25 C C = 0.2 C_C=0.2 CC=0.2 F B = 5 F_B=5 FB=5 F C = 4 F_C=4 FC=4 p 1 = 200 p_1=200 p1=200 p 3 = 38 p_3=38 p3=38 Δ x = 2 \Delta{x}=2 Δx=2
    以上各单位都是协调的(也就是说不用管单位,直接带入数值计算就得了),试采用 S I M P L E SIMPLE SIMPLE算法确定 p 2 、 u B 、 u C p_2、u_B、u_C p2uBuC的值?
    在这里插入图片描述

    之前对于 ∣ u ∣ u |u|u uu的写法很不理解,20210803突然想到,这样的写法不就体现出方向性了吗?如果写成 u 2 u^2 u2就不能体现方向性了。所以这样的写法本人猜测其核心是为了体现方向性

    求解思想

    首先,从图1我们可以看出,本算例采用的是交错网格(staggered grid):速度 u u u和压力 p p p分别存储于两套不同的网格系统;速度 u u u存储于压力 p p p控制容积的东、西界面上。
    其次,我们主要是体会 S I M P L E SIMPLE SIMPLE算法的精髓:压力修正+预估校正。
    最后,在程序中体现迭代校正。自拟初始估计值。

    压力修正方法的基本思想

    1、假设一个压力场,记为 p ∗ p^{*} p;
    2、利用 p ∗ p^{*} p,求解动量离散方程,得到相应的速度 u ∗ u^{*} u;
    3、利用质量守恒方程来改进压力场,要求与改进后的压力场相对应的速度场能满足连续方程;
    4、用 p ′ p^{'} p表示压力修正值, u ′ u^{'} u表示速度修正值,以 p ∗ + p ′ p^{*}+p^{'} p+p u ∗ + u ′ u^{*}+u^{'} u+u作为本层次的解,并据此开始下一层次的迭代计算。

    两个关键问题

    1、如何获得压力修正值 p ′ p^{'} p,使( p ∗ + p ′ p^{*}+p^{'} p+p)满足连续性方程?
    2、获得 p ′ p^{'} p后,如何确定 u ′ u^{'} u?

    求解步骤及说明

    求解步骤与压力修正方法的基本思想不同,一个是先假定速度,一个是先假定压力,哪个是对的?
    答:没有矛盾。先假定速度是为了求动量离散方程的系数及源项

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    疑惑: a e a_e ae表示动量方程离散系数?怎么求解?

    在这里插入图片描述

    本案例求解分析过程

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

    本案例计算Python代码

    # 已知参数
    CB = 0.25
    CC = 0.2
    FB = 5
    FC = 4
    p1 = 200
    p3 = 38
    dx = 2
    
    # B、C处的速度试算值和初始压力假定
    uB_0 = 8
    uC_0 = 10
    p2_star = 100
    
    # 设置统计迭代次数变量
    count = 1
    
    # 迭代开始
    while 1:
        # 求速度预估值
        uB_star = 0.5 * uB_0 - (p2_star - p1) / (2 * CB * uB_0 * dx)
        uC_star = 0.5 * uC_0 - (p3 - p2_star) / (2 * CC * uC_0 * dx)
    
        # 求压力修正值
        p2_correction = 2 * dx * (uB_star * FB - uC_star * FC) / ((FB / (CB * uB_0)) + (FC / (CC * uC_0)))
    
        # 求速度修正值
        uB_correction = - p2_correction / (2 * CB * uB_0 * dx)
        uC_correction = p2_correction / (2 * CC * uC_0 * dx)
    
        # 迭代结果
        p2 = p2_star + p2_correction
        uB = uB_star + uB_correction
        uC = uC_star + uC_correction
        
        # 打印第count次的迭代结果(第count次的真值,可能不是最终真正的真值)
        print('第',count,'次迭代结果:')
        print('p2 = ',p2) 
        print('uB = ',uB)
        print('uC = ',uC)
        print('\n')
        # 退出循环判断
       # if abs(p2_correction) < 0.000001 and abs(uB_correction) < 0.000001 and abs(uC_correction) <0.000001:#退出循环条件
        if abs(p2 - p2_star) < 0.000001 and abs(uB - uB_0) < 0.000001 and abs(uC - uC_0) <0.000001:
            # 条件判断中的p2_star、uB_0、uC_0是上一迭代层的真值(从下面else中赋值语句可以看出)
            # 条件判断中的p2、uB、uC是最新迭代层(本迭代层)的真值
            # 以上一迭代层的真值与本迭代层的真值的差值小于10的负6次方作为收敛条件(收敛后每一迭代层的真值理应相等)
            break
        else:
            p2_star = p2
            uB_0 = uB
            uC_0 = uC
            count += 1
            
    print('只需迭代',count,'次使得残差值都小于10的负6次方')
    print('最终求得的结果为:')
    print('p2 = ',p2) 
    print('uB = ',uB)
    print('uC = ',uC)
      
    
    

    计算结果如下:
    在这里插入图片描述

    展开全文
  • 采用宏、细观相连接的损伤模型,描述了多孔介质的骨架刚度和流体通道演变过程,在连续介质力学守恒定律的基础上,推出了非饱和土的动力方程.进而,可方便地退化为饱和多孔介质的动力方程和一维固结方程.图3,参11.
  • 针对多孔介质建模的复杂性,从随机动力学的朗之万方程出发,引入概率的影响,建立构造随机多孔介质的微观生长模型.通过不同的初始固相分布,该模型可以演化出各种不同性质的多孔介质微观结构.由于考虑了颗粒之间的...
  • 针对岩体是一种可变形的多孔介质,...在研究渗流的过程中,考虑岩体蠕变的影响,基于岩体流变力学和渗流力学理论,应用对应性原理,建立了含有蠕变参数的动态应力平衡方程,推导了蠕变作用下多孔介质渗流的固流耦合数学模型。
  • 多孔介质理论的基础上,基于非线性达西定律并假设水分迁移过程为单向、可逆及水分迁移过程中无溶质迁移,推导得出了无相变以及考虑相变的水分迁移方程;引入土体传热方程和土骨架质量密度变化方程,过程中考虑了冻土中...
  • 基于多孔介质微观孔隙结构的分形特征、毛管模型和Poiseuille方程建立了计算分形多孔介质宏观物性参数渗透率与孔隙度的理论模型,给出了二者之间新的理论关系式。分形多孔介质宏观物性参数及其关系式是多孔介质微观...
  • 基于分形方法的多孔介质有效应力模型研究,闫铁,李玮,泰尔扎吉(Terzaghi)的有效应力方程是描述多孔介质应力关系的经典模型,后来的学者相继开展了很多研究,建立了相关模型,但基本都�
  • 其次,在Gurson理论已有工作的基础上,推出了一个规则多孔介质基体服从 Drucker-Prager屈服准则的近似屈服面方程,进而得到了相关联的弹塑本构模型.最后,以裂纹密度为细观损伤变量,修正裂纹密度细观损伤对屈服面...
  • 应用连续介质力学的基本原理,针对流体饱和两相多孔介质的特点,建立起增量形式的两相多孔介质弹塑性波动方程组,以实现对两相多孔介质弹塑性动力反应的描述.运用伽辽金方法对该波动方程组进行空间离散,得到两相...
  • Kozeny-Carman (KC) 方程多孔介质渗流领域最著名的半经验公式,长期以来,KC方程及其推广形式被广泛用于估算多孔介质的渗透率。但是,方程中的KC常数是一个没有确切物理意义的经验常数,且被证明并非一个常数值。...
  • 写出了多孔介质内不可压缩流体流动的基本方程组;利用流线迭代法对充满多孔介质的变截面圆管内的流体的流动进行数值计算,并通过算例对充满多孔介质管道中的流动情况进行分析.
  • 多孔介质条件 多孔介质模型可以应用于很多问题,如通过充满介质的流动、通过过滤纸、...当你使用这一模型时,你就定义了一个具有多孔介质的单元区域,而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定...

    转载自:http://zhengjun0228.blog.163.com/blog/static/71377014200971895419613/

    多孔介质条件

    多孔介质模型可以应用于很多问题,如通过充满介质的流动、通过过滤纸、穿孔圆盘、流量分配器以及管道堆的流动。当你使用这一模型时,你就定义了一个具有多孔介质的单元区域,而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定。通过介质的热传导问题也可以得到描述,它服从介质和流体流动之间的热平衡假设,具体内容可以参考多孔介质中能量方程的处理一节。

       

    多孔介质的一维化简模型,被称为多孔跳跃,可用于模拟具有已知速度/压降特征的薄膜。多孔跳跃模型应用于表面区域而不是单元区域,并且在尽可能的情况下被使用(而不是完全的多孔介质模型),这是因为它具有更好的鲁棒性,并具有更好的收敛性。详细内容请参阅多孔跳跃边界条件。

       

    多孔介质模型的限制

       

    如下面各节所述,多孔介质模型结合模型区域所具有的阻力的经验公式被定义为"多孔"。事实上多孔介质不过是在动量方程中具有了附加的动量损失而已。因此,下面模型的限制就可以很容易的理解了。

      流体通过介质时不会加速,因为事实上出现的体积的阻塞并没有在模型中出现。这对于过渡流是有很大的影响的,因为它意味着FLUENT不会正确的描述通过介质的过渡时间。

      多孔介质对于湍流的影响只是近似的。详细内容可以参阅湍流多孔介质的处理一节。

       

    多孔介质的动量方程

       

    多孔介质的动量方程具有附加的动量源项。源项由两部分组成,一部分是粘性损失项(Darcy),另一个是内部损失项:

    其中S_ii(x, y, or z)动量源项,DC是规定的矩阵。在多孔介质单元中,动量损失对于压力梯度有贡献,压降和流体速度(或速度方阵)成比例。

    对于简单的均匀多孔介质:

    其中a是渗透性,C_2时内部阻力因子,简单的指定DC分别为对角阵1/a C_2其它项为零。

    FLUENT还允许模拟的源项为速度的幂率:

    其中C_0C_1为自定义经验系数。

    注意:在幂律模型中,压降是各向同性的,C_0的单位为国际标准单位。

       

       

    多孔介质的Darcy定律

       

    通过多孔介质的层流流动中,压降和速度成比例,常数C_2可以考虑为零。忽略对流加速以及扩散,多孔介质模型简化为Darcy定律:

    在多孔介质区域三个坐标方向的压降为:

    其中 为多孔介质动量方程1中矩阵D的元素vj为三个方向上的分速度,D n_x D n_y、以及D n_z为三个方向上的介质厚度。

    在这里介质厚度其实就是模型区域内的多孔区域的厚度。因此如果模型的厚度和实际厚度不同,你必须调节1/a_ij的输入。.

       

    多孔介质的内部损失

    在高速流动中,多孔介质动量方程1中的常数C_2提供了多孔介质内部损失的矫正。这一常数可以看成沿着流动方向每一单位长度的损失系数,因此允许压降指定为动压头的函数。

       

    如果你模拟的是穿孔板或者管道堆,有时你可以消除渗透项而只是用内部损失项,从而得到下面的多孔介质简化方程:

    写成坐标形式为:

       

       

    多孔介质中能量方程的处理

       

    对于多孔介质流动,FLUENT仍然解标准能量输运方程,只是修改了传导流量和过度项。在多孔介质中,传导流量使用有效传导系数,过渡项包括了介质固体区域的热惯量:

    其中:

    h_f=流体的焓

    h_s=固体介质的焓

    f=介质的多孔性

    k_eff=介质的有效热传导系数

    S^h_f=流体焓的源项

    S^h_s=固体焓的源项

       

    多孔介质的有效传导率

       

    多孔区域的有效热传导率k_eff是由流体的热传导率和固体的热传导率的体积平均值计算得到:

    其中:

    f=介质的多孔性

    k_f=流体状态热传导率(包括湍流的贡献k_t

    k_s=固体介质热传导率

    如果得不到简单的体积平均,可能是因为介质几何外形的影响。有效传导率可以用自定义函数来计算。然而,在所有的算例中,有效传导率被看成介质的各向同性性质。

       

    多孔介质中的湍流处理

       

    在多孔介质中,默认的情况下FLUENT会解湍流量的标准守恒防城。因此,在这种默认的方法中,介质中的湍流被这样处理:固体介质对湍流的生成和耗散速度没有影响。如果介质的渗透性足够大,而且介质的几何尺度和湍流涡的尺度没有相互作用,这样的假设是合情合理的。但是在其它的一些例子中,你会压制了介质中湍流的影响。

    如果你使用k-e模型或者Spalart-Allmaras模型,你如果设定湍流对粘性的贡献m_t为零,你可能会压制了湍流对介质的影响。当你选择这一选项时,FLUENT会将入口湍流的性质传输到介质中,但是它对流动混合和动量的影响被忽略了。除此之外,在介质中湍流的生成也被设定为零。要实现这一解策略,请在流体面板中打开层流选项 。激活这个选项就意味着多孔介质中的m_t为零,湍流的生成也为零。如果去掉该选项(默认)则意味着多孔介质中的湍流会像大体积流体流动一样被计算。。

       

    概述

       

    模拟多孔介质流动时,对于问题设定需要的附加输入如下:

    1.         定义多孔区域

    2.         确定流过多孔区域的流体材料

    3.         设定粘性系数(多孔介质动量方程3中的1/a_ij)以及内部阻力系数(多孔介质动量方程3中的C_2_ij),并定义应用它们的方向矢量。幂率模型的系数也可以选择指定。

    4.         定义多孔介质包含的材料属性和多孔性

    5.         设定多孔区域的固体部分的体积热生成速度(或任何其它源项,如质量、动量)(此项可选)。

    6.         如果合适的话,限制多孔区域的湍流粘性。

    7.         如果相关的话,指定旋转轴和/或区域运动。

    在定义粘性和内部阻力系数中描述了决定阻力系数和/或渗透性的方法。如果你使用多孔动量源项的幂律近似,你需要输入多孔介质动量方程5中的C_0C_1来取代阻力系数和流动方向。

    在流体面板中(下图)你需要设定多孔介质的所有参数,该面板是从边界条件菜单中打开的(详细内容请参阅边界条件的设定一节)

    Figure 1:多孔区域的流体面板

       

    定义多孔区域

       

                正如定义边界条件概述中所提到的,多孔区域是作为特定类型的流体区域来模拟的。亚表明流体区域是多孔区域,请在流体面板中激活多孔区域选项。面板会自动扩展到多孔介质输入状态。

       

    定义穿越多孔介质的流体

       

                在材料名字下拉菜单中选择适当的流体就可以定义通过多孔介质的流体了。如果你模拟组分输运或者多相流,流体面板中就不会出现材料名字下拉菜单了。对于组分计算,所有流体和/或多孔区域的混合材料就是你在组分模型面板中指定的材料。对于多相流模型,所有流体和/或多孔区域的混合材料就是你在多相流模型面板中指定的材料。

       

    定义粘性和内部阻力系数

       

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

       

    在三维问题中也可能会使用圆锥(或圆柱)坐标系来定义系数,具体如下:

    定义阻力系数的过程如下:

    1.         定义方向矢量。

      使用笛卡尔坐标系,简单指定方向1矢量,如果是三维问题,指定方向2矢量。每一个方向都应该是从(0,0)或者(0,0,0)到指定的(X,Y)(X,Y,Z)矢量。(如果方向不正确请按上面的方法解决)

      对于有些问题,多孔介质的主轴和区域的坐标轴不在一条直线上,你不必知道多孔介质先前的方向矢量。在这种情况下,三维中的平面工具或者二维中的线工具可以帮你确定这些方向矢量。

    1.      捕捉"Snap"平面工具(或者线工具)到多孔区域的边界。(请遵循使用面工具和线工具中的说明,它在已存在的表面上为工具初始化了位置)。

    2.      适当的旋转坐标轴直到它们和多孔介质区域成一条线。

    3.      当成一条线之后,在流体面板中点击从平面工具更新或者从线工具更新按钮。FLUENT会自动将方向1矢量指向为工具的红(三维)或绿(二维)箭头所指的方向。

      要使用圆锥坐标系(比方说环状、锥状顾虑单元),请遵循下面步骤(这一选项只用于三维问题):

    1.         打开圆锥选项

    2.         指定圆锥轴矢量和在锥轴上的点。圆锥轴矢量的方向将会是从(0,0,0)到指定的(X,Y,Z)方向的矢量。FLUENT将会使用圆锥轴上的点将阻力转换到笛卡尔坐标系。

    3.         设定锥半角(锥轴和锥表面之间的角度,如下图),使用柱坐标系,锥半角为0.

    Figure 1:锥半角

      对于有些问题,锥形过滤单元的主轴和区域的坐标轴不在一条直线上,你不必知道锥轴先前的方向矢量以及锥轴上的点。在这种情况下,三维中的平面工具或者二维中的线工具可以帮你确定这些方向矢量。一种方法如下:

    1.         在点击捕捉到区域按钮之前,你可以在下拉菜单中选择垂直于锥轴矢量的轴过滤单元的边界区域。

    2.         点击捕捉到区域按钮,FLUENT会自动将平面工具捕捉到边界。它也会设定锥轴矢量和锥轴上的点(需注意的是你还要自己设定锥半角)。

      另一种方法为:

    1.      捕捉"Snap"平面工具到多孔区域的边界。(请遵循使用面工具和线工具中的说明,它在已存在的表面上为工具初始化了位置)。

    2.      旋转和平移工具坐标轴,直到工具的红箭头指向锥的轴向。工具的起点在轴上。

    3.      当轴和工具的起点成一条线时,在流体面板中点击从平面工具更新按钮。FLUENT会自动设定轴向矢量以及在轴上的点(注意:你还是要自己设定锥的半角)。

    2.         在粘性阻力中指定每个方向的粘性阻力系数1/a,在内部阻力中指定每一个方向上的内部阻力系数C_2(你可能需要将滚动条向下滚动来查看这些输入)。如果你使用锥指定方法,方向1为锥轴方向,方向2为垂直于锥表面(对于圆柱就是径向)方向,方向3圆周(q)方向。

    在三维问题中可能有三种可能的系数,在二维问题中有两种:

      在各向同性算例中,所有方向上的阻力系数都是相等的(如海绵)。在各向同性算例中你必须将每个方向上的阻力系数设定为相等。

      在三维问题中只有两个方向上的系数相等,第三个方向上的阻力系数和前两个不等,或者在二维问题中两个方向上的系数不等,你必须准确的指定每一个方向上的系数。例如,如果你得多孔区域是由具有小洞的细管组成,细管平行于流动方向,流动会很容易的通过细管,但是流动在其它两个方向上(通过小洞)会很小。如果你有一个平的盘子垂直于流动方向,流动根本就不会穿过它而只在其它两个方向上。

      在三维问题中还有一种可能就是三个系数各不相同。例如,如果多孔区域是由不规则间隔的物体(如针脚)组成的平面,那么阻碍物之间的流动在每个方向上都不同。此时你就需要在每个方向上指定不同的系数(请注意指定各向同性系数时,多孔介质的解策略的注解)。

       

    推导粘性和内部损失系数的方法在定义粘性和内部阻力系数一节中介绍。

       

    当你使用多孔介质模型时,你必须记住FLUENT中的多孔单元是100%打开的,而且你所指定1/a_ij/C_2_ij的值必须是基于这个假设的。然而,假如你知道通过真实装置压降和速度之间的的变化,它只是部分地对流动开放。下面的练习会告诉你如何对FLUENT模型计算适当的C_2值。

       

    假定穿孔圆盘只有25%对流动开放。已知通过圆盘的压降为0.5。在圆盘内真实流体速度基础上,即通过%开放区域的的基础上,损失系数由下式定义的损失系数K_L0.5

    要计算适当的C_2值,请注意在FLUENT模型中:

    1.         通过穿孔圆盘的速度假定圆盘为100%开放的。

    2.         损失系数必须转化为多孔区域每个单位长度的动压头损失。

    对于第一条,第一步是计算并调节损失因子K_L',它应该是在100%开放区域的速度基础上的:

    或者注意对于相同的流速,v_25% open = 4 v_100% open,

    调节之后的损失系数为8。对于第二条,你必须将它转换为穿孔圆盘每个单位厚度的损失系数。假定圆盘的厚度为1.0 mm。内部损失系数为(国际标准单位):

    注意,对于各向异性介质,这些信息必须分别从每一个坐标方向上计算。

       

    第二个例子,考虑模拟充满介质的流动。在湍流流动中,充满介质的流动用渗透性和内部损失系数来模拟。推导适当常数的方法包括了Ergun方程[49]的使用,对于在很大范围雷诺数内和许多类型的充满形式,有一个半经验的关系式:

    当模拟充满介质的层流流动时,上面方程中的第二项可能是个小量,从而得到Blake-Kozeny方程[49]

    在这些方程中,m是粘性,D_p是平均粒子直径,e空间所占的分数(即空间的体积除以总体积)。比较多孔介质中Darcy定律的方程1和内部损失系数为9的方程1,则每一方向上的渗透性和内部损失系数定义为:

    第三个例子我们会考虑Van Winkle等人[146][121]的方程,并表明如何通过具有方孔圆盘的多孔介质输入来计算压力损失。

       

    作者所声明的应用在通过在等边三角形上的方洞圆盘的湍流中的表达式为:

    其中:

    m(dot)=通过圆盘的质量流速

    A_f=剩下的面积或者洞的总面积

    A_p=圆盘的面积(固体和洞)

    C=对于不同D/t的不同雷诺数范围被列成不同的表的系数

    D/t=洞的直径和圆盘厚度的比例

    对于t/D > 1.6Re > 4000,系数C近似为0.98,其中雷诺数是基于洞的直径与速度的

    使用下式整理方程17

    除以圆盘的厚度D x = t有:

    其中v是表面速度而不是洞内的速度。与多孔介质内部损失系数中的方程1比较可以看出,对于垂直于圆盘方向,常数C_2可由下式计算:

    考虑通过由随机方向的纤维或者玻璃材料组成的垫子或者过滤器的层流。对于可以二选一的方程Blake-Kozeny(方程11),我们可能会选择将实验数据列成表。很多类型的纤维都由这一类相关的数据[70]

       

    固体体积分数f  玻璃丝织品的无量纲渗透性Q

    0.262         0.25

    0.258         0.26

    0.221         0.40

    0.218         0.41

    0.172         0.80

    其中Q = a为纤维直径。使用多孔介质的Darcy定律中的方程1可以很容易从给定的纤维直径和体积分数种计算出 

       

    使用幂律模型

       

                对于多孔介质动量源项(多孔介质动量方程中的方程5),如果你使用幂律模型近似,你只要在流体面板的幂律模型中输入系数C_0C_1就可以了。如果C_0C_1为非零值,解算器会忽略面板中除了多孔介质幂律模型之外的所有输入。

       

    定义热传导

       

                如果你选择在多孔介质中模拟热传导,你必须指定多孔介质中的材料以及多孔性。要定义多孔介质的材料,向下拉流体面板中阻力输入下面的滚动条,然后在多孔热传导的固体材料下拉列表中选中适当的固体。

                然后在多孔热传导下设定多孔性。多孔性f是多孔介质中流体的体积分数(即介质的开放体积分数)。多孔性用于介质中的热传导预测,处理方法请参阅多孔介质能量方程的处理一节。它还对介质中的反应源项和体力的计算有影响。这个源项和介质中流体的体积成比例。如果你想要模拟完全开放的介质(固体介质没有影响),你应该设定多孔性为1.0。当多孔性为1.0时,介质的固体部分对于热传导和(或)热源项/反应源项没有影响。注意:多孔性永远不会影响介质中的流体速度,这已经在多孔介质的动量方程一节中介绍了。不管你将多孔性设定为何值,,FLUENT所预测的速度都是介质中的表面速度。

       

    定义源项

       

                如果你想在多孔流动的能量方程中包括热的影响,请激活源项选项并设定非零的能量源项。FLUENT会计算多孔区域所生成的能量,该能量为能量源项值乘以组成多孔区域的单元所有体积值。你也可以定义质量、动量、湍流、组分或者其它标量的源项,详细内容请参阅、质量、动量、能量和其它源项的定义。

       

    在多孔区域内压制湍流源项

       

                如多孔介质的湍流处理中所讨论的,湍流在多孔介质中的计算和大量(bulk)流体流动是一样的。如果你使用k-e模型或者Spalart-Allmaras模型,你想要压制湍流在多孔区域的影响可以打开流体区域面板中的层流区域选项(从而使得多孔区域的湍流生成为零)。

       

    指定旋转轴并定义区域运动

                 

                旋转轴和区域运动的输入和标准流体区域的输入是相同的,详细情况可以参阅流体区域的输入一节。

       

    多孔介质的解策略

                一般说来,在模拟多孔介质时,你可以使用标准的解算步骤以及解参数的设置。然而你会发现如果多孔区域在流动方向上压降相当大(比如:渗透性a很低或者内部因子C_2很大)的话,解的收敛速度就会变慢。这就表明由于动量源项中出现了多孔介质的压降(方程的矩阵不再是对角占优了),收敛性问题就出现了。解决多孔介质区域收敛性差最好的补救办法就是对于通过介质的流向压降有一个很好初始预测。猜测的办法之一就是,在介质流体单元的上游或者下游补偿一个压力值,详细内容请参阅所选单元的补偿值一节。必须记住的是,当补偿压力时,你所输入的压力可以定义为解算器所使用的gauge压力(即在操作条件面板中定义的相对于操作压力的压力)。

                另一个处理收敛性差的方法是临时取消多孔介质模型(在流体面板中关闭多孔区域)然后获取一个不受多孔区域影响的初始流场。取消多孔区域后,FLUENT会将多孔区域处理为流体区域并按相应的流体区域来计算。一旦获取了初始解,或者计算很容易收敛,你就可以激活多孔模型继续计算包含多孔区域的流场(对于大阻力多孔介质不推荐使用该方法)。

                对于高度各向异性的多孔介质,有时会造成收敛性的麻烦。对于这些问题你可以将多孔介质的各向异性系数(1/a_ijC_2_i,j)限制在二阶或者三阶的量级。即使在某一方向上介质的阻力为无穷大,你也不需要将它设定超过初始流动方向上的1000倍。

       

    多孔介质的后处理

       

                可以通过检查速度分量和压力值来确定多孔区域对于流场的影响。你可能对下列变量或函数的图形(XY图,等值线图或者矢量图)或者文档报告感兴趣:

      X,Y,Z速度(在速度类别中)

      静压(在压力类别中)

                这些变量会在后处理面板的变量选择下拉菜单制定类别中出现。

                需要注意的是多孔区域的热报告不影响固体介质的属性。所报告的多孔区域内的热容、传导率以及焓是流体的属性不包括固体介质的影响。

    转载于:https://www.cnblogs.com/liusuanyatong/p/11259785.html

    展开全文
  • 多孔介质内部大圆球绕流问题进行理论分析。首先根据基本假设写出简化的基本方程组,然后求出多孔介质内部大圆球绕流问题的解析解,最后对解析解进行分析和讨论。
  • 结果表明,Elovich模型和零级动力学方程分别能很好描述多孔介质中汽油和柴油挥发动力学曲线,且汽油在粗砂中的挥发速率系数较在亚黏土中大,而柴油在亚黏土中的挥发速率系数是在粗砂中的2.0~2.5倍。采用线性方程和...
  • 三维MHD模型在多孔介质中的整体吸引子问题,钟婷茹,,在本篇论文中,作者主要研究了三维MHD模型在多孔介质中的整体吸引子问题,详细验证了相应的半群有吸收集.通过方程组相应的半群具�
  • 文中以水为工质,对多孔介质中Forchheimer区和湍流区下阻力压降进行研究,并将其与经典Ergun型方程的计算值进行了比较和分析.结果表明,湍流区中Ergun型方程的值与实验值出现很大偏差,但二者之间有比较相似的总体...
  • 根据城市空间结构与多孔介质的相似性,在城市热岛效应的研究中引入了多孔介质理论,采用体积平均技术和雷诺时均方程法建立了适用于城市湍流流动与换热的多孔介质数学模型,并利用微尺度数值模拟计算结果的空间体积...
  • 文章目录一、多孔介质(porous media)二、使用步骤1.达西定律2.Darcy-Forchheimer3.OpenFOAM fvOptions4.paraview动画每帧显示时间(Annotate Time)5.参考内容总结 一、多孔介质(porous media) 在流体力学中,...


    一、多孔介质(porous media)

    在流体力学中,流体流经多孔介质是指流体流经多孔介质时的行为方式,如海绵或木材,或使用沙子或其他多孔材料过滤水。
    我们将展示如何设置多孔层,并使用OpenFOAM模拟流体通过该介质。

    二、使用步骤

    1.达西定律

    达西定律最初是由达西通过实验确定的,但后来通过均匀化方法从Navier-Stokes方程推导出来。热传导领域类似于傅里叶定律,电网络领域类似于欧姆定律,扩散理论类似于菲克定律。
    达西定律的一个应用是用于分析通过含水层的水流;达西定律和质量守恒方程简化了水文地质基本关系之一的地下水流动方程。

    2.Darcy-Forchheimer

    Darcy Forchheimer模型用于模拟多孔介质,允许我们在流体域内简单地添加一个孔隙带,而无需任何费用。

    在动量方程中,Sm作为源项
    在这里插入图片描述
    D和F在fvOption中修改
    在这里插入图片描述

    3.OpenFOAM fvOptions

    为了使用这个模型,你必须把fvOptions文件放到 0/constant 常量文件夹中,其中包括以下内容:

    达西公式
    0/constant/fvOptions

    porosity1
    {
        type            explicitPorositySource;//模拟多孔介质区域(显示)
        active          yes;
    
        explicitPorositySourceCoeffs
        {
            selectionMode   cellZone;//在toposet中配置
            cellZone        fixedWalls;//在toposet中配置
    
            type            DarcyForchheimer;
    
    	    f 4050; 
            d 0; //本案例不考虑粘性损失
            DarcyForchheimerCoeffs
            {
                d   ($d $d $d);
                f   ($f $f $f);
                           
            coordinateSystem
            {
                type    cartesian;
                origin  (0 0 0);
                coordinateRotation
                {
                    type    axesRotation;
                    e1  (1 0 0);
                    e2  (0 1 0);
                }
            }
        }
      }
    }
    

    多孔介质模型的位置
    0/system/topoSetDict

    actions
    (
       {
    	name	porousWalls;
            type	cellSet;//faceSet;
    	action	new;
    	source  boxToCell;
    	sourceInfo
    	{
    	   box (0 0 -1) (0.2 1 1); 
    	}
        } 
       //faceZone:leftAndRightWall
       {
    	name	fixedWalls;
    	type	cellZoneSet;//faceZoneSet;
    	action	new;
     	source  setToCellZone;
    	sourceInfo
    	{	
    	   //faceSet	  porousWalls;
    	   set	  porousWalls;
    	}
       } 
    
        
    );
    
    

    4.paraview动画每帧显示时间(Annotate Time)

    paraview动画每帧显示时间(Annotate Time)

    1. 点击Filters->Alphabetical->ProgrammableFilter,在Properties->Script中设置以下两段代码中任意一段,按照自己的需求选择并修改
    pdo =  self.GetOutput()
    sexaTime = vtk.vtkFloatArray()
    sexaTime.SetName("SexaTime")
    t = inputs[0].GetInformation().Get(vtk.vtkDataObject.DATA_TIME_STEP())
    #convertedTime = t+18.5 #compute your time here in python
    convertedTime = t//12+18+(t*5)%60/100
    sexaTime.InsertNextValue(convertedTime)
    pdo.GetFieldData().AddArray(sexaTime)
    
    #set as string
    pdo =  self.GetOutput()
    sexaTime = vtk.vtkStringArray()
    sexaTime.SetName("SexaTime")
    t = inputs[0].GetInformation().Get(vtk.vtkDataObject.DATA_TIME_STEP())
    #timeAsAString = f(t) #compute your time here in python in a string format
    s1=t//12+18
    s2=(t*5)%60
    timeAsAString = str(t//12+18)+':'+str((t*5)%60)
    sexaTime.InsertNextValue(timeAsAString)
    pdo.GetFieldData().AddArray(sexaTime)
    
    1. 在Pipeline Browser中选中ProgrammableFilter然后点击Filters->Annotation->AnnotateGlobalData,然后在Properties中设置AnnotateGlobalData的显示格式
      结果如下所示:
      在这里插入图片描述

    5.参考内容

    Paraview设置动画显示时间(Annotate Time)格式:https://www.jianshu.com/p/404f183e7a77
    DarcyForchheimer-openfoamwiki:
    https://openfoamwiki.net/index.php/DarcyForchheimer
    fvOptions----运行时可选的物理操作:https://blog.csdn.net/hanbingchegu/article/details/107417054
    DarcyForchheimer-openfoamwiki:
    https://openfoamwiki.net/index.php/DarcyForchheimer
    (详细)Flow through a porous media: using OpenFOAM:https://www.hyperlyceum.com/product/flow-through-a-porous-media-using-openfoam/


    总结

    OpenFOAM : fvOptions, topoSetDict; Paraview : Annotate Time
    CFD小白的进阶之路。

    展开全文
  • 为研究流体在石墨多孔介质中的流动规律,通过水的单相饱和渗流实验,揭示了水在不同孔隙率的石墨多孔介质中的非Darcy渗流规律.实验结果表明:流动偏离达西定律,呈现非线性,可以用一元二次方程的形式拟合,而且...
  • 本文基于Luikov热湿耦合传递方程,提出了一种计算一维多孔介质瞬态热湿传递的理论分析方法.通过引入无量纲温度、湿度对Luikov方程进行拉普拉斯变换,采用传递函数方法(TFM)得到了一维瞬态热湿传递方程的解析解,该解析...
  • 采用随机漫步法生成了分形孔模型,根据分子运动论建立了能应用于分形孔中的气体扩散控制方程,在所建立的扩散模型基础上,研究了分形多孔介质中气体稳态扩散规律并得到了稳态扩散系数表达式。结果表明:扩散系数正比...
  • 方法通过建立水平土柱模型,应用一维非饱和砂土热量运移定解方程,在室内测取不同回灌温度条件下不同多孔介质的热弥散度.结果相同温度回灌条件下,中砂热弥散度最大,粗砂次之,细砂最小;同种多孔介质回灌条件下,...
  • 多孔介质中的流动、传热与化学反应-研究现状及展望 ,姜元勇,,多孔介质分布广泛,研究其中的流动、传热与化学反应问题具有重要的实际应用和科研价值。评述了多孔介质中流动、传热与化学反应的

空空如也

空空如也

1 2 3 4 5 ... 17
收藏数 323
精华内容 129
关键字:

多孔介质方程