1、实验研究实验研究理论分析理论分析数值模拟数值模拟研究流体运动规律研究流体运动规律的三种基本方法的三种基本方法连续性方程为：三种基本方法的比较三种基本方法的比较方法方法优优 点点缺缺 点点实验实验研究研究1.结果真实可信结果真实可信2.理论研究和数值方法的基理论研究和数值方法的基础础1.受模型尺寸、人身安全、受模型尺寸、人身安全、测量精度等限制测量精度等限制2.耗费人力、物力、财力耗费人力、物力、财力3.周期长周期长理论理论分析分析1.结果具有普遍性结果具有普遍性2.指导实验、验证数值计算指导实验、验证数值计算1.抽象和简化计算对象抽象和简化计算对象2.非线性方程，很难给出非线性方程，很难给出解
2、析结果解析结果数值数值模拟模拟1.克服了实验研究和理论分克服了实验研究和理论分析的缺点析的缺点2.形象地再现流动情景形象地再现流动情景程序的编制及资料的收集、程序的编制及资料的收集、整理，很大程度上依赖于整理，很大程度上依赖于经验与技巧经验与技巧数值模拟数值模拟iijijijjijigxpuuxuxuuxtu)()(连续性方程为：连续性方程：连续性方程：任何流体运动规律的基础：任何流体运动规律的基础：动量方程：动量方程：能量方程：能量方程：v一组非线性偏微分方程一组非线性偏微分方程v获得这些偏微分方程的精确解析解：很难获得这些偏微分方程的精确解析解：很难miiSxutRfjpjjjpjpSST
3、ucxTxTucxtTc)()(v计算流体力学（计算流体力学（CFDCFD）看作是流动基本方程控制下对流看作是流动基本方程控制下对流动的数值模拟。通过计算机数值计算和图形显示，对包动的数值模拟。通过计算机数值计算和图形显示，对包含有流体运动和热传导等相关物理现象的系统所作的分含有流体运动和热传导等相关物理现象的系统所作的分析。析。vCFDCFD的基本思想：的基本思想：把原来在时间、空间坐标中连续的物把原来在时间、空间坐标中连续的物理量的场（如速度场、温度场、浓度场等），用有限个理量的场（如速度场、温度场、浓度场等），用有限个离散点上的值的集合来代替，按一定方式建立起关于这离散点上的值的集合来代
4、替，按一定方式建立起关于这些值的代数方程（称为离散方程）并求解之，以获得物些值的代数方程（称为离散方程）并求解之，以获得物理量场的近似值。理量场的近似值。v 数值解法很好地利用了数值解法很好地利用了离散离散的概念。的概念。v常用的离散化方法常用的离散化方法：有限差分法、有限元法和有限体积有限差分法、有限元法和有限体积法。法。vFluent是基于有限体积法的，是基于有限体积法的， 是是目前目前CFD应用最广的一种应用最广的一种 方法，方法，将计算区域划分为网将计算区域划分为网 格，并使每个网格点周围有格，并使每个网格点周围有 一个互不重复的控制体积，一个互不重复的控制体积， 将待解的微分方程对每
5、个控将待解的微分方程对每个控 制体积积分，从而得到一组制体积积分，从而得到一组 离散方程。离散方程。 离散化方法离散化方法数值解法的基本思路与步骤数值解法的基本思路与步骤 数值解法的思路与步骤数值解法的思路与步骤 通俗理解通俗理解 1. 1.建立反映工程问题或物理问题本质的数学模型，建立建立反映工程问题或物理问题本质的数学模型，建立针对控制方程的数值离散化方法。针对控制方程的数值离散化方法。 2.2.前处理，建立坐标，构建几何体，网格划分，确定边前处理，建立坐标，构建几何体，网格划分，确定边界条件的类型。界条件的类型。 构体和网格划分技巧构体和网格划分技巧 3. 3.程序计算，选取模型，输入初
6、始条件，确定边界条件、程序计算，选取模型，输入初始条件，确定边界条件、设定控制参数。设定控制参数。 经验取值经验取值 ；与；与2 2有反复修正的过程有反复修正的过程 4. 4.后处理，显示计算结果，以图表展示，判断分析结果后处理，显示计算结果，以图表展示，判断分析结果的准确性和质量。的准确性和质量。vPHOENICSPHOENICS：世界上第一套计算流体动力学与传热学的商世界上第一套计算流体动力学与传热学的商用软件。用软件。vCFXCFX：第一个通过第一个通过ISO9001ISO9001质量认证的商业质量认证的商业CFDCFD软件，英国软件，英国AEAAEA TechnologyTechnol
7、ogy公司开发，公司开发，20032003年被年被ANSYSANSYS收购。收购。vFluentFluent：继继PHOENICSPHOENICS之后的第二个投放市场的基于有限之后的第二个投放市场的基于有限体积法的软件，美国体积法的软件，美国FluentFluent公司于公司于19831983年推出，年推出，20062006年年被被ANSYSANSYS收购。目前功能最全面、适用性最广、国内使用收购。目前功能最全面、适用性最广、国内使用最广泛的最广泛的CFDCFD软件之一，用软件之一，用C C语言编写。语言编写。vSTAR-CDSTAR-CD，FIDAPFIDAP等软件。等软件。 常用的常用的C
8、FD软件软件v前处理器：前处理器：主要负责网格的生成。主要负责网格的生成。 FLUENTFLUENT划分网格的途径有两种：划分网格的途径有两种： 用用FLUENTFLUENT提供的专用网格软件提供的专用网格软件GAMBITGAMBIT进行网格划分进行网格划分 由其他的由其他的CADCAD软件完成造型工作，再导入软件完成造型工作，再导入GAMBITGAMBIT中生成网中生成网格，如格，如Pro/EPro/E、SolidWorksSolidWorks、SolidedgeSolidedge等。等。v求解器：求解器：流体计算的核心，所有计算在此完成。流体计算的核心，所有计算在此完成。v后处理器：后处理
9、器：FluentFluent求解器本身就附带有比较强大的后处求解器本身就附带有比较强大的后处理功能。理功能。OriginOrigin、TecplotTecplot等也是比较专业的后处理器，等也是比较专业的后处理器，可以把一些数据可视化。可以把一些数据可视化。 Fluent软件包软件包v安装安装ExceedExceed：ExceedExceed是在是在WindowsWindows环境下模拟的环境下模拟的UNIXUNIX软件，因为软件，因为GambitGambit必须在必须在UNIXUNIX环境下才可以运行。环境下才可以运行。v安装安装GambitGambit，拷贝，拷贝LicenceLicenc
10、e到到Gambit.IncGambit.Inc/license/license目目录下。录下。v安装安装FluentFluent，拷贝，拷贝LicenceLicence到到Fluent.IncFluent.Inc/license/license目目录下。录下。 如需要卸载，以相反的顺序卸载即可。如需要卸载，以相反的顺序卸载即可。 Fluent软件包的安装软件包的安装v .dbsdbs文件：文件：GambitGambit工作文件，若想修改网格，可以打开工作文件，若想修改网格，可以打开这个文件进行再编辑。这个文件进行再编辑。v . mshmsh文件：文件：GambitGambit输出的网格文件。输
11、出的网格文件。v . cascas文件：文件：. .mshmsh文件经过文件经过FluentFluent处理以后得到的文件。处理以后得到的文件。v . datdat文件：文件：FluentFluent计算数据结果的数据文件。计算数据结果的数据文件。v . joujou文件：文件：日志文档，可以编辑运行。日志文档，可以编辑运行。 Fluent软件包的文件形式软件包的文件形式利用利用Fluent软件进行求解的步骤：软件进行求解的步骤：v（1）运行）运行Gambit，确定几何形状，构造几何体，输出，确定几何形状，构造几何体，输出 .dbs文件文件v（2）生成计算网格，输出）生成计算网格，输出.msh
12、文件文件v（3）运行）运行Fluent软件，选择合适的解算器：软件，选择合适的解算器：2D（二维单（二维单精度）、精度）、3D（三维单精度）、（三维单精度）、2DDP（二维双精度）、（二维双精度）、3DDP（三维双精度）。（三维双精度）。v（4）输入网格，即读入）输入网格，即读入.msh文件文件v（5）检查网格）检查网格 网格检查是最容易出的问题是网格体积为负数。如网格检查是最容易出的问题是网格体积为负数。如果出现负体积，计算就不可以进行。果出现负体积，计算就不可以进行。利用利用Fluent软件进行求解的步骤：软件进行求解的步骤：v（6）选择解算器的格式）选择解算器的格式 Fluent提供提供
13、3种不同的解格式：分离解、隐式耦合解、种不同的解格式：分离解、隐式耦合解、显式耦合解。分离解是按顺序解连续性方程、动量方程显式耦合解。分离解是按顺序解连续性方程、动量方程、能量方程以及组分方程，耦合解是同时解。隐式和显、能量方程以及组分方程，耦合解是同时解。隐式和显式解法的区别在于线性耦合方程的方式不同。式解法的区别在于线性耦合方程的方式不同。 （7）选择需要解的基本方程）选择需要解的基本方程 层流还是湍流（或者无粘流）、化学组分还是化学层流还是湍流（或者无粘流）、化学组分还是化学反应、热传导模型等。反应、热传导模型等。 （8）指定材料物理性质）指定材料物理性质 （9）指定边界条件）指定边界条
14、件利用利用Fluent软件进行求解的步骤：软件进行求解的步骤：v（10）调节解的控制参数（松弛因子、限值、一、二阶）调节解的控制参数（松弛因子、限值、一、二阶、迭代误差等）、迭代误差等） 设定松弛因子：设定松弛因子： 计算流体力学中要求解非线性方程，求解过程中控计算流体力学中要求解非线性方程，求解过程中控制变量的变化可以通过松弛因子来实现，来控制变量在制变量的变化可以通过松弛因子来实现，来控制变量在每次迭代中的变化，从而控制收敛的速度并改善收敛的每次迭代中的变化，从而控制收敛的速度并改善收敛的状况。等于状况。等于1，相当于不用松弛因子；大于，相当于不用松弛因子；大于1，为超松弛，为超松弛因子，
15、加快收敛速度；小于因子，加快收敛速度；小于1，为欠从弛因子，改善收敛，为欠从弛因子，改善收敛的条件。的条件。Fluent里面用的是欠松弛，在里面用的是欠松弛，在01之间，越小之间，越小表示两次迭代值之间变化越小，也就越稳定，但收敛也表示两次迭代值之间变化越小，也就越稳定，但收敛也就越慢。就越慢。利用利用Fluent软件进行求解的步骤：软件进行求解的步骤：v（11）保存设置，输出）保存设置，输出.cas文件文件v（12）初始化流场）初始化流场 需要初始化流场提供一个初始解。可以从一个或者多个需要初始化流场提供一个初始解。可以从一个或者多个边界条件算出初始解，也可以分别输入流场的数值。初始化边界条
16、件算出初始解，也可以分别输入流场的数值。初始化流场有时候对于后面计算的影响很大，需要依靠经验选择最流场有时候对于后面计算的影响很大，需要依靠经验选择最佳的边界，设定合理的数值，进行初始化。佳的边界，设定合理的数值，进行初始化。v（13）开始计算）开始计算 输入迭代步数，报告步长等。输入迭代步数，报告步长等。 计算过程中，有选择性地观察残差，调整相关系数，使计算过程中，有选择性地观察残差，调整相关系数，使得残差趋于平稳，并达到规定的残差值，保证收敛。收敛是得残差趋于平稳，并达到规定的残差值，保证收敛。收敛是一个很重要的问题，一般达到收敛，才可以认为计算结束。一个很重要的问题，一般达到收敛，才可以
17、认为计算结束。利用利用Fluent软件进行求解的步骤：软件进行求解的步骤：v（14)检查结果检查结果 速度、温度、浓度分布的等值线、矢量图速度、温度、浓度分布的等值线、矢量图 (15)保存结果，输出保存结果，输出.dat文件文件 (16)根据需要进行后处理，利用根据需要进行后处理，利用Oigin、Tecplot等软件等软件完善温度分布图、速度分布场、烟气成分场等等。必要完善温度分布图、速度分布场、烟气成分场等等。必要的话，可以细化网格，改变数值以及物理模型。的话，可以细化网格，改变数值以及物理模型。网格化分网格化分 v GAMBIT 软件既可以采用结构化网格，也可以采用非结软件既可以采用结构化
18、网格，也可以采用非结构化网格。对于不同问题，可以划分二维的三角形和四边构化网格。对于不同问题，可以划分二维的三角形和四边形网格，三维的四面体网格、六面体网格、金字塔型网格形网格，三维的四面体网格、六面体网格、金字塔型网格、楔型网格，以及由上述网格类型构成的混合型网格。、楔型网格，以及由上述网格类型构成的混合型网格。v 结构化网格生成速度快、质量好、结构简单，但只适用于结构化网格生成速度快、质量好、结构简单，但只适用于形状规则的图形，因此使用的范围比较窄。对于结构复杂形状规则的图形，因此使用的范围比较窄。对于结构复杂的几何结构，可以选用非结构化网格，降低了划分的难度的几何结构，可以选用非结构化网
19、格，降低了划分的难度，但是会减慢收敛速度、降低准确性。，但是会减慢收敛速度、降低准确性。v 在在Gambit中，只有中，只有map和和submap生成的是结构化网格，生成的是结构化网格，其余均为非结构化网格。其余均为非结构化网格。网格化分网格化分v GAMBIT提供了以下面网格划分提供了以下面网格划分Type选项。选项。 选项选项 详细说明详细说明Map 生成一个普通的结构化网格单元网格生成一个普通的结构化网格单元网格 Submap 将一个不可图示的面分成可图示区域并在每将一个不可图示的面分成可图示区域并在每个区域生成结构化网格单元网格个区域生成结构化网格单元网格 Pave 生成非结构化网格单
20、元网格生成非结构化网格单元网格 Tri Primitive 将一个三侧面分成三个四边形区域并在每个将一个三侧面分成三个四边形区域并在每个区域生成可图示的网格区域生成可图示的网格 Wedge Primitive 在楔形面的尖部生成三角形网格单元并从尖在楔形面的尖部生成三角形网格单元并从尖部向外生成放射状网格部向外生成放射状网格 网格化分网格化分v 如上所述，每个如上所述，每个Elements选项与上面列举的一个或者多个选项与上面列举的一个或者多个Type选项相关联。下表中显示了每个面网格划分选项相关联。下表中显示了每个面网格划分Elements和和Type选项之间的关系。（选项之间的关系。（注意
21、注意：由：由“X”标记得阴影单元代表允标记得阴影单元代表允许选项的组合。）许选项的组合。）ElementsType QuadTriQuad/TriMap X XSubmap X Pave XXXTri Primitive X Wedge Primitive XFluent模拟效果图 煤粉炉煤粉炉 低低NOx燃烧器燃烧器 算例算例v 1.冷、热水混合器内部二维流动冷、热水混合器内部二维流动 一个冷、热水混合器的内部流动与热量交换的问题。一个冷、热水混合器的内部流动与热量交换的问题。温度为温度为T=350K的热水自上部的热水小管嘴流入，与自下的热水自上部的热水小管嘴流入，与自下部右侧小管嘴流入的温
22、度为部右侧小管嘴流入的温度为290K的冷水在混合器内进行的冷水在混合器内进行热量与动量的交换后，自下部左侧的小管嘴流出。混合器热量与动量的交换后，自下部左侧的小管嘴流出。混合器结构如图所示。结构如图所示。算例算例冷、热水混合后的出流口冷、热水混合后的出流口冷水入口，冷水入口，T=290K, V=10m/s热水入口，热水入口，T=350K, V=10m/s算例算例v 冷、热水混合器内的三维流动与换热冷、热水混合器内的三维流动与换热 冷水和热水分别自混合器的两侧沿水平切向方向流冷水和热水分别自混合器的两侧沿水平切向方向流入，在容器内混合后经过下部渐缩通道流入等径的出流管入，在容器内混合后经过下部渐
23、缩通道流入等径的出流管，最后流入大气。这是一个三维流动问题，所研究的内容，最后流入大气。这是一个三维流动问题，所研究的内容是混合器内的流场、压力分布和温度场。混合器结构如图是混合器内的流场、压力分布和温度场。混合器结构如图所示：所示：算例算例热水入口热水入口速度：速度：1m/s温度：温度：320K冷水入口冷水入口速度：速度：1m/s温度：温度：280K出水口出水口相对压强：相对压强：0PavFluent入门与进阶教程入门与进阶教程：新手入门：新手入门vFluent流体计算应用教程流体计算应用教程：全面：全面v计算流体动力学分析计算流体动力学分析CFD软件原理与应用软件原理与应用：通俗易懂通俗易懂vFluent技术基础与应用实例技术基础与应用实例：含实例：含实例vFluent中英文帮助中英文帮助v流体力学专业论坛流体力学专业论坛 推荐书籍推荐书籍LOGO
