二十六、多孔介质模型(二)-催化器
1.概念介绍
许多工业应用，如过滤器、催化剂床和填料，都涉及到模拟通过多孔介质的流动。本案例为气体通过催化转换器的流动。催化转换器通常用于净化汽油和柴油发动机的所产生的废气，将有害环境的废气转化为可接受的物质。
2.模型描述
这里的催化转化器模型如图所示，用于流动建模的催化转化器几何形状。氮气以125 m/s、800K的匀速通过进口，通过一对方形通道的陶瓷整体基板，从出口排出。
当入口和出口段中的流动为湍流时，通过基板的流动为层流，并且具有沿入口轴的惯性和粘性损失系数。基底在其他方向不可渗透。使用比主流方向高三个数量级的损失系数对该特性进行建模。
3.基本设置
1) 导入网格：
使用Fluent软件打开Chapter26.msh.gz网格文件，文件在本文末尾链接资源内。
2) 修改网格单位
Domain → Mesh → Units...
选中length，将长度单位修改为mm。
注意：此处只是将Fluent的默认长度单位由m改为了mm，并没有修改模型的尺寸。
3) 求解器设置
保持默认即可
4.设置计算模型
本例考虑多孔介质对温度的影响，因此打开能量方程
由于N2气体的流速125m/s非常大，流动达到了湍流，因此使用湍流模型。
SST k-ω模型湍流粘度考虑到了湍流剪应力的传播。相较于标准k-ω模型，SST k-ω模型在广泛的流动领域中有更高的精度和可信度。
5.材料设置
在Materials下，单击Create/Edit...，弹出如下面板。
Materials → Create/Edit...
单击Fluent Database，选择nitrogen（N2），单击copy，然后单击close，N2材料添加成功
6.设置计算域
单击信息树Cell Zone Conditions，
选择计算域fluid:0，双击fluid：0，弹出流体域设置对话框，将Material Name改为nitrogen，即将工作流体改为N2而不再是空气，其他设置保持默认。
对fluid:1和fluid:3进行同样的设置，将工作流体改为N2。
接下来就是多孔介质的设置了，多孔介质设置只需要在计算域进行设置即可，我们进行较为详细的讲解
双击fluid:substrate:1，弹出流体域设置对话框，勾选Porous Zone表明这部分计算域为多孔介质区域。
将Material Name改为nitrogen，表示多孔介质中的流体为N2，单击Porous Zone栏，对多孔介质进行设置。
1) Laminar Zone：
勾选Laminar Zone表明在多孔介质区域中抑制湍流的产生：如果不勾选，多孔介质区域的固体边界仿佛对湍流的产生和耗散没有影响，这个假定只有在渗透性很大并且多孔介质的几何尺度对湍流涡旋的几何尺度没有影响的情况下可行
2) Direction -1 Vector和Direction -1 Vector：
将多孔介质的主方向分别设置为z轴方向和y轴方向
多孔介质分为各向同性多孔介质，指的是多孔介质各个方向的阻力相同。各向异性指各个方向阻力不同，有的方向流体容易通过，有的方向流体很难通过。
若多孔介质为各向同性，此设置无意义。但若为各向异性，则方向1矢量表示多孔介质的第一个主方向为x方向，方向2矢量表示多孔介质的第二个主方向为y方向，第三个方向与这两个方向垂直，不必指定。
3) Viscous Resistance：
分别按下图设置三个方向的粘性阻力系数和惯性阻力系数
多孔介质模型主要是对流动阻力产生影响，因此其会在动量方程中添加动量源项来充当阻力。
x方向上的动量源项为下图，第一项为粘性损失项，第二项为惯性损失项。y方向和z方向同理，当为各向异性时，各方向上的粘性阻力系数可互不相同，惯性阻力系数也可不同。
Viscous Resistance表示粘性阻力系数，方向1表示Du，方向2表示Dv(y方向粘性阻力系数)，方向3表示Dw(z方向粘性阻力系数)
4) Inertial Resistance：
表示惯性阻力系数，方向1表示Cu，方向2表示Cv(y方向惯性阻力系数)，方向3表示Cw(z方向惯性阻力系数)
对计算域fluid:substrate:2进行同样的设置。
这里使用copy会更加方便
单击copy按钮，把刚才设置的计算域fluid:substrate:1复制给fluid:substrate:2。
From：被复制的计算域
To：要进行复制参数的计算域
单击copy后，弹出下框，点击ok即可。
Copy在进行边界条件设置时同样适用，且使用方法相同。
其余设置保持默认，可参照文章二十。
7.设置边界条件
inlet：采用velocity-inlet边界，速度为125m/s，在 Turbulence中的 Specification Method中选择 Intensity and Hydraulic Diameter，在 Turbulent Intensity中填入5，在 Hydraulic Diameter中填入500，单击OK按钮确认。
Thermal中流体进口温度设置为800K
outlet：设置为pressure-outlet，在 Turbulence中的 Specification Method中选择 Intensity and Hydraulic Diameter，在 Turbulent Intensity中填入5，在 Hydraulic Diameter中填入500，回流温度设置为800K，单击OK按钮确认。
其它边界，保持默认设置
8.求解方法
选择Coupled，动量选择二阶迎风格式，勾选伪瞬态选项
监视器设置
Solution → Reports → Definitions → New → Surface Report → Mass Flow Rate
创建一个监视器，用于监控出口流量
Name名称任意，勾选Report File、Report Plot和Print to Console。Report File表示生成一个用于存储出口流量数据的文件，Report Plot表示窗口显示流量曲线图，Print to Console表示流量数据会在控制面板打印出来。
Solution Controls设置保持默认
9.初始化
选择标准初始化，Compute From选择inlet，单击Initialize，完成初始化
对于多孔介质问题，推荐使用Standard标准化，若使用混合标准化，可能会产生不合实际的初始速度场，难以收敛。
10.计算设置
设置计算步长150步
11.后处理
1) 出口流量曲线图
2) 速度矢量图
创建面y=-425，在该面上显示速度矢量图
3) 各面的速度云图
在XY平面上创建各个面，显示各面的速度云图如下。
4) Wall面上的静压分布云图
文章创作不易，从建模到画网格，再到编写代码、画示意图都是我一点点摸索出来的，一篇原创文章可能要花费很久的时间。希望对大家有用的话，能够分享、在看、点赞支持一下。
为了方便大家学习交流，我创建了一个QQ群：975634031，大家可以在其中讨论相关的问题，同时我会将所有文章的源文件都放到QQ群中，还会放一些其它的学习文件。欢迎大家加入。
源文件链接：
链接：https://pan.baidu.com/s/1Z5dHlOAKkFvtTCNpTdDcnQ
提取码：siwg
