fluent多孔介质参数设定

Chapter 10: Modeling Flow Through Porous Media

This tutorial is divided into the following sections:

10.1. Introduction

10.2. Prerequisites

10.3. Problem Description

10.4. Setup and Solution

10.5. Summary

10.6. Further Improvements

10.1. Introduction

Many industrial applications such as filters, catalyst beds and packing, involve modeling the flow through porous media.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, ANSYS 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.

10.2. Prerequisites

This tutorial is written with the assumption that you have completed one or more of the introductory tutorials found in this manual:

?Introduction to Using ANSYS FLUENT in ANSYS Workbench: Fluid Flow and Heat Transfer in a Mixing Elbow (p.1)

?Parametric Analysis in ANSYS Workbench Using ANSYS FLUENT (p.77)

Chapter 10: Modeling Flow Through Porous Media

?Introduction to Using ANSYS FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow (p.131)

and that you are familiar with the ANSYS FLUENT navigation pane and menu structure. Some steps in the setup and solution procedure will not be shown explicitly.

10.3. Problem Description

The catalytic converter modeled here is shown in Figure 10.1 (p.448).The nitrogen flows 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.

Figure 10.1 Catalytic Converter Geometry for Flow Modeling

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 along the inlet axis.The substrate is imper-meable in other directions.This characteristic is modeled using loss coefficients that are three orders of magnitude higher than in the main flow direction.

10.4. Setup and Solution

The following sections describe the setup and solution steps for this tutorial:

10.4.1. Preparation

10.4.2. Step 1: Mesh

10.4.3. Step 2: General Settings

10.4.4. Step 3: Models

10.4.5. Step 4: Materials

10.4.6. Step 5: Cell Zone Conditions

10.4.7. Step 6: Boundary Conditions

10.4.8. Step 7: Solution

10.4.9. Step 8: Postprocessing

Setup and Solution 10.4.1. Preparation

1.Extract the file porous.zip from the ANSYS_Fluid_Dynamics_Tutorial_Inputs.zip archive

which is available from the Customer Portal.

Note

For detailed instructions on how to obtain the ANSYS_Fluid_Dynamics_Tutori-

al_Inputs.zip file, please refer to Preparation (p.3) in Introduction to Using ANSYS

FLUENT in ANSYS Workbench: Fluid Flow and Heat Transfer in a Mixing Elbow (p.1).

2.Unzip porous.zip to your working folder.

The mesh file catalytic_converter.msh can be found in the porous directory created after

unzipping the file.

https://www.wendangku.net/doc/b93254445.html,e the FLUENT Launcher to start the 3D version of ANSYS FLUENT.

For more information about FLUENT Launcher, see Starting ANSYS FLUENT Using FLUENT

Launcher in the User’s Guide.

4.Enable Double-Precision.

Note

The Display Options are enabled by default.Therefore, once you read in the mesh, it

will be displayed in the embedded graphics window.

10.4.2. Step 1: Mesh

1.Read the mesh file (catalytic_converter.msh).

File?Read?Mesh...

2.Check the mesh.

General?Check

ANSYS FLUENT will perform various checks on the mesh and report the progress in the console. Make

sure that the reported minimum volume is a positive number.

3.Scale the mesh.

General?Scale...

a.Select mm from the Mesh Was Created In drop-down list.

b.Click Scale .

c.

Select mm from the View Length Unit In drop-down list.All dimensions will now be shown in millimeters.

d.Close the Scale Mesh dialog box.

4.

Check the mesh.General ? Check

Note

It is a good idea to check the mesh after you manipulate it (i.e., scale, convert to poly-hedra, merge, separate, fuse, add zones, or smooth and swap.) This will ensure that the quality of the mesh has not been compromised.

5.

Examine the mesh.

Rotate the view and zoom in to get the display shown in Figure 10.2 (p.451).The hex mesh on the geometry contains a total of 34,580 cells.

Chapter 10: Modeling Flow Through Porous Media

Figure 10.2 Mesh for the Catalytic Converter Geometry

10.4.3. Step 2: General Settings

General Setup and Solution

Chapter 10: Modeling Flow Through Porous Media

1.Retain the default solver settings.

10.4.4. Step 3: Models

Models

1.Select the standard - turbulence model.

Models?Viscous?Edit...

Setup and Solution

a.Select k-epsilon (2eqn) in the Model list.

The original Viscous Model dialog box will now expand.

b.Retain the default settings for k-epsilon Model and Near-Wall Treatment and click OK to close

the Viscous Model dialog box.

10.4.5. Step 4: Materials

Materials

1.Add nitrogen to the list of fluid materials by copying it from the FLUENT Database of materials.

Materials?air?Create/Edit...

Chapter 10: Modeling Flow Through Porous Media

a.Click the FLUENT Database... button to open the FLUENT Database Materials dialog box.

Setup and Solution

i.Select nitrogen (n2) in the FLUENT Fluid Materials selection list.

ii.Click Copy to copy the information for nitrogen to your list of fluid materials.

iii.Close the FLUENT Database Materials dialog box.

b.Click Change/Create and close the Create/Edit Materials dialog box.

10.4.6. Step 5: Cell Zone Conditions

Cell Zone Conditions

Chapter 10: Modeling Flow Through Porous Media

1.Set the cell zone conditions for the fluid (fluid ).

Cell Zone Conditions?fluid?Edit...

a.Select nitrogen from the Material Name drop-down list.

b.Click OK to close the Fluid dialog box.

2.Set the cell zone conditions for the substrate (substrate ).

Cell Zone Conditions?substrate?Edit...Setup and Solution

Chapter 10: Modeling Flow Through Porous Media

a.Select nitrogen from the Material Name drop-down list.

b.Enable Porous Zone to activate the porous zone model.

c.Enable Laminar Zone 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 Table 10.1:Values for the

Principle Direction Vectors (p.459).

