完整word版,fluent中多孔介质设置问题和算例
FLUENT多孔介质数值模拟设置
FLUENT多孔介质数值模拟设置FLUEN■多孔介质数值模拟设置多孔介质条件多孔介质模型可以应用于很多问题，如通过充满介质的流动、通过过滤纸、穿孔圆盘、流量分配器以及管道堆的流动。
当你使用这一模型时，你就定义了一个具有多孔介质的单元区域，而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定。
通过介质的热传导问题也可以得到描述，它服从介质和流体流动之间的热平衡假设，具体内容可以参考多孔介质中能量方程的处理一节。
多孔介质的一维化简模型，被称为多孔跳跃，可用于模拟具有已知速度/压降特征的薄膜。
多孔跳跃模型应用于表面区域而不是单元区域，并且在尽可能的情况下被使用(而不是完全的多孔介质模型)，这是因为它具有更好的鲁棒性，并具有更好的收敛性。
详细内容请参阅多孔跳跃边界条件。
多孔介质模型的限制如下面各节所述，多孔介质模型结合模型区域所具有的阻力的经验公式被定义为“多孔”。
事实上多孔介质不过是在动量方程中具有了附加的动量损失而已。
因此，下面模型的限制就可以很容易的理解了。
—I 二流体通过介质时不会加速，因为事实上出现的体积的阻塞并没有在模型中出现。
这对于过渡流是有很大的影响的，因为它意味着FLUENT S 会正确的描述通过介质的过渡时间。
多孔介质对于湍流的影响只是近似的。
详细内容可以参阅湍流多孔介质的处理一节。
多孔介质的动量方程多孔介质的动量方程具有附加的动量源项。
源项由两部分组成，一部分是粘性损失项(Darcy)，另一个是内部损失项：其中S_i是i向(x, y, or z) 动量源项，D和C是规定的矩阵。
在多孔介质单元中，动量损失对于压力梯度有贡献，压降和流体速度(或速度方阵)成比例。
对于简单的均匀多孔介质：其中a是渗透性，C_2时内部阻力因子，简单的指定D和C分别为对角阵1/a 和C_2其它项为零。
FLUENT?允许模拟的源项为速度的幕率：其中C_0和C_1为自定义经验系数。
注意：在幕律模型中，压降是各向同性的，C_0的单位为国际标准单位。
fluent多孔介质简单操作
[转]fluent中多孔介质porous media设置问题
经过痛苦的一段经历，终于将局部问题真相大白，为了使保位同仁不再经过我之痛苦，现在将本人多孔介质经验公布如下，希望各位能加精：
1。
划分网格之后，定义需要做为多孔介质的区域为fluid,与缺省的fluid分别开来，再定义其名称，我习惯将名称定义为porous;
2。
选中porous zone与laminar复选框,再点击porous zone标签即出现一个带有滚动条的界面；
3。
porous zone设置方法：
1）定义矢量：二维定义一个矢量，第二个矢量方向不用定义，是与第一个矢量方向正交的；
三维定义二个矢量，第三个矢量方向不用定义，是与第一、二个矢量方向正交的；
（如何知道矢量的方向：打开grid图，看看X,Y,Z的方向，如果是X向，矢量为1,0,0,同理Y向为0,1,0，Z向为0,0,1,如果所需要的方向与坐标轴正向相反，则定义矢量为负）
圆锥坐标与球坐标请参考fluent帮助。
2）定义粘性阻力1/a与内部阻力C2：请参看本人上一篇博文“终于搞清fluent中多孔粘性阻力与内部阻力的计算方法”，此处不赘述；
3）如果了定义粘性阻力1/a与内部阻力C2,就不用定义C1与C0，因为这是两种不同的定义方法，C1与C0只在幂率模型中出现，该处保持默认就行了；
4）定义孔隙率porousity，默认值1表示全开放，此值按实验测值填写即可。
完了，其他设置与普通k-e或RSM相同。
总结一下，与君共享!。
多孔介质系数计算方法(含模型文件和Excel表格下载)
多孔介质系数的计算方法1.STAR-CCM+中的多孔介质系数（文中所述的Excel文档以及模型源文件均可通过文中图片获取下载地址）在STAR-CCM+中主要需要设置多孔惯性阻力、多孔粘性阻力、孔隙率等参数。
孔隙率为流体在多孔介质中所占的分数，由标量分布法指定，取值介于0-1之间。
通过多孔介质的流量可根据达西定律描述，它基于渗透力度量以使流速与压力梯度相关联。
单相流体压降Δp(压力梯度)指定如下：Δp=−ρ(α|v n|+β)v n式中v n为表面法向速度，官方文档里写的α和β分被为多孔惯性阻力和多孔粘性阻力，但通过量纲分析，在软件参数输入区域α∙ρ为多孔惯性阻力[kg/m4]，β∙ρ为多孔粘性阻力[kg/m3·s]，ρ为交界面处流体密度。
在软件中，多孔介质需要设定单独的区域，并在类型属性中设定为多孔区域，多孔区域与流体区域通过内部交界面传递数据。
展开多孔区域的物理值菜单，需要输入多孔惯性阻力和多孔粘性阻力，分别表示单位长度的惯性量和粘性量值。
下图以正交各项异性参数为例展示需要输入参数的位置。
多孔惯性阻力和多孔粘性阻力的值通过试验数据的拟合求得，以散热器为例，对散热器芯体进行单体性能试验，得到所通过的气体(密度为1.18415kg/m3)流量和总压降的关系，转化为速度和总压降的关系如下所示。
在Excel中对数据进行拟合，将总压力转换为压力梯度，通过计算可得多孔惯性阻力和多孔粘性阻力。
得到的多孔惯性阻力和多孔粘性阻力是沿流动方向的阻力系数，另外两个方向通常取为比主流方向阻力系数大2-3个数量级的值，以阻止流体在这两个方向上的流动，但取值过大可能会影响收敛性。
下图为计算的Excel文件，填入试验数据即可得两个阻力系数。
下载本文档，将扩展名docx改为rar，得到一个压缩包文件，双击打开，再依次打开word、media文件夹。
然后将media文件夹下的image4.jpg拖拽至桌面或某个文件夹，将其扩展名jpg改为rar，解压缩可得一个包含模型和Excel文件下载地址的txt文本文档download.txt。
Fluent计算多孔介质模型资料
广东省深圳市宝安区沙井辛养社区西部工业园 TEL:+86-755-3366-8888 FAX:+86-755-3366-0612Fluent计算多孔介质模型资料这是一个多孔介质例子，进口速度为0.01m/s,组份为液态水和氧气，其中氧气从多孔介质porous jump 渗透过去，如何看氧气在tissue中扩散的。
porous jump的face permeability1 a=e-8 m_2thickness 设为0.0001pressure jump coefficient为默认porous zone设置如下：direction vector 1, 1,viscous resistance 100 eachinertial resistance 100 eachporosity 0.1边界条件设置如下：Ab – wall - defaultBc – wall – defaultBe – porous jump – face permeability 1e-8, porous medium thickness0.0001Cd – outflow rating – 0.5De – wall – defaultDefault interior – interiorDefault interior001 – interiorDefault interior019 – interiorEf – wall - defaultFg – outflow rating – 1Fluid - porous zone - direction vector 1, 1, viscous resistance 100 each,inertial resistance 100 each, porosity 0.1Gh- wall - defaultHi – wall - defaultHk - porous jump same conditions as otherIj – outflow – 0.5Jk – wall – defaultKl – wall – defaultLa – velocity inlet – 0.01 m/s, temperature 300K, 0.5 mass fraction O2 Lfluid – porous zone - direction vector 1, 1, viscous resistance 100 each,inertial resistance 100 each, porosity 0.1Pipefluid – fluid – default (no porous zone)Models – species transport – water and oxygen mixtureVariations – different boundary conditions at top and bottom (outflow, wall ect)注意,其中porous zone在gambit中设置为fluid，在fluent中设置为porous zone边界条件设置如下：Ab – wall - defaultBc – wall – defaultBe – porous jump – face permeability 1e-8, porous medium thickness0.0001Cd – outflow rating – 0.5De – wall – defaultDefault interior – interiorDefault interior001 – interiorDefault interior019 – interiorEf – wall - defaultFg – outflow rating – 1Fluid - porous zone - direction vector 1, 1, viscous resistance 100 each,inertial resistance 100 each, porosity 0.1Gh- wall - defaultHi – wall - defaultHk - porous jump same conditions as otherIj – outflow – 0.5Jk – wall – defaultKl – wall – defaultLa – velocity inlet – 0.01 m/s, temperature 300K, 0.5 mass fraction O2 Lfluid – porous zone - direction vector 1, 1, viscous resistance 100 each,inertial resistance 100 each, porosity 0.1Pipefluid – fluid – default (no porous zone)Models – species transport – water and oxygen mixtureVariations – different boundary conditions at top and bottom (outflow, wall ect) 注意,其中porous zone在gambit中设置为fluid，在fluent中设置为porous zone。
FLUENT多孔介质数值模拟设置
FLUENT多孔介质数值模拟设置多孔介质条件多孔介质模型可以应用于很多问题，如通过充满介质的流动、通过过滤纸、穿孔圆盘、流量分配器以及管道堆的流动。
当你使用这一模型时，你就定义了一个具有多孔介质的单元区域，而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定。
通过介质的热传导问题也可以得到描述，它服从介质和流体流动之间的热平衡假设，具体内容可以参考多孔介质中能量方程的处理一节。
多孔介质的一维化简模型，被称为多孔跳跃，可用于模拟具有已知速度/压降特征的薄膜。
多孔跳跃模型应用于表面区域而不是单元区域，并且在尽可能的情况下被使用（而不是完全的多孔介质模型），这是因为它具有更好的鲁棒性，并具有更好的收敛性。
详细内容请参阅多孔跳跃边界条件。
多孔介质模型的限制如下面各节所述，多孔介质模型结合模型区域所具有的阻力的经验公式被定义为“多孔”。
事实上多孔介质不过是在动量方程中具有了附加的动量损失而已。
因此，下面模型的限制就可以很容易的理解了。
流体通过介质时不会加速，因为事实上出现的体积的阻塞并没有在模型中出现。
这对于过渡流是有很大的影响的，因为它意味着FLUENT不会正确的描述通过介质的过渡时间。
多孔介质对于湍流的影响只是近似的。
详细内容可以参阅湍流多孔介质的处理一节。
多孔介质的动量方程多孔介质的动量方程具有附加的动量源项。
源项由两部分组成，一部分是粘性损失项 (Darcy)，另一个是内部损失项：其中S_i是i向(x, y, or z)动量源项，D和C是规定的矩阵。
在多孔介质单元中，动量损失对于压力梯度有贡献，压降和流体速度（或速度方阵）成比例。
对于简单的均匀多孔介质：其中a是渗透性，C_2时内部阻力因子，简单的指定D和C分别为对角阵1/a 和C_2其它项为零。
FLUENT还允许模拟的源项为速度的幂率：其中C_0和C_1为自定义经验系数。
注意：在幂律模型中，压降是各向同性的，C_0的单位为国际标准单位。
多孔介质的Darcy定律通过多孔介质的层流流动中，压降和速度成比例，常数C_2可以考虑为零。
FLUENT多孔介质数值模拟设置
FLUENT多孔介质数值模拟设置C=对于不同D/t的不同雷诺数范围被列成不同的表的系数A_p=圆盘的面积（固体和洞)如果你选择在多孔介质中模拟热传导，你必须指定多孔介质中的材料以及多孔性。
要定义多孔介质的材料，向下拉流体面板中阻力输入底下的滚动条，然后在多孔热传导的固体材料下拉列表中选中适当的固体。
另一个处理收敛性差的要领是临时取消多孔介质模型（在流体面板中关闭多孔区域）然后获取一个不受多孔区域影响的初始流场。
取消多孔区域后，FLUEＮT会将多孔区域处理为流体区域并按响应的流体区域来计算。
一旦获取了初始解，或者计算很容易收敛，你就可以激活多孔模型继续计算包罗多孔区域的流场（对于大阻力多孔介质不保举使用该要领）。
这些变量会在后处理面板的变量选择下拉菜谱制定类别中出现。
然后在多孔热传导下设定多孔性。
多孔性f是多孔介质中流体的体积分数（即介质的开放体积分数）。
多孔性用于介质中的热传导预测，处理要领请参阅多孔介质能量方程的处理一节。
它还对介质中的反应源项和体力的计算有影响。
这个源项和介质中流体的体积成比例。
如果你想要模拟完全开放的介质（固体介质没有影响），你应该设定多孔性为1.0。
当多孔性为1.0时，介质的固体部门对于热传导和（或）热源项/反应源项没有影响。
注意：多孔性永远不会影响介质中的流体速率，这已经在多孔介质的动量方程一节中介绍了。
不管你将多孔性设定为何值，，FLUEＮT所预测的速率都是介质中的外貌速率。
对于多孔介质动量源项（多孔介质动量方程中的方程5），如果你使用幂律模型近似，你只要在流体面板的幂律模型中输入系数C_0和C_1就可以了。
如果C_0或C_1为非零值，解算器会忽略面板中除了多孔介质幂律模型之外的所有输入。
定义源项一般说来，在模拟多孔介质时，你可以使用标准的解算步骤以及解参数的设置。
然而你会发现如果多孔区域在流动方向上压降至关大（比如：渗透性a很低或者内部因数C_2很大）的话，解的收敛速率就会变慢。
FLUENT中定义多孔介质的方法
FLUENT中定义多孔介质的方法多孔介质计算中需要输入的项目如下：（1）定义多孔介质区域。
（2）定义多孔介质速度函数形式。
（3）定义流过多孔介质区的流体属性。
（4）设定多孔区的化学反应。
（5）设定粘性阻力系数。
（6）设定多孔介质的多孔率。
（7）在计算热交换的过程中选择多孔介质的材料。
（8）设定多孔介质固体部分的体热生成率。
（9）设定流动区域上的任意固定值的流动参数。
（10）需要的话，将多孔区流动设为层流，或取消湍流计算。
（11）定义旋转轴或区域的运动。
经过痛苦的一段经历，终于将局部问题真相大白，为了使保位同仁不再经过我之痛苦，现在将本人多孔介质经验公布如下，希望各位能加精：1。
Gambit中划分网格之后，定义需要做为多孔介质的区域为fluid,与缺省的fluid 分别开来，再定义其名称，我习惯将名称定义为porous;2。
在fluent中定义边界条件define-boundary condition-porous(刚定义的名称)，将其设置边界条件为fluid,点击set按钮即弹出与fluid边界条件一样的对话框，选中porous zone与laminar复选框,再点击porous zone标签即出现一个带有滚动条的界面；3。
porous zone设置方法：1）定义矢量：二维定义一个矢量，第二个矢量方向不用定义，是与第一个矢量方向正交的；三维定义二个矢量，第三个矢量方向不用定义，是与第一、二个矢量方向正交的；（如何知道矢量的方向：打开grid图，看看X,Y,Z的方向，如果是X向，矢量为1,0,0,同理Y向为0,1,0，Z向为0,0,1,如果所需要的方向与坐标轴正向相反，则定义矢量为负）圆锥坐标与球坐标请参考fluent帮助。
2）定义粘性阻力1/a与内部阻力C2：见附fluent中多孔粘性阻力系数与惯性阻力系数的方法：对比Darcy 定律v K J =⋅J 为水力坡度，即流经路径长度为L 的水头损失，m/m ；K 为比例系数，也称渗流系数，m/sg p v Kρ∆=- p μνα∆=- 这样可以得到，g K ρμα=，所以1g K ραμ= 其中： 1α为粘性阻力系数，1/m 2； K 为比例系数，也称渗流系数，m/s ；μ为动力粘度，Pa.s ；ρ为流体密度，kg/m 3；g 为重力加速度，m/s 2；采用20℃饱和水的物性参数：密度为998.2kg/m 3 ，运动粘度1.004×10-3Pa.s ；渗流速度取1.0×10-6m/s ，K 取1.0×10-4m/s ，水力坡度约为0.01m/m ； 10461998.29.89.74101.010100410g K ραμ--⨯===⨯⨯⨯⨯ 3）如果了定义粘性阻力1/a 与内部阻力C2,就不用定义C1与C0，因为这是两种不同的定义方法，C1与C0只在幂率模型中出现，该处保持默认就行了；4）定义孔隙率porousity ，默认值1表示全开放，此值按实验测值填写即可。
(完整word版)fluent中多孔介质设置问题和算例
经过痛苦的一段经历，终于将局部问题真相大白，为了使保位同仁不再经过我之痛苦，现在将本人多孔介质经验公布如下，希望各位能加精：1。
Gambit中划分网格之后，定义需要做为多孔介质的区域为fluid,与缺省的fluid分别开来，再定义其名称，我习惯将名称定义为porous;2。
在fluent中定义边界条件define-boundary condition-porous(刚定义的名称)，将其设置边界条件为fluid,点击set按钮即弹出与fluid边界条件一样的对话框，选中porous zone与laminar复选框,再点击porous zone标签即出现一个带有滚动条的界面；3。
porous zone设置方法：1）定义矢量：二维定义一个矢量，第二个矢量方向不用定义，是与第一个矢量方向正交的；三维定义二个矢量，第三个矢量方向不用定义，是与第一、二个矢量方向正交的；（如何知道矢量的方向：打开grid图，看看X,Y,Z的方向，如果是X向，矢量为1,0,0,同理Y向为0,1,0，Z向为0,0,1,如果所需要的方向与坐标轴正向相反，则定义矢量为负）圆锥坐标与球坐标请参考fluent帮助。
2）定义粘性阻力1/a与内部阻力C2：请参看本人上一篇博文“终于搞清fluent中多孔粘性阻力与内部阻力的计算方法”，此处不赘述；3）如果了定义粘性阻力1/a与内部阻力C2,就不用定义C1与C0，因为这是两种不同的定义方法，C1与C0只在幂率模型中出现，该处保持默认就行了；4）定义孔隙率porousity，默认值1表示全开放，此值按实验测值填写即可。
完了，其他设置与普通k-e或RSM相同。
总结一下，与君共享!Tutorial 7. Modeling Flow Through Porous MediaIntroductionMany industrial applications involve the modeling of flow through porous media, such as filters, catalyst beds, and packing. This tutorial illustrates how to set up and solve a problem involving gas flow through porous media.The industrial problem solved here involves gas flow through a catalytic converter. Catalytic converters are commonly used to purify emissions from gasoline and diesel engines by converting environmentally hazardous exhaust emissions to acceptable substances.Examples of such emissions include carbon monoxide (CO), nitrogen oxides (NOx), and unburned hydrocarbon fuels. These exhaust gas emissions are forced through a substrate, which is a ceramic structure coated with a metal catalyst such as platinum or palladium.The nature of the exhaust gas flow is a very important factor in determining the performance of the catalytic converter. Of particular importance is the pressure gradient and velocity distribution through the substrate. Hence CFD analysis is used to design efficient catalytic converters: by modeling the exhaust gas flow, the pressure drop and the uniformity of flow through the substrate can be determined. In this tutorial, FLUENT is used to model the flow of nitrogen gas through a catalytic converter geometry, so that the flow field structure may be analyzed.This tutorial demonstrates how to do the following:_ Set up a porous zone for the substrate with appropriate resistances._ Calculate a solution for gas flow through the catalytic converter using the pressure based solver. _ Plot pressure and velocity distribution on specified planes of the geometry._ Determine the pressure drop through the substrate and the degree of non-uniformity of flow through cross sections of the geometry using X-Y plots and numerical reports.Problem DescriptionThe catalytic converter modeled here is shown in Figure 7.1. The nitrogen flows in through the inlet with a uniform velocity of 22.6 m/s, passes through a ceramic monolith substrate with square shaped channels, and then exits through the outlet.While the flow in the inlet and outlet sections is turbulent, the flow through the substrate is laminar and is characterized by inertial and viscous loss coefficients in the flow (X) direction. The substrate is impermeable in other directions, which is modeled using loss coefficients whose values are three orders of magnitude higher than in the X direction.Setup and SolutionStep 1: Grid1. Read the mesh file (catalytic converter.msh).File /Read /Case...2. Check the grid. Grid /CheckFLUENT will perform various checks on the mesh and report the progress in the console. Make sure that the minimum volume reported is a positive number.3. Scale the grid.Grid! Scale...(a) Select mm from the Grid Was Created In drop-down list.(b) Click the Change Length Units button. All dimensions will now be shown in millimeters.(c) Click Scale and close the Scale Grid panel.4. Display the mesh. Display /Grid...(a) Make sure that inlet, outlet, substrate-wall, and wall are selected in the Surfaces selection list.(b) Click Display.(c) Rotate the view and zoom in to get the display shown in Figure 7.2.(d) Close the Grid Display panel.The hex mesh on the geometry contains a total of 34,580 cells.Step 2: Models1. Retain the default solver settings. Define /Models /Solver...2. Select the standard k-ε turbulence model. Define/ Models /Viscous...Step 3: Materials1. Add nitrogen to the list of fluid materials by copying it from the Fluent Database for materials.Define /Materials...(a) Click the Fluent Database... button to open the Fluent Database Materials panel.i. Select nitrogen (n2) from the list of Fluent Fluid Materials.ii. Click Copy to copy the information for nitrogen to your list of fluid materials. iii. Close the Fluent Database Materials panel.(b) Close the Materials panel.Step 4: Boundary Conditions. Define /Boundary Conditions...1. Set the boundary conditions for the fluid (fluid).(a) Select nitrogen from the Material Name drop-down list.(b) Click OK to close the Fluid panel.2. Set the boundary conditions for the substrate (substrate).(a) Select nitrogen from the Material Name drop-down list.(b) Enable the Porous Zone option to activate the porous zone model.(c) Enable the Laminar Zone option to solve the flow in the porous zone without turbulence.(d) Click the Porous Zone tab.i. Make sure that the principal direction vectors are set as shown in Table7.1. Use the scroll bar to access the fields that are not initially visible in the panel.ii. Enter the values in Table 7.2 for the Viscous Resistance and Inertial Resistance. Scroll down to access the fields that are not initially visible in the panel.(e) Click OK to close the Fluid panel.3. Set the velocity and turbulence boundary conditions at the inlet (inlet).(a) Enter 22.6 m/s for the Velocity Magnitude.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdown list in the Turbulence group box.(c) Retain the default value of 10% for the Turbulent Intensity.(d) Enter 42 mm for the Hydraulic Diameter.(e) Click OK to close the Velocity Inlet panel.4. Set the boundary conditions at the outlet (outlet).(a) Retain the default setting of 0 for Gauge Pressure.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdown list in the Turbulence group box.(c) Enter 5% for the Backflow Turbulent Intensity.(d) Enter 42 mm for the Backflow Hydraulic Diameter.(e) Click OK to close the Pressure Outlet panel.5. Retain the default boundary conditions for the walls (substrate-wall and wall) and close the Boundary Conditions panel.Step 5: Solution1. Set the solution parameters. Solve /Controls /Solution...(a) Retain the default settings for Under-Relaxation Factors.(b) Select Second Order Upwind from the Momentum drop-down list in the Discretization group box.(c) Click OK to close the Solution Controls panel.2. Enable the plotting of residuals during the calculation. Solve/Monitors /Residual...(a) Enable Plot in the Options group box.(b) Click OK to close the Residual Monitors panel.3. Enable the plotting of the mass flow rate at the outlet.Solve / Monitors /Surface...(a) Set the Surface Monitors to 1.(b) Enable the Plot and Write options for monitor-1, and click the Define... button to open the Define Surface Monitor panel.i. Select Mass Flow Rate from the Report Type drop-down list.ii. Select outlet from the Surfaces selection list.iii. Click OK to close the Define Surface Monitors panel.(c) Click OK to close the Surface Monitors panel.4. Initialize the solution from the inlet. Solve /Initialize /Initialize...(a) Select inlet from the Compute From drop-down list.(b) Click Init and close the Solution Initialization panel.5. Save the case file (catalytic converter.cas). File /Write /Case...6. Run the calculation by requesting 100 iterations. Solve /Iterate...(a) Enter 100 for the Number of Iterations.(b) Click Iterate.The FLUENT calculation will converge in approximately 70 iterations. By this point the mass flow rate monitor has attended out, as seen in Figure 7.3.(c) Close the Iterate panel.7. Save the case and data files (catalytic converter.cas and catalytic converter.dat).File /Write /Case & Data...Note: If you choose a file name that already exists in the current folder, FLUENTwill prompt you for confirmation to overwrite the file.Step 6: Post-processing1. Create a surface passing through the centerline for post-processing purposes.Surface/Iso-Surface...(a) Select Grid... and Y-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Retain the default value of 0 for the Iso-Values.(d) Enter y=0 for the New Surface Name.(e) Click Create.2. Create cross-sectional surfaces at locations on either side of the substrate, as well as at its center.Surface /Iso-Surface...(a) Select Grid... and X-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Enter 95 for Iso-Values.(d) Enter x=95 for the New Surface Name.(e) Click Create.(f) In a similar manner, create surfaces named x=130 and x=165 with Iso-Values of 130 and 165, respectively. Close the Iso-Surface panel after all the surfaces have been created.3. Create a line surface for the centerline of the porous media.Surface /Line/Rake...(a) Enter the coordinates of the line under End Points, using the starting coordinate of (95, 0, 0) and an ending coordinate of (165, 0, 0), as shown.(b) Enter porous-cl for the New Surface Name.(c) Click Create to create the surface.(d) Close the Line/Rake Surface panel.4. Display the two wall zones (substrate-wall and wall). Display /Grid...(a) Disable the Edges option.(b) Enable the Faces option.(c) Deselect inlet and outlet in the list under Surfaces, and make sure that only substrate-wall and wall are selected.(d) Click Display and close the Grid Display panel.(e) Rotate the view and zoom so that the display is similar to Figure 7.2.5. Set the lighting for the display. Display /Options...(a) Enable the Lights On option in the Lighting Attributes group box.(b) Retain the default selection of Gourand in the Lighting drop-down list.(c) Click Apply and close the Display Options panel.6. Set the transparency parameter for the wall zones (substrate-wall and wall).Display/Scene...(a) Select substrate-wall and wall in the Names selection list.(b) Click the Display... button under Geometry Attributes to open the Display Properties panel.i. Set the Transparency slider to 70.ii. Click Apply and close the Display Properties panel.(c) Click Apply and then close the Scene Description panel.7. Display velocity vectors on the y=0 surface.Display /Vectors...(a) Enable the Draw Grid option. The Grid Display panel will open.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(b) Enter 5 for the Scale.(c) Set Skip to 1.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Vectors panel.The flow pattern shows that the flow enters the catalytic converter as a jet, with recirculation on either side of the jet. As it passes through the porous substrate, it decelerates and straightens out, and exhibits a more uniform velocity distribution.This allows the metal catalyst present in the substrate to be more effective.Figure 7.4: Velocity Vectors on the y=0 Plane8. Display filled contours of static pressure on the y=0 plane.Display /Contours...(a) Enable the Filled option.(b) Enable the Draw Grid option to open the Display Grid panel.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(c) Make sure that Pressure... and Static Pressure are selected from the Contours of drop-down lists.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Contours panel.Figure 7.5: Contours of the Static Pressure on the y=0 planeThe pressure changes rapidly in the middle section, where the fluid velocity changes as it passes through the porous substrate. The pressure drop can be high, due to the inertial and viscous resistance of the porous media. Determining this pressure drop is a goal of CFD analysis. In the next step, you will learn how to plot the pressure drop along the centerline of the substrate.9. Plot the static pressure across the line surface porous-cl.Plot /XY Plot...(a) Make sure that the Pressure... and Static Pressure are selected from the Y Axis Function drop-down lists.(b) Select porous-cl from the Surfaces selection list.(c) Click Plot and close the Solution XY Plot panel.Figure 7.6: Plot of the Static Pressure on the porous-cl Line SurfaceIn Figure 7.6, the pressure drop across the porous substrate can be seen to be roughly 300 Pa.10. Display filled contours of the velocity in the X direction on the x=95, x=130 and x=165 surfaces.Display /Contours...(a) Disable the Global Range option.(b) Select Velocity... and X Velocity from the Contours of drop-down lists.(c) Select x=130, x=165, and x=95 from the Surfaces selection list, and deselect y=0.(d) Click Display and close the Contours panel.The velocity profile becomes more uniform as the fluid passes through the porous media. The velocity is very high at the center (the area in red) just before the nitrogen enters the substrate and then decreases as it passes through and exits the substrate. The area in green, which corresponds to a moderate velocity, increases in extent.Figure 7.7: Contours of the X Velocity on the x=95, x=130, and x=165 Surfaces11. Use numerical reports to determine the average, minimum, and maximum of the velocity distribution before and after the porous substrate.Report /Surface Integrals...(a) Select Mass-Weighted Average from the Report Type drop-down list.(b) Select Velocity and X Velocity from the Field Variable drop-down lists.(c) Select x=165 and x=95 from the Surfaces selection list.(d) Click Compute.(e) Select Facet Minimum from the Report Type drop-down list and click Compute again.(f) Select Facet Maximum from the Report Type drop-down list and click Compute again.(g) Close the Surface Integrals panel.The numerical report of average, maximum and minimum velocity can be seen in the main FLUENT console, as shown in the following example:The spread between the average, maximum, and minimum values for X velocity gives the degree to which the velocity distribution is non-uniform. You can also use these numbers to calculate the velocity ratio (i.e., the maximum velocity divided by the mean velocity) and the space velocity (i.e., the product of the mean velocity and the substrate length).Custom field functions and UDFs can be also used to calculate more complex measures ofnon-uniformity, such as the standard deviation and the gamma uniformity index.SummaryIn this tutorial, you learned how to set up and solve a problem involving gas flow through porous media in FLUENT. You also learned how to perform appropriate post-processing to investigate the flow field, determine the pressure drop across the porous media and non-uniformity of the velocity distribution as the fluid goes through the porous media.Further ImprovementsThis tutorial guides you through the steps to reach an initial solution. You may be able to obtain a more accurate solution by using an appropriate higher-order discretization scheme and by adapting the grid. Grid adaption can also ensure that the solution is independent of the grid. These steps aredemonstrated in Tutorial 1.。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
经过痛苦的一段经历，终于将局部问题真相大白，为了使保位同仁不再经过我之痛苦，现在将本人多孔介质经验公布如下，希望各位能加精：1。
Gambit中划分网格之后，定义需要做为多孔介质的区域为fluid,与缺省的fluid分别开来，再定义其名称，我习惯将名称定义为porous;2。
在fluent中定义边界条件define-boundary condition-porous(刚定义的名称)，将其设置边界条件为fluid,点击set按钮即弹出与fluid边界条件一样的对话框，选中porous zone与laminar复选框,再点击porous zone标签即出现一个带有滚动条的界面；3。
porous zone设置方法：1）定义矢量：二维定义一个矢量，第二个矢量方向不用定义，是与第一个矢量方向正交的；三维定义二个矢量，第三个矢量方向不用定义，是与第一、二个矢量方向正交的；（如何知道矢量的方向：打开grid图，看看X,Y,Z的方向，如果是X向，矢量为1,0,0,同理Y向为0,1,0，Z向为0,0,1,如果所需要的方向与坐标轴正向相反，则定义矢量为负）圆锥坐标与球坐标请参考fluent帮助。
2）定义粘性阻力1/a与内部阻力C2：请参看本人上一篇博文“终于搞清fluent中多孔粘性阻力与内部阻力的计算方法”，此处不赘述；3）如果了定义粘性阻力1/a与内部阻力C2,就不用定义C1与C0，因为这是两种不同的定义方法，C1与C0只在幂率模型中出现，该处保持默认就行了；4）定义孔隙率porousity，默认值1表示全开放，此值按实验测值填写即可。
完了，其他设置与普通k-e或RSM相同。
总结一下，与君共享!Tutorial 7. Modeling Flow Through Porous MediaIntroductionMany industrial applications involve the modeling of flow through porous media, such as filters, catalyst beds, and packing. This tutorial illustrates how to set up and solve a problem involving gas flow through porous media.The industrial problem solved here involves gas flow through a catalytic converter. Catalytic converters are commonly used to purify emissions from gasoline and diesel engines by converting environmentally hazardous exhaust emissions to acceptable substances.Examples of such emissions include carbon monoxide (CO), nitrogen oxides (NOx), and unburned hydrocarbon fuels. These exhaust gas emissions are forced through a substrate, which is a ceramic structure coated with a metal catalyst such as platinum or palladium.The nature of the exhaust gas flow is a very important factor in determining the performance of the catalytic converter. Of particular importance is the pressure gradient and velocity distribution through the substrate. Hence CFD analysis is used to design efficient catalytic converters: by modeling the exhaust gas flow, the pressure drop and the uniformity of flow through the substrate can be determined. In this tutorial, FLUENT is used to model the flow of nitrogen gas through a catalytic converter geometry, so that the flow field structure may be analyzed.This tutorial demonstrates how to do the following:_ Set up a porous zone for the substrate with appropriate resistances._ Calculate a solution for gas flow through the catalytic converter using the pressure based solver. _ Plot pressure and velocity distribution on specified planes of the geometry._ Determine the pressure drop through the substrate and the degree of non-uniformity of flow through cross sections of the geometry using X-Y plots and numerical reports.Problem DescriptionThe catalytic converter modeled here is shown in Figure 7.1. The nitrogen flows in through the inlet with a uniform velocity of 22.6 m/s, passes through a ceramic monolith substrate with square shaped channels, and then exits through the outlet.While the flow in the inlet and outlet sections is turbulent, the flow through the substrate is laminar and is characterized by inertial and viscous loss coefficients in the flow (X) direction. The substrate is impermeable in other directions, which is modeled using loss coefficients whose values are three orders of magnitude higher than in the X direction.Setup and SolutionStep 1: Grid1. Read the mesh file (catalytic converter.msh).File /Read /Case...2. Check the grid. Grid /CheckFLUENT will perform various checks on the mesh and report the progress in the console. Make sure that the minimum volume reported is a positive number.3. Scale the grid.Grid! Scale...(a) Select mm from the Grid Was Created In drop-down list.(b) Click the Change Length Units button. All dimensions will now be shown in millimeters.(c) Click Scale and close the Scale Grid panel.4. Display the mesh. Display /Grid...(a) Make sure that inlet, outlet, substrate-wall, and wall are selected in the Surfaces selection list.(b) Click Display.(c) Rotate the view and zoom in to get the display shown in Figure 7.2.(d) Close the Grid Display panel.The hex mesh on the geometry contains a total of 34,580 cells.Step 2: Models1. Retain the default solver settings. Define /Models /Solver...2. Select the standard k-ε turbulence model. Define/ Models /Viscous...Step 3: Materials1. Add nitrogen to the list of fluid materials by copying it from the Fluent Database for materials.Define /Materials...(a) Click the Fluent Database... button to open the Fluent Database Materials panel.i. Select nitrogen (n2) from the list of Fluent Fluid Materials.ii. Click Copy to copy the information for nitrogen to your list of fluid materials. iii. Close the Fluent Database Materials panel.(b) Close the Materials panel.Step 4: Boundary Conditions. Define /Boundary Conditions...1. Set the boundary conditions for the fluid (fluid).(a) Select nitrogen from the Material Name drop-down list.(b) Click OK to close the Fluid panel.2. Set the boundary conditions for the substrate (substrate).(a) Select nitrogen from the Material Name drop-down list.(b) Enable the Porous Zone option to activate the porous zone model.(c) Enable the Laminar Zone option to solve the flow in the porous zone without turbulence.(d) Click the Porous Zone tab.i. Make sure that the principal direction vectors are set as shown in Table7.1. Use the scroll bar to access the fields that are not initially visible in the panel.ii. Enter the values in Table 7.2 for the Viscous Resistance and Inertial Resistance. Scroll down to access the fields that are not initially visible in the panel.(e) Click OK to close the Fluid panel.3. Set the velocity and turbulence boundary conditions at the inlet (inlet).(a) Enter 22.6 m/s for the Velocity Magnitude.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdown list in the Turbulence group box.(c) Retain the default value of 10% for the Turbulent Intensity.(d) Enter 42 mm for the Hydraulic Diameter.(e) Click OK to close the Velocity Inlet panel.4. Set the boundary conditions at the outlet (outlet).(a) Retain the default setting of 0 for Gauge Pressure.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdown list in the Turbulence group box.(c) Enter 5% for the Backflow Turbulent Intensity.(d) Enter 42 mm for the Backflow Hydraulic Diameter.(e) Click OK to close the Pressure Outlet panel.5. Retain the default boundary conditions for the walls (substrate-wall and wall) and close the Boundary Conditions panel.Step 5: Solution1. Set the solution parameters. Solve /Controls /Solution...(a) Retain the default settings for Under-Relaxation Factors.(b) Select Second Order Upwind from the Momentum drop-down list in the Discretization group box.(c) Click OK to close the Solution Controls panel.2. Enable the plotting of residuals during the calculation. Solve/Monitors /Residual...(a) Enable Plot in the Options group box.(b) Click OK to close the Residual Monitors panel.3. Enable the plotting of the mass flow rate at the outlet.Solve / Monitors /Surface...(a) Set the Surface Monitors to 1.(b) Enable the Plot and Write options for monitor-1, and click the Define... button to open the Define Surface Monitor panel.i. Select Mass Flow Rate from the Report Type drop-down list.ii. Select outlet from the Surfaces selection list.iii. Click OK to close the Define Surface Monitors panel.(c) Click OK to close the Surface Monitors panel.4. Initialize the solution from the inlet. Solve /Initialize /Initialize...(a) Select inlet from the Compute From drop-down list.(b) Click Init and close the Solution Initialization panel.5. Save the case file (catalytic converter.cas). File /Write /Case...6. Run the calculation by requesting 100 iterations. Solve /Iterate...(a) Enter 100 for the Number of Iterations.(b) Click Iterate.The FLUENT calculation will converge in approximately 70 iterations. By this point the mass flow rate monitor has attended out, as seen in Figure 7.3.(c) Close the Iterate panel.7. Save the case and data files (catalytic converter.cas and catalytic converter.dat).File /Write /Case & Data...Note: If you choose a file name that already exists in the current folder, FLUENTwill prompt you for confirmation to overwrite the file.Step 6: Post-processing1. Create a surface passing through the centerline for post-processing purposes.Surface/Iso-Surface...(a) Select Grid... and Y-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Retain the default value of 0 for the Iso-Values.(d) Enter y=0 for the New Surface Name.(e) Click Create.2. Create cross-sectional surfaces at locations on either side of the substrate, as well as at its center.Surface /Iso-Surface...(a) Select Grid... and X-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Enter 95 for Iso-Values.(d) Enter x=95 for the New Surface Name.(e) Click Create.(f) In a similar manner, create surfaces named x=130 and x=165 with Iso-Values of 130 and 165, respectively. Close the Iso-Surface panel after all the surfaces have been created.3. Create a line surface for the centerline of the porous media.Surface /Line/Rake...(a) Enter the coordinates of the line under End Points, using the starting coordinate of (95, 0, 0) and an ending coordinate of (165, 0, 0), as shown.(b) Enter porous-cl for the New Surface Name.(c) Click Create to create the surface.(d) Close the Line/Rake Surface panel.4. Display the two wall zones (substrate-wall and wall). Display /Grid...(a) Disable the Edges option.(b) Enable the Faces option.(c) Deselect inlet and outlet in the list under Surfaces, and make sure that only substrate-wall and wall are selected.(d) Click Display and close the Grid Display panel.(e) Rotate the view and zoom so that the display is similar to Figure 7.2.5. Set the lighting for the display. Display /Options...(a) Enable the Lights On option in the Lighting Attributes group box.(b) Retain the default selection of Gourand in the Lighting drop-down list.(c) Click Apply and close the Display Options panel.6. Set the transparency parameter for the wall zones (substrate-wall and wall).Display/Scene...(a) Select substrate-wall and wall in the Names selection list.(b) Click the Display... button under Geometry Attributes to open the Display Properties panel.i. Set the Transparency slider to 70.ii. Click Apply and close the Display Properties panel.(c) Click Apply and then close the Scene Description panel.7. Display velocity vectors on the y=0 surface.Display /Vectors...(a) Enable the Draw Grid option. The Grid Display panel will open.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(b) Enter 5 for the Scale.(c) Set Skip to 1.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Vectors panel.The flow pattern shows that the flow enters the catalytic converter as a jet, with recirculation on either side of the jet. As it passes through the porous substrate, it decelerates and straightens out, and exhibits a more uniform velocity distribution.This allows the metal catalyst present in the substrate to be more effective.Figure 7.4: Velocity Vectors on the y=0 Plane8. Display filled contours of static pressure on the y=0 plane.Display /Contours...(a) Enable the Filled option.(b) Enable the Draw Grid option to open the Display Grid panel.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(c) Make sure that Pressure... and Static Pressure are selected from the Contours of drop-down lists.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Contours panel.Figure 7.5: Contours of the Static Pressure on the y=0 planeThe pressure changes rapidly in the middle section, where the fluid velocity changes as it passes through the porous substrate. The pressure drop can be high, due to the inertial and viscous resistance of the porous media. Determining this pressure drop is a goal of CFD analysis. In the next step, you will learn how to plot the pressure drop along the centerline of the substrate.9. Plot the static pressure across the line surface porous-cl.Plot /XY Plot...(a) Make sure that the Pressure... and Static Pressure are selected from the Y Axis Function drop-down lists.(b) Select porous-cl from the Surfaces selection list.(c) Click Plot and close the Solution XY Plot panel.Figure 7.6: Plot of the Static Pressure on the porous-cl Line SurfaceIn Figure 7.6, the pressure drop across the porous substrate can be seen to be roughly 300 Pa.10. Display filled contours of the velocity in the X direction on the x=95, x=130 and x=165 surfaces.Display /Contours...(a) Disable the Global Range option.(b) Select Velocity... and X Velocity from the Contours of drop-down lists.(c) Select x=130, x=165, and x=95 from the Surfaces selection list, and deselect y=0.(d) Click Display and close the Contours panel.The velocity profile becomes more uniform as the fluid passes through the porous media. The velocity is very high at the center (the area in red) just before the nitrogen enters the substrate and then decreases as it passes through and exits the substrate. The area in green, which corresponds to a moderate velocity, increases in extent.Figure 7.7: Contours of the X Velocity on the x=95, x=130, and x=165 Surfaces11. Use numerical reports to determine the average, minimum, and maximum of the velocity distribution before and after the porous substrate.Report /Surface Integrals...(a) Select Mass-Weighted Average from the Report Type drop-down list.(b) Select Velocity and X Velocity from the Field Variable drop-down lists.(c) Select x=165 and x=95 from the Surfaces selection list.(d) Click Compute.(e) Select Facet Minimum from the Report Type drop-down list and click Compute again.(f) Select Facet Maximum from the Report Type drop-down list and click Compute again.(g) Close the Surface Integrals panel.The numerical report of average, maximum and minimum velocity can be seen in the main FLUENT console, as shown in the following example:The spread between the average, maximum, and minimum values for X velocity gives the degree to which the velocity distribution is non-uniform. You can also use these numbers to calculate the velocity ratio (i.e., the maximum velocity divided by the mean velocity) and the space velocity (i.e., the product of the mean velocity and the substrate length).Custom field functions and UDFs can be also used to calculate more complex measures ofnon-uniformity, such as the standard deviation and the gamma uniformity index.SummaryIn this tutorial, you learned how to set up and solve a problem involving gas flow through porous media in FLUENT. You also learned how to perform appropriate post-processing to investigate the flow field, determine the pressure drop across the porous media and non-uniformity of the velocity distribution as the fluid goes through the porous media.Further ImprovementsThis tutorial guides you through the steps to reach an initial solution. You may be able to obtain a more accurate solution by using an appropriate higher-order discretization scheme and by adapting the grid. Grid adaption can also ensure that the solution is independent of the grid. These steps aredemonstrated in Tutorial 1.。