1、fluent中多孔介质设置问题和算例之欧阳育创编经过痛苦的一段经历，终于将局部问题真相大白，为了使保位同仁不再经过我之痛苦，现在将本人多孔介质经验公布如下，希望各位能加精：时间：2021.02.04创作：欧阳育1。Gambit中划分网格之后，定义需要做为多孔介质的区域为fluid,与缺省的fluid分别开来，再定义其名称，我习惯将名称定义为porous;2。在fluent中定义边界条件define-boundary condition-porous(刚定义的名称)，将其设置边界条件为fluid,点击set按钮即弹出与fluid边界条件一样的对话框，选中porous zone与laminar复选
2、框,再点击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、述；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, suchas filters, cat
4、alyst beds, and packing. This tutorial illustrates how to set up and solve aproblem involving gas flow through porous media.The industrial problem solved here involves gas flow through a catalytic converter. Catalyticconverters are commonly used to purify emissions from gasoline and diesel enginesby
5、 converting environmentally hazardous exhaust emissions to acceptable substances.Examples of such emissions include carbon monoxide (CO), nitrogen oxides (NOx), andunburned hydrocarbon fuels. These exhaust gas emissions are forced through a substrate,which is a ceramic structure coated with a metal
6、catalyst such as platinum or palladium.The nature of the exhaust gas flow is a very important factor in determining the performanceof the catalytic converter. Of particular importance is the pressure gradientand velocity distribution through the substrate. Hence CFD analysis is used to designefficie
7、nt catalytic converters: by modeling the exhaust gas flow, the pressure drop andthe uniformity of flow through the substrate can be determined. In this tutorial, FLUENTis used to model the flow of nitrogen gas through a catalytic converter geometry, so thatthe flow field structure may be analyzed.Th
8、is 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 pressurebasedsolver._ Plot pressure and velocity distribution on specified planes of the geometry._ Det
9、ermine the pressure drop through the substrate and the degree of non-uniformityof 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 inthrough the inlet with a uniform ve
10、locity of 22.6 m/s, passes through a ceramic monolithsubstrate 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 substrateis laminar and is characterized by inertial and viscous loss coefficients in the f
11、low (X)direction. The substrate is impermeable in other directions, which is modeled using losscoefficients 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
12、/CheckFLUENT will perform various checks on the mesh and report the progress in theconsole. 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 w
13、ill 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 Surfacesselection list.(b) Click Display.(c) Rotate the view and zoom in to get the display shown in Figure
14、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 mater
15、ials by copying it from the Fluent Databasefor materials.Define /Materials.(a) Click the Fluent Database. button to open the Fluent Database Materialspanel.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 materi
16、als.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 bou
17、ndary 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 withoutturbulence.(d) Click the Porous Zone tab.i. Make
18、 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 thepanel.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 visi
19、ble 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 dropdownlist in the Turbulence group box.(c) Reta
20、in 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 f
21、rom the Specification Method dropdownlist 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
22、) andclose 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 Discretizationgroup box.(c) Click OK to close the Solu
23、tion 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
24、Monitors to 1.(b) Enable the Plot and Write options for monitor-1, and click the Define. buttonto 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 p
25、anel.(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 /
26、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 thispoint the mass flow rate monitor has attended out, as seen in Figure 7.3.(c) Close the Iter
27、ate 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
28、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
29、Name.(e) Click Create.2. Create cross-sectional surfaces at locations on either side of the substrate, as wellas 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 Is
30、o-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-Valuesof 130 and 165, respectively. Close the Iso-Surface panel after all the surfaceshave been created.3. Create a line surface for the centerline of the porous
31、media.Surface /Line/Rake.(a) Enter the coordinates of the line under End Points, using the starting coordinateof (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. D
32、isplay 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 onlysubstrate-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.