ANSYS FLUENT automatically calculates the third (z-direction) vector based on your inputs for the first two vectors.The direction vectors determine which axis the viscous and inertial resistance coefficients act upon.

Table 10.1 Values for the Principle Direction Vectors

Direction-2 Vector Direction-1 Vector Axis 01X 10Y 0

Z

Use the scroll bar to access the fields that are not initially visible in the dialog box.ii.

Enter the values in Table 10.2:Values for the Viscous and Inertial Resistance (p.459)Viscous Resistance and Inertial Resistance .

Direction-2 and Direction-3 are set to arbitrary large numbers.These values are several orders of magnitude greater than that of the Direction-1 flow and will make any radial flow insignificant.

Scroll down to access the fields that are not initially visible in the panel.

Table 10.2 Values for the Viscous and Inertial Resistance

Inertial Resistance (1/m)Viscous Resistance (1/m2)Direction 20.4143.846e+07Direction-1204143.846e+10Direction-220414

3.846e+10

Direction-3

e.

Click OK to close the Fluid dialog box.

10.4.7. Step 6: Boundary Conditions

Boundary Conditions

Setup and Solution

1.Set the velocity and turbulence boundary conditions at the inlet (inlet ).

Boundary Conditions?

inlet?Edit...

Chapter 10: Modeling Flow Through Porous Media

a.Enter 22.6 m/s for Velocity Magnitude.

b.Select Intensity and Hydraulic Diameter from the Specification Method drop-down 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 dialog box.

2.Set the boundary conditions at the outlet (outlet).

Boundary Conditions?

outlet?Edit...

Setup and Solution

Chapter 10: Modeling Flow Through Porous Media

a.Retain the default setting of 0 for Gauge Pressure.

b.Select Intensity and Hydraulic Diameter from the Specification Method drop-down 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 dialog box.

3.Retain the default boundary conditions for the walls (substrate-wall and wall).

10.4.8. Step 7: Solution

1.Set the solution parameters.

Solution Methods

Setup and Solution

a.Select Coupled from the Scheme drop-down list.

b.Retain the default selection of Least Squares Cell Based from the Gradient drop-down list in

the Spatial Discretization group box.

c.Retain the default selection of Second Order Upwind from the Momentum drop-down list.

d.Enable Pseudo Transient.

2.Enable the plotting of residuals during the calculation.

Monitors?Residuals?Edit...

a.Retain the default settings.

b.Click OK to close the Residual Monitors dialog box.

3.

Enable the plotting of the mass flow rate at the outlet.Monitors (Surface Monitors) ?

Create...

Chapter 10: Modeling Flow Through Porous Media

Setup and Solution

a.Enable Plot and Write.

b.Select Mass Flow Rate from the Report Type drop-down list.

c.Select outlet in the Surfaces selection list.

d.Click OK to close the Surface Monitor dialog box.

4.Initialize the solution from the inlet.

Solution Initialization

a.Retain the default selection of Hybrid Initialization from the Initialization Methods group box.

b.Click Initialize.

Note

A warning is displayed in the console stating that the convergence tolerance of

1.000000e-06 not reached during Hybrid Initialization.This means that the default

number of iterations is not enough.You will increase the number of iterations and

re-initialize the flow. For more information refer to Hybrid Initialization in the User's

Guide.

c.Click More Settings....

Chapter 10: Modeling Flow Through Porous Media

i.Increase the Number of Iterations to 15.

ii.Click OK to close the Hybrid Initialization dialog box.

d.Click Initialize once mor

e.

Note

Click OK in the Question dialog box, where it asks to discard the current data.The

console displays that hybrid initialization is done.

Note

For flows in complex topologies, hybrid initialization will provide better initial velocity

and pressure fields than standard initialization.This will improve the convergence be-

havior of the solver.

5.Save the case file (catalytic_converter.cas.gz).

File?Write?Case...

6.Run the calculation by requesting 100 iterations.

Run Calculation

多孔介质在fluent中的操作方法 网络上传版本

如何在Fluent中实现多孔介质双能量方程(LNTE) How to use Non-equilibrium Thermal equation (LNTE) model for Porous media in Fluent Software ●请参照本人发表的文章： ●Please refer to the following papers： 1)Wang Fu–Qiang*，Shuai Yong*，Wang Zhi–Q iang，Leng Yu，Tan He–Ping. Thermal and chemical reaction performance analyses of steam methane reforming in porous media solar thermochemical reactor，International Journal of Hydrogen Energy，39(2)：718-730，2014 关键词：Porous, Solar, Hydrogen, Methane, Reforming, P1 approximation, radiative heat transfer 2)Wang Fu–Qiang*，Shuai Yong*，Tan He–Ping，Zhang Xiao-Feng，Mao Qian-Jun，Heat transfer analyses of porous media receiver with multi–dish collector by coupling MCRT and FVM method，Solar Energy，93：158–168，2013 关键词：Solar, Porous, dish concentrator, Receiver, Monte Carlo 3)Wang Fu–Qiang*，Shuai Yong*，Tan He–Ping，Yu Chun–Liang，Thermal Performance Analysis of Porous Media Receiver with Concentrated Solar Irradiation，International Journal of Heat and Mass Transfer，62：247–254，2013 关键词：Solar, Porous, dish concentrator, Receiver, Monte Carlo

学习fluent(流体常识及软件计算参数设置)

