文档分类：
在线文档经过高度压缩，下载原文更清晰。
Fluent经典问题答疑 Fluent 经典问题及答疑1.在 gambit 中对一体积成功的进行了体网格,网格进行了 examine mesh,也没有什么问题,可当要进行边界类型(boundary type)的设定时,却发现 type 只有 node, element_side 两项,没有什么 wall,pressure_outlet 等。为何无法定义边界?===& 因为没有选择求解器为 fluent 5/62.在 FLUENT 模拟以后用 display 下的操作都无法显示,不过刚开始用的是好的,然后就不行了,为什么?===& DirectX 控制面板中的“加速”功能禁用即可3.把带网格的几个 volume,copy 到另一处,但原来 split 的界面,现在都变成了 wall,怎么才能把 wall 变成内部流体呢?===〉直接边界面定义为 interior 即可第 3 题:在数值模拟过程中,离散化的目的是什么?如何对计算区域进行离散化?离散化时通常使用哪些网格?如何对控制方程进行离散?离散化常用的方法有哪些?它们有什么不同?注:我将原题目的提问顺序进行了修改调整,这样更利于回答。4.FLUENT 中常用的文件格式类型:dbs,msh,cas,dat,trn,jou,profile 等有什么用处?在 Gambit 目录中,有三个文件,分别是 default_id.dbs,jou,trn 文件,对 Gambit 运行save,将会在工作目录下保存这三个文件:default_id.dbs,default_id.jou,default_id.trn。jou 文件是 gambit 命令记录文件,可以通过运行 jou 文件来批处理 gambit 命令;dbs 文件是 gambit 默认的储存几何体和网格数据的文件;trn 文件是记录 gambit 命令显示窗(transcript)信息的文件;msh 文件可以在 gambit 划分网格和设置好边界条件之后 export 中选择 msh 文件输出格式,该文件可以被 fluent 求解器读取。Case 文件包括网格,边界条件,解的参数,用户界面和图形环境。Data 文件包含每个网格单元的流动值以及收敛的历史纪录(残差值)。Fluent 自动保存文件类型,默认为 date 和 case 文件Profile 文件边界轮廓用于指定求解域的边界区域的流动条件。例如,它们可以用于指定入口平面的速度场。读入轮廓文件,点击菜单 File/Read/Profile...弹出选择文件对话框,你就可以读入边界轮廓文件了。写入轮廓文件,你也可以在指定边界或者表面的条件上创建轮廓文件。例如:你可以在一个算例的出口条件中创建一个轮廓文件,然后在其它算例中读入该轮廓文件,并使用出口轮廓作为新算例的入口轮廓。要写一个轮廓文件,你需要使用 Write Profile 面板(Figure 1),菜单:File/Write/Profile...是关于 Tecplot 软件使用的;在这里给 Tecplot 新手推荐一个学习的方法:如果你对 Tecplot 一点都不熟悉的话,别紧张,没关系的,你可以直接看 Tecplot 的动画 Demo很快就能入门。Demo 网页:在“开始”-&&所有程序“-&&Tecplot 10&-&&Tutorials&-&&Tutorial&打开网页,里面有 Tecplot 使用的大部分功能的动画教程,花上一二十分钟,就能对 Tecplot有所了解,并可以作初步的使用了。这八道问题的答案基本都在里面了。以下是 Demo 的内容:Let's take care of the warning &turbulent viscosity limited to viscosity ratio****& which is notphysical. This problem is mainly due to one of the following:1)Poor mesh quality(i.e.,skewness & 0.85 for Quad/Hex, or skewness & 0.9 for Tri/Tetra elements).{what values do you have?}2)Use of improper turbulent boudary conditions.3)Not supplying good initial values for turbulent quantities.出现这个警告,一般来讲,最可能的就是网格质量的问题,尤其是 Y+值的问题;在划分网格的时候要注意,第一层网格高度非常重要,可以使用 NASA 的 Viscous Grid Space Calculator来计算第一层网格高度;如果这方面已经注意了,那就可能是边界条件中有关湍流量的设置问题,关于这个,本版中已经有专门的帖子进行了讨论,Fluent 培训的教程中也有讲到,请大家参考。回答问题:4 常见离散格式的性能的对比(稳定性、精度和经济性)。请参考王福军的书《计算流体动力学分析—CFD 理论与应用》离散格式稳定性及稳定条件精度与经济性中心差分条件稳定 Peclet 小于等于 2 在不发生振荡的参数范围内,可以获得校准确的结果。一阶迎风绝对稳定虽然可以获得物理上可接受的解,但当 Peclet 数较大时,假扩散较严重。为避免此问题,常需要加密计算网格。二阶迎风绝对稳定精度较一阶迎风高,但仍有假扩散问题。混合格式绝对稳定当 Peclet 小于等于 2 时,性能与中心差分格式相同。当Peclet 大于 2 时,性能与一阶迎风格式相同。指数格式、乘方格式绝对稳定主要适用于无源项的对流扩散问题,对有非常数源项的场合,当 Peclet 数较高时有较大误差。QUICK 格式条件稳定 Peclet 小于等于8/3可以减少假扩散误差,精度较高,应用较广泛,但主要用于六面体和四边形网格。改进的 QUICK 格式绝对稳定性能同标准 QUICK 格式,只是不存在稳定性问题。5.在利用有限体积法建立离散方程时,必须遵守哪几个基本原则?具体参考王福军的书《计算流体动力学—CFD 软件原理与应用》的第 52-54 页,这里只作简短介绍。在利用有限体积法建立离散方程时,必须遵守如下四条基本原则:1.控制体积界面上的连续性原则;2.正系数原则;3.源项的负斜率线性化原则;4.主系数等于相邻节点系数之和原则。有限体积法的四条基本原则.pdf6 .流场数值计算的目的是什么?主要方法有哪些?其基本思路是什么?各自的适用范围是什么?答:这个问题的范畴好大啊。简要的说一下个人的理解吧:流场数值求解的目的就是为了得到某个流动状态下的相关参数,这样可以节省实验经费,节约实验时间,并且可以模拟一些不可能做实验的流动状态。主要方法有有限差分,有限元和有限体积法,好像最近还有无网格法和波尔兹曼法(格子法)。基本思路都是将复杂的非线性差分/积分方程简化成简单的代数方程。相对来说,有限差分法对网格的要求较高,而其他的方法就要灵活的多第 7 题:可压缩流动和不可压缩流动,在数值解法上各有何特点?为何不可压缩流动在求解时反而比可压缩流动有更多的困难?注:这个问题不是一句两句话就能说清楚的,大家还是看下面的两篇小文章吧,摘自《计算流体力学应用》,读完之后自有体会。8 什么叫边界条件?有何物理意义?它与初始条件有什么关系?边界条件与初始条件是控制方程有确定解的前提。边界条件是在求解区域的边界上所求解的变量或其导数随时间和地点的变化规律。对于任何问题,都需要给定边界条件。初始条件是所研究对象在过程开始时刻各个求解变量的空间分布情况,对于瞬态问题,必须给定初始条件,稳态问题,则不用给定。对于边界条件与初始条件的处理,直接影响计算结果的精度。在瞬态问题中,给定初始条件时要注意的是:要针对所有计算变量,给定整个计算域内各单元的初始条件;初始条件一定是物理上合理的,要靠经验或实测结果。10 在数值计算中,偏微分方程的双曲型方程、椭圆型方程、抛物型方程有什么区别?PS:这个问题想来应该是比较基础的问题,既然没人回答,我就插几句吧;嘿嘿。我们知道很多描述物理问题的控制方程最终就可以归结为偏微分方程,描述流动的控制方程也不例外。从数学角度,一般将偏微分方程分为椭圆型(影响域是椭圆的,与时间无关,且是空间内的闭区域,故又称为边值问题),双曲型(步进问题,但依赖域仅在两条特征区域之间),抛物型(影响域以特征线为分界线,与主流方向垂直;具体来说,解的分布与瞬时以前的情况和边界条件相关,下游的变化仅与上游的变化相关;也称为初边值问题);从物理角度,一般将方程分为平衡问题(或稳态问题),时间步进问题。两种角度,有这样的关系:椭圆型方程描述的一般是平衡问题(或稳态问题),双曲型和抛物型方程描述的一般是步进问题。至于具体的分类方法,大家可以参考一般的偏微分方程专著,里面都有介绍。关于各种不同近似水平的流体控制方程的分类,大家可以参考张涵信院士编写《计算流体力学—差分方法的原理与应用》里面讲的相当详细。三种类型偏微分方程的基本差别如下:1)三种类型偏微分方程解的适定性(即解存在且唯一,并且解稳定)要求对定解条件有不同的提法;2)三种类型偏微分方程解的光滑性不同,对定解条件的光滑性要求也不同;椭圆型和抛物型方程的解是充分光滑的,因此对定解条件的光滑性要求不高。而双曲型方程允许有所谓的弱解存在(如流场中的激波),即解的一阶导数可以不连续,所以对定解条件的光滑性要求很高,这也正是采用有限元法求解双曲型方程困难较多的原因之一。3)三种类型偏微分方程的影响区域和依赖区域不一样。在双曲型和抛物型方程所控制的流场中,某一点的影响区域是有界的,可采用步进求解。如对双曲型方程求解时,为了与影响区域的特征一致,采用上风格式比较适宜。而椭圆型方程的影响范围遍及全场,必须全场求解,所采用的差分格式也要采用相应的中心格式。以上只是一些较为肤浅的概念,如想深入,可参考相关的偏微分方程及数值计算等书籍。可压缩 Euler 及 Navier-Stokes 方程数值解描述无粘流动的基本方程组是 Euler 方程组,描述粘性流动的基本方程组是 Navier-Stokes方程组。用数值方法通过求解 Euler 方程和 Navier-Stokes 方程模拟流场是计算流体动力学的重要内容之一。由于飞行器设计实际问题中的绝大多数流态都具有较高的雷诺数,这些流动粘性区域很小,由对流作用主控,因此针对 Euler 方程发展的计算方法,在大多数情况下对Navier-Stokes 方程也是有效的,只需针对粘性项用中心差分离散。用数值方法求解无粘 Euler 方程组的历史可追溯到 20 世纪 50 年代,具有代表性的方法是 1952 年 Courant 等人以及 1954 年 Lax 和 Friedrichs 提出的一阶方法。从那时开始,人们发展了大量的差分格式。Lax 和 Wendroff 的开创性工作是非定常 Euler(可压缩 Navier-Stokes)方程组数值求解方法发展的里程碑。二阶精度 Lax-Wendroff 格式应用于非线性方程组派生出了一类格式,其共同特点是格式空间对称,即在空间上对一维问题是三点中心格式,在时间上是显式格式,并且该类格式是从时间空间混合离散中导出的。ormack 格式。采用时空混合离散方法,其数值解趋近于定常时依赖于计算中采用的时间步长。尽管由时间步长项引起的误差与截断误差在数量级上相同,但这却体现了一个概念上的缺陷,因为在计算得到的定常解中引进了一个数值参数。将时间积分从空间离散中分离出来就避免了上述缺陷。常用的时空分别离散格式有中心型格式和迎风型格式。空间二阶精度的中心型格式(一维问题是三点格式)就属于上述范畴。该类格式最具代表性的是 Beam-Warming 隐式格式和 Jameson 等人采用的 Runge-Kutta 时间积分方法发展的显式格式。迎风型差分格式共同特点是所建立起的特征传播特性与差分空间离散方向选择的关系是与无粘流动的物理特性一致的。第一个显式迎风差分格式是由 Courant 等人构造的,并推广为二阶精度和隐式时间积分方法。基于通量方向性离散的 Steger-Warming 和 Van Leer 矢通量分裂方法可以认为是这类格式的一种。该类格式的第二个分支是 Godunov 方法,该方法在每个网格步求解描述相邻间断(Riemann 问题)的当地一维 Euler 方程。根据这一方法 Engquist、Osher 和 Roe 等人构造了一系列引入近似 Riemann 算子的格式,这就是著名的通量差分方法。对于没有大梯度的定常光滑流动,所有求解 Euler 方程格式的计算结果都是令人满意的,但当出现诸如激波这样的间断时,其表现确有很大差异。绝大多数最初发展起来的格式,如Lax-Wendroff 格式中心型格式,在激波附近会产生波动。人们通过引入人工粘性构造了各种方法来控制和限制这些波动。在一个时期里,这类格式在复杂流场计算中得到了应用。然而,由于格式中含有自由参数,对不同问题要进行调整,不仅给使用上带来了诸多不便,而且格式对激波分辨率受到影响,因而其在复杂流动计算中的应用受到了一定限制。另外一种方法是力图阻止数值波动的产生,而不是在其产生后再进行抑制。这种方法是建立在非线性限制器的概念上,这一概念最初由 Boris 和 Book 及 Van Leer 提出,并且通过Harten 发展的总变差减小(TVD, Total Variation Diminishing)的重要概念得以实现。通过这一途径,数值解的变化以非线性的方式得以控制。这一类格式的研究和应用,在 20 世纪 80年代形成了一股发展浪潮。1988 年,张涵信和庄逢甘利用热力学熵增原理,通过对差分格式修正方程式的分析,构造了满足熵增条件能够捕捉激波的无波动、无自由参数的耗散格式(NND 格式)。该类格式在航空航天飞行器气动数值模拟方面得到了广泛应用。1987 年,Harten 和 Osher 指出,TVD 格式最多能达到二阶精度。为了突破这一精度上的限制引入了实质上无波动(ENO)格式的概念。该类格式“几乎是 TVD”的,Harten 因此推断这些格式产生的数值解是一致有界的。继 Harten 和 Osher 之后,Shu 和 Osher 将 ENO 格式从一维推广到多维。J.Y.Yang 在三阶精度 ENO 差分格式上也做了不少工作。1992 年,张涵信另辟蹊径,在 NND 格式的基础上,发展了一种能捕捉激波的实质上无波动、无自由参数的三阶精度差分格式(简称 ENN 格式)。1994 年,Liu、Osher 和 Chan 发展了 WENO(WeightedEssentially Non-Oscillatory)格式。WENO 格式是基于 ENO 格式构造的高阶混合格式,它在保持了 ENO 格式优点的同时,计算流场中虚假波动明显减少。此后,Jiang 提出了一种新的网格模板光滑程度的度量方法。目前高阶精度格式的研究与应用是计算流体力学的热点问题之一。不可压缩 Navier-Stokes 方程求解不可压缩流体力学数值解法有非常广泛的需求。从求解低速空气动力学问题,推进器内部流动,到水动力相关的液体流动以及生物流体力学等。满足这么广泛问题的研究,要求有与之相应的较好的物理问题的数学模型以及鲁棒的数值算法。相对于可压缩流动,不可压缩流动的数值求解困难在于,不可压缩流体介质的密度保持常数,而状态方程不再成立,连续方程退化为速度的散度为零的方程。由此,在可压缩流动的计算中可用于求解密度和压力的连续方程在不可压缩流动求解中仅是动量方程的一个约束条件,由此求解不可压缩流动的压力称为一个困难。求解不可压缩流动的各种方法主要在于求解不同的压力过程。目前,主要有两类求解不可压缩流体力学的方法,原始变量方法和非原始变量方法。求解不可压缩流动的原始变量方法是将 Navier-Stokes 方程写成压力和速度的形式,进行直接求解,这种形式已被广为应用。非原始变量方法主要有 Fasel 提出的流函数-涡函数法、Aziz和 Hellums 提出的势函数-涡函数方法。在求解三维流动问题时,上述每一个方法都需要反复求解三个 Possion 方程,非常耗时。原始变量方法可以分为三类:第一种方法是 Harlow和 Welch 首先提出的压力 Possion 方程方法。该方法首先将动量方程推进求得速度场,然后利用 Possion 方程求解压力,这一种方法由于每一时间步上需要求解 Possion 方程,求解非常耗时。第二种方法是 Patanker 和 Spalding 的 SIMPLE(Semi-Implicit Method forPressure-Linked Equation)法,它是通过动量方程求得压力修正项对速度的影响,使其满足速度散度等于零的条件作为压力控制方程。第三种方法是虚拟压缩方法,这一方法是 Chorin于 1967 年提出的。该方法的核心就是通过在连续方程中引入一个虚拟压缩因子,再附加一项压力的虚拟时间导数,使压力显式地与速度联系起来,同时方程也变成了双曲型方程。这样,方程的形式就与求解可压缩流动的方程相似,因此,许多求解可压缩流动的成熟方法都可用于不可压缩流动的求解。目前,由于基于求解压力 Possion 方程的方法非常复杂及耗时,难于求解具体的工程实际问题,因此用此方法解决工程问题的工作越来越少。工程上常用的主要是 SIMPLE 方法和虚拟压缩方法。10. 椭圆形方程求解区中各点上的值是相互影响的,因而各节点上的代数方程必须联立求解,而不能先解得区域中某一部分上的值后再去确定其余地区上的值抛物型方程描写物理上的一类步进问题,这类问题中因变量与时间有关或问题中有类似时间的变量,因而又称初值问题。其求解区域是个开区间,计算时从已知的初值出发,逐步向前推进,依次获得适合于给定边界条件的解.双曲型方程数值求解也是一个步进过程。物理学中的波动过程,空气动力学中的无粘流体稳态超音速流动及无粘流体的非稳态流动,都是双曲型问题。第 11 题:关于贴体网格的定义,随便找本 CFD 的书上面基本上都有介绍,在此就不多言。在此谈谈网格独立性的问题。可以参考帖子:/dvbbs/dispbbs.asp?boardID=61&ID=1208&page=5数值计算的与实验值之间的误差来源只要有这几个:物理模型近似误差(无粘或有粘,定常与非定常,二维或三维等等),差分方程的截断误差及求解区域的离散误差(这两种误差通常统称为离散误差),迭代误差(离散后的代数方程组的求解方法以及迭代次数所产生的误差),舍入误差(计算机只能用有限位存储计算物理量所产生的误差)等等。在通常的计算中,离散误差随网格变细而减小,但由于网格变细时,离散点数增多,舍入误差也随之加大。第 13 题:具体内容可以参考 Gambit Documentation 中的 Quality Type Definitions 章节。判断网格质量的方面有:Area 单元面积,适用于 2D 单元,较为基本的单元质量特征。Aspect Ratio 长宽比,不同的网格单元有不同的计算方法,等于 1 是最好的单元,如正三角形,正四边形,正四面体,正六面体等;一般情况下不要超过 5:1.Diagonal Ratio 对角线之比,仅适用于四边形和六面体单元,默认是大于或等于 1 的,该值越高,说明单元越不规则,最好等于 1,也就是正四边形或正六面体。Edge Ratio 长边与最短边长度之比,大于或等于 1,最好等于 1,解释同上。EquiAngle Skew 通过单元夹角计算的歪斜度,在 0 到 1 之间,0 为质量最好,1 为质量最差。最好是要控制在 0 到 0.4 之间。EquiSize Skew 通过单元大小计算的歪斜度,在 0 到 1 之间,0 为质量最好,1 为质量最差。2D 质量好的单元该值最好在 0.1 以内,3D 单元在 0.4 以内。MidAngle Skew 通过单元边中点连线夹角计算的歪斜度,仅适用于四边形和六面体单元,在0 到 1 之间,0 为质量最好,1 为质量最差。Size Change 相邻单元大小之比,仅适用于 3D 单元,最好控制在 2 以内。Stretch 伸展度。通过单元的对角线长度与边长计算出来的,仅适用于四边形和六面体单元,在 0 到 1 之间,0 为质量最好,1 为质量最差。Taper 锥度。仅适用于四边形和六面体单元,在 0 到 1 之间,0 为质量最好,1 为质量最差。Volume 单元体积,仅适用于 3D 单元,划分网格时应避免出现负体积。Warpage 翘曲。仅适用于四边形和六面体单元,在 0 到 1 之间,0 为质量最好,1 为质量最差。以上只是针对 Gambit 帮助文件的简单归纳,不同的软件有不同的评价单元质量的指标,使用时最好仔细阅读帮助文件。另外,在 Fluent 中的窗口键入:grid quality 然后回车,Fluent 能检查网格的质量,主要有以下三个指标:1.Maxium cell squish: 如果该值等于 1,表示得到了很坏的单元;2.Maxium cell skewness: 该值在 0 到 1 之间,0 表示最好,1 表示最坏;3.Maxium 'aspect-ratio': 1 表示最好。以上的一些只是简略提要,具体的请参考相关资料。15 对于自己的模型,大多数人有这样的想法:我的模型如何来画网格?用什么样的方法最简单?这样做网格到底对不对?我一般首先关注的是我的模型中那一部分的流态情况最复杂,需要加密网格;那一部分得流动情况是我最关注的,需要加密网格。之后根据上述情况开始考虑怎样对模型进行分块:分块的原则,一方面要保证流态复杂得部位要独立进行网格加密,流动平稳的区域要相应加大网格以节省资源,另一方面,要进行保证分块后的各个子模块结构比较规整,尽可能采用结构网格,提高计算精度。最后,进行网格划分,开始关注网格划分得顺序,要先对那些线和面进行网格设定约束,保证各个区域满足不同得网格疏密,之后再对各个不同的体积分块进行划分,我一般采用从内到外,从最关注得体到非关注得体的顺序进行划分16 在两个面的交界线上如果出现网格间距不同的情况时,即两块网格不连续时,怎么样克服这种情况呢?这个问题就是非连续性网格的设置,一般来说就是把两个交接面设置为一对 interface。另外,作此操作可能出现的问题及可供参考的解决方法为:问题:把两个面(其中一个实际是由若干小面组成,将若干小面定义为了 group 了)拼接在一起,也就是说两者之间有流体通过,两个面个属不同的体,网格导入到 fluent 时,使用 interface时出现网格 check 的错误,将 interface 的边界条件删除,就不会发生网格检查的错误,如何将两个面的网格相连?原因:interface 后的两个体的交接面,fluent 以将其作为内部流体处理(非重叠部分默认为wall,合并后网格会在某些地方发生畸变,导致合并失败,也可能准备合并的两个面几何位置有误差,应该准确的在同一几何位置(合并的面大小相等时),在合并之前要合理分块。解决方法:为了避免网格发生畸变(可能一个面上的网格跑到另外的面上了),可以一面网格粗,一面网格细避免; 再者就是通过将一个面的网格直接映射到另一面上的,两个面默认为interior.也可以将网格拼接一起.17 依据实体在 GAMBIT 建模之前简化时,必须遵循哪几个原则?答:最根本的原则就是简化后对实际流动影响不大我觉得在建模前首先要考虑你模型的结构,物理意义上模型是否为轴对称结构或对称结构,如果是的话看能否简化为二维问题,因为二维问题不管从建模上还是求解上都远远方便与三维模型,而且也能达到相应的精度。其次,在有些梯度比较大的地方这些问题不能简化,像有很多拐角的地方往往存在一些集中,这些不能忽略。18 在设置 GAMBIT 边界层类型时需要注意的几个问题:a、没有定义的边界线如何处理?b、计算域内的内部边界如何处理(2D)?答:gambit 默认为 wall,一般情况下可以到 fluent 再修改边界类型。内部边界如果是 split产生的,那么就不需再设定了,如果不是,那么就需要设定为 interface 或者是 internal19 为何在划分网格后,还要指定边界类型和区域类型?常用的边界类型和区域类型有哪些?答:要得到一个问题的定解就需要有定解条件,而边界条件就属于定解条件。也就是说边界条件确定了结果。不同的流体介质有不同的物理属性,也就会得到不同的结果,所以必须指定区域类型。对于 gambit 来说,默认的区域类型是 fluid,所以一般情况下不需要再指定了。第 20 题:何为流体区域(fluid zone)和固体区域(solid zone)?为什么要使用区域的概念?FLUENT 是怎样使用区域的?Fluid Zone 是一个单元组,是求解域内所有流体单元的综合。所激活的方程都要在这些单元上进行求解。向流体区域输入的信息只是流体介质(材料)的类型。对于当前材料列表中没有的材料,需要用户自行定义。注意,多孔介质也当作流体区域对待。Solid Zone 也是一个单元组,只不过这组单元仅用来进行传热计算,不进行任何的流动计算。作为固体处理的材料可能事实上是流体,但是假定其中没有对流发生,固体区域仅需要输入材料类型。Fluent 中使用 Zone 的概念,主要是为了区分分块网格生成,边界条件的定义等等;21 如何监视 FLUENT 的计算结果?如何判断计算是否收敛?在 FLUENT 中收敛准则是如何定义的?分析计算收敛性的各控制参数,并说明如何选择和设置这些参数?解决不收敛问题通常的几个解决方法是什么?可以采用残差控制面板来显示;或者采用通过某面的流量控制;如监控出口上流量的变化;采用某点或者面上受力的监视;涡街中计算达到收敛时,绕流体的面上受的升力为周期交变,而阻力为平缓的直线。怎样判断计算结果是否收敛?1、观察点处的值不再随计算步骤的增加而变化;2、各个参数的残差随计算步数的增加而降低,最后趋于平缓;3、要满足质量守恒(计算中不牵涉到能量)或者是质量与能量守恒(计算中牵涉到能量)。特别要指出的是,即使前两个判据都已经满足了,也并不表示已经得到合理的收敛解了,因为,如果松弛因子设置得太紧,各参数在每步计算的变化都不是太大,也会使前两个判据得到满足。此时就要再看第三个判据了。还需要说明的就是,一般我们都希望在收敛的情况下,残差越小越好,但是残差曲线是全场求平均的结果,有时其大小并不一定代表计算结果的好坏,有时即使计算的残差很大,但结果也许是好的,关键是要看计算结果是否符合物理事实,即残差的大小与模拟的物理现象本身的复杂性有关,必须从实际物理现象上看计算结果。比如说一个全机模型,在大攻角情况下,解震荡得非常厉害,而且残差的量级也总下不去,但这仍然是正确的,为什么呢,因为大攻角下实际流动情形就是这样的,不断有涡的周期性脱落,流场本身就是非定常的,所以解也是波动的,处理的时候取平均就可以呢:)22 什么叫松弛因子?松弛因子对计算结果有什么样的影响?它对计算的收敛情况又有什么样的影响?1、亚松驰(Under Relaxation):所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减,以避免由于差值过大而引起非线性迭代过程的发散。用通用变量来写出时,为松驰因子(Relaxation Factors)。《数值传热学-214》2、FLUENT 中的亚松驰:由于 FLUENT 所解方程组的非线性,我们有必要控制的变化。一般用亚松驰方法来实现控制,该方法在每一部迭代中减少了的变化量。亚松驰最简单的形式为:单元内变量等于原来的值加上亚松驰因子 a 与变化的积, 分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程,包括耦合解算器所解的非耦合方程(湍流和其他标量)都会有一个相关的亚松驰因子。在 FLUENT 中,所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题,但是对于一些特殊的非线性问题(如:某些湍流或者高 Rayleigh 数自然对流问题),在计算开始时要慎重减小亚松驰因子。使用默认的亚松驰因子开始计算是很好的习惯。如果经过 4 到 5 步的迭代残差仍然增长,你就需要减小亚松驰因子。有时候,如果发现残差开始增加,你可以改变亚松驰因子重新计算。在亚松驰因子过大时通常会出现这种情况。最为安全的方法就是在对亚松驰因子做任何修改之前先保存数据文件,并对解的算法做几步迭代以调节到新的参数。最典型的情况是,亚松驰因子的增加会使残差有少量的增加,但是随着解的进行残差的增加又消失了。如果残差变化有几个量级你就需要考虑停止计算并回到最后保存的较好的数据文件。注意:粘性和密度的亚松驰是在每一次迭代之间的。而且,如果直接解焓方程而不是温度方程(即:对 PDF 计算),基于焓的温度的更新是要进行亚松驰的。要查看默认的亚松弛因子的值,你可以在解控制面板点击默认按钮。对于大多数流动,不需要修改默认亚松弛因子。但是,如果出现不稳定或者发散你就需要减小默认的亚松弛因子了,其中压力、动量、k 和 e 的亚松弛因子默认值分别为 0.2,0.5,0.5 和 0.5。对于 SIMPLEC 格式一般不需要减小压力的亚松弛因子。在密度和温度强烈耦合的问题中,如相当高的 Rayleigh 数的自然或混合对流流动,应该对温度和/或密度(所用的亚松弛因子小于 1.0)进行亚松弛。相反,当温度和动量方程没有耦合或者耦合较弱时,流动密度是常数,温度的亚松弛因子可以设为1.0。对于其它的标量方程,如漩涡,组分,PDF 变量,对于某些问题默认的亚松弛可能过大,尤其是对于初始计算。你可以将松弛因子设为 0.8 以使得收敛更容易。SIMPLE 与 SIMPLEC 比较在 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 邻近校正中所用的方法。24 这个问题的意思是出现了回流,这个问题相对于湍流粘性比的警告要宽松一些,有些 case可能只在计算的开始阶段出现这个警告,随着迭代的计算,可能会消失,如果计算一段时间之后,警告消失了,那么对计算结果是没有什么影响的,如果这个警告一直存在,可能需要作以下处理:1.如果是模拟外部绕流,出现这个警告的原因可能是边界条件取得距离物体不够远,如果边界条件取的足够远,该处可能在计算的过程中的确存在回流现象;对于可压缩流动,边界最好取在 10 倍的物体特征长度之处;对于不可压缩流动,边界最好取在 4 倍的物体特征长度之处。2.如果出现了这个警告,不论对于外部绕流还是内部流动,可以使用 pressure-outlet 边界条件代替 outflow 边界条件改善这个问题。25 燃烧过程中经常遇到一个“头疼”问题是计算后温度场没什么变化?即点火问题,解决计算过程中点火的方法有哪些?什么原因引起点火困难的问题?这个问题确实比较常见,大概说一下自己的解决方法。1) 点火问题。一般来说,对简单的问题我采取初始化时给一个较高温度的办法,当然,PATCH 也是很好的办法。对那些希望局部燃烧的问题,我是画网格的时候予以区分,然后关闭反应项。2)不着火的问题。原因比较复杂,可能是点火不成功,也可能是其他原因。对一个甲烷燃烧问题,或许可以修改活化能来解决。另外就是对具体问题需要选合适的燃烧模型。计算后温度场没有变化,说明化学反应根本没有发生!fluent 解决计算过程中的点火问题,是采用&patch&的方法在已经初始化的流场中设定一个高温的点火区,该选项中设定的温度值一定要高于燃料的最低点火温度(例如甲烷最低点火温度为 800K,这里就可以将其设定为 2000K 或 3000K),其余参数都可以和初始化时一样。不要担心点火区的高温会对整个计算的影响!在燃料点火和燃烧中,最关键的是燃烧反应机理,fluent 中可以采用一步或两步总体反应(fluent 自带)或通过 import 导入正确的详细反应机理。在具体设置过程中,一定要注意:播放器加载中，请稍候...
该用户其他文档
文档介绍：
Fluent经典问题答疑 Fluent 经典问题及答疑1.在 gambit 中对一体积成功的进行了体网格,网格进行了 examine mesh,也没有什么问题,可当要进行边界类型(boundary type)的设定时,却发现 type 只有 node, element_side 两项,没有什么 wall,pressure_outlet 等。为何无法定义边界?===& 因为没有选择求解器为 fluent 5/62.在 FL...
内容来自淘豆网转载请标明出处.
时间： 23:58