FLUENT模拟残差曲线震荡与收敛问题Word文档格式.docx
《FLUENT模拟残差曲线震荡与收敛问题Word文档格式.docx》由会员分享，可在线阅读，更多相关《FLUENT模拟残差曲线震荡与收敛问题Word文档格式.docx（14页珍藏版）》请在冰点文库上搜索。
23
在FLUENT运行过程中，经常会出现“turbulenceviscousrate”超过了极限值，此时如何解决？
而这里的极限值指的是什么值？
修正后它对计算结果有何影响？
Let'
stakecareofthewarning"
turbulentviscositylimitedtoviscosityratio****"
whichisnotphysical.Thisproblemismainlyduetooneofthefollowing:
1）Poormeshquality（i.e.,skewness>
0.85forQuad/Hex,orskewness>
0.9forTri/Tetraelements）.{whatvaluesdoyouhave?
}
2）Useofimproperturbulentboudaryconditions.
3）Notsupplyinggoodinitialvaluesforturbulentquantities.
出现这个警告，一般来讲，最可能的就是网格质量的问题，尤其是Y值的问题；
在划分网格的时候要注意，第一层网格高度非常重要，可以使用NASA的ViscousGridSpaceCalculator来计算第一层网格高度；
如果这方面已经注意了，那就可能是边界条件中有关湍流量的设置问题，关于这个，本版中已经有专门的帖子进行了讨论，Fluent培训的教程中也有讲到，请大家参考。
24
在FLUENT运行计算时，为什么有时候总是出现“reversedflow”？
其具体意义是什么？
有没有办法避免？
如果一直这样显示，它对最终的计算结果有什么样的影响？
这个问题的意思是出现了回流，这个问题相对于湍流粘性比的警告要宽松一些，有些case可能只在计算的开始阶段出现这个警告，随着迭代的计算，可能会消失，如果计算一段时间之后，警告消失了，那么对计算结果是没有什么影响的，如果这个警告一直存在，可能需要作以下处理：
1.如果是模拟外部绕流，出现这个警告的原因可能是边界条件取得距离物体不够远，如果边界条件取的足够远，该处可能在计算的过程中的确存在回流现象；
对于可压缩流动，边界最好取在10倍的物体特征长度之处；
对于不可压缩流动，边界最好取在4倍的物体特征长度之处。
2.如果出现了这个警告，不论对于外部绕流还是部流动，可以使用pressure-outlet边界条件代替outflow边界条件改善这个问题。
22什么叫松弛因子？
松弛因子对计算结果有什么样的影响？
它对计算的收敛情况又有什么样的影响？
1、亚松驰（UnderRelaxation）：
所谓亚松驰就是将本层次计算结果与上一层次结果的差值作适当缩减，以避免由于差值过大而引起非线性迭代过程的发散。
用通用变量来写出时，为松驰因子（RelaxationFactors）。
《数值传热学-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邻近校正中所用的方法。
1对于刚接触到FLUENT新手来说，面对铺天盖地的学习资料和令人难读的FLUENThelp，如何学习才能在最短的时间入门并掌握基本学习方法呢？
学习任何一个软件，对于每一个人来说，都存在入门的时期。
认真勤学是必须的，什么是最好的学习方法，我也不能妄加定论，在此，我愿意将我三年前入门FLUENT心得介绍一下，希望能给学习FLUENT的新手一点帮助。
由于当时我需要学习FLUENT来做毕业设计，老师给了我一本书，占忠的《FLUENT流体工程仿真计算实例与应用》，当然，学这本书之前必须要有两个条件，第一，具有流体力学的基础，第二，有FLUENT安装软件可以应用。
然后就照着书上二维的计算例子，一个例子，一个步骤地去学习，然后学习三维，再针对具体你所遇到的项目进行针对性的计算。
不能急于求成，从前处理器GAMBIT，到通过FLUENT进行仿真，再到后处理，如TECPLOT，进行循序渐进的学习，坚持，效果是非常显著的。
如果身边有懂得FLUENT的老师，那么遇到问题向老师请教是最有效的方法，碰到不懂的问题也可以上网或者查找相关书籍来得到答案。
另外我还有本《计算流体动力学分析》王福军的，两者结合起来学习效果更好。
2CFD计算中涉及到的流体及流动的基本概念和术语：
理想流体和粘性流体；
牛顿流体和非牛顿流体；
可压缩流体和不可压缩流体；
层流和湍流；
定常流动和非定常流动；
亚音速与超音速流动；
热传导和扩散等。
A.理想流体（IdealFluid）和粘性流体（ViscousFluid）：
流体在静止时虽不能承受切应力，但在运动时，对相邻的两层流体间的相对运动，即相对滑动速度却是有抵抗的，这种抵抗力称为粘性应力。
流体所具备的这种抵抗两层流体相对滑动速度，或普遍说来抵抗变形的性质称为粘性。
粘性的大小依赖于流体的性质，并显著地随温度变化。
实验表明，粘性应力的大小与粘性及相对速度成正比。
当流体的粘性较小（实际上最重要的流体如空气、水等的粘性都是很小的），运动的相对速度也不大时，所产生的粘性应力比起其他类型的力如惯性力可忽略不计。
此时我们可以近似地把流体看成无粘性的，这样的流体称为理想流体。
十分明显，理想流体对于切向变形没有任何抗拒能力。
这样对于粘性而言，我们可以将流体分为理想流体和粘性流体两大类。
应该强调指出，真正的理想流体在客观实际中是不存在的，它只是实际流体在某些条件下的一种近似模型。
B.牛顿流体（NewtonianFluid）和非牛顿流体（non-NewtonianFluid）：
日常生活和工程实践中最常遇到的流体其切应力与剪切变形速率符合下式的线性关系，称为牛顿流体。
而切应力与变形速率不成线性关系者称为非牛顿流体。
图2-1（a）中绘出了切应力与变形速率的关系曲线。
其中符合上式的线性关系者为牛顿流体。
其他为非牛顿流体，非牛顿流体中又因其切应力与变形速率关系特点分为膨胀性流体（Dilalant），拟塑性流体（Pseudoplastic），具有屈服应力的理想宾厄流体（IdealBinghamFluid）和塑性流体（PlasticFluid）等。
通常油脂、油漆、牛奶、牙膏、血液、泥浆等均为非牛顿流体。
非牛顿流体的研究在化纤、塑料、石油、化工、食品及很多轻工业中有着广泛的应用。
图2-1（b）还显示出对于有些非牛顿流体，其粘滞特性具有时间效应，即剪切应力不仅与变形速率有关而且与作用时间有关。
当变形速率保持常量，切应力随时间增大，这种非牛顿流体称为震凝性流体（RheopecticFluid）。
当变形速率保持常量而切应力随时间减小的非牛顿流体则称为触变性流体（ThixotropicFluid）。
C.可压缩流体（CompressibleFluid）和不可压缩流体（IncompressibleFluid）：
在流体的运动过程中，由于压力、温度等因素的改变，流体质点的体积（或密度，因质点的质量一定），或多或少有所改变。
流体质点的体积或密度在受到一定压力差或温度差的条件下可以改变的这个性质称为压缩性。
真实流体都是可以压缩的。
它的压缩程度依赖于流体的性质及外界的条件。
例如水在100个大气压下，容积缩小0.5%，温度从20°
变化到100°
，容积降低4%。
因此在一般情况下液体可以近似地看成不可压的。
但是在某些特殊问题中，例如水中爆炸或水击等问题，则必须把液体看作是可压缩的。
气体的压缩性比液体大得多，所以在一般情形下应该当作可压缩流体处理。
但是如果压力差较小，运动速度较小，并且没有很大的温度差，则实际上气体所产生的体积变化也不大。
此时，也可以近似地将气体视为不可压缩的。
在可压缩流体的连续方程中含密度，因而可把密度视为连续方程中的独立变量进行求解，再根据气体的状态方程求出压力。
不可压流体的压力场是通过连续方程间接规定的。
由于没有直接求解压力的方程，不可压流体的流动方程的求解具有其特殊的困难。
D.层流（LaminarFlow）和湍流（TurbulentFlow）：
实验表明，粘性流体运动有两种形态，即层流和湍流。
这两种形态的性质截然不同。
层流是流体运动规则，各部分分层流动互不掺混，质点的轨线是光滑的，而且流动稳定。
湍流的特征则完全相反，流体运动极不规则，各部分激烈掺混，质点的轨线杂乱无章，而且流场极不稳定。
这两种截然不同的运动形态在一定条件下可以相互转化。
E.定常流动（SteadyFlow）和非定常流动（UnsteadyFlow）：
以时间为标准，根据流体流动的物理量（如速度、压力、温度等）是否随时间变化，将流动分为定常与非定常两大类。
当流动的物理量不随时间变化，为定常流动；
反之称为非定常流动。
定常流动也称为恒定流动，或者稳态流动；
非定常流动也称为非恒定流动、非稳态流动。
许多流体机械在起动或关机时的流体流动一般是非定常流动，而正常运转时可看作是定常流动。
F.亚音速流动（Subsonic）与超音速流动（Supersonic）：
当气流速度很大，或者流场压力变化很大时，流体就受到了压速性的影响。
马赫数定义为当地速度与当地音速之比。
当马赫数小于1时，流动为亚音速流动；
当马赫数远远小于1（如M<
0.1）时，流体的可压速性及压力脉动对密度变化影响都可以忽略。
当马赫数接近1时候（跨音速），可压速性影响就显得十分重要了。
如果马赫数大于1，流体就变为超音速流动。
FLUENT对于亚音速，跨音速以及超音速等可压流动都有模拟能力。
G.热传导（HeatTransfer）及扩散（Diffusion）：
除了粘性外，流体还有热传导及扩散等性质。
当流体中存在温度差时，温度高的地方将向温度低的地方传送热量，这种现象称为热传导。
同样地，当流体混合物中存在组元的浓度差时，浓度高的地方将向浓度低的地方输送该组元的物质，这种现象称为扩散。
流体的宏观性质，如扩散、粘性和热传导等，是分子输运性质的统计平均。
由于分子的不规则运动，在各层流体间交换着质量、动量和能量，使不同流体层的平均物理量均匀化，这种性质称为分子运动的输运性质。
质量输运宏观上表现为扩散现象，动量输运表现为粘性现象，能量输运表象为热传导现象。
理想流体忽略了粘性，即忽略了分子运动的动量输运性质，因此在理想流体中也不应考虑质量和能量输运性质——扩散和热传导，因为它们具有相同的微观机制
3在数值模拟过程中，离散化的目的是什么？
如何对计算区域进行离散化？
离散化时通常使用哪些网格？
如何对控制方程进行离散？
离散化常用的方法有哪些？
它们有什么不同？
首先说一下CFD的基本思想：
把原来在时间域及空间域上连续的物理量的场，如速度场，压力场等，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。
然后，我们再讨论下这些题目。
离散化的目的：
我们知道描述流体流动及传热等物理问题的基本方程为偏微分方程，想要得它们的解析解或者近似解析解，在绝大多数情况下都是非常困难的，甚至是不可能的，就拿我们熟知的Navier-Stokes方程来说，现在能得到的解析的特解也就70个左右；
但为了对这些问题进行研究，我们可以借助于我们已经相当成熟的代数方程组求解方法，因此，离散化的目的简而言之，就是将连续的偏微分方程组及其定解条件按照某种方法遵循特定的规则在计算区域的离散网格上转化为代数方程组，以得到连续系统的离散数值逼近解。
计算区域的离散及通常使用的网格：
在对控制方程进行离散之前，我们需要选择与控制方程离散方法相适应的计算区域离散方法。
网格是离散的基础，网格节点是离散化的物理量的存储位置，网格在离散过程中起着关键的作用。
网格的形式和密度等，对数值计算结果有着重要的影响。
一般情况下，二维问题，有三角形单元和四边形，三位问题中，有四面体，六面体，棱锥体，楔形体及多面体单元。
网格按照常用的分类方法可以分为：
结构网格，非结构网格，混合网格；
也可以分为：
单块网格，分块网格，重叠网格；
等等。
上面提到的计算区域的离散方法要考虑到控制方程的离散方法，比如说：
有限差分法只能使用结构网格，有限元和有限体积法可以使用结构网格也可以使用非结构网格。
控制方程的离散及其方法：
上面已经提到了离散化的目的，控制方程的离散就是将主控的偏微分方程组在计算网格上按照特定的方法离散成代数方程组，用以进行数值计算。
按照应变量在计算网格节点之间的分布假设及推到离散方程的方法不同，控制方程的离散方法主要有：
有限差分法，有限元法，有限体积法，边界元法，谱方法等等。
这里主要介绍最常用的有限差分法，有限元法及有限体积法。
（1）有限差分法（FiniteDifferenceMethod，简称FDM）是数值方法中最经典的方法。
它是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程（控制方程）的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。
求差分方程组（代数方程组）的解，就是微分方程定解问题的数值近似解，这是一种直接将微分问题变为代数问题的近似数值解法。
这种方法发展较早，比较成熟，较多用于求解双曲型和抛物型问题（发展型问题）。
用它求解边界条件复杂，尤其是椭圆型问题不如有限元法或有限体积法方便。
（2）有限元法（FiniteElementMethod，简称FEM）与有限差分法都是广泛应用的流体力学数值计算方法。
有限元法是将一个连续的求解域任意分成适当形状的许多微小单元，并于各小单元分片构造插值函数，然后根据极值原理（变分或加权余量法），将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为个单元极值之和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组就得到各节点上待求的函数值。
有限元法的基础是极值原理和划分插值，它吸收了有限差分法中离散处理的核，又采用了变分计算中选择逼近函数并对区域积分的合理方法，是这两类方法相互结合，取长补短发展的结果。
它具有广泛的适应性，特别适用于几何及物理条件比较复杂的问题，而且便于程序的标准化。
对椭圆型问题（平衡态问题）有更好的适应性。
有限元法因求解速度较有限差分法和有限体积法满，因此，在商用CFD软件中应用并不普遍，目前的商用CFD软件中，FIDAP采用的是有限元法。
而有限元法目前在固体力学分析中占绝对比例，几乎所有的固体力学分析软件都是采用有限元法。
（3）有限体积法（FiniteVolumeMethod，简称FVM）是近年发展非常迅速的一种离散化方法，其特点是计算效率高。
目前在CFD领域得到了广泛的应用。
其基本思路是：
将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积；
将待解的微分方程（控制方程）对每一个控制体积分，从而得到一组离散方程。
其中的未知数是网格点上的因变量，为了求出控制体的积分，必须假定因变量值在网格点之间的变化规律。
从积分区域的选取方法看来，有限体积法属于加权余量法中的子域法，从未知解的近似方法看来，有限体积法属于采用局部近似的离散方法。
简言之，子域法加离散，就是有限体积法的基本方法。
各种离散化方法的区别：
简短而言，有限元法，将物理量存储在真实的网格节点上，将单元看成由周边节点及型函数构成的统一体；
有限体积法往往是将物理量存储在网格单元的中心点上，而将单元看成围绕中心点的控制体积，或者在真实网格节点上定义和存储物理量，而在节点周围构造控制题。
4常见离散格式的性能的对比（稳定性、精度和经济性）
请参考王福军的书《计算流体动力学分析—CFD理论与应用》
离散格式
稳定性及稳定条件
精度与经济性
中心差分
条件稳定Peclet小于等于2
在不发生振荡的参数围，可以获得校准确的结果。
一阶迎风
绝对稳定
虽然可以获得物理上可接受的解，但当Peclet数较大时，假扩散较严重。
为避免此问题，常需要加密计算网格。
二阶迎风
精度较一阶迎风高，但仍有假扩散问题。
混合格式
当Peclet小于等于2时，性能与中心差分格式相同。
当Peclet大于2时，性能与一阶迎风格式相同。
指数格式、乘方格式
主要适用于无源项的对流扩散问题，对有非常数源项的场合，当Peclet数较高时有较大误差。
QUICK格式
条件稳定Peclet小于等于8/3
可以减少假扩散误差，精度较高，应用较广泛，但主要用于六面体和四边形网格。
改进的QUICK格式
性能同标准QUICK格式，只是不存在稳定性问题。
5在利用有限体积法建立离散方程时，必须遵守哪几个基本原则？
1.控制体积界面上的连续性原则；
2.正系数原则；
3.源项的负斜率线性化原则；
4.主系数等于相邻节点系数之和原则。
6流场数值计算的目的是什么？
主要方法有哪些？
其基本思路是什么？
各自的适用围是什么？
这个问题的畴好大啊。
简要的说一下个人的理解吧：
流场数值求解的目的就是为了得到某个流动状态下的相关参数，这样可以节省实验经费，节约实验时间，并且可以模拟一些不可能做实验的流动状态。
主要方法有有限差分，有限元和有限体积法，好像最近还有无网格法和波尔兹曼法（格子法）。
基本思路都是将复杂的非线性差分/积分方程简化成简单的代数方程。
相对来说，有限差分法对网格的要求较高，而其他的方法就要灵活的多
7可压缩流动和不可压缩流动，在数值解法上各有何特点？
为何不可压缩流动在求解时反而比可压缩流动有更多的困难？
可压缩Euler及Navier-Stokes方程数值解
描述无粘流动的基本方程组是Euler方程组，描述粘性流动的基本方程组是Navier-Stokes方程组。
用数值方法通过求解Euler方程和Navier-Stokes方程模拟流场是计算流体动力学的重要容之一。
由于飞行器设计实际问题中的绝大多数流态都具有较高的雷诺数，这些流动粘性区域很小，由对流作用主控，因此针对Euler方程发展的计算方法，在大多数情况下对Navier-Stokes方程也是有效的，只需针对粘性项用中心差分离散。
用数值方法求解无粘Euler方程组的
升级会员 