luent 中一些问题 ( 目录 ) 离散化的目的 计算区域的离散及通常使用的网格 控制方程的离散及其方法 各种离散化方法的区别 8 9 10在GAMBIT 中显示的“check 主要通过哪几种来判断其网格的质量？及其在做网格时大 致注意到哪些细节？ 11 在两个面的交界线上如果出现网格间距不同的情况时，即两块网格不连续时，怎么样克 服这种情况呢？ 12在设置GAMBIT 边界层类型时需要注意的几个问题： a 、没有定义的边界线如何处理？ b 、计算域内的内部边界如何处理( 2D )？ 13 为何在划分网格后，还要指定边界类型和区域类型？常用的边界类型和区域类型有哪 些？ 14 20 何为流体区域( fluid zone )和固体区域( solid zone )？为什么要使用区域的概念？ FLUENT 是怎样使用区域的？ 15 21 如何监视 FLUENT 的计算结果？如何判断计算是否收敛？在 FLUENT 中收敛准则是 如何定义的？分析计算收敛性的各控制参数，并说明如何选择和设置这些 参数？解决不收1 如何入门 2 CFD 2.1 2.2 2. 3 2.4 2.5 2.6 计算中涉及到的流体及流动的基本概念和术语 理想流体( Ideal Fluid )和粘性流体( Viscous Fluid ) 牛顿流体( Newtonian Fluid )和非牛顿流体( non-Newtonian Fluid ) 可压缩流体 ( Compressible Fluid )和不可压缩流体( Incompressible Fluid ) 层流( Laminar Flow )和湍流( Turbulent Flow ) 定常流动( Steady Flow )和非定常流动( Unsteady Flow ) 亚音 速流动 (Subsonic) 与超音速流动( Supersonic ) 热传导( Heat Transfer )及扩散 ( Diffusion ) 2.7 3 在数值模拟过程中，离散化的目的是什么？如何对计算区域进行离散化？离散化时通常 使用哪些网格？如何对控制方程进行离散？离散化常用的方法有哪些？它们有 什么不 同？ 3.1 3.2 3.3 3.4 4 常见离散格式的性能的对比(稳定性、精度和经济性) 5 流场数值计算的目的是什么？主要方法有哪些？其基本思路是什么？各自的适用范围是 什 么？ 6 可压缩流动和不可压缩流动，在数值解法上各有何特点？为何不可压缩流动在求解时反 而比 可压缩流动有更多的困难？ 6.1 可压缩 Euler 及 Navier-Stokes 方程数值解 6.2 不可压缩 Navier-Stokes 方程求解 什么叫边界条件？有何物理意义？它与初始条件有什么关系？ 在数值计算中，偏微分方程的 双曲型方程、椭圆型方程、抛物型方程有什么区别？ 在网格生成技术中，什么叫贴体坐标 系？什么叫网格独立解？

学习fluent (流体常识及软件计算参数设置)

luent中一些问题----(目录) 1 如何入门 2 CFD计算中涉及到的流体及流动的基本概念和术语 2.1 理想流体（Ideal Fluid）和粘性流体（Viscous Fluid） 2.2 牛顿流体（Newtonian Fluid）和非牛顿流体（non-Newtonian Fluid） 2.3 可压缩流体（Compressible Fluid）和不可压缩流体（Incompressible Fluid） 2.4 层流（Laminar Flow）和湍流（Turbulent Flow） 2.5 定常流动（Steady Flow）和非定常流动（Unsteady Flow） 2.6 亚音速流动(Subsonic)与超音速流动（Supersonic） 2.7 热传导（Heat Transfer）及扩散（Diffusion） 3 在数值模拟过程中，离散化的目的是什么？如何对计算区域进行离散化？离散化时通常使用哪些网格？如何对控制方程进行离散？离散化常用的方法有哪些？它们有什么不 同？ 3.1 离散化的目的 3.2 计算区域的离散及通常使用的网格 3.3 控制方程的离散及其方法 3.4 各种离散化方法的区别 4 常见离散格式的性能的对比（稳定性、精度和经济性） 5 流场数值计算的目的是什么？主要方法有哪些？其基本思路是什么？各自的适用范围是什么？ 6 可压缩流动和不可压缩流动，在数值解法上各有何特点？为何不可压缩流动在求解时反而比可压缩流动有更多的困难？ 6.1 可压缩Euler及Navier-Stokes方程数值解 6.2 不可压缩Navier-Stokes方程求解 7 什么叫边界条件？有何物理意义？它与初始条件有什么关系？ 8 在数值计算中，偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别？ 9 在网格生成技术中，什么叫贴体坐标系？什么叫网格独立解？ 10 在GAMBIT中显示的“check”主要通过哪几种来判断其网格的质量？及其在做网格时大致注意到哪些细节？ 11 在两个面的交界线上如果出现网格间距不同的情况时，即两块网格不连续时，怎么样克服这种情况呢？ 12 在设置GAMBIT边界层类型时需要注意的几个问题：a、没有定义的边界线如何处理？ b、计算域内的内部边界如何处理（2D）？ 13 为何在划分网格后，还要指定边界类型和区域类型？常用的边界类型和区域类型有哪些？ 14 20 何为流体区域（fluid zone）和固体区域（solid zone）？为什么要使用区域的概念？FLUENT是怎样使用区域的？ 15 21 如何监视FLUENT的计算结果？如何判断计算是否收敛？在FLUENT中收敛准则是如何定义的？分析计算收敛性的各控制参数，并说明如何选择和设置这些参数？解决不收

FLUENT中的求解器、算法和离散方法

FLUENT中的求解器、算法和离散方法 作为一个非科班出身的CFD工程师，一开始常常被CFD软件里各种概念搞的晕头转向。最近终于静下心来看了看CFD理论的书，理清了一些概念。就此写一遍博文，顺便整理一下所学内容。 I 求解器： FLUENT中求解器的选择在如下图所示界面中设置： FLUENT中的求解器主要是按照是否联立求解各控制方程来区分的，详见下图：

II 算法： 算法是求解时的策略，即按照什么样的方式和步骤进行求解。FLUENT中算法的选择在如下图所示的界面中设置：

