精华内容
下载资源
问答
  • 2021-04-24 02:03:05

    【实例简介】

    好文献大家一定要顶,不要错过啊!确定多孔介质流动参数的格子Boltzmann方法

    4期

    钱吉裕等:确定多孔介质流动参数的格子 Boltzmann方法

    657

    3多孔介质的构造及边界条件

    30

    孔隙率为0.5058

    250

    孔隙率为04112

    本文中的多孔介质是用实际泡沫金属材料的电

    ▲孔隙率为04563

    20U

    镜照片(如图3(a))经过图像处理方法得到的,图像

    处理过程主要包括灰度拉伸、梯度锐化、噪声去除

    100

    阈值变换以及中值滤波等。图3(b)和3(c)是处理得

    到的两种不同孔隙率的多孔介质,图像中的每个像

    素对应一个格子,利用格子上的灰度值,可以判定每

    00000.020.00400060.008

    个格子点上的材料属性。根据位图中个像素上的灰

    速度/m/s

    度值,可统计得到孔隙率,有效孔径等结构参数。

    图5压力梯度与速度的关系

    本义的数值计算中,流体工质选用水,根据式5结论

    c=1,T=1.0;格子数150×300。边界条件设定

    本文釆用D2Q9不可压模型计算了多孔介质中

    如下:入口边界为un=Uhm,un=0;出口边界为的流场,针对多孔介质中的LB模型给出了计算多孔

    尸=Pm;上下边界取周期边界;流体通道与固体介质渗透率的一种简便方法,验证∫ Darcy定律。运

    骨架之间采用无滑移反弹边界条件

    用该方法,可以研究孔隙中介观层次的流动机制

    还可以计算一般各向同性或非同性多孔介质的宏观

    4数值结果及讨论

    渗透率。

    本文运用图象处理方法构造的多孔介质结构,

    图4是针对上述方法构造的多孔介质,采用LB能较好地体现实际多孔介质的儿何特性,该方法是

    方法计算得到的速度矢量图,流体从左边界流进,

    种构造二维多孔介质的快速有效的方法。

    从图中可以看出横向孔道中的流速普遍大于纵向孔

    道的流速,这是因为在该边界条件下,压力梯度主参考文獻

    要在横向

    林瑞泰多孔介质传热传质引论.北京:科学出版社,1995

    在层流条件下,不考虑重力的作用

    arc

    y定

    [2]郭照立,李青,郑楚光双扩散自然对流的格子 Boltzmann

    律可以写成阿

    模拟.计算物理,2002,19(6):483487

    ]许友生,刘慈群,俞慧丹多孔介质中两相驱离的格子 Boltz

    man模型新研究.应用数学和力学,2002,23(4):353-359

    yp

    (10)

    [4 Bruce J Palmer, David R. Rector Lattice Boltzmann Al-

    gorithm for Simulating Thermal Flow in Compressible

    对于非规则多孔介质,如果宏观上是各向同性材

    Fluids. J. Computational Physics, 2000, 161: 1-20

    料,选择足够大的计算区域就可以认为K是一个标

    5] Zhaoli Guo, Baochang Shi, Chuguang Zheng, A Couple

    Lattice BGK Model for the Boussinesq Equations. IntJ

    量,根据速度与压差的关系就可以计算出渗透率(如

    Numer. Meth. Fluids. 2002. 39:325-342

    图5),图中,孔隙率为04112,04563,0.5058的多刘雅仁,陈若航李华兵等二维对流扩散方程的格子Bo

    孔介质的渗透率分别为548×10-11mn2,666×10-11

    amn方法.物理学报,1999,48(10):1800-1803

    7]孔样言.高等渗流力学.合肥中国科学技术大学出版社,

    m2,8.19×10-1m2。

    1999

    图4多孔介质内的速度矢量图(右图是左图指定区域的局部放大图〕

    【实例截图】

    【核心代码】

    更多相关内容
  • 关于matlab在多孔介质中的数值模拟应用及介绍
  • 多孔介质相关求解问题,以及关于如何利用多孔介质理论去简化模型。
  • 介绍随机生长 四参数生成法, 构造黏土微观结构的基本原理和详细的生成步骤, 联合应用 Ma tlab 和 AutoCAD VBA 编制相应的 程序, 讨论构造的黏土微观结构对地下水渗流的影响.程序生成的微观结构可以直接导入有限元...
  • 在实验研究的基础上,应用流体力学分析软件Fluent 6.3对多孔介质孔道内流体流动进行了模拟。研究了进口气速、多孔介质长度、多孔介质孔隙率对孔道进出口压降的影响规律。研究结果表明,多孔介质长度较小时,压降与...
  • 我国褐煤资源丰富,但是由于其水分高,导致了褐煤燃烧发电效率低下,不宜...分析了多孔介质的干燥过程特点,阐述了目前国内外的多孔介质干燥机理,总结了相关的多孔介质干燥模型,为多孔介质褐煤热力脱水模型化提供理论支持.
  • 为了减小采空区渗透率分布对流场数值模拟精确程度的影响,在岩梁理论和“O”型圈理论研究的基础上,分析了采空区非均质多孔介质孔隙分布规律,建立了采空区孔隙率及渗透率三维分布数学模型;以姚桥煤矿7271综放工作...
  • 多孔介质中的流动、传热与化学反应-研究现状及展望 ,姜元勇,,多孔介质分布广泛,研究其中的流动、传热与化学反应问题具有重要的实际应用和科研价值。评述了多孔介质中流动、传热与化学反应的
  • 讲述COMSOL在多孔介质中的应用案例及功能,有助于了解comsol在此方面的应用实例。
  • 多孔介质渗透率的准确确定是多孔介质传热传质研究的一个重要组成部分,而传统的Taylor公式对于不均匀多孔介质的预测效果很不理想。为此,利用实际多孔介质剖面图,系统研究了颗粒粒径分布分形维数乙、多孔介质孔隙面积...
  • 多孔介质在全腔-肺吻合术计算流体力学模型中模拟肺血管阻力应用探讨,孙彦隽,刘金龙,【目的】 探讨多孔介质在全腔-肺吻合术(Total Cavo-pulmonary Connection,TCPC)的计算流体力学(Computational Fluid ...
  • Micro-CT技术在多孔介质两相流动实验研究中的应用,赵越超,陈俊霖,可视化及量化分析技术是开展多孔介质多相渗流实验研究的重要手段,为了探索微焦点X射线CT(Micro-CT)技术在多孔介质两相流动实验研�
  • 针对岩体是一种可变形的多孔介质,...在研究渗流的过程中,考虑岩体蠕变的影响,基于岩体流变力学和渗流力学理论,应用对应性原理,建立了含有蠕变参数的动态应力平衡方程,推导了蠕变作用下多孔介质渗流的固流耦合数学模型。
  • CFD多孔介质

    2022-02-02 21:33:42
    多孔介质

    多孔介质

    1. 应用场景

    在实际的工程计算中经常会遇到如:土壤,海绵,岩石,过滤器,滤网这一类的问题计算。

    对于这一类介质的计算来说,其几何空隙非常多,建立真实的几何非常麻烦,且会产生非常多的网格,工作量和计算量都非常大。因此,可以简化为多孔介质进行计算。

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

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

    2. 基本原理

    将流体域设置为多孔介质后,可以人为定义各个方向上的阻力系数,来代替多孔介质中的固体对流体的阻力。即添加一个与速度相关的动量汇,其表达公式为
    S i = − ( ∑ j = 1 3 D i j μ v j + ∑ j = 1 3 C i j 1 2 ρ ∣ v ∣ v j ) S_{i}=-\left(\sum_{j=1}^{3} D_{i j} \mu v_{j}+\sum_{j=1}^{3} C_{i j} \frac{1}{2} \rho|v| v_{j}\right) Si=(j=13Dijμvj+j=13Cij21ρvvj)
    其中, S i S_{i} Si 为第i(x,y,z)方向上的动量方程的源项;v为速度值;

    D 与C 为指定系数矩阵,D系数这一项代表粘性损失项,C系数这一项代表惯性损失项

    对于均匀的多孔介质来说,可改写为下面这个式子 S i = − ( D μ v i + C 2 1 2 ρ ∣ v ∣ v i ) S_{i}=-\left(D \mu v_{i}+ C_2 \frac{1}{2} \rho|v| v_{i}\right) Si=(Dμvi+C221ρvvi)

    其中,

    D为黏性阻力系数,又可以被改写为1/α,α为渗透率;

    C 2 C_2 C2为惯性阻力系数;v为速度

    动量汇 S i S_i Si作用于流体产生了压力梯度 D p D_p Dp,即 D p = − S i ∗ l D_p = -S_i * l Dp=Sil , l l l多孔介质域的厚度

    压降和流速的关系是一个一元二次方程的形式 Y = A x 2 + B x Y=Ax^{2} + Bx Y=Ax2+Bx

    我们在计算多孔介质区域时最重要的两个参数 D D D C 2 C_2 C2 可以通过拟合压降和流速的关系式来求得。

    例如有压降与速度关系式 △ p = 1398.8 u 2 + 1785.1 u \triangle p = 1398.8 u^{2} + 1785.1u p=1398.8u2+1785.1u

    于是有:

    1398.8 = C 2 1 2 ρ ∗ l 1398.8 = C_2\frac {1}{2}\rho*l 1398.8=C221ρl, 根据流体密度滤网(多孔介质)的厚度可以计算得到惯性阻力系数 C 2 C_2 C2

    1785.1 = D ∗ μ ∗ l 1785.1 = D*\mu*l 1785.1=Dμl, 根据流体的动力黏度滤网厚度可以计算得到黏性阻力系数 D D D,也就是 1 α \frac{1}{\alpha} α1

    3. 用CFD方法获取压降与速度的关系式

    建立滤网几何模型的思路

    下图是一个滤网的示意图,从图中可以看到小孔的布置方式为六边形布置,即每个小孔和其周围的六个小孔都是等距的,因此可以将整个滤网两侧的流域简化为计算单个小孔,并将周围的流体域取成六棱柱的形状,边长就是两个小孔之间的距离。

    这样就可以使用对称边界来将滤网前后的整个流体域简化为六棱柱型的小块流体域。

    img

    3D几何和网格

    img

    4. fluent设置

    csdn

    其中,relative velocity resistance formulation 勾选表示粘性阻力系数用的是相对速度

    alternative formulation 勾选的作用是,当各向异性差异十分明显时(比如10000:1),用以增加收敛性

    若流体在某主流方向(如x方向)流通,而其他两方向(y,z方向)不流通,则需要设置其他两个方向(y,z方向)的粘性阻力系数和惯性阻力系数比该主流方向(x方向)的阻力系数大三个数量级

    csdn多孔介质

    孔隙率为流体所占据的体积分数

    thermal model中热平衡指多孔介质基体温度约等于孔间流体温度,非热平衡指基体与孔间流体温度差距较大

    relative viscosity(相对粘度)用于计算effective viscosity(有效粘度, which is introduced to account for the effect of the porous medium on the diffusion term in the momentum equations), 计算公式如下
    μ e = μ r ∗ μ \mu_e = \mu_r*\mu μe=μrμ

    其中相对粘度 μ r \mu_r μr可由以下模型计算得到
    在这里插入图片描述

    其中
    γ \gamma γ是 porosity of the medium(the volume fraction of fluid within the porous region (that is, the open volume fraction of the medium)

    展开全文
  • 随着研究的不断深入,超声波探测技术不受压力温度限制、不改变流体成分、检测成本低、对人体无害等优点逐渐得到重视,并在多孔介质内相态研究领域逐渐得到了应用。在参考国内外相关文献资料的基础上,简述了超声波的...
  • 为准确描述石墨多孔介质的孔隙...压汞实验的结果表明,石墨试样体积分形维数在2.68~2.92之间,与Menger海绵的分形维数2.73较为接近,实验值与理论计算值基本吻合,可以应用分形理论描述石墨多孔介质的孔隙度与比表面积。
  • 通过试验验证,应用Fluent软件对多孔介质结构的微燃烧室进行数值模拟,分析孔隙率、氢氧体积混合比和流量比等因素对燃烧室性能的影响.模拟结果表明,采用多孔介质结构可以改善燃烧室的燃烧性能,燃烧室在多孔介质...
  • 地下的各种现象的特点是变形结构之间的相互作用,如裂缝和耦合的热、水力和机械过程。地下动力学的模拟可以提供有价值...关键词:热流体力学、断裂、断裂变形、多孔介质、多点有限体积、剪切膨胀、离散断裂矩阵、混合维
  • 二十、多孔介质模型案例

    千次阅读 2022-04-26 16:41:54
    Fluent自带多孔介质模型,对于多孔介质的模拟,不考虑流体在多孔介质内部的流动,只考虑多孔介质对于流动阻力及能量方程产生的影响。 2 模型描述 本例的模型如图所示。模型有一个进口和一个出口,入口气流流速...

    1  概念介绍

    多孔介质就是固体物质内部和表面有许多孔隙,如海绵等,由固体物质组成的骨架和由骨架分隔成大量密集成群的微小空隙所构成的物质。多孔介质内的流体以渗流方式运动。

    Fluent自带多孔介质模型,对于多孔介质的模拟,不考虑流体在多孔介质内部的流动,只考虑多孔介质对于流动阻力及能量方程产生的影响。

                                                  

    2  模型描述

    本例的模型如图所示。模型有一个进口和一个出口,入口气流流速为22.6m/s,中间经过多孔介质区域,出口压力为0,最后经OUTLET流出。

     

    3  导入网格

    使用Workbench打开工程文件,文件在本文末尾链接资源内。

    4  Scale网格尺寸

    Scale修改网格尺寸。如图所示。

    确保计算域尺寸是我们所需要的。

    5  设置求解器

    选择压力基(pressure-based)求解器,同时选择稳态模拟。

     

    6  设置计算模型

    本例空气进口流速较大为22.6m/s,因此湍流模型选择为Standard k-e湍流模型,标准壁面函数。

    Standard k-e湍流模型具有高稳定性、经济性及精确性,应用最为广泛,但对于旋流等问题精确性较差。

    不考虑能量方程,因此能量方程保持关闭

     

    7  材料设置

    在Materials下,双击air材料,弹出如下面板。

    单击Fluent Database,选择nitrogen(N2),单击copy,然后单击close,N2材料添加成功

     

    8 设置计算域

    单击信息树Cell Zone Conditions,出现两个计算域分别为fluid和substrate(想要出现多个计算域,需要在建模和画网格时进行一定的设置)

    双击Fluid,弹出流体域设置对话框,将Material Name改为nitrogen,即将工作流体改为N2而不再是空气,其他设置保持默认。

     

    接下来就是多孔介质的设置了,多孔介质设置只需要在计算域进行设置即可,我们进行较为详细的讲解

     

    双击substrate,弹出流体域设置对话框,勾选Porous Zone表明这部分计算域为多孔介质区域。

    将Material Name改为nitrogen,表示多孔介质中的流体也为N2,由于不考虑能量方向,因此不必指定固体域材料

    Laminar Zone勾选Laminar Zone表明多孔介质区域中抑制湍流的产生:

    如果不勾选,多孔介质区域的固体边界仿佛对湍流的产生和耗散没有影响,这个假定只有在渗透性很大并且多孔介质的几何尺度对湍流涡旋的几何尺度没有影响的情况下可行

    选择Porous Zone,出现许多设置,接下来我们一一说明

    Direction -1 Vector和Direction -1 Vector:

    多孔介质分为各向同性多孔介质,指的是多孔介质各个方向的阻力相同。各向异性指各个方向阻力不同,有的方向流体容易通过,有的方向流体很难通过。

     

    若多孔介质为各向同性,此设置无意义。但若为各向异性,则方向1矢量表示多孔介质的第一个主方向为x方向,方向2矢量表示多孔介质的第二个主方向为y方向,第三个方向与这两个方向垂直,不必指定。

     

    Relative Velocity Resistance Formulation:

    当坐标系未旋转坐标系或者动网格时,此选项勾选才有意义。表示计算阻力是否按照相对速度计算。若不是动网格,则勾选无意义。

     

    Viscous Resistance:

    多孔介质模型主要是对流动阻力产生影响,因此其会在动量方程中添加动量源项来充当阻力。

    x方向上的动量源项为下图,第一项为粘性损失项,第二项为惯性损失项。y方向和z方向同理,当为各向异性时,各方向上的粘性阻力系数可互不相同,惯性阻力系数也可不同。

    Viscous Resistance表示粘性阻力系数,方向1表示Du,方向2表示Dv(y方向粘性阻力系数),方向3表示Dw(z方向粘性阻力系数)

    Inertial Resistance:

    表示惯性阻力系数,方向1表示Cu,方向2表示Cv(y方向惯性阻力系数),方向3表示Cw(z方向惯性阻力系数)

    Alternative Formulation:

    当多孔介质为各向异性时,第一个方向的阻力系数与第二个方向阻力系数数量级上差距较大时,会引起收敛问题,勾选此选项可加快收敛。

    Fluid Porosity:

    孔隙率,表示流体所占据的体积分数,影响动量和能量方程

    9 设置边界条件

    inlet:采用velocity-inlet边界,速度为22.6m/s,在 Turbulence中的 Specification Method中选择 Intensity and Hydraulic Diameter,在 Turbulent Intensity中填入10,在 Hydraulic Diameter中填入42,单击OK按钮确认。

    outlet:设置为pressure-outlet,在 Turbulence中的 Specification Method中选择 Intensity and Hydraulic Diameter,在 Turbulent Intensity中填入5,在 Hydraulic Diameter中填入42,单击OK按钮确认。

    wall和substrate-wall:wall边界,保持默认设置

    Porous-in和porous-out:interior

    10 求解方法

    选择Coupled,Gradient栏选择Green-Gauss Cell Based,动量选择二阶迎风格式,勾选伪瞬态选项

    Solution Controls设置保持默认

    11 初始化

    选择标准初始化,Compute From选择inlet,单击Initialize,完成初始化

    12 计算设置

    设置计算步长100步

    13 后处理

    速度云图,选择y=0面显示速度云图

    通过速度云图,我们发现,当流体运动到多孔介质时,流体速度明显降低,这是因为多孔介质产生阻力的原因

    选择x=130、165、95等面显示速度云图(关于这些面的创建,我们以后会专门讲解)

    若不勾选多孔介质区域,则速度云图如下

     

    以上案例的cas和dat文件均可以免费获取,需要的朋友,只需要点赞关注收藏一键三连后私信我即可哦

    展开全文
  • 平面多孔介质的分数阶理论描述,张亚鹏,高峰,本文利用分数阶微积分的展开式建立多孔介质的本构模型。多孔介质由于在其形状上具有统计自相似性,因此其可以应用分数维进行描述
  • 多孔介质中的超绝热燃烧是一种先进的燃烧技术,具有高效低污染的特点。将这一技术应用于发动机领域,有可能引起内燃机技术和产业的一场重大革新。介绍超绝热燃烧的概念,讨论了多孔介质中往复流动下超绝热燃烧的特点...
  • 使用Fluent模拟多孔介质,进而模拟固体区域

    简介

    Fluent多孔介质仿真的具体步骤可以参考【简单学微课堂】Fluent多孔介质应用
    需要注意多孔介质区与流体介质区都属于“流体”材料。二者的交界面一定要设成interior

    测试对象

    在这里插入图片描述

    如图为以旋转对称的圆管,内部有一固定圆柱体状的“多孔介质”。空气从左侧流入,右侧流出。

    仿真结果

    1 不设置多孔介质区

    如果全设置为空气区,则流速场如下图
    在这里插入图片描述
    静压场如下
    在这里插入图片描述

    2 设置多孔介质区为0阻尼

    在这里插入图片描述
    计算得到的速度场
    在这里插入图片描述
    计算得到的压力场
    在这里插入图片描述
    可见计算结果与不设置多孔介质区完全一致

    3 设置多孔介质区为低阻尼

    取阻尼系数Viscous Resistance为 1 0 7 10^7 107
    速度场流线
    在这里插入图片描述
    压力场
    在这里插入图片描述

    4 设置多孔介质区为高阻尼

    取阻尼系数Viscous Resistance为 1 0 10 10^{10} 1010
    流速场如下
    在这里插入图片描述
    流线如下
    在这里插入图片描述
    压力场如下
    在这里插入图片描述

    5 设置多孔介质区为极高阻尼

    取阻尼系数Viscous Resistance为 1 0 20 10^{20} 1020(收敛性变差)
    流速场如下
    在这里插入图片描述
    流线如下
    在这里插入图片描述
    局部速度矢量
    在这里插入图片描述
    压力场如下
    在这里插入图片描述
    阻尼系数Viscous Resistance为 1 0 20 10^{20} 1020 1 0 10 10^{10} 1010结果相差不大。多孔介质区已经成为事实上纯固体区

    总结

    使用多孔介质是可以模拟固体的。多孔介质模型通过修改NS方程的动量方程(源项)将固体视为了一种特殊的流体,流体、多孔介质和固体之间的区别消失了可以在统一的网格上求解

    多孔介质模型的具体算法可以参考Fluent手册中help/flu_ug/flu_ug_sec_bc_porous_media.html这一章

    展开全文
  • 多孔介质应用在不同行业中有无数的应用:汽车,航空航天,民用(商业,住宅),环境噪声控制和生物医学。 过去,涉及多孔材料的设计问题是通过“坐垫决定”解决的,这导致了多个/迭代的原型和实验,这既昂贵又费时...
  • 采用局部非热平衡模型,在方腔左侧壁面温度正弦波变化的边界条件下,应用SIMPLER算法数值模拟研究了固体骨架发热多孔介质方腔内的稳态非达西自然对流,主要探讨了不同正弦波波动参数Ⅳ、振幅4及方腔的高宽比M/L对方...
  • 多孔介质的建模具有许多应用程序,其技术可用于汽车,航空航天,石油勘探和生物医学领域。 这项工作集中在汽车内饰的噪声和振动(NV)开发上,但是这些思想可以转化为上述领域。 NV开发需要将NV目标设置为不同的级别...
  • 应用Fourier展开和Hankel变换,求解了弹性饱和多孔介质非轴对称Biot波动方程,得到半空间问题的一般解.引入了位移分量和应力分量组合量,建立了这些组合量Hankel变换之间的关系.用传递矩阵方法处理了弹性饱和多孔介质半...
  • 多孔介质模型可以应用于很多问题,如通过充满介质的流动、通过过滤纸、穿孔圆盘、流量分配器以及管道堆的流动。当你使用这一模型时,你就定义了一个具有多孔介质的单元区域,而且流动的压力损失由多孔介质的动量方程...
  • 大数据-算法-预混气体在多孔介质内燃烧的工业应用和数值模拟研究.pdf
  • 自然储层孔隙结构复杂,孔隙和孔喉共存且往往会表现出分形特征。
  • OpenFOAM-多孔介质模型(porous media)

    千次阅读 2021-07-21 13:40:05
    文章目录一、多孔介质(porous media)二、使用步骤1.达西定律2.Darcy-Forchheimer3.OpenFOAM fvOptions4.paraview动画每帧显示时间(Annotate Time)5.参考内容总结 一、多孔介质(porous media) 在流体力学中,...

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 551
精华内容 220
关键字:

多孔介质的应用