fluent第一章概述

中国科学技术大学FLUENT讲义：第一章概论--1--FLUENT简介FLUENT是用于计算复杂几何条件下流动和传热问题的程序。它提供的无结构网格生成程序，把计算相对复杂的几何结构问题变得容易和轻松。可以生成的网格包括二维的三角形和四边形网格；三维的四面体、六面体及混合网格。并且，可以根据计算结果调整网格。这种网格的自适应能力对于精确求解有较大梯度的流场如自由剪切流和边界层问题有很实际的作用。同时，网格自适应和调整只是在需要加密的流动区域里实施，而非整个流动场，因此可以节约计算时间。程序的结构FLUENT程序软件包应该包括以下几个部分：1，FLUENT解法器2，prePDF，用于模拟PDF燃烧过程3，GAMBIT，网格生成4，TGrid，额外的处理器，用于从现有的边界网格生成体网格。5，Filters(Translators)，转换其它程序生成的网格，用于FLUENT计算。可以接口的程序包括：ANSYS,I-DEAS,NASTRAN,PATRAN等。基本程序结构示意图。GAMBIT设置几何形状生成2D或3D网格其它软件包，如CAD，CAE等FLUENT网格输入及调整物理模型边界条件流体物性确定计算结果后处理TGrid2D三角网格3D四面体网格2D和3D混合网格prePDFPDF查表2D或3D网格几何形状或网格PDF程序网格边界和（或）体网格边界网格中国科学技术大学FLUENT讲义：第一章概论--2--FLUENT程序的用途1，采用三角形、四边形、四面体、六面体及其混合网格计算二维和三位流动问题；计算过程中，网格可以自适应。2，可压缩与不可压缩流动问题；3，稳态和瞬态流动问题；4，无粘流，层流及湍流问题；5，牛顿流体及非牛顿流体；6，对流换热问题（包括自然对流和混合对流）；7，导热与对流换热耦合问题；8，辐射换热；9，惯性坐标系和非惯性坐标系下的流动问题模拟；10，多运动坐标系下的流动问题；11，化学组分混合与反应；12，可以处理热量、质量、动量和化学组分的源项；13，用Lagrangian轨道模型模拟稀疏相（颗粒，水滴，气泡等）；14，多孔介质流动；15，一维风扇、热交换器性能计算；16，两相流问题；17，复杂表面形状下的自由面流动；二维网格：三维网格：图1－1，FLUENT的基本控制体形状用FLUENT程序求解问题的步骤1，确定几何形状，生成计算网格（用GAMBIT，也可以读入其它指定程序生成的网格）；2，选择2D或3D来模拟计算；3，输入网格；tetrahedrontrianglequadrilateralhexahedronpyramidprismorwedge中国科学技术大学FLUENT讲义：第一章概论--3--4，检查网格；5，选择解法器；6，选择求解的方程：层流或湍流（或无粘流），化学组分或化学反应，传热模型等。确定其它需要的模型如：风扇、热交换器、多孔介质等模型。7，确定流体物性；8，指定边界条件；9，条件计算控制参数；10，流场初始化；11，计算；12，检查结果；13，保存结果，后处理等。关于FLUENT解法器的说明1，FLUENT2D，二维单精度解法器2，FLUENT3D，三维单精度解法器3，FLUENT2ddp，二维双精度解法器4，FLUENT3ddp，三维双精度解法器FLUENT程序启动方法1，WINDOEWSNT下，点击FLUENT5。2，在MS－DOS下，键入命令。FLUENT2D/3D/2ddp/3ddp。平行计算命令为：FLUENT2D/3D/2ddp/3ddp-tx。x是处理器编号。如，我们用3号处理器计算三维双精度问题，命令为：FLUENT3ddp–t3FLUENT命令的一般形式为：FLUENT[version][-help][options]FLUENT求解方法的选择1．非耦合求解(Segregated)2．耦合隐式求解(CoupledImplicit)3．耦合显式求解(CoupledExplicit)非耦合求解方法主要用于不可压缩或压缩性不强的流体流动。耦合求解则可以用在高速可压缩流动。FLUENT默认设置是非耦合求解，但对于高速可压流动，有强的体积力（浮力或离心力）的流动，求解问题时网格要比较密，建议采用耦合隐式求解方法，可以耦合求解能量和动量方程，能比较快地得到收敛解。缺点是需要的内存比较大（是非耦合求解迭代时间的1.5-2倍）。如果必须要耦合求解，但是你的机器内存不够，这时候可以考虑用耦合显式解法器求解问题。该解法器也耦合了动量，能量及组分方程，但内存却比隐式求解方法小。缺点是收敛时间比较长。这里需要指出的是非耦合求解的一些模型在耦合求解解法器里并不都有。耦合解法器没有的模型包括：多相流模型，混合分数/PDF燃烧模型，预混燃烧模型，污染物生成模型，相变模型，Rosseland辐射模型，确定质量流率的周期性流动模型及周期性换热模型等。中国科学技术大学FLUENT讲义：第一章概论--4--隐式(Implicit)：对于给定变量，单元内的未知值用邻近单元的已知和未知值计算得出。因此，每一个未知值会在不止一个方程中出现，这些方程必须同时解来给出未知量。显式(Explicit)：对于给定变量，每一个单元内的未知量用只包含已知量的关系式计算得到。因此未知量只在一个方程中出现，而且每一个单元内的未知量的方程只需解一次就可以给出未知量的值。一阶迎风格式(FirstOrderUpwind)：当需要一阶精度时，我们假定描述单元内变量平均值的单元中心变量就是整个单元内各个变量的值，而且单元表面的量等于单元内的量。因此，当选择一阶迎风格式时，表面值被设定等于迎风单元的单元中心值。二阶迎风格式(SecondOrderUpwind)：当需要二阶精度时，使用多维线性重建方法来计算单元表面处的值。在这种方法中，通过单元中心解在单元中心处的泰勒展开来实现单元表面的二阶精度值。因此，当使用二阶迎风格式时，用下面的方程来计算表面值；QUICK格式：对于四边形和六面体网格，我们可以确定它们唯一的上游和下游表面以及单元。FLUENT还提供了计算对流变量在表面处高阶值的QUICK格式。QUICK类型的格式是通过变量的二阶迎风与中心插值加上适当的权因子得到的；亚松驰(Under-Relaxation)：由于FLUENT所解方程组的非线性，我们有必要控制的变化。一般用亚松驰方法来实现控制，该方法在每一部迭代中减少了f的变化量。亚松驰最简单的形式为：单元内变量f等于原来的值f_old加上亚松驰因子a与f变化的积：oldSIMPLE：SIMPLE算法使用压力和速度之间的相互校正关系来强制质量守恒并获取压力场。一阶与二阶的比较当流动和网格成一条线时（如：矩形网格或者六面体网格模拟矩形导管的层流流动），可以使用一阶迎风离散格式。但是，当流动和网格不在一条线上时（即：流动斜穿网格线）一阶对流离散增加了对流离散的误差（数值耗散）。对于三角形和四面体网格，流动从来就不会和网格成一条线，此时一般要使用二阶离散来获取更高精度的结果。对于四边形或者六面体网格，如果使用二阶离散格式，尤其是对于复杂流动来说，你可以获取更好的结果。总而言之，一阶离散一般会比二阶离散收敛得好，但是精度要差，尤其是对于三角形或者四面体网格精度更差。对于大多数情况，你可以在计算的开始使用二阶格式。对于有些情况，你应该以一阶离散开始计算，在进行了初步迭代之后再转到二阶格式。例如，如果你解高马赫数流动问题，初始解科所预期的解相差较大，你就应该先用一阶格式迭代几步然后打开二阶格式继续计算直至收敛。对于与网格成一条线的简单流动（如：划分为矩形网格或者六面体网格的矩形导管的层流流动），数值耗散自然会很低，所以一般使用一阶格式替代二阶格式而不损失精度。最后，如果你使用二阶格式遇到收敛性问题，你就应该尝试使用一阶格式。选择压力插值格式如压力插值格式所述，当使用分离求解器时我们可以采用很多压力插值格式。对于大多数情况，标准格式已经足够了，但是对于特定的某些模型使用其它格式可能会更好：对于具有较大体积力的问题，推荐使用体积力加权格式。对于具有高涡流数，高Rayleigh数自然对流，高速旋转流动，包含多孔介质的流动和高度扭曲区域的流动，使用PRESTO!格式。注意：PRESTO！只能用于四边形或者六面体网格。中国科学技术大学FLUENT讲义：第一章概论--5--对于可压流动推荐使用二阶格式。当其它格式不适用时，使用二阶格式来提高精度（如：对于流过具有非六面体或者非四边形网格的曲面边界的流动。）选择压力速度耦合方法在分离求解器中，FLUENT提供了压力速度耦合的三种方法：SIMPLE，SIMPLEC以及PISO。定常状态计算一般使用SIMPLE或者SIMPLEC方法，对于过渡计算推荐使用PISO方法。PISO方法还可以用于高度倾斜网格的定常状态计算和过渡计算。需要注意的是压力速度耦合只用于分离求解器，对于耦合求解器你不可以使用它。SIMPLE与SIMPLEC比较在FLUENT中，可以使用标准SIMPLE算法和SIMPLEC（SIMPLE-Consistent）算法，默认是SIMPLE算法，但是对于许多问题如果使用SIMPLEC可能会得到更好的结果，尤其是可以应用增加的亚松驰迭代时，具体介绍如下：对于相对简单的问题（如：没有附加模型激活的层流流动），其收敛性已经被压力速度耦合所限制，你通常可以用SIMPLEC算法很快得到收敛解。在SIMPLEC中，压力校正亚松驰因子通常设为1.0，它有助于收敛。但是，在有些问题中，将压力校正松弛因子增加到1.0可能会导致不稳定。对于这种情况，你需要使用更为保守的亚松驰或者使用SIMPLE算法。对于包含湍流和/或附加物理模型的复杂流动，只要用压力速度耦合做限制，SIMPLEC会提高收敛性。它通常是一种限制收敛性的附加模拟参数，在这种情况下，SIMPLE和SIMPLEC会给出相似的收敛速度。设定亚松驰因子分离求解器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着，使用分离求解器解的方程，包括耦合求解器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。在FLUENT中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是对于一些特殊的非线性问题（如：某些湍流或者高Rayleigh数自然对流问题），在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如果经过4到5步的迭代残差仍然增长，你就需要减小亚松驰因子。有时候，如果发现残差开始增加，你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现这种情况。最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件，并对解的算法做几步迭代以调节到新的参数。最典型的情况是，亚松驰因子的增加会使残差有少量的增加，但是随着解的进行残差的增加又消失了。如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。对于大多数流动，不需要修改默认亚松弛因子。但是，如果出现不稳定或者发散你就需要减小默认的亚松弛因子了，其中压力、动量、k和e的亚松弛因子默认值分别为0.2，0.5，0.5和0.5。对于SIMPLEC格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合的问题中，如相当高的Rayleigh数的自然或混合对流流动，应该对温度和/或密度（所用的亚松弛因子小于1.0）进行亚松弛。相反，当温度和动量方程没有耦合或者耦合较弱时，流动密度是常数，温度的亚松弛因子可以设为1.0。对于其它的标量方程，如漩涡，组分，PDF变量，对于某些问题默认的亚松弛可能过大，尤其是对于初始计算。你可以将松弛因子设为0.8以使得收敛更容易。中国科学技术大学FLUENT讲义：第一章概论--6--FLUENT求解过程FLUENT的两种解法都可以解守恒型积分方程，其中包括动量、能量、质量以及其他标量如湍流和化学组分的守恒。在两种情况下都应用了控制体技术，它包括：使用计算网格对流体区域进行划分；离散方程的线化以及获取线性方程结果以更新相关变量的值。