这里简单介绍一下SIMPLE、SIMPLEC、PISO等算法的基本思想和适用范围。 SIMPLE算法：基本思想如前面讲求解器的那张图中解释分离式求解器的例子所示的一样，这里再贴一遍： 1.假设初始压力场分布。 2.利用压力场求解动量方程，得到速度场。 3.利用速度场求解连续性方程，使压力场得到修正。 4.根据需要，求解湍流方程及其他方程 5.判断但前计算是否收敛。若不收敛，返回第二步。 简单说来，SIMPLE算法就是分两步走：第一步预测，第二步修正，即预测－修正。 SIMPLC算法：是对SIMPLE算法的一种改进，其计算步骤与SIMPLE算法相同，只是压力修正项中的一些系数不同，可以加快迭代过程的收敛。 PISO算法：比SIMPLE算法增加了一个修正步，即分三步：第一步预测，第二步修正得到一个修正的场分布，第三步在第二步基础上在进行一侧修正。即预测－修正－修正。PISO算法在求解瞬态问题时有明显优势。对于稳态问题可能SIMPLE 或SIMPLEC更合适。 如果你实在不知道该如何选择，就保持FLUENT的默认选项好了。因为默认选项可以很好解决70％以上的问题，而且对于大部分出了问题的计算来说，也很少是因为算法选择不恰当所致。 III 离散方法： 离散方法是指按照什么样的方式将控制方程在网格节点离散，即将偏微分格式的控制方程转化为各节点上的代数方程组。FLUENT中离散方法的选择在如下图所示的界面中设置：

FLUENT多孔介质中平面面板(plane surface)工具的使用

1、输出grid图形 2、选择surface---plane,打开plane surface面板 3、通过确定三个点来确定平面位置。单击slect point，出现提示，不点选cancel.在grid 图形的多孔介质区域任意位置右键点选3个点。 4、回到plane surface面板，勾选plane tool，则在grid图形的多孔介质区域出现一个平面。

若出现的平面与我们的预期相差比较大的话，可以单击reset points,可以获得一个特殊位置的平面。 5、打开多孔介质的控制面板，选择porou zone标签，点击update from plane tool按钮，获得方向矢量1，和方向矢量2的原始值，并与左下角的坐标系统比较，确定我们大概的旋转方向。 6、对比grid图形左下角的坐标系统，红线和红色箭头代表的是方向矢量1，绿线和绿色箭头代表的是方向矢量2 应该使红线和X正方向平行，绿线和Y正方向平行。具体的操作应该是： 一：先单击白线的蓝色箭头，固定了该方向在旋转过程中不变，可以保证在旋转的过程比较有规律，然后右键点选白线的红色箭头旋转红线的红色箭头到X的正轴； 二: 接下来应该是单击白线的红色箭头，固定该方向不变，单击白线的蓝色箭头，旋转绿线的绿色箭头指向Y的正轴。（所以多孔介质区域我们一般是设置在坐标系统里面，轴线等 与坐标系统无非直角角度关系）。把平面移动到图形外有利于旋转，比较清楚。平面

法线方向的移动是用鼠标右键单击平面阴影部分并拖动，横向移动则需按下shift并进行如上操作。 7、旋转到适当的位置后（鼠标右键拖动箭头），再次点击update from plane tool按钮，获得方向矢量1，和方向矢量2。 得到的数值很可能不是整数，这个时候我们可以把他简化为整数。例如：0.9123可以简化为1，0.01245可以简化为0，以此类推。

辐射和对流模型Fluent参数设置

辐射和对流模型Fluent参数设置 1.读入***.mesh文件，并对网格文件进行进行检查，Grid→cheek，主要看最小体积和最小面积不能为负，之后进行刻度转换，Grid→scale，在Gmbit 里面建模默认尺寸为米，与实际尺寸之间要进行转化，如下图： 2.选择求解器，Define→Models→sover……根据情况选择，如上图：接着选择辐射模型，Define→Models→Radiation，如下图，当Radiation Model面板上 点击ok时，会出现一个信息提示框，告诉你新 的材料物性被添加了，你将在后面设置物性参 数，因此现在只需单击ok确认这个信息即可， 如下图： 注意：当你激活辐射模型后，Fluent会自动打开能量求解器，如下图： 不用再Define→Models→Energy……

3.设置流体粘性，由于模型中空气流速比较大，设成双方程模型：如下图： 4.设置操作条件，此模型此有流体，属有重力情况，Define→Operating Conditions，选中 Gravity.Y方向加速度设置为-9.8 2 m,击OK确定。 /s 设置工作温度，在后面要激活的Boussinesq model要用到，（Boussinesq model：

考虑温度变化而忽略压强变化引起的密度变化叫做Boussinesq 假设） 5. 定义材料并设置其物理属性 Define →Material …… 先定义空气物性，要定义成有浮力的，取Boussinesq 选项。 Density=1.1653/m kg ，()k kg j C p ?=/1005 Thermal Conductivity=0.0267()k m w ?/，Material Type ：fluid ； Thermal Expansion Coefficient =0.0033()k /1。 通过滚动条使先前面板中不可见的物性显示出来。在Scattering Coefficient 和Scattering Phase Function 中保持默认值，在要解决的问题中不涉及到散射问题；设定热扩散系数（用boussinesq 模型时）为1e-5K -1。单击Change/Create ，关闭Materials 面板。 6.设置边界条件Define → Boundary Conditions ……

多孔介质条件多孔介质模型可以应用于很多问题,如通过充满介质的流动

