Tutorial 7. Modeling Flow Through Porous Media Introduction
Many industrial applications involve the modeling of ow through porous media, such as _lters, catalyst beds, and packing. This tutorial illustrates how to set up and solve a problem involving gas ow through porous media.
The industrial problem solved here involves gas ow 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 ow 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 designe_cient catalytic converters: by modeling the exhaust gas ow, the pressure drop and the uniformity of ow through the substrate can be determined. In this tutorial, FLUENT
is used to model the ow of nitrogen gas through a catalytic converter geometry, so that the ow _eld 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 ow through the catalytic converter using the pressurebased solver.
_ Plot pressure and velocity distribution on speci_ed planes of the geometry.
_ Determine the pressure drop through the substrate and the degree of non-uniformity
of ow through cross sections of the geometry using X-Y plots and numerical reports.
许多工业应用都涉及通过多孔介质（如过滤器，催化剂床和填料）的流动模型。本教程说明如何建立和解决涉及气体通过多孔介质的问题。
这里解决的工业问题涉及通过催化转换器的气体流量。催化转化器通常用于通过将对环境有害的废气排放物转化为可接受的物质来净化汽油和柴油发动机的排放物。
这种排放的例子包括一氧化碳（CO），氮氧化物（NOx）和未燃烧的碳氢化合物燃料。这些废气排放物被迫通过衬底，该衬底是涂覆有诸如铂或钯的金属催化剂的陶瓷结构。
排气流量的性质是决定催化转化器性能的一个非常重要的因素。特别重要的是通过基底的压力梯度和速度分布。因此，使用CFD分析来设计催化转换器：通过对排气流量进行建模，可以确定通过基板的流量的压降和流量的均匀性。在本教程中，FLUENT
用于模拟通过催化转化器几何形状的氮气流量，从而可以分析流量结构。
本教程演示了如何执行以下操作：
_设置具有适当阻力的基材的多孔区域。
_使用基于压力的解算器计算通过催化转化器的气体流量的解决方案。
_绘制几何体特定平面上的压力和速度分布。
_确定通过基材的压降和不均匀的程度
通过使用X-Y图和数字报告的几何横截面的流量。
Prerequisites
This tutorial assumes that you are familiar with the menu structure in FLUENT and that
you have completed Tutorial 1. Some steps in the setup and solution procedure will not
be shown explicitly.
本教程假设您熟悉FLUENT中的菜单结构您已完成教程1.设置和解决方案过程中的某些步骤不会明确显示。Problem Description
The catalytic converter modeled here is shown in Figure 7.1. The nitrogen ows 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.
图7.1所示为这里建模的催化转化器。氮流入以22.6m / s的均匀速度通过入口，穿过陶瓷整体具有方形通道的基底，然后通过出口离开。
While the ow in the inlet and outlet sections is turbulent, the ow through the substrate is laminar and is characterized by inertial and viscous loss coe_cients in the ow (X) direction. The substrate is impermeable in other directions, which is modeled using losscoe_cients whose values are three orders of magnitude higher than in the X direction.
当入口和出口段的湍流是湍流时，通过基底的流动是层流的，并且在流动（X）方向上具有惯性和粘性损失系数。基材在其他方向上是不透水的，其使用的损失系数比在X方向上高三个数量级。
Setup and Solution
Preparation
1. Download porous.zip from the Fluent Inc. User Services Center or copy it from
the FLUENT documentation CD to your working folder (as described in Tutorial 1).
2. Unzip porous.zip.
catalytic converter.msh can be found in the porous folder created after unzippingthe _le.
3. Start the 3D (3d) version of FLUENT.
1.从Fluent Inc.用户服务中心下载porous.zip或从中复制
FLUENT文档光盘放到您的工作文件夹中（如教程1所述）。
2.解压多孔.zip。
催化转换器.msh可以在解压后形成的多孔文件夹中找到。
3.启动FLUENT的3D（3D）版本。
Step 1: Grid
1. Read the mesh _le (catalytic converter.msh).
File ! Read !Case...
2. Check the grid.
Grid !Check
FLUENT 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: Models
1. Retain the default solver settings. De_ne ! Models !Solver...
2. Select the standard k-_ turbulence model.
De_ne ! Models !Viscous...
Step 3: Materials
1. Add nitrogen to the list of uid materials by copying it from the Fluent Database for materials. De_ne !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 uid materials. iii. Close the Fluent Database Materials panel.
(b) Close the Materials panel.
Step 4: Boundary Conditions De_ne !Boundary Conditions...
1. Set the boundary conditions for the uid (uid).
(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 ow 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 _elds 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 _elds 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 Speci_cation 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 Speci_cation Method dropdown list in the Turbulence group box.
(c) Enter 5% for the Backow Turbulent Intensity.
(d) Enter 42 mm for the Backow 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: Solution
1. 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 ow 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 De_ne... button to open the De_ne 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 De_ne 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 _le (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 ow rate monitor has attened out, as seen in Figure 7.3.
(c) Close the Iterate panel.
7. Save the case and data _les (catalytic converter.cas and catalytic converter.dat).
File ! Write !Case & Data...
Note: If you choose a _le name that already exists in the current folder, FLUENT
will prompt you for con_rmation to overwrite the _le.
Step 6: Postprocessing
1. Create a surface passing through the centerline for postprocessing 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 ow pattern shows that the ow enters the catalytic converter as a jet, withrecirculation on either side of the jet. As it passes through the porous substrate, itdecelerates 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 Plane
8. Display _lled 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 plane
The pressure changes rapidly in the middle section, where the uid 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 Surface
In Figure 7.6, the pressure drop across the porous substrate can be seen to be roughly 300 Pa.
10. Display _lled 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 deselecty=0.
(d) Click Display and close the Contours panel.
The velocity pro_le becomes more uniform as the uid 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
thesubstrate. 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 Surfaces
11. Use numerical reports to determine the average, minimum, and maximum of thevelocity 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 inthe 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 _eld functions and UDFs can be also used to calculate more complex measures of non-uniformity, such as the
standard deviation and the gamma uniformityindex.
Summary
In this tutorial, you learned how to set up and solve a problem involving gas ow through porous media in FLUENT. You also learned how to perform appropriate postprocessing to investigate the ow _eld, determine the pressure drop across the porous media and non-uniformity of the velocity distribution as the uid goes through the porous media.See Section 7.19 of the User's Guide for additional details about modeling ow through porous media (including heat transfer and reaction modeling).
Further Improvements
This 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 are demonstrated in Tutorial 1.
Tutorial 8. Using a Single Rotating Reference Frame Introduction
This tutorial considers the ow within a 2D, axisymmetric, co-rotating disk cavity system.
Understanding the behavior of such ows is important in the design of secondary air
passages for turbine disk cooling.
This tutorial demonstrates how to do the following:
_ Set up a 2D axisymmetric model with swirl, using a rotating reference frame.
_ Use the standard k-_ and RNG k-_ turbulence models with the enhanced near-wall
treatment.
_ Calculate a solution using the pressure-based solver.
_ Display velocity vectors and contours of pressure.
_ Set up and display XY plots of radial velocity.
_ Restart the solver from an existing solution.
Prerequisites
This tutorial assumes that you are familiar with the menu structure in FLUENT and that
you have completed Tutorial 1. Some steps in the setup and solution procedure will not
be shown explicitly.
Problem Description
The problem to be considered is shown schematically in Figure 8.1. This case is similar
to a disk cavity con_guration that was extensively studied by Pincombe [1].
Air enters the cavity between two co-rotating disks. The disks are 88.6 cm in diameter
and the air enters at 1.146 m/s through a circular bore 8.86 cm in diameter. The disks,
which are 6.2 cm apart, are spinning at 71.08 rpm, and the air enters with no swirl. As
the ow is diverted radially, the rotation of the disk has a signi_cant e_ect on the viscous
ow developing along the surface of the disk.
2004-06 FLUENT流体工程仿真计算实例与应用韩占忠王敬兰小平北京理工大学出版社 第一章流体力学基础与fluent简介 第二章二维流动与传热的数值计算 第一节冷、热水混合器内部二维流动 第二节喷管内二维非定常流动 第三节三角翼的可压缩外部绕流 第四节三角翼不可压缩的外部绕流(空化模型应用) 第五节vof模型的应用 第六节组分传输与气体燃烧 第三章三维流动与传热的数值计算 第一节冷、热水混合器内的三维流动与换热 第二节粘性流体通过圆管弯头段的三维流动 第三节三维稳态热传导问题 第四节动网格问题 第五节叶轮机械的mixing plane模型 2004-09 计算流体动力学分析CFD软件原理与应用王福军清华大学出版社（偏重理论） 第1章计算流动力学基础知识 第2章基于有限体积法的控制方程离散 第3章基于SIMPLE算法的流场数值计算
第4章三维流模型及其在CFD中的应用 第5章边界条件的应用 第6章网格的生成 第7章FLUENT软件的基本用法 第8章CFD综合应用实例 2007-02 FLUENT技术基础与应用实例王瑞金张凯王刚清华大学出版社 第1章Fluent概述 第2章流体力学基础知识 第3章流体力学数值模拟基础 第4章Fluent软件介绍 第5章速度场的计算 第6章温度场的计算 第7章多相流模型 第8章凝固和融化模型 第9章可动区域中流动问题的模拟 第10章动网格模型 第11章UDF和UDS 第12章Fluent并行计算 第13章Tecplot软件 2008-07 Fluent高级应用与实例分析江帆，黄鹏清华大学出版社第1章 CFD基础 第2章Fluent基本介绍
FLUENT 教程 赵玉新 I、目录 第一章、开始 第二章、操作界面 第三章、文件的读写 第四章、单位系统 第五章、读入和操作网格 第六章、边界条件 第七章、物理特性 第八章、基本物理模型 第九章、湍流模型 第十章、辐射模型 第十一章、化学输运与反应流 第十二章、污染形成模型 第十三章、相变模拟 第十四章、多相流模型 第十五章、动坐标系下的流动 第十六章、解算器的使用 第十七章、网格适应 第十八章、数据显示与报告界面的产生 第十九章、图形与可视化 第二十章、Alphanumeric Reporting 第二十一章、流场函数定义 第二十二章、并行处理 第二十三章、自定义函数 第二十四章、参考向导 第二十五章、索引（Bibliography） 第二十六章、命令索引 II、如何使用该教程 概述 本教程主要介绍了FLUENT 的使用，其中附带了相关的算例，从而能够使每一位使用 者在学习的同时积累相关的经验。本教程大致分以下四个部分：第一部分包括介绍信息、用户界面信息、文件输入输出、单位系统、网格、边界条件以及物理特性。第二和第三部分包含物理模型，解以及网格适应的信息。第四部分包括界面的生成、后处理、图形报告、并行处理、自定义函数以及FLUENT 所使用的流场函数与变量的定义。 下面是各章的简略概括 第一部分： z开始使用：本章描述了FLUENT 的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式，并且概述了问题解决的大致步骤。在本章中，我们给出
了一个可以在你自己计算机上运行的简单的算例。 z使用界面：本章描述了用户界面、文本界面以及在线帮助的使用方法。同时也提供了远程处理与批处理的一些方法。（请参考关于特定的文本界面命令的在线帮助） z读写文件：本章描述了FLUENT 可以读写的文件以及硬拷贝文件。 z单位系统：本章描述了如何使用FLUENT 所提供的标准与自定义单位系统。 z读和操纵网格：本章描述了各种各样的计算网格来源，并解释了如何获取关于网格的诊断信息，以及通过尺度化（scale）、分区（partition）等方法对网格的修改。本章还描述了非一致（nonconformal）网格的使用. z边界条件：本章描述了FLUENT 所提供的各种类型边界条件，如何使用它们，如何定义它们and how to define boundary profiles and volumetric sources. z物理特性：本章描述了如何定义流体的物理特性与方程。FLUENT 采用这些信息来处理你的输入信息。 第二部分： z基本物理模型：本章描述了FLUENT 计算流体流动和热传导所使用的物理模型（包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流）。以及在使用这些模型时你需要输入的数据，本章也包含了自定义标量的信息。 z湍流模型：本章描述了FLUENT 的湍流模型以及使用条件。 z辐射模型：本章描述了FLUENT 的热辐射模型以及使用条件。 z化学组分输运和反应流：本章描述了化学组分输运和反应流的模型及其使用方法。本章详细的叙述了prePDF 的使用方法。 z污染形成模型：本章描述了NOx 和烟尘的形成的模型，以及这些模型的使用方法。 第三部分： z相变模拟：本章描述了FLUENT 的相变模型及其使用方法。 z离散相变模型：本章描述了FLUENT 的离散相变模型及其使用方法。 z多相流模型：本章描述了FLUENT 的多相流模型及其使用方法。 z Flows in Moving Zones（移动坐标系下的流动）：本章描述了FLUENT 中单一旋转坐标系，多重移动坐标系，以及滑动网格的使用方法。 z Solver 的使用：本章描述了如何使用FLUENT 的解法器（solver）。 z网格适应：本章描述了explains the solution-adaptive mesh refinement feature in FLUENT and how to use it 第四部分： z显示和报告数据界面的创建：本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data z图形和可视化：本章描述了检验FLUENT 解的图形工具 z Alphanumeric Reporting：本章描述了如何获取流动、力、表面积分以及其它解的数据。 z流场函数的定义：本章描述了如何定义FLUENT 面板内出现的变量选择下拉菜单中的流动变量，并且告诉我们如何创建自己的自定义流场函数。 z并行处理：本章描述了FLUENT 的并行处理特点以及使用方法 z自定义函数：本章描述了如何通过用户定义边界条件，物理性质函数来形成自己的FLUENT 软件。 如何使用该手册 z根据你对CFD 以及FLUENT 公司的熟悉，你可以通过各种途径使用该手册 对于初学者，建议如下：
Tutorial 7. Modeling Flow Through Porous Media Introduction Many industrial applications involve the modeling of ow through porous media, such as _lters, catalyst beds, and packing. This tutorial illustrates how to set up and solve a problem involving gas ow through porous media. The industrial problem solved here involves gas ow 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 ow 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 designe_cient catalytic converters: by modeling the exhaust gas ow, the pressure drop and the uniformity of ow through the substrate can be determined. In this tutorial, FLUENT is used to model the ow of nitrogen gas through a catalytic converter geometry, so that the ow _eld 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 ow through the catalytic converter using the pressurebased solver. _ Plot pressure and velocity distribution on speci_ed planes of the geometry. _ Determine the pressure drop through the substrate and the degree of non-uniformity of ow through cross sections of the geometry using X-Y plots and numerical reports. 许多工业应用都涉及通过多孔介质（如过滤器，催化剂床和填料）的流动模型。本教程说明如何建立和解决涉及气体通过多孔介质的问题。 这里解决的工业问题涉及通过催化转换器的气体流量。催化转化器通常用于通过将对环境有害的废气排放物转化为可接受的物质来净化汽油和柴油发动机的排放物。 这种排放的例子包括一氧化碳（CO），氮氧化物（NOx）和未燃烧的碳氢化合物燃料。这些废气排放物被迫通过衬底，该衬底是涂覆有诸如铂或钯的金属催化剂的瓷结构。 排气流量的性质是决定催化转化器性能的一个非常重要的因素。特别重要的是通过基底的压力梯度和速度分布。因此，使用CFD分析来设计催化转换器：通过对排气流量进行建模，可以确定通过基板的流量的压降和流量的均匀性。在本教程中，FLUENT 用于模拟通过催化转化器几何形状的氮气流量，从而可以分析流量结构。 本教程演示了如何执行以下操作： _设置具有适当阻力的基材的多孔区域。 _使用基于压力的解算器计算通过催化转化器的气体流量的解决方案。 _绘制几何体特定平面上的压力和速度分布。 _确定通过基材的压降和不均匀的程度 通过使用X-Y图和数字报告的几何横截面的流量。 Prerequisites This tutorial assumes that you are familiar with the menu structure in FLUENT and that you have completed Tutorial 1. Some steps in the setup and solution procedure will not be shown explicitly. 本教程假设您熟悉FLUENT中的菜单结构您已完成教程1.设置和解决方案过程中的某些步骤不会明确显示。
1.pressure based 和density based Coupled会同时求解所有的方程（质量守恒方程、动量守恒方程和能量守恒方程）而不是单个方程求解（方程互相分离）。当速度和压力高度耦合（高压和高速）时应该使用耦合求解，但这样会需要较长的计算时间。 在耦合求解中，能量方程中总是包含组分扩散（Species Diffusion Term）项。 当使用segregated求解时，fluent允许指定固体材料的各项异性传导性。 求解方法主要根据要求解的模型来选择。Segregated方法是基于压力，而coupled求解是基于密度的。这样就使得segregated求解低速流动较好而coupled求解音速/超音速问题较好。我不推荐使用coupled求解所有低于马赫数4的流动（直到基于压力的coupled求解方法出现在下一个fluent版本中）。我曾经用segregated方法求解直到1.5马赫的问题，并且结果很好。但是速度越高，需要的网格就越多（因为segregated趋向于“平滑”波动），所以必须多加注意划分网格。 Coupled方法使用默认设置时往往是比较稳定的。Segregated方法常常对容许极根很敏感。当使用segregated方法求解时，不要提高turbulent viscosity ration limit（除非你根据过去的经验或者你的物理模型有很好的理由超过这个极限，但我从没有听说这样是比较理想的）。不要给压力和温度极限限定的合理的范围（例如Plimits=Pstatic+/-(2*dynamic pressure)）来计算适当的温度。 1，pressure based 求解方法在求解不可压流体时，如果我们联立求解从动量方程和连续性方程离散得到的代数方程组，可以直接得到各速度分量及相应的压力值，但是要占用大量的计算内存，这一方法已可以在Fluent6.3中实现，所需内存为分离算法的1.5-2倍，同时Fluent6.3中的压力基耦合求解器也很适合求解带有激波的高速空气动力问题（可压流体），这是一个新变化。本人也在尝试用这个模型模拟一些噪声问题。 2，density based求解方法是针对可压流体设计的，因而更适合于可压流场的计算。以速度分量、密度（密度基）作为基本变量，压力则由状态方程求解。 Pressure-Based Solver是Fluent的优势，它是基于压力法的求解器，使用的是压力修正算法，求解的控制方程是标量形式的，擅长求解不可压缩流动，对于可压流动也可以求解；Fluent 6.3以前的版本求解器，只有Segregated Solver和Coupled Solver，其实也就是Pressure-Based Solver 的两种处理方法； Density-Based Solver是Fluent 6.3新发展出来的，它是基于密度法的求解器，求解的控制方程是矢量形式的，主要离散格式有Roe，AUSM+，该方法的初衷是让Fluent具有比较好的求解可压缩流动能力，但目前格式没有添加任何限制器，因此还不太完善；它只有Coupled的算法；对于低速问题，他们是使用Preconditioning方法来处理，使之也能够计算低速问题。 Density-Based Solver下肯定是没有SIMPLEC，PISO这些选项的，因为这些都是压力修正算法，不会在这种类型的求解器中出现的；一般还是使用Pressure-Based Solver解决问题。 2.连续性方程不收敛是怎么回事？ 在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事 这和Fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。 你可以试验SIMPLEC方法，应该会收敛快些。 3.边界条件对应的一般设定方法 边界条件对应的一般设定方法：
Fluent经典问题及解答 1 对于刚接触到FLUENT新手来说，面对铺天盖地的学习资料和令人难读的FLUENT help，如何学习才能在最短的时间内入门并掌握基本学习方法呢？(#61) 2 CFD计算中涉及到的流体及流动的基本概念和术语：理想流体和粘性流体；牛顿流体和非牛顿流体；可压缩流体和不可压缩流体；层流和湍流；定常流动和非定常流动；亚音速与超音速流动；热传导和扩散等。（13楼） 3 在数值模拟过程中，离散化的目的是什么？如何对计算区域进行离散化？离散化时通常使用哪些网格？如何对控制方程进行离散？离散化常用的方法有哪些？它们有什么不同？（#80） 4 常见离散格式的性能的对比（稳定性、精度和经济性）(#62) 5 在利用有限体积法建立离散方程时，必须遵守哪几个基本原则？（#81） 6 流场数值计算的目的是什么？主要方法有哪些？其基本思路是什么？各自的适用范围是什么？(#130) 7 可压缩流动和不可压缩流动，在数值解法上各有何特点？为何不可压缩流动在求解时反而比可压缩流动有更多的困难？（#55） 8 什么叫边界条件？有何物理意义？它与初始条件有什么关系？（#56） 9 在一个物理问题的多个边界上，如何协调各边界上的不同边界条件？在边界条件的组合问题上，有什么原则？ 10 在数值计算中，偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别？（#143） 11 在网格生成技术中，什么叫贴体坐标系？什么叫网格独立解？（#35） 12 在GAMBIT的foreground和background中，真实体和虚实体、实操作和虚操作四个之间是什么关系？ 13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量？及其在做网格时大致注意到哪些细节？（#38） 14 画网格时，网格类型和网格方法如何配合使用？各种方法有什么样的应用范围及做网格时需注意的问题？(#169) 15 对于自己的模型，大多数人有这样的想法：我的模型如何来画网格？用什么样的方法最简单？这样做网格到底对不对？(#154) 16 在两个面的交界线上如果出现网格间距不同的情况时，即两块网格不连续时，怎么样克服这种情况呢？（#40） 17 依据实体在GAMBIT建模之前简化时，必须遵循哪几个原则？(#170) 18 在设置GAMBIT边界层类型时需要注意的几个问题：a、没有定义的边界线如何处理？b、计算域内的内部边界如何处理（2D）？（#128） 19 为何在划分网格后，还要指定边界类型和区域类型？常用的边界类型和区域类型有哪些？（#127） 20 何为流体区域（fluid zone）和固体区域（solid zone）？为什么要使用区域的概念？FLUENT是怎样使用区域的？(#41) 21 如何监视FLUENT的计算结果？如何判断计算是否收敛？在FLUENT中收敛准则是如何定义的？分析计算收敛性的各控制参数，并说明如何选择和设置这些参数？解决不收敛问题通常的几个解决方法是什么？（9楼） 22 什么叫松弛因子？松弛因子对计算结果有什么样的影响？它对计算的收敛情况又有什么样的影响？（7楼）
QUICK格式可能产生比二阶精度更好的结果。但是，一般情况下，用二阶精度就已足够，即使使用QUICK格式，结果也不一定好。乘方格式(Power-law Scheme)一般产生与一阶精度格式相同精度的结果。中心差分格式一般只用于大涡模拟，而且要求网格很细的情况。 53 对于FLUENT的耦合解算器，对时间步进格式的主要控制是Courant数（CFL），那么Courant 数对计算结果有何影响？ courant number实际上是指时间步长和空间步长的相对关系，系统自动减小courant数，这种情况一般出现在存在尖锐外形的计算域，当局部的流速过大或者压差过大时出错，把局部的网格加密再试一下。 在Fluent中，用courant number来调节计算的稳定性与收敛性。一般来说，随着courant number的从小到大的变化，收敛速度逐渐加快，但是稳定性逐渐降低。所以具体的问题，在计算的过程中，最好是把courant number从小开始设置，看看迭代残差的收敛情况，如果收敛速度较慢而且比较稳定的话，可以适当的增加courant number的大小，根据自己具体的问题，找出一个比较合适的courant number，让收敛速度能够足够的快，而且能够保持它的稳定性。 54 在分离求解器中，FLUENT提供了压力速度耦和的三种方法：SIMPLE，SIMPLEC及PISO，它们的应用有什么不同？ 在FLUENT中，可以使用标准SIMPLE算法和SIMPLEC（SIMPLE-Consistent）算法，默认是SIMPLE算法，但是对于许多问题如果使用SIMPLEC可能会得到更好的结果，尤其是可以应用增加的亚松驰迭代时，具体介绍如下： 对于相对简单的问题（如：没有附加模型激活的层流流动），其收敛性已经被压力速度耦合所限制，你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中，压力校正亚松驰因子通常设为1.0，它有助于收敛。但是，在有些问题中，将压力校正松弛因子增加到1.0可能会导致不稳定。 对于所有的过渡流动计算，强烈推荐使用PISO算法邻近校正。它允许你使用大的时间步，而且对于动量和压力都可以使用亚松驰因子1.0。对于定常状态问题，具有邻近校正的PISO 并不会比具有较好的亚松驰因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。 当你使用PISO邻近校正时，对所有方程都推荐使用亚松驰因子为1.0或者接近1.0。如果你只对高度扭曲的网格使用PISO倾斜校正，请设定动量和压力的亚松驰因子之和为1.0比如：压力亚松驰因子0.3，动量亚松驰因子0.7）。如果你同时使用PISO的两种校正方法，推荐参阅PISO邻近校正中所用的方法 55 对于大多数情况，在选择选择压力插值格式时，标准格式已经足够了，但是对于特定的某些模型使用其它格式有什么特别的要求？ 压力插值方式的列表只在使用Pressure-based求解器中出现。一般情况下可选择Standard；对于含有高回旋数的流动，高 Rayleigh数的自然对流，高速旋转流动，多孔介质流动，高曲率计算区域等流动情况，选择PRESTO格式；对于可压缩流动，选择Second Order；当然也可以选择Second Order以提高精度；对于含有大体力的流动，选择Body Force Weighted。 注意：Second Order格式不可以用于多孔介质；在使用VOF和Mixture多相流模型时，只能
FLUENT中文手册（简化版） 本手册介绍FLUENT的使用方法，并附带了相关的算例。下面是本教程各部分各章节的简略概括。 第一部分： ?开始使用：描述了FLUENT的计算能力以及它与其它程序的接口。介绍了如何对具体的应用选择适当的解形式，并且概述了问题解决的大致步骤。在本章中给出了一个简单的算例。 ?使用界面：描述用户界面、文本界面以及在线帮助的使用方法，还有远程处理与批处理的一些方法。?读写文件：描述了FLUENT可以读写的文件以及硬拷贝文件。 ?单位系统：描述了如何使用FLUENT所提供的标准与自定义单位系统。 ?使用网格：描述了各种计算网格来源，并解释了如何获取关于网格的诊断信息，以及通过尺度化（scale）、分区（partition）等方法对网格的修改。还描述了非一致（nonconformal）网格的使用. ?边界条件：描述了FLUENT所提供的各种类型边界条件和源项，如何使用它们，如何定义它们等 ?物理特性：描述了如何定义流体的物理特性与方程。FLUENT采用这些信息来处理你的输入信息。 第二部分： ?基本物理模型：描述了计算流动和传热所用的物理模型（包括自然对流、周期流、热传导、swirling、旋转流、可压流、无粘流以及时间相关流）及其使用方法，还有自定义标量的信息。 ?湍流模型：描述了FLUENT的湍流模型以及使用条件。 ?辐射模型：描述了FLUENT的热辐射模型以及使用条件。 ?化学组分输运和反应流：描述了化学组分输运和反应流的模型及其使用方法，并详细叙述了prePDF 的使用方法。 ?污染形成模型：描述了NOx和烟尘的形成的模型，以及这些模型的使用方法。 第三部分： ?相变模拟：描述了FLUENT的相变模型及其使用方法。 ?离散相变模型：描述了FLUENT的离散相变模型及其使用方法。 ?多相流模型：描述了FLUENT的多相流模型及其使用方法。 ?移动坐标系下的流动：描述单一旋转坐标系、多重移动坐标系、以及滑动网格的使用方法。 ?解法器（solver）的使用：描述了如何使用FLUENT的解法器。 ?网格适应：描述了如何优化网格以适应计算需求。 第四部分： ?显示和报告数据界面的创建：本章描述了explains how to create surfaces in the domain on which you can examine FLUENT solution data ?图形和可视化：本章描述了检验FLUENT解的图形工具 ?Alphanumeric Reporting：本章描述了如何获取流动、力、表面积分以及其它解的数据。 ?流场函数的定义：本章描述了如何定义FLUENT面板内出现的变量选择下拉菜单中的流动变量，并且告诉我们如何创建自己的自定义流场函数。 ?并行处理：本章描述了FLUENT的并行处理特点以及使用方法 ?自定义函数：本章描述了如何通过用户定义边界条件，物理性质函数来形成自己的FLUENT软件。 如何使用该手册 对于初学者，建议从阅读“开始”这一章起步。 对于有经验的使用者，有三种不同的方法供你使用该手册：按照特定程序的步骤从按程序顺序排列的目录列表和主题列表中查找相关资料；从命令索引查找特定的面板和文本命令的使用方法；从分类索引查找特定类别信息（在线帮助中没有此类索引，只能在印刷手册中找到它）。 什么时候使用Support Engineer：Support Engineer能帮你计划CFD模拟工程并解决在使用FLUENT 中所遇到的困难。在遇到困难时我们建议你使用Support Engineer。但是在使用之前有以下几个注意事项：●仔细阅读手册中关于你使用并产生问题的命令的信息 ●回忆导致你产生问题的每一步 ●如果可能的话，请记下所出现的错误信息 ●对于特别困难的问题，保存FLUENT出现问题时的日志以及手稿。在解决问题时，它是最好的资源。
计算流体力学作业FLUENT 模拟燃烧 问题描述：长为2m、直径为0.45m的圆筒形燃烧器结构如图1所示，燃烧筒壁上嵌有三块厚为0.0005 m，高0.05 m的薄板，以利于甲烷与空气的混合。燃烧火焰为湍流扩散火焰。在燃烧器中心有一个直径为0.01 m、长为0.01 m、壁厚为0.002 m的小喷嘴，甲烷以60 m/s的速度从小喷嘴注入燃烧器。空气从喷嘴周围以0.5 m/s的速度进入燃烧器。总当量比大约是0.76(甲烷含量超过空气约28%)，甲烷气体在燃烧器中高速流动，并与低速流动的空气混合，基于甲烷喷嘴直径的雷诺数约为5.7×103。 假定燃料完全燃烧并转换为:CH4+2O2→CO2+2H2O 反应过程是通过化学计量系数、形成焓和控制化学反应率的相应参数来定义的。利用FLUENT的finite-rate化学反应模型对一个圆筒形燃烧器内的甲烷和空气的混合物的流动和燃烧过程进行研究。 1、建立物理模型，选择材料属性，定义带化学组分混合与反应的湍流流动边界条件 2、使用非耦合求解器求解燃烧问题 3、对燃烧组分的比热分别为常量和变量的情况进行计算，并比较其结果 4、利用分布云图检查反应流的计算结果 5、预测热力型和快速型的NO X含量 6、使用场函数计算器进行NO含量计算 一、利用GAMBIT建立计算模型 第1步启动GAMBIT，建立基本结构 分析：圆筒燃烧器是一个轴对称的结构，可简化为二维流动，故只要建立轴对称面上的
二维结构就可以了，几何结构如图2所示。 (1)建立新文件夹 在F盘根目录下建立一个名为combustion的文件夹。 (2)启动GAMBIT (3)创建对称轴 ①创建两端点。A(0,0,0)，B(2,0,0) ②将两端点连成线 (4)创建小喷嘴及空气进口边界 ①创建C、D、E、F、G点
网格质量与那些因素有关？ 网格质量本身与具体问题的具体几何特性、流动特性及流场求解算法有关。因此，网格质量最终要由计算结果来评判，但是误差分析以及经验表明，CFD计算对计算网格有一些一般性的要求，例如光滑性、正交性、网格单元的正则性以及在流动变化剧烈的区域分布足够多的网格点等。对于复杂几何外形的网格生成，这些要求往往并不可能同时完全满足。例如，给定边界网格点分布，采用Laplace 方程生成的网格是最光滑的，但是最光滑的网格不一定满足物面边界正交性条件，其网格点分布也很有可能不能捕捉流动特征，因此，最光滑的网格不一定是最好的网格。对计算网格的一个最基本的要求当然是所有网格点的Jacobian必须为正值，即网格体积必须为正，其他一些最常用的网格质量度量参数包括扭角（skew angle）、纵横比（aspect ratio、Laplacian）、以及弧长（arc length）等。通过计算、检查这些参数，可以定性的甚至从某种程度上定量的对网格质量进行评判。Parmley等给出了更多的基于网格元素和网格节点的网格质量度量参数。有限元素法关于插值逼近误差估计的理论，实际上也对网格单元的品质给出了基本的规定：即每个单元的内切球半径与外切球半径之，应该是一个适当的，与网格疏密无关的常数。 实体与虚体的区别 在建模中，经常会遇到实体、实面与虚体、虚面，虚体的计算域也可以进行计算并得到所需的结果。那么它们的区别是什么呢？ 对于求解是没有任何区别的，只要你能在虚体或者实体上划分你需要的网格。关键是看你网格生成的质量如何，与实体虚体无关。 gambit的实体和虚体在生成网格和计算的时候对于结果没有任何影响，实体和虚体的主要区别有以下几点： 1.实体可以进行布尔运算但是虚体不能，虽然不能进行布尔运算，但是虚体存在merge，split 等功能。 2.实体运算在很多cad软件里面都有，但是虚体是gambit的一大特色，有了虚体以后，gambit 的建模和网格生成的灵活性增加了很多。 3.在网格生成的过程中，如果有几个相对比较平坦的面，你可以把它们通过merge合成一个，这样，作网格的时候，可以节省步骤，对于曲率比较大的面，可能生成的网格质量不好，这时候，你可以采取用split的方式把它划分成几个小面以提高网格质量。 在Fluent中进行非稳态（unsteady）计算时如何设置步长？
湍流与黏性有什么关系？ 湍流和粘性都是客观存在的流动性质。 湍流的形成需要一定的条件，粘性是一切流动都具有的。 流体流动方程本身就是具非线性的。 NS方程中的粘性项就是非线性项，当然无粘的欧拉方程也是非线性的。 粘性是分子无规则运动引起的，湍流相对于层流的特性是由涡体混掺运动引起的。 湍流粘性是基于湍流体的parcel湍流混掺是类比于层流体中的分子无规则运动，只是分子无规则运动遥远弱些吧了。不过，这只是类比于，要注意他们可是具有不同的属性。 粘性是耗散的根源，实际流体总是有耗散的。 而粘性是制约湍流的。 LANDAU说，粘性的存在制约了湍流的自由度。 湍流粘性系数和层流的是不一样的,层流的粘性系数基本可认为是常数,可湍流中层流底层中粘性系数很小,远小于层流时的粘性系数;而在过渡区,与之相当,在一个数量级;在充分发展的湍流区,又远大于层流时的粘性系数.这是鮑辛内斯克1987年提出的。 1 FLUENT的初始化面板中有一项是设置从哪个地方开始计算(compute from),选择从不同的边界开始计算有很大的区别吗?该怎样根据具体问题选择从哪里计算呢?比如有两个速度入口A和B,还有压力出口等等，是选速度入口还是压力出口?如果选速度入口,有两个,该选哪个呀?有没有什么原则标准之类的东西？ 一般是选取ALL ZONE，即所有区域的平均处理，通常也可选择有代表性的进口(如多个进口时)进行初始化。对于一般流动问题，初始值的设定并不重要，因为计算容易收敛。但当几何条件复杂，而且流动速度高变化快(如音速流动)，初始条件要仔细选择。如果不收敛，还应试验不同的初始条件，甚至逐次改变边界条件最后达到所要求的条件。 2 要判断自己模拟的结果是否是正确的，似乎解的收敛性要比那些初始条件和边界条件更重要，可以这样理解吗？也就是说，对于一个具体的问题，初始条件和边界条件的设定并不是唯一的，为了使解收敛，需要不断调整初始条件和边界条件直到解收敛为止，是吗？如果解收敛了，是不是就可以基本确定模拟的结果是正确的呢？ 对于一个具体的问题，边界条件的设定当然是唯一的，只不过初始化时可以选择不同的初始条件（指定常流），为了使解的收敛比较好，我一般是逐渐的调节边界条件到额定值（"额定值"是指你题目中要求的入口或出口条件，例如计算一个管内流动，要求入口压力和温度为10MPa和3000K，那么我开始叠代时选择入口压力和温度为1MPa和500K（假设，这看你自己问题了），等流场计算的初具规模、收敛的较好了，再逐渐调高压力和温度，经过好几次调节后最终到达额定值10MPa和3000K，这样比一开始就设为10MPa 和3000K收敛的要好些）这样每次叠代可以比较容易收敛，每次调节后不用再初始化即自动调用上次的解为这次的初始解，然后继续叠代。即使解收敛了，这并不意味着就可以基本确定模拟的结果是正确的，还需要和实验的结果以及理论分析结果进行对比分析。 连续性方程不收敛是怎么回事？ 在计算过程中其它指数都收敛了，就continuity不收敛是怎么回事 这和Fluent程序的求解方法SIMPLE有关。SIMPLE根据连续方程推导出压力修正方法求解压力。由于连续方程中流场耦合项被过渡简化，使得压力修正方程不能准确反映流场的变化，从而导致该方程收敛缓慢。你可以试验SIMPLEC方法，应该会收敛快些。 边界条件对应的一般设定方法 边界条件对应的一般设定方法： *Genaeral--- pressure inlet;pressure outlet *Compressible flows---mass flow inlet;pressure far-field *Incompressible ---velocity inlet;outflow
Fluent动网格----layering个一个简单实例我这几天看了点动网格技术方面的东西，在学习过程中发现这方面的例子很少，自己也走了一些弯路。现在还好，弄明白了一些，能够应付现在我的工作。为了让更多学习者快速了解动网格，我打算尽量把我学习心得在这里和大家分享，这里给出一个layering的一个简单例子。 1.Gambit画网格 本例很简单，在Gambit里画一个10*10的矩形，网格间隔为1，也就是有100个网格，具体见下图。都学动网格的人了，不至于这个不会做！ 这里需要注意一个问题：设置边界条件的时候，一定要把要移动的边单独设定，本例中一右边界作为移动的边，设成wall就可以，这里再后面需要制定。 2.编写UDF #include "udf.h" #include "unsteady.h" #include "stdio.h" #include "stdlib.h" /************************************************************/ real current_time = 0.0 ; Domain * domain ; Thread * thread ; real NV_VEC( origin ),NV_VEC( force ),NV_VEC( moment ) ; /************************************************************/ DEFINE_CG_MOTION(throttle,dt,vel,omega,time,dtime) { current_time = CURRENT_TIME ; vel[0] = 30; Message("time=%f omega=%f\n",current_time) ; }
fluent 经典问题请问双CPU并行计算的效率问题.txt27信念的力量在于即使身处逆境，亦能帮助你鼓起前进的船帆；信念的魅力在于即使遇到险运，亦能召唤你鼓起生活的勇气；信念的伟大在于即使遭遇不幸，亦能促使你保持崇高的心灵。发信人: rao (绕绕), 信区: NumComp 标题: [合集] 请问双CPU并行计算的效率问题 发信站: BBS 水木清华站 (Mon Jul 7 03:32:43 2003), 站内 ☆─────────────────────────────────────☆ xuzheng (天使暂时离开@_@反方向的钟) 于 (Fri Jul 4 11:03:44 2003) 提到: 大致上只有一个CPU在工作，或者两CPU占有率相当于一个CPU mpich1.2.5+fortran 怎么配置可以使两个CPU同时工作？？ BOW ☆─────────────────────────────────────☆ luxz (panda--在热死和冻死边缘挣扎) 于 (Fri Jul 4 11:04:57 2003) 提到: mpirun -np 2 *.exe 【在 xuzheng (天使暂时离开@_@反方向的钟) 的大作中提到: 】 : 大致上只有一个CPU在工作，或者两CPU占有率相当于一个CPU : mpich1.2.5+fortran : 怎么配置可以使两个CPU同时工作？？ : BOW ☆─────────────────────────────────────☆ xuzheng (天使暂时离开@_@反方向的钟) 于 (Fri Jul 4 11:06:27 2003) 提到: 不是，你误解了我的意思 再具体点说就是16个节点双CPU的集群，并行计算过程中 每个节点的CPU效率大概只有50％
Introduction:This tutorial demonstrates how to model the2D turbu-lent?ow across a circular cylinder using LES(Large Eddy Simula-tion),and compute?ow-induced noise(aero-noise)using FLUENT’s acoustics model. In this tutorial you will learn how to: ?Perform2D Large Eddy Simulation(LES) ?Set parameters for an aero-noise calculation ?Save surface pressure data for an aero-noise calculation ?Calculate aero-noise quantities ?Postprocess an aero-noise solution Prerequisites:This tutorial assumes that you are familiar with the menu structure in FLUENT,and that you have solved or read Tu-torial1.Some steps in the setup and solution procedure will not be shown explicitly. Problem Description:The problem considers turbulent air?ow over a2D circular cylinder at a free stream velocity U of69.19m/s. The cylinder diameter D is1.9cm.The Reynolds number based on the?ow parameters is about90000.The computational do-main(Figure3.0.1)extends5D upstream and20D downstream of the cylinder,and5D on both sides of it.If the computational domain is not taken wide enough on the downstream side,so that no reversed?ow occurs,the accuracy of the aero-noise prediction may be a?ected.The rule of thumb is to take at least20D on the downstream side of the obstacle. c Fluent Inc.June20,20023-1
Fluent经典问题及答疑1 1 对于刚接触到FLUENT新手来说，面对铺天盖地的学习资料和令人难读的FLUENT help，如何学习才能在最短的时间内入门并掌握基本学习方法呢？(#61) 2 CFD计算中涉及到的流体及流动的基本概念和术语：理想流体和粘性流体；牛顿流体和非牛顿流体；可压缩流体和不可压缩流体；层流和湍流；定常流动和非定常流动；亚音速与超音速流动；热传导和扩散等。（13楼） 3 在数值模拟过程中，离散化的目的是什么？如何对计算区域进行离散化？离散化时通常使用哪些网格？如何对控制方程进行离散？离散化常用的方法有哪些？它们有什么不同？（#80） 4 常见离散格式的性能的对比（稳定性、精度和经济性）(#62) 5 在利用有限体积法建立离散方程时，必须遵守哪几个基本原则？（#81） 6 流场数值计算的目的是什么？主要方法有哪些？其基本思路是什么？各自的适用范围是什么？(#130) 7 可压缩流动和不可压缩流动，在数值解法上各有何特点？为何不可压缩流动在求解时反而比可压缩流动有更多的困难？（#55） 8 什么叫边界条件？有何物理意义？它与初始条件有什么关系？（#56） 9 在一个物理问题的多个边界上，如何协调各边界上的不同边界条件？在边界条件的组合问题上，有什么原则？ 10 在数值计算中，偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别？（#143） 11 在网格生成技术中，什么叫贴体坐标系？什么叫网格独立解？（#35） 12 在GAMBIT的foreground和background中，真实体和虚实体、实操作和虚操作四个之间是什么关系？ 13 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量？及其在做网格时大致注意到哪些细节？（#38） 14 画网格时，网格类型和网格方法如何配合使用？各种方法有什么样的应用范围及做网格时需注意的问题？(#169) 15 对于自己的模型，大多数人有这样的想法：我的模型如何来画网格？用什么样的方法最简单？这样做网格到底对不对？(#154) 16 在两个面的交界线上如果出现网格间距不同的情况时，即两块网格不连续时，怎么样克服这种情况呢？（#40） 17 依据实体在GAMBIT建模之前简化时，必须遵循哪几个原则？(#170) 18 在设置GAMBIT边界层类型时需要注意的几个问题：a、没有定义的边界线如何处理？b、计算域内的内部边界如何处理（2D）？（#128） 19 为何在划分网格后，还要指定边界类型和区域类型？常用的边界类型和区域类型有哪些？（#127） 20 何为流体区域（fluid zone）和固体区域（solid zone）？为什么要使用区域的概念？FLUENT是怎样使用区域的？(#41) 21 如何监视FLUENT的计算结果？如何判断计算是否收敛？在FLUENT中收敛准则是如何定义的？分析计算收敛性的各控制参数，并说明如何选择和设置这些参数？解决不收敛问题通常的几个解决方法是什么？（9楼） 22 什么叫松弛因子？松弛因子对计算结果有什么样的影响？它对计算的收敛情况又有什么样的影响？（7楼） 23 在FLUENT运行过程中，经常会出现“turbulence viscous rate”超过了极限值，此时如何解决？而这里的极限值指的是什么值？修正后它对计算结果有何影响？(#28) 24 在FLUENT运行计算时，为什么有时候总是出现“reversed flow”？其具体意义是什么？有没有