多孔介质条件 多孔介质模型可以应用于很多问题，如通过充满介质的流动、通过过滤纸、穿孔圆盘、流量分配器以及管道堆的流动。当你使用这一模型时，你就定义了一个具有多孔介质的单元区域，而且流动的压力损失由多孔介质的动量方程中所输入的内容来决定。通过介质的热传导问题也可以得到描述，它服从介质和流体流动之间的热平衡假设，具体内容可以参考多孔介质中能量方程的处理一节。 多孔介质的一维化简模型，被称为多孔跳跃，可用于模拟具有已知速度/压降特征的薄膜。多孔跳跃模型应用于表面区域而不是单元区域，并且在尽可能的情况下被使用（而不是完全的多孔介质模型），这是因为它具有更好的鲁棒性，并具有更好的收敛性。详细内容请参阅多孔跳跃边界条件。 1、多孔介质模型的限制 如下面各节所述，多孔介质模型结合模型区域所具有的阻力的经验公式被定义为“多孔”。事实上多孔介质不过是在动量方程中具有了附加的动量损失而已。因此，下面模型的限制就可以很容易的理解了。 ● 流体通过介质时不会加速，因为事实上出现的体积的阻塞并没有在模型中出现。这对于过渡流是有很大的影响的，因为它意味着FLUENT 不会正确的描述通过介质的过渡时间。 ● 多孔介质对于湍流的影响只是近似的。详细内容可以参阅湍流多孔介质的处理一节。 2、多孔介质的动量方程 多孔介质的动量方程具有附加的动量源项。源项由两部分组成，一部分是粘性损失项 (Darcy)，另一个是内部损失项： ∑∑==+=31312 1j j j j ij j ij i v v C v D S ρμ 其中S_i 是i 向(x, y, or z)动量源项，D 和C 是规定的矩阵。在多孔介质单元中，动量损失对于压力梯度有贡献，压降和流体速度（或速度方阵）成比例。 对于简单的均匀多孔介质： j j i i v v C v S ραμ2 12+= 其中a 是渗透性，C2是内部阻力因子，简单的指定D 和C 分别为对角阵1/a 和C2，其它项为零。 FLUENT 还允许模拟的源项为速度的幂率： ()i C C j i v v C v C S 10011-== 其中C_0和C_1为自定义经验系数。 注意：在幂律模型中，压降是各向同性的，C_0的单位为国际标准单位。

Fluent求解参数设置

求解参数设置（Solution Methods/Solution Controls）： 在设置完计算模型和边界条件后，即可开始求解计算了，因为常会出现求解不收敛或者收敛速度很慢的情况，所以就要根据具体的模型制定具体的求解策略，主要通过修改求解参数来完成。在求解参数中主要设置求解的控制方程、选择压力速度耦合方法、松弛因子、离散格式等。 在VOF模型中，PISO比较适合于不复杂的流体，SIMPLE和SIMPLEC适合于可压缩的流体或者处于封闭域中的流体。 ? 求解的控制方程: 在求解参数设置中，可以选择所需要求 解的控制方程。可选择的方程包括Flow(流动方 程)、Turbulence(湍流方程)、Energy(能量方 程)、Volume Fraction(体积分数方程)等。在 求解过程中，有时为了得到收敛的解，先关闭 一些方程，等一些简单的方程收敛后，再开启 复杂的方程一起计算。 ? 选择压力速度耦合方法: 在基于压力求解器中，FLUENT提供了压 力速度耦合的4种方法，即SIMPLE、 SIMPLEC(SIMPLE．Consistent)、PISO以及 Coupled。定常状态计算一般使用SIMPLE或者 SIMPLEC方法，对于过渡计算推荐使用PISO方 法。PISO方法还可以用于高度倾斜网格的定常 状态计算和过渡计算。需要注意的是压力速度 耦合只用于分离求解器，在耦合求解器中不可 以使用。 在FLUENT中，可以使用标准SIMPLE算法和SIMPLEC算法，默认是SIMPLE算法，但对于许多问题如果使用SIMPLEC可能会得到更好的结果，尤其是可以应用增加的亚松弛迭代时。 对于相对简单的问题(如没有附加模型激活的层流流动)，其收敛性可以被压力速度耦合所限制，用户通常可以使用SIMPLEC算法很快得到收敛解。在SIMPLEC算法中，压力校正亚松弛因子通常设为1.0，它有助于收敛，但是，在有些问题中，将压力校正松弛因子增加到1.0可能会导致流动不稳定，对于这种情况，则需要使用更为保守的亚松弛或者使用SIMPLE算法。对于包含湍流或附加物理模型的复杂流动，只要用压力速度耦合做限制，SIMPLEC就会提高收敛性，它通常是一种限制收敛性的附加模拟参数，在这种情况下，SIMPLE和SIMPLEC 会给出相似的收敛速度。 对于所有的过渡流动计算，推荐使用PISO算法邻近校正。它允许用户使用大的时间步，而且对于动量和压力都可以使用亚松弛因子1.0。对于定常状态问题，具有邻近校正的PISO并不会比具有较好的亚松弛因子的SIMPLE或SIMPLEC好。对于具有较大扭曲网格上的定常状态和过渡计算推荐使用PISO倾斜校正。 当使用PISO邻近校正时，对所有方程都推荐使用亚松弛因子为1.0或者接近1.0。如果只对高度扭曲的网格使用PISO倾斜校正，则要设定动量和压力的亚松弛因子之和为1.0(例如，压力亚松弛因子0.3，动量亚松弛因子0.7)。

多孔介质边界条件

7.19多孔介质边界条件 多孔介质模型适用的范围非常广泛，包括填充床，过滤纸，多孔板，流量分配器，还有管群，管束系统。当使用这个模型的时候，多孔介质将运用于网格区域，流场中的压降将由输入的条件有关，见Section 7.19.2.同样也可以计算热传导，基于介质和流场热量守恒的假设，见Section 7.19.3. 通过一个薄膜后的已知速度/压力降低特性可以简化为一维多孔介质模型，简称为“多孔跳跃”。多孔跳跃模型被运用于一个面区域而不是网格区域，而且也可以代替完全多孔介质模型在任何可能的时候，因为它更加稳定而且能够很好地收敛。见Section 7.22. 7.19.1 多孔介质模型的限制和假设 多孔介质模型就是在定义为多孔介质的区域结合了一个根据经验假设为主的流动阻力。本质上，多孔介质模型仅仅是在动量方程上叠加了一个动量源项。这种情况下，以下模型方面的假设和限制就可以很容易得到： ?因为没有表示多孔介质区域的实际存在的体，所以fluent默认是计算基于连续性方程的虚假速度。做为一个做精确的选项，你可以适用fluent 中的真是速度，见section7.19.7。 ?多孔介质对湍流流场的影响，是近似的，见7.19.4。 ?当在移动坐标系中使用多孔介质模型的时候，fluent既有相对坐标系也可以使用绝对坐标系，当激活相对速度阻力方程。这将得到更精确的源项。 相关信息见section7.19.5和7.19.6。 ?当需要定义比热容的时候，必须是常数。 7.19.2 多孔介质模型动量方程 多孔介质模型的动量方程是在标准动量方程的后面加上动量方程源项。源项包含两个部分：粘性损失项（达西公式项，方程7.19-1右边第一项），和惯性损失项（方程7.19-1右边第二项） (7.19-1)

详细FLUENT实例讲座翼型计算

详细FLUENT实例讲座翼型计算 部门： xxx 时间： xxx 整理范文，仅供参考，可下载自行编辑

CAE联盟论坛精品讲座系列 详细FLUENT实例讲座-翼型计算 主讲人：流沙 CAE联盟论坛总版主 1.1 问题描述 翼型升阻力计算是CFD最常规的应用之一。本例计算的翼型为 RAE2822，其几何参数可以查看翼型数据库。本例计算在来流速度0.75马赫，攻角3.19°情况下，翼型的升阻系数及流场分布，并将计算结果与实验数据进行对比。模型示意图如图1所示。 b5E2RGbCAP 1.p ng(12.13 K>2018/7/29 23:41:251.2 FLUENT前处理设置Step 1：导入计算模型 以3D，双精度方式启动FLUENT14.5。 利用菜单【File】>【Read】>【Mesh…】，在弹出的文件选择对话框中选择网格文件rae2822_coarse.msh，点击OK按钮选择文件。如图2所示。p1EanqFDPw

点击FLUENT模型树按钮General，在右侧设置面板中点击按钮Display…，在弹出的设置对话框中保持默认设置，点击Display按钮，显示网格。如图3所示。DXDiTa9E3d 2.png(11.51 K>2018/7/29 23:41:25

3.png(33.41 K>2018/7/29 23:41:253-2.png(52.04 K>2018/7/29 23:41:25Step 2：检查网格 采用如图4所示步骤进行网格的检查与显示。点击FLUENT模型树节点General节点，在右侧面板中通过按钮Scale…、Check及 Report Quality实现网格检查。 4.png(12. 10 K>RTCrpUDGiT2018/7/29 23:41:25点击按钮Check，在命令输出按钮出现如图5所示网格统计信息。从图中可以看出，网格尺寸分布： x轴：-48.97~50m

fluent求解器

Model slover slover:求解器 Pressuere based: 基于压力 Density based：基于密度 Formulation:算法 implicit: 隐式算法 explicit：显式算法 space：选择空间属性 2D：二维空间 Axisymmetric：轴对称空间 Axisymmetric swirl：轴对称旋转空间 3D：三维空间 time：时间 steady:稳态 unsteady:非稳态 (~~~~~~~~~~~~~~~~~) velocity formulation：速度属性 absolute：绝对速度 relative：相对速度 Gradient option:梯度选项 Green-Gauss Cell-Based：格林－高斯基于单元体－默认方法；解有伪扩散（求解域的拖尾现象）。伪扩散是指在平流扩散方程数值解中因平流项有限差分的截断误差引起的虚假扩散。这是解方程欧拉型模式所特有的。其大小与所用的有限差分格式有关，有时甚至完全掩盖方程中其他扩散项的作用。为克服伪扩散，须采取特殊的技术措施和各种不同的差分格式。 Green-Gauss Node-Based：格林－高斯基于节点－更精确；最小化伪扩散；推荐用在三角网格上。 least-quares cell based：基于单元体的最小二乘法插值。推荐用于多面体网格，与基于节点的格林-高斯格式具有相同的 精度和格式。 porous formulation：选择多孔算法 superficial velocity：表面速度算法 physical velocity：物理速度算法 Multiphase Model:多相流模型 Volume of Fluid：VOF模型 Mixture：混合模型 Eulerian：欧拉模型 Energy：能量方程

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可以考虑为零。忽略对流加速以及扩散，多孔介质模型简化为Darcy定律： 在多孔介质区域三个坐标方向的压降为：

FLUENT求解器设置

FLUENT求解器设置主要包括：1、压力-速度耦合方程格式选择2、对流插值3、梯度插值4、压力插值 下面对这几种设置做详细说明。 一、压力-速度耦合方程求解算法 FLUENT中主要有四种算法：SIMPLE，SIMPLEC，PISO，FSM (1)SIMPLE(semi-implicit method for pressure-linked equations)半隐式连接压力方程方法，是FLUENT的默认格式。 (2)SIMPLEC(SIMPLE-consistent)。对于简单的问题收敛非常快速，不对压力进行修正，所以压力松弛因子可以设置为1 (3)Pressure-Implicit with Splitting of Operators (PISO)。对非定常流动问题或者包含比平均网格倾斜度更高的网格适用 (4)Fractional Step Method (FSM)对非定常流的分步方法。用于NITA格式，与PISO具有相同的特性。 二、对流插值（动量方程） FLUENT有五种方法：一阶迎风格式、幂率格式、二阶迎风格式、MUSL三阶格式、QUICK格式 （1）FLUENT默认采用一阶格式。容易收敛，但精度较差，主要用于初值计算。 （2）Power Lar.幂率格式，当雷诺数低于5时，计算精度比一阶格式要高。 （3）二阶迎风格式。二阶迎风格式相对于一阶格式来说，使用更小的截断误差，适用于三角形、四面体网格或流动与网格不在同一直线上；二阶格式收敛可能比较慢。 （4）MUSL(monotone upstream-centered schemes for conservation laws).当地3阶离散格式。主要用于非结构网格，在预测二次流，漩涡，力等时更精确。 （5）QUICK（Quadratic upwind interpolation）格式。此格式用于四边形/六面体时具有三阶精度，用于杂交网格或三角形/四面体时只具有二阶精度。 三、梯度插值梯度插值主要是针对扩散项。 FLUENT有三种梯度插值方案：green-gauss cell-based，Green-gauss node-based，least-quares cell based. （1）格林-高斯基于单元体。求解方法可能会出现伪扩散。 （2）格林-高斯基于节点。求解更精确，最小化伪扩散，推荐用于三角形网格上 （3）基于单元体的最小二乘法插值。推荐用于多面体网格，与基于节点的格林-高斯格式具有相同的精度和格式。 四、压力插值压力基分离求解器主要有五种压力插值算法。 （1）标准格式（Standard)。为FLUENT缺省格式，对大表妹边界层附近的曲线发现压力梯度流动求解精度会降低（但不能用于流动中压力急剧变化的地方——此时应该使用PRESTO!格式代替） （2）PRESTO!主要用于高旋流，压力急剧变化流（如多孔介质、风扇模型等），或剧烈弯曲的区域。 （3）Linear（线性格式）。当其他选项导致收敛困难或出现非物理解时使用此格式。 （4）second order（二阶格式）。用于可压缩流动，不能用于多孔介质、阶跃、风扇、VOF/MIXTURE多相流。 （5）Body Force Weighted体积力。当体积力很大时，如高雷诺数自然对流或高回旋流动中采用此格式。

完整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相同。总结一下，与君共享!

Fluent计算多孔介质模型资料

广东省深圳市宝安区沙井辛养社区西部工业园 TEL:+86-755-3366-8888 FAX:+86-755-3366-0612 Fluent计算多孔介质模型资料 这是一个多孔介质例子，进口速度为0.01m/s,组份为液态水和氧气，其中氧气从多孔介质porous jump 渗透过去，如何看氧气在tissue中扩散的。 porous jump的face permeability1 a=e-8 m_2 thickness 设为0.0001 pressure jump coefficient为默认 porous zone设置如下： direction vector 1, 1, viscous resistance 100 each inertial resistance 100 each porosity 0.1 边界条件设置如下： Ab – wall - default Bc – wall – default Be – porous jump – face permeability 1e-8, porous medium thickness 0.0001 Cd – outflow rating – 0.5 De – wall – default Default interior – interior Default interior001 – interior Default interior019 – interior Ef – wall - default Fg – outflow rating – 1 Fluid - porous zone - direction vector 1, 1, viscous resistance 100 each, inertial resistance 100 each, porosity 0.1 Gh- wall - default Hi – wall - default Hk - porous jump same conditions as other Ij – outflow – 0.5 Jk – wall – default Kl – wall – default

多孔介质的网络模型构建-宫法明

2012年齐鲁大学生软件设计大赛命题 多孔介质的网络模型构建 （中国石油大学宫法明） 一、课题背景简介 多孔介质是指内部含有大量空隙（void）的固体，固体骨架遍及多孔介质所占据的体积空间。多孔介质内部的空隙极其微小。储集石油和天然气的砂岩地层的空隙直径大多在不足1微米到500微米之间；毛细血管内径一般为5～15微米；肺泡-微细支气管系统的空隙直径一般为200微米左右或更小；植物体内输送水分和糖分的空隙直径一般不大于40微米。 一般多孔介质的空隙都是相通的，也可能是部分连通、部分不连通的。由于多孔介质本身的不均匀性、随机性和几何拓扑结构的复杂性，其内部渗透特性、流体传递过程等难以实测。因此，利用计算机对多孔介质进行微观建模，通过计算获取多孔介质的相关构造参数具有重要的研究价值。 注：本竞赛题目来自目前在研的一项国家科技重大专项课题，是其中的一部分，属于比较关键的基础研究，选报本题目的参赛选手在锻炼自己的同时，取得的任何一点成果，都很有可能会为国家做出重要的贡献。 二、课题研究的基本思路及环节 用计算机对多孔介质进行相关研究的基本思路及环节是： ①借助工业用微焦点CT 系统（目前已在使用纳米测量技术，数据更丰富，精度更准确），获取一系列能够真实描述多孔介质的微观空隙结构的CT 切片图像；图1所示为其中一张：

图1：CT切片图像 ②对每幅CT图像进行分割，找到空隙和固体骨架之间的边界，从而可以将固体部分剔除，只留下空隙部分所占据的平面区域；图2所示为分割结果（一个矩形的部分区域）中的一张（黑色 部分为空隙）： 图2：分割结果

③将一系列CT图像中空隙部分所占据的区域叠加在一起，便构成了整个体积空间中所有空隙构成的一个三维体，从而可以用三维显示技术将空隙空间显示出来；如图3所示： 图3：空隙的三维体数据 ④上述步骤产生的空隙体数据一般数据量较大，影响显示的实时性，且大量空隙相互遮挡，不利观察，也不利于后续的各种参数计算，因此需要构建空隙空间的几何模型，通俗的说，就是在空隙体数据外围包上一层皮（一般是网状的，如四边形网格或三角网格），对这个“皮”进行材质、光照等设置之后显示出来，效果就有了较真实的展示。如图4所示：

FLUENT操作过程及全参数选择

振动流化床仿真操作过程及参数选择 1创建流化床模型。 根据靳海波论文提供的试验机参数，创建流化床模型。流化床直148mm 高1m开孔率9%孔径2mm在筛板上铺两层帆布保证气流均布。 因为实验机为一个圆形的流化床，所以可简化为仅二维模型。而实际实验中流化高度远小于1m甚至500mm所以为提高计算时间，可将模型高度缩为500mm由于筛板上铺设两层帆布以达到气流均分的目的，所以认为沿整个筛板的进口风速为均匀的。最终简化模型如下图所示： 上图为流化后的流化床模型，可以看出流化床下端的网格相对上端较密，因为流化行为主要发生的流化床下端，为了加快计算时间，所以采用这种下密上疏的划分方式。其中进口设置为velocity inlet ;出口设置为outflow ；左右两边分为设置为wall。在GAMBIT中设置完毕后，输出二维模型vfb.msh。 outflow 边界条件不需要给定任何入口的物理条件，但是应用也会有限制，大致为以下四点： 1.只能用于不可压缩流动

2.出口处流动充分发展 3.不能与任何压力边界条件搭配使用（压力入口、压力出口） 4.不能用于计算流量分配问题（比如有多个出口的问题） 2 打开FLUENT 6326，导入模型vfb.msh 点击GRID—CHECK检查网格信息及模型中设置的信息，核对是否正确，尤其查看是否出现负体积和负面积，如出现马上修改。核对完毕后，点击GRID-SCAL弹出SCALEGRID窗口，设置单位为mm 并点击change length unit 按钮。具体设置如下： 3设置求解器 保持其他设置为默认，更改TIME为unsteady，因为实际流化的过程是随时间变化的。 （1）pressure based 求解方法在求解不可压流体时，如果我们联立求解 从动量方程和连续性方程离散得到的代数方程组，可以直接得到各速

学习fluent(流体常识及软件计算参数设置)

luent 中一些问题( 目录) 1 如何入门 2 CFD 计算中涉及到的流体及流动的基本概念和术语 2.1 理想流体( Ideal Fluid )和粘性流体( Viscous Fluid ) 2.2 牛顿流体( Newtonian Fluid )和非牛顿流体( non-Newtonian Fluid ) 2.3 可压缩流体( Compressible Fluid )和不可压缩流体( Incompressible Fluid ) 2.4 层流( Laminar Flow )和湍流( Turbulent Flow ) 2.5 定常流动( Steady Flow )和非定常流动( Unsteady Flow ) 2.6 亚音速流动(Subsonic) 与超音速流动( Supersonic ) 2.7 热传导( Heat Transfer )及扩散( Diffusion ) 3 在数值模拟过程中，离散化的目的是什么？如何对计算区域进行离散化？离散化时通常使用哪些网格？如何对控制方程进行离散？离散化常用的方法有哪些？它们有什么不同？ 3.1 离散化的目的 3.2 计算区域的离散及通常使用的网格 3.3 控制方程的离散及其方法 3.4 各种离散化方法的区别 4 常见离散格式的性能的对比(稳定性、精度和经济性) 5 流场数值计算的目的是什么？主要方法有哪些？其基本思路是什么？各自的适用范围是什么？ 6 可压缩流动和不可压缩流动，在数值解法上各有何特点？为何不可压缩流动在求解时反而比可压缩流动有更多的困难？ 6.1 可压缩Euler 及Navier-Stokes 方程数值解 6.2 不可压缩Navier-Stokes 方程求解 7 什么叫边界条件？有何物理意义？它与初始条件有什么关系？ 8 在数值计算中，偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别？ 9 在网格生成技术中，什么叫贴体坐标系？什么叫网格独立解？ 10在GAMBIT中显示的“check主要通过哪几种来判断其网格的质量？及其在做网格时大致注意到哪些细节？ 11 在两个面的交界线上如果出现网格间距不同的情况时，即两块网格不连续时，怎么样克服这种情况呢？ 12 在设置GAMBIT 边界层类型时需要注意的几个问题：a 、没有定义的边界线如何处理？ b、计算域内的内部边界如何处理( 2D)？ 13 为何在划分网格后，还要指定边界类型和区域类型？常用的边界类型和区域类型有哪 些？ 14 20 何为流体区域( fluid zone )和固体区域( solid zone )？为什么要使用区域的概念？FLUENT 是怎样使用区域的？ 15 21 如何监视FLUENT 的计算结果？如何判断计算是否收敛？在FLUENT 中收敛准则是 如何定义的？分析计算收敛性的各控制参数，并说明如何选择和设置这些参数？解决不收

fluent求解器资料

压力基求解器 在压力基求解器中，控制方程是依次求解的。压力基求解器是从原来的分离式求解器发展来的，按顺序仪次求解动量方程、压力修正方程、能量方程和组分方程及其他标量方程，如湍流方程等，和之前不同的是，压力基求解器还增加了耦合算法，可以自由在分离求解和耦合求解之间转换，需要注意的是，在压力基求解器中提供的几个物理模型，在密度基求解器中是没有的。这些物理模型包括：流体体积模型(VOF)，多项混合模型，欧拉混合模型，PDF燃烧模型，预混合燃烧模型，部分预混合燃烧模型，烟灰和NOx模型，Rosseland辐射模型，熔化和凝固等相变模型，指定质量流量的周期流动模型，周期性热传导模型和壳传导模型等。 与密度基求解器的区别： 区别1：压力基求解器主要用于低速不可压缩流动的求解，而密度基求解器则主要针对高速可压缩流动而设计，但是现在两种方法都已经拓展成为可以求解很大流动速度范围的求解方法。两种求解方法的共同点是都使用有限容积的离散方法，但线性化和求解离散方程的方法不同。 区别2：密度基求解器从原来的耦合求解器发展来的，同时求解连续性方程、动量方程、能量方程和组分方程。然后依次再求解标量方程。（注：密度基求解器不求解压力修正方程，因为其压力是由状态方程得出的）。密度基求解器收敛速度快，需要内存和计算量比压力基求解器要大！ 特点： 适用于压力基但不适用于密度基的模型： （1）空化模型 （2） VOF模型 （3） Mixture多相流模型 （4） Eulerian多相流模型 （5）非预混燃烧模型 （6）预混燃烧模型 （7）部分预混燃烧模型 （8）组合PDF传输模型 密度基求解器(Coupled Sover）是同时fluent求解连续方程、动量方程、能量方程及组分输运方程的耦合方程组，然后逐一地求解湍流标量方程。由于控制方程是非线性的，且相互之间是耦合的，因此，在得到收敛解之前，要经过多轮迭代: 1)根据当前的解的结果，更新所有流动变量。如果计算刚刚开始，则用初始值来更新。 2)同时求解连续方程、动量方程、能量方程及组分输运方程的耦合方程组